DOI QR코드

DOI QR Code

Molecular Authentication of Magnoliae Flos Using Robust SNP Marker Base on trnL-F and ndhF Region

  • Kim, Min-Kyeoung (KM Fundamental Research Division, Korea Institute of Oriental Medicine) ;
  • Noh, Jong-Hun (Korean Ginseng Center for Most Valuable Products & Ginseng Genetic Resource Bank, Kyung Hee University) ;
  • Yan, Deok-Chun (Korean Ginseng Center for Most Valuable Products & Ginseng Genetic Resource Bank, Kyung Hee University) ;
  • Lee, Sanghun (KM Fundamental Research Division, Korea Institute of Oriental Medicine) ;
  • Lee, Hee-Nyeong (Korean Ginseng Center for Most Valuable Products & Ginseng Genetic Resource Bank, Kyung Hee University) ;
  • Jin, Chi-Gyu (Korean Ginseng Center for Most Valuable Products & Ginseng Genetic Resource Bank, Kyung Hee University)
  • 투고 : 2015.04.07
  • 심사 : 2015.06.17
  • 발행 : 2015.06.30

초록

Magnoliae Flos (Sini in Korean) is one of the most important oriental medicinal plants. In the Korean Herbal Pharmacopeia, the bud of the all species in Manolia denudate and Manolia genus were regarded as the botanical sources for ‘Sini’. Most the dried bud of Manolia denudata, Manolia biondii and Manolia sprengeri were used as ‘Xin-yi’ in China. Therefore, the purpose of this study was to determine and compare the ‘Magnolia’ species, four species including Manolia denudata, M. biondii, M. liliiflora and M. Kobus were analysis of sequencing data revealed DNA polymorphisms. The based on tRNA coding leucine/phenylalanine (trnL-F) and NADH-plastoquinone oxidoreductase subunit 5 (ndhF) sequences in chloroplast DNA. For the identification of ‘Magnolia’ species, polymerase chain reaction (PCR) analysis of chloroplast DNA regions such as ndhF have proven an appropriate method. A single nucleotide polymorphism (SNP) has been identified between genuine “Sini” and their fraudulent and misuse. Specific PCR primers were designed from this polymorphic site within the sequence data, and were used to detect true plants via multiplex PCR.

키워드

Introduction

Magnoliae Flos (Sini in Korean), the dried bud of Manolia denudate including Manolia genus is one of the most widely used herbal medicinal plants in Korean Herbal Pharmacopeia (KHP, 4th). Among the “Magnolia” species, four species, Manolia denudata, M. biondii, M. liliiflora and M. kobus were most used as publicly certified medicinal materials in Korean Herbal Pharmacopeia (KHP, 5th).

Magnoliae Flos has an astringent effect, which can improve local blood circulation and promote the absorption of secretions, thereby diminishing inflammation, clearing nasal passage, and relieving or eliminating symptoms (Wang et al., 1996), it also has an anti-allergic effect, which can effectively fight against the nasal itching, sneezing, runny nose, and other symptoms caused by allergic rhinitis (Kim et al., 2002). However, it is the absence of the test component for each of the four species. Also, the dried bud of four species are always misidentified as Magnoliae Flos due to their morphological similarities, most of medicinal plant products in the markets are packaged of powders of slices and Processed things that no longer bear the original morphological features of the plants (Park et al., 2006; Jigden et al., 2010; Jin et al., 2013).

Traditional authentication methods, which have relied on morphological and histological differences, are limited and quite often unreliable. Chemical analysis are significantly affected by environmental growth conditions as well as storage conditions that in comparison, DNA analysis by molecular techniques is not influenced by growth stage and environmental conditions of plants (Zhu et al., 2004).

Several molecular methods, such as arbitrarily primed polymerase chain reaction (AP-PCR) (Cheung et al., 1994), random amplified polymorphic DNA (RAPD) (Cui et al., 2003; Shim et al., 2003; Shaw and But, 1995), loop-mediated isothermal amplification (LAMP) (Sasaki et al., 2008), sequence characterized-amplified region (SCAR) (Choi et al., 2008; Wang et al., 2001), restriction fragment length polymorphism (RFLP) (Ngan et al., 1999), amplification fragment length polymorphism (AFLP) (Ha et al., 2002), and DNA microarray (Zhu et al., 2008) have previously been described. The main drawback of RAPD are lack of reliability and reproducibility. RFLP and AFLP need stringent reaction conditions and tedious operation, and they are not suitable for identification of processed sample due to the degradation of genomic DNA.

