Search

The complete mitochondrial genome of Alboglossiphonia lata (basionym: Glossiphonia lata), sourced from a biodiversity hotspot of China, has been determined and reported in this study. It was 15,236 bp in length and consisted of 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes and three control regions. The mitogenome was deposited GenBank under the accession number PP165800. A. lata and other species within the Glossiphoniidae family were clustered together with high bootstrap values. The mitochondrial genome of A. lata provides valuable molecular data for further phylogenetic research on the Glossiphoniidae family.

Rhodiola L. is a genus exhibiting rapid radiation and represents a typical case for studying plastid gene adaptation in species that spread from high altitudes to low altitudes. In this study, 23 samples of 18 Rhodiola species were collected from the Qinghai-Tibetan Plateau and five scattered alpine areas, and the plastid genomes (plastomes) of these species were sequenced, annotated, and compared between high-altitude and widely distributed groups. The plastomes of Rhodiola were found to be highly conserved in terms of gene size, content, and order but highly variable in several lineage-specific features, such as codon usage bias, IR boundary shifting, and distinct repeat sequence structures binding to SSRs. Codon usage in the genes of photosystem II exhibited an obvious preference, reflecting significant environmental adaptation pressures. In this study, three repeat regions compounded with trinucleotide and mononucleotide repeats were found for the first time in R. forrestii, R. himalensis, and R. yunnanensis. High-variability regions such as ndhF, ycf1, trnH-psbA, and rpoC1-rpoB were screened, laying the foundation for the precise identification of these species. The phylogenetic analysis revealed the occurrence of cyto-nuclear discordance, likely originating from the frequent interspecific hybridization events observed within Rhodiola species during rapid radiation. Dioecious and hermaphrodite species can be broadly categorized into two subclades, probably they have different environmental adaptation strategies in response to climate change. In addition, the phylogenetic tree supported the monophyly of R. forrestii and R. yunnanensis, which compose R. Sect. Pseudorhodiola. In conclusion, plastome data enrich the genetic information available for the Rhodiola genus and may provide insight into species migration events during climate change.

Accurate identification the species composition in mixtures poses a significant challenge, especially in processed mixtures comprising multiple species, such as those found in food and pharmaceuticals. Therefore, we have attempted to utilize shotgun metabarcoding technology to tackle this issue. In this study, the method was initially established using two mock samples of the Mongolian compound preparation Gurigumu-7 (G-7), which was then applied to three pharmaceutical products and 12 hospital-made preparations. A total of 119.72 Gb of raw data sets were obtained through shotgun metagenomic sequencing. By combining ITS2, matK, and rbcL, all the labeled bio-ingredients specified in the G-7 prescription can be detected, although some species may not be detectable in all samples. The prevalent substitution of Akebia quinata can be found in all the pharmaceutical and hospital samples, except for YN02 and YN12. The toxic alternative to Akebia quinata, Aristolochia manshuriensis, was exclusively identified in the YN02 sample. To further confirm this result, we validated it in YN02 using HPLC and real-time PCR with TaqMan probes. The results showed that aristolochic acid A (AAA) was detected in YN02 using HPLC, and the ITS2 sequence of Aristolochia manshuriensis has been validated in YN02 through qPCR and the use of a TaqMan probe. This study confirms that shotgun metabarcoding can effectively identify the biological components in Mongolian medicine compound preparation G-7. It also demonstrates the method’s potential to be utilized as a general identification technique for mixtures containing a variety of plants.

Abstract The freshwater leech Whitmania pigra (W. pigra) Whitman (Annelida phylum) is a model organism for neurodevelopmental studies. However, molecular biology research on its embryonic development is still scarce. Here, we described a series of developmental stages of the W. pigra embryos and defined five broad stages of embryogenesis: cleavage stages, blastocyst stage, gastrula stage, organogenesis and refinement, juvenile. We obtained a total of 239.64 Gb transcriptome data of eight representative developmental phases of embryos (from blastocyst stage to maturity), which was then assembled into 21,482 unigenes according to our reference genome sequenced by single-molecule real-time (SMRT) long-read sequencing. We found 3114 genes differentially expressed during the eight phases with phase-specific expression pattern. Using a comprehensive transcriptome dataset, we demonstrated that 57, 49 and 77 DEGs were respectively related to morphogenesis, signal pathways and neurogenesis. 49 DEGs related to signal pathways included 30 wnt genes, 14 notch genes, and 5 hedgehog genes. In particular, we found a cluster consisting of 7 genes related to signal pathways as well as synapses, which were essential for regulating embryonic development. Eight genes cooperatively participated in regulating neurogenesis. Our results reveal the whole picture of W. pigra development mechanism from the perspective of transcriptome and provide new clues for organogenesis and neurodevelopmental studies of Annelida species.

Pseudostellaria davidii (Franch.) Pax belongs to subseries distancs of Pseudostellaria (Caryophyllaceae), and is mainly distributed in north-eastern Asia. The complete chloroplast (cp) genome of P. davidii was assembled and annotated for the first time in this study. The cp genome of P. davidii is 149,732 bp in length with the GC content of 36.57%, and it consists of four subregions: a large single-copy (LSC) region of 81,156 bp, a small single-copy (SSC) region of 16,894 bp and two inverted repeats (IR) regions of 25,841 bp each. The cp genome of P. davidii encodes a total of 111 unique genes, which are 77 protein-coding genes, four rRNA genes, and 30 tRNA genes. The results of phylogenetic analysis strongly suggested that Pseudostellaria was a monophyletic group and P. davidii forms an independent sister clade to other species of Pseudostellaria.

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.

Atractylodes species are widely distributed across East Asia and are cultivated as medicinal herbs in China, Japan, and Korea. Their unclear morphological characteristics and low levels of genetic divergence obscure the taxonomic relationships among these species. In this study, 24 plant samples were collected representing five species of Atractylodes located in China; of these, 23 belonged to members of the A. lancea complex. High-throughput sequencing was used to obtain the concatenated nrDNA sequences (18S-ITS1-5.8S-ITS2-28S) and plastid genomes. The concatenated nrDNA sequence lengths for all the Atractylodes species were 5,849 bp, and the GC content was 55%. The lengths of the whole plastid genome sequences ranged from 152,138 bp (A. chinensis) to 153,268 bp (A. lancea), while their insertion/deletion sites were mainly distributed in the intergenic regions. Furthermore, 33, 34, 36, 31, and 32 tandem repeat sequences, as well as 30, 30, 29, 30, and 30 SSR loci, were detected in A. chinensis, A. koreana, A. lancea, A. japonica, and A. macrocephala, respectively. In addition to these findings, a considerable number of heteroplasmic variations were detected in the plastid genomes, implying a complicated phylogenetic history for Atractylodes. The results of the phylogenetic analysis involving concatenated nrDNA sequences showed that A. lancea and A. japonica formed two separate clades, with A. chinensis and A. koreana constituting their sister clade, while A. lancea, A. koreana, A. chinensis, and A. japonica were found based on plastid datasets to represent a mixed clade on the phylogenetic tree. Phylogenetic network analysis suggested that A. lancea may have hybridized with the common ancestor of A. chinensis and A. japonica, while ABBA–BABA tests of SNPs in the plastid genomes showed that A. chinensis was more closely related to A. japonica than to A. lancea. This study reveals the extensive discordance and complexity of the relationships across the members of the A. lancea complex (A. lancea, A. chinensis, A. koreana, and A. japonica) according to cytonuclear genomic data; this may be caused by interspecific hybridization or gene introgression.

Isodon rubescens (Hemsley) H. Hara is the source of Donglingcao under the monograph Rabdosiae Rubescentis Herba in Chinese Pharmacopoeia. In the local marketplace, this medicine can be accidentally contaminated, deliberately substituted, or mixed with other related species. The contaminants of herbal products are a threat to consumer safety. Due to the scarcity of genetic information on Isodon plants, more molecular markers are needed to avoid misidentification. In the present study, the complete chloroplast (cp) genome of seven species of Isodon was sequenced, de novo assembled and characterized. The cp genomes of these species universally exhibited a conserved quadripartite structure, i.e., two inverted repeats (IRs) containing most of the ribosomal RNA genes and two unique regions (large single copy and small single copy). Moreover, the genome structure, codon usage, and repeat sequences were highly conserved and showed similarities among the seven species. Five highly variable regions (trnS-GCU-trnT-CGU, atpH-atpI, trnE-UUC-trnT-GGU, ndhC-trnM-CAU, and rps15-ycf1) might be potential molecular markers for identifying I. rubescens and its contaminants. These findings provide valuable information for further species identification, evolution, and phylogenetic research of Isodon.

Briggsia chienii W. Y. Chun 1946 is an endemic herbaceous perennial species distributed in southern China. In this study, we firstly characterized the complete chloroplast genome sequence of B. chienii and provided new molecular resources for promoting its conservation and taxonomic assignment. Its complete chloroplast genome is 154,082 bp in length and contains the typical quadripartite structure of angiosperm plastome, including two inverted repeat (IR) regions of 25,447 bp, a large single-copy (LSC) region of 85,035 bp, and a small single-copy (SSC) region of 18,153 bp. The plastome contains 114 genes, consisting of 80 protein-coding genes, 30 tRNA gene, and 4 rRNA genes. The overall GC content in the plastome of B. chienii is 37.4%, which is lower than lots of angiosperm plastome. The phylogenetic result indicated that B. chienii exhibited the closest relationship with Oreocharis cotinifolia W. T. Wang 1983, and provided new information for the phylogeny relationship of genus Briggsia.

ABSTRACT This study determined the composition of fungal communities and characterized the enriched fungal species in raw and roasted malts via the third-generation PacBio-based full-length single-molecule real-time (SMRT) sequencing of the full-length amplicon of the internal transcribed spacer (ITS) region. In total, one kingdom, six phyla, 23 classes, 56 orders, 120 families, 188 genera, 333 species, and 780 operational taxonomic units (OTUs) were detected with satisfactory sequencing depth and sample size. Wickerhamomyces (56%), Cyberlindnera (15%), Dipodascus (12%), and Candida (6.1%) were characterized as the dominant genera in the raw malts, and Aspergillus (35%), Dipodascus (21%), Wickerhamomyces (11%), and Candida (3.5%) in the roasted malts. Aspergillus proliferans, Aspergillus penicillioides, and Wickerhamomyces anomalus represented the crucial biomarkers causing intergroup differences. Correlation analysis regarding environmental factors indicated that the water activity (aw) of the samples affected the composition of the fungal communities in the malts. In practice, special attention should be paid to the mycotoxin-producing fungi, as well as other fungal genera that are inversely correlated with their growth, to ensure the safe use of malt and its end products. IMPORTANCE Fungal contamination and secondary metabolite accumulation in agricultural products represent a global food safety challenge. Although high-throughput sequencing (HTS) is beneficial for explaining fungal communities, it presents disadvantages, such as short reads, species-level resolution, and uncertain identification. This work represents the first attempt to characterize the fungal community diversity, with a particular focus on mycotoxin-producing fungi, in malt via the third-generation PacBio-based full-length SMRT sequencing of the ITS region, aiming to explore and compare the differences between the fungal communities of raw and roasted malts. The research is beneficial for developing effective biological control and conservation measures, including improving the roasting conditions, monitoring the environmental humidity and aw, and effectively eliminating and degrading fungi in the industry chain according to the diverse fungal communities determined, for the safe use of malts and their end products, such as beers. In addition, the third-generation SMRT sequencing technology allows highly efficient analysis of fungal community diversity in complex matrices, yielding fast, high-resolution long reads at the species level. It can be extended to different research fields, updating modern molecular methodology and bioinformatics databases.

The complete mitochondrial genome of Haemadipsa tianmushana Song 1977 from China has been determined and reported for the first time in this study. It was 14,625 bp in length and consisted of 22 tRNA genes, 2 rRNA genes, 13 protein-coding genes (PCGs), and 3 control regions. The nucleotide base content of the complete mitogenome for this species was 35.1% A, 10.5% C, 11.6% G, and 42.8% T. The tRNA genes were ranged from 57 bp (SerTCT) to 66 bp (GlnTTG) in length. The phylogenetic analyses indicated that Hirudinea is a mono-phyletic clade. And it includes Whitmania acranulata, Whitmania pigra, Whitmania laevis, Zeylanicobdella arugamensis, Ozobranchus jantseanus and Placobdella lamothei. In Hirudiniformes, H. tianmushana and three species of Haemopidae were obviously clustered into two independent branches. This result is consistent with a taxonomy that they all belong to the same suborder. This study adds to the genetic resources currently available for the species.

A “signal off” fluorescent aptasensor based on graphene oxide (GO) nanosheet and double-stranded DNA structure was fabricated for OTA detection. In the absence of OTA, the aptamer and its complementary DNA (cDNA) formed double-stranded conjugates that could coexist with GO, presenting fluorescence responses. Then, the presented OTA was captured by the aptamers, causing the release of cDNA-FAM probes. The free probes were adsorbed by GO, leading to an OTA concentration-dependent fluorescence quenching through fluorescence resonance energy transfer. Under optimum conditions, the fluorescent aptasensor exhibited outstanding sensitivity with a LOD of 11 pg/mL and a wide dynamic range of 0.04–30 ng/mL. The selectivity of the aptasensor was confirmed against other four mycotoxins, and the feasibility and reliability were verified by determining OTA in the spiked malt samples with satisfactory recovery of 95.50%-112.05%. This aptasensing platform could be adapted to measure other mycotoxins by replacing the aptamer sequences for food safety.

Affected by some external conditions and internal factors, pesticides can be transferred from tea into its infusion, causing subsequent damage to humans as tea infusion is generally consumed. This study aimed to explore the inherent regularity in transfer behavior of 23 pesticides belonging to different classes from honeysuckle to its tea infusion, and to understand the effects of external brewing conditions and internal physicochemical parameters of the pesticides on their transfer rates. Results indicated that the transfer rates (Rt) of pesticides from honeysuckle into tea solutions increased with prolonged brewing time, or adding a cover on a container, but decreased with increasing the times of infusion. In addition, the transfer potential of these pesticides greatly depended on their physicochemical properties but not their type. The pesticides with high water solubility and low water partition coefficient (LogKow, e.g., omethoate) were more easily transferred into tea infusions than those with low water solubility and high LogKow (e.g., chlorpyrifos). Compared the tea brewing in a covered container, the empirical models obtained in an uncovered cup predicted the transfer behavior and drinking risk of pesticides potentially introduced into honeysuckle and its tea infusion. The linear equation was as follow: Rt = 10.756 LogWS + 7.517, R = 0.8771. In practice, honeysuckle should be brewed in an uncovered cup within a short brewing time, and the first tea infusion should be abandoned to reduce the transfer percentage of pesticides. This study provided beneficial references for pesticide application in honeysuckle plantation to establish realistic maximum residue limits of multi-pesticides in honeysuckle tea and related products.

Bupleuri Radix is the dry root of certain species of the genus Bupleurum and is commonly used in traditional Chinese medicine. The increasing global demand for Bupleuri Radix cannot be fulfilled with wild populations only. Therefore, cultivated Bupleurum is now the main commercial source of this medicinal product. Different species of Bupleurum show different medicinal properties and clinical effects, making reliable authentication and assignment of correct botanical origin for medicinal species critical. However, accurate identification of the cultivated Bupleurum species is difficult due to dramatic morphological variations resulting from cultivation. In this study, we sampled 56 cultivated Bupleurum populations of six different morphotypes (Types A-F) from the main production areas of China, and 10 wild populations of four species were used as reference materials. Conventional DNA barcoding was conducted to identify cultivated Bupleurum species. Additionally, verification based on complete chloroplast genomes was performed and new chloroplast markers were developed and evaluated. The combination of these methods resulted in the successful identification of all cultivated Bupleurum individuals. Three chloroplast regions are recommended as additional barcodes for the genus: ycf4_cemA, psaJ_rpl33, and ndhE_ndhG. This is a reliable and promising strategy that can be applied to the authentication of natural products and the identification of other medicinal plant species with similar taxonomic problems.

Atractylodes japonica Koidz. ex Kitam. is a perennial herbal plant, and its dried rhizomes have been widely used as traditional medicine in China and Japan. In this study, we assembled and annotated the complete chloroplast (cp) genome sequence of A. japonica using the high-throughput sequencing approach. The cp genome of A. japonica is 153,208 bp in length with the overall GC content of 37.7%, including two inverted repeat (IR) regions of 25,147 bp, which was separated by a large single-copy (LSC) region of 84,255 bp and a small single-copy (SSC) region of 18,659 bp. 113 unique genes were annotated in the genome, including 80 protein-coding genes, 29 represented tRNA genes, and four denoted rRNA genes. A maximum-likelihood phylogenetic analysis with 38 complete cp sequences showed that Atractylodes formed a monophyletic clade, and A. japonica and A. koreana formed a subclade in Atractylodes. This study provides the chloroplast genome structure features and phylogenetic relationship of A. japonica.

Atractylodes koreana (Nakai) Kitam is a perennial herb of Asteraceae, mainly distributed in China and Korea, which is the main adulterant of traditional herbal medicine ‘Cangzhu’. In the present study, we reported the complete chloroplast (cp) genome of A. koreana with the total length of 153,232 bp, which is consisted of four regions, including one large single copy (LSC) region of 84,250 bp, one small single copy (SSC) region of 18,690 bp, and two inverted repeat regions (IRa and IRb) of 25,146 bp. The GC content of the complete cp genome is 37.7%. A total of 110 unique genes were annotated, comprising 79 protein-coding genes, 27 transfer RNA (tRNA) genes and four ribosome RNA (rRNA) genes. Moreover, nine protein-coding genes contained one intron and three protein-coding genes (clpP, ycf3, and rps12) contained two introns. The phylogenetic analysis indicated that A. koreana is a sister group of A. chinensis and A. lancea.

Bupleurum yinchowense Shan & Yin Li was first described as a new Bupleurum species in 1974, but its classification status has always been disputed. Here, its complete chloroplast genome was provided to resolve this issue. The length of the B. yinchowense chloroplast genome is 155,851 bp and composed of two inverted repeats (IR: 26,307 bp), a large single-copy region (LSC: 85,625 bp), and a small single-copy region (SSC: 17,612 bp). The overall GC content is 37.6%. The chloroplast genome consists of 113 genes, including 79 protein-coding genes, four rRNA genes, and 30 tRNA genes. Phylogenetic analysis suggested that Bupleurum yinchowense holds a distinct phylogenetic position and can be considered as an accepted species.

The root of Bupleurum marginatum var. stenophyllum (H. Wolff) Shan & Yin Li (Apiaceae), a new substitution for the popular Chinese medicinal material, Bupleuri Radix (Chai hu), is not easily distinguishable via traditional methods. The complete chloroplast genome sequence of B. marginatum var. stenophyllum was characterized using next-generation sequencing and the de novo assembly method. The complete genome was 155,576 bp in length and contained two inverted repeat (IR) regions of 26,311 bp, a large single-copy (LSC) region of 85,351 bp, and a small single-copy (SSC) region of 17,603 bp. It encoded 113 unique genes consisting of 79 protein-coding genes (PCGs), 30 transfer RNA genes, and four ribosomal RNA genes. Importantly, three genes (petB, petD and rps16) with small exon, and one trans-splicing gene (rps12) were correctly annotated. The overall GC content of the B. marginatum var. stenophyllum chloroplast genome is 37.7%. The phylogenetic analyses indicated that B. marginatum var. stenophyllum was closely related to B. marginatum. Moreover, many genetic information sites were available for distinguishing B. marginatum var. stenophyllum from the official ‘Chai hu’ plant sources, B. scorzonerifolium Willd. and B. chinense DC.

With the widespread use of traditional medicine around the world, the safety and efficacy of traditional herbal patent medicine have become an increasing concern to the public. However, it is difficult to supervise the authenticity of herbal materials in mixed herbal products according to the current quality standards, especially for traditional herbal patent medicine, with a distinct variance in the dosage of herbal materials. This study utilized the shotgun metabarcoding approach to analyze the biological ingredients of Fuke Desheng Wan (FKDSW), which is an effective traditional herbal product for the treatment of dysmenorrhea. Six herbal materials were collected, and a lab-made mock FKDSW sample was produced to establish a method for the authentication assessment of biological ingredients in traditional herbal patent medicine based on shotgun metabarcoding. Furthermore, four commercial FKDSW samples were collected to verify the practicality of the shotgun metabarcoding approach. Then, a total of 52.16 Gb raw data for 174 million paired-end reads was generated using the Illumina NovaSeq sequencing platform. Meanwhile, 228, 23, and 14 operational taxonomic units (OTUs) were obtained for the ITS2, matK, and rbcL regions, respectively, after bioinformatic analysis. Moreover, no differences were evident between the assembly sequences obtained via shotgun metabarcoding and their corresponding reference sequences of the same species obtained via Sanger sequencing, except for part of the ITS2 and matK assembly sequences of Paeonia lactiflora Pall., Saussurea costus (Falc.) Lipsch. and Bupleurum chinense DC. with 1–6 different bases. The identification results showed that all six prescribed ingredients were successfully detected and that the non-authentic ingredient of Bupleuri Radix (Chaihu, Bupleurum chinense DC. or Bupleurum scorzonerifolium Willd.) was found in all the commercial samples, namely Bupleurum falcatum L. Here, 25 weed species representing 16 genera of ten families were detected. Moreover, 26 fungal genera belonging to 17 families were found in both lab-made and commercial FKDSW samples. This study demonstrated that the shotgun metabarcoding approach could overcome the biased PCR amplification and authenticate the biological ingredients of traditional herbal patent medicine with a distinct variance in the dosage of the herbal materials. Therefore, this provides an appropriate evaluation method for improving the safety and efficacy of traditional herbal patent medicine.

As one of the high-incidence diseases in the world, pharyngitis seriously affects the lives of those with the condition. Qingguo Wan is a herbal medicine used for treating pharyngitis, and its quality evaluation is currently only accomplished via traditional identification. However, precise identification becomes challenging with fake products on the market or fungal contamination during the production process. This study used the Illumina NovaSeq platform for targeting the ITS2, psbA-trnH, matK, and rbcL sequences to survey the species composition of lab-made and commercial samples. The results showed that a total of 34.56 Gb of raw data that was obtained represented more than 0.23 billion reads. After assembly, annotation, and operational taxonomic unit clustering, 103, 12, 10, and 12 OTUs were obtained, which belonged to the ITS2, psbA-trnH, matK, and rbcL sequences of the mock lab-made and commercial samples. The analytical results indicated that the sequences of all the prescription ingredients were successfully obtained in the two lab-made samples. The positive control medicinal Panax quinquefolius L. sequence was obtained in HSZY175, while Scutellaria baicalensis Georgi, Lonicera japonica Thunb. Menispermum dauricum DC. and Paeonia lactiflora Pall. were detected in the three commercial samples. The detection results of the other four herbs in different fragments were not all the same. In addition, a total of 28 fungi OTUs, representing 19 families and 20 genera, were obtained from both the commercial and mock lab-made samples. Aspergillus, Cladosporium, and Penicillium dominated among the 20 genera. This study demonstrated that the shotgun metabarcoding method is a powerful tool for the molecular identification of the biological ingredients in Qingguo Wan. It can be used to effectively supplement traditional methods while providing a new technique for the quality evaluation of Qingguo Wan.

Traditional herbal patent medicine typically consists of multiple ingredients, making it challenging to supervise contamination by impurities and the improper use of raw materials. This study employed shotgun metabarcoding for the species identification of biological ingredients in traditional herbal patent medicine, Wuhu San. The five prescribed herbal materials found in Wuhu San were collected, and their reference sequences were obtained by traditional DNA barcoding using Sanger sequencing. Two lab-made and three commercial Wuhu San samples were collected, and a total of 37.14 Gb of shotgun sequencing data was obtained for these five samples using the Illumina sequencing platform. A total of 1,421,013 paired-end reads were enriched for the Internal Transcribed Spacer 2 (ITS2), psbA and trnH intergenic spacer region (psbA-trnH), maturase k (matK), and ribulose-1, 5-bisphosphate carboxylase (rbcL) regions. Furthermore, 80, 11, 9, and 8 operational taxonomic units were obtained for the ITS2, psbA-trnH, matK, and rbcL regions, respectively, after metagenomic assembly, annotation, and chimeric detection. In the two lab-made mock samples, all labeled ingredients in the Wuhu San prescription were successfully detected, and the positive control, Panax quinquefolius L., was detected in the HSZY172 mock sample. Three species, namely Angelica sinensis (Oliv.) Diels, Saposhnikovia divaricata (Turcz. ex Ledeb.) Schischk., and Carthamus tinctorius L., belonging to three labeled ingredients, Angelicae Sinensis Radix (Danggui), Saposhnikoviae Radix (Fangfeng), and Carthami Flos (Honghua), were detected in the three commercial samples. Angelica dahurica (Hoffm.) Benth. & Hook. f. ex Franch. & Sav., the original Angelicae Dahuricae Radix (Baizhi) species, was only detected in WHS003. Arisaema erubescens (Wall.) Schott, Arisaema heterophyllum Blume, or Arisaema amurense Maxim., the original Arisaematis Rhizoma (Tiannanxing) species, were not detected in any of the commercial samples, which could be attributed to the fact that this medicinal material underwent extensive processing. In addition, the Saposhnikovia divaricata adulterant was detected in all the commercial samples, while 24 fungal genera, including Aspergillus, were identified in both the lab-made and commercial samples. This study showed that shotgun metabarcoding provided alternative strategy and technical means for identifying prescribed ingredients in traditional herbal patent medicine and displayed the potential to effectively complement traditional methods.

Increased public awareness of nutritional and health issues has resulted in the increasing consumption of food and herbal products made from the root of Pueraria montana var. lobata (Willd.) Maesen & S. M. Almeida ex Sanjappa & Predeep (kudzu vine) and P. montana var. thomsonii (Benth.) M. R. Almeida. The famous herbal medicine Yufeng Ningxin, which is used to treat cardiovascular diseases, can be legally produced only using P. montana var. lobata. However, precise identification at the subspecies level is usually challenging when these products’ ingredients lose their morphological characteristics after deep processing. Here, six herbarium specimens, 21 expert-identified original plant samples, 30 raw material samples, 10 food products and 12 herbal products were collected to test the subspecies-level authentication abilities of ITS2 sequences. The results showed that ITS2 sequences can distinguish P. montana var. lobata from P. montana var. thomsonii with stable single nucleotide polymorphism (SNP) sites. A total of 93.3% of raw material samples were consistent with the markings on their labels, but only 50% of Gegen Powder samples were made from P. montana var. lobata. High-quality ITS2 sequences were successfully obtained from nine of the 12 herbal products using Sanger sequencing. Substitution and fungal contamination were detected in 3 herbal products by further DNA metabarcoding, even though thin-layer chromatography (TLC) and high-performance liquid chromatography (HPLC) tests verified that the products met existing quality standards. This study demonstrated that DNA barcoding is a powerful tool for the identification of P. montana var. lobata and P. montana var. thomsonii at the subspecies level, and we conclude that DNA barcodes can be broadly applied to trace the raw materials of food and herbal products.

Seed is not only the main reproductive organ of most of herbal plants but also an important part of Traditional Chinese Medicine (TCM). Seed TCMs possess important medicinal properties and have been widely used as components of pharmaceutical products. In parallel with the increasing popularity and accessibility of seeds as medicinal products in recent years, numerous substitutes and adulterants have also appeared on the market. Due to the small volume and similar appearances of many seed TCMs, they are very difficult to accurately identify the constituent plant species through organoleptic methods. Usage of the wrong herb may be ineffective or may worsen the condition and even cause death. Correct identification of seed herbal medicines is therefore essential for their safe use. Here, we acquired 177 ITS2 sequences and 15 psbA-trnH sequences from 51 kinds of seed TCMs belonging to 64 species that have been described in the Chinese Pharmacopoeia. Tree-building analysis showed that the ITS2 sequences of 48 seed TCMs can be differentiated from each other, and they formed distinct, non-overlapping groups in the maximum-likelihood tree. Furthermore, all of the sequences acquired in this study have been submitted to the public DNA barcoding system for herbal medicine, and this integrated database was used to identify 400 seed TCM samples purchased from medicinal markets, drug stores, and the Internet, enabling the identification of 7.5% of the samples as containing non-declared species. This study provides a brief operating procedure for the identification of seed TCMs found in herbal medicine. In the future, researchers and traditional herbal medicine enterprises can use this system to test their herbal materials.

Herbal drug authentication is an important task in traditional medicine; however, it is challenged by the limitations of traditional authentication methods and the lack of trained experts. DNA barcoding is conspicuous in almost all areas of the biological sciences and has already been added to the British pharmacopeia and Chinese pharmacopeia for routine herbal drug authentication. However, DNA barcoding for the Korean pharmacopeia still requires significant improvements. Here, we present a DNA barcode reference library for herbal drugs in the Korean pharmacopeia and developed a species identification engine named KP-IDE to facilitate the adoption of this DNA reference library for the herbal drug authentication. Using taxonomy records, specimen records, sequence records, and reference records, KP-IDE can identify an unknown specimen. Currently, there are 6,777 taxonomy records, 1,054 specimen records, 30,744 sequence records (ITS2 and psbA-trnH) and 285 reference records. Moreover, 27 herbal drug materials were collected from the Seoul Yangnyeongsi herbal medicine market to give an example for real herbal drugs authentications. Our study demonstrates the prospects of the DNA barcode reference library for the Korean pharmacopeia and provides future directions for the use of DNA barcoding for authenticating herbal drugs listed in other modern pharmacopeias.

Substandard traditional patent medicines may lead to global safety-related issues. Protecting consumers from the health risks associated with the integrity and authenticity of herbal preparations is of great concern. Of particular concern is quality control for traditional patent medicines. Here, we establish an effective approach for verifying the biological composition of traditional patent medicines based on single-molecule real-time (SMRT) sequencing and DNA barcoding. Yimu Wan (YMW), a classical herbal prescription recorded in the Chinese Pharmacopoeia, was chosen to test the method. Two reference YMW samples were used to establish a standard method for analysis, which was then applied to three different batches of commercial YMW samples. A total of 3703 and 4810 circular-consensus sequencing (CCS) reads from two reference and three commercial YMW samples were mapped to the ITS2 and psbA-trnH regions, respectively. Moreover, comparison of intraspecific genetic distances based on SMRT sequencing data with reference data from Sanger sequencing revealed an ITS2 and psbA-trnH intergenic spacer that exhibited high intraspecific divergence, with the sites of variation showing significant differences within species. Using the CCS strategy for SMRT sequencing analysis was adequate to guarantee the accuracy of identification. This study demonstrates the application of SMRT sequencing to detect the biological ingredients of herbal preparations. SMRT sequencing provides an affordable way to monitor the legality and safety of traditional patent medicines.

Scutellaria baicalensis Georgi has long been used in traditional medicine to treat various such widely varying diseases and has been listed in the Chinese Pharmacopeia, the Japanese Pharmacopeia, the Korean Pharmacopoeia and the European Pharmacopoeia. Flavonoids, especially wogonin, wogonoside, baicalin, and baicalein, are its main functional ingredients with various pharmacological activities. Although pharmaological studies for these flavonoid components have been well conducted, the molecular mechanism of their biosynthesis remains unclear in S. baicalensis. In this study, Illumina/Solexa deep sequencing generated more than 91 million paired-end reads and 49,507 unigenes from S. baicalensis roots, stems, leaves and flowers. More than 70% unigenes were annotated in at least one of the five public databases and 13,627 unigenes were assigned to 3,810 KEGG genes involved in 579 different pathways. 54 unigenes that encode 12 key enzymes involved in the pathway of flavonoid biosynthesis were discovered. One baicalinase and three baicalein 7-O-glucuronosyltransferases genes potentially involved in the transformation between baicalin/wogonoside and baicalein/wogonin were identified. Four candidate 6-hydroxylase genes for the formation of baicalin/baicalein and one candidate 8-O-methyltransferase gene for the biosynthesis of wogonoside/wogonin were also recognized. Our results further support the conclusion that, in S. baicalensis, 3,5,7-trihydroxyflavone was the precursor of the four above compounds. Then, the differential expression models and simple sequence repeats associated with these genes were carefully analyzed. All of these results not only enrich the gene resource but also benefit research into the molecular genetics and functional genomics in S. baicalensis.

The DNA barcoding technology uses a standard region of DNA sequence for species identification and discovery. At present, "DNA barcode" actually refers to DNA sequences, which are not amenable to information storage, recognition, and retrieval. Our aim is to identify the best symbology that can represent DNA barcode sequences in practical applications. A comprehensive set of sequences for five DNA barcode markers ITS2, rbcL, matK, psbA-trnH, and CO1 was used as the test data. Fifty-three different types of one-dimensional and ten two-dimensional barcode symbologies were compared based on different criteria, such as coding capacity, compression efficiency, and error detection ability. The quick response (QR) code was found to have the largest coding capacity and relatively high compression ratio. To facilitate the further usage of QR code-based DNA barcodes, a web server was developed and is accessible at http://qrfordna.dnsalias.org. The web server allows users to retrieve the QR code for a species of interests, convert a DNA sequence to and from a QR code, and perform species identification based on local and global sequence similarities. In summary, the first comprehensive evaluation of various barcode symbologies has been carried out. The QR code has been found to be the most appropriate symbology for DNA barcode sequences. A web server has also been constructed to allow biologists to utilize QR codes in practical DNA barcoding applications.

BackgroundInternal transcribed spacer of nuclear ribosomal DNA (nrDNA) is already one of the most popular phylogenetic and DNA barcoding markers. However, the existence of its multiple copies has complicated such usage and a detailed characterization of intra-genomic variations is critical to address such concerns.Methodology/principal findingsIn this study, we used sequence-tagged pyrosequencing and genome-wide analyses to characterize intra-genomic variations of internal transcribed spacer 2 (ITS2) regions from 178 plant species. We discovered that mutation of ITS2 is frequent, with a mean of 35 variants per species. And on average, three of the most abundant variants make up 91% of all ITS2 copies. Moreover, we found different congeneric species share identical variants in 13 genera. Interestingly, different species across different genera also share identical variants. In particular, one minor variant of ITS2 in Eleutherococcus giraldii was found identical to the ITS2 major variant of Panax ginseng, both from Araliaceae family. In addition, DNA barcoding gap analysis showed that the intra-genomic distances were markedly smaller than those of the intra-specific or inter-specific variants. When each of 5543 variants were examined for its species discrimination efficiency, a 97% success rate was obtained at the species level.ConclusionsIdentification of identical ITS2 variants across intra-generic or inter-generic species revealed complex species evolutionary history, possibly, horizontal gene transfer and ancestral hybridization. Although intra-genomic multiple variants are frequently found within each genome, the usage of the major variants alone is sufficient for phylogeny construction and species determination in most cases. Furthermore, the inclusion of minor variants further improves the resolution of species identification.

BACKGROUND: The trnH-psbA intergenic spacer region has been used in many DNA barcoding studies. However, a comprehensive evaluation with rigorous sequence preprocessing and statistical testing on the utility of trnH-psbA and its combinations as DNA barcodes is lacking. METHODOLOGY/PRINCIPAL FINDINGS: Sequences were searched from GenBank for a meta-analysis on the usefulness of trnH-psbA and its combinations as DNA barcodes. After preprocessing, we constructed full and matching data sets that contained 17 983 trnH-psbA sequences and 2190 sets of trnH-psbA, matK, rbcL, and ITS2 sequences from the same sample, repectively. These datasets were used to analyze the ability of trnH-psbA and its combinations to discriminate species by the BLAST and BLAST+P methods. The Fisher's exact test was used to evaluate the significance of performance differences. For the full data set, the identification success rates of trnH-psbA exceeded 70% in 18 families and 12 genera, respectively. For the matching data set, the identification rates of trnH-psbA were significantly higher than those of the other loci in two families and four genera. Similarly, the identification rates of trnH-psbA+ITS2 were significantly higher than those of matK+rbcL in 18 families and 21 genera. CONCLUSION/SIGNIFICANE: This study provides valuable information on the higher utility of trnH-psbA and its combinations. We found that trnH-psbA+ITS2 combination performs better or equally well compared with other combinations in most taxonomic groups investigated. This information will guide the optimal usage of trnH-psbA and its combinations for species identification.

BackgroundThe plant working group of the Consortium for the Barcode of Life recommended the two-locus combination of rbcL+matK as the plant barcode, yet the combination was shown to successfully discriminate among 907 samples from 550 species at the species level with a probability of 72%. The group admits that the two-locus barcode is far from perfect due to the low identification rate, and the search is not over.Methodology/principal findingsHere, we compared seven candidate DNA barcodes (psbA-trnH, matK, rbcL, rpoC1, ycf5, ITS2, and ITS) from medicinal plant species. Our ranking criteria included PCR amplification efficiency, differential intra- and inter-specific divergences, and the DNA barcoding gap. Our data suggest that the second internal transcribed spacer (ITS2) of nuclear ribosomal DNA represents the most suitable region for DNA barcoding applications. Furthermore, we tested the discrimination ability of ITS2 in more than 6600 plant samples belonging to 4800 species from 753 distinct genera and found that the rate of successful identification with the ITS2 was 92.7% at the species level.ConclusionsThe ITS2 region can be potentially used as a standard DNA barcode to identify medicinal plants and their closely related species. We also propose that ITS2 can serve as a novel universal barcode for the identification of a broader range of plant taxa.