- Open Access
Molecular identification of Uncaria (Gouteng) through DNA barcoding
Chinese Medicine volume 11, Article number: 3 (2016)
While DNA barcoding is an important technology for the authentication of the botanical origins of Chinese medicines, the suitable markers for DNA barcoding of the genus Uncaria have not been reported yet. This study aims to determine suitable markers for DNA barcoding of the genus Uncaria (Gouteng).
Genomic DNA was extracted from the freshly dried leaves of Uncaria plants by a Bioteke’s Plant Genomic DNA Extraction Kit. Five candidate DNA barcode sites (ITS2, rbcL, psbA–trnH, ITS, and matK) were amplified by PCR with established primers. The purified PCR products were bidirectionally sequenced with appropriate amplification primers in an ABI-PRISM3730 instrument. The candidate DNA barcodes of 257 accessions of Uncaria in GenBank were aligned by ClustalW. Sequence assembly and consensus sequence generation were performed with CodonCode Aligner 3.7.1. The identification efficiency of the candidate DNA barcodes was evaluated with BLAST and nearest distance methods. The interspecific divergence and intraspecific variation were assessed by the Kimura 2-Parameter model. Genetic distances were computed with Molecular Evolutionary Genetics Analysis 6.0.
The accessions of the five candidate DNA barcodes from 11 of 12 species of Uncaria in China and four species from other countries were included in the analysis, while 54 of total accessions were submitted to GenBank. In a comparison of the interspecific genetic distances of the five candidate barcodes, psbA–trnH exhibited the highest interspecific divergence based on interspecific distance, theta prime, and minimum interspecific distance, followed by ITS2. The distribution of the interspecific distance of ITS2 and psbA–trnH was higher than the corresponding intraspecific distance. Additionally, psbA–trnH showed 95.9 % identification efficiency by both the BLAST and nearest distance methods regardless of species or genus level. ITS2 exhibited 92.2 % identification efficiency by the nearest distance method, but 87 % by the BLAST method.
While psbA–trnH and ITS2 (used alone) were applicable barcodes for species authentication of Uncaria, psbA–trnH was a more suitable barcode for authentication of Uncaria macrophylla.
Uncaria rhynchophylla (Miq.) Jacks is used to treat convulsion, hypertension, epilepsy, eclampsia, migraine, and cerebral diseases [1–3]. Rhynchophylline, isorhynchophylline, corynoxeine, and isocorynoxeine are the major components of U. rhynchophylla . Oleanane and ursane-type triterpenes, (including uncarinic acids, ursolic acid, 3-hydroxyurs-12-en-27,28-dioic acid, hyperin, and catechin) were found in Uncaria [1, 5]. Uncaria comprises 34 species , 10 of which are found in the Guangxi Zhuang Autonomous Region. Among the 10 species of Uncaria in Guangxi, U. rhynchophylla and Uncaria macrophylla are the most widely and abundantly distributed . Stems with hooks from several species of Uncaria, including U. rhynchophylla, U. macrophylla, Uncaria hirsuta, Uncaria sinensis, and Uncaria sessilifructus, have been used in Chinese medicine (CM) preparations, Gouteng in Chinese. Only the above five species plants of the genus Uncaria can serve as the botanical origins of Gouteng according to the Chinese Pharmacopeia (10th edition) . Adulterants of Gouteng include Uncaria laevigata, Uncaria lancifolia, Uncaria scandens, Uncaria rhynchophylloides, and Uncaria homomalla [7, 9], due to similar organoleptic characteristics to those of U. rhynchophylla. But their chemical constituents and therapeutic effects are distinct from those of U. rhynchophylla [2, 10, 11].
DNA barcoding can accurately identify species on the basis of short standardized genes or DNA regions [12, 13], without confounding factors such as environmental influence, growth phase, and morphological diversity within species [14–16]. The mitochondrial gene encoding cytochrome c oxidase subunit 1 (co1) is a potential DNA barcode in most animal species as well as some fungal species. However, the co1 gene and other mitochondrial genes from plants have limited use in identifying plant species across a wide range of taxa, due to their low genetic variations and variable mitochondrial genomes . Several DNA regions, such as ITS2, psbA–trnH, matK, rbcL, ITS, ycf5, and rpoC1 [14, 18–21] have been evaluated as potential DNA barcodes in medicinal plants. Among these candidate barcoding loci, the ITS2 locus not only had the highest identification efficiency among all tested regions, but also discriminated a wide range of plant taxa [14, 22]. By contrast, ITS1 was a useful barcode for identifying Salvia species . The psbA–trnH intergenic region was a suitable DNA marker for identification of flowering plants [17, 18], pteridophytes , Lonicera japonica Thunb from Caprifoliaceae , and aquatic plant species .
The authentication of the botanical origins of Gouteng is based on the morphological characteristics, microscopic structures, or chemical components of specimens . The accuracy is often affected by environmental and subjective factors, especially for dry medicinal materials from different origins . Chemical analysis methods, such as high-performance liquid chromatography (HPLC) and HPLC coupled with quadrupole time-of-flight mass spectrometry, have also been studied . Multiple genetic molecular markers have been used to screen Uncaria, such as random amplified polymorphic DNA (RAPD) and rDNAs (including 5.8S rDNA, ITS1, and ITS2) .
This study aims to determine suitable markers for DNA barcoding of the genus Uncaria. In this study, five candidate loci (ITS2, rbcL, psbA–trnH, ITS, and matK) were tested for their potential as DNA barcodes for Uncaria.
Fifty-four sequences from our laboratory (all submitted to GenBank), among which 12 samples of six species of Uncaria (U. rhynchophylla, U. macrophylla, U. sessilifructus, U. hirsuta, U. lancifolia and U. homomalla) are used as Gouteng in CM markets, were collected from areas in Guangxi Province, including Rongshui, Sanjiang, Shanglin, Ningming, and Jinxi county, Nanning Sitang town, and Guangxi Medicinal Botanical Garden, in 2009 and 2010 by Professor Ruisong Huang. The plant species were identified by Shouyang Liu, Yiling Zhu, and Kejian Yan through morphological characteristics and analysis of microscopic structures [7, 10]. All of the voucher specimens (all the voucher numbers can be seen in Table 1) were deposited in the Key Laboratory of Biological Molecular Medicine Research of Guangxi Higher Education, Guangxi Medical University.
In total, 257 accessions related to the five candidate DNA barcoding sites (ITS2, rbcL, psbA–trnH, ITS, and matK) from 89 samples belonging to 15 species of Uncaria were analyzed in this study. All accession data were downloaded from GenBank, except for the above 54 sequences, which were amplified and sequenced in our laboratory. All datasets of Uncaria species used in the study contained more than two samples, except for Uncaria africana, Uncaria guianensis, and Uncaria lanosa. Some accessions in which the sequences contained undetermined bases or were from sp. species (taxa of species unclear or unnamed) were not selected. In this study, the correctness of the accessions downloaded from GenBank was tested through blasting against those of congener plants. Only the sequences with both a similarity ratio and query cover ratio higher than 90 % in the same species were suitable for selection. However, some accessions containing inversion sequences were collected in this dataset because they could influence the sequence divergence and supply some important genetic characters . The total data and sample information used in this study are shown in Table 1.
DNA extraction, PCR amplification, and sequencing
In this study, genomic DNA was extracted from the freshly dried leaves of Uncaria plants by the improved protocol of a new rapid Plant Genomic DNA Extraction Kit (centrifugal column type, DP3112; Bioteke Corporation, Beijing, China). The Uncaria leaves were ground in liquid nitrogen, and the cell nuclear separation solution (3 ml for 0.5 g sample) was immediately added to the samples to remove impurities from the cytoplasm before the cell nuclei were lysed . PCR amplification of the five candidate DNA barcode sites was performed in a Tprofessional Gradient 96 Type (Biometra, Göttingen, Germany) with approximately 30 ng of genomic DNA as a template in a 25-µL reaction mixture. Each reaction contained 1 × PCR buffer (2.0 mM MgCl2, 0.2 mM each dNTP, 0.1 µM each primer; synthesized by Sangon Biotech, Co., Ltd., Shanghai, China), and 1.0 U Taq DNA polymerase (TaKaRa Biotechnology Co., Ltd., Dalian, China). The primers and reaction conditions used were the same as those used by Chen et al. . The PCR products were electrophoresed in a 1.5 % agarose gel in 1 × TAE buffer, then purified with a TIANGel Midi Purification Kit (Tiangen Biotech Co. Ltd, Beijing, China). The purified PCR products were bidirectionally sequenced with appropriate amplification primers (Additional file 1) in an ABI-PRISM3730 instrument (Thermo Fisher Scientific, MA, USA) by Sangon Biotech, Co., Ltd., Shanghai, China.
Sequence alignment and data analysis
Sequence assembly and consensus sequence generation were performed by CodonCode Aligner 3.7.1 (CodonCode Co., MA, USA) by trimming the low quality sequence and primer areas. The matK and rbcL regions were delimited by alignment with known sequences in databases by CodonCode Aligner. After removal of the psbA and trnH genes at the ends of psbA–trnH, the boundary of the psbA–trnH intergenic spacer was determined according to the annotations of similar sequences in GenBank. The five candidate DNA barcodes were aligned by ClustalW (EMBL-EBI, Heidelberg, German). Kimura 2-Parameter (K2P) genetic distances were computed with Molecular Evolutionary Genetics Analysis 6.0 (The Biodesign Institute, AZ, USA) . All interspecific and intraspecific distances, including theta prime, minimum interspecific distance, theta, and coalescent depth for all accessions of each locus, were calculated and compared to evaluate the interspecific divergence and intraspecific variation by the K2P model. Meanwhile, statistical analysis of the distribution divergency of the genetic distance between different sequences was performed through the Wilcoxon signed-rank test to assess the barcoding gap for different candidate loci with SPSS software (SPSS 16.0: International Business Machines Corporation Statistical Product and Service Solutions, Armonk, New York, USA), which the test statistical W+ and W− were calculated for two side test, as described previously [14, 22]. The BLAST1 and nearest distance methods were used to evaluate the species identification efficiency [32, 33].
PCR amplification and base composition of the five loci of Uncaria
The sequence length and GC content of the five candidate loci (ITS2, rbcL, psbA–trnH, ITS, and matK) were obtained from the CodonCode Aligner and Clustal W alignment results (Table 2). The GC content of psbA–trnH was the lowest, while that of ITS2 was the highest. The variability of the length range of the psbA–trnH intergenic spacer was greater than that of the other candidates. The psbA–trnH region of U. macrophylla was more divergent than that of the other Uncaria plants.
Genetic interspecific divergence and intraspecific variation
Six parameters (Table 3) represented the genetic divergences of species in Uncaria. In a comparison of the intraspecific distances of the five candidate barcodes among Uncaria species, the intraspecific distance of psbA–trnH was higher than that of the other loci at the species level. Meanwhile, the interspecific genetic distance of the psbA–trnH intergenic spacer exhibited the highest divergence according to the interspecific distance, theta prime, and minimum interspecific distance. The interspecific distance of ITS2 was the second highest after psbA–trnH. All interspecific divergences of ITS2, psbA–trnH, and ITS were greatly higher than the corresponding intraspecific divergences. Furthermore, the overall mean distance of psbA–trnH was the highest among the five loci (Fig. 1).
The psbA–trnH intergenic spacer had the highest interspecific divergence among all the loci based on the Wilcoxon signed-rank test. The second highest interspecific divergence was shown by ITS2. The scale of the interspecific divergence of psbA–trnH was higher than ITS2, ITS, matK and rbcL, respectively (all P < 0.001), that of ITS2 was higher than ITS, matK and rbcL, respectively (all P < 0.001, Table 4). Furthermore, the intraspecific divergences between ITS and matK, rbcL and matK, ITS2 and matK, psbA–trnH and matK, and ITS and rbcL did not exhibit any significant differences (P > 0.05, Table 5).
Analysis of barcoding gaps
As a barcode for identifying botanical species, the divergence between species should be higher than the variation within species . Although the histogram of the K2P genetic distance analysis revealed a partial overlap “barcoding gap” between the intraspecific and interspecific divergence of ITS2 or psbA–trnH (Fig. 2), the intraspecific variation of psbA–trnH and ITS2 was considerably lower than the distribution of their interspecific divergence. The genetic divergence distribution of ITS was similar to that of ITS2. No clear “barcoding gap” corresponding to the rbcL or matK loci was observed, wherein the genetic distance distribution of more than 90 % of accessions was less than 0.020. However, the distribution of the interspecific divergence of ITS2 and psbA–trnH provided a better resolution than that of rbcL and matK.
Identification efficiency and characteristics of Clustal W alignment
The BLAST and nearest distance methods were employed to test the applicability of the five loci for species identification of Uncaria. psbA–trnH presented 95.9 % identification efficiency with both the BLAST and nearest distance methods at the species or genus level. ITS2 exhibited 92.2 % identification efficiency by the nearest distance method, but 87 % by the BLAST method, whereas rbcL showed only 76.2 % by the nearest distance method and 42.9 % by the BLAST method (Table 2). Meanwhile, psbA–trnH of U. macrophylla exhibited more obvious characteristics than U. rhynchophylla and the other species tested (Figs. 3, 4, 5). Two insertion fragments existed in the psbA–trnH sequence of U. macrophylla, including a serial seven A fragment at 171–177 bp, and another double repeat “ATTAAA” at 234–247 bp. The psbA–trnH intergenic spacer can be used as a barcode for the identification of Uncaria plants. The phylogeny of Uncaria ITS2 (computed model: Maximum Composite Likelihood)  showed that only four accessions (4/77 accessions) were in the incorrect taxonomic category (Fig. 6), which was less than the other loci tested. Thus, ITS2 could be another suitable DNA barcode for Uncaria.
Significance of authentication of Uncaria by DNA barcoding
Gouteng is commonly exploited as the major ingredient herb of CM prescriptions for hypertension or migraine treatment [2, 35]. The amount of stems with hooks of U. rhynchophylla (Gouteng) required in traditional clinic and pharmaceutical production, has been increased; while the natural growth of U. rhynchophylla, U. hirsuta and U. macrophylla which could serve as the botanical origins of Gouteng was limited with the rising of collection. Some other species of the genus Uncaria are often collected to adulterate Guoteng, such as U. laevigata, U. lancifolia, U. scandens . Therefore, the correct genotypic identification of Uncaria plant material is essential in order to protect public health and for industrial production.
Although some methods have been developed to distinguish Uncaria plants based on morphotype, microcharacter, or physical and chemical reactions [8, 9], these are dependent on taxonomy experts. Currently, the genetic molecular markers for the genus Uncaria were related to RAPD, rDNA, and ITS, while DNA barcoding assays have not yet been reported. This study included 11 of 12 species of Uncaria in China, with U. rhynchophylloides missing in the screen for suitable DNA barcodes for Uncaria.
In the present study, psbA–trnH presented 95.9 % identification efficiency for Uncaria accessions tested with both BLAST and nearest distance methods at the species or genus level. ITS2 also exhibited high identification efficiency at 92.2 or 87 % with the nearest distance or BLAST method, respectively.
Quality and amplification efficiency of DNA from Uncaria
The DNA of Uncaria was not extracted efficiently, due to the large amounts of polysaccharides, polyphenols, and alkaloids present in the samples. A cell nuclear separation solution was used to remove the impurities from genomic DNA . The quality of the DNA extracted from the Uncaria plants satisfied the requirements for PCR amplification and sequencing. The efficiency of both PCR amplification and sequencing for psbA–trnH was the highest among the five candidate loci. Specifically, PCR amplification showed 96.7 % efficiency, while sequencing showed 100 % efficiency. Because the average GC content of ITS2 was 66.3 %, which was higher than that of the other loci, the resulting DNA extract was slightly difficult to amplify.
Selection of candidate DNA barcodes
In this study, the length of psbA–trnH of Uncaria ranged from 235 to 315 bp (mean 287 bp), which was longer than that of ITS2, but shorter than that of rbcL, ITS, and matK. Additionally, psbA–trnH of Uncaria exhibited the highest interspecific divergence among the five loci tested, based on the results of six parameters of the K2P model or Wilcoxon signed-rank test of interspecific divergence. The interspecies divergence of psbA–trnH was higher than the relevant intraspecies variation. Furthermore, psbA–trnH of U. macrophylla was significantly distinct from that of U. rhynchophylla and the other species because of two insertion fragments: one was a seven A repeat inserted at 171–177 bp and the other was two cis-repeats of ATTAAA at 233–247 bp (Figs. 3, 4, 5). Although one TAAAAAA repeat was observed at 171–177 bp in psbA–trnH from Uncaria yunnanensis, no double cis-repeats of ATTAAA were observed at 233–247 bp. Meanwhile, one inversion sequence of length 73–74 bp with identity ratios of more than 98 % in psbA–trnH of Uncaria was found in this study (Additional file 2). The intragenic variation of the genus Uncaria was large because of this inversion phenomenon existing in psbA–trnH. This situation was also observed in psbA–trnH of Aconitum L. . The characteristics of the insertion sequences in psbA–trnH could effectively authenticate Uncaria species.
ITS2 was another suitable locus for distinguishing different species of Uncaria. The length range of ITS2 was 210–221 bp (mean 219.9 bp), which was the shortest among the five loci. Consequently, 95.8 % efficiency could be reached by PCR amplification. In a comparison of the interspecific genetic distances of the five candidate barcodes among Uncaria species, the mean interspecific distance of ITS2 was higher than its mean intraspecific divergence, and the values were second only to those of psbA–trnH (Table 3). Based on the phylogenetic analysis of ITS2 by the neighbor-joining method and the evolutionary distances computed by the Maximum Composite Likelihood model, more than 93 % of Uncaria at the species level in this study were divided into monophyla as recognized species. Among 77 accessions of ITS2, comprising 14 species of Uncaria, only four accessions were in an incorrect taxonomic category, according to the construction of a phylogenetic tree for ITS2 (Fig. 6). Uncaria manifested complex morphological features and genetic backgrounds, and even some specimens with obvious differences in appearance possessed similar ITS sequences . This could explain the existence of some accessions that appeared in different monophyla from their original morphological taxa. Some species submitted to GenBank may have been wrongly categorized. Sequences with lengths of less than 100 bp, those with ambiguous bases containing more than one “N”, or those belonging to unnamed species (such as those with spp. and aff. in the species name) were excluded  from this study to guarantee the reliability of the selected sequences.
A better “barcoding gap” was observed between the interspecific divergence and intraspecific variation of ITS2 compared with the other loci. ITS, which contained three fragments (ITS1, 5.8S rDNA, ITS2), exhibited a similar identification efficiency to that of ITS2. Both rbcL and matK were unsuitable genetic loci for authentication of the botanical origins of Gouteng, because of the absence of a clear barcoding gap between the interspecific divergence and intraspecific variation by the K2P model. The overall mean distance of rbcL was only 0.002 and that for matK was 0.005, as computed by the Maximum Composite Likelihood model (Fig. 1). Moreover, we found that the combination of psbA–trnH with ITS2 would provide a better result for the authentication of Uncaria plants, and could even distinguish between incorrect and correct taxa or identify some cryptic species. Currently, a preliminary system for DNA barcoding of herbal materials has been established based on a two-locus combination of ITS2 and psbA–trnH barcodes . Recently, ITS2 was successfully exploited in a survey involving commercial Rhodiola products, including decoction pieces .
psbA–trnH and ITS2 also exhibited high authentication power for different species of Uncaria. Both psbA–trnH and ITS2 revealed the distinct divergence of U. macrophylla from U. rhynchophylla and the other species at the species level.
While psbA–trnH and ITS2 (used alone) were applicable barcodes for species authentication of Uncaria, psbA–trnH was a more suitable barcode for authentication of U. macrophylla.
internal transcribed spacer
- psbA–trnH :
gene spacer between psbA and trnH in chloroplast DNA
high-performance liquid chromatography
random amplified polymorphic DNA
Zhang YB, Yang WZ, Yao CL, Feng RH, Yang M, Guo DA, Wu WY. New triterpenic acids from Uncaria rhynchophylla: chemistry, NO-inhibitory activity, and tandem mass spectrometric analysis. Fitoterapia. 2014;96:39–47.
Ndagijimana A, Wang X, Pan G, Zhang F, Feng H, Olaleye O. A review on indole alkaloids isolated from Uncaria rhynchophylla and their pharmacological studies. Fitoterapia. 2013;86:35–47.
Xian YF, Lin ZX, Mao QQ, Hu Z, Zhao M, Che CT, Ip SP. Bioassay-Guided Isolation of Neuroprotective Compounds from Uncaria rhynchophylla against Beta-Amyloid-Induced Neurotoxicity. Evid Based Complement Alternat Med. 2012;2012:802625.
Qu J, Gong T, Ma B, Zhang L, Kano Y, Yuan D. Comparative study of fourteen alkaloids from Uncaria rhynchophylla hooks and leaves using HPLC-diode array detection-atmospheric pressure chemical ionization/MS method. Chem Pharm Bull (Tokyo). 2012;60(1):23–30.
Shimada Y, Goto H, Kogure T, Shibahara N, Sakakibara I, Sasaki H, Terasawa K. Protective effect of phenolic compounds isolated from the hooks and stems of Uncaria sinensis on glutamate-induced neuronal death. Am J Chin Med. 2001;29(1):173–80.
Flora of China (Rubiaceae) [http://www.efloras.org/florataxon.aspx?flora_id=3&taxon_id=134135].
Huang RS, Zhu YL, Yan KJ. Investigation on Plant Resources of Uncaria Schreb in Guangxi and Folk Application in Zhuang and Yao. Guangxi YiKe DaXue XueBao. 2012;29(2):306–8.
Commission CP. The Pharmacopoeia of the People’s Republic of China, vol. 1. 10th ed. Beijing: China Medical Science Press; 2010. p. 240.
He Q, Zeng X. The authentication about the appearance characters of Gouteng and its adulterants. ShiYong ZhongYiYao ZaZhi. 2011;29(4):279.
Huang RS, Zhang P, Zhu YL, Yan KJ, Qin DJ. The analysis of Rhynchophylline contents in part of Uncaria Schreb plant breeds. ZhongGuo ShiYan FangJiXue ZaZhi. 2012;18(8):125–8.
Gattuso M, Di Sapio O, Gattuso S, Pereyra EL. Morphoanatomical studies of Uncaria tomentosa and Uncaria guianensis bark and leaves. Phytomedicine. 2004;11(2–3):213–23.
Hebert PD, Gregory TR. The promise of DNA barcoding for taxonomy. Syst Biol. 2005;54(5):852–9.
Schindel DE, Miller SE. DNA barcoding a useful tool for taxonomists. Nature. 2005;435(7038):17.
Chen S, Yao H, Han J, Liu C, Song J, Shi L, Zhu Y, Ma X, Gao T, Pang X, et al. Validation of the ITS2 region as a novel DNA barcode for identifying medicinal plant species. PLoS One. 2010;5(1):e8613.
Pang X, Song J, Zhu Y, Xie C, Chen S. Using DNA Barcoding to Identify Species within Euphorbiaceae. Planta Med. 2010;76(15):1784–6.
Techen N, Parveen I, Pan Z, Khan IA. DNA barcoding of medicinal plant material for identification. Curr Opin Biotechnol. 2014;25C:103–10.
Kress WJ, Wurdack KJ, Zimmer EA, Weigt LA, Janzen DH. Use of DNA barcodes to identify flowering plants. Proc Natl Acad Sci USA. 2005;102(23):8369–74.
Kress WJ, Erickson DL. A two-locus global DNA barcode for land plants: the coding rbcL gene complements the non-coding trnH–psbA spacer region. PLoS One. 2007;2(6):e508.
Newmaster SG, Fazekas AJ, Steeves RA, Janovec J. Testing candidate plant barcode regions in the Myristicaceae. Mol Ecol Resour. 2008;8(3):480–90.
Yao H, Song J, Liu C, Luo K, Han J, Li Y, Pang X, Xu H, Zhu Y, Xiao P, et al. Use of ITS2 region as the universal DNA barcode for plants and animals. PLoS One. 2010;5(10):e13102.
Sun Z, Gao T, Yao H, Shi L, Zhu Y, Chen S. Identification of Lonicera japonica and its Related Species Using the DNA Barcoding Method. Planta Med. 2011;77:301–6.
Gao T, Yao H, Song J, Liu C, Zhu Y, Ma X, Pang X, Xu H, Chen S. Identification of medicinal plants in the family Fabaceae using a potential DNA barcode ITS2. J Ethnopharmacol. 2010;130(1):116–21.
Wang M, Zhao HX, Wang L, Wang T, Yang RW, Wang XL, Zhou YH, Ding CB, Zhang L. Potential use of DNA barcoding for the identification of Salvia based on cpDNA and nrDNA sequences. Gene. 2013;528(2):206–15.
Ma XY, Xie CX, Liu C, Song JY, Yao H, Luo K, Zhu YJ, Gao T, Pang XH, Qian J, et al. Species identification of medicinal pteridophytes by a DNA barcode marker, the chloroplast psbA–trnH intergenic region. Biol Pharm Bull. 2010;33(11):1919–24.
Ghahramanzadeh R, Esselink G, Kodde LP, Duistermaat H, van Valkenburg JL, Marashi SH, Smulders MJ, van de Wiel CC. Efficient distinction of invasive aquatic plant species from non-invasive related species using DNA barcoding. Mol Ecol Resour. 2013;13(1):21–31.
Zhu S, Zhou L, Pang H, Huang H, Gao X, Zeng C. Molecular identification of one Uncaria plant. Zhongguo ZhongYao ZaZhi. 2011;36(5):535–7.
Xie S, Shi Y, Wang Y, Wu C, Liu W, Feng F, Xie N. Systematic identification and quantification of tetracyclic monoterpenoid oxindole alkaloids in Uncaria rhynchophylla and their fragmentations in Q-TOF-MS spectra. J Pharm Biomed Anal. 2013;81–82:56–64.
Gang T, Liu T, Zhu Y, Liu ZY. Molecular identification of medicinal plant genus Uncaria in Guizhou. Zhong Yao Cai. 2008;31(6):825–8.
He J, Wong KL, Shaw PC, Wang H, Li DZ. Identification of the medicinal plants in Aconitum L. by DNA barcoding technique. Planta Med. 2010;76(14):1622–8.
Xu P, Wu YS, Huang RS, Zhou J, Luo Y, Li KZ, Tang YL, Liu Y. Polymorphism of Three Kinds of Medicinal Plants Uncaria from Guangxi by RAPD. ShiZhen GuoYi GuoYao. 2012;23(3):3.
Tamura K, Dudley J, Nei M, Kumar S. MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007;24(8):1596–9.
Meyer CP, Paulay G. DNA barcoding: error rates based on comprehensive sampling. PLoS Biol. 2005;3(12):e422.
Lahaye R, van der Bank M, Bogarin D, Warner J, Pupulin F, Gigot G, Maurin O, Duthoit S, Barraclough TG, Savolainen V. DNA barcoding the floras of biodiversity hotspots. Proc Natl Acad Sci USA. 2008;105(8):2923–8.
Hebert PD, Ratnasingham S, deWaard JR. Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species. Proc Biol Sci. 2003;270(Suppl 1):S96–9.
Chang YY, Tsai YT, Lai JN, Yeh CH, Lin SK. The traditional Chinese medicine prescription patterns for migraine patients in Taiwan: a population-based study. J Ethnopharmacol. 2014;151(3):1209–17.
Chen S, Pang X, Song J, Shi L, Yao H, Han J, Leon C. A renaissance in herbal medicine identification: from morphology to DNA. Biotechnol Adv. 2014;32(7):1237–44.
Xin T, Li X, Yao H, Lin Y, Ma X, Cheng R, Song J, Ni L, Fan C, Chen S. Survey of commercial Rhodiola products revealed species diversity and potential safety issues. Sci Rep. 2015;5:8337.
YSW conceived and designed the study. YLT, RSH, PX, KZL extracted DNA and performed PCR; YLT, DZC and YLuo checked the quality of PCR products and analyzed the sequencing results; YLT, NXC and YLiu analyzed the data; YSW wrote and revised the manuscript. All authors read and approved the final manuscript.
We are grateful to Mr. Shou-yang Liu and Mr. Yi-ling Zhu, Mr. Ke-jian Yan for their help with the collection and identification of the Uncaria specimens. We appreciate the essential comments provided by Dr. Hui Yao and all others that much improved the manuscript. We thank the Institute of Medicinal Plant Development, Chinese Academy of Medical Sciences, P.R.China, for providing the calculation and analysis platform. The work was supported by Guangxi Natural Science Foundation (2010GXNSFA013145); Scientific Research Project Foundation for Higher Education School from Guangxi Provincial Department of Education (201203YB038); National Natural Scientific Foundation of P.R. China (No.31260069).
The authors declare that they have no competing interests.