In this study, we describe general ways to differentiate medicinally important plants in a more reproducible and robust approach by analyzing Single nucleotide polymorphisms (SNP). For the SNP analysis, used regions include ribosomal genes, chloroplast genes and mitochondrial genes.

In this case, in order to determine the base sequence of trnL (tRNA-Leu) and ndhF intron region were targeted for molecular analysis and these regions were proved to be useful for discrimination of M. denudata, M. biondii, M. liliiflora and M. kobus.

Chloroplast DNA (cpDNA) intron and intergenic spacer sequence is used for evolutionary and phylogenetic study (Olmstead and Palmer 1994). The cpDNA trnL (tRNA-Leu) and ndhF intron region in land plants comprising the trnL (UAA) intron and trnL (UAA)- trnF (GAA) noncoding intergenic spacer is one of most widely used chloroplast markers for phylogenetic analysis in plants (Quandt and Stech, 2004). It has been used successfully to infer phylogenetic relationships within and among angiosperm families (Olmstead et al., 1992)

It was possible to test paternity of medicinal herbs, as it is possible to identify true Bos taurus coreanae and young antlers of deers, using genetic base sequencing (An et al., 2006). To identify plants using DNA base sequencing, chloroplast DNA or trnL-F (tRNA-Leu/Phe) /ndhF (NADH dehydrogenase subunit F) is frequently used as it is found only in plants (Lang et al., 2006).

A marker that recognizes four species of the genus Magnolia (Manolia denudata, M. biondii, M. liliiflora and M. kobus) for herbal use was developed after determining the trnL-F region and base sequencing ndhF gene and verifying its SNP.

Although traditional herbal medicine industry enrooted in medicinal plants as resources is one of the major medical organizations partnering with Western medicine, it currently imports more than 90% of raw materials from China, and there are cases of adverse reactions because of mixture of unsorted ingredients, improper applications and misunderstandings due to lack of manufacture and distribution management system. Therefore, the aim of this research is to trace the origin of Magnoliae Flos through analysis of DNA base sequence, and to develop identifying technology such as SNP markers so that there will be standardization of distribution of the herbal Magnoliae Flos.

This research promotes to verify authenticity of Magnoliae Flos and to develop identification methods, as demanded by Korea Food and Drug Administration-linked Center for Herbal Medicine.

 

Materials and Methods

Plant

The collected samples of Magnoliae Flos are shown in Table 1. Ten samples of herbal Magnoliae Flos currently distributed in Korea and China were purchased. And the standard ten samples were collected from the following places, The Ok-cheon medicinal plants authentication center of Korea Food and Drug Administration (http://agri.oc.go.kr), Kyung Hee University in Su-won, Po-cheon, and Young-cheon. Other eight samples were collected from substance separation and analysis team at Andong National University. All samples were correctly identified by the sequence of trnL-F and ndhF of NCBI-registered Magnolia plants were screened and confirmed.

Table 1.List of ‘Magnoliae Flos (Sini)’ samples used in this study

DNA extraction

A total of twenty eight collected samples were frozen using liquid nitrogen, ground using porcelain mortar and pestle, and the genomic DNA were extracted using G-spinTM Genomic DNA Extraction Kit (iNtRON, Seongnam, Koera) and each dried sample was isolated using a modified cetyltrimethylammonium bromid (CTAB) method (Murray et al., 1980).

PCR amplification of the trnL (tRNA-Leu) and ndhF (ndhF) intron region

In order to determine the base sequence of trnL-F, universal primers trnLF-c (forward) and trnLF-f (reverse) of the trnL-F region were used in PCR amplification. The base sequence of the primers were trnLF-c (5'-CGA AAT CGG TAG ACG CTA-3') and trnLF-f (5'-ATT TGA ACT GGT GAC ACG AG-3') (Taberlet et al., 1991).

The conditions of running PCR were set as the following: pre-denaturation at 96℃ for 2 minutes; denaturation at 96℃ for 30 seconds; annealing at 57℃ for 30 seconds; and extension at 72℃ for 2 minutes, for 36 cycles. Universal primers ndhF-F1 (forward) and ndhF-R1 (reverse) of the ndhF region were used in PCR amplification to determine the base sequence of ndhF. The base sequence of the primers were ndhF-F1 (5'-TTG GGA ATT GGT GGG AAT GTG-3') and ndhF-R1 (5'-TTC CTA TGG ACC CAA CGA AC-3') (Zhang et al., 2003), and the conditions of running PCR were set as the following: pre-denaturation at 95℃ for 3 minutes; denaturation at 95℃ for 30 seconds; annealing at 61℃ for 30 seconds; and extension at 72℃ for 2 minutes, for 36 cycles. Each PCR amplification was performed in a volume of 20㎕, and the reaction mixture consisted of each of the primers at a concentration of 0.5 ㎛, 20 ng of template DNA. The PCR products were migrated on a 1.0% agarose gel electrophoresis and detected by ethidium bromide staining under UV.

Sequencing and DNA sequence analysis

The PCR products were purified using a GENEALL PCR SV Purification Kit (GeneAll Biotechnology, Seoul, Koera) per the manufacturer’s instructions and then sequenced by Genotech, Inc (Genotech, Daejeon, Korea). The DNA sequence of the trnL (tRNA-Leu) and ndhF (ndhF) intron regions obtained in sequencing experiments were compiled using SeqMan software, and the sequence were edited with BioEdit program (Hall et al., 1999). Multiplex sequence alignments were performed using online Clustal W2 program (http://www.ebi.ac.uk/Tools/clustalw2/).

Design of specific primers

Specific primer were designed for M. denudata, M. liliiflora, M. Kobus and M. biondii, respectively, on the basis of the DNA polymorphisms detected (Table 2). In order to design specific primers to recognize the four species of Magnolia (M. denudata, M. liliiflora, M. kobus, M. biondii) through DNA obtained from plants and herbal Magnoliae Flos used in this research, SNP unique to ndhF of each species were determined.

Table 2.Specific primer design adjacent to SNP position using the mismatching technology

There was no SNP unique to M. kobus, and therefore “SIkbnF” specific primer was designed to recognize M. kobus and M. biondii, a specific primer “SIbinF” was designed only to identify M. biondii. Likewise, no unique SNP could be verified in M. denudata that after designing a specific primer “SIdlnF” to recognize M. denudata and M. liliiflora, a specific primer “SIlinF” was designed only to identify M. liliiflora. Specific primers were designed through “mismatching” technology (Kim et al., 2005; Shiokai et al., 2008; Hayashi et al., 2004).

Multiplex PCR

Molecular authentication of M. denudata, M. liliiflora, M. Kobus and M. biondii was performed using multiplex-PCR. Six primers (ndhF-F1, SIbinF, SIkbnF, SIdlnF, SIlinF and ndhF-R1) were used simultaneously in multiplex PCR amplification. The reaction mixture was identical to the one described earlier except the concentrations of ndhF-F1, SIkbnF, SIdlnF, SIlinF and ndhF-R1 were 0.63 ㎛, 0.31 ㎛, 0.63 ㎛, 1.13 ㎛, 0.5 ㎛ and 0.81 ㎛ respectively. PCR amplification was performed in a total volume of 20㎕. The conditions for DNA amplification were set as the following: pre-denaturation at 95℃ for 3 minutes; denaturation at 95℃ for 30 seconds; annealing at 61℃ for 30 seconds; and extension at 72℃ for 2 minutes, for 33 cycles. The PCR products were visualized on a 1% agarose gel.

Results and Discussion

Chloroplast DNA (cpDNA) sequence variations are now widely used to investigate interspecific relationships among angiosperms and other plants (palmer et al., 1988).

In this study, each cpDNA was separated and the regions of trnL-F and ndhF were amplified using PCR, in order to verify the information on the base sequence of the collected 28 samples of Magnoliae Flos. As a result, each sample yielded 950 bp and 1,248 bp of trnL-F band and ndhF band, respectively (Fig. 1, Fig. 2). As a result of analysis of base sequencing of trnL-F region, the samples were divided into three groups by the different base of SNP position (333 bp and 835 bp).

Fig. 1.Amplification of four species ‘Magnoliae Flos (Sini)’ DNA sequence (trnL-F).

Fig. 2.Amplification of four species ‘Magnoliae Flos (Sini)’ DNA sequence (ndhF).

Twenty eight plants samples of Magnoliae Flos trnL-F region DNA sequences were collected form NCBI, and aligned in Bioedit program with ClustalX program and a phylogenetic tree was built using Mega4 program. The A group, B group and most group were classified cousin relationship with M. denudate, M. Kobus and M. biondii is classified distant relationship (Fig. 3).

Fig. 3.The Neighbor-joining tree of trnL-F sequences of species belonging to the genus Magnolia.

In contrast, the analysis of base sequencing of ndhF domain divided the samples into four groups by the different base of SNP position (213 bp, 395 bp, 517 bp and 926 bp, respectively) (Table 3).

Table 3.The base of SNP position and classification of Magnoliae Flos (Sini)

The tRNA coding sequences (trnL-F), the regions cannot be used for authentication of Magnoliae Flos due to their less polymorphism. In comparison, the ndhF (NADH dehydrogenase subunit F noncoding region was targeted for molecular analysis and differentiation among Magnoliae Flos. These regions can be used for authentication of M. denudata, M. liliiflora, M. Kobus and M. biondii.

The ndhF noncoding regions of four species were PCR-amplified using the ndhF-F1 and ndhF-R1 universal primer set has 1,248 bp of band. All the herbal Magnoliae Flos that are distributed currently in Korea and China are verified as M. biondii. M. denudata and M. liliiflora are related closely in their familial lines, and this could be verified through the analysis of divergence between the Magnolia species. The difference came out to be only 1 bp out of 1,143 bp of the aligned sequence for divergence analysis (Table 4).

Table 4.Nucleotide divergences of ndhF regions

Twenty eight plants samples of Magnoliae Flos ndhF region DNA sequences were collected form NCBI, and aligned in Bioedit program with ClustalX program and a phylogenetic tree was built using Mega4 program. The A group , B group, C group and most group were classified cousin relationship with M. denudate, M. Kobus, M. liliiflora and M. biondii is classified distant relationship (Fig. 4).

Fig. 4.The Neighbor-joining tree of ndhF sequences of species belonging to the genus Magnolia.

The ndhF universal band at 1,248 bp was used as control in order to verify the presence or absence of PCR amplification. Their sequences were deposited in GeneBank (KP742958-KP74296).

Multiplex alignment of ndhF noncoding region of M. denudata, M. liliiflora, M. Kobus and M. biondii were performed with CLUSTALX program. Specific primer were designed for M. denudata, M. liliiflora, M. Kobus and M. biondii based on the detected SNP sites. Molecular discrimination of Magnoliae Flos was conducted using multiplex PCR with the six primers described. The combination of six specific primers described, as shown in Fig. 5, yield expected amplicons for different species. The specific primer “SIbinF” developed a specific band at 1,032 bp, which is amplified only in M. biondii; “SIkbnF” showed a specific band at 860 bp, amplified only in M. kobus / M. biondii; “SIdlnF” expressed at 738 bp, amplified uniquely in M. denudata / M. liliiflora; and “SIlinF” band appeared at 328 bp, amplified only in M. liliiflora (Fig. 5).

Fig. 5.Multiplex PCR for identification of four species of Magnoliae Flos (Sini).

The specific primer “SIbinF” developed a specific band at 1,032 bp, which is amplified only in M. biondii; “SIkbnF” showed a specific band at 860 bp, amplified only in M. kobus / M. biondii; “SIdlnF” expressed at 738 bp, amplified uniquely in M. denudata / M. liliiflora; and “SIlinF” band appeared at 328 bp, amplified only in M. liliiflora.

Therefore, authentic Magnoliae Flos can be identified through its specific band from the influx of commercially available herbal plants, distinguished by its SNP markers (specific primers) using multiplex-PCR. The results obtained through this research can be expected not only to verify the origin of Magnoliae Flos that is circulated in today’s market but also to recognize its species, as well as to separate from inauthentic and misused plants.

The described method has important implications in both the production and sale of this medicinal products, allowing for the Prevention of fraud and misuse, and also revealing the possible presence of other with cheaper plant material. This method is reliable, efficient, and can be used for numerous repeated tests of many medicinal plants, and the methodology presented in this study can be adapted for authentication of other medicinal materials.

참고문헌

  1. An, S.M., J.A. Ryuk, Y.H. Kim, B.C. Chae, H.J. Kim, K.H. Kim, K.K. Kang, B.S. Ko and M.Y. Lee. 2006. Genetic analysis of polygonati rhizoma and polygonati odorati rhizoma using random amplified microsatlite polymorphism. Korean J Medicinal Crop Sci. 14:125-129.
  2. Cheung, K.S., H.S. Kwan, P.P. But and P.C. Shaw. 1994. Pharmaceutical identification of American and oriental ginseng roots by genomic fingerprinting using arbitrarily primed polymerase chain reaction (AP-PCR). J Ethnopharmacol. 42:67-69. https://doi.org/10.1016/0378-8741(94)90025-6
  3. Choi, Y.E., C.H. Ahn, B.B. Kim and E.S. Yoon. 2008. Development of species specific AFLP-derived SCAR marker for authentication of Panax japonicas C. A. Meyer. Biological and Pharmaceutical Bulletin 31:135-138. https://doi.org/10.1248/bpb.31.135
  4. Cui, X.M., C.K. Lo, K.L. Yip, K.W. Dong and K.T. Sim. 2003. Authentication of Panax notoginseng by 5S-rRNA spacer domain and random amplified polymorphic DNA (RAPD) Analysis. Planta Medica 69:584-586. https://doi.org/10.1055/s-2003-40632
  5. Jigden, B., H. Wang, Y.J. Kim, N. Samdan, J.G. In and D.C. Yang. 2010. Molecular identification of oriental medicinal plant Schizonepeta tenuifolia bunge (‘Hyung-Gae’) by multiplex PCR. Plant Biotechnol Res. 4:223-228. https://doi.org/10.1007/s11816-010-0142-3
  6. Jin, C.G., M.K. Kim, J.Y. Kim, M.S. Sun, W.S. Kwon and D.C. Yang. 2013. Molecular authentication of Morus folium using mitochondrial nad7 intron 2 region. Korean J. Plant Res. 26(3):397-402. https://doi.org/10.7732/kjpr.2013.26.3.397
  7. Ha, W.Y., P.C. Shaw, J. Liu, F.C. Yau and J. Wang. 2002. Authentication of Panax ginseng and Panax quinquefolius using amplified fragment length polymorphism (AFLP) and directed amplification of minisatellite region DNA (DAMD). Journal of Agricultural and Food Chemistry 50:1871-1875. https://doi.org/10.1021/jf011365l
  8. Hall, T.A. 1999. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. Ser. 41:95-98.
  9. Hayashi, K., N. Hashimoto and M. Daigen. 2004. Development of PCR-based SNP markers for rice blast resistance genes at the Piz locus. Theor. Appl. Genet. 108(7):1212-1220. https://doi.org/10.1007/s00122-003-1553-0
  10. Kim, H.M., J.M. Yi and K.S. Lim. 2002. Magnoliae Flos inhibits mast cell-dependent immediate-type allergic reactions. Pharmacological Res. 32(2):107-111.
  11. Kim, M.Y., K. Van and P. Lestari. 2005. SNP identification and SNP marker development for a GmNARK gene controlling supernodulation in soybean. Theor. Appl. Genet. 110(6):1003-1010. https://doi.org/10.1007/s00122-004-1887-2
  12. Kim, M.K., J. Baigalmaa, H. Sun, J.H. Noh, S.Y. Kim and D.C. Yang. 2008. Phylogenetic analysis of Ji-Mo (Anemarrhena asphodeloides) on the basis of chloroplast DNA sequences. Korean J Medicinal Crop Sci. 16(1):20-26.
  13. Kim, M.K., G.H. Jang, D.C. Yang, S.H. Lee, H.N. Lee and C.G. Jin. 2014. Molecular authentication of Acanthopanacis cortex by multiplex-PCR analysis Tools. Korean J. Plant Res. 27(6):680-686. https://doi.org/10.7732/kjpr.2014.27.6.680
  14. Lang, P., F. Dane and T.L. Kubisiak. 2006. Phylogeny of Castanea (Fagaceae) based on chloroplast trnL-F sequence data. Tree Genetics & Genomes 2:132-139. https://doi.org/10.1007/s11295-006-0036-2
  15. Murray, M.G and W.F. Thompson. 1980. Rapid isolation of high molecular weight plant DNA. Nucleic acids Research 8:4321-4325. https://doi.org/10.1093/nar/8.19.4321
  16. Ngan, F., P. Shaw, P.P. But and J. Wang. 1999. Molecular authentication of Panax species. Phytochemistry 50:787-791. https://doi.org/10.1016/S0031-9422(98)00606-2
  17. Olmstead, R.G., R. Scotland, S. Wagstaff, J. Sweere and P. Reeves. 1992. Application of the chloroplast gene ndhF to angiosperm phylogenetic studies. Plant Molecular Evolution Newsletter 2: 27-30.
  18. Olmstead, R.G and J.D. Palmer. 1994. Chloroplast DNA systematics: a review of methods and data analysis. American Journal of Botany 81:1205-1224. https://doi.org/10.2307/2445483
  19. Palmer, J.D., R.K. Jansen, H.J. Michaes, M.W. Chase and J.R. Manhart. 1988. Chloroplast DNA variation and plant phylogeny. Ann Missouri Bot Garden 75:1180-1206. https://doi.org/10.2307/2399279
  20. Park, M.J., M.K. Kim, J.G. In and D.C.Yang. 2006. Molecular identification of Korean ginseng by amplification refractory mutation system-PCR. Food Res. Int. 39:568-574. https://doi.org/10.1016/j.foodres.2005.11.004
  21. Quadt, D and M. Stech. 2004. Molecular evolution of the trnTUGU trnFGAA region in Bryophtes. Plant Biology 6(5):545-554. https://doi.org/10.1055/s-2004-821144
  22. Sasaki, Y., K. Komatsu and S. Nagumo. 2008. Rapid detection of Panax ginseng by loop-mediated isothermal amplification and its application to authentication of ginseng. Biological and Pharmaceutical Bulletin 31:1806-1808. https://doi.org/10.1248/bpb.31.1806
  23. Shiokai, S., H. Kitashiba and K. Shirasawa. 2008. Leaf-punch method to prepare a large number of PCR templates from plants for SNP analysis. Molecular Breeding 23(2):329-336. https://doi.org/10.1007/s11032-008-9244-9
  24. Shaw, P.C and P.P. But. 1995. Authentication of Panax species and their adulterants by random-primed polymerase chain reaction. Planta Medica 61:466-469. https://doi.org/10.1055/s-2006-958138
  25. Shim, Y.H., J.H. Choi, C.D. Park, C.J. Lim, J.H. Cho and H.J. Kim. 2003. Molecular differentiation of Panax species by RAPD analysis. Archives of Pharmacal Research 26:601-605. https://doi.org/10.1007/BF02976708
  26. Taberlet, P., L. Gielly, G. Pautou and J. Bouvet. 1991. Universal primers for amplification of three non-coding region of chloroplast DNA. Plant Molecular Biology 17:1105-1109. https://doi.org/10.1007/BF00037152
  27. Wang, B., C.G. Zhang, P.F. Gao, X.M. Wu and Y. Zhao. 1996. Research progress on Nidus vespae, a traditional Chinese medicine derived from insents. Journal of Pharmaceutical and Scientific Innovation 2(6):1-9.
  28. Wang, J., W.Y. Ha, F.N. Ngan, P.P. But and P.C. Shaw. 2001. Application of sequence characterized amplified region (SCAR) analysis to authenticate Panax species and their adulterants. Planta Medica 67:781-783. https://doi.org/10.1055/s-2001-18340
  29. Zhang, W.H., Z.D. Chen, J.H. Li, H.B. Chen and Y.C. Tang. 2003. Phylogeny of the Dipsacales s.I. based on chloroplast trnL-F and ndhF sequences. Molecular Phylogenetics and Evolution 26:176-189. https://doi.org/10.1016/S1055-7903(02)00303-2
  30. Zhu. S., H. Fushimi, S. Cai and K. Komatsu. 2004. Species identification from ginseng drugs by multiplex amplification refractory mutation system (MARMS). Planta Medica 70:189-192. https://doi.org/10.1055/s-2004-815502
  31. Zhu. S., H. Fushimi and K. Komatsu. 2008. Development of a DNA microarray for authentication of ginseng drugs based on 18S rRNA gene sequence. Journal of Agricultural and Food Chem. 56:3953-3959. https://doi.org/10.1021/jf0732814