Verification of Thai ethnobotanical medicine “Kamlang Suea Khrong” driven by multiplex PCR and powerful TLC techniques

Kamlang Suea Khrong (KSK) crude drug, a traditional Thai medicine used for oral tonic and analgesic purposes, is obtained from three origins: the inner stem bark of Betula alnoides (BA) or the stems of Strychnos axillaris (SA) or Ziziphus attopensis (ZA). According to the previous reports, SA contains strychnine-type alkaloids that probably cause poisoning; however, only organoleptic approaches are insufficient to differentiate SA from the other plant materials. To ensure the botanical origin of KSK crude drug, powerful and reliable tools are desperately needed. Therefore, molecular and chemical identification methods, DNA barcoding and thin-layer chromatography (TLC), were investigated. Reference databases, i.e., the ITS region and phytochemical profile of the authentic plant species, were conducted. In case of molecular analysis, multiplex polymerase chain reaction (PCR) based on species-specific primers was applied. Regarding species-specific primers designation, the suitability of three candidate barcode regions (ITS, ITS1, and ITS2) was evaluated by genetic distance using K2P model. ITS2 presented the highest interspecific variability was verified its discrimination power by tree topology. Accordingly, ITS2 was used to create primers that successfully specified plant species of authentic samples. For chemical analysis, TLC with toluene:ethyl acetate:ammonia (1:9:0.025) and hierarchical clustering were operated to identify the authentic crude drugs. The developed multiplex PCR and TLC methods were then applied to identify five commercial KSK crude drugs (CK1-CK5). Both methods correspondingly indicated that CK1-CK2 and CK3-CK5 were originated from BA and ZA, respectively. Molecular and chemical approaches are convenient and effective identification methods that can be performed for the routine quality-control of the KSK crude drugs for consumer reliance. According to chemical analysis, the results indicated BA, SA, and ZA have distinct chemical profiles, leading to differences in pharmacological activities. Consequently, further scientific investigations are required to ensure the quality and safety of Thai ethnobotanical medicine known as KSK.


Introduction
In recent decades, there has been a revival of interest in traditional medicine since it is available, accessible, affordable, and acceptable to the local population [1]. However, taking traditional medicines cause efficacy and safety concerns due to poor quality and adulterated or counterfeit products. Accordingly, WHO strategies intend to ensure the quality, safety, proper use, and effectiveness of traditional medicines [2]. In general, medicinal herbs have been traded as crude drugs in the form of dried slices or powders of the parts used, such as the leaves, bark, and roots; therefore, identification using only organoleptic methods is impossible. Moreover, their vernacular names are often homonymic, i.e., local names refer to multiple scientific plant species. For example, Curcuma spp., Rinorea spp., Thunbergia laurifolia, and Crotalaria spectabilis are known in Thailand as "Rang Chuet". However, C. spectabilis should be avoided because it contains poisonous pyrrolizidine alkaloids, which cause hepatotoxicity in mammals [3]. Although several medicinal herbs on the market have the same vernacular name and traditional uses, their phytochemical constituents, pharmacological activities, and toxicities might differ. Hence, reliable plant identification of commercial crude drugs is required to ensure quality and safety. The focus of this study was on "Kamlang Suea Khrong" (KSK), a homonymic Thai medicinal herb that means "the power of tiger".
KSK, a Thai ethnobotanical medicine, refers to at least three medicinal plants, i.e., Betula alnoides Buch.-Ham. ex D.Don (Betulaceae), Strychnos axillaris Colebr. (Strychnaceae), and Ziziphus attopensis Pierre (Rhamnaceae) [4][5][6][7][8]. B. alnoides is a tree with exfoliating aromatic stem bark native to northern Thailand [4][5][6]. S. axillaris is a liana distributed in northeastern Thailand [4]. While, the scandent shrub Z. attopensis is predominantly found in northern and northeastern regions of Thailand [7,8]. Furthermore, these three plants species are mainly located in other countries, including Cambodia, China, Laos, and Myanmar [9]. KSK crude drug, a material intended for medicinal use, is obtained from the inner stem bark of B. alnoides (BA) or the stem of S. axillaris (SA) or the stem of Z. attopensis (ZA) [4][5][6][7][8]. KSK crude drug has been traditionally used as an herbal tonic to strengthen the body in Thailand [10]. It has also been used for centuries by hard-working farmers to boost their energy and relieve their aches and pains [11]. KSK crude drug (BA or SA or ZA) has been combined with other herbs such as stems of Cryptolepis buchananii Roem. & Schult. (Asclepiadaceae) and Salacia chinensis L. (Celastraceae) to prepare oral longevity recipes, including decoction (Ya Tom) and spirits (Ya Dong Sura) [10][11][12]. Ya Tom is generally prepared by boiling plant materials in water (1 :3 v/v) at medium heat, then simmering for 15 min in a clay pot before serving. It has been suggested to sip for one glass (250 mL) twice a day regularly [10][11][12][13]. Ya Dong Sura is typically made by soaking dried plant materials in alcohol (~1:1 v/v) for 15-30 days. This preparation should not be consumed more than two small cups (30 mL) once daily [10][11][12][13]. KSK has also been certified by the Thai Ministry of Health to be a composition of analgesic traditional household remedy [14]. As a result, herbal pharmaceutical products contained KSK are now allowed to be registered and marketed freely throughout Thailand. According to our survey on commercial KSK crude drugs (raw materials) on Thai herbal markets, they have been sold as a single herb or an ingredient mixed with other plant materials to formulate recipes, including Ya Tom and Ya Dong Sura (S1 Fig). Commercial KSK crude drug has also been claimed to improve sexual performance, in addition to the aforementioned properties. Because it has been recommended for long-term use, the toxicity of KSK crude drug should therefore be considered.
Only clinical toxicity studies for BA and ZA aqueous extracts have been reported. The BA extract revealed no acute toxicity in mice [15]. Similarly, the extract of ZA had neither acute nor chronic toxicity in Sprague-Dawley rats [16]. However, no clinical toxicity studies of SA have been described. Interestingly, it was reported that SA has been used in the preparation of arrow poison in the Malay peninsula [17,18]. Moreover, it contains strychnine-type alkaloids, i.e., spermostrychnine, strychnospermine, and their deacetyl derivatives that probably induce clonic convulsions [19,20]. These data suggest that prolonged intake of SA should be thoughtfully considered because it may cause toxicity. Since the appearances of commercial BA, SA, and ZA are similar, the identification of commercial KSK crude drugs by organoleptic analysis, which is common to persons with experience and expertise, could be misleading. Official identification methods of the origin of commercial KSK crude drugs have not previously been addressed; therefore, the development of effective methods is essential.
To provide more trustworthy and accurate identification, integrated molecular and chemical techniques, e.g., DNA barcoding and chromatography, are being applied [21]. DNA barcoding is a rapid and robust technique for species identification based on DNA sequences. It has proven to be a reliable and appropriate tool applied for identifying herbal medicines [22]. The barcode loci that can differentiate between closely related plant species such as matK, rbcL, ITS, trnH-psbA have been used [23]. To succeed in DNA barcoding, the quality of genomic DNA, the affinity of primer for amplification, sequencing methodologies, and accurate reference library are necessary to consider [22][23][24]. TLC is a powerful tool that is widely adopted for qualitative and quantitative herbal analyses [25]. It is a simple, cost-and timeeffective technique because various samples and also reference compounds can be analyzed simultaneously [26]. These combination methods enable accurate identification of several plant species with the same vernacular name, e.g., "Baimaoteng" (Aristolochia mollissima Thunb. and Solanum lyratum Hance) [27] and "Krai-Krue" (Aristolochia pierrei, A. tagala, A. pothieri, Raphistemma pulchellum, Gymnopetalum integrifollum, and Jasminum spp.) [28]. This study aims to develop and apply DNA barcoding and thin-layer chromatography (TLC) to afford reliable identification of BA, SA, and ZA, which share the identical Thai name as "KSK" crude drugs.

Collection of KSK
As presented in Table 1, under the permission of Department of National Parks, Wildlife and Plant Conservation (DNPWPC) and Queen Sirikit Botanic Garden (QSBG), Betula alnoides Buch.-Ham. ex D.Don (Betulaceae), Strychnos axillaris Colebr. (Strychnaceae), and Ziziphus attopensis Pierre (Rhamnaceae) were collected from different locations in Thailand. These plants were authenticated by Ms. Wannaree Charoensup, a botanist at the Faculty of Pharmacy, Chiang Mai University, Thailand. The authentic plant specimens (S1 Appendix) obtained from BA1-BA13, SA1-SA8, and ZA1-ZA4 were prepared, preserved, and labeled voucher specimens in the Herbarium of Faculty of Pharmacy, Chiang Mai University. The commercial KSK crude drugs, CK1-CK5, which were used for botanical origin identification, were purchased from Thai herbal markets as informed in Table 1 [29], which were ITS5S (5 0 -CCT TATCATTTAGAGGAAGGAG-3 0 ) and ITS4 (5 0 -TCCTCCGCTTATTGATATGC-3 0 ), 100 ng of genomic plant DNA, and 0.5 μL of 0.5 units of KOD FX Neo DNA polymerase (Toyobo, Japan). Deionized water was added to a final volume of 25 μL.
PCR amplification was carried out under the following cycling parameters: initial denaturation at 94˚C for 2 min followed by 30 cycles of denaturation at 94˚C for 15 s, annealing at 53˚C for 30 s, and elongation at 68˚C for 45 s. A final elongation step was performed at 68˚C for 5 min. The PCR products were visualized by gel electrophoresis using a 2.2% agarose gel. They were then subjected to purification with a MEGAquick-spin TM Plus Total Fragment DNA Purification Kit (Intron Biotechnology, Korea). A purified PCR fragment was directly sequenced using an ABI PRISM 3730 XL sequencer (Applied Biosystems, USA).
Data analysis. The obtained nucleotide sequences were examined, aligned, and manually edited using BioEdit version 7.0.5.3 and MEGA-X version 10.0.5 software. All sequences were then submitted to the GenBank database under accession number LC509575-99 (Table 1). Single nucleotide polymorphisms (SNPs) of each plant species were explored using the latter software. The discrimination power among three candidate barcodes (ITS, ITS1, and ITS2) (Fig 1) was compared using Kimura 2-parameter (K2P) distance matrices. The selected region was then tested for reliability via tree topology analysis based on the neighbor-joining (NJ) method with 1000 bootstrap replicates. Finally, DNA sequences from the selected barcode region were used as a reference database for species-specific primer design to identify commercial KSK crude drugs. Species identification of commercial KSK crude drugs by multiplex PCR. Molecular identification of CK1-CK5 was carried out using multiplex PCR. The nucleotide sequences of the selected region of the authentic medicinal plants in this study together with other medicinal plants in the genera Betula, Strychnos, and Ziziphus from the GenBank database were aligned (S1 Spreadsheet) to search the different sites in order to design species-specific primers. Three sets of species-specific primers ( Table 2), i.e., KSKITSF/BeAITS2R, KSKITSF/ StAITS2R, and KSKITSF/ZiAITS2R, as well as an internal control primer, KSKITSF/ITS4, were utilized simultaneously in multiplex PCR amplification. Strychnos nux-blanda A.W. Hill (Strychnaceae), which was authenticated by the botanist and kept in the herbarium (voucher specimen no. YS19-StN1), was used as a negative control. The reaction mixture and PCR amplification conditions used were the same as those used in the aforementioned protocol. The amplicons were finally determined by 2.2% agarose gel electrophoresis and captured on a UV tray using an imager (Gel Doc TM EZ, Bio-Rad Laboratories, USA) and Image Lab version 3.0 software.

Chemical identification
Sample extraction. The stem bark of BA1-BA5, stems of SA1-SA2, stems of ZA1-ZA2, and plant materials CK1-CK5 were cut into small pieces and dried in a hot air oven at 50˚C for 72 h. The dried samples were ground into fine powder and extracted with 95% ethanol on a shaker (SHR-2D, Daihan Scientific, Korea). After 24 h of extraction, the samples were filtered, and the filtrates were concentrated under reduced pressure. For each sample, this extraction process was repeated three times. The crude extracts of all samples were kept in a refrigerator until use.  TLC condition development. TLC analysis was carried out using activated silica gel 60 F 254 TLC plates (Merck, Germany) as the stationary phase. Betulinic acid (Chemfaces, China) and lupeol (Chemfaces, China) were utilized as standard references. Five microliters of each 1 mg/mL standard reference solution and each 15 mg/mL sample extract were loaded with a 4 mm bandwidth on TLC plates using a sample applicator (Linomat V, Camag, Switzerland) equipped with a 100 μL syringe. Then the plates were developed in the optimized mobile phase to a distance of 80 mm. After that, the dried plates were sprayed with anisaldehyde-sulfuric acid reagent before drying in an oven at 105˚C for 20 min. The TLC chromatograms were viewed and photo-recorded under UV light at 366 nm using a TLC viewer (TLC visualizer 2, Camag, Switzerland). The results from this method were orderly rearranged using visionCATS version 2.3 software.
Data analysis. Hierarchical cluster analysis (HCA), a simple and rapid clustering method, was applied for qualitative analysis to differentiate the three plant species. The TLC data were used to generate row and column data matrices. The data consisted of sample codes as rows, R f values as columns, and the number 0 (absent) and 1 (available) as the input variables (S2 Appendix). These data were clustered as a dendrogram using R Studio version 1.2.5033 software.
Species identification of commercial KSK crude drugs by TLC. Chemical identification of CK1-CK5 was analyzed using TLC chromatograms and HCA. The authentic samples of each plant species, which showed the clearest chromatograms in the prior experiment, were chosen as references for the species identification of CK1-CK5.

Generation of the DNA barcode database and DNA nucleotide analysis
DNA sequences in the ITS region obtained from the leaf extracts of the authentic plants (BA1-BA13, SA1-SA8, and ZA1-ZA4) were chosen as the DNA barcode database for KSK in this study (S2 Spreadsheet). They were submitted to the GenBank database under accession number LC509575-99 (Table 1). The lengths of the ITS regions from B. alnoides, S. axillaris and, Z. attopensis were 601 bp, 616-617 bp, and 645 bp, respectively. Moreover, there was no intraspecific variation for B. alnoides, as BA1-BA13 shared the same haplotype. Additionally, the DNA sequences of S. axillaris and Z. attopensis demonstrated single nucleotide polymorphisms (SNPs) with four haplotypes for each species (S3 Appendix). As shown in Table 3 Since commercial KSK crude drugs are available as their dried stem bark or stems form, their DNA quality and quantity might not be sufficient for amplification. To solve this problem, species-specific primers providing short amplicons were applied. Three candidate DNA barcodes, ITS, ITS1, and ITS2, were investigated to obtain a suitable database for species-specific primer design. Table 4 demonstrates the properties of ITS, ITS1, and ITS2 region alignment from authentic samples (S3 Spreadsheet). The results showed that the ITS possessed the most extended length range and alignment range of DNA, followed by the ITS1 and then the ITS2 regions. The ITS2 region expressed the most variable sites and parsimony-informative sites, whereas ITS data indicated the most conserved sites. The discrimination power of the candidate regions was then compared using genetic distance analysis based on the K2P model (S4 Spreadsheet). The percentage of K2P pairwise distances was calculated as the relative distribution of interspecific divergence (Fig 2). The K2P pairwise distances (%) of the ITS, ITS1, and ITS2 regions among B. alnoides, S. axillaris, and Z. attopensis were 0.00-0.40, 0.00-0.60, and 0.00-0.70, respectively. These results indicated that the ITS2 region exhibited the highest variation in interspecific distances, followed by the ITS1 and ITS regions. ITS2 was therefore considered to be the most suitable DNA region used for commercial KSK crude drug identification.
Because the development of actual protocols for species identification by DNA barcoding remains ambiguous, method performance should be evaluated [31]. Thus, tree topology analysis using the neighbor-joining (NJ) method was carried out. DNA sequences at the ITS2 region of BA1-BA13, SA1-SA8, and ZA1-ZA4 of this work, as well as those of other medicinal plant species in the genera Betula, Strychnos, and Ziziphus registered in the GenBank database were selected as the in-groups (S1 Spreadsheet) while Cocos nucifera retrieved from the GenBank database (GenBank accession no. HQ265515) was employed as an out-group. The NJ tree demonstrated that the obtained DNA sequences from the ITS2 region of B. alnoides (LC509575-87), S. axillaris (LC509588-95), and Z. attopensis (LC509596-99) clustered into their own clades and were separated from the other species (Fig 3). These results suggested that their nucleotide sequences in the ITS2 region were species-specific; subsequently, they were employed as reference databases for identifying the origin of commercial KSK crude drugs.

Chemical profiling of KSK
The mobile phase system that gave the best resolution and separation between each species was the combination of toluene:ethyl acetate:ammonia (1:9:0.025). The chromatograms (Fig 4)

PLOS ONE
Verification of Thai ethnobotanical medicine "Kamlang Suea Khrong" SA1-SA2, and ZA1-ZA2 (S2 Appendix), was clearly clustered among different species (Fig 4). The TLC chromatograms of BA1, SA1, and ZA1 gave the clearest separation and resolution; therefore, they were chosen as references for identifying the origin of the commercial KSK crude drugs.

Species identification of the commercial KSK crude drugs by multiplex PCR and TLC
In this study, CK1-CK5, which were unidentified by the organoleptic method and purchased from herbal markets, were identified by both molecular and chemical methods. The results from authentic samples, i.e., BA1, SA1, and ZA1, were used as references. For molecular analysis, multiplex PCR, a molecular technique that amplifies DNA using species-specific primers under a single set of PCRs [32], was established. DNA nucleotides of selected region ITS2 of the authentic samples were applied for manually designing and creating species-specific primers. As presented in Fig 5, reverse species-specific primers based on the ITS2 region, i.e., BeAITS2R, StAITS2R, and ZiAITS2R, were designed for B. alnoides, S. axillaris, and Z. attopensis, respectively. Furthermore, a forward primer, KSKITSF, was designed and used to amplify the 5.8S rDNA region, while ITS4, a universal reverse primer, was also used to amplify 26S rDNA. KSKITSF and ITS4 were used as internal controls. Then, multiplex PCR, concurrently utilizing all of these primers in a single PCR, was established to identify KSK crude drugs.
To verify the efficiency of species-specific primers, the selected authentic samples (BA1, SA1, and ZA1) and negative control sample (S. nux-blanda) were examined using multiplex PCR. As shown in Fig 6A and S3 Fig, the   of BA1, 416 and 311 bp of SA1, 413 and 240 bp of ZA1 as internal controls and species-specific amplicons, respectively. For the negative control, only the internal control amplicon was detected. These obtained results corresponding to the expected results (Fig 5) indicated that all designed primers were proper to identify CK1-CK5. The amplicons from the multiplex PCRs of CK1-CK5 in Fig 6A and S3 Fig also comprised 2 fragments, the upper and lower bands. All  upper bands gave fragment sizes similar to those of the internal control. CK1-CK2 possessed a lower band corresponding to BA1, whereas CK3-CK5 provided a lower band consistent with ZA1. These results demonstrated that CK1-CK2 originated from B. alnoides and that CK3-CK5 were derived from Z. attopensis.
Chemical identification was applied to CK1-CK5 using the TLC technique followed by HCA analysis. Separation was performed using toluene:ethyl acetate:ammonia (1:9:0.025) as the mobile phase. TLC chromatograms were observed under UV light at 366 nm after derivatization with anisaldehyde-sulfuric acid. The cluster dendrogram (S2 Appendix) and TLC chromatograms in Fig 6B show that CK1-CK2 categorized as BA1 represented B. alnoides and CK3-CK5 as ZA1 represented Z. attopensis, which corresponded to the molecular analysis-based results.

Discussion
DNA barcoding is particularly useful for differentiating closely related plant species with similar morphologies and chemical profiles. However, dried or processed materials or even plant specimens with high numbers of secondary metabolites might cause a problem in terms of the quality and quantity of DNA [33]. Although chromatography, especially TLC, is one of the official identification methods available in various herbal Pharmacopoeias, certain limitations exist because of the variation in chemical constituents in plant materials due to age, growth, and storage conditions [25,26]. In addition, closely related species with similar chemical profiles cannot be distinguished from each other [21]. To alleviate these challenges, DNA barcoding and TLC were incorporated to identify commercial KSK crude drugs in this study. This integrated approach has also previously been applied to successfully identify and control the quality of ethnobotanical medicines, for example, Marsdenia tenacissima (traditional Chinese medicine) and Galphimia glauca (Mexican medicinal species) [34,35].
Generally, DNA barcodes of unknown samples would typically be compared with reference DNA barcodes to find the matching species [32]. However, crude drugs derived from the bark or stems of plants are DNA-poor sources in living trees. Moreover, their DNA nucleotides are complicated to isolate and often highly degraded or sometimes contaminated with endophytic fungi [36]. These problems might affect DNA sequencing; therefore, species-specific primers were considered for certain successful identification [33]. To produce a DNA reference database of KSK, ITS, the commonly used region for DNA barcoding of several plant species, was selected [37]. With regard to this study, the ITS2 region expressed the highest discrimination power when performing genetic distance analysis using the K2P model. It was then subjected to a performance test for its reliability using tree topology analysis via the NJ method. The findings proved that the ITS2 region provided adequate information to identify the plant species of KSK and could distinguish between the closely related medicinal plant species. ITS2 is an ideal barcode due to its short length, ease of amplification, and highly variable discrimination of closely related species [38]. Several previous reports have indicated that the use of the ITS2 region, a mini-barcode, could successfully discriminate plants with the same vernacular name, for example, traditional Chinese medicines named Gao-ben (Ligusticum sinense and L. jeholense) [39] and Luoshiteng (Trachelospermum jasminoides, Ficus tikoua, F. pumila, and Euonymus fortune [40]. Therefore, the ITS2 region was suitable for generating species-specific primers, which were employed to identify the origin of the commercial KSK crude drugs. When designing specific primers, many aspects are taken into account including primer length, primer length difference, GC content, nucleotide composition at the primer 3 0 end, melting temperature (Tm), and melting temperature difference [41]. Furthermore, the DNA analysis of closely related species was also considered [33]. Upon careful species-specific primer design, a successful multiplex PCR method to identify the botanical origins of the commercial KSK crude drugs was obtained. This technique has also been used to authenticate other medicinal plant species such as Aralia continentalis, Angelica biserrata [42]. Moreover, it has been practically applied for identification of Atractylodes macrocephala, Glycyrrhiza uralensis, and Zingiber officinale in ginseng decoctions [43].
The chemical profiles of the crude ethanolic extracts from the bark samples of BA1-BA5, stems of SA1-SA2, and stems of ZA1-ZA2 were obtained using TLC, a time-saving, economical and uncomplicated technique [26]. The extraction method applied in this study imitated the traditional Thai preparation "Ya Dong Sura", which was formulated by soaking sliced or powdered medicinal plants in liquor. Therefore, the obtained chemical constituents would correspond to cultural consumption. The TLC fingerprints (Fig 4) of the authentic plant extracts (BA1-BA5, SA1-SA2, and ZA1-ZA2) were compared within the same species. Lupeol and betulinic acid were used as reference markers. The chromatographic pattern of BA3 was very weak when compared to BA1-BA5. This result indicated that BA3 demonstrated fewer chemical constituents than the other samples. This variation could be due to differences in environmental conditions and the age of the collected plants [44]. BA1-BA5 were gathered from northern areas of Thailand, at elevations ranging from 1,167 to 1,660 meters above sea level. The collected area of BA3 was different than BA1, BA2, and BA4, but similar to BA5. Meanwhile, BA3 was obtained from smaller and younger tree than others. Consequently, the difference in age of this plant species is a promising factor in the chemical component differences. These findings are in agreement with several prior investigations [45][46][47]. Due to the variation of the constituents, the standardization which is used to ensure the quality, including efficacy and safety, of herbal medicines is necessary. A standardization involves pre-determining one or more biochemical constituents as either active or as marker compounds [48]. The chemical markers including lupeol and betulinic acid are concentrated in this work. The comparison of chemical profiles among BA, SA, and ZA revealed that betulinic acid was the major substance in ZA, while lupeol was the main component of BA. However, betulinic acid and lupeol have also been found in other plants such as Arbutus menziesii, Lycospersicon esculentum, and Morus alba [49,50]. Therefore, these chemical substances could be used as specific markers to differentiate the origin of KSK crude drug, including BA, SA, and ZA. This chemical analysis was not addressed whether these markers could discriminate KSK crude drug from other plant species. In addition, the dissimilarity of chemical substances among these three species may result in different pharmacological effects. Additional scientific studies are needed to confirm this issue.
The plant origins of five commercial KSK crude drugs (CK1-CK5) obtained from the north and northeastern Thailand as raw materials, the dried fragmented bark or stems of the individual plants, were determined. The results of both the multiplex PCR and TLC-based techniques indicated that CK1-CK2 were BA, CK3-CK5 were ZA, and none of the commercial KSK crude drugs in this experiment were SA. The concerned SA was not observed as commercial KSK crude drugs because SA may be generally used in only specific areas, or it may be preferable to collect for use rather than for sale. Because KSK crude drug has been officially accepted by the Thai Ministry of Health to be an ingredient of analgesic traditional household remedy [14], it could be found in herbal medicine products as various dosage forms including tincture, bolus, and capsule, all of which are freely sold in Thailand despite being marketed as raw materials. The developed chemical and molecular methods could also be adapted to standardize KSK in order to produce safe and effective herbal medicinal products as well.

Conclusion
Because KSK crude drug has been recommended for long-term usage, its toxicity should be considered. To avoid the intake of SA, which may induce toxicity after prolonged use, identification of KSK crude drug is desperately needed. The findings of this study suggested simple, fast, and effective multiplex PCR and TLC chromatographic techniques that were successfully adopted for the identification of commercial KSK crude drugs distributed in Thailand. The proposed molecular and chemical identification methods could be applied for routine quality-control of KSK crude drugs for consumer reliance. Despite the fact that BA, SA, and ZA have the same name as KSK crude drug, their chemical profiles differ, potentially resulting in a wide range of pharmacological activities. Accordingly, further studies should provide insight into this issue.