Plant DNA barcoding: from gene to genome

DNA metabarcoding diet analysis in a generalist omnivore: feeding trials reveal the efficacy of extraction kits and a multi-locus approach for identifying diverse diets

Metabarcoding-based diet analysis is a valuable tool for understanding the feeding behavior of a wide range of species. However, many studies using these methods for wild animals assume accuracy and precision without experimental evaluation with known positive control food items. Here, we conducted a feeding trial experiment with a positive control community in pasture-raised chickens and assessed the efficacy of several commonly used DNA extraction kits and primer sets. We hand-fed 22 known food items, including insects and plants, to six backyard laying hens and collected their excreta for eight h. We evaluated the efficacy of three DNA extraction kits, three primer sets for plant identification (targeting rbcL, trnL, and internal transcribed spacer 2 [ITS2]), and three primer sets for arthropod identification (targeting cytochrome oxidase subunit I [COI]). The detection success rate of our positive control food items was highly variable, ranging from 2.04% to 93.88% for all kit/primer combinations and averaging 37.35% and 43.57% for the most effective kit/primer combination for plants and insects, respectively. Extraction kits using bead-based homogenization positively affected the recovery proportion of plant and insect DNA in excreta samples. The minimum time to detect known food items was 44 min post-feeding. Two COI primer sets significantly outperformed the third, and both recovery proportion and taxonomic resolution from ITS2 were significantly higher than those from rbcL and trnL. Taken together, these results display the potential variability that can be inherently present in DNA-based diet analyses and highlight the utility of experimental feeding trials in validating such approaches, particularly for omnivores with diverse diets.

Phylogeny, adaptive evolution, and taxonomy of Acronema (Apiaceae): evidence from plastid phylogenomics and morphological data

Introduction

The genus Acronema, belonging to Apiaceae, includes approximately 25 species distributed in the high-altitude Sino-Himalayan region from E Nepal to SW China. This genus is a taxonomically complex genus with often indistinct species boundaries and problematic generic delimitation with Sinocarum and other close genera, largely due to the varied morphological characteristics.

Methods

To explore the phylogenetic relationships and clarify the limits of the genus Acronema and its related genera, we reconstructed a reliable phylogenetic framework with high support and resolution based on two molecular datasets (plastome data and ITS sequences) and performed morphological analyses.

Results

Both phylogenetic analyses robustly supported that Acronema was a non-monophyletic group that fell into two clades: Acronema Clade and East-Asia Clade. We also newly sequenced and assembled sixteen Acronema complete plastomes and performed comprehensively comparative analyses for this genus. The comparative results showed that the plastome structure, gene number, GC content, codon bias patterns were high similarity, but varied in borders of SC/IR and we identified six different types of SC/IR border. The SC/IR boundaries of Acronema chienii were significantly different from the other Acronema members which was consistent with the type VI pattern in the genus Tongoloa. We also identified twelve potential DNA barcode regions (ccsA, matK, ndhF, ndhG, psaI, psbI, rpl32, rps15, ycf1, ycf3, psaI-ycf4 and psbM-trnD) for species identification in Acronema. The molecular evolution of Acronema was relatively conservative that only one gene (petG) was found to be under positive selection (ω = 1.02489).

Discussion

The gene petG is one of the genes involved in the transmission of photosynthetic electron chains during photosynthesis, which plays a crucial role in the process of photosynthesis in plants. This is also a manifestation of the adaptive evolution of plants in high-altitude areas to the environment. In conclusion, our study provides novel insights into the plastome adaptive evolution, phylogeny, and taxonomy of genus Acronema.

The Complete Chloroplast Genome Sequence of the Medicinal Moss Rhodobryum giganteum (Bryaceae, Bryophyta): Comparative Genomics and Phylogenetic Analyses

Rhodobryum giganteum (Bryaceae, Bryophyta), a rare medicinal bryophyte, is valued for its cardiovascular therapeutic properties in traditional Chinese medicine. This study presents the first complete chloroplast genome sequence of R. giganteum, including its assembly and annotation. The circular chloroplast genome of R. giganteum is 124,315 bp in length, displaying a typical quadripartite structure with 128 genes: 83 protein-coding genes, 37 tRNAs, and 8 rRNAs. Analyses of codon usage bias, repetitive sequences, and simple sequence repeats (SSRs) revealed an A/U-ending codon preference, 96 repetitive sequences, and 385 SSRs in the R. giganteum chloroplast genome. Nucleotide diversity analysis identified 10 high mutational hotspots. Ka/Ks ratio analysis suggested potential positive selection in rpl20, rps18, petG, and psbM genes. Phylogenetic analysis of whole chloroplast genomes from 38 moss species positioned R. giganteum within Bryales, closely related to Rhodobryum laxelimbatum. This study augments the chloroplast genomic data for Bryales and provides a foundation for molecular marker development and genetic diversity analyses in medicinal bryophytes.

Repetitive Sequences, Codon Usage Bias and Phylogenetic Analysis of the Plastome of Miliusa glochidioides

Annonaceae is the largest family in Magnoliales, exhibiting the greatest diversity among and within genera. In this study, we conducted an analysis of repetitive sequences and codon usage bias in the previously acquired plastome of Miliusa glochidioides. Using a concatenated dataset of shared genes, we constructed the phylogenetic relationships among 27 Annonaceae species. The results showed that the size of the plastomes in the Annonaceae ranged from 159 to 202 kb, with the size of the inverted repeat region ranging from 40 to 65 kb. Within the plastome of M. glochidioides, we identified 42 SSRs, 36 tandem repeats, and 9 dispersed repeats. These SSRs consist of three nucleotide types and eight motif types, with a preference for A/T bases, primarily located in the large single-copy regions and intergenic spacers. Tandem and dispersed repeat sequences were predominantly detected in the IR region. Through codon usage bias analysis, we identified 30 high-frequency codons and 11 optimal codons. The plastome of M. glochidioides demonstrated relatively weak codon usage bias, favoring codons with A/T endings, primarily influenced by natural selection. Phylogenetic analysis revealed that all four subfamilies formed monophyletic groups, with Cananga odorata (Ambavioideae) and Anaxagorea javanica (Anaxagoreoideae) successively nested outside Annonoideae + Malmeoideae. These findings improve our understanding of the plastome of M. glochidioides and provide additional insights for studying plastome evolution in Annonaceae.

Advances in Physalis molecular research: applications in authentication, genetic diversity, phylogenetics, functional genes, and omics

The plants of the genus Physalis L. have been extensively utilized in traditional and indigenous Chinese medicinal practices for treating a variety of ailments, including dermatitis, malaria, asthma, hepatitis, and liver disorders. The present review aims to achieve a comprehensive and up-to-date investigation of the genus Physalis, a new model crop, to understand plant diversity and fruit development. Several chloroplast DNA-, nuclear ribosomal DNA-, and genomic DNA-based markers, such as psbA-trnH, internal-transcribed spacer (ITS), simple sequence repeat (SSR), random amplified microsatellites (RAMS), sequence-characterized amplified region (SCAR), and single nucleotide polymorphism (SNP), were developed for molecular identification, genetic diversity, and phylogenetic studies of Physalis species. A large number of functional genes involved in inflated calyx syndrome development (AP2-L, MPF2, MPF3, and MAGO), organ growth (AG1, AG2, POS1, and CNR1), and active ingredient metabolism (24ISO, DHCRT, P450-CPL, SR, DUF538, TAS14, and 3β-HSB) were identified contributing to the breeding of novel Physalis varieties. Various omic studies revealed and functionally identified a series of reproductive organ development-related factors, environmental stress-responsive genes, and active component biosynthesis-related enzymes. The chromosome-level genomes of Physalis floridana Rydb., Physalis grisea (Waterf.) M. Martínez, and Physalis pruinosa L. have been recently published providing a valuable resource for genome editing in Physalis crops. Our review summarizes the recent progress in genetic diversity, molecular identification, phylogenetics, functional genes, and the application of omics in the genus Physalis and accelerates efficient utilization of this traditional herb.

Characterization of the complete chloroplast genome of the rare medicinal plant: Mandragora caulescens (Solanaceae)

In this study, we assembled high-quality chloroplast genomes of Mandragora caulescens through a reference-guided approach using high-throughput Illumina sequencing reads. The resulting chloroplast genome assembly displayed a typical quadripartite structural organization, comprising a large single-copy (LSC) region of 85,233 bp, two inverted repeat (IR) regions of 25,685 bp each, and a small single-copy (SSC) region of 18,207 bp. The chloroplast genome harbored 141 complete genes, and its overall GC content was 38.0%. In maximum-likelihood (ML) and Bayesian inference (BI) trees, the 19 Solanaceae species formed a monophyletic group, dividing into two main clades. M. caulescens and Nicandra physalodes formed a monophyletic group, suggesting a close relationship between the two species. The M. caulescens cp genome presented in this study lays a good foundation for further genetic and genomic studies of the Solanaceae.

The chloroplast genome sequences of Ipomoea alba and I. obscura (Convolvulaceae): genome comparison and phylogenetic analysis

Ipomoea species have diverse uses as ornamentals, food, and medicine. However, their genomic information is limited; I. alba and I. obscura were sequenced and assembled. Their chloroplast genomes were 161,353 bp and 159,691 bp, respectively. Both genomes exhibited a quadripartite structure, consisting of a pair of inverted repeat (IR) regions, which are separated by the large single-copy (LSC) and small single-copy (SSC) regions. The overall GC content was 37.5% for both genomes. A total of 104 and 93 simple sequence repeats, 50 large repeats, and 30 and 22 short tandem repeats were identified in the two chloroplast genomes, respectively. G and T were more preferred than C and A at the third base position based on the Parity Rule 2 plot analysis, and the neutrality plot revealed correlation coefficients of 0.126 and 0.105, indicating the influence of natural selection in shaping the codon usage bias in most protein-coding genes (CDS). Genome comparative analyses using 31 selected Ipomoea taxa from Thailand showed that their chloroplast genomes are rather conserved, but the presence of expansion or contraction of the IR region was identified in some of these Ipomoea taxa. A total of five highly divergent regions were identified, including the CDS genes accD, ndhA, and ndhF, as well as the intergenic spacer regions psbI-atpA and rpl32-ccsA. Phylogenetic analysis based on both the complete chloroplast genome sequence and CDS datasets of 31 Ipomoea taxa showed that I. alba is resolved as a group member for series (ser.) Quamoclit, which contains seven other taxa, including I. hederacea, I. imperati, I. indica, I. nil, I. purpurea, I. quamoclit, and I. × sloteri, while I. obscura is grouped with I. tiliifolia, both of which are under ser. Obscura, and is closely related to I. biflora of ser. Pes-tigridis. Divergence time estimation using the complete chloroplast genome sequence dataset indicated that the mean age of the divergence for Ipomoeeae, Argyreiinae, and Astripomoeinae, was approximately 29.99 Mya, 19.81 Mya, and 13.40 Mya, respectively. The node indicating the divergence of I. alba from the other members of Ipomoea was around 10.06 Mya, and the split between I. obscura and I. tiliifolia is thought to have happened around 17.13 Mya. The split between the I. obscura accessions from Thailand and Taiwan is thought to have taken place around 0.86 Mya.

Structure, gene composition, divergence time and phylogeny analysis of the woody desert species Neltuma alba, Neltuma chilensis and Strombocarpa strombulifera

Neltuma alba (Algarrobo blanco), Neltuma chilensis (Algarrobo Chileno) and Strombocarpa strombulifera (Fortuna) are some of the few drought resistant trees and shrubs found in small highly fragmented populations, throughout the Atacama Desert. We reconstructed their plastid genomes using de novo assembly of paired-end reads from total genomic DNA. We found that the complete plastid genomes of N. alba and N. chilensis are larger in size compared to species of the Strombocarpa genus. The Strombocarpa species presented slightly more GC content than the Neltuma species. Therefore, we assume that Strombocarpa species have been exposed to stronger natural selection than Neltuma species. We observed high variation values in the number of cpSSRs (chloroplast simple sequence repeats) and repeated elements among Neltuma and Strombocarpa species. The p-distance results showed a low evolutionary divergence within the genus Neltuma, whereas a high evolutionary divergence was observed between Strombocarpa species. The molecular divergence time found in Neltuma and Strombocarpa show that these genera diverged in the late Oligocene. With this study we provide valuable information about tree species that provide important ecosystem services in hostile environments which can be used to determine these species in the geographically isolated communities, and keep the highly fragmented populations genetically healthy.

Comparative and phylogenetic analysis of the complete chloroplast genomes of 10 Artemisia selengensis resources based on high-throughput sequencing

BACKGROUND

Artemisia selengensis, classified within the genus Artemisia of the Asteraceae family, is a perennial herb recognized for its dual utility in culinary and medicinal domains. There are few studies on the chloroplast genome of A. selengensis, and the phylogeographic classification is vague, which makes phylogenetic analysis and evolutionary studies very difficult.

RESULTS

The chloroplast genomes of 10 A. selengensis in this study were highly conserved in terms of gene content, gene order, and gene intron number. The genome lengths ranged from 151,148 to 151,257 bp and were typical of a quadripartite structure with a total GC content of approximately 37.5%. The chloroplast genomes of all species encode 133 genes, including 88 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Due to the contraction and expansion of the inverted repeats (IR), the overlap of ycf1 and ndhF genes occurred at the inverted repeats B (IRB) and short single copy sequence (SSC) boundaries. According to a codon use study, the frequent base in the chloroplast genome of A. selengensis' third codon position was A/T. The number of SSR repeats was 42-44, most of which were single nucleotide A/T repeats. Sequence alignment analysis of the chloroplast genome showed that variable regions were mainly distributed in single copy regions, nucleotide diversity values of 0 to 0.009 were calculated by sliding window analysis, 8 mutation hotspot regions were detected, and coding regions were more conserved than non-coding regions. Analysis of non-synonymous substitution (Ka) and synonymous substitution (Ks) revealed that accD, rps12, petB, and atpF genes were affected by positive selection and no genes were affected by neutral selection. Based on the findings of the phylogenetic analysis, Artemisia selengensis was sister to the genus Artemisia Chrysanthemum and formed a monophyletic group with other Artemisia genera.

CONCLUSIONS

In this research, the present study systematically compared the chloroplast genomic features of A. selengensis and provided important information for the study of the chloroplast genome of A. selengensis and the evolutionary relationships among Asteraceae species.

Population genomics of the invasive Northern Giant Hornet Vespa mandarinia in North America and across its native range

The northern giant hornet Vespa mandarinia (NGH) is a voracious predator of other insect species, including honey bees. NGH's native range spans subtropical and temperate regions across much of east and southeast Asia and, in 2019, exotic populations of the species were discovered in North America. Despite this broad range and invasive potential, investigation of the population genomic structure of NGH across its native and introduced ranges has thus far been limited to a small number of mitochondrial samples. Here, we present analyses of genomic data from NGH individuals collected across the species' native range and from exotic individuals collected in North America. We provide the first survey of whole-genome population variation for any hornet species, covering this species' native and invasive ranges, and in doing so confirm likely origins in Japan and South Korea for the two introductions. We additionally show that, while this introduced population exhibited strongly elevated levels of inbreeding, these signatures of inbreeding are also present in some long-standing native populations, which may indicate that inbreeding depression alone is insufficient to prevent the persistence of NGH populations. As well as highlighting the importance of ongoing monitoring and eradication efforts to limit the spread of this species outside of its natural range, our data will serve as a foundational database for future genomic studies into introduced hornet populations.

The complete chloroplast genome assembly of Solidago altissima (Lineaus, 1753) (Astereae; Asteraceae)

This study aims to report the complete chloroplast genome of Solidago altissima L., a globally recognized invasive plant. The complete genome length of S. altissima is 152,961 bp; S. altissima has a typical quadripartite structure (including a large single copy of 84,829 bp, a small single copy of 18,084 bp, and two inverted repeat regions of 25,024 bp), which is commonly found in angiosperms. The genome contains 129 genes, consisting of 85 coding sequences, 36 tRNA genes, and 8 rRNA genes. To understand the phylogenetic relationship between S. altissima and its related species, maximum likelihood analysis was performed. The results revealed that S. altissima is closely related to Symphyotrichum subulatum. The findings of the present study could provide fundamental data for the future phylogenetic and evolutionary studies, while also research on species invasion and resolving complexity of S. altissima.

Chemical Relationship among Genetically Authenticated Medicinal Species of Genus Angelica

The genus Angelica comprises various species utilized for diverse medicinal purposes, with differences attributed to the varying levels or types of inherent chemical components in each species. This study employed DNA barcode analysis and HPLC analysis to genetically authenticate and chemically classify eight medicinal Angelica species (n = 106) as well as two non-medicinal species (n = 14) that have been misused. Nucleotide sequence analysis of the nuclear internal transcribed spacer (ITS) region revealed differences ranging from 11 to 117 bp, while psbA-trnH showed variances of 3 to 95 bp, respectively. Phylogenetic analysis grouped all samples except Angelica sinensis into the same cluster, with some counterfeits forming separate clusters. Verification using the NCBI database confirmed the feasibility of species identification. For chemical identification, a robust quantitative HPLC analysis method was developed for 46 marker compounds. Subsequently, two A. reflexa-specific and seven A. biserrata-specific marker compounds were identified, alongside non-specific markers. Moreover, chemometric clustering analysis reflecting differences in chemical content between species revealed that most samples formed distinct clusters according to the plant species. However, some samples formed mixed clusters containing different species. These findings offer crucial insights for the standardization and quality control of medicinal Angelica species.

Solanum aculeatissimum and Solanum torvum chloroplast genome sequences: a comparative analysis with other Solanum chloroplast genomes

BACKGROUND

Solanum aculeatissimum and Solanum torvum belong to the Solanum species, and they are essential plants known for their high resistance to diseases and adverse conditions. They are frequently used as rootstocks for grafting and are often crossbred with other Solanum species to leverage their resistance traits. However, the phylogenetic relationship between S. aculeatissimum and S. torvum within the Solanum genus remains unclear. Therefore, this paper aims to sequence the complete chloroplast genomes of S. aculeatissimum and S. torvum and analyze them in comparison with 29 other previously published chloroplast genomes of Solanum species.

RESULTS

We observed that the chloroplast genomes of S. aculeatissimum and S. torvum possess typical tetrameric structures, consisting of one Large Single Copy (LSC) region, two reverse-symmetric Inverted Repeats (IRs), and one Small Single Copy (SSC) region. The total length of these chloroplast genomes ranged from 154,942 to 156,004 bp, with minimal variation. The highest GC content was found in the IR region, while the lowest was in the SSC region. Regarding gene content, the total number of chloroplast genes and CDS genes remained relatively consistent, ranging from 128 to 134 and 83 to 91, respectively. Nevertheless, there was notable variability in the number of tRNA genes and rRNAs. Relative synonymous codon usage (RSCU) analysis revealed that both S. aculeatissimum and S. torvum preferred codons that utilized A and U bases. Analysis of the IR boundary regions indicated that contraction and expansion primarily occurred at the junction between SSC and IR regions. Nucleotide polymorphism analysis and structural variation analysis demonstrated that chloroplast variation in Solanum species mainly occurred in the LSC and SSC regions. Repeat sequence analysis revealed that A/T was the most frequent base pair in simple repeat sequences (SSR), while Palindromic and Forward repeats were more common in long sequence repeats (LSR), with Reverse and Complement repeats being less frequent. Phylogenetic analysis indicated that S. aculeatissimum and S. torvum belonged to the same meristem and were more closely related to Cultivated Eggplant.

CONCLUSION

These findings enhance our comprehension of chloroplast genomes within the Solanum genus, offering valuable insights for plant classification, evolutionary studies, and potential molecular markers for species identification.

The complete chloroplast genome assembly of Commelina caroliniana Walter (Commelinaceae)

We sequenced and published the chloroplast genome of Commelina caroliniana Walter, which was previously misidentified as C. diffusa owing to their morphological similarities until 1989. The genome of C. caroliniana is 160,857 bp long [large single copy region: 88,064 bp; a small single copy region: 18,549 bp; two inverted repeat regions: 27,122 bp] and has a GC content of 35.7%. The genome comprises 133 genes, including 87 coding sequences (CDSs), 38 tRNAs, and eight rRNAs. The phylogenetic relationship between C. caroliniana and related species was analyzed using the maximum likelihood method based on the 79 CDSs of the chloroplast genome. Phylogenetic analysis revealed that C. caroliniana is closely related to C. communis. Our findings will contribute to studies on species identification and phylogenetic and evolutionary research. These results enhance our understanding of the Commelina genus.

Intraspecific and Intrageneric Genomic Variation across Three Sedum Species (Crassulaceae): A Plastomic Perspective

Sedum is the largest succulent genus in Crassulaceae. Because of predominant maternal inheritance, little recombination, and slow evolution, plastomes can serve as powerful super barcodes for inter- or intra-species phylogenetic analyses. While previous research has focused on plastomes between Sedum species, intra-species studies are scarce. Here, we sequenced plastomes from three Sedum species (Sedum alfredii, Sedum plumbizincicola, and Sedum japonicum) to understand their evolutionary relationships and plastome structural evolution. Our analyses revealed minimal size and GC content variation across species. However, gene distribution at IR boundaries, repeat structures, and codon usage patterns showed diversity at both inter-specific and intra-specific levels. Notably, an rps19 gene expansion and a bias toward A/T-ending codons were observed. Codon aversion motifs also varied, potentially serving as markers for future studies. Phylogenetic analyses confirmed the non-monophyly of Sedum and divided the Acre clade into two groups. Individuals from the same species clustered together, with strong support for the relationships between S. alfredii, S. tricarpum, and S. plumbizincicola. Additionally, S. japonicum clearly affiliates with the Acre clade. This study provides valuable insights into both intra-specific and intra-generic plastome variation in Sedum, as well as overall plastome evolution within the genus.

Three complete chloroplast genomes from two north American Rhus species and phylogenomics of Anacardiaceae

BACKGROUND

The suamc genus Rhus (sensu stricto) includes two subgenera, Lobadium (ca. 25 spp.) and Rhus (ca. 10 spp.). Their members, R. glabra and R. typhina (Rosanae: Sapindales: Anacardiaceae), are two economic important species. Chloroplast genome information is of great significance for the study of plant phylogeny and taxonomy.

RESULTS

The three complete chloroplast genomes from two Rhus glabra and one R. typhina accessions were obtained with a total of each about 159k bp in length including a large single-copy region (LSC, about 88k bp), a small single-copy regions (SSC, about 19k bp) and a pair of inverted repeats regions (IRa/IRb, about 26k bp), to form a canonical quadripartite structure. Each genome contained 88 protein-coding genes, 37 transfer RNA genes, eight ribosomal RNA genes and two pseudogenes. The overall GC content of the three genomes all were same (37.8%), and RSCU values showed that they all had the same codon prefers, i.e., to use codon ended with A/U (93%) except termination codon. Three variable hotspots, i.e., ycf4-cemA, ndhF-rpl32-trnL and ccsA-ndhD, and a total of 152-156 simple sequence repeats (SSR) were identified. The nonsynonymous (Ka)/synonymous (Ks) ratio was calculated, and cemA and ycf2 genes are important indicators of gene evolution. The phylogenetic analyses of the family Anacardiaceae showed that the eight genera were grouped into three clusters, and supported the monophyly of the subfamilies and all the genera. The accessions of five Rhus species formed four clusters, while, one individual of R. typhina grouped with the R. glabra accessions instead of clustering into the two other individuals of R. typhina in the subgenus Rhus, which showed a paraphyletic relationship.

CONCLUSIONS

Comparing the complete chloroplast genomes of the Rhus species, it was found that most SSRs were A/T rich and located in the intergenic spacer, and the nucleotide divergence exhibited higher levels in the non-coding region than in the coding region. The Ka/Ks ratio of cemA gene was > 1 for species collected in America, while it was < 1 for other species in China, which dedicated that the Rhus species from North America and East Asia have different evolutionary pressure. The phylogenetic analysis of the complete chloroplast genome clarified the Rhus placement and relationship. The results obtained in this study are expected to provide valuable genetic resources to perform species identification, molecular breeding, and intraspecific diversity of the Rhus species.

Transcriptome and comparative chloroplast genome analysis of Taxus yunnanensis individuals with high and low paclitaxel yield

Paclitaxel is a potent anti-cancer drug that is mainly produced through semi-synthesis, which still requires plant materials as precursors. The content of paclitaxel and 10-deacetyl baccatin III (10-DAB) in Taxus yunnanensis has been found to differ from that of other Taxus species, but there is little research on the mechanism underlying the variation in paclitaxel content in T. yunnanensis of different provenances. In this experiment, the contents of taxoids and precursors in twigs between a high paclitaxel-yielding individual (TG) and a low paclitaxel-yielding individual (TD) of T. yunnanensis were compared, and comparative analyses of transcriptomes as well as chloroplast genomes were performed. High-performance liquid chromatography (HPLC) detection showed that 10-DAB and baccatin III contents in TG were 18 and 47 times those in TD, respectively. Transcriptomic analysis results indicated that genes encoding key enzymes in the paclitaxel biosynthesis pathway, such as taxane 10-β-hydroxylase (T10βH), 10-deacetylbaccatin III 10-O-acetyltransferase (DBAT), and debenzoyl paclitaxel N-benzoyl transferase (DBTNBT), exhibited higher expression levels in TG. Additionally, qRT-PCR showed that the relative expression level of T10βH and DBAT in TG were 29 and 13 times those in TD, respectively. In addition, six putative transcription factors were identified that may be involved in paclitaxel biosynthesis from transcriptome data. Comparative analysis of plastid genomes showed that the TD chloroplast contained a duplicate of rps12, leading to a longer plastid genome length in TD relative to TG. Fifteen mutation hotspot regions were identified between the two plastid genomes that can serve as candidate DNA barcodes for identifying high-paclitaxel-yield individuals. This experiment provides insight into the difference in paclitaxel accumulation among different provenances of T. yunnanensis individuals.

Complete chloroplast genome and comparison of herbicides toxicity on Aeschynomene indica (Leguminosae) in upland direct-seeding paddy field

BACKGROUND

Indian jointvetch (Aeschynomene indica) is a common and pernicious weed found in the upland direct-seeding rice fields in the lower reaches of the Yangtze River in China. However, there are few reports on the degree of harm, genetic characteristics, and management methods of this weed. The purpose of this study is to clarify the harm of Indian jointvetch to upland direct-seeding rice, analyze the genetic characteristics of this weed based on chloroplast genomics and identify its related species, and screen herbicides that are effective in managing this weed in upland direct-seeding rice fields.

RESULTS

In a field investigation in upland direct-seeding rice paddies in Shanghai and Jiangsu, we determined that the plant height and maximum lateral distance of Indian jointvetch reached approximately 134.2 cm and 57.9 cm, respectively. With Indian jointvetch present at a density of 4/m2 and 8/m2, the yield of rice decreased by approximately 50% and 70%, respectively. We further obtained the first assembly of the complete chloroplast (cp.) genome sequence of Indian jointvetch (163,613 bp). There were 161 simple sequence repeats, 166 long repeats, and 83 protein-encoding genes. The phylogenetic tree and inverted repeat region expansion and contraction analysis based on cp. genomes demonstrated that species with closer affinity to A. indica included Glycine soja, Glycine max, and Sesbania cannabina. Moreover, a total of 3281, 3840, and 3838 single nucleotide polymorphisms were detected in the coding sequence regions of the cp. genomes of S. cannabina voucher IBSC, G. soja, and G. max compared with the A. indica sequence, respectively. A greenhouse pot experiment indicated that two pre-emergence herbicides, saflufenacil and oxyfluorfen, and two post-emergence herbicides, florpyrauxifen-benzyl and penoxsulam, can more effectively manage Indian jointvetch than other common herbicides in paddy fields. The combination of these two types of herbicides is recommended for managing Indian jointvetch throughout the entire growth period of upland direct-seeding rice.

CONCLUSIONS

This study provides molecular resources for future research focusing on the identification of the infrageneric taxa, phylogenetic resolution, and biodiversity of Leguminosae plants, along with recommendations for reliable management methods to control Indian jointvetch.

Characterization of the complete chloroplast genome of Wolffia arrhiza and comparative genomic analysis with relative Wolffia species

Lemnoideae, commonly referred to as the duckweed, are aquatic plants found worldwide. Wolffia species are known for their extreme reduction in size and complexity, lacking both roots and leaves, and they hold the distinction of being the smallest plants among angiosperms. Interestingly, it belongs to the Araceae family, despite its apparent morphological differences from land plants in the same family. Traditional morphological methods have limitations in classifying these plants, making molecular-level information essential. The chloroplast genome of Wolffia arrhiza is revealed that a total length of 169,602 bp and a total GC content of 35.78%. It follows the typical quadripartite structure, which includes a large single copy (LSC, 92,172 bp) region, a small single copy (SSC, 13,686 bp) region, and a pair of inverted repeat (IR, 31,872 bp each) regions. There are 131 genes characterized, comprising 86 Protein-Coding Genes, 37 Transfer RNA (tRNA) genes, and 8 ribosomal RNA (rRNA) genes. Moreover, 48 simple sequence repeats and 32 long repeat sequences were detected. Comparative analysis between W. arrhiza and six other Lemnoideae species identified 12 hotspots of high nucleotide diversity. In addition, a phylogenetic analysis was performed using 14 species belonging to the Araceae family and one external species as an outgroup. This analysis unveiled W. arrhiza and Wolffia globosa as closely related sister species. Therefore, this research has revealed the complete chloroplast genome data of W. arrhiza, offering a more detailed understanding of its evolutionary position and phylogenetic categorization within the Lemnoideae subfamily.

Species delimitation of tea plants (Camellia sect. Thea) based on super-barcodes

BACKGROUND

The era of high throughput sequencing offers new paths to identifying species boundaries that are complementary to traditional morphology-based delimitations. De novo species delimitation using traditional or DNA super-barcodes serve as efficient approaches to recognizing putative species (molecular operational taxonomic units, MOTUs). Tea plants (Camellia sect. Thea) form a group of morphologically similar species with significant economic value, providing the raw material for tea, which is the most popular nonalcoholic caffeine-containing beverage in the world. Taxonomic challenges have arisen from vague species boundaries in this group.

RESULTS

Based on the most comprehensive sampling of C. sect. Thea by far (165 individuals of 39 morphospecies), we applied three de novo species delimitation methods (ASAP, PTP, and mPTP) using plastome data to provide an independent evaluation of morphology-based species boundaries in tea plants. Comparing MOTU partitions with morphospecies, we particularly tested the congruence of MOTUs resulting from different methods. We recognized 28 consensus MOTUs within C. sect. Thea, while tentatively suggesting that 11 morphospecies be discarded. Ten of the 28 consensus MOTUs were uncovered as morphospecies complexes in need of further study integrating other evidence. Our results also showed a strong imbalance among the analyzed MOTUs in terms of the number of molecular diagnostic characters.

CONCLUSION

This study serves as a solid step forward for recognizing the underlying species boundaries of tea plants, providing a needed evidence-based framework for the utilization and conservation of this economically important plant group.

Comparative Analysis of Chloroplast Genome of Meconopsis (Papaveraceae) Provides Insights into Their Genomic Evolution and Adaptation to High Elevation

The Meconopsis species are widely distributed in the Qinghai-Tibet Plateau, Himalayas, and Hengduan Mountains in China, and have high medicinal and ornamental value. The high diversity of plant morphology in this genus poses significant challenges for species identification, given their propensity for highland dwelling, which makes it a question worth exploring how they cope with the harsh surroundings. In this study, we recently generated chloroplast (cp) genomes of two Meconopsis species, Meconopsis paniculata (M. paniculata) and M. pinnatifolia, and compared them with those of ten Meconopsis cp genomes to comprehend cp genomic features, their phylogenetic relationships, and what part they might play in plateau adaptation. These cp genomes shared a great deal of similarities in terms of genome size, structure, gene content, GC content, and codon usage patterns. The cp genomes were between 151,864 bp and 154,997 bp in length, and contain 133 predictive genes. Through sequence divergence analysis, we identified three highly variable regions (trnD-psbD, ccsA-ndhD, and ycf1 genes), which could be used as potential markers or DNA barcodes for phylogenetic analysis. Between 22 and 38 SSRs and some long repeat sequences were identified from 12 Meconopsis species. Our phylogenetic analysis confirmed that 12 species of Meconopsis clustered into a monophyletic clade in Papaveraceae, which corroborated their intrageneric relationships. The results indicated that M. pinnatifolia and M. paniculata are sister species in the phylogenetic tree. In addition, the atpA and ycf2 genes were positively selected in high-altitude species. The functions of these two genes might be involved in adaptation to the extreme environment in the cold and low CO2 concentration conditions at the plateau.

Whole-genome sequencing in medicinal plants: current progress and prospect

Advancements in genomics have dramatically accelerated the research on medicinal plants, and the development of herbgenomics has promoted the "Project of 1K Medicinal Plant Genome" to decipher their genetic code. However, it is difficult to obtain their high-quality whole genomes because of the prevalence of polyploidy and/or high genomic heterozygosity. Whole genomes of 123 medicinal plants were published until September 2022. These published genome sequences were investigated in this review, covering their classification, research teams, ploidy, medicinal functions, and sequencing strategies. More than 1,000 institutes or universities around the world and 50 countries are conducting research on medicinal plant genomes. Diploid species account for a majority of sequenced medicinal plants. The whole genomes of plants in the Poaceae family are the most studied. Almost 40% of the published papers studied species with tonifying, replenishing, and heat-cleaning medicinal effects. Medicinal plants are still in the process of domestication as compared with crops, thereby resulting in unclear genetic backgrounds and the lack of pure lines, thus making their genomes more difficult to complete. In addition, there is still no clear routine framework for a medicinal plant to obtain a high-quality whole genome. Herein, a clear and complete strategy has been originally proposed for creating a high-quality whole genome of medicinal plants. Moreover, whole genome-based biological studies of medicinal plants, including breeding and biosynthesis, were reviewed. We also advocate that a research platform of model medicinal plants should be established to promote the genomics research of medicinal plants.

Semi-automated sequence curation for reliable reference datasets in ITS2 vascular plant DNA (meta-)barcoding

One of the most critical steps for accurate taxonomic identification in DNA (meta)-barcoding is to have an accurate DNA reference sequence dataset for the marker of choice. Therefore, developing such a dataset has been a long-term ambition, especially in the Viridiplantae kingdom. Typically, reference datasets are constructed with sequences downloaded from general public databases, which can carry taxonomic and other relevant errors. Herein, we constructed a curated (i) global dataset, (ii) European crop dataset, and (iii) 27 datasets for the EU countries for the ITS2 barcoding marker of vascular plants. To that end, we first developed a pipeline script that entails (i) an automated curation stage comprising five filters, (ii) manual taxonomic correction for misclassified taxa, and (iii) manual addition of newly sequenced species. The pipeline allows easy updating of the curated datasets. With this approach, 13% of the sequences, corresponding to 7% of species originally imported from GenBank, were discarded. Further, 259 sequences were manually added to the curated global dataset, which now comprises 307,977 sequences of 111,382 plant species.

Phylogeny and diversification of genus Sanicula L. (Apiaceae): novel insights from plastid phylogenomic analyses

BACKGROUND

The genus Sanicula L. is a unique perennial herb that holds important medicinal values. Although the previous studies on Sanicula provided us with a good research basis, its taxonomic system and interspecific relationships have not been satisfactorily resolved, especially for those endemic to China. Moreover, the evolutionary history of this genus also remains inadequately understood. The plastid genomes possessing highly conserved structure and limited evolutionary rate have proved to be an effective tool for studying plant phylogeny and evolution.

RESULTS

In the current study, we newly sequenced and assembled fifteen Sanicula complete plastomes. Combined with two previously reported plastomes, we performed comprehensively plastid phylogenomics analyses to gain novel insights into the evolutionary history of this genus. The comparative results indicated that the seventeen plastomes exhibited a high degree of conservation and similarity in terms of their structure, size, GC content, gene order, IR borders, codon bias patterns and SSRs profiles. Such as all of them displayed a typical quadripartite structure, including a large single copy region (LSC: 85,074-86,197 bp), a small single copy region (SSC: 17,047-17,132 bp) separated by a pair of inverted repeat regions (IRs: 26,176-26,334 bp). And the seventeen plastomes had similar IR boundaries and the adjacent genes were identical. The rps19 gene was located at the junction of the LSC/IRa, the IRa/SSC junction region was located between the trnN gene and ndhF gene, the ycf1 gene appeared in the SSC/IRb junction and the IRb/LSC boundary was located between rpl12 gene and trnH gene. Twelve specific mutation hotspots (atpF, cemA, accD, rpl22, rbcL, matK, ycf1, trnH-psbA, ycf4-cemA, rbcL-accD, trnE-trnT and trnG-trnR) were identified that can serve as potential DNA barcodes for species identification within the genus Sanicula. Furthermore, the plastomes data and Internal Transcribed Spacer (ITS) sequences were performed to reconstruct the phylogeny of Sanicula. Although the tree topologies of them were incongruent, both provided strong evidence supporting the monophyly of Saniculoideae and Apioideae. In addition, the sister groups between Saniculoideae and Apioideae were strongly suggested. The Sanicula species involved in this study were clustered into a clade, and the Eryngium species were also clustered together. However, it was clearly observed that the sections of Sanicula involved in the current study were not respectively recovered as monophyletic group. Molecular dating analysis explored that the origin of this genus was occurred during the late Eocene period, approximately 37.84 Ma (95% HPD: 20.33-52.21 Ma) years ago and the diversification of the genus was occurred in early Miocene 18.38 Ma (95% HPD: 10.68-25.28 Ma).

CONCLUSION

The plastome-based tree and ITS-based tree generated incongruences, which may be attributed to the event of hybridization/introgression, incomplete lineage sorting (ILS) and chloroplast capture. Our study highlighted the power of plastome data to significantly improve the phylogenetic supports and resolutions, and to efficiently explore the evolutionary history of this genus. Molecular dating analysis explored that the diversification of the genus occurred in the early Miocene, which was largely influenced by the prevalence of the East Asian monsoon and the uplift of the Hengduan Mountains (HDM). In summary, our study provides novel insights into the plastome evolution, phylogenetic relationships, taxonomic framework and evolution of genus Sanicula.

CMAUP database update 2024: extended functional and association information of useful plants for biomedical research

Knowledge of the collective activities of individual plants together with the derived clinical effects and targeted disease associations is useful for plant-based biomedical research. To provide the information in complement to the established databases, we introduced a major update of CMAUP database, previously featured in NAR. This update includes (i) human transcriptomic changes overlapping with 1152 targets of 5765 individual plants, covering 74 diseases from 20 027 patient samples; (ii) clinical information for 185 individual plants in 691 clinical trials; (iii) drug development information for 4694 drug-producing plants with metabolites developed into approved or clinical trial drugs; (iv) plant and human disease associations (428 737 associations by target, 220 935 reversion of transcriptomic changes, 764 and 154121 associations by clinical trials of individual plants and plant ingredients); (v) the location of individual plants in the phylogenetic tree for navigating taxonomic neighbors, (vi) DNA barcodes of 3949 plants, (vii) predicted human oral bioavailability of plant ingredients by the established SwissADME and HobPre algorithm, (viii) 21-107% increase of CMAUP data over the previous version to cover 60 222 chemical ingredients, 7865 plants, 758 targets, 1399 diseases, 238 KEGG human pathways, 3013 gene ontologies and 1203 disease ontologies. CMAUP update version is freely accessible at https://bidd.group/CMAUP/index.html.

Exploring DNA barcode for accurate identification of threatened Aconitum L. species from Western Himalaya

BACKGROUND

Aconitum species, belonging to Ranunculaceae, have high medicinal importance but due to their overexploitation come under IUCN (International Union for Conservation of Nature) red list. The precise identification of the Aconitum species is equally important because they are used in herbal formulations. The present study aimed to develop an efficient DNA barcode system for the authentic identification of Aconitum species.

METHODS AND RESULTS

A set of 92 barcode gene sequences (including 12 developed during the present study and 80 retrieved from NCBI) of 5 Aconitum species (A. heterophyllum, A. vialoceum, A. japonicum, A. napellus, and A. stapfianum) were analyzed using three methods (tree-based, distance-based, and similarity-based) for species discrimination. The PWG-distance method was found most effective for species discrimination. The discrimination rate of PWG- distance ranged from 33.3% (rbcL + trnH-psbA) to 100% (ITS, rbcL + ITS, ITS + trnH-psbA and rbcL + ITS + trnH-psbA). Among DNA barcodes and their combinations, the ITS marker had the highest degree of species discrimination (NJ-40%, PWG-100% and BLAST-40%), followed by trnH-psbA (NJ-20%, PWG-60% and BLAST-20%). ITS also had higher barcoding gap as compared to other individual barcodes and their combinations. Further, we also analyzed six Aconitum species (A. balfourii, A. ferox, A. heterophyllum, A. rotundifolium, A. soongaricum and A. violaceum) existing in Western Himalaya. These species were distinguished clearly through tree-based method using the ITS barcode gene with 100% species resolution.

CONCLUSION

ITS showed the best species discrimination power and was used to develop species-specific barcodes for Aconitum species. DNA barcodes developed during the present study can be used to identify Aconitum species.

DNA Barcoding and Metabarcoding Protocols for Species Identification

The article presents the several steps to be performed on a plant, fungal, insect, or soil sample to obtain DNA sequences for DNA barcode analysis. The chapter begins with a description of sample preparation including procedures for cleaning and proceeds to DNA extraction with methods adapted for the specific type of sample. Next, DNA quantification is described so the proper amount is used for the amplification of the selected barcode regions. Information is provided for reaction mixes and amplification conditions for several referenced barcode primer pairs tuned for the individual sample of interest. This is followed by a description of procedures to access the success of amplification, cleanup, and quantification of the product ready for either Sanger sequencing or library preparation for massive parallel sequencing (MPS). Finally, procedures are provided for Sanger sequencing, library preparation, and MPS sequencing. The chapter provides several references of barcode regions for different sample types.

Stepping towards pollen DNA metabarcoding: A breakthrough in forensic sciences

This review is engaged in determining the capability of plant pollen as a significant source of evidence for the linkage between suspects and crime location in forensic sciences. Research and review articles were collected from Google Scholar, the Web of Science, and PubMed. Articles were searched using specific keywords such as "Forensic Palynology," "Pollen metabarcoding," "Plant forensics," and "Pollen" AND "criminal investigation." Boolean logic was also utilized to narrow the articles to be included in this review article. Through the literature and exploratory research, it has been observed in the current study that with advancements in technology, forensic palynology has found its application in creating an association between the crime scene and suspected individuals to have a link to it, as pollen DNA is a long-lasting investigative tool that can effectively help forensic investigations. Moreover, the literature shows that the DNA of pollen and spores has helped forensic scientists link suspects to crime scenes, and the introduction of pollen DNA metabarcoding tools has eased the efforts of palynologists to analyze pollen DNA. The introduction of DNA metabarcoding techniques to analyze pollen from plants has helped identify the geological locations of the plants and ultimately identify the culprit.

Insights into the phylogenetic relationships and species boundaries of the Myricaria squamosa complex (Tamaricaceae) based on the complete chloroplast genome

Myricaria plants are widely distributed in Eurasia and are helpful for windbreak and embankment protection. Current molecular evidence has led to controversy regarding species boundaries within the Myricaria genus and interspecific phylogenetic relationships between three specific species-M. bracteata, M. paniculata and M. squamosa-which have remained unresolved. This study treated these three unresolved taxa as a species complex, named the M. squamosa complex. The genome skimming approach was used to determine 35 complete plastome sequences and nuclear ribosomal DNA sequences for the said complex and other closely related species, followed by de novo assembly. Comparative analyses were conducted across Myricaria to identify the genome size, gene content, repeat type and number, SSR (simple sequence repeat) abundance, and codon usage bias of chloroplast genomes. Tree-based species delimitation results indicated that M. bracteata, M. paniculata and M. squamosa could not be distinguished and formed two monophyletic lineages (P1 and P2) that were clustered together. Compared to plastome-based species delimitation, the standard nuclear DNA barcode had the lowest species resolution, and the standard chloroplast DNA barcode and group-specific barcodes delimitated a maximum of four out of the five species. Plastid phylogenomics analyses indicated that the monophyletic M. squamosa complex is comprised of two evolutionarily significant units: one in the western Tarim Basin and the other in the eastern Qinghai-Tibet Plateau. This finding contradicts previous species discrimination and promotes the urgent need for taxonomic revision of the threatened genus Myricaria. Dense sampling and plastid genomes will be essential in this effort. The super-barcodes and specific barcode candidates outlined in this study will aid in further studies of evolutionary history.

The low level of plastome differentiation observed in some lineages of Poales hinders molecular species identification

Chloroplast genomes are a source of information successfully used in various fields of plant genetics, including molecular species identification. However, recent studies indicate an extremely low level of interspecific variability in the plastomes of some taxonomic groups of plants, including the genus Stipa L., which is a representative of the grass family. In this study we aimed to analyze the level of chloroplast genome diversity within particular genera as well as the effectiveness of identifying plant species in the Poaceae family and the other representatives of Poales order. Analysis of complete plastid genome alignments created for 96 genera comprising 793 species and 1707 specimens obtained from the GenBank database allowed defining and categorizing molecular diagnostic characters distinguishing the analyzed species from the other representatives of the genus. The results also demonstrate which species do not have any species-specific mutations, thereby they cannot be identified on the basis of differences between the complete chloroplast genomes. Our research showed a huge diversity of the analyzed species in terms of the number of molecular diagnostic characters and indicated which genera pose a particular challenge in terms of molecular species identification. The results show that a very low level of genetic diversity between plastomes is not uncommon in Poales. This is the first extensive research on super-barcoding that tests this method on a large data set and illustrates its effectiveness against the background of phylogenetic relationships.

Comparative and phylogenetic analysis of Chiloschista (Orchidaceae) species and DNA barcoding investigation based on plastid genomes

BACKGROUND

Chiloschista (Orchidaceae, Aeridinae) is an epiphytic leafless orchid that is mainly distributed in tropical or subtropical forest canopies. This rare and threatened orchid lacks molecular resources for phylogenetic and barcoding analysis. Therefore, we sequenced and assembled seven complete plastomes of Chiloschista to analyse the plastome characteristics and phylogenetic relationships and conduct a barcoding investigation.

RESULTS

We are the first to publish seven Chiloschista plastomes, which possessed the typical quadripartite structure and ranged from 143,233 bp to 145,463 bp in size. The plastomes all contained 120 genes, consisting of 74 protein-coding genes, 38 tRNA genes and eight rRNA genes. The ndh genes were pseudogenes or lost in the genus, and the genes petG and psbF were under positive selection. The seven Chiloschista plastomes displayed stable plastome structures with no large inversions or rearrangements. A total of 14 small inversions (SIs) were identified in the seven Chiloschista plastomes but were all similar within the genus. Six noncoding mutational hotspots (trnNGUU-rpl32 > rpoB-trnCGCA > psbK-psbI > psaC-rps15 > trnEUUC-trnTGGU > accD-psaI) and five coding sequences (ycf1 > rps15 > matK > psbK > ccsA) were selected as potential barcodes based on nucleotide diversity and species discrimination analysis, which suggested that the potential barcode ycf1 was most suitable for species discrimination. A total of 47-56 SSRs and 11-14 long repeats (> 20 bp) were identified in Chiloschista plastomes, and they were mostly located in the large single copy intergenic region. Phylogenetic analysis indicated that Chiloschista was monophyletic. It was clustered with Phalaenopsis and formed the basic clade of the subtribe Aeridinae with a moderate support value. The results also showed that seven Chiloschista species were divided into three major clades with full support.

CONCLUSION

This study was the first to analyse the plastome characteristics of the genus Chiloschista in Orchidaceae, and the results showed that Chiloschista plastomes have conserved plastome structures. Based on the plastome hotspots of nucleotide diversity, several genes and noncoding regions are suitable for phylogenetic and population studies. Chiloschista may provide an ideal system to investigate the dynamics of plastome evolution and DNA barcoding investigation for orchid studies.

Assessment of genetic diversity and phylogenetic relationship of local coffee populations in southwestern Saudi Arabia using DNA barcoding

The genetic diversity of local coffee populations is crucial to breed new varieties better adapted to the increasingly stressful environment due to climate change and evolving consumer preferences. Unfortunately, local coffee germplasm conservation and genetic assessment have not received much attention. Molecular tools offer substantial benefits in identifying and selecting new cultivars or clones suitable for sustainable commercial utilization. New annotation methods, such as chloroplast barcoding, are necessary to produce accurate and high-quality phylogenetic analyses. This study used DNA barcoding techniques to examine the genetic relationships among fifty-six accessions collected from the southwestern part of Saudi Arabia. PCR amplification and sequence characterization were used to investigate the effectiveness of four barcoding loci: atpB-rbcl, trnL-trnF, trnT-trnL, and trnL. The maximum nucleotide sites, nucleotide diversity, and an average number of nucleotide differences were recorded for atpB-rbcl, while trnT-trnL had the highest variable polymorphic sites, segregating sites, and haploid diversity. Among the four barcode loci, trnT-trnL recorded the highest singleton variable sites, while trnL recorded the highest parsimony information sites. Furthermore, the phylogenetic analysis clustered the Coffea arabica genotypes into four different groups, with three genotypes (KSA31, KSA38, and KSA46) found to be the most divergent genotypes standing alone in the cluster and remained apart during the analysis. The study demonstrates the presence of considerable diversity among coffee populations in Saudi Arabia. Furthermore, it also shows that DNA barcoding is an effective technique for identifying local coffee genotypes, with potential applications in coffee conservation and breeding efforts.

A DNA barcode library for woody plants in tropical and subtropical China

The application of DNA barcoding has been significantly limited by the scarcity of reliable specimens and inadequate coverage and replication across all species. The deficiency of DNA barcode reference coverage is particularly striking for highly biodiverse subtropical and tropical regions. In this study, we present a comprehensive barcode library for woody plants in tropical and subtropical China. Our dataset includes a standard barcode library comprising the four most widely used barcodes (rbcL, matK, ITS, and ITS2) for 2,520 species from 4,654 samples across 49 orders, 144 families, and 693 genera, along with 79 samples identified at the genus level. This dataset also provides a super-barcode library consisting of 1,239 samples from 1,139 species, 411 genera, 113 families, and 40 orders. This newly developed library will serve as a valuable resource for DNA barcoding research in tropical and subtropical China and bordering countries, enable more accurate species identification, and contribute to the conservation and management of tropical and subtropical forests.

Authentication of Allium ulleungense, A. microdictyon and A. ochotense based on super-barcoding of plastid genome and 45S nrDNA

Allium ulleungense (AU) and A. microdictyon (AM) are valuable medicinal and edible vegetables, referred to as mountain garlic in Korea. The identification of AU, AM and a neighboring species A. ochotense (AO) is difficult because of their morphological similarities. We collected samples from three species and 46 cultivated collections to understand the genetic diversity of these valuable Allium species. Among them, we sequenced six collections, including three species and three cultivating collections to obtain data from the plastid genome (plastome) and nuclear 45S ribosomal DNA (nrDNA) for super-barcoding. The AM and AO showed around 60 single nucleotide polymorphisms (SNPs) and 39 Insertion/Deletion (InDels) in the plastome but no variations in the nrDNA sequences. Conversely, the AU and AM showed more than 170 SNPs and 80 InDels in the plastomes, and 20 SNPs and 1 InDel were found in the 45S nrDNA sequences. Among the three cultivating collections, one TB collection was determined to be the AU type in both plastome and nrDNA sequences. However, the other two collections, JB and SA, showed the AM type plastome but were heterozygous in the 45S nrDNA sequences, indicating both AU and AM types (putative AM x AU hybrid). Ten molecular markers were developed based on sequence variations to identify these three species and assess their genetic diversity. A total of 49 collections were genotyped using the ten developed markers and classified into five groups: 14 AU, 22 AM, 1 AO, 3 putative AM x AU hybrids, and 9 putative AU x AM hybrid collections. Super-barcoding with plastomes and nrDNAs revealed the genetic diversity of the three Allium species and putative hybrids between species. The newly developed markers will facilitate species and hybrid identification, thereby benefiting marker-assisted molecular breeding of Allium species.

Creating, curating and evaluating a mitogenomic reference database to improve regional species identification using environmental DNA

Species detection using eDNA is revolutionizing global capacity to monitor biodiversity. However, the lack of regional, vouchered, genomic sequence information-especially sequence information that includes intraspecific variation-creates a bottleneck for management agencies wanting to harness the complete power of eDNA to monitor taxa and implement eDNA analyses. eDNA studies depend upon regional databases of mitogenomic sequence information to evaluate the effectiveness of such data to detect and identify taxa. We created the Oregon Biodiversity Genome Project to create a database of complete, nearly error-free mitogenomic sequences for all of Oregon's fishes. We have successfully assembled the complete mitogenomes of 313 specimens of freshwater, anadromous and estuarine fishes representing 24 families, 55 genera and 129 species and lineages. Comparative analyses of these sequences illustrate that many regions of the mitogenome are taxonomically informative, that the short (~150 bp) mitochondrial 'barcode' regions typically used for eDNA assays do not consistently diagnose for species and that complete single or multiple genes of the mitogenome are preferable for identifying Oregon's fishes. This project provides a blueprint for other researchers to follow as they build regional databases, illustrates the taxonomic value and limits of complete mitogenomic sequences and offers clues as to how current eDNA assays and environmental genomics methods of the future can best leverage this information.

Comparison of chloroplast genomes and phylogenetic analysis of four species in Quercus section Cyclobalanopsis

The identification in Quercus L. species was considered to be difficult all the time. The fundamental phylogenies of Quercus have already been discussed by morphological and molecular means. However, the morphological characteristics of some Quercus groups may not be consistent with the molecular results (such as the group Helferiana), which may lead to blurring of species relationships and prevent further evolutionary researches. To understand the interspecific relationships and phylogenetic positions, we sequenced and assembled the CPGs (160,715 bp-160842 bp) of four Quercus section Cyclobalanopsis species by Illumina pair-end sequencing. The genomic structure, GC content, and IR/SC boundaries exhibited significant conservatism. Six highly variable hotspots were detected in comparison analysis, among which rpoC1, clpP and ycf1 could be used as molecular markers. Besides, two genes (petA, ycf2) were detected to be under positive selection pressure. The phylogenetic analysis showed: Trigonobalanus genus and Fagus genus located at the base of the phylogeny tree; The Quercus genus species were distincted to two clades, including five sections. All Compound Trichome Base species clustered into a single branch, which was in accordance with the results of the morphological studies. But neither of group Gilva nor group Helferiana had formed a monophyly. Six Compound Trichome Base species gathered together in pairs to form three branch respectively (Quercus kerrii and Quercus chungii; Quercus austrocochinchinensis with Quercus gilva; Quercus helferiana and Quercus rex). Due to a low support rate (0.338) in the phylogeny tree, the interspecies relationship between the two branches differentiated by this node remained unclear. We believe that Q. helferiana and Q. kerrii can exist as independent species due to their distance in the phylogeny tree. Our study provided genetic information in Quercus genus, which could be applied to further studies in taxonomy and phylogenetics.

Comparative chloroplast genomes provided insights into the evolution and species identification on the Datureae plants

Generally, chloroplast genomes of angiosperms are always highly conserved but carry a certain number of variation among species. In this study, chloroplast genomes of 13 species from Datureae tribe that are of importance both in ornamental gardening and medicinal usage were studied. In addition, seven chloroplast genomes from Datureae together with two from Solanaceae species retrieved from the National Center for Biotechnology Information (NCBI) were integrated into this study. The chloroplast genomes ranged in size from 154,686 to 155,979 and from 155,497 to 155,919 bp for species of Datura and Brugmansia, respectively. As to Datura and Brugmansia, a total of 128 and 132 genes were identified, in which 83 and 87 protein coding genes were identified, respectively; Furthermore, 37 tRNA genes and 8 rRNA genes were both identified in Datura and Brugmansia. Repeats analysis indicated that the number and type varied among species for Simple sequence repeat (SSR), long repeats, and tandem repeats ranged in number from 53 to 59, 98 to 99, and 22 to 30, respectively. Phylogenetic analysis based on the plastid genomes supported the monophyletic relationship among Datura and Brugmansia and Trompettia, and a refined phylogenic relationships among each individual was resolved. In addition, a species-specific marker was designed based on variation spot that resulted from a comparative analysis of chloroplast genomes and verified as effective maker for identification of D. stramonium and D. stramonium var. inermis. Interestingly, we found that 31 genes were likely to be under positive selection, including genes encoding ATP protein subunits, photosystem protein subunit, ribosome protein subunits, NAD(P)H dehydrogenase complex subunits, and clpP, petB, rbcL, rpoCl, ycf4, and cemA genes. These genes may function as key roles in the adaption to diverse environment during evolution. The diversification of Datureae members was dated back to the late Oligocene periods. These chloroplast genomes are useful genetic resources for taxonomy, phylogeny, and evolution for Datureae.

The complete chloroplast genome of Cicer reticulatum and comparative analysis against relative Cicer species

The chloroplast (cp) genome is an adequate genomic resource to investigate evolutionary relationships among plant species and it carries marker genes available for species identification. The Cicer reticulatum is one of perennial species as the progenitor of cultivated chickpeas. Although a large part of the land plants has a quadruple chloroplast genome organization, the cp genome of C. reticulatum consists of one LSC (Large Single Copy Region), one SSC (Small Single Copy Region), and one IR (Inverted Repeat) region, which indicates that it has an untypical and unique structure. This type of chloroplast genome belongs to the IR-lacking clade. Chloroplast DNA (cpDNA) was extracted from fresh leaves using a high salt-based protocol and sequencing was performed using DNA Nanoball Sequencing technology. The comparative analysis employed between the species to examine genomic differences and gene homology. The study also included codon usage frequency analysis, hotspot divergence analysis, and phylogenetic analysis using various bioinformatics tools. The cp genome of C. reticulatum was found 125,794 bp in length, with an overall GC content of 33.9%. With a total of 79 protein-coding genes, 34 tRNA genes, and 4 rRNA genes. Comparative genomic analysis revealed 99.93% similarity between C. reticulatum and C. arietinum. Phylogenetic analysis further indicated that the closest evolutionary relative to C. arietinum was C. reticulatum, whereas the previously sequenced wild Cicer species displayed slight distinctions across their entire coding regions. Several genomic regions, such as clpP and ycf1, were found to exhibit high nucleotide diversity, suggesting their potential utility as markers for investigating the evolutionary relationships within the Cicer genus. The first complete cp genome sequence of C. reticulatum will provide novel insights for future genetic research on Cicer crops.

Comparative and phylogenetic analyses of the chloroplast genomes of Filipendula species (Rosoideae, Rosaceae)

Genus Filipendula (Rosoideae, Rosaceae) comprises about 15 species and mainly distributed in Northern Hemisphere. The phylogenetic relationships based on the nrITS marker are not consistent with the traditional taxonomic systems of the genus. Here, we first analysed the complete chloroplast (cp) genomes of seven Filipendula species (including two varieties of F. palmate). Our results indicated that the cp genomes of Filipendula species had few changes in size, ranging from 154,205 bp to 154,633 bp and the average of 36.63% GC content. A total of 126 annotated genes had the identical order and orientation, implying that the cp genome structure of Filipendula species was rather conserved. However, the cp genomes of Filipendula species exhibited structural differences, including gene loss, transposition and inversion when compared to those of other genera of Rosoideae. Moreover, SSRs with the different number were observed in the cp genome of each Filipendula species and sequence divergence mainly occurred in noncoding regions, in which four mutational hotspots were identified. In contrast, only two positive selection genes (matK and rps8) were found. Phylogenetic and molecular-dating analysis indicated that Filipendula species were divergent from other genera of Rosoideae at about 82.88 Ma. Additionally, Filipendula species from East Asia were split at about 9.64 Ma into two major clades. These results provide a basis for further studying the infrageneric classification of Filipendula.

Comparative analyses of Linderniaceae plastomes, with implications for its phylogeny and evolution

Introduction

The recently established Linderniaceae, separated from the traditionally defined Scrophulariaceae, is a taxonomically complicated family. Although previous phylogenetic studies based on a few short DNA markers have made great contributions to the taxonomy of Linderniaceae, limited sampling and low resolution of the phylogenetic tree have failed to resolve controversies between some generic circumscriptions. The plastid genome exhibits a powerful ability to solve phylogenetic relationships ranging from shallow to deep taxonomic levels. To date, no plastid phylogenomic studies have been carried out in Linderniaceae.

Methods

In this study, we newly sequenced 26 plastid genomes of Linderniaceae, including eight genera and 25 species, to explore the phylogenetic relationships and genome evolution of the family through plastid phylogenomic and comparative genomic analyses.

Results

The plastid genome size of Linderniaceae ranged from 152,386 bp to 154,402 bp, exhibiting a typical quartile structure. All plastomes encoded 114 unique genes, comprising 80 protein-coding genes, 30 tRNA genes, and four rRNA genes. The inverted repeat regions were more conserved compared with the single-copy regions. A total of 1803 microsatellites and 1909 long sequence repeats were identified, and five hypervariable regions (petN-psbM, rps16-trnQ, rpl32-trnL, rpl32, and ycf1) were screened out. Most protein-coding genes were relatively conserved, with only the ycf2 gene found under positive selection in a few species. Phylogenomic analyses confirmed that Linderniaceae was a distinctive lineage and revealed that the presently circumscribed Vandellia and Torenia were non-monophyletic.

Discussion

Comparative analyses showed the Linderniaceae plastomes were highly conservative in terms of structure, gene order, and gene content. Combining morphological and molecular evidence, we supported the newly established Yamazakia separating from Vandellia and the monotypic Picria as a separate genus. These findings provide further evidence to recognize the phylogenetic relationships among Linderniaceae and new insights into the evolution of the plastid genomes.

Advancing herbal medicine: enhancing product quality and safety through robust quality control practices

This manuscript provides an in-depth review of the significance of quality control in herbal medication products, focusing on its role in maintaining efficiency and safety. With a historical foundation in traditional medicine systems, herbal remedies have gained widespread popularity as natural alternatives to conventional treatments. However, the increasing demand for these products necessitates stringent quality control measures to ensure consistency and safety. This comprehensive review explores the importance of quality control methods in monitoring various aspects of herbal product development, manufacturing, and distribution. Emphasizing the need for standardized processes, the manuscript delves into the detection and prevention of contaminants, the authentication of herbal ingredients, and the adherence to regulatory standards. Additionally, it highlights the integration of traditional knowledge and modern scientific approaches in achieving optimal quality control outcomes. By emphasizing the role of quality control in herbal medicine, this manuscript contributes to promoting consumer trust, safeguarding public health, and fostering the responsible use of herbal medication products.

Plastid phylogenomics and species discrimination in the "Chinese" clade of Roscoea (Zingiberaceae)

Roscoea is an alpine or subalpine genus from the pan-tropical family Zingiberaceae, which consists of two disjunct groups in geography, namely the "Chinese" clade and the "Himalayan" clade. Despite extensive research on the genus, Roscoea species remain poorly defined and relationships between these species are not well resolved. In this study, we used plastid genomes of nine species and one variety to resolve phylogenetic relationships within the "Chinese" clade of Roscoea and as DNA super barcodes for species discrimination. We found that Roscoea plastid genomes ranged in length from 163,063 to 163,796 bp, and encoded 113 genes, including 79 protein-coding genes, 30 tRNA genes, four rRNA genes. In addition, expansion and contraction of the IR regions showed obvious infraspecific conservatism and interspecific differentiation. Plastid phylogenomics revealed that species belonging to the "Chinese" clade of Roscoea can be divided into four distinct subclades. Furthermore, our analysis supported the independence of R. cautleoides var. pubescens, the recovery of Roscoea pubescens Z.Y. Zhu, and a close relationship between R. humeana and R. cautloides. When we used the plastid genome as a super barcode, we found that it possessed strong discriminatory power (90%) with high support values. Intergenic regions provided similar resolution, which was much better than that of protein-coding regions, hypervariable regions, and DNA universal barcodes. However, plastid genomes could not completely resolve Roscoea phylogeny or definitively discriminate species. These limitations are likely related to the complex history of Roscoea speciation, poorly defined species within the genus, and the maternal inheritance of plastid genomes.

Comparative and phylogenetic analysis of complete chloroplast genome sequences of Salvia regarding its worldwide distribution

Salvia is widely used as medicine, food, and ornamental plants all over the world, with three main distribution centers, the Central and western Asia/Mediterranean (CAM), the East Aisa (EA), and the Central and South America (CASA). Along with its large number of species and world-wide distribution, Salvia is paraphyletic with multiple diversity. Chloroplast genomes (CPs) are useful tools for analyzing the phylogeny of plants at lower taxonomic levels. In this study, we reported chloroplast genomes of five species of Salvia and performed phylogenetic analysis with current available CPs of Salvia. Repeated sequence analysis and comparative analysis of Salvia CPs were also performed with representative species from different distribution centers. The results showed that the genetic characters of the CPs are related to the geographic distribution of plants. Species from CAM diverged first to form a separate group, followed by species from EA, and finally species from CASA. Larger variations of CPs were observed in species from CAM, whereas more deficient sequences and less repeated sequences in the CPs were observed in species from CASA. These results provide valuable information on the development and utilization of the worldwide genetic resources of Salvia.

Phylogeny of Leontopodium (Asteraceae) in China-with a reference to plastid genome and nuclear ribosomal DNA

The infrageneric taxonomy system, species delimitation, and interspecies systematic relationships of Leontopodium remain controversial and complex. However, only a few studies have focused on the molecular phylogeny of this genus. In this study, the characteristics of 43 chloroplast genomes of Leontopodium and its closely related genera were analyzed. Phylogenetic relationships were inferred based on chloroplast genomes and nuclear ribosomal DNA (nrDNA). Finally, together with the morphological characteristics, the relationships within Leontopodium were identified and discussed. The results showed that the chloroplast genomes of Filago, Gamochaeta, and Leontopodium were well-conserved in terms of gene number, gene order, and GC content. The most remarkable differences among the three genera were the length of the complete chloroplast genome, large single-copy region, small single-copy region, and inverted repeat region. In addition, the chloroplast genome structure of Leontopodium exhibited high consistency and was obviously different from that of Filago and Gamochaeta in some regions, such as matk, trnK (UUU)-rps16, petN-psbM, and trnE (UUC)-rpoB. All the phylogenetic trees indicated that Leontopodium was monophyletic. Except for the subgeneric level, our molecular phylogenetic results were inconsistent with the previous taxonomic system, which was based on morphological characteristics. Nevertheless, we found that the characteristics of the leaf base, stem types, and carpopodium base were phylogenetically correlated and may have potential value in the taxonomic study of Leontopodium. In the phylogenetic trees inferred using complete chloroplast genomes, the subgen. Leontopodium was divided into two clades (Clades 1 and 2), with most species in Clade 1 having herbaceous stems, amplexicaul, or sheathed leaves, and constricted carpopodium; most species in Clade 2 had woody stems, not amplexicaul and sheathed leaves, and not constricted carpopodium.

Validation of the shotgun metabarcoding approach for comprehensively identifying herbal products containing plant, fungal, and animal ingredients

Identifying plant, fungal, and animal ingredients in a specific mixture remains challenging during the limitation of PCR amplification and low specificity of traditional methods. Genomic DNA was extracted from mock and pharmaceutical samples. Four type of DNA barcodes were generated from shotgun sequencing dataset with the help of a local bioinformatic pipeline. Taxa of each barcode was assigned by blast to TCM-BOL, BOLD, and GenBank. Traditional methods including microscopy, thin layer chromatography (TLC), and high-performance liquid chromatography (HPLC) were carried out according to Chinese pharmacopoeia. On average, 6.8 Gb shotgun reads were sequenced from genomic DNA of each sample. Then, 97, 11, 10, 14, and one operational taxonomic unit (OTU) were generated for ITS2, psbA-trnH, rbcL, matK, and COI, respectively. All the labeled ingredients including eight plant, one fungal, and one animal species were successfully detected in both the mock and pharmaceutical samples, in which Chebulae Fructus, Poria, and Fritilariae Thunbergia Bulbus were identified via mapping reads to organelle genomes. In addition, four unlabeled plant species were detected from pharmaceutical samples, while 30 genera of fungi, such as Schwanniomyces, Diaporthe, Fusarium were detected from mock and pharmaceutical samples. Furthermore, the microscopic, TLC, and HPLC analysis were all in accordance with the standards stipulated by Chinese Pharmacopoeia. This study indicated that shotgun metabarcoding could simultaneously identified plant, fungal, and animal ingredients in herbal products, which has the ability to serve as a valuable complement to traditional methods.

Phylogeny of Trigonotis in China-with a special reference to its nutlet morphology and plastid genome

The genus Trigonotis comprises nearly 60 species mainly distributed in East and Southeast Asia. China has the largest number of Trigonotis species in the world, with a total of 44 species, of which 38 are endemic. Nutlet morphology is useful for the taxonomic delimitation of Trigonotis. However, there are still controversial circumscriptions of nutlet shape in some species. In previous studies, interspecies phylogenetic relationships were inferred using few DNA markers and very few taxa, which possibly led to erroneous or incomplete conclusions. In this study, the nutlet morphology of 39 Trigonotis taxa and the characteristics of 34 complete chloroplast genomes (29 taxa) were investigated and analyzed. Then, the phylogenetic relationships were discussed within this genus based on complete chloroplast genomes. To the best of our knowledge, this study is the first comprehensive analysis of nutlet morphology and complete chloroplast genome of Trigonotis. Based on nutlet morphology, Trigonotis can be divided into two groups: Group 1, hemispherical or oblique tetrahedron with carpopodiums, and Group 2, inverted tetrahedron without carpopodiums. The chloroplast genome of Trigonotis exhibited a typical quadripartite structure, including 84-86 protein-coding, 37 transfer RNA, and 8 ribosomal RNA genes, with a total length of 147,247-148,986 bp. Genes in the junctions were well conserved in Trigonotis, similar to those in other Boraginaceae s.str. species. Furthermore, Trigonotis chloroplast genomes showed relatively high diversity, with more conserved genic regions than intergenic regions; in addition, we detected 14 hot spots (Pi > 0.005) in non-coding regions. Phylogenetic analyses based on chloroplast genome data identified highly resolved relationships between Trigonotis species. Specifically, Trigonotis was divided into two clades with strong support: one clade included species with hemispherical or oblique tetrahedron nutlets with carpopodiums and bracts, whereas the other clade included species with inverted tetrahedron nutlets without carpopodiums or bracts. Our results may inform future taxonomic, phylogenetic, and evolutionary studies on Boraginaceae.

Comparison of Biological and Genetic Characteristics between Two Most Common Broad-Leaved Weeds in Paddy Fields: Ammannia arenaria and A. multiflora (Lythraceae)

Ammannia arenaria and A. multifloras, morphologically similar at the seedling stage, are the most common broad-leaved weeds in paddy fields. Our study showed that A. arenaria occupied more space than A. multifloras when competing with rice. However, A. multifloras germination has lower temperature adaptability. No difference in sensitivity to common herbicides between two Ammannia species was observed. Chloroplast (cp) genomes could be conducive to clarify their genetic relationship. The complete cp genome sequences of A. arenaria (158,401 bp) and A. multiflora (157,900 bp) were assembled for the first time. In A. arenaria, there were 91 simple sequence repeats, 115 long repeats, and 86 protein-encoding genes, one, sixteen, and thirty more than those in A. multiflora. Inverted repeats regions expansion and contraction and the phylogenetic tree based on cp genomes demonstrated the closely relationship between the two species. However, in A. arenaria, 20 single nucleotide polymorphisms in the CDS region were detected compared to A. multiflora, which can be used to distinguish the two species. Moreover, there was one unique gene, infA, only in A. arenaria. This study provides reliable molecular resources for future research focusing on the infrageneric taxa identification, phylogenetic resolution, population structure, and biodiversity of Ammannia species.

Genome skimming as an efficient tool for authenticating commercial products of the pharmaceutically important Paris yunnanensis (Melanthiaceae)

BACKGROUND

Paris yunnanensis (Melanthiaceae) is a traditional Chinese medicinal plant of significant pharmaceutical importance. Due to previous taxonomic confusion, a congeneric species, Paris liiana, has been mistaken for P. yunnanensis and cultivated on a large scale, leading to the mixing of commercial products (i.e., seedlings and processed rhizomes) of P. yunnanensis with those of P. liiana. This may have adverse effects on quality control in the standardization of P. yunnanensis productions. As the lack of PCR amplifiable genomic DNA within processed rhizomes is an intractable obstacle to the authentication of P. yunnanensis products using PCR-based diagnostic tools, this study aimed to develop a PCR-free method to authenticate commercial P. yunnanensis products, by applying genome skimming to generate complete plastomes and nrDNA arrays for use as the molecular tags.

RESULTS

Based on a dense intraspecies sampling of P. liiana and P. yunnanensis, the robustness of the proposed authentication systems was evaluated by phylogenetic inferences and experimental authentication of commercial seedling and processed rhizome samples. The results indicate that the genetic criteria of both complete plastomes and nrDNA arrays were consistent with the species boundaries to achieve accurate discrimination of P. yunnanensis and P. liinna. Owing to its desirable accuracy and sensitivity, genome skimming can serve as an effective and sensitive tool for monitoring and controlling the trade of P. yunnanensis products.

CONCLUSION

This study provides a new way to solve the long-standing problem of the molecular authentication of processed plant products due to the lack of PCR amplifiable genomic DNA. The proposed authentication system will support quality control in the standardization of P. yunnanensis products in cultivation and drug production. This study also provides molecular evidence to clarify the long-standing taxonomic confusion regarding the species delimitation of P. yunnanensis, which will contribute to the rational exploration and conservation of the species.

Comparitive Analysis of the Chloroplast Genomes of Three Houpoea Plants

The genus Houpoea belongs to the family Magnoliaceae, and the species in this genus have important medicinal values. However, the investigation of the correlation between the evolution of the genus and its phylogeny has been severely hampered by the unknown range of species within the genus and the paucity of research on its chloroplast genome. Thus, we selected three species of Houpoea: Houpoea officinalis var officinalis (OO), Houpoea officinalis var. biloba (OB), and Houpoea rostrata (R). With lengths of 160,153 bp (OO), 160,011 bp (OB), and 160,070 bp (R), respectively, the whole chloroplast genomes (CPGs) of these three Houpoea plants were acquired via Illumina sequencing technology, and the findings were annotated and evaluated. These three chloroplast genomes were revealed by the annotation findings to be typical tetrads. A total of 131, 132, and 120 different genes were annotated. The CPGs of the three species had 52, 47, and 56 repeat sequences, which were primarily found in the ycf2 gene. A useful tool for identifying species is the approximately 170 simple sequence repeats (SSRs) that have been found. The border area of the reverse repetition region (IR) was studied, and it was shown that across the three Houpoea plants, it is highly conservative, with only changes between H. rostrata and the other two plants observed. Numerous highly variable areas (rps3-rps19, rpl32-trnL, ycf1, ccsA, etc.) have the potential to serve as the barcode label for Houpoea, according to an examination of mVISTA and nucleotide diversity (Pi). Phylogenetic relation indicates that Houpoea is a monophyletic taxon, and its genus range and systematic position are consistent with the Magnoliaceae system of Sima Yongkang-Lu Shugang, including five species and varieties of H. officinalis var. officinalis, H. rostrata, H. officinalis var. biloba, Houpoea obovate, and Houpoea tripetala, which evolved and differentiated from the ancestors of Houpoea to the present Houpoea in the above order. This study provides valuable information on the genus Houpoea, enriches the CPG information on Houpoea genus, and provides genetic resources for the further classification of and phylogenetic research on Houpoea.

Uncovering the first complete chloroplast genomics, comparative analysis, and phylogenetic relationships of the medicinal plants Rhamnus cathartica and Frangula alnus (Rhamnaceae)

Rhamnus cathartica and Frangula alnus are economically valuable medicinal plants from the Rhamnaceae family. However, their chloroplast genome structure, phylogenetic position, relationships, and evolution remain poorly understood. Herein, the complete chloroplast genome resources of R. cathartica and F. alnus have been added. The first comparative analysis of the Rhamnus and Frangula species based on complete chloroplast genomes was provided. The chloroplast genomes of R. cathartica and F. alnus exhibited a quadripartite structure, with total lengths of 161,149 bp and 161,255 bp, respectively. The lack of the infA and psbL genes does not negatively impact the normal functioning of Rhamnus and Frangula species. The rpl20 and rpl33 genes are undergoing rapid evolution. Rhamnus and Frangula species prefer amino acids with A/U-terminal codons. There were between 100 and 126 simple sequence repeats and between 38 and 100 long repeats. Several highly divergent intergenic regions (trnK-UUU-trnQ-UUG, atpH-atpI, trnY-GUA-trnE-UUC, trnG-GCC-trnfM-CAU, trnT-UGU-trnF-GAA, rpl20-rps12, and rpl22-rps19) and highly divergent genes (ycf3, ndhA, rpl32, and ycf1) were identified, which could serve as potential phylogenetic markers due to their variability. We reconstructed the phylogenetic relationships among Rhamnus species and F. alnus using complete chloroplast genomes. There is no significant correlation between the medicinal value of the species analyzed and their phylogenetic relationships. These results provide valuable insights for understanding the phylogenetic relationship and evolution of Rhamnus and Frangula species. These findings could serve as a foundation for future studies on the Rhamnaceae.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01331-7.