Next Article in Journal
Ontogeny and Anatomy of the Dimorphic Pitchers of Nepenthes rafflesiana Jack
Previous Article in Journal
Verticillium Wilt of Mint in the United States of America
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

DNA-Based Authentication and Metabolomics Analysis of Medicinal Plants Samples by DNA Barcoding and Ultra-High-Performance Liquid Chromatography/Triple Quadrupole Mass Spectrometry (UHPLC-MS)

by
Marta Sánchez
,
Elena González-Burgos
,
Pradeep Kumar Divakar
and
M. Pilar Gómez-Serranillos
*
Department of Pharmacology, Pharmacognosy and Botany, Faculty of Pharmacy, Universidad Complutense de Madrid (UCM), 28040 Madrid, Spain
*
Author to whom correspondence should be addressed.
Plants 2020, 9(11), 1601; https://doi.org/10.3390/plants9111601
Submission received: 28 October 2020 / Revised: 14 November 2020 / Accepted: 17 November 2020 / Published: 18 November 2020

Abstract

:
There is growing interest for medicinal plants in the world drug market. Particularly, Matricaria recutita L., Valeriana officinalis L., Tilia spp., and Camellia sinensis (L.) Kuntze are some of the most consumed medicinal plants for treatment of minor health problems. Medicinal plants are seen as natural and safe; however, they can cause interactions and produce adverse reactions. Moreover, there is lack of consensus in medicinal plants regulation worldwide. DNA barcoding and UHPLC-MS technique are increasingly used to correctly identify medicinal plants and guarantee their quality and therapeutic safety. We analyzed 33 samples of valerian, linden, tea, and chamomile acquired in pharmacies, supermarkets, and herbal shops by DNA barcoding and UHPLC-MS. DNA barcoding, using matk as a barcode marker, revealed that CH1 sold as Camellia sinensis was Blepharocalyx tweediei, and sample TS2 sold as linden belong to Malvales. On the other hand, UHPLC-MS analysis revealed the presence of bioactive compounds (apigenin-7-glucoside, acetoxy valerenic acid, valerenic acid, epigallocatechin, and tiliroside). However, none of samples met minimum content of these active principles (except for valerenic acid in VF3) according to the European Medicines Agency (EMA) and Real Spanish Pharmacopeia. In conclusion, this study revealed the need to incorporate DNA barcoding and HPLC-MS techniques in quality controls of medicinal plants.

Graphical Abstract

1. Introduction

Medicinal plants constitute the basis of traditional and modern primary healthcare. Over 80% of the population, mainly of developing countries, depend on traditional and herbal medicine. Moreover, at least 25% of drugs in the modern pharmacopeia are derived from plants [1]. Furthermore, it is estimated that pharmacological activity has been evaluated in only 15% of all 300,000 plant species identified [2]. The consumption of medicinal plants for disease prevention and health promotion has increased significantly in the last two decades [3]. The reasons that explain the rise in therapeutic use of medicinal plants are several. They include natural tendency in population, erroneous perception of its safety, lower economic cost compared to conventional medicines, and polypharmacy [4,5]. The high demand for medicinal plants is reflected in the economic data of world market. The global trade in herbs was over USD 83 billion in 2012, being especially high in India, China, and Germany [6,7]. These medicinal plants are available in pharmacies, supermarkets, and herbal shops [4,8].
These medicinal plants must meet standards of quality, safety, and efficacy. In this context, pharmacopeia monographs include tests and acceptance criteria ranging from botanical identification, pharmacognostic evaluation, and chemical characterization with chromatographic methods to evaluate the quality of herbal medicines [9]. However, one of the biggest difficulties related to quality is that commercial medicinal plants are crushed or powdered, being problematic to identify phenotype or part of the plant [10]. The DNA barcoding technique is postulated as an effective tool to overcome limitations in the quality controls of commercial medicinal plants. DNA-based techniques consist of using short DNA sequences from standardized gene regions (rbcL, matK, and ITS2) [11,12,13]. One of the great advantages of DNA barcoding technique is that the result is not influenced by harvesting period, growth condition, environmental factors, and sample age, among other factors [14,15,16]. Moreover, high-performance liquid chromatography (HPLC) and mass chromatography (MS) techniques are widely useful to identify and quantify bioactive compounds found in medicinal plants [14]. Therefore, the correct botanical identification by DNA barcoding and the precise bioactive compounds determination by HPLC-MS constitute an integrated approach to guarantee quality and safety of market medicinal plants [17,18,19].
Matricaria recutita L. (chamomile), Valeriana officinalis L. (valerian), Tilia spp. (linden), and Camellia sinensis (L.) Kuntze (tea) are found among the most consumed medicinal plants. They are commonly acquired in pharmacies, herbal shops, and supermarkets. Chamomile flowers as infusions have a beneficial effect on digestion; valerian root in capsules are used for reducing anxiety and a nervous state and improving sleep; linden leaves as infusions reduce anxiety symptoms; tea leaves as capsule or infusion bring relief from digestive problems [20].
The aim of the present work is to apply DNA barcoding and UHPLC-MS methods as a tool to evaluate the quality of market samples of Matricaria recutita, Valeriana officinalis, Tilia spp., and Camellia sinensis.

2. Results and Discussion

2.1. DNA Barcoding Analysis

Commercial herbal products can be adulterer, replaced, or suffer some kind of contamination [20,21]. Herbal products are sold as ground or powdered form of a raw herb, which makes correct botanical identification difficult [22]. DNA barcoding constitutes a very useful tool for quality control and, consequently, for clinical safety [23,24]. Molecular analyses is crucial for accurate and fast identification of medicinal plants, since the plant fragments sold in the market is difficult to identify using traditional methods especially due to lack of morphological features. Moreover, they often mixed with other plant materials, and in such cases, molecular analysis is one of the best approaches for accurate sample identification [25,26].
In this study, the DNA barcode marker matk was used to identify four of the most consumed medicinal plants, which were M. recutita, V. officinalis, Tilia spp., and C. sinensis [27,28,29]. A total of 33 market samples (pharmacies, herbal shops, and supermarkets) were investigated (nine samples from M. recutita flowers, seven from V. officinalis root, nine from Tilia spp. Leaves, and eight from C. sinensis leaves). DNA was successfully extracted from 23 of 33 markets samples. Raw material conditions (i.e., very dry, storage quality) could explain why DNA could not be extracted in 10 of the samples. A total of eighteen new matk sequences were generated for 23 samples isolated from chamomile, linden, tea, and valeria acquired in different distribution channels. These matk sequences were aligned with 38 sequences downloaded from GenBank. A total of 10 sequences were in C. sinensis data matrix, 20 in M. recutita, nine in V. officinalis, and 17 in Tilia spp. The length of all newly generated matk sequences was above 700 bp. All new sequences generated for this study have been deposited in GenBank.
The results of the molecular phylogenetic and DNA barcoding analyses were largely congruent. In the phylogenetic analysis, all samples were grouped together forming monophyletic clades except for CH1, TS2, and VH2 samples (Figures S1–S4). The authenticity of these samples was confirmed with DNA barcoding analysis using the Barcode of Life Data (BOLD) system. The results of the BOLD system blast are depicted in Figures S5–S8. The sample codes MH1 to MH3, MF2, and MS1 to MS3 were identified as M. recutita; CH2 and CF3 as C. sinensis; VH1, VS1, and VF2 as Valeriana hirtella Kunth.; TH1 to TH3 resulted Tilia cordata Mill. (Table 1). The sample CH1 sold as C. sinensis did not correspond to this species but to Blepharocalyx tweediei [30] Berg. Blepharocalyx tweediei is a tree native to Argentina and Uruguay that is traditionally used as infusion for cough, bronchospasm, diarrhea, and other intestinal disorders. Further, TS2 sample sold as Linden matched with Malvales; no closely related species was detected. Remarkably, samples sold as V. officinalis and Tilia platyphyllos/europea Scop. were identified as Valeriana hirtella Kunth. and T. cordata, respectively. On the other hand, all samples sold as chamomile were correctly identified as Matricaria recutita. Therefore, our study has identified alterations in Camellia sinensis and Tilia spp. in one of the three and one of the four samples, respectively, of these two genera belonged to entirely different species/order. Figure 1 shows the DNA final barcode identification of the analyzed marked medicinal plants. A percentage of 36.4 of plant samples were identified as right species, 48.5% were incomplete samples, and 15.2% were consider as replacements. This work reveals the need for a correct botanical identification for valerian and linden to ensure the precise labeling. DNA-based sample authentication is necessary for market medicinal plants for correct identification [13,31,32].

2.2. UHPLC/MS Analysis

UHPLC-MS based metabolomics approach allows to qualitatively and quantitatively analyze all metabolites in medicinal plant species with high sensitivity and precision [33]. The first of these techniques (molecular assays) can only authenticate the medicinal plant, while the second assesses (chromatographic assays) its quality and can provide information on the presence and concentration of compounds with pharmacological activity, but can never determine the identity of the species in question (providing no direct evidence of fraud). The simultaneous use of both techniques is therefore an additional advantage for the evaluation of the quality of medicinal plants and thus for their efficacy and safety. The 33 market samples of the medicinal plants M. recutita, V. officinalis, Tilia spp. and C. sinensis were analyzed by HPLC-MS to identify and quantify the main bioactive compounds responsible for their pharmacological activity [34] (Table 2, Table 3, Table 4 and Table 5, Figure 2). Particularly, the compound apigenin-7-glucoside was identified in M. recutita samples. The concentration of apigenin-7-glucoside ranged from 0.001 to 0.035%. Most of the samples acquired in the different commercial establishments had an average content of this active principle of 0.003%. Only two commercial samples of chamomile had a higher content in apigenin-7-glucoside, specifically MH1 sample acquired in an herbal shop (0.035%) and MS3 sample acquired in a supermarket (0.016%). The European Medicines Agency and the Real Spanish Pharmacopeia establish that apigenin-7-glucoside content should be at least 0.25% of dried drug [35,36], thus being the content in this bioactive compound lower in all analyzed samples. Moreover, previous studies have identified that the content for apigenin-7-glucoside in dry material of Matricaria chamomilla varied from 210 to 1110 mg/100 g using UPLC-UV method [37].
For V. officinalis samples, the bioactive compounds identified were acetoxy valerenic acid and valerenic acid. Acetoxy valerenic acid concentration ranged from 0.020 to 0.053%. The average content for this compound for all samples (except for VF3) was 0.025%. The content in acetoxy valerenic acid for samples VF3 acquired in pharmacy was 0.053%. On the other hand, a variable content of valerenic acid was identified among the different valerian samples (ranging from 0.048 to 0.167%). Hence, the content of this bioactive compound was similar in those samples from herbal shops (from 0.101 to 0.115%) and very different among pharmacy samples (0.048, 0.078, and 0.167%). On the other hand, the content of valerenic acid in the valerian sample from supermarket was 0.084%. According to the European Medicines Agency and the Real Spanish Pharmacopeia, the content in sesquiterpeneic acids should be not less than 0.17% expressed as valerenic acid [36,38]. The only sample that meets this requirement is VF3; this sample acquired in pharmacy could not be identified in the DNA barcoding study. Moreover, Navarrete et al. (2013) analyzed valerenic acid content in valerian species using liquid chromatography with ultraviolet detection. Particularly, the content of valerenic acid was 0.88% [39].
The compound epigallocatechin was identified and quantified in C. sinensis samples. Its content varied between 1.51 and 4.72%. The lowest content corresponded to a pharmacy sample and the highest content to an herbal shop sample. The most common average content identified in four of six samples from pharmacies and supermarkets was 2.39%. On the other hand, the epigallocatechin content of samples from supermarkets was very similar (3.29 and 3.71%). The content in flavan-3-ols including epigallocatechin should range from 10–25% of dried drug according to European Medicines Agency [40,41]. Previous works reported that the content in epigallocatechin was 4.62% in Camellia sinensis samples [42].
Finally, tiliroside was identified in tea samples with a content variation between 0.008 and 0.043%. The low content of 0.008% was quantified in an herbal shop sample. The rest of linden samples had an average of 0.034%. There are no data on tiliroside content neither in the European Medicines Agency nor in the Real Spanish Pharmacopeia [43]. The content of tiliroside has been identified to be higher in inflorescences that in leaves in Tilia cordata (49.2 µg/g versus 16.1 µg/g) [44].
The analysis of the secondary metabolites revealed that all the plant species analyzed, regardless of whether or not they had been correctly identified botanically, contained to a greater or lesser extent the bioactive compounds that were sought. This shows that these compounds are presented in other species of the same genus (i.e., Valeriana hirtella) and of other orders (i.e., Malvales). Specifically, the tiliroside has been identified in many different families including Malvaceae and Tiliaceae [44].

3. Materials and Methods

3.1. Reagents

Methanol HPLC grade was obtained from Panreac Química (Barcelona, Spain). Acetonitrile HPLC grade, formic acid HPLC grade, SYBR safe DNA gel stain, and the primers MatK-1RKIM-f and MatK-3FKIM-r were from Thermo Fisher Scientific (Runcorn, Cheshire, UK). Agarose MB 250 was purchased from Biotools Biotechnological and Medical Laboratories SA (Madrid, Spain). Apigenin-7-glucoside, epigallocatechin, and tiliroside were acquired from Extrasynthese (Genay Cedex, France). Acetoxy acid was from Sigma-Aldrich (St. Louis, MO, USA), and valerenic acid was from Chromadex (Irvine, CA, USA).

3.2. Herbal Products

The 33 samples of the four medicinal plants Matricaria recutita L., Valeriana officinalis L., Tilia spp., and Camellia sinensis (L.) Kuntze were acquired from supermarkets, herbal shops, and pharmacies located in the Autonomous Community of Madrid (Spain). Particularly, nine of these samples were from M. recutita flowers (3 samples from pharmacies, 3 samples from herbal shops, and 3 samples from supermarkets), seven from V. officinalis root (3 samples from pharmacies, 3 samples from herbal shops, and 1 sample from supermarkets), nine from Tilia spp. leaves (3 samples from pharmacies, 3 samples from herbal shops, and 3 samples from supermarkets), and eight from C. sinensis leaves (3 samples from pharmacies, 3 samples from herbal shops, and 2 samples from supermarkets) (Table 6). They were stored in conditions of temperature and ambient humidity.

3.3. DNA Barcoding Analysis

3.3.1. DNA Extraction

DNA was extracted from market medicinal plants samples of roots (V. officinalis), leaves (Tilia spp. and C. sinensis), and flowers (M. recutita) using the Speed Tools tissue DNA extraction kit Biotools Biotechnological and Medical Laboratories following manufacturer’s instructions. Approximately, 100 mg of samples were pulverized with a sterile mortar in liquid nitrogen at room temperature and secondary metabolites were eliminated with methanol before starting DNA extraction. Then, samples were soaked in methanol for 2 h to remove potential secondary metabolites and dried overnight. Later, samples were incubated into lysis buffer initially at 65 °C for 30 min. In addition, root material samples were kept at room temperature overnight. After DNA extraction, samples were revealed in a 1% agarose gel stained with SYBR safe to check DNA quality [18].

3.3.2. PCR and Sequencing

Matk was chosen in this study, since it is one of the universal DNA barcode markers for land plants [25]. The matK gene of chloroplast is 1500 bp long, located within the intron of the trnK. Since this gene is larger in length, a fragment of this is used for DNA barcoding analysis. Further, the gene contains high substitution rates within the species and is a potential candidate for DNA barcoding studies [11]. PCR amplifications of matk were performed using specific primers MatK-1RKIM-f and MatK-3FKIM-r (Ki-Joong Kim, pers. comm.). The reaction mixture (25 µL final volume) contained 5 μL of DNA (1:10), 4.5 μL sterile water, 12.5 μL REDTaq ReadyMix PCR Reaction M, and 1.5 μL of each primer (forward and reverse) at 10 μM. PCR amplifications were carried out in a Techne R TC-3000 thermal cycler with the following conditions: one initial heating step of 45 s at 98 °C, followed by 35 cycles of 10 s at 98 °C, 30 s at 54 °C, and 40 s at 72 °C. A final extension step of 10 min at 72 °C was added, after which samples were kept at 4 °C.
Once PCR amplification is completed, a 1% agarose gel stained with SYBR safe was used to visualize DNA. PCR products were purified by Speed-Tools PCR Clean Up (Biotools Biotechnological and Medical Laboratories SA, Madrid, Spain) kit following the manufacturer’s instructions, and sequencing was performed with labelling using BigDye Terminator v. 3.1 Kit (Applied Biosystems, Madrid, Spain) as follows: 35 cycles of 20 sec at 94 °C, 20 sec at 48 °C, and 4 min at 60 °C. Sequences were obtained in an ABI PRISM 3130 Genetic Analyzer (Life Technologies, Alcobendas, Madrid, Spain). This is a standard approach to analyze PCR products.

3.3.3. Data Analysis

First, DNA sequences were assembled and manually adjusted in BioEdit sequence alignment editor software (v 7.2). Then, a second edition and assembly of the sequence fragments was made with the program SeqMan v.7 (Lasergene R, DNASTAR, Madison, WI, USA). Sequence identity was assessed using the mega-BLAST search function in GenBank [45,46]. Sequences with equal and above 95% similarity were downloaded from the GenBank and aligned with the newly generated sequences for this study. Separate dataset for each taxa viz. M. recutita, V. officinalis, Tilia spp., and C. sinensis were prepared for phylogenetic analysis. Each dataset was aligned using MAFFT v.7 [47] implementing the G-INS-I alignment algorithm, “1PAM/K = 2” scoring matrix with an offset value of 0.0, and the remaining parameters set to default values. The ambiguous regions in the matk alignment were assessed and removed using the least stringent option in Gblocks v.0.91b [48]. The alignments were analyzed using the maximum likelihood approach with RAxML v.8.2.6 program [49] as implemented on the CIPRES Web Portal, with the GTRGAMMA model. Nodal support was evaluated using 1000 bootstrap pseudoreplicates. The maximum likelihood (ML) trees were mid-point rooted. Phylogenetic trees were drawn using FigTree v.1.4.2 [50]. Moreover, sample identification was also assessed using a genetic distance-based blast option on the Barcode of Life Data System (BOLD Systems v3) [51]. The default setting was used for BLAST algorithm of the standard BOLD identification engine for matK sequences.

3.4. UHPLC-MS/MS Analysis

Samples were previously pulverized with a mortar and pestle sterile. Thirty milligrams of each medicinal plant samples were dissolved in 1 mL of ethanol/water 70/30% (v/v). Then, chamomile, linden, and tea samples were diluted 1:5 or 1:500, and valerian samples were diluted 1:4. HPLC standards were prepared at a concentration of 20 mg/L in methanol HPLC grade. Dilutions were prepared in a range of 0.05 to 1 mg/L in ethanol/water 70/30% (v/v) [52,53].
An ultra-high-performance liquid chromatography coupled with triple quadrupole mass spectrometry technique (UHPLC-QqQ-MS/MS) was developed using a LC-QQQ 8030 equipment (Shimazdu, Tokyo, Japan). The column was Phenomenex Gemini 5u C18 110A, 150 × 2 mm (Phenomenex, Alcobendas, Spain). The gradient mode was 7 min 5–95% Phase B; 8 min 95% Phase B; 8.5 min 5% Phase B using acetonitrile; in Phase A 0.1% formic acid in water. The flow rate was 0.5 mL/min and the infection volume was 10 µl for all medicinal plant samples, except for valerian, which was an injection volume of 20 μL.
Regarding LC-MS analysis, the mass spectrometer electrospray capillary voltage was maintained at 4.0 kV and the drying gas temperature at 250 °C with a flow rate of 15 mL/min. Nebulizer working flow was set at 1.5 mL/min. Nitrogen was used as both nebulizing and drying gas. Detection was carried out in Multiple reaction monitoring (MRM) mode, with a dwell time of 100 ms by monitoring three selective transitions for each parent compound. MRM transitions and their collision energy (CE) are shown in Table 7. The sample injection volume was 10 μL.
Data are expressed as means ± standard deviations of triplicate independent analyses.

4. Conclusions

DNA barcoding is revealed as an effective and necessary tool in the identification and authentication of plant species that are used in therapeutics. Since DNA barcoding does not provide qualitative or quantitative information on the metabolites of the plant raw material, its use together with chromatographic techniques such as HPLC-MS allows us to determine the precise identification of the species and the metabolic profile with great sensitivity and precision. The HPLC-MS technique combines the high resolution of ultra-high-performance liquid chromatography with high-resolution MS for molecular quantification. This technique has shown a high efficiency and resolution, a short time analysis, and increased sensitivity. The current work showed a high success rate in obtaining PCR amplification and sequencing of matk locus for the accurate sample identification of the market samples of medicinal plants. DNA barcoding demonstrated that the labeling of some medicinal plants acquired in different distribution channels is incorrect, demonstrating the need to apply DNA barcoding methods in the quality control of herbal products to ensure correct botanical identification and therefore ensure product quality. Moreover, the UHPLC-MS analysis displayed that all the analyzed samples presented the bioactive compounds (apigenin-7-glucoside, acetoxy valerenic acid, valerenic acid, epigallocatechin, and tiliroside) responsible for the pharmacological activity. This content was very variable for some plant species, as is the case of valerenic acid in valerian. In addition, none of the analyzed samples (except for valerenic acid in VF3) met the minimum content of these active principles according to the European Medicines Agency and the Real Spanish Pharmacopeia. This study shows the need to incorporate leading-edge molecular and analytical techniques for the quality control of plant species that are used in therapeutics to guarantee patient safety.
Future trends in the field of medicinal plants should be aimed at conducting more quality control studies in other widely consumed species and incorporating the combination of these molecular and analytical techniques in the plant industry for therapeutic purposes and even incorporating them into official documents (Pharmacopoeia, EMA).

Supplementary Materials

The following are available online at https://www.mdpi.com/2223-7747/9/11/1601/s1, Figure S1: Phylogenetic tress for M. recutita. Figure S2: Phylogenetic tress for V. officinalis. Figure S3: Phylogenetic tress for Tilia spp. Figure S4: Phylogenetic tress for C. sinensis. Figure S5: BOLD System for M. recutita. Figure S6: BOLD System for V. officinalis. Figure S7: BOLD System for Tilia spp. Figure S8: BOLD System for C. sinensis.

Author Contributions

Conceptualization, P.K.D. and M.P.G.-S.; methodology, P.K.D.; software, P.K.D.; validation, P.K.D. and M.S.; formal analysis, P.K.D. and E.G.-B.; investigation, M.S. and P.K.D.; writing—original draft preparation, P.K.D., M.S., and E.G.-B.; writing—review and editing, P.K.D. and E.G.-B.; supervision, M.P.G.-S. and P.K.D.; funding acquisition, P.K.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. WHO (World Health Organization). The World Traditional Medicines Situation. In Traditional Medicines: Global Situation, Issues and Challenges; WHO (World Health Organization): Geneva, Switzerland, 2011; pp. 1–14. [Google Scholar]
  2. De Luca, V.; Salim, V.; Atsumi, S.M.; Yu, F. Mining the biodiversity of plants: A revolution in the making. Science 2012, 336, 1658–1661. [Google Scholar] [CrossRef]
  3. Debbie, S.; Ladds, G.; Duez, P.; Williamson, E.; Chan, K. Pharmacovigilance of herbal medicine. J. Ethnopharmacol. 2012, 140, 513–518. [Google Scholar]
  4. Ekor, M. The growing use of herbal medicines: Issues relating to adverse reactions and challenges in monitoring safety. Front. Pharmacol. 2014, 4, 177. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Güney, O.I. Consumption attributes and preferences on medicinal and aromatic plants: A consumer segmentation analysis. Cienc. Rural 2019, 49, 5. [Google Scholar] [CrossRef]
  6. Sofowora, A.; Ogunbodede, E.; Onayade, A. The role and place of medicinal plants in the strategies for disease prevention. Afr. J. Tradit. Complement. Altern. Med. 2013, 10, 210–229. [Google Scholar] [CrossRef] [PubMed]
  7. Allkin, B. Useful Plants–Medicines: At Least 28,187 Plant Species are Currently Recorded as Being of Medicinal Use. In State of the World’s Plants; Willis, K.J., Ed.; Royal Botanic Gardens: London, UK, 2017. [Google Scholar]
  8. Wei, S.; Luo, Z.; Cui, S.; Qiao, J.; Zhang, Z.; Zhang, L.; Fu, J.; Ma, X. Molecular Identification and Targeted Quantitative Analysis of Medicinal Materials from Uncaria Species by DNA Barcoding and LC-MS/MS. Molecules 2019, 24, 175. [Google Scholar] [CrossRef] [Green Version]
  9. Neves, E.O.; de Sales, P.M.; Silveira, D. Pharmacopeial specifications and analytical data from post-marketing quality sampling and testing programs: A perspective beyond out-of-specification results. J. Pharm. Biomed. Anal. 2020, 178, 112935. [Google Scholar] [CrossRef]
  10. Palhares, R.M.; Gonçalves, M.; Dos Santos, B.; Pereira, G.; das Graças, M.; Oliveira, G. Medicinal plants recommended by the world health organization: DNA barcode identification associated with chemical analyses guarantees their quality. PLoS ONE 2015, 10, e0127866. [Google Scholar] [CrossRef] [Green Version]
  11. CBOL Plant Working Group. A DNA barcode for land plants. Proc. Natl. Acad. Sci. USA 2009, 106, 12794–12797. [Google Scholar] [CrossRef] [Green Version]
  12. Tungmunnithum, D.; Renouard, S.; Drouet, S.; Blondeau, J.P.; Hano, C. A Critical Cross-Species Comparison of Pollen from Nelumbo nucifera Gaertn. vs. Nymphaea lotus L. for Authentication of Thai Medicinal Herbal Tea. Plants 2020, 9, 921. [Google Scholar]
  13. Pawar, R.S.; Handy, S.M.; Cheng, R.; Shyong, N.; Grundel, E. Assessment of the Authenticity of Herbal Dietary Supplements: Comparison of Chemical and DNA Barcoding Methods. Planta Med. 2017, 83, 921–936. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Abubakar, M.B.; Salleh, M.F.; Shamsir Omar, M.S.; Wagiran, A. DNA Barcoding and Chromatography Fingerprints for the Authentication of Botanicals in Herbal Medicinal Products. Evid. Based Complement Altern. Med. 2017, 1352948. [Google Scholar] [CrossRef]
  15. Mili Bhargava, M.; Sharma, A. DNA barcoding in plants: Evolution and applications of in silico approaches and resources. Mol. Phylogenet. Evol. 2013, 67, 631–641. [Google Scholar] [CrossRef] [PubMed]
  16. Ichim, M.C. The DNA-Based Authentication of Commercial Herbal Products Reveals Their Globally Widespread Adulteration. Front. Pharm. 2019, 10, 1227. [Google Scholar] [CrossRef] [PubMed]
  17. Belščak-Cvitanović, A.; Valinger, D.; Benković, M.; Tušek, A.J.; Jurina, T.; Komes, D.; Kljusurić, J.G. Integrated approach for bioactive quality evaluation of medicinal plant extracts using HPLC-DAD, spectrophotometric, near infrared spectroscopy and chemometric techniques. Int. J. Food Prop. 2017, 20, 2463–2480. [Google Scholar] [CrossRef]
  18. Zuo, Y.; Chen, Z.; Kondo, K.; Funamoto, T.; Wen, J.; Zhou, S. DNA barcoding of Panax species. Planta Med. 2011, 77, 182–187. [Google Scholar] [CrossRef]
  19. Jianping, H.; Xiaohui Pang, X.; Baosheng, L.; Hui, Y.; Jingyuan, S.; Shilin, C. An authenticity survey of herbal medicines from markets in China using DNA barcoding. Sci. Rep. 2016, 6, 18723. [Google Scholar]
  20. Sánchez, M.; González-Burgos, E.; Iglesias, I.; Lozano, R.; Gómez-Serranillos, M.P. Current uses and knowledge of medicinal plants in the Autonomous Community of Madrid (Spain): A descriptive cross-sectional study. BMC Complement. Med. Ther. 2020, 20, 306. [Google Scholar] [CrossRef]
  21. Newmaster, S.G.; Grguric, M.; Shanmughanandhan, D.; Ramalingam, S.; Ragupathy, S. DNA barcoding detects contamination and substitution in North American herbal products. BMC Med. 2013, 11, 222. [Google Scholar] [CrossRef] [Green Version]
  22. Wachtel-Galor, S.; Benzie, I.F.F. Herbal Medicine: An Introduction to Its History, Usage, Regulation, Current Trends, and Research Needs. In Herbal Medicine: Biomolecular and Clinical Aspects, 2nd ed.; Benzie, I.F.F., Wachtel-Galor, S., Eds.; CRC Press/Taylor & Francis: Boca Raton, FL, USA, 2011; Chapter 1. [Google Scholar]
  23. Howard, C.; Lockie-Williams, C.; Slater, A. Applied Barcoding: The Practicalities of DNA Testing for Herbals. Plants 2020, 9, 1150. [Google Scholar] [CrossRef]
  24. Upton, R.; David, B.; Gafner, S.; Glasl, S. Botanical ingredient identification and quality assessment: Strengths and limitations of analytical techniques. Phytochem. Rev. 2019, 19, 1157–1177. [Google Scholar] [CrossRef] [Green Version]
  25. Hausmann, A.; Haszprunar, G.; Hebert, P.D.N. DNA Barcoding the Geometrid Fauna of Bavaria (Lepidoptera): Successes, Surprises, and Questions. PLoS ONE 2011, 6, e17134. [Google Scholar] [CrossRef] [PubMed]
  26. Janjua, S.; Fakhar-i-Abbas; William, K. DNA Mini-barcoding for wildlife trade control: A case study on identification of highly processed animal materials. Mitochondrial DNA 2016, 27, 544–546. [Google Scholar] [CrossRef] [PubMed]
  27. Ford, C.S.; Ayres, K.L.; Haider, N.; Toomey, N.; van-Alpen-Stahl, J.; Kelly, L.J.; Wikström, N.; Hollingsworth, P.M.; Duff, R.J.; Hoot, S.B.; et al. Selection of candidate DNA barcoding regions for use on land plants. Bot. J. Linn. Soc. 2009, 159, 1–11. [Google Scholar] [CrossRef] [Green Version]
  28. Raclariu, A.C.; Paltinean, R.; Vlase, L.; Labarre, A.; Manzanilla, V.; Ichim, M.C.; Crisan, G.; Brysting, A.K.; Boer, H. Comparative authentication of Hypericum perforatum herbal products using DNA metabarcoding, TLC and HPLC-MS. Sci. Rep. 2017, 7, 1291. [Google Scholar] [CrossRef] [PubMed]
  29. Raclariu, A.C.; Ţebrencu, C.E.; Ichim, M.C.; Ciupercǎ, O.T.; Brysting, A.K.; Boer, H. What’s in the box? Authentication of Echinacea herbal products using DNA metabarcoding and HPTLC. Phytomedicine 2018, 15, 32–38. [Google Scholar] [CrossRef]
  30. Hernández, J.; Raggone, M.; Bonazzola, P.; Bandoni, A.; Consolini, A. Antitussive, antispasmodic, bronchodilating and cardiac inotropic effects of the essential oil from Blepharocalyx salicifolius leaves. J. Ethnopharmacol. 2018, 210, 107–117. [Google Scholar] [CrossRef] [Green Version]
  31. Ivanova, N.V.; Kuzmina, M.L.; Braukmann, T.W.A.; Borisenko, A.V.; Zakharov, E.V. Authentication of Herbal Supplements Using Next- Generation Sequencing. PLoS ONE 2016, 11, e0156426. [Google Scholar] [CrossRef]
  32. Techen, N.; Parveen, I.; Pan, Z.; Khan, I.A. DNA barcoding of medicinal plant material for identification. Curr. Opin. Biotechnol. 2014, 25, 103–110. [Google Scholar] [CrossRef]
  33. Li, X.; Zhang, X.; Ye, L.; Kang, Z.; Jia, D.; Yang, L.; Zhang, B. LC-MS-Based Metabolomic Approach Revealed the Significantly Different Metabolic Profiles of Five Commercial Truffle Species. Front. Microbiol. 2019, 10, 2227. [Google Scholar] [CrossRef]
  34. Saleem, H.; Htar, T.T.; Naidu, R.; Anwar, S.; Zengin, G.; Locatelli, M.; Ahemad, N. HPLC-PDA Polyphenolic Quantification, UHPLC-MS Secondary Metabolite Composition, and In Vitro Enzyme Inhibition Potential of Bougainvillea glabra. Plants 2020, 9, 388. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. European Medicines Agency (EMA). Available online: https://www.ema.europa.eu/en/documents/herbal-report/final-assessment-report-matricaria-recutita-l-flos-matricaria-recutita-l-aetheroleum_en.pdf (accessed on 30 September 2020).
  36. Real Farmacopea Española. Available online: https://extranet.boe.es/farmacopea/ (accessed on 1 October 2020).
  37. Haghi, G.; Hatami, A.; Safaei, A.; Mehran, M. Analysis of phenolic compounds in Matricaria chamomilla and its extracts by UPLC-UV. RPS 2014, 9, 31–37. [Google Scholar] [PubMed]
  38. European Medicines Agency (EMA). Available online: https://www.ema.europa.eu/en/documents/herbal-report/final-assessment-report-valeriana-officinalis-l-radix-valeriana-officinalis-l-aetheroleum_en.pdf (accessed on 30 September 2020).
  39. Navarrete, A.; Avula, B.; Choi, Y.W.; Khan, I.A. Chemical fingerprinting of valeriana species: Simultaneous determination of valerenic acids, flavonoids, and phenylpropanoids using liquid chromatography with ultraviolet detection. J. AOAC Int. 2006, 89, 8–15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. European Medicines Agency (EMA). Available online: https://www.ema.europa.eu/en/documents/herbal-report/draft-assessment-report-camellia-sinensis-l-kuntze-non-fermentatum-folium_en.pdf (accessed on 1 October 2020).
  41. European Medicines Agency. Available online: https://www.ema.europa.eu/en/documents/scientific-guideline/guideline-specifications-test-procedures-acceptance-criteria-herbal-substances-herbal-preparations/traditional-herbal-medicinal-products-revision-1_en.pdf (accessed on 25 September 2020).
  42. Raghavendra, H.L.; Kumar, S.V.P.; Kekuda, T.R.P.; Ejeta, E.; Molla, B.; Anilakumar, K.R.; Khanum, F. HPLC method for chemical composition and in vitro antioxidant activity of Camellia sinensis Linn. Anal. Chem. Lett. 2011, 1, 361–369. [Google Scholar] [CrossRef]
  43. Pieczykolan, A.; Pietrzak, W.; Nowak, R.; Pielczyk, J.; Lamacz, K. Optimization of Extraction Conditions for Determination of Tiliroside in Tilia L. Flowers Using an LC-ESI-MS/MS Method. J. Anal. Methods Chem. 2019, 9052425. [Google Scholar] [CrossRef] [Green Version]
  44. Nowak, R. Separation and quantification of tiliroside from plant extracts by SPE/RP-HPLC. Pharm. Biol. 2003, 41, 627–630. [Google Scholar] [CrossRef]
  45. Sayers, E.W.; Barrett, T.; Benson, D.A.; Bolton, E.; Bryant, S.H.; Canese, K.; Chetvernin, V.; Church, D.M.; Dicuccio, M.; Federhen, S.; et al. Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2011, 39, D38–D51. [Google Scholar] [CrossRef] [Green Version]
  46. Pedales, R.D.; Damatac, A.M.; Limbo, C.A.; Marquez, C.M.; Navarro, A.I.B.; Molina, J. DNA Barcoding of Philippine Herbal Medicinal Products. J. AOAC Int. 2016, 99, 1479–1489. [Google Scholar] [CrossRef]
  47. Katoh, K.; Toh, H. Recent developments in the MAFFT multiple sequence alignment program. Brief. Bioinform. 2008, 9, 286–298. [Google Scholar] [CrossRef] [Green Version]
  48. Talavera, G.; Castresana, J. Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst. Biol. 2007, 56, 564–577. [Google Scholar] [CrossRef] [Green Version]
  49. Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 2014, 30, 1312–1313. [Google Scholar] [CrossRef] [PubMed]
  50. Rambaut, A. “FigTree 1.2.2.”. Available online: http://tree.bio.ed.ac.uk/software/figtree/ (accessed on 13 August 2020).
  51. Ratnasingham, S. Hebert PDN. BOLD: The Barcode of Life Data System. Mol. Ecol. Notes 2007, 7, 355–364. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Wang, C.; Zheng, Z.; Deng, X.; Ma, X.; Wang, S.; Liu, J.; Liu, Y.; Shi, J. Flexible and powerful strategy for qualitative and quantitative analysis of valepotriates in Valeriana jatamansi Jones using high-performance liquid chromatography with linear ion trap Orbitrap mass spectrometry. J. Sep. Sci. 2017, 40, 1906–1919. [Google Scholar] [CrossRef] [PubMed]
  53. Scoparo, C.T.; de Souza, L.M.; Dartora, N.; Sassaki, G.L.; Gorin, P.A.; Iacomini, M. Analysis of Camellia sinensis green and black teas via ultra high performance liquid chromatography assisted by liquid-liquid partition and two-dimensional liquid chromatography (size exclusion × reversed phase). J. Chromatogr. A 2012, 1222, 29–37. [Google Scholar] [CrossRef] [PubMed]
Figure 1. DNA-sequence-based identification of the analyzed medicinal plants marked samples.
Figure 1. DNA-sequence-based identification of the analyzed medicinal plants marked samples.
Plants 09 01601 g001
Figure 2. Representative UHPLC-ESI-QqQ-MS/MS chromatograms for (A) acetoxyvalerenic and valerenic acids presented in valerian samples, (B) epigallocatechin presented in tea samples, (C) apigenin presented in chamomile, and (D) tiliroside presented in linden.
Figure 2. Representative UHPLC-ESI-QqQ-MS/MS chromatograms for (A) acetoxyvalerenic and valerenic acids presented in valerian samples, (B) epigallocatechin presented in tea samples, (C) apigenin presented in chamomile, and (D) tiliroside presented in linden.
Plants 09 01601 g002aPlants 09 01601 g002b
Table 1. Sample identification using the Barcode of Life Data (BOLD) systems.
Table 1. Sample identification using the Barcode of Life Data (BOLD) systems.
Sample CodeLabelledBOLD System Identification
MH1MatricariaMatricaria recutita
MH2ChamomileMatricaria recutita
MH3Matricaria chamomillaMatricaria recutita
MF2ChamomileMatricaria recutita
MS1ChamomileMatricaria recutita
MS2Matricaria recutitaMatricaria recutita
MS3ChamomileMatricaria recutita
CH1Camellia sinensisBlepharocalyx tweediei
CH2TheaCamellia sinensis
CF3TheaCamellia sinensis
VH1Valeriana officinalisValeriana hirtella
VS1Valeriana officinalisValeriana hirtella
VF2Valeriana officinalisValeriana hirtella
TH1Tilia platyphyllosTilia cordata
TH2LindenTilia cordata
TH3Tilia EuropeaTilia cordata
TS2LindenMalvales
Table 2. Content of apigenin-7-glucoside in Matricaria recutita in samples labeled as acquired in pharmacies, herbal shops, and pharmacies. Data are expressed as means ± standard deviations of triplicate independent analyses.
Table 2. Content of apigenin-7-glucoside in Matricaria recutita in samples labeled as acquired in pharmacies, herbal shops, and pharmacies. Data are expressed as means ± standard deviations of triplicate independent analyses.
Matricaria recutita L.
SampleApigenin-7-glucoside (mg/g) (mean ± SD)
MH10.35 ± 0.015
MH20.02 ± 0.004
MH30.05 ± 0.005
MS10.07 ± 0.004
MS20.06 ± 0.003
MS30.16 ± 0.009
MF10.03 ± 0.002
MF20.01 ± 0.004
MF30.03 ± 0.005
Table 3. Content of valerenic acid and acetoxyvalerenic acid in Valeriana in samples labeled as acquired in pharmacies, herbal shops, and pharmacies. Data are expressed as means ± standard deviations of triplicate independent analyses.
Table 3. Content of valerenic acid and acetoxyvalerenic acid in Valeriana in samples labeled as acquired in pharmacies, herbal shops, and pharmacies. Data are expressed as means ± standard deviations of triplicate independent analyses.
Valeriana Officinalis L.
SampleAcetoxyvalerenic Acid (mg/g) (mean ± SD)Valerenic Acid (mg/g) (mean ± SD)
VH10.26 ± 0.030 1.09 ± 0.017
VH20.24 ± 0.0111.15 ± 0.017
VH30.24 ± 0.0141.01 ± 0.020
VF10.20 ± 0.0200.48 ± 0.015
VF20.29 ± 0.0170.78 ± 0.026
VF30.53 ± 0.0141.67 ± 0.030
VS10.25 ± 0.0550.84 ± 0.015
Table 4. Content of epigallocatechin in Camellia sinensis in samples labeled as acquired in pharmacies, herbal shops, and pharmacies. Data are expressed as means ± standard deviations of triplicate independent analyses.
Table 4. Content of epigallocatechin in Camellia sinensis in samples labeled as acquired in pharmacies, herbal shops, and pharmacies. Data are expressed as means ± standard deviations of triplicate independent analyses.
Camellia Sinensis (L.) Kuntze
SampleEpigallocatechin (mg/g) (mean ± SD)
CH121.2 ± 0.025
CH247.2 ± 0.015
CH323.2 ± 0.060
CS132.9 ± 0.025
CS237.1 ± 0.011
CF115.1 ± 0.020
CF225.9 ± 0.011
CF325.5 ± 0.030
Table 5. Content of tiliroside in Tilia in samples labeled as acquired in pharmacies, herbal shops, and pharmacies. Data are expressed as means ± standard deviations of triplicate independent analyses.
Table 5. Content of tiliroside in Tilia in samples labeled as acquired in pharmacies, herbal shops, and pharmacies. Data are expressed as means ± standard deviations of triplicate independent analyses.
Tilia spp.
SampleTiliroside (mg/g) (mean ± SD)
TH10.32 ± 0.056
TH20.39 ± 0.045
TH30.08 ± 0.020
TS10.27 ± 0.026
TS20.42 ± 0.032
TS30.43 ± 0.025
TF10.23 ± 0.020
TF20.36 ± 0.010
TF30.35 ± 0.011
Table 6. Samples of the medicinal plants chamomile, valerian, tea, and linden acquired from supermarkets, herbal shops, and pharmacies located in the Autonomous Community of Madrid.
Table 6. Samples of the medicinal plants chamomile, valerian, tea, and linden acquired from supermarkets, herbal shops, and pharmacies located in the Autonomous Community of Madrid.
ChamomileValerian
PharmaciesHerbal ShopsSupermarketsPharmaciesHerbal ShopsSupermarkets
MF1MH1MS1VF1VH1VS1
MF2MH2MS2VF2VH2
MF3MH3MS3VF3VH3
LindenTea
PharmaciesHerbal ShopsSupermarketsPharmaciesHerbal ShopsSupermarkets
TF1TH1TS1CF1CH1CS1
TF2TH2TS2CF2CH2CS2
TF3TH3TS3CF3CH3
Table 7. Multiple reaction monitoring (MRM) transitions and t collision energy conditions.
Table 7. Multiple reaction monitoring (MRM) transitions and t collision energy conditions.
Active CompoundTransitionCEDwell (msec)
Acetoxyalerenic acid291.3 > 59.0027 V100
291.3 > 249.0517 V100
Valerenic acid233.2 > 40.9222 V100
233.2 > 84.2025 V100
Epigallocatechin457 > 169.0518 V100
457 > 125.0045 V100
Tiliroside592.8 > 284.9533 V100
592.8 > 254.9555V100
Apigenin-7-O-glucoside431.3 > 267.9035 V100
431.3 > 150.9552 V100
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sánchez, M.; González-Burgos, E.; Divakar, P.K.; Gómez-Serranillos, M.P. DNA-Based Authentication and Metabolomics Analysis of Medicinal Plants Samples by DNA Barcoding and Ultra-High-Performance Liquid Chromatography/Triple Quadrupole Mass Spectrometry (UHPLC-MS). Plants 2020, 9, 1601. https://doi.org/10.3390/plants9111601

AMA Style

Sánchez M, González-Burgos E, Divakar PK, Gómez-Serranillos MP. DNA-Based Authentication and Metabolomics Analysis of Medicinal Plants Samples by DNA Barcoding and Ultra-High-Performance Liquid Chromatography/Triple Quadrupole Mass Spectrometry (UHPLC-MS). Plants. 2020; 9(11):1601. https://doi.org/10.3390/plants9111601

Chicago/Turabian Style

Sánchez, Marta, Elena González-Burgos, Pradeep Kumar Divakar, and M. Pilar Gómez-Serranillos. 2020. "DNA-Based Authentication and Metabolomics Analysis of Medicinal Plants Samples by DNA Barcoding and Ultra-High-Performance Liquid Chromatography/Triple Quadrupole Mass Spectrometry (UHPLC-MS)" Plants 9, no. 11: 1601. https://doi.org/10.3390/plants9111601

APA Style

Sánchez, M., González-Burgos, E., Divakar, P. K., & Gómez-Serranillos, M. P. (2020). DNA-Based Authentication and Metabolomics Analysis of Medicinal Plants Samples by DNA Barcoding and Ultra-High-Performance Liquid Chromatography/Triple Quadrupole Mass Spectrometry (UHPLC-MS). Plants, 9(11), 1601. https://doi.org/10.3390/plants9111601

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop