Phylogenetic analysis of Fritillaria cirrhosa D. Don and its closely related species based on complete chloroplast genomes

Chloroplast Genome Comparison and Phylogenetic Analysis of the Commercial Variety Actinidia chinensis 'Hongyang'

Actinidia chinensis 'Hongyang', also known as red yangtao (red heart kiwifruit), is a vine fruit tree native to China possessing significant nutritional and economic value. However, information on its genetic diversity and phylogeny is still very limited. The first chloroplast (cp) genome of A. chinensis 'Hongyang' cultivated in China was sequenced using de novo technology in this study. A. chinensis 'Hongyang' possesses a cp genome that spans 156,267 base pairs (bp), exhibiting an overall GC content of 37.20%. There were 132 genes that were annotated, with 85 of them being protein-coding genes, 39 transfer RNA (tRNA) genes, and 8 ribosomal RNA (rRNA) genes. A total of 49 microsatellite sequences (SSRs) were detected, mainly single nucleotide repeats, mostly consisting of A or T base repeats. Compared with 14 other species, the cp genomes of A. chinensis 'Hongyang' were biased towards the use of codons containing A/U, and the non-protein coding regions in the A. chinensis 'Hongyang' cpDNA showed greater variation than the coding regions. The nucleotide polymorphism analysis (Pi) yielded nine highly variable region hotspots, most in the large single copy (LSC) region. The cp genome boundary analysis revealed a conservative order of gene arrangement in the inverted repeats (IRs) region of the cp genomes of 15 Actinidia plants, with small expansions and contractions of the boundaries. Furthermore, phylogenetic tree indicated that A. chinensis 'Hongyang' was the closest relative to A. indochinensis. This research provides a useful basis for future genetic and evolutionary studies of A. chinensis 'Hongyang', and enriches the biological information of Actinidia species.

Complete chloroplast genome sequence of Camellia sinensis: genome structure, adaptive evolution, and phylogenetic relationships

The chloroplast (cp) genome holds immense potential for a variety of applications including species identification, phylogenetic analysis, and evolutionary studies. In this study, we utilized Illumina NovaSeq 6000 to sequence the DNA of Camellia sinensis L. cultivar 'Zhuyeqi', followed by the assembly of its chloroplast genome using SPAdes v3.10.1, with subsequent analysis of its features and phylogenetic placement. The results showed that the cp genome of 'Zhuyeqi' was 157,072 bp, with a large single-copy region (LSC, 86,628 bp), a small single-copy region (SSC,18,282 bp), and two inverted repeat regions (IR, 26,081 bp). The total AT and GC contents of the cp genome of 'Zhuyeqi' were observed to be 62.21% and 37.29%, respectively. The cp genome encoded 135 unique genes, including 90 protein-coding genes (CDS), 37 tRNA genes, and 8 rRNA genes. Moreover, 31 codons and 247 simple sequence repeats (SSRs) were identified. The cp genomes of 'Zhuyeqi' were found to be relatively conserved, with conservation observed in the IR region, which showed no evidence of inversions or rearrangements. The five regions with the largest variations were identified, with four regions (rps12, rps19, rps16, and rpl33) located in the LSC region and one divergent region (trnI-GAU) in the IR region. Phylogenetic analysis revealed that Camellia sinensis (KJ996106.1) was closely related to 'Zhuyeqi', indicating a close phylogenetic relationship between these two species. These findings could provide important genetic information for further research into breeding of tea tree, phylogeny, and evolution of Camellia sinensis.

Developing population identification tool based on polymorphism of rDNA for traditional Chinese medicine: Artemisia annua L

BACKGROUND

Artemisia annua, a well-known traditional Chinese medicine, is the main source for production of artemisinin, an anti-malaria drug. A. annua is distributed globally, with great diversity of morphological characteristics and artemisinin contents. Diverse traits among A. annua populations impeded the stable production of artemisinin, which needs an efficient tool to identify strains and assess population genetic homogeneity.

PURPOSE

In this study, ribosomal DNA (rDNA), were characterized for A. annua for strains identification and population genetic homogeneity assessment.

METHODS

The ribosomal RNA (rRNA) genes were identified using cmscan and assembled using rDNA unit of LQ-9 as a reference. rDNA among Asteraceae species were compared performing with 45S rDNA. The rDNA copy number was calculated based on sequencing depth. The polymorphisms of rDNA sequences were identified with bam-readcount, and confirmed by Sanger sequencing and restriction enzyme experiment. The ITS2 amplicon sequencing was used to verify the stability of ITS2 haplotype analysis.

RESULTS

Different from other Asteraceae species, 45S and 5S linked-type rDNA was only found in Artemisia genus. Rich polymorphisms of copy number and sequence of rDNA were identified in A. annua population. The haplotype composition of internal transcribed spacer 2 (ITS2) region which had moderate sequence polymorphism and relative short size was significantly different among A. annua strains. A population discrimination method was developed based on ITS2 haplotype analysis with high-throughput sequencing.

CONCLUSION

This study provides comprehensive characteristics of rDNA and suggests that ITS2 haplotype analysis is ideal tool for A. annua strain identification and population genetic homogeneity assessment.

Comparative plastome genomics, taxonomic delimitation and evolutionary divergences of Tetraena hamiensis var. qatarensis and Tetraena simplex (Zygophyllaceae)

The Zygophyllum and Tetraena genera are intriguingly important ecologically and medicinally. Based on morphological characteristics, T. hamiensis var. qatarensis, and T. simplex were transferred from Zygophyllum to Tetraena with the least genomic datasets available. Hence, we sequenced the T. hamiensis and T. simplex and performed in-depth comparative genomics, phylogenetic analysis, and estimated time divergences. The complete plastomes ranged between 106,720 and 106,446 bp-typically smaller than angiosperms plastomes. The plastome circular genomes are divided into large single-copy regions (~ 80,964 bp), small single-copy regions (~ 17,416 bp), and two inverted repeats regions (~ 4170 bp) in both Tetraena species. An unusual shrinkage of IR regions 16-24 kb was identified. This resulted in the loss of 16 genes, including 11 ndh genes which encode the NADH dehydrogenase subunits, and a significant size reduction of Tetraena plastomes compared to other angiosperms. The inter-species variations and similarities were identified using genome-wide comparisons. Phylogenetic trees generated by analyzing the whole plastomes, protein-coding genes, matK, rbcL, and cssA genes exhibited identical topologies, indicating that both species are sisters to the genus Tetraena and may not belong to Zygophyllum. Similarly, based on the entire plastome and proteins coding genes datasets, the time divergence of Zygophyllum and Tetraena was 36.6 Ma and 34.4 Ma, respectively. Tetraena stem ages were 31.7 and 18.2 Ma based on full plastome and protein-coding genes. The current study presents the plastome as a distinguishing and identification feature among the closely related Tetraena and Zygophyllum species. It can be potentially used as a universal super-barcode for identifying plants.

Identification of three cultivated varieties of Scutellaria baicalensis using the complete chloroplast genome as a super-barcode

Scutellaria baicalensis has been one of the most commonly used traditional Chinese medicinal plants in China for more than 2000 years. The three new varieties cultivated could not be distinguished by morphology before flowering. It will hinder the promotion of later varieties. Chloroplast DNA has been widely used in species identification. Moreover, previous studies have shown that complete chloroplast genome sequences have been suggested as super barcodes for identifying plants. Therefore, we sequenced and annotated the complete chloroplast genomes of three cultivated varieties. The chloroplast genomes of SBW, SBR, and SBP were 151,702 bp, 151,799 bp, and 151,876 bp, which contained 85 protein-coding genes, 36 tRNA genes, and 8 rRNA genes. The analysis of the repeat sequences, codon usage, and comparison of chloroplast genomes shared a high degree of conservation. However, the sliding window results show significant differences among the three cultivated varieties in matK-rps16 and petA-psbJ. And we found that the matK-rps16 sequence can be used as a barcode for the identification of three varieties. In addition, the complete chloroplast genome contains more variations and can be used as a super-barcode to identify these three cultivated varieties. Based on the protein-coding genes, the phylogenetic tree demonstrated that SBP was more closely related to SBW, in the three cultivated varieties. Interestingly, we found that S. baicalensis and S. rehderiana are closely related, which provides new ideas for the development of S. baicalensis. The divergence time analysis showed that the three cultivated varieties diverged at about 0.10 Mya. Overall, this study showed that the complete chloroplast genome could be used as a super-barcode to identify three cultivated varieties of S. baicalensis and provide biological information, and it also contributes to bioprospecting.

Comparison of plastid genomes and ITS of two sister species in Gentiana and a discussion on potential threats for the endangered species from hybridization

BACKGROUND

Gentiana rigescens Franchet is an endangered medicinal herb from the family Gentianaceae with medicinal values. Gentiana cephalantha Franchet is a sister species to G. rigescens possessing similar morphology and wider distribution. To explore the phylogeny of the two species and reveal potential hybridization, we adopted next-generation sequencing technology to acquire their complete chloroplast genomes from sympatric and allopatric distributions, as along with Sanger sequencing to produce the nrDNA ITS sequences.

RESULTS

The plastid genomes were highly similar between G. rigescens and G. cephalantha. The lengths of the genomes ranged from 146,795 to 147,001 bp in G. rigescens and from 146,856 to 147,016 bp in G. cephalantha. All genomes consisted of 116 genes, including 78 protein-coding genes, 30 tRNA genes, four rRNA genes and four pseudogenes. The total length of the ITS sequence was 626 bp, including six informative sites. Heterozygotes occurred intensively in individuals from sympatric distribution. Phylogenetic analysis was performed based on chloroplast genomes, coding sequences (CDS), hypervariable sequences (HVR), and nrDNA ITS. Analysis based on all the datasets showed that G. rigescens and G. cephalantha formed a monophyly. The two species were well separated in phylogenetic trees using ITS, except for potential hybrids, but were mixed based on plastid genomes. This study supports that G. rigescens and G. cephalantha are closely related, but independent species. However, hybridization was confirmed to occur frequently between G. rigescens and G. cephalantha in sympatric distribution owing to the lack of stable reproductive barriers. Asymmetric introgression, along with hybridization and backcrossing, may probably lead to genetic swamping and even extinction of G. rigescens.

CONCLUSION

G. rigescens and G. cephalantha are recently diverged species which might not have undergone stable post-zygotic isolation. Though plastid genome shows obvious advantage in exploring phylogenetic relationships of some complicated genera, the intrinsic phylogeny was not revealed because of matrilineal inheritance here; nuclear genomes or regions are hence crucial for uncovering the truth. As an endangered species, G. rigescens faces serious threats from both natural hybridization and human activities; therefore, a balance between conservation and utilization of the species is extremely critical in formulating conservation strategies.

Morphological variability of wild-growing crown imperial (Fritillaria imperialis L.) germplasm in central region of Iran-implications for in-situ conservation initiatives

BACKGROUND

Crown imperial (Fritillaria imperialis L.) is a threatened bulbous plant which has great ornamental and medicinal values and importance. In the present study, a total of 100 specimens of wild-growing F. imperialis from 10 natural areas of Markazi province, Iran, representing one of the main centers of genetic diversity of this species, were evaluated using 37 phenotypic attributes during April 2021.

RESULTS

High level of genetic variation within populations (75%) and low levels of genetic variation among populations (25%) was revealed. The highest coefficient of variation (CV) was found in leaf trichome (82.00%) and then margin of crown leaves (80.44%). In addition, flower color (CV = 50.86%), flower number (CV = 44.61%), peduncle diameter (CV = 33.44%), and plant length (CV = 32.55%)-all important from an ornamental point of view- showed relatively high CV values. The CV was the lowest for flower shape, filament color, bulb shape, bulblet number, and floral scent. Ward cluster analysis identified two main clusters, containing 14 and 86 specimens, respectively. The first group consisted mainly of specimens from the adjacent Shahbaz and Rasvand populations. According to the principal component analysis (PCA), the first six components of data accounted for 88.36% of total variance. The Shahbaz-1, Shahbaz-2, Shahbaz-6, Shahbaz-7, Shahbaz-9, and Bolagh-8 specimens showed the highest variation and were separated from others, which they can be used further in breeding programs, while Sarchal-2, Bolagh-3, and Chepeqli-4 specimens showed the lowest variability. Moreover, the studied populations were clustered into four distinct groups, each including populations that were geographically close to one another.

CONCLUSIONS

Although the examined specimens revealed high genetic diversity herein, the results indicated that wild-growing populations of F. imperialis are still at risk suffering from overcollection in the most of studied areas, especially in Deh-Sad and Tureh.

Exploring the potential mechanism of Fritiliariae Irrhosae Bulbus on ischemic stroke based on network pharmacology and experimental validation

Objective: To study the potential targets and molecular mechanisms of Fritiliariae Irrhosae Bulbus (FIB) in the treatment of ischemic strokes based on a network pharmacology strategy, with a combination of molecular docking and animal experiments.

Methods: The active components and targets of FIB were screened by TCMSP database and TCMIP database, and the related targets of ischemic strokes were screened by GeneCards, OMIM, CTD, and DrugBank, then the intersection targets of the two were taken. The protein interaction network was constructed by STRING, the PPI network diagram was drawn by using Cytoscape software, and the key targets of FIB treatment of ischemic strokes were analyzed by MCODE. The DAVID database was used for GO and KEGG enrichment analysis, and the potential pathway of FIB against ischemic strokes was obtained. Molecular docking was performed by using AutoDock Tools 1.5.6 software. Finally, a mouse model of ischemic stroke was established, and the results of network pharmacology were verified by in vivo experiments. Realtime Polymerase Chain Reaction was used to detect the expression levels of relevant mRNAs in the mouse brain tissue. Western blot was used to detect the expression levels of related proteins in the mouse brain tissue.

Results: 13 kinds of active components of FIB were screened, 31 targets were found in the intersection of FIB and ischemic strokes, 10 key targets were obtained by MCODE analysis, 236 biological processes were involved in GO enrichment analysis, and key targets of KEGG enrichment analysis were mainly concentrated in Neuroactive light receptor interaction, Calcium signaling pathway, Cholinergic synapse, Hepatitis B, Apoptosis—multiple specifications, Pathways in cancer and other significantly related pathways. There was good binding activity between the screened main active components and target proteins when molecular docking was performed. Animal experiments showed that the infarct volume of brain tissue in the FIB treatment group was considerably reduced. RT-qPCR and the results of Western Blot showed that FIB could inhibit the expression of active-Caspase3, HSP90AA1, phosphorylated C-JUN, and COX2.

Conclusion: Based on network pharmacology, the effect of FIB in the treatment of ischemic strokes was discussed through the multi-component-multi-target-multi-pathway. The therapeutic effect and potential mechanisms of FIB on ischemic strokes were preliminarily explored, which provided a ground work for further researches on the pharmacodynamic material basis, mechanism of action and clinical application.

Application of Analytical Technologies in the Discrimination and Authentication of Herbs from Fritillaria: A Review

Medicinal plants of Fritillaria are widely distributed in numerous countries around the world and possess excellent antitussive and expectorant effects. In particular, Fritillariae Bulbus (FB) as a precious traditional medicine has thousands of years of medical history in China. Herbs of Fritillaria have a high market value and demand while limited by harsh growing circumstances and scarce wild resources. As a consequence, fraudulent behaviors are regularly engaged by the unscrupulous merchants in an attempt to reap greater profits. It is of an urgent need to evaluate the quality of Fritillaria herbs and their products using various analytical instruments and techniques. This review has scrutinized approximately 160 articles from 1995 to 2022 published on the investigation of Fritillaria herbs and related herbal products. The botanical classification of genus Fritillaria, types of counterfeits, technologies applied for differentiating Fritillaria species were comprehensively summarized and discussed in the current review. Molecular and chromatographic identification were the dominant technologies in the authentication of Fritillaria herbs. Additionally, we brought some potential and promising technologies and analytical strategies into attention, which are worthy attempting in the future researches. This review could conduce to excellent reference value for further investigations of the authenticity assessment of Fritillaria species.

Structural Characterization of the Acer ukurunduense Chloroplast Genome Relative to Related Species in the Acer Genus

Acer ukurunduense refers to a deciduous tree distributed in Northeast Asia and is a widely used landscaping tree species. Although several studies have been conducted on the species’ ecological and economic significance, limited information is available on its phylo-genomics. Our study newly constitutes the complete chloroplast genome of A. ukurunduense into a 156,645-bp circular DNA, which displayed a typical quadripartite structure. In addition, 133 genes were identified, containing 88 protein-coding genes, 37 tRNA genes, and eight rRNA genes. In total, 107 simple sequence repeats and 49 repetitive sequences were observed. Thirty-two codons indicated that biased usages were estimated across 20 protein-coding genes (CDS) in A. ukurunduense. Four hotspot regions (trnK-UUU/rps16, ndhF/rpl32, rpl32/trnL-UAG, and ycf1) were detected among the five analyzed Acer species. Those hotspot regions may be useful molecular markers and contribute to future population genetics studies. The phylogenetic analysis demonstrated that A. ukurunduense is most closely associated with the species of Sect. Palmata. A. ukurunduense and A. pubipetiolatum var. pingpienense diverged in 22.11 Mya. We selected one of the hypervariable regions (trnK-UUU/rps16) to develop a new molecular marker and designed primers and confirmed that the molecular markers could accurately discriminate five Acer species through Sanger sequencing. By sequencing the cp genome of A. ukurunduense and comparing it with the relative species of Acer, we can effectively address the phylogenetic problems of Acer at the species level and provide insights into future research on population genetics and genetic diversity.

Complete chloroplast genome of Lilium ledebourii (Baker) Boiss and its comparative analysis: lights into selective pressure and adaptive evolution

Lilium ledebourii (Baker) Boiss is a rare species, which exhibits valuable traits. However, before its genetic diversity and evolutionary were uncovered, its wild resources were jeopardized. Moreover, some ambiguities in phylogenetic relationships of this genus remain unresolved. Therefore, obtaining the whole chloroplast sequences of L. ledebourii and its comparative analysis along with other Lilium species is crucial and pivotal to understanding the evolution of this genus as well as the genetic populations. A multi-scale genome-level analysis, especially selection pressure, was conducted. Detailed third‑generation sequencing and analysis revealed a whole chloroplast genome of 151,884 bp, with an ordinary quadripartite and protected structure comprising 37.0% GC. Overall, 113 different genes were recognized in the chloroplast genome, consisting of 30 distinct tRNA genes, four distinct ribosomal RNAs genes, and 79 unique protein-encoding genes. Here, 3234 SSRs and 2053 complex repeats were identified, and a comprehensive analysis was performed for IR expansion and contraction, and codon usage bias. Moreover, genome-wide sliding window analysis revealed the variability of rpl32-trnL-ccsA, petD-rpoA, ycf1, psbI-trnS-trnG, rps15-ycf1, trnR, trnT-trnL, and trnP-psaJ-rpl33 were higher among the 48 Lilium cp genomes, displaying higher variability of nucleotide in SC regions. Following 1128 pairwise comparisons, ndhB, psbJ, psbZ, and ycf2 exhibit zero synonymous substitution, revealing divergence or genetic restriction. Furthermore, out of 78 protein-coding genes, we found that accD and rpl36 under positive selection: however, at the entire-chloroplast protein scale, the Lilium species have gone through a purifying selection. Also, a new phylogenetic tree for Lilium was rebuilt, and we believe that the Lilium classification is clearer than before. The genetic resources provided here will aid future studies in species identification, population genetics, and Lilium conservation.

A Comparative Analysis of the Chloroplast Genomes of Four Polygonum Medicinal Plants

Polygonum is a generalized genus of the Polygonaceae family that includes various herbaceous plants. In order to provide aid in understanding the evolutionary and phylogenetic relationship in Polygonum at the chloroplast (cp) genome-scale level, we sequenced and annotated the complete chloroplast genomes of four Polygonum species using next-generation sequencing technology and CpGAVAS. Then, repeat sequences, IR contractions, and expansion and transformation sites of chloroplast genomes of four Polygonum species were studied, and a phylogenetic tree was built using the chloroplast genomes of Polygonum. The results indicated that the chloroplast genome construction of Polygonum also displayed characteristic four types of results, comparable to the published chloroplast genome of recorded angiosperms. The chloroplast genomes of the four Polygonum plants are highly consistent in genome size (159,015 bp-163,461 bp), number of genes (112 genes, including 78 protein-coding genes, 30 tRNA genes, and 4 rRNA genes), gene types, gene order, codon usage, and repeat sequence distribution, which identifies the high preservation among the Polygonum chloroplast genomes. The Polygonum phylogenetic tree was recreated by a full sequence of the chloroplast genome, which illustrates that the P. bistorta, P. orientale, and P. perfoliatum are divided into the same branch, and P. aviculare belongs to Fallopia. The precise system site of lots base parts requires further verification, but the study would provide a basis for developing the available genetic resources and evolutionary relationships of Polygonum.

The first complete chloroplast genome of Vicatia thibetica de Boiss.: genome features, comparative analysis, and phylogenetic relationships

Vicatia thibetica de Boiss.: a herb in the family Apiaceae, has been used for over a hundred years as an essential medicinal and edible plant in the Bai ethnic group of Dali City. However, due to the lack of study on plastid genomes of V. thibetica, studies of comparison and phylogeny with other related species remain scarce. In the current study, we assembled, annotated, and characterized the entire chloroplast (cp) genome of V. thibetica through high-throughput sequencing for the first time, compared with published whole chloroplast genomes from the same family. A phylogenetic analysis of the chloroplast genome has also been performed. The whole chloroplast genome of V. thibetica was 145,796 in size and consisted of a large single-copy region (LSC; 92,186 bp), a small single-copy region (SSC; 17,452 bp), and a pair of inverted repeat regions (IRs; 18,079 bp) forming a circular quadripartite structure. Annotation resulted in 128 genes, including 84 protein-coding genes (PCGs), 35 transfer RNA genes (tRNAs), eight ribosomal genes (rRNAs), and one pseudogene. Repeat sequence analysis displayed V. thibetica plastid genome contains 75 simple repeats, 37 long repeats, and 29 tandem repeats. Compared with the cp genome of other Apiaceae species, a common feature was that the IR regions of the genome were more conservative compared to the LSC and SSC regions. Highly variable hotspots included rps16, ndhC-trnV-UAC, clpP, ycf1, and ndhB in the genomes, which supply valuable molecular markers for phylogeny, identification, and classification in the Apiaceae family. The results of phylogenetic analysis strongly supported the genus Vicatia as an independent genus in the family Apiaceae, in which the closest affinities to the related species of Angelica, Peucedanum, and Ligusticum were observed. In conclusion, the first chloroplast genome of Vicatia reported in this study may  improve our understanding of phylogenetic relationship of different genera of Apiaceae. In addition, the current data will be valuable as chloroplast genomic resource for species identification and population genetics.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01154-y.

Bulbus Fritillariae Cirrhosae as a Respiratory Medicine: Is There a Potential Drug in the Treatment of COVID-19?

Bulbus fritillariae cirrhosae (BFC) is one of the most used Chinese medicines for lung disease, and exerts antitussive, expectorant, anti-inflammatory, anti-asthmatic, and antioxidant effects, which is an ideal therapeutic drug for respiratory diseases such as ARDS, COPD, asthma, lung cancer, and pulmonary tuberculosis. Through this review, it is found that the therapeutic mechanism of BFC on respiratory diseases exhibits the characteristics of multi-components, multi-targets, and multi-signaling pathways. In particular, the therapeutic potential of BFC in terms of intervention of “cytokine storm”, STAT, NF-κB, and MAPK signaling pathways, as well as the renin-angiotensin system (RAS) that ACE is involved in. In the “cytokine storm” of SARS-CoV-2 infection there is an intense inflammatory response. ACE2 regulates the RAS by degradation of Ang II produced by ACE, which is associated with SARS-CoV-2. For COVID-19, may it be a potential drug? This review summarized the research progress of BFC in the respiratory diseases, discussed the development potentiality of BFC for the treatment of COVID-19, explained the chemical diversity and biological significance of the alkaloids in BFC, and clarified the material basis, molecular targets, and signaling pathways of BFC for the respiratory diseases. We hope this review can provide insights on the drug discovery of anti-COVID-19.

Identification of evolutionary relationships and DNA markers in the medicinally important genus Fritillaria based on chloroplast genomics

The genus Fritillaria has attracted great attention because of its medicinal and ornamental values. At least three reasons, including the accurate discrimination between various Fritillaria species, protection and sustainable development of rare Fritillaria resources as well as understanding of relationship of some perplexing species, have prompted phylogenetic analyses and development of molecular markers for Fritillaria species. Here we determined the complete chloroplast (CP) genomes for F. unibracteata, F. przewalskii, F. delavayi, and F. sinica through Illumina sequencing, followed by de novo assembly. The lengths of the genomes ranged from 151,076 in F. unibracteata to 152,043 in F. przewalskii. Those CP genomes displayed a typical quadripartite structure, all including a pair of inverted repeats (26,078 to 26,355 bp) separated by the large single-copy (81,383 to 81,804 bp) and small single-copy (17,537 to 17,569 bp) regions. Fritillaria przewalskii, F. delavayi, and F. sinica equivalently encoded 133 unique genes consisting of 38 transfer RNA genes, eight ribosomal RNA genes, and 87 protein coding genes, whereas F. unibracteata contained 132 unique genes due to absence of the rps16 gene. Subsequently, comparative analysis of the complete CP genomes revealed that ycf1, trnL, trnF, ndhD, trnN-trnR, trnE-trnT, trnN, psbM-trnD, atpI, and rps19 to be useful molecular markers in taxonomic studies owning to their interspecies variations. Based on the comprehensive CP genome data collected from 53 species in Fritillaria and Lilium genera, a phylogenomic study was carried out with three Cardiocrinum species and five Amana species as outgroups. The results of the phylogenetic analysis showed that Fritillaria was a sister to Lilium, and the interspecies relationships within subgenus Fritillaria were well resolved. Furthermore, phylogenetic analysis based on the CP genome was proved to be a promising method in selecting potential novel medicinal resources to substitute current medicinal species that are on the verge of extinction.

Comparative Analysis of Chloroplast Genome and New Insights Into Phylogenetic Relationships of Polygonatum and Tribe Polygonateae

Members of Polygonatum are perennial herbs that have been widely used in traditional Chinese medicine to invigorate Qi, moisten the lung, and benefit the kidney and spleen among patients. However, the phylogenetic relationships and intrageneric taxonomy within Polygonatum have long been controversial because of the complexity of their morphological variations and lack of high-resolution molecular markers. The chloroplast (cp) genome is an optimal model for deciphering phylogenetic relationships in related families. In the present study, the complete cp genome of 26 species of Trib. Polygonateae were de novo assembled and characterized; all species exhibited a conserved quadripartite structure, that is, two inverted repeats (IR) containing most of the ribosomal RNA genes, and two unique regions, large single sequence (LSC) and small single sequence (SSC). A total of 8 highly variable regions (rps16-trnQ-UUG, trnS-GCU-trnG-UCC, rpl32-trnL-UAG, matK-rps16, petA-psbJ, trnT-UGU-trnL-UAA, accD-psaI, and trnC-GCA-petN) that might be useful as potential molecular markers for identifying Polygonatum species were identified. The molecular clock analysis results showed that the divergence time of Polygonatum might occur at ∼14.71 Ma, and the verticillate leaf might be the ancestral state of this genus. Moreover, phylogenetic analysis based on 88 cp genomes strongly supported the monophyly of Polygonatum. The phylogenetic analysis also suggested that Heteropolygonatum may be the sister group of the Polygonatum, but the Disporopsis, Maianthemum, and Disporum may have diverged earlier. This study provides valuable information for further species identification, evolution, and phylogenetic research of Polygonatum.

DNA Barcoding and Phylogenomic Analysis of the Genus Fritillaria in China Based on Complete Chloroplast Genomes

The Fritillaria is an extremely complicated genus in taxonomy and phylogeny, which contains numerous medicinal species in China. Both traditional characteristic-based taxonomy and universal DNA barcodes (ITS, trnH-psbA, and rbcL) are difficult to effectively identify the species. Here, we generated a large dataset of chloroplast genomes from multiple accessions per species of Fritillaria to evaluate their effectiveness in species discrimination. Moreover, phylogeny of species in China was explored based on the complete chloroplast genomes, and then divergence times of each node were estimated. The results showed that all 21 species in Fritillaria here (including two suspicious species) could be correctly discriminated using cpDNA genomes except F. cirrhosa, which suggested that DNA super-barcode could greatly enhance species discriminatory resolution for complicated genera. Furthermore, four regions (ycf1, matK-trnG-GCC, rpoC1, and matK) gained remarkably higher resolution than that of other plastid regions, but only matK might be suitable to identify Fritillaria species in consideration of its lengths. Phylogenomic analysis showed that the subgenus Fritillaria in China was divided into four major clades with obvious geographic structure. Among them, Clade I, mainly distributed in southwest China, was a young and complicated group. Moreover, according to the analysis, taxonomic treatments of the two suspicious species, namely "F. omeiensis" and "F. hupehensis" in Flora of China (2000) are questionable and might need further revision. Molecular dating revealed that both origin and divergence of subgenus Fritillaria, as well as its four major clades, were significantly associated with geological and climatic fluctuations during the Middle to Late Miocene. This study would enrich case studies of DNA super-barcode and provide new insights on speciation, lineage diversification, and biogeography of the Fritillaria in China.

Complete Chloroplast Genome Sequence of Sonchus brachyotus Helps to Elucidate Evolutionary Relationships with Related Species of Asteraceae

Sonchus brachyotus DC. possesses both edible and medicinal properties and is widely distributed throughout China. In this study, the complete cp genome of S. brachyotus was sequenced and assembled. The total length of the complete S. brachyotus cp genome was 151,977 bp, including an LSC region of 84,553 bp, SSC region of 18,138 bp, and IR region of 24,643 bp. Sequence analyses revealed that the cp genome encoded 132 genes, including 87 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. The GC content was 37.6%. One hundred mononucleotide microsatellites, 4 dinucleotide microsatellites, 67 trinucleotide microsatellites, 4 tetranucleotide microsatellites, and 1 long repeat were identified. The SSR frequency of the LSC region was significantly greater than that of the IR and SSC regions. In total, 175 SSRs and highly variable regions were recognized as potential cp markers. By analyzing the IR/LSC and IR/SSC boundaries, structural differences between S. brachyotus and 6 other species were detected. According to phylogenetic analyses, S. brachyotus was most closely related to S. arvensis and S. oleraceus. Overall, this study provides complete cp genome resources for S. brachyotus that will be beneficial for identifying potential molecular markers and evolutionary patterns of S. brachyotus and its closely related species.

Chloroplast genome sequencing based on genome skimming for identification of Eriobotryae Folium

Background

Whole chloroplast genome (cpDNA) sequence is becoming widely used in the phylogenetic studies of plant and species identification, but in most cases the cpDNA were acquired from silica gel dried fresh leaves. So far few reports have been available to describe cpDNA acquisition from crude drugs derived from plant materials, the DNA of which usually was seriously damaged during their processing. In this study, we retrieved cpDNA from the commonly used crude drug Eriobotryae Folium (Pipaye in Chinese, which is the dried leaves of Eriobotrya japonica, PPY) using genome skimming technique.

Results

We successfully recovered cpDNA sequences and rDNA sequences from the crude drug PPY, and bioinformatics analysis showed a high overall consistency between the cpDNA obtained from the crude drugs and fresh samples. In the ML tree, each species formed distinct monophyletic clades based on cpDNA sequence data, while the phylogenetic relationships between Eriobotrya species were poorly resolved based on ITS and ITS2.

Conclusion

Our results demonstrate that both cpDNA and ITS/ITS2 are effective for identifying PPY and its counterfeits derived from distantly related species (i.e. Dillenia turbinata and Magnolia grandiflora), but cpDNA is more effective for distinguishing the counterfeits derived from the close relatives of Eriobotrya japonica, suggesting the potential of genome skimming for retrieving cpDNA from crude drugs used in Traditional Chinese Medicine for their identification.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12896-021-00728-0.

Plant super-barcode: a case study on genome-based identification for closely related species of Fritillaria

Background

Although molecular analysis offers a wide range of options for species identification, a universal methodology for classifying and distinguishing closely related species remains elusive. This study validated the effectiveness of utilizing the entire chloroplast (cp) genome as a super-barcode to help identify and classify closely related species.

Methods

We here compared 26 complete cp genomes of ten Fritillaria species including 18 new sequences sequenced in this study. Each species had repeats and the cp genomes were used as a whole DNA barcode to test whether they can distinguish Fritillaria species.

Results

The cp genomes of Fritillaria medicinal plants were conserved in genome structure, gene type, and gene content. Comparison analysis of the Fritillaria cp genomes revealed that the intergenic spacer regions were highly divergent compared with other regions. By constructing the phylogenetic tree by the maximum likelihood and maximum parsimony methods, we found that the entire cp genome showed a high discrimination power for Fritillaria species with individuals of each species in a monophyletic clade. These results indicate that cp genome can be used to effectively differentiate medicinal plants from the genus Fritillaria at the species level.

Conclusions

This study implies that cp genome can provide distinguishing differences to help identify closely related Fritillaria species, and has the potential to be served as a universal super-barcode for plant identification.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13020-021-00460-z.

Complete chloroplast genome of Salvia plebeia: organization, specific barcode and phylogenetic analysis

Salvia plebeia has been in use as traditional Chinese medicine (TCM) for more than 500 years. In this study, the complete chloroplast (cp) genome of S. plebeia was sequenced, assembled and compared to those of other five published Salvia cp genomes. It was found that the cp genome structure of S. plebeia was well conserved and had a total size of 151 062 bp. Four parameters were used to display the usage conditions of the codons of the amino acids in Salvia genus. Although the number of protein-coding genes in each species was the same, the total number of codons was different. Except for amino acids Trp and Met whose Relative Synonymous Codon Usage (RSCU) value of one condon was equal to 1, the remaining 19 amino acids had 1-3 preferred codons. The preferred codon names of each amino acid were coincident. The period size for the tandem repeats of six species ranged from 9 to 410 bp. Salvia cp genomes mainly possessed tandem repeats with a copy number less than or equal to 3. The sequence length of tandem repeats of the six species ranged from 25 to 824 bp. Highly viarable regions including four intergenic spacers and six partial genes were discovered as potential specific barcodes for Salvia species through cp genome-wide comparison. Finally, we performed phylogenetic analyses based on the complete cp genome and coding sequences respectively. These results provide information to help construct the cp genome library for Salvia, which may support studies of phylogenetics, DNA barcoding, population and transplastomics.

The Complete Chloroplast Genome of the Vulnerable Oreocharis esquirolii (Gesneriaceae): Structural Features, Comparative and Phylogenetic Analysis

Oreocharis esquirolii, a member of Gesneriaceae, is known as Thamnocharis esquirolii, which has been regarded a synonym of the former. The species is endemic to Guizhou, southwestern China, and is evaluated as vulnerable (VU) under the International Union for Conservation of Nature (IUCN) criteria. Until now, the sequence and genome information of O. esquirolii remains unknown. In this study, we assembled and characterized the complete chloroplast (cp) genome of O. esquirolii using Illumina sequencing data for the first time. The total length of the cp genome was 154,069 bp with a typical quadripartite structure consisting of a pair of inverted repeats (IRs) of 25,392 bp separated by a large single copy region (LSC) of 85,156 bp and a small single copy region (SSC) of18,129 bp. The genome comprised 114 unique genes with 80 protein-coding genes, 30 tRNA genes, and four rRNA genes. Thirty-one repeat sequences and 74 simple sequence repeats (SSRs) were identified. Genome alignment across five plastid genomes of Gesneriaceae indicated a high sequence similarity. Four highly variable sites (rps16-trnQ, trnS-trnG, ndhF-rpl32, and ycf 1) were identified. Phylogenetic analysis indicated that O. esquirolii grouped together with O. mileensis, supporting resurrection of the name Oreocharis esquirolii from Thamnocharisesquirolii. The complete cp genome sequence will contribute to further studies in molecular identification, genetic diversity, and phylogeny.

Comparative Genomics and Phylogenetic Analyses of Christia vespertilionis and Urariopsis brevissima in the Tribe Desmodieae (Fabaceae: Papilionoideae) Based on Complete Chloroplast Genomes

Taxonomic and phylogenetic relationships of Christia, Urariopsis, Uraria and related genera within the tribe Desmodieae (Fabaceae: Papilionoideae) have long been controversial. Here, we report the complete chloroplast (cp) genomes of Christia vespertilionis and Urariopsis brevissima and perform comparative and phylogenetic analyses with Uraria lagopodioides and other relatives in the Desmodieae. The cp genomes of C. vespertilionis and U. brevissima are 149,656 and 149,930 bp long, with 128 unique genes (83 protein-coding genes, 37 tRNA genes and 8 rRNA genes), respectively. Comparative analyses revealed 95-129 simple sequence repeats (SSRs) and eleven highly variable regions (trnK-rbcL, rbcL-atpB, ndhJ-trnF, trnL-trnT, psbD-rpoB, accD-cemA, petA-psbL, psbE-petL, rps11-rps19, ndhF-ccsA, and rps15-ycf1) among six Desmodieae species. Phylogenetic analyses clearly resolved two subtribes (Desmodiinae and Lespedezinae) of Desmodieae as monophyletic, and the newly reported C. vespertilionis and U. brevissima clustered in subtribe Desmodiinae. A sister relationship of C. vespertilionis to U. lagopodioides was supported. Evidence was presented to support the treatment of Urariopsis as a distinct genus rather than in synonymy with Uraria. The results provide valuable information for further studies on species delimitation, phylogenetics, population genetics, and the evolutionary process of speciation in the Desmodieae.

Adaptation Evolution and Phylogenetic Analyses of Species in Chinese Allium Section Pallasia and Related Species Based on Complete Chloroplast Genome Sequences

The section Pallasia is one of the components of the genus Allium subgenus Allium (Amaryllidaceae), and species relationship in this section is still not resolved very well, which hinders further evolutionary and adaptive studies. Here, the complete chloroplast genomes of five sect. Pallasia species were reported, and a comparative analysis was performed with other three related Allium species. The genome size of the eight species ranged from 151,672 bp to 153,339 bp in length, GC content changed from 36.7% to 36.8%, and 130 genes (except Allium pallasii), 37 tRNA, and 8 rRNA were identified in each genome. By analyzing the IR/LSC and IR/SSC boundary, A. pallasii exhibited differences compared with other seven species. Phylogenetic analysis achieved high supports in each branch, seven of the eight Allium species cluster into a group, and A. pallasii exhibit a close relationship with A. obliquum. Higher pairwise Ka/Ks ratios were found in A. schoenoprasoides compared to A. caeruleum and A. macrostemon while a lower value of Ka/Ks ratios was detected between A. caeruleum and A. macrostemon. This study will be a great contribution to the future phylogenetic and adaptive research in Allium.