Skip to main content

Authentication of milk thistle commercial products using UHPLC-QTOF-ESI + MS metabolomics and DNA metabarcoding



Milk thistle is one of the most popular hepatoprotectants, and is often sold in combination with other ingredients. Botanical supplements are known to be vulnerable to contamination and adulteration, and emerging technologies show promise to improve their quality control.


Untargeted and semi-targeted metabolomics based on UHPLC-QTOF-ESI+MS techniques, UV spectrometry, and DNA metabarcoding using Illumina MiSeq were used to authenticate eighteen milk thistle botanical formulations (teas, capsules, tablets, emulsion).


Untargeted metabolomics separated 217 molecules and by multivariate analysis the discrimination between the different preparations was established. The semi-targeted metabolomics focused on 63 phytochemicals, mainly silymarin flavonolignans and flavonoids, that may be considered as putative biomarkers of authenticity. All formulations contained molecules from silymarin complexes at different levels. The quantitative evaluation of silybins was done using in parallel UV spectrometry and UHPLC-QTOF-ESI+MS and their correlations were compared. DNA metabarcoding detected milk thistle in eleven out of sixteen retained preparations, whereas two others had incomplete evidence of milk thistle despite metabolomics validating specific metabolites, e.g., silymarin complex, identified and quantified in all samples. Meanwhile, the DNA metabarcoding provided insights into the total species composition allowing the interpretation of the results in a broad context.


Our study emphasizes that combining spectroscopic, chromatographic, and genetic techniques bring complementary information to guarantee the quality of the botanical formulations.

Peer Review reports


Silybum marianum (L.) Gaertn. (Asteraceae, milk thistle, MT, Fig. 1) preparations are among the most commonly used botanical-based hepatoprotectants in complementary and alternative medicine [1,2,3,4,5,6]. Milk thistle has been purported to have also other health-promoting effects being used for the treatment of dyspeptic complaints, alcohol or drug-induced hepatic cirrhosis and fibrosis, and support treatment in hepatitis and other chronic inflammatory liver conditions [2, 7, 8], for stimulation of milk production in lactating mothers [9,10,11], and has been investigated for oncological indications and metabolic syndrome [2, 12, 13]. In Europe, intravenous silibinin, a flavonolignan isolated from milk thistle, has been approved as an antidote in patients intoxicated with Amanita phalloides, a mushroom that causes fatal poisoning [14].

Fig. 1
figure 1

Silybum marianum (L.) Gaertn. (A) Red purple flower head with spiny bracts; (B) Variegated leaf with lobed margins; (C) Mature flower head with seeds; (D) Fruits. (Photos: A.C. Raclariu-Manolică; M. Naie)

The therapeutic activity of milk thistle is associated to a great extent with a mixture of flavonolignans, known as silymarin complex including mainly silibinin (or silybin) A and B, isosilibinin (or isosilybin) A and B, silychristin and silydianin [15,16,17,18,19]. These major flavonolignans together with the flavonoid taxifolin are considered the marker compounds for milk thistle identification [20]. Silymarin is present in seeds, fruits, and leaves of milk thistle, but mature seeds are being reported as having the maximum concentration [19, 21]. Milk thistle contains also other flavonoids (e.g. kaempferol, quercetin, rutin, luteolin, naringin, kaempferol, apigenin), proteins, sugars (arabinose, rhamnose, xylose, glucose), tocopherol, sterols (cholesterol, campesterol, stigmasterol, sitosterol), and lipids in the form of triglycerides (linoleic, oleic and palmitic acids) [22]. The putative mechanisms of action of the main bioactive compounds of milk thistle have been discussed by numerous pharmacological, pharmacokinetic, and toxicological studies [23,24,25,26].

Regarded to be “natural” and thus “safe”, thousands of botanical preparations are advertised, marketed, and sold via various channels, and are often preferred over synthetic pharmaceuticals by consumers [27, 28]. The claimed therapeutic properties, the high popularity among consumers, folk traditions, and the increased global market demand, have spurred studies on various aspects of botanicals and their derived preparations [29,30,31]. Ensuring their quality and safety and reducing the potential risks related to their intake are key priorities for these commodities, with a paramount focus on consumer health [32,33,34,35,36]. These matters have raised the interest in finding novel testing and quality monitoring strategies, including emerging technologies applied for authentication purposes [36, 37]. But, finding a single comprehensive analytical approach for the authentication of botanicals and their derived preparation is a complicated task hampered by the complexity of these products and by the lack of harmonization regarding regulations, definitions, and quality standards that vary between countries and continents [38,39,40]. Botanicals are inherent chemical mixtures prone to variability under natural conditions that are often reflected in the batch-to-batch composition variation of the final preparations [41, 42]. Moreover, they often have long and complex supply chains, where numerous ingredients are extracted and processed differently, and key aspects such as identity and authenticity remain challenging to assess, hindering accurate monitoring and quality control processes [42, 43]. These are only a few aspects that make botanical preparations some of the most vulnerable commodities worldwide, especially to accidental contamination and fraud through adulteration [44].

Treatments involving milk thistle are generally well tolerated in recommended doses, with a low incidence of adverse drug reactions and mild side effects when reported [23, 45, 46]. However, incorrectly used nomenclature for milk thistle-based material has generated a strong degree of skepticism regarding the safe use and efficacy of this botanical and its derived preparations [47, 48]. Regarding this, strong arguments have been made concerning the description of milk thistle compounds, particularly “silymarin” and “silibinin”, terms that were reported to be often used interchangeably [48]. Furthermore, the proportion of these marker compounds was reported as being prone to variability under natural conditions and significantly affected by the production and processing steps [47, 49,50,51,52,53,54]. Nevertheless, the chemical composition of the milk thistle used material was rarely determined in most of the studies focusing on its biological activity [52]. Thus, notwithstanding the purported beneficial role of milk thistle, the vague description of the chemical composition of the derived preparations represents probably, for several studies, the main pitfall in proving its clinical efficacy [52, 55, 56].

Moreover, widespread contamination with fungi, microbes, and pesticides has been reported in milk thistle-based dietary supplements raising serious safety issues for human health, as botanicals-induced hepatotoxicity may occur [55, 57, 58]. Studies have shown alarming discrepancies between declared and detected chemical content between brands of marketed milk thistle preparations, as well as within batches of the same preparations and manufacturers. These differences can critically alter the expected therapeutic effects [51, 55, 56, 59, 60]. Some of the studies focusing on the quantitative analysis of silymarin showed that a large number of investigated marketed preparations contained a lower amount than declared and some were even completely missing this marker compound of milk thistle [51, 53, 59, 61]. On top of this, some of the studies provided evidence of the presence of foreign matters in the products, reported as probably belonging to undeclared adulterants, but the identity of these adulterants was not determined [60, 62].

Considering the reviewed challenges, one may wonder if the consumption of commercial milk thistle preparations may rather be harmful than beneficial with regard to possible accidental contaminants and/or adulterants. Pharmacovigilance of these preparations remains difficult since they are sold over-the-counter with no medical prescription, and limited legislative framework to trace or monitor adverse reactions [63, 64]. In addition, the standard quality control analytical methods do not always have sufficient resolution for the identification of target plant species within complex preparations, and often are not able to detect non-targeted plant ingredients that may be present as contaminants or adulterants [39, 65]. Thus, new technologies and fit-for-purpose methodologies need to be adopted for the quality control of botanicals and their derived complex preparations [6, 36].

In this study, we propose and evaluate a novel analytical approach to investigate the authenticity of commercial botanical preparations labeled as containing Silybum marianum (L.) Gaertn. (milk thistle), either as unique ingredient or in combination with other plant-based ingredients. Using untargeted and semi-targeted metabolomics based on ultra-high-performance liquid chromatography coupled with quadrupole-time of flight mass spectrometry (UHPLC-QTOF-ESI+MS) data and ultraviolet spectroscopy (UV), alongside high-throughput DNA metabarcoding, we aimed to answer the following research questions: (1) Can UHPLC-QTOF-ESI+MS untargeted metabolomics unveil potential molecular markers that differentiate between unique and multiple ingredient milk thistle-based preparations?; (2) Can semi-targeted metabolomics identify key molecules to assess the authenticity of milk thistle preparations, and to detect any deviation compared to other plant ingredients stated on the label of the commercial preparation?; (3) Can UV spectrometry be useful as a fast method comparable to UHPLC-QTOF-ESI+MS for authentication of milk-thistle in botanical formulations?; (4) Can DNA metabarcoding be used to test for the presence of milk thistle in botanical preparations, and to detect the presence of off-label plant species? Ultimately, this study aims to provide a new fit-for-purpose complementary analytical approach to assess the quality of complex botanicals and derived formulations, to enable more rapid advances in the regulatory context.


Botanical formulations and reference material

Eighteen herbal preparations that included Silybum marianum (L.) Gaertn and or other derived compounds (i.e., silymarin) according to the label were randomly purchased from Romania (15) and Germany (3) in the autumn of 2021. The samples were bought from herbal shops (8), via e-commerce (6), retail stores (3), and pharmacies (1), and were sold as herbal teas (7), tablets (6), capsules (4), and one emulsion (Tables 1 and Additional file 1.A and 1.B). According to the label information, there were 7 unique ingredients (U) and 11 multi-plant ingredient preparations (M), as presented in Table 1. These products for scientific analysis were imported into Norway under Norwegian Medicines Agency license no. 18/13493–2. Each sample was given a specific ID number ranging from PA1 to PA18 and a code referring to the type and pharmaceutical form of the preparations, as follows: “U” unique ingredient, “M” multiple ingredients, “T” teas, “Tb” tablets, “C” capsules, and “E” emulsions. An overview of the samples including label information, but not the producer/importer name, lot number, expiration date, or any other information that could lead to the identification of that specific product can be found in Additional file 1.A. The five genuine MT herbal materials used as references for the identification and quantification of the main target compounds in the metabolomics analysis were kindly provided by the Agricultural Research and Development Station Secuieni (Neamt County, RO), Vegetable Research and Development Station Bacău (Bacău County, RO), and a local farmer from Fundu Tutovei (Bacău County, RO), and they have ID collection codes ranging from ACM24 to 27, and ACM29, and the code “Genuine U0” (their description can be found in Additional file 1.C.). Ancuța Cristina Raclariu-Manolică undertook the formal identification of the plant material used as a reference in this study. Voucher specimens are deposited at the National Institute of Research and Development for Biological Sciences,” Stejarul” Biological Research Centre (Romania), having deposition numbers ranging from PlantCheck_ACM24 to PlantCheck_ ACM27, and PlantCheck_ACM29, and available on request.

Table 1 Categories of herbal formulations (teas, capsules, tablets and emulsion) submitted to different analysis, and their codes applied for unique (U) or multiple (M) ingredients. Abbreviations: C-capsules; T-teas; Tb-Tablets; E-emulsion

Extraction of phytochemicals

The same quantity of 5 g from each sample (5 genuine powdered MT and 18 herbal supplements) was suspended in 100 ml ethanol 70%, mixed 3 min by vortex and kept in an ultrasonic bath for 3 × 20 min at 50 °C. After the storage, 24 h at room temperature, each extract was centrifuged at 12,500 rpm and the supernatant was collected and filtered through 0.25 mm membrane filter. All extractions were made in triplicate.

UV spectroscopy

The UV spectra (200–340 nm) were recorded using a UV/VIS Lambda 25 (Perkin Elmer Inc, Waltham, Massachusetts, USA) spectrometer and the measurements were done in quartz cuvettes, comparative to a blank sample (ethanol 70%). Each extract was filtered through a 0.4 μm nylon membrane and diluted with ethanol 70% in different proportions to fit in the spectral absorbance scale. The specific absorbances located in the region 286–288 nm were recorded to evaluate the levels of flavonolignans (FL). In parallel, a calibration curve was built with pure silybins A + B (25 to 75 micrograms/ml) having the following equation: y = 0.0246x–0.2363 (R2 = 0.9968), as presented in Additional file 2.A. Considering the calibration curve, the results for each formulation were expressed in mg silybin equivalents per g dry matter (d.m.). This method offered a preliminary information and a rough evaluation of the silymarin flavonolignans found in the herbal preparations that claimed their presence on the label.

UHPLC-QTOF-ESI+MS metabolomics

Solvents, reagents, and analytical standards

HPLC grade pure ethanol, acetonitrile, and methanol were purchased from Merck (Darmstadt, Germany), and formic acid (99.99%) was purchased from Sigma-Aldrich (St. Louis, Missouri, United States). Deionized water was produced by a Milli-Q system (Millipore, Bedford, MA, USA). The analytical standard of Silybin (a mixture of silybin A and B) was purchased from Sigma Aldrich (CAS nr. 802918-57-6, St. Louis, Missouri, United States).

Untargeted and semi-targeted metabolomics

The metabolomic fingerprints of all ethanolic extracts were performed using the ultra-high-performance liquid chromatography coupled with electrospray -quadrupole-time of flight-mass spectrometry using the positive ionization (UHPLC-QTOF-ESI+MS) on a UltiMate 3000 UHPLC system equipped with a quaternary pump Dionex delivery system (Thermo Fisher Scientific Inc., Waltham, Massachusetts, USA), and mass spectrometry (MS) detection by a QqTOF MaXis Impact (Bruker Daltonics GmbH, Bremen, Germany). The metabolites were separated using a Kinetex column (Phenomenex Inc, Torrance, USA) (5 μm, 150 × 2.1 mm, 100 Å) at 25 °C. The flow rate was set at 0.8 ml·min− 1 and the volume of each injected extract was 8 µl. The mobile phase consisted of 0.1% formic acid in water (A) and 0.1% formic acid in acetonitrile (B). The gradient was: 20 to 40% B (0–5 min), 40–60% B (5–8 min), 60–70% B (8–10 min), 70–20% B (10–16 min), and 20% B isocratic until 24 min. Several quality control (QC) samples obtained from a pool of extracts were used in parallel to calibrate the separations. The chromatograms were processed using Chromeleon software (Dionex, Thermo Fisher Scientific Inc, Waltham, Massachusetts, USA). The MS parameters were: ionization ESI+, calibrated with sodium formate, capillary voltage 3500 V, nebulizing gas pressure of 2.8 bar, drying gas flow 12 l/min, drying temperature 300 °C. The control of the instrument and the data processing were done using the specific software TofControl 3.2, HyStar 3.2, Data Analysis 4.2 (Bruker Daltonics GmbH, Bremen, Germany).

Metabolomic data processing and statistical analysis

The Base Peak chromatograms and all MS spectra were recorded and processed by Compass DataAnalysis 4.2 (Bruker Daltonics, GmbH, Bremen, Germany) using the find molecular feature (FMF) algorithm. The time alignment, spectral background extraction, normalization by the median values, of the bucket values in analysis, and an 80% bucket filter were the used parameters. From the initial metabolic matrix, including retention time, MS peak intensity, signal/noise ratios (> 10), and mass-to-charge ratio (m/z) values of separated molecules, after retention of 60% common molecules, a total of 217 molecules, having m/z values from 270 to 615 Dalton, were selected.

The statistical analysis was done by the Metaboanalyst v5.0 online software [66] and algorithms ( From the matrices representing the MS peak intensity versus mass-to-charge ratio (m/z) values of each molecule from each sample the most relevant statistical parameters were tested to reflect the discrimination between sample groups, the prediction, and the correlation maps. Therefore, the Principal Component Analysis (PCA) and Sparse Partial Least Square Discriminant Analysis (sPLSDA), the Heatmaps and the Random Forest-based prediction were used to evaluate the similarities between samples, the identification of the putative biomarkers. According to the statistical analysis, molecules, which may explain the discriminations between samples and some specific putative biomarkers for authenticity, were selected and identified using international databases, e.g. Human Metabolome Database [67], Lipid Maps [68], Phenol-Explorer (version 3.6) [69], and PubChem [70].

As mentioned before, the untargeted analysis revealed 217 molecules to be considered for the discrimination between the different categories of formulations. In the second step, the semi-targeted analysis was performed using the same Metaboanalyst 5.0. software, One way ANOVA algorithm [66]. Sixty-three molecules were considered as potential authenticity markers, from phytochemical classes characteristic of MT seeds. Such semi-targeted analysis focused on silymarin complex, including taxifolin, but also were lignan precursors (coumaric acid and coniferyl derivatives), phytosterols, phenolic acids, flavonoids, fatty acids, and derivatives, as mentioned in Additional file 3.

For a quantitative evaluation of silybins in each formulation, using UHPLC-QTOF-ESI+MS technique, the calibration curve was also built using a stock solution of 4 mg/ml pure silybins A + B. Five different volumes of stock solution (from 1.25 to 7.5 µl were injected (corresponding to 5, 10, 15, 20, and 30 micrograms). The equation: y = 337240x + 108,540 (R2 = 0.9785) was considered to calculate the concentration of silybins in every sample, expressed in mg silybin equivalents per g dry matter curve (see Additional file 2.B).

DNA metabarcoding

DNA Extraction

Each sample was already made into powder, and total DNA from all samples was extracted from the homogenized contents using the E.Z.N.A.®SP plant DNA kit (Omega Biotek Inc, Norcross Georgia) following the manufacturer’s instructions. DNA extracts were then quantified using a Qubit 2.0 Fluorometer with dsDNA Broad-Range assay kit (Invitrogen, USA). In the case of non-successful DNA extraction using the E.Z.N.A. Plant DNA Mini Kit, subsamples were extracted following a modified CTAB extraction method as described by Doyle and Doyle [71], and adapted by Raclariu et al. [72]. The final elution volume was 100 µl.

DNA libraries preparation and sequencing

All amplicon libraries were prepared in three technical replicates on 96-well polymerase chain reaction (PCR) plates. On each plate, we also included negative controls consisting of extraction blanks (created by performing all steps of the DNA extraction on “empty” samples) alongside the DNA extraction of other materials, and PCR controls (created by replacing the template DNA with ddH2O at the PCR step). This resulted in a total of 21 negative controls across the project (i.e., six extraction blanks analyzed in triplicate and four PCR controls). The amplicon libraries for the nuclear ribosomal target sequences, internal transcribed spacer nrITS2, were performed using indexed ITS-3p62plF1 and ITS-4unR1 primers designed in [73] following the indexing strategy as in [74].

PCR was conducted with the following conditions: 1X Q5 hot start high fidelity mastermix (New England Biolabs Inc, UK), 1X Q5 enhancer (New England Biolabs Inc, UK), 0.5 µM of each indexed ITS-3p62plF1 and ITS-4unR1 primer [73] and 3 µl of extracted DNA in a final volume of 25 µl. Unique dual-index primer combinations were used for each subsample as in [75], and the thermocycling protocol was as described in [73].

Amplicons were visualized on agarose gels and quantified using ImageLab Software v6.0 (Bio-Rad Laboratories, Inc., USA). Following quantification, uniform amounts of each amplicon were merged using a Biomek4000 liquid handling robot (Beckman Coulter, USA). The DNA library was then cleaned using 1.0X AMpure beads (Beckman Coulter, USA), size selected using BluePippin (Sage Science, USA), and quantified on a Fragment Analyzer (Advanced Analytical Technologies, Inc., USA) using the High Sensitivity Genomic DNA Kit (Agilent). The library was finally sequenced on MiSeq platform v3 (Illumina, San Diego, CA, USA), alongside samples from other projects.

Bioinformatics data analysis

Bioinformatic processes for the metabarcoding analysis were conducted as in the annotated scripts provided on the following GitHub page: In brief, forward and reverse raw sequencing files obtained following MiSeq sequencing were merged using PEAR 0.9.3 [76] and demultiplexed using the ngsfilter command from the OBITools software suite [77]. The obigrep command was used to select fragments > 400 bp, to check, and further quality filtering was conducted to remove sequences < 100 bp using the fastq_filter command from the USEARCH algorithm [78]. Sequences were then dereplicated using the fastx_uniques command from the USEARCH algorithm [78], and sequences with less than 10 occurrences in the dataset were removed. Our dataset was then denoised using UNOISE algorithm (i.e. unoise3 command from USEARCH) [79], and ZOTUs (e.g., Zero-radius Operational Taxonomic Units) were retrieved. Finally, the taxonomic assignment was performed using the blastn command from the BLAST + application [80]. We applied strict filtering control to remove any false positive detection. For each sample, we first selected all ZOTUs corresponding to a unique species and discarded the ZOTUs that didn’t have at least two reads in at least two PCR replicates. Then, for each retained ZOTUs, we subtracted the highest number of reads which could be found in the corresponding ZOTU in any of all negative controls (extraction blanks and PCR controls). This conservative approach was applied in all PCR replicates of the sample. This was manually done for each analyzed sample of our study. We chose this approach to ensure that potential contamination or “tag-jump” will not lead to potential false positive results. Finally, ZOTUs were manually checked and all ZOTUS corresponding to a unique species were pooled together in a unique species identifier, and the read numbers were added (see Additional file 4.A. and B.) to avoid overinflation of the species diversity detected in this study.



Phytochemical fingerprinting by untargeted UHPLC-QTOF-ESI+MS metabolomics

Multivariate analysis was first applied to discriminate between the five MT genuine samples (code ACM) versus herbal formulations commercialized as teas (T), tablets (Tb), or capsules (C) containing MT as a unique ingredient (U), or in combination with other plant-based ingredients (M), in agreement with the information stated on their labels, as presented in Table 1 and Additional file 1.A. A total number of 217 molecules were separated and included in the matrix for multivariate analysis using Metaboanalyst v5.0 software. Figure 2.A., 3B., and 3 C. show the sPLSDA score plots for the different categories of herbal formulations. First, a comparison between the fingerprints specific to genuine MT seeds (ACM), teas (T), capsules (C), and tablets (Tb) claiming to contain MT ingredients are presented in Fig. 2.A. A significant discrimination was observed between the genuine MT seeds and formulations, with a co-variance of 24.8% for the first 2 components. The Tb group was significantly different in this case, also showing a higher heterogeneity among the six formulations. The C and T groups were partly superposed since capsules and teas contain powders of raw plant tissues. These two groups were discriminated against the genuine ACM group (MT seeds) since they contained other added ingredients and excipients.

Fig. 2
figure 2

The sparse PLSDA score plots showing the discrimination between the different categories of herbal preparations. (A) Comparison of the fingerprints specific to genuine MT seeds (ACM), teas (T), capsules (C), and tablets (Tb) claiming to contain MT ingredients, as unique or multiple combinations. (B) Discrimination between the fingerprints of capsules (C), teas (T), and tablets (Tb). (C) Discrimination between the fingerprints of the herbal preparation containing unique (U) and multiple (M) ingredients

Further, it was seen the discrimination between the fingerprints of capsules, teas, and tablets, therefore similar score plots were obtained, with a co-variance of 29.7% (Fig. 2.B.). Here, the difference between the composition of tablets vs. teas or capsules was more visible. According to Additional file 1.A., one can see that tablets, more than capsules, include non-herbal ingredients (standardized extracts, dextrins, flavonoid pigments) which may explain this discrimination. The influence of unique (U) versus multiple (M) ingredients upon the discrimination between capsules, teas, and tablets was also plotted, the co-variance being 22.7% for the first two components (Fig. 2.C.). Here, also a clear discrimination between TU and TM subgroups of teas, containing unique and multiple ingredients was noticed. This shows that the addition of other plant phytochemicals in tea can be easily identified by untargeted metabolomics. Additionally, the capsules and teas with multiple components have similar fingerprints, different from tablets.

Semi-targeted metabolite profiles

The semi-targeted analysis focused on specific groups of molecules that were previously identified in MT seeds and formulations, namely flavonolignans including silymarin complex, taxifolin, lignan precursors (coumaric acid and coniferyl derivatives), phytosterols, flavonoids, phenolic acids, fatty acids and polar lipid derivatives, a number of 63 molecules being selected and identified (see Additional file 3) using the match of m/z values with HMDB and other databases ( as mentioned in Materials and Methods).

The multivariate analysis focused on the most common molecules that may discriminate and characterize the profile of individual teas, capsules, or tablets. Figure 3.A. presents the heatmap of sample clusters (T and ACM-green, C-red, Tb-Blue) vs. the main 25 molecules responsible for the discrimination, as selected by Metaboanalyst algorithm. Specific MT molecules like silybins A + B and silyhermin are readily identified, as well as phenolic acids, flavonoids, and polar lipids as putative biomarkers for discrimination. Some capsules (C1/M2, C2-C4/M3-M5) and tablets (Tb1/U1 and Tb2/U2) showed specifically higher levels of such molecules. The Random Forest (RF) analysis (Fig. 3.B.) classified the top 15 molecules to be considered most significant as putative biomarkers, according to Mean Decrease Accuracy (MDA) values > 0.002.

Fig. 3
figure 3

(A) The Heatmap showing the clusters of the three groups of samples C, T, and Tb vs. the top of 25 molecules selected as most relevant for the discrimination between the teas (T/U vs. T/M), four multiple ingredient capsules (C1-C4/M2-M5), and tablets T/U, T/M. (B) The RF analysis plot showing the top of 15 molecules to be considered as potential biomarkers, according to the Mean Decrease Accuracy value

Considering the mean values per formulation, one can see that, for example in the Tb group, Silybins A + B had significantly higher levels, followed by capsules (C) and teas (T). Considering all 15 molecules, the highest levels of phytochemicals were found in capsules (C) with 9 out of 15 molecules. Focusing on these 15 to 25 molecules, out of the 63 separated and identified by UHPLC-QTOF-MS, would be an effective approach for developing qualitative or quantitative evaluation of these herbal supplements.

Profiles of silymarin complex

Since silymarin flavonolignans are specifically related to genuine MT products, a more targeted analysis focused on this subclass of molecules and compared their levels in the individual products with unique (U) or multiple ingredients (M). Figure 4.A. and 4.B. show the MS peak intensities of the silymarin class of molecules, which may be considered as relevant biomarkers of product quality and authenticity. The targeted silymarin flavonolignans were silybins A + B, silychristin and silydianin, taxifolin, dehydrosilybin, silyhermin. The total intensity, as the sum of all silymarin flavolignans was also calculated, as presented in green columns.

Fig. 4
figure 4

(A) Comparative values of the mean MS peak intensities recorded for the genuine MT seeds (ACM) comparative to herbal preparations (teas – T; tablets- Tb), which declared to have a unique ingredient (U). (B) Comparative values of MS peak intensities for the herbal preparations (teas – T; tablets- Tb; capsules-C; emulsion-E) with multiple ingredients (M). The error bars (± SD) from triplicate measurements represented 20–30% of the mean values represented in the graphic

Silybins were identified in all samples followed by silychristin and silydianin. In the group U (unique ingredient) of samples (Fig. 4.A.), the richest preparations were Tb2/U2, Tb1/U1, Tb3/U3, Tb6/U7 followed by T3/U6, and T2/U4. The sample Tb4/U5 had the smallest content, but still with an acceptable content of silybins, almost 50% of the genuine samples (the mean value from the ACM group was considered). Taxifolin, silychristin and silydianin were found in higher levels especially in samples Tb2/U2, Tb1/U1, and Tb3/U3.

In the group M (multiple ingredients) of samples (Fig. 4.B.), the richest preparations were C4/M5, C3/M4, and C1/M2 followed by T1/M1 and T4/M6, Tb5/M10, E/M9, T7/M11. The samples C2/M3 and T6/M8 had the smallest content of silybins. Silychristin and silydianin were found in C3/M4, C4/M5, C1/M2, and T1/M1, at levels representing around 60% of the levels found in the samples with unique ingredients. Taxifolin was also found, but at lower levels.

Generally, tea products with multiple ingredients contained lower levels of silymarin complex. Meanwhile, the capsules and especially tablets contained higher levels of silymarin complex probably due to the use of standardized, concentrated MT extracts as ingredients, as can be seen in samples Tb1-Tb3, C3, and C4. Taxifolin was identified at higher levels in tablets Tb3-Tb6, but also in some teas and capsules.

As expected, the ratio between the silymarin complex in the U vs. M group was around 2, as presented in Additional file 5. A significant variability was noticed, explained by the different ingredients used and claimed on the label and the type of formulation (teas vs. tablets vs. capsules).

To authenticate by accurate analysis is still difficult since, as can be seen in Additional file 1.A., the producers of these herbal formulations did not report the concentration of active compounds belonging to the silymarin complex on the product labels. Although semi-targeted analysis could confirm the presence of silymarins as well as their relative levels, an accurate comparison with the claimed composition mentioned on the label was not possible.

Quantitative evaluation of silybins using UV-spectrometry and UHPLC-QTOF-ESI+MS

The semi-targeted analyses showed that the silymarin complex compounds were especially suitable as quality indicators. In order to further evaluate these compounds for authenticity assessment, calibrations curves were built with analytical pure standards of silybins A + B, using UHPLC-QTOF-ESI+MS for accurate calculation of silybins and UV spectrometry (as a fast, less accurate but indicative method of silymarin complex-absorbing molecules in ethanol at 288 nm) as presented in Additional file 2.A.

Based on the calibration curves, a comparative evaluation of silybins concentrations (mg/g d.m.) in all herbal preparations, as determined both, by UV spectrometry and UHPLC-QTOF-ESI+MS analysis is presented in Fig. 5.

Fig. 5
figure 5

Comparative evaluation of silybins concentrations (mg silybins/g sample) in all herbal preparations, as determined by UV spectrometry and UHPLC-QTOF-ESI+MS analysis. For codes and abbreviations see Table 1 and Additional file 1.A. The error bars (± SD) from triplicate measurements represented 20–30% of the mean values represented in graphic

With only a few exceptions (samples T3/U6 and C4/M5), the silybins concentrations determined by UHPLC-QTOF-ESI+MS were around 2–3 times lower compared to data released by UV spectrometry in unique ingredient samples (group U) and up to 10 times lower in multiple ingredient samples (group M), respectively. This is explained by the UV absorbance of other phenolic derivatives, besides silybins, at 288 nm. Therefore, the UV analysis overestimates the silybins concentration and is indicating mostly the pool of flavonoids, including flavonolignans. Meanwhile the UHPLC-QTOF-ESI+MS, gives a much more accurate information, regarding the types and levels of the different molecules, offering a more complete picture of the identity of the products. Nevertheless, these results show that UV-spectrometry can be applied as a preliminary, rough evaluation of the formulation, whereas UHPLC-QTOF-ESI+MS may target more precisely the molecules to be authenticated and also other components, in a semi-targeted or quantitative way.

However, in this study, the information provided on the label by the producers was not precisely mentioned, and in many cases lacking important details. For instance, it was not always clearly defined what the proportions were between an ingredient such as MT standardized extracts and additional powdered MT seeds, or the concentration of key-compounds responsible for the claimed effect, e.g., other ingredients added in the preparation.

DNA metabarcoding

Qubit fluorometer quantitation showed large differences in total DNA concentrations among the eighteen analyzed MT-based botanical preparations (see Additional file 6). The sequencing success rate was 100% (18/18 samples). A dataset consisting of 1,163,356 reads fulfilling our initial trimming and filtering quality criteria was obtained, with an average of 77,557 reads per sample. Zero-radius operational taxonomic units (ZOTUs) were obtained for all preparations (100%) (see Additional file 4.A.). Sixteen preparations (89%) had ZOTUs that passed bioinformatics trimming and filtering quality criteria that require ZOTUs to have at least 10 reads in the whole dataset and at the samples level, more than one read and being detected in at least 2 out of 3 replicates in order to be retained for further analysis. Two samples (tablets PA2 and PA5) did not fulfill the imposed criteria and were excluded (i.e., there was no read remaining for any ZOTU following the various filtering steps). ZOTUs and their read numbers for the same species were merged for further analysis. Across all sixteen retained samples, a total of 59 different species (declared and non-declared on the label), were identified using the basic local alignment search tool (BLAST) from the retained ZOTUs.

The main targeted plant ingredient - S. marianum (milk thistle), was detected in eleven out of sixteen retained preparations. Out of the five single ingredient samples - those containing only S. marianum according to the label, S. marianum was detected in four samples (PA3, PA7, PA9, PA10) and in one not (PA17). Out of eleven multiple ingredient samples - those containing S. marianum together with other species according to the label, S. marianum was detected in seven samples (PA1, PA6, PA8, PA11, PA12, PA14, PA15) and in four not (PA4, PA13, PA16, PA18). The fidelity for S. marianum in single-ingredient products was 80% (4 out of 5), and for multi-ingredient products, 64% (7 out of 11) (see Additional file 7.A and 7.B).

All five retained single-ingredient samples contained species not mentioned on the label. Two of the multi-ingredient samples contained all species listed on the label (the capsules PA6 and PA8), but both also contained off-label species, and seven contained fewer species than listed on the label (PA1, PA4, PA11, PA12, PA14, PA15) and apart from PA14 they all contained additional off label species. Two samples contained none of the species from the label (PA13, PA16, PA18) but instead contained off-label species. The overall ingredient fidelity (detected species from product label/total number of species on the label) for multi-ingredient products was 45% and for all products 56% (see Additional file 8).

A total of 47 species non-listed on the label were detected. The most abundant species by sequence reads were Hordeum vulgare L., Urtica radicans Wight, Viola sp. (Viola arcuata Blume; Viola arvensis Murray), Avena sativa L., and Helianthus sp. (Helianthus divaricatus L.; H. giganteus L.; Helianthus grosseserratus M.Martens). The plant taxa detected in the samples are presented in Fig. 6.

Fig. 6
figure 6

Sankey diagram summarizing detected species (declared and non-declared on the label), from the retained samples, using DNA metabarcoding. Only species represented by ZOTUs detected in ≥2 replicates and with ≥2 reads are shown. Sizes of flows denote proportions of reads at the species level

Comparative results

The results from the UHPLC-QTOF-ESI+MS untargeted analysis gave initial valuable information about the general fingerprint of the different categories of formulations, identifying a large number of molecules (217) that can be used as specific indicators of genuine ingredients (MT seeds and other plant components) and subgroups of herbal formulations, such as tablets that may have different fingerprints due to the more complex non-herbal ingredients. Nevertheless, the untargeted analysis does not offer enough information to make a more precise authentication of individual products and to find the relevant biomarkers for their identity.

The semi-targeted analysis, which focused on 63 molecules belonging to relevant phytochemicals for MT identification, gave better indications for the key molecules to be considered as authenticity biomarkers. The silymarin complex compounds (silybins and taxifolin) showed to be relevant as biomarkers of authenticity. The quantitative evaluation applied comparatively, using UV-spectrometry (based on absorbance at 288 nm and expressed in silybin-equivalents) and the accurate determination of silybins by UHPLC-QTOF-ESI+MS showed different contributions of MT-based ingredients in unique-ingredient products and multiple-ingredient products.

DNA metabarcoding had a good resolution in detecting MT at the species level and provided insights into the total species composition of herbal preparations labeled as containing unique (U) and multi-ingredients (M). The DNA metabarcoding and LC-MS semi-targeted metabolomics results were in accordance for eleven samples (61%). LC-MS could validate the presence of S. marianum in another five samples (100%) that have not passed the filtering criteria of DNA metabarcoding analysis.


It is well-accepted that botanicals and their derived herbal supplements are susceptible to various issues that raise serious quality and safety concerns [33]. Challenges may occur throughout the value chain, from cultivation or wild harvesting of the medicinal plants as sources of raw material to the final marketed product [35, 81].

In spite of the less exigent legislative regulations regarding herbal supplements, many laboratories apply new analytical approaches for the analysis of milk thistle content in raw materials and/or derived herbal formulations, including high-performance liquid chromatography (HPLC) with different types of detectors, thin layer chromatography (TLC), high-performance thin layer chromatography (HPTLC), hyphenated mass spectrometry, but also UV spectrometry [17, 52, 55, 82, 83]. However, milk thistle preparations are often highly processed and usually mixed with other plant ingredients, limiting the accuracy of traditional analytical methods in identifying the targeted plant species, and making it even more challenging to detect non-target species. Hence, applying new fit-for-purpose technologies and methodologies will perhaps enable a more accurate quality assessment of milk thistle-derived preparations [36, 65]. In this study, we combined two emerging technologies - metabolomics and DNA metabarcoding for the authentication of milk thistle-based preparations.

Generally, metabolomics is defined as the holistic qualitative and quantitative measurement of the complete set of small metabolites in a biological system at a given time [84,85,86]. As one of the most rapidly evolving fields, metabolomics found its applications in a plethora of basic and applied studies of the life sciences [87, 88]. Recent innovation and progress in metabolomics technologies have been used also to address a wide range of biological questions within the field of natural products [89,90,91,92]. These developments have opened new perspectives in the field of botanicals and derived herbal preparations by providing, among others, powerful tools for their authentication and quality assessment [37, 93]. Metabolomics comprises methods and advanced analytical platforms, with a high degree of sensitivity, selectivity, and reproducibility that enable a broader insight into the highly diverse metabolome complexity [94, 95], as already shown in several studies focusing on metabolic profiling in complex mixtures [96, 97]. Nevertheless, each analytical method has its own advantages and disadvantages, and the choice of a certain method is typically driven by the focus of the study, followed inter alia by the nature of the samples, costs, or accessibility [96,97,98,99]. UHPLC-QTOF-MS using ESI+ fragmentation is a versatile technique that imparts great promise for the comprehensive authentication of botanicals and botanical preparations and was the method of choice in this study. Here, untargeted and semi-targeted metabolomics as well as a quantitative evaluation of silybins as biomarkers of MT presence in different herbal formulations have been used and compared for their analytical efficacy in the context of authenticity and quality control of milk thistle-derived commercial preparations.

While targeted metabolomics has limited coverage of the metabolome as it aims to measure a predefined set of known metabolites, untargeted metabolomics focuses on a rather wider coverage, or ideally, complete measurement, of the relative levels of metabolites in a sample [100], enabling the simultaneous comparison of several samples without having a priori information about their content or suspicion of contamination or adulteration [84]. Here, by untargeted approach, there were identified molecules that can compare unique (including genuine samples) ingredient samples against multi-ingredient MT preparations and various subgroups of herbal formulations, in an attempt to find possible metabolites to be used as key markers in the authentication process. However, with the untargeted approach, it was not possible to find a single precise authentication of the products. Meanwhile, semi-targeted metabolomics is a promising, alternative, enabling the measurement of predefined metabolites known to reflect the presence of MT [101] e.g. flavonolignans. In our study, the semi-targeted analysis focused on 63 relevant phytochemicals and gave good indications regarding the key molecules to be further considered as authenticity markers, the silymarin complex (silybins and taxifolin) being the most relevant. MT was found to be present in all preparations, in agreement with the genuine MT seed composition, but at different levels. The comparative evaluation using both UV spectrometry and UHPLC-QTOF-ESI+MS showed different contributions of MT-based ingredients. These findings corroborate previous results showing various degrees of substitution and adulteration of the botanicals [102,103,104,105,106,107,108].

However, the usefulness of metabolomics based on specific phytochemicals is still a challenge for herbal product authentication and has a limited capacity to detect other botanical ingredients e.g. contaminants or adulterants [37].

While metabarcoding is focused on ingredient authentication by DNA recognition (qualitative), the metabolomic approach is looking at molecules that can be found either in the key ingredient (to check its presence and quantity in the product) or in other ingredients (which are usually found in plant mixtures of teas, powders of standardized extracts in capsules or tablets, etc.) or even non-declared excipients. Here we tried to see the capability of the metabolomic approach to identify similarities and differences between these products compared to genuine plant metabolites. Also, we followed in parallel the classical “phytochemical analysis” identifying and quantifying just the molecules belonging to the silymarin complex, e.g., silybins A + B and taxifolin.

DNA metabarcoding brings together the innovation of high-throughput sequencing (HTS) technologies and the DNA barcoding concept, enabling simultaneous multi-taxa identification from a pool of genetic material containing DNA from different origins [109, 110]. This approach generated an emerging area of research with practical applicability in the analysis of species composition of a wide range of multi-ingredient and highly processed samples, being used today in regulatory, conservation, and commercial contexts [38, 111,112,113,114,115]. This has lately emerged as a cost-effective and reliable method to improve the authentication and quality control process of botanicals and derived preparations [38, 65, 116].

Here, DNA metabarcoding results indicate a high level of inconsistencies between the identified species and those listed on the labels of the sixteen retained preparations. The targeted plant ingredient, milk thistle, was detected in 11 (68,75%) out of sixteen retained preparations, including three unique and eight multi-ingredient preparations (five herbal teas, two tablets, three capsules, and one emulsion). However, we emphasize that in four other preparations (PA2, PA4, PA5, and PA13) MT was detected – but we couldn’t validate its positive identification since the samples did not fulfill the imposed bioinformatics quality requirements. Either way, semi-targeted metabolomics confirmed the presence of MT in all four products. Six preparations contained all the listed plant-based ingredients (P3, P6, P7, P8, P9, P10), but additional plant species were detected in all of them. According to the label, this included four unique and two multiple ingredients (two tablets, two capsules and two herbal teas) preparations. Two samples (PA2 and PA5) did not fulfill the trimming and filtering quality criteria and they were not considered in the results and discussion. The findings corroborate previous studies, showing significant incongruences between the detected species and those listed on the labels of some marketed botanical preparations [72, 117,118,119,120,121,122].

Considering the nature of sourced botanicals and the long value chain to the final preparation, the discrepancies between the species detected using DNA metabarcoding and those listed on the product labels require a careful evaluation regarding the possible sources of contamination or adulteration. In this study, we used the information found on the label/leaflet for each preparation to define some variables that can impact the interpretation of the results. Thus, the evaluation of the authentication results was made in line with a priori information such as the origin and cultivation conditions, and a posteriori information such as the taxonomic identification of ZOTUs. In this regard, we highlight also that DNA metabarcoding is a very sensitive method, and even traces of another species or a pollen grain will give a positive identification. For instance, in this study, various wind-pollinated plant species (anemophilous) were detected and their presence in our results can be expected and considered as normal trace contamination – if they are in quantities that do not pose any quality issues of the preparation and if they are within the permitted contamination range. However, quantifying the relative species abundance within a sample is beyond the technical characteristics of DNA metabarcoding since many potentially confounding factors can affect read numbers, thereby, in these cases, appropriate phytochemistry methods should be used on a punctual-based evaluation when contamination or adulteration is suspected. The off-label detected species can be interpreted as contamination or adulteration, but also as being generated by amplification bias (i.e. PCR chimeras), sequencing errors, or false-positive taxonomic identifications due to errors in the barcode sequences reference databases [123,124,125]. To mitigate false positives and to increase the overall reliability of results, in this study we used technical replicates based on three independent PCR amplified products from the same preparation. Further, very strict filtering and trimming thresholds for sequence reads were applied to overcome sequencing errors, followed by very conservative selection standards in order to validate a positive identification. The varying degrees of success in identifying the species listed on the labels can be due to a failure to detect them with metabarcoding. This can be explained by false negatives that are often due to a combination of highly degraded DNA resulting from harvesting, drying, storage, transportation, and processing [119, 126, 127], the inability of recovering DNA due to the presence of pharmaceutical excipients affecting DNA extraction [128], or as result of poor primer fit and amplification biases [129], stochasticity due to low DNA concentrations [130], or incomplete reference databases.

Proper analytical validation of DNA metabarcoding is necessary before this can be implemented for molecular diagnostics, both in quality monitoring programs in a regulatory context, and in supply chain management systems by the industry sector. Important steps have been taken toward validating and standardizing DNA metabarcoding for quality control in commercial applications and regulatory contexts. A very good practical example is the study commissioned by the Federal Office of Consumer Protection and Food Safety (BVL) in Germany [131]. In this study within an inter-laboratory ring trial including 15 laboratories, the reproducibility, robustness, and measurement of DNA metabarcoding uncertainties, have been analyzed using meat-based multi-ingredient samples. The study concluded that DNA metabarcoding is a robust authentication tool and can be used in routine analysis by official food control laboratories [131]. While some DNA barcoding methods are validated and standardized for quality control in commercial applications and regulatory contexts [132], so far, no similar large inter-laboratory DNA metabarcoding protocols were performed for its validation as an authentication tool in the field of botanicals and their derived preparations. Even if DNA metabarcoding addresses a number of limitations when using plant-based samples, we expect that a common effort for a validation study will be performed and propose a DNA metabarcoding protocol applicable to the quality control systems of botanical preparations.

The application of emerging and innovative techniques and fit-for-purpose methodologies to advance the evaluation of botanical preparations in the context of quality assessment is strongly advocated today [36, 133]. Each analytical technique has its benefits and limitations, and interdisciplinary approaches have been shown to improve the quality assessment process of botanical preparations [65, 116, 133, 134]. The results of this study corroborate previous results confirming the advantages of combining analytical approaches for the quality assessment of botanical preparations [72, 121, 135].


This study used untargeted and semi-targeted metabolomics analysis based on UHPLC-QTOF-ESI+MS data and UV spectrometry, alongside high-throughput DNA metabarcoding using Illumina MiSeq to authenticate eighteen botanical preparations labeled as containing Silybum marianum (L.) Gaertn. (milk thistle) either as a unique ingredient or in combination with other plant-based ingredients. The results confirm that DNA metabarcoding using Illumina MiSeq can be used to test for the presence of S. marianum and simultaneously to detect other plant ingredients within complex herbal preparations with results to be interpreted in a broad context. It should be emphasized however that DNA metabarcoding detected milk thistle in only eleven out of sixteen retained preparations, and the other two had incomplete evidence of milk thistle despite metabolomics validating its presence, challenging its use as a stand-alone approach for routine screening. Moreover, the high sensitivity of DNA metabarcoding requires careful consideration of the total species composition detected by interpreting the results in a broad context, particularly concerning the detection of false positives versus possible contaminants and adulterants. Further, DNA metabarcoding does not provide information on the active metabolites of the botanical preparations, and this narrows its analytical capabilities to the identification of target species and confirmation of presence, but not the absence of other species. The clear advantage of semi-targeted metabolomics based on UHPLC-QTOF-ESI+MS consisted in the analytical ability to detect the quantity of the predefined set of phytochemical markers compounds and showing clearly that all investigated milk thistle preparations contained molecules from silymarin complexes at different concentrations. Moreover, metabolomics realized a wider coverage of the relative levels of other metabolites enabling the comparison and discrimination between the different groups of formulations without having a priori information about their content. This study shows that the combination of complementary methods offers a robust analytical approach to advance authentication and quality control of botanical preparations.

Data Availability

The DNA sequencing read data generated and analyzed during the current study are available in Zenodo



Milk thistle


ultra-high-performance liquid chromatography coupled with electrospray -quadrupole-time of flight-mass spectrometry using the positive ionization


Ulraviolt visible


Dry matter


Quality control


Mass-to-charge ratio


Principal Component Analysis


Sparse Partial Least Square Discriminant Analysis


Random Forest


Mean Decrease Accuracy


Polymerase chain reaction


Nuclear ribosomal internal transcribed spacer


Zero-radius operational taxonomic unit


  1. Schuppan D, Jia J-D, Brinkhaus B, Hahn EG. Herbal products for liver diseases: a therapeutic challenge for the new millennium. Hepatology. 1999;30:1099–104.

    Article  CAS  PubMed  Google Scholar 

  2. Post-White J, Ladas EJ, Kelly KM. Advances in the use of milk thistle (Silybum marianum). Integr Cancer Ther. 2007;6:104–9.

    Article  CAS  PubMed  Google Scholar 

  3. Ross SM. Milk thistle (Silybum marianum): an ancient botanical medicine for modern times. Holist Nurs Pract. 2008;22:299–300.

    Article  PubMed  Google Scholar 

  4. Siegel AB, Stebbing J. Milk thistle: early seeds of potential. Lancet Oncol. 2013;14:929–30.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Krepkova LV, Babenko AN, Saybel’ OL, Lupanova IA, Kuzina OS, Job KM et al. Valuable hepatoprotective plants - how can we pptimize waste free uses of such highly versatile resources? Front Pharmacol 2021;12.

  6. Raclariu-Manolică AC, Socaciu C. Detecting and profiling of milk thistle metabolites in food supplements: a safety-oriented approach by advanced analytics. Metabolites. 2023;13:440.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Blumenthal M, Bussmann RW, Goldberg A, Gruenwald J, Hall T, Riggins CW, et al. The complete german commission E monographs: therapeutic guide to herbal medicines. Integrative medicine communications. Austin, TX/Boston, MA,: American Botanical Council; 1998.

    Google Scholar 

  8. Gillessen A, Schmidt HH-J. Silymarin as supportive treatment in liver diseases: a narrative review. Adv Ther. 2020;37:1279–301.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Zecca E, Zuppa AA, D’Antuono A, Tiberi E, Giordano L, Pianini T, et al. Efficacy of a galactogogue containing silymarin-phosphatidylserine and galega in mothers of preterm infants: a randomized controlled trial. Eur J Clin Nutr. 2016;70:1151–4.

    Article  CAS  PubMed  Google Scholar 

  10. Bazzano AN, Cenac L, Brandt AJ, Barnett J, Thibeau S, Theall KP. Maternal experiences with and sources of information on galactagogues to support lactation: a cross-sectional study. Int J Womens Health. 2017;9:105–13.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Lawrence RM, Lawrence RA. 11 - medications, herbal preparations, and natural products in breast milk. In: Lawrence RA, Lawrence RM, editors. Breastfeeding (ninth edition). Philadelphia: Elsevier; 2022. pp. 326–92.

    Chapter  Google Scholar 

  12. Greenlee H, Abascal K, Yarnell E, Ladas E. Clinical applications of Silybum marianum in Oncology. Integr Cancer Ther. 2007;6:158–65.

    Article  CAS  PubMed  Google Scholar 

  13. Tajmohammadi A, Razavi BM, Hosseinzadeh H. Silybum marianum (milk thistle) and its main constituent, silymarin, as a potential therapeutic plant in metabolic syndrome: a review. Phytother Res. 2018;32:1933–49.

    Article  CAS  PubMed  Google Scholar 

  14. Mengs U, Pohl R-T, Mitchell T, Legalon® SIL. The antidote of choice in patients with acute hepatotoxicity from amatoxin poisoning. Curr Pharm Biotechnol. 2012;13:1964–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Lee DY-W, Liu Y. Molecular structure and stereochemistry of silybin A, silybin B, isosilybin A, and isosilybin B, isolated from Silybum marianum (milk thistle). J Nat Prod. 2003;66:1171–4.

    Article  CAS  PubMed  Google Scholar 

  16. Smith WA, Lauren DR, Burgess EJ, Perry NB, Martin RJ. A silychristin isomer and variation of flavonolignan levels in milk thistle (Silybum marianum) fruits. Planta Med. 2005;71:877–80.

    Article  CAS  PubMed  Google Scholar 

  17. Csupor D, Csorba A, Hohmann J. Recent advances in the analysis of flavonolignans of Silybum marianum. J Pharm Biomed Anal. 2016;130:301–17.

    Article  CAS  PubMed  Google Scholar 

  18. Bijak M. Silybin, a major bioactive component of milk thistle (Silybum marianum L. Gaernt.)—chemistry, bioavailability, and metabolism. Molecules 2017;22:1942.

  19. Aziz M, Saeed F, Ahmad N, Ahmad A, Afzaal M, Hussain S, et al. Biochemical profile of milk thistle (Silybum Marianum L.) with special reference to silymarin content. Food Sci Nutr. 2021;9:244–50.

    Article  CAS  PubMed  Google Scholar 

  20. AbouZid SF, Chen S-N, Pauli GF. Silymarin content in Silybum marianum populations growing in Egypt. Ind Crops Prod. 2016;83:729–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Javeed A, Ahmed M, Sajid AR, Sikandar A, Aslam M, Hassan T, ul, et al. Comparative assessment of phytoconstituents, antioxidant activity and chemical analysis of different parts of milk thistle Silybum marianum L. Molecules. 2022;27:2641.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Abenavoli L, Izzo AA, Milić N, Cicala C, Santini A, Capasso R. Milk thistle (Silybum marianum): a concise overview on its chemistry, pharmacological, and nutraceutical uses in liver diseases. Phytother Res. 2018;32:2202–13.

    Article  PubMed  Google Scholar 

  23. Rainone F. Milk thistle. Am Fam Physician. 2005;72:1285–8.

    PubMed  Google Scholar 

  24. Calani L, Brighenti F, Bruni R, Del Rio D. Absorption and metabolism of milk thistle flavanolignans in humans. Phytomedicine. 2012;20:40–6.

    Article  CAS  PubMed  Google Scholar 

  25. Zhu H-J, Brinda BJ, Chavin KD, Bernstein HJ, Patrick KS, Markowitz JS. An assessment of pharmacokinetics and antioxidant activity of free silymarin flavonolignans in healthy volunteers: a dose escalation study. Drug Metab Dispos. 2013;41:1679–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Khazaei R, Seidavi A, Bouyeh M. A review on the mechanisms of the effect of silymarin in milk thistle (Silybum marianum) on some laboratory animals. Veterinary Med Sci. 2022;8:289–301.

    Article  CAS  Google Scholar 

  27. Garcia-Alvarez A, Egan B, de Klein S, Dima L, Maggi FM, Isoniemi M, et al. Usage of plant food supplements across six european countries: findings from the PlantLIBRA consumer survey. PLoS ONE. 2014;9:e92265.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Smith T, Majid F, Eckl V, Reynolds CM. Herbal Supplement sales in US increase by record-breaking 17.3% in 2020 2017:14.

  29. Ekor M. The growing use of herbal medicines: issues relating to adverse reactions and challenges in monitoring safety. Front Pharmacol. 2014;4:177.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Binns CW, Lee MK, Lee AH. Problems and prospects: public health regulation of dietary supplements. Annu Rev Public Health. 2018;39:403–20.

    Article  PubMed  Google Scholar 

  31. Fitzgerald M, Heinrich M, Booker A. Medicinal plant analysis: a historical and regional discussion of emergent complex techniques. Front Pharmacol. 2020;10.

  32. World Health Organization (WHO). (2007). World Health Organization guidelines for assessing quality of herbal medicines with reference to contaminants and residues.

  33. World Health Organization (WHO). World Health Organization traditional medicine strategy: 2014–2023. WHO; 2013.

  34. Rietjens IMCM, Slob W, Galli C, Silano V. Risk assessment of botanicals and botanical preparations intended for use in food and food supplements: emerging issues. Toxicol Lett. 2008;180:131–6.

    Article  CAS  PubMed  Google Scholar 

  35. Heinrich M. Quality and safety of herbal medical products: regulation and the need for quality assurance along the value chains. Br J Clin Pharmacol. 2015;80:62–6.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Thakkar S, Anklam E, Xu A, Ulberth F, Li J, Li B, et al. Regulatory landscape of dietary supplements and herbal medicines from a global perspective. Regul Toxicol Pharmacol. 2020;114:104647.

    Article  PubMed  Google Scholar 

  37. Simmler C, Graham JG, Chen S-N, Pauli GF. Integrated analytical assets aid botanical authenticity and adulteration management. Fitoterapia. 2018;129:401–14.

    Article  PubMed  Google Scholar 

  38. De Boer HJ, Ichim MC, Newmaster SG. DNA barcoding and pharmacovigilance of herbal medicines. Drug Saf. 2015;38:611–20.

    Article  CAS  PubMed  Google Scholar 

  39. Dwyer JT, Coates PM, Smith MJ. Dietary supplements: regulatory challenges and research resources. Nutrients. 2018;10:41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Heinrich M, Appendino G, Efferth T, Fürst R, Izzo AA, Kayser O, et al. Best practice in research – overcoming common challenges in phytopharmacological research. J Ethnopharmacol. 2020;246:112230.

    Article  CAS  PubMed  Google Scholar 

  41. Kunle. Standardization of herbal medicines - a review. Int J Biodvers Conserv. 2012;4.

  42. Zhang J, Wider B, Shang H, Li X, Ernst E. Quality of herbal medicines: challenges and solutions. Complement Ther Med. 2012;20:100–6.

    Article  PubMed  Google Scholar 

  43. Bilia AR. Science meets regulation. J Ethnopharmacol. 2014;158(Pt B):487–94.

    Article  Google Scholar 

  44. Ichim MC. The DNA-Based authentication of commercial herbal products reveals their globally widespread adulteration. Front Pharmacol 2019;10.

  45. Saller R, Meier R, Brignoli R. The use of silymarin in the treatment of liver diseases. Drugs. 2001;61:2035–63.

    Article  CAS  PubMed  Google Scholar 

  46. Xue Y, Sheng Y, Wang J, Huang Q, Zhang F, Wen Y et al. Fast screening and identification of illegal adulterated glucocorticoids in dietary supplements and herbal products using UHPLC-QTOF-MS with All-Ion Fragmentation Acquisition combined with characteristic fragment Ion List classification. Frontiers in Chemistry 2021;9.

  47. Šimánek V, Kren V, Ulrichová J, Vicar J, Cvak L. Silymarin: What is in the name ...? An appeal for a change of editorial policy. Hepatology 2000;32:442–4.

  48. Kroll DJ, Shaw HS, Oberlies NH. Milk thistle nomenclature: why it matters in cancer research and pharmacokinetic studies. Integr Cancer Ther. 2007;6:110–9.

    Article  CAS  PubMed  Google Scholar 

  49. Martin RJ, Lauren DR, Smith WA, Jensen DJ, Deo B, Douglas JA. Factors influencing silymarin content and composition in variegated thistle (Silybum marianum). Null. 2006;34:239–45.

    Article  CAS  Google Scholar 

  50. Habán M, Habánová M, Otepka P, Kobida Ľ. Milk thistle (Silybum marianum [L.] Gaertn. Cultivated in polifunctional crop rotation and its evaluation 2010:7.

  51. Anthony K, Saleh MA. Chemical profiling and antioxidant activity of commercial milk thistle food supplements. J Chem Pharm Res 2012;4.

  52. Chambers CS, Holečková V, Petrásková L, Biedermann D, Valentová K, Buchta M, et al. The silymarin composition? and why does it matter??? Food Research International 2017;100:339–53.

  53. Pendry BA, Kemp V, Hughes MJ, Freeman J, Nuhu HK, Sanchez-Medina A, et al. Silymarin content in Silybum marianum extracts as a biomarker for the quality of commercial tinctures. J Herb Med. 2017;10:31–6.

    Article  Google Scholar 

  54. Marceddu R, Dinolfo L, Carrubba A, Sarno M, Di Miceli G. Milk thistle (Silybum marianum L.) as a novel multipurpose crop for agriculture in marginal environments: a review. Agronomy. 2022;12:729.

    Article  CAS  Google Scholar 

  55. Fenclova M, Novakova A, Viktorova J, Jonatova P, Dzuman Z, Ruml T, et al. Poor chemical and microbiological quality of the commercial milk thistle-based dietary supplements may account for their reported unsatisfactory and non-reproducible clinical outcomes. Sci Rep. 2019;9:11118.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Viktorova J, Stranska-Zachariasova M, Fenclova M, Vitek L, Hajslova J, Kren V, et al. Complex evaluation of antioxidant capacity of milk thistle dietary supplements. Antioxidants. 2019;8:317.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Stickel F, Egerer G, Seitz HK. Hepatotoxicity of botanicals. Public Health Nutr. 2000;3:113–24.

    Article  CAS  PubMed  Google Scholar 

  58. Tournas VH, Rivera Calo J, Sapp C. Fungal profiles in various milk thistle botanicals from US retail. Int J Food Microbiol. 2013;164:87–91.

    Article  CAS  PubMed  Google Scholar 

  59. Frommenwiler DA, Sharaf MHM, Reich E. The truth behind herbal products: how HPTLC can help herbal industry detect adulteration? Planta Med. 2019;85:ISL–EA.

    Google Scholar 

  60. McCutcheon A. Adulteration of milk thistle (Silybum marianum). Botanical Adulterants Prevention Bulletin (Austin, TX: ABC-AHPNCNPR Botanical Adulterants Prevention Program) 2020.

    Google Scholar 

  61. Eklund L, Simon JP, Ballenger J. High performance liquid chromatography of flavonolignans in commercial milk thistle supplements. Bios. 2009;80:164–9.

    Article  CAS  Google Scholar 

  62. Booker A, Heinrich M. Value chains of botanical and herbal medicinal products: a european perspective. HerbalGram. 2016;112:40–5.

    Google Scholar 

  63. De Boer HJ, Ichim MC, Newmaster SG. DNA barcoding and pharmacovigilance of herbal medicines. Drug Saf. 2015;38:611–20.

    Article  CAS  PubMed  Google Scholar 

  64. Lüde S, Vecchio S, Sinno-Tellier S, Dopter A, Mustonen H, Vucinic S, et al. Adverse effects of plant food supplements and plants consumed as food: results from the poisons centres-based PlantLIBRA Study. Phytother Res. 2016;30:988–96.

    Article  PubMed  Google Scholar 

  65. Raclariu AC, Heinrich M, Ichim MC, de Boer H. Benefits and limitations of DNA barcoding and metabarcoding in herbal product authentication. Phytochem Anal. 2018;29:123–8.

    Article  CAS  PubMed  Google Scholar 

  66. Pang Z, Zhou G, Ewald J, Chang L, Hacariz O, Basu N, et al. Using MetaboAnalyst 5.0 for LC-HRMS spectra processing, multi-omics integration and covariate adjustment of global metabolomics data. Nat Protoc. 2022;17:1735–61.

    Article  CAS  PubMed  Google Scholar 

  67. Wishart DS, Guo A, Oler E, Wang F, Anjum A, Peters H, et al. HMDB 5.0: the human metabolome database for 2022. Nucleic Acids Res. 2022;50:D622–31.

    Article  CAS  PubMed  Google Scholar 

  68. Fahy E, Subramaniam S, Murphy RC, Nishijima M, Raetz CRH, Shimizu T, et al. Update of the LIPID MAPS comprehensive classification system for lipids. J Lipid Res. 2009;50:9–14.

    Article  CAS  Google Scholar 

  69. Rothwell JA, Perez-Jimenez J, Neveu V, Medina-Remón A, M’Hiri N, García-Lobato P, et al. Phenol-explorer 3.0: a major update of the phenol-explorer database to incorporate data on the effects of food processing on polyphenol content. Database. 2013;2013:bat070.

    Article  PubMed  PubMed Central  Google Scholar 

  70. Kim S, Chen J, Cheng T, Gindulyte A, He J, He S, et al. PubChem in 2021: new data content and improved web interfaces. Nucleic Acids Res. 2021;49:D1388–95.

    Article  CAS  PubMed  Google Scholar 

  71. Doyle JJ, Doyle JL. Isolation of plant DNA from fresh tissue. Focus. 1987;12:13–5.

    Google Scholar 

  72. Raclariu AC, Paltinean R, Vlase L, Labarre A, Manzanilla V, Ichim MC, et al. Comparative authentication of Hypericum perforatum herbal products using DNA metabarcoding, TLC and HPLC-MS. Sci Rep. 2017;7:1291.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Kolter A, Gemeinholzer B. Internal transcribed spacer primer evaluation for vascular plant metabarcoding. Metabarcoding and Metagenomics. 2021;5:e68155.

    Article  Google Scholar 

  74. Fadrosh DW, Ma B, Gajer P, Sengamalay N, Ott S, Brotman RM, et al. An improved dual-indexing approach for multiplexed 16S rRNA gene sequencing on the Illumina MiSeq platform. Microbiome. 2014;2:6.

    Article  PubMed  PubMed Central  Google Scholar 

  75. Ribas TFA, Sales JB, de Boer L, Anmarkrud H, Oliveira JA, Laux RRM. Unexpected diversity in the diet of Doryteuthis sanpaulensis (Brakoniecki, 1984) (Mollusca: Cephalopoda) from the southern brazilian sardine fishery identified by metabarcoding. Fish Res. 2021;239:105936.

    Article  Google Scholar 

  76. Zhang J, Kobert K, Flouri T, Stamatakis A. PEAR: a fast and accurate Illumina paired-end reAd mergeR. Bioinformatics. 2014;30:614–20.

    Article  CAS  PubMed  Google Scholar 

  77. Boyer F, Mercier C, Bonin A, Le Bras Y, Taberlet P, Coissac E. Obitools: a unix-inspired software package for DNA metabarcoding. Mol Ecol Resour. 2016;16:176–82.

    Article  CAS  PubMed  Google Scholar 

  78. Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1.

    Article  CAS  PubMed  Google Scholar 

  79. Edgar RC. UNOISE2: improved error-correction for Illumina 16S and ITS amplicon sequencing 2016:081257.

  80. Camacho C, Madden T, Tao T, Agarwala R, Morgulis A. BLAST® Command Line Applications user Manual. National Center for Biotechnology Information (US); 2008.

  81. Booker A, Johnston D, Heinrich M. Value chains of herbal medicines—Research needs and key challenges in the context of ethnopharmacology. J Ethnopharmacol. 2012;140:624–33.

    Article  PubMed  Google Scholar 

  82. Wallace S, Carrier DJ, Beitle RR, Clausen EC, Griffis CL, HPLC-UV, LC-MS-MS. Characterization of silymarin in milk thistle seeds and corresponding products. J Nutraceuticals Funct Med Foods. 2003;4:37–48.

    Article  Google Scholar 

  83. Petrásková L, Káňová K, Biedermann D, Křen V, Valentová K. Simple and Rapid HPLC separation and quantification of flavonoid, flavonolignans, and 2,3-dehydroflavonolignans in silymarin. Foods. 2020;9:116.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  84. Fiehn O. Metabolomics–the link between genotypes and phenotypes. Plant Mol Biol. 2002;48:155–71.

    Article  CAS  PubMed  Google Scholar 

  85. Goodacre R, Vaidyanathan S, Dunn WB, Harrigan GG, Kell DB. Metabolomics by numbers: acquiring and understanding global metabolite data. Trends Biotechnol. 2004;22:245–52.

    Article  CAS  PubMed  Google Scholar 

  86. Villas-Bôas SG, Rasmussen S, Lane GA. Metabolomics or metabolite profiles? Trends Biotechnol. 2005;23:385–6.

    Article  CAS  PubMed  Google Scholar 

  87. Alonso A, Marsal S, Julià A. Analytical methods in untargeted metabolomics: state of the art in 2015. Front Bioeng Biotechnol 2015;3.

  88. Yang Q, Zhang A, Miao J, Sun H, Han Y, Yan G, et al. Metabolomics biotechnology, applications, and future trends: a systematic review. RSC Adv. 2019;9:37245–57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  89. Song Q, Zhang A, Yan G, Liu L, Wang X. Technological advances in current metabolomics and its application in tradition chinese medicine. RSC Adv. 2017;7:53516–24.

    Article  CAS  Google Scholar 

  90. Atanasov A, Zotchev S, Dirsch V, Orhan I, Banach M, Rollinger J, et al. Natural products in drug discovery: advances and opportunities. Nat Rev Drug Discovery. 2021;20:1–17.

    Article  CAS  Google Scholar 

  91. Ahmad S, Katiyar CK, Ulrich-Merzenich GS, Mukherjee PK. Editorial: Metabolomics and ethnopharmacology in the development of herbal and traditional medicine. Front Pharmacol 2022;13.

  92. Zhao J, Wang M, Saroja SG, Khan IA. NMR technique and methodology in botanical health product analysis and quality control. J Pharm Biomed Anal. 2022;207:114376.

    Article  CAS  PubMed  Google Scholar 

  93. Abraham EJ, Kellogg JJ. Chemometric-guided approaches for profiling and authenticating botanical materials. Front Nutr. 2021;0.

  94. Socaciu C. From phytochemistry to metabolomics: eight decades of research in plant and food science. Studia UBB Chemia. 2019;64:205–24.

    Article  CAS  Google Scholar 

  95. Harrieder E-M, Kretschmer F, Böcker S, Witting M. Current state-of-the-art of separation methods used in LC-MS based metabolomics and lipidomics. J Chromatogr B. 2022;1188:123069.

    Article  CAS  Google Scholar 

  96. Dunn WB, Ellis DavidI. Metabolomics: current analytical platforms and methodologies. TRAC Trends Anal Chem. 2005;24:285–94.

    Article  CAS  Google Scholar 

  97. Lee K-M, Jeon J-Y, Lee B-J, Lee H, Choi H-K. Application of metabolomics to quality control of natural product derived medicines. Biomol Ther (Seoul). 2017;25:559–68.

    Article  CAS  PubMed  Google Scholar 

  98. Alseekh S, Fernie AR. Metabolomics 20 years on: what have we learned and what hurdles remain? Plant J. 2018;94:933–42.

    Article  CAS  PubMed  Google Scholar 

  99. Emwas A-H, Roy R, McKay RT, Tenori L, Saccenti E, Gowda GAN et al. NMR spectroscopy for metabolomics research. Metabolites 2019;9.

  100. Roberts LD, Souza AL, Gerszten RE, Clish CB. Targeted metabolomics. Curr Protoc Mol Biol. 2012;CHAPTER. :Unit30.2.

  101. Billet K, Malinowska MA, Munsch T, Unlubayir M, Adler S, Delanoue G, et al. Semi-targeted metabolomics to validate biomarkers of grape downy mildew infection under field conditions. Plants. (Basel). 2020;9:1008.

    Article  CAS  Google Scholar 

  102. Choi YH, Choi H-K, Hazekamp A, Bermejo P, Schilder Y, Erkelens C, et al. Quantitative analysis of bilobalide and ginkgolides from Ginkgo biloba Leaves and Ginkgo products using. Chem Pharm Bull. 2003;51:158–61.

    Article  CAS  Google Scholar 

  103. Yang SY, Kim HK, Lefeber AW, Erkelens C, Angelova N, Choi YH, et al. Application of two-dimensional nuclear magnetic resonance spectroscopy to quality control of ginseng commercial products. Planta Med. 2006;72:364–9.

    Article  CAS  PubMed  Google Scholar 

  104. Booker A, Zhai L, Gkouva C, Li S, Heinrich M. From traditional resource to global commodities: a comparison of Rhodiola species using NMR spectroscopy-metabolomics and HPTLC. Front Pharmacol 2016;7.

  105. Booker A, Jalil B, Frommenwiler D, Reich E, Zhai L, Kulic Z, et al. The authenticity and quality of Rhodiola rosea products. Phytomedicine. 2016;23:754–62.

    Article  CAS  PubMed  Google Scholar 

  106. Seethapathy GS, Tadesse M, Urumarudappa SKJ, Gunaga V, Vasudeva S, Malterud R. Authentication of Garcinia fruits and food supplements using DNA barcoding and NMR spectroscopy. Sci Rep. 2018;8:10561.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  107. Windarsih A, Rohman A, Swasono RT. Application of 1H-NMR based metabolite fingerprinting and chemometrics for authentication of Curcuma longa adulterated with C. heyneana. J Appl Res Med Aromatic Plants. 2019;13:100203.

    Article  Google Scholar 

  108. Marchev AS, Koycheva IK, Aneva IY, Georgiev MI. Authenticity and quality evaluation of different Rhodiola species and commercial products based on NMR-spectroscopy and HPLC. Phytochem Anal. 2020;31:756–69.

    Article  CAS  PubMed  Google Scholar 

  109. Hebert PDN, Cywinska A, Ball SL, deWaard JR. Biological identifications through DNA barcodes. Proceedings of the Royal Society of London B: Biological Sciences 2003;270:313–21.

  110. Taberlet P, Coissac E, Pompanon F, Brochmann C, Willerslev E. Towards next-generation biodiversity assessment using DNA metabarcoding. Mol Ecol. 2012;21:2045–50.

    Article  CAS  PubMed  Google Scholar 

  111. Epp LS, Boessenkool S, Bellemain EP, Haile J, Esposito A, Riaz T, et al. New environmental metabarcodes for analysing soil DNA: potential for studying past and present ecosystems. Mol Ecol. 2012;21:1821–33.

    Article  CAS  PubMed  Google Scholar 

  112. Coghlan ML, White NE, Murray DC, Houston J, Rutherford W, Bellgard MI, et al. Metabarcoding avian diets at airports: implications for birdstrike hazard management planning. Invest Genet. 2013;4:27.

    Article  Google Scholar 

  113. Nielsen UN, Wall DH. The future of soil invertebrate communities in polar regions: different climate change responses in the Arctic and Antarctic? Ecol Lett. 2013;16:409–19.

    Article  PubMed  Google Scholar 

  114. Soininen EM, Zinger L, Gielly L, Bellemain E, Bråthen KA, Brochmann C, et al. Shedding new light on the diet of norwegian lemmings: DNA metabarcoding of stomach content. Polar Biol. 2013;36:1069–76.

    Article  Google Scholar 

  115. Dormontt EE, van Dijk K, Bell KL, Biffin E, Breed MF, Byrne M et al. Advancing DNA barcoding and metabarcoding applications for plants requires systematic analysis of herbarium collections—an australian perspective. Front Ecol Evol 2018;6.

  116. Raclariu-Manolică AC, de Boer HJ. Chapter 8 - DNA barcoding and metabarcoding for quality control of botanicals and derived herbal products. In: Mukherjee PK, editor. Evidence-Based Validation of Herbal Medicine (Second Edition), Elsevier; 2022, p. 223–38.

  117. Coghlan ML, Haile J, Houston J, Murray DC, White NE, Moolhuijzen P, et al. Deep sequencing of plant and animal DNA contained within traditional chinese medicines reveals legality issues and health safety oncerns. PLOS Genet. 2012;8:e1002657.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  118. Coghlan ML, Maker G, Crighton E, Haile J, Murray DC, White NE, et al. Combined DNA, toxicological and heavy metal analyses provides an auditing toolkit to improve pharmacovigilance of traditional chinese medicine (TCM). Sci Rep. 2015;5:17457.

    Article  CAS  Google Scholar 

  119. Cheng X, Su X, Chen X, Zhao H, Bo C, Xu J, et al. Biological ingredient analysis of traditional chinese medicine preparation based on high-throughput sequencing: the story for Liuwei Dihuang Wan. Sci Rep. 2014;4:5147.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  120. Ivanova NV, Kuzmina ML, Braukmann TWA, Borisenko AV, Zakharov EV. Authentication of herbal supplements using next-generation sequencing. PLoS ONE. 2016;11:e0156426.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  121. Raclariu AC, Mocan A, Popa MO, Vlase L, Ichim MC, Crisan G, et al. Veronica officinalis product authentication using DNA metabarcoding and HPLC-MS reveals widespread adulteration with Veronica chamaedrys. Front Pharmacol. 2017;8.

  122. Seethapathy GS, Raclariu-Manolica A-C, Anmarkrud JA, Wangensteen H, de Boer HJ. DNA metabarcoding authentication of ayurvedic herbal products on the european market raises concerns of quality and fidelity. Front Plant Sci. 2019;10.

  123. Robasky K, Lewis NE, Church GM. The role of replicates for error mitigation in next-generation sequencing. Nat Rev Genet. 2014;15:56–62.

    Article  CAS  PubMed  Google Scholar 

  124. Ficetola GF, Pansu J, Bonin A, Coissac E, Giguet-Covex C, De Barba M, et al. Replication levels, false presences and the estimation of the presence/absence from eDNA metabarcoding data. Mol Ecol Resour. 2015;15:543–56.

    Article  CAS  PubMed  Google Scholar 

  125. Pawluczyk M, Weiss J, Links MG, Egaña Aranguren M, Wilkinson MD, Egea-Cortines M. Quantitative evaluation of bias in PCR amplification and next-generation sequencing derived from metabarcoding samples. Anal Bioanal Chem. 2015;407:1841–8.

    Article  CAS  PubMed  Google Scholar 

  126. Novak J, Grausgruber-Gröger S, Lukas B. DNA-based authentication of plant extracts. Food Res Int. 2007;40:388–92.

    Article  CAS  Google Scholar 

  127. Särkinen T, Staats M, Richardson JE, Cowan RS, Bakker FT. How to open the treasure chest? Optimising DNA extraction from herbarium specimens. PLoS ONE. 2012;7:e43808.

    Article  PubMed  PubMed Central  Google Scholar 

  128. Costa J, Amaral JS, Fernandes TJR, Batista A, Oliveira MBPP, Mafra I. DNA extraction from plant food supplements: influence of different pharmaceutical excipients. Mol Cell Probes. 2015;29:473–8.

    Article  CAS  PubMed  Google Scholar 

  129. Piñol J, Mir G, Gomez-Polo P, Agustí N. Universal and blocking primer mismatches limit the use of high-throughput DNA sequencing for the quantitative metabarcoding of arthropods. Mol Ecol Resour. 2015;15:819–30.

    Article  CAS  PubMed  Google Scholar 

  130. Giguet-Covex C, Pansu J, Arnaud F, Rey P, Griggo C, Gielly L, et al. Long livestock farming history and human landscape shaping revealed by lake sediment DNA. Nat Communications; Lond. 2014;5:3211.

    Article  CAS  Google Scholar 

  131. Dobrovolny S, Uhlig S, Frost K, Schlierf A, Nichani K, Simon K, et al. Interlaboratory validation of a DNA metabarcoding assay for mammalian and poultry species to detect food adulteration. Foods. 2022;11:1108.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  132. Sgamma T, Lockie-Williams C, Kreuzer M, Williams S, Scheyhing U, Koch E, et al. DNA barcoding for industrial quality assurance. Planta Med. 2017;83:1117–29.

    Article  CAS  PubMed  Google Scholar 

  133. Raclariu-Manolica AC, Mauvisseau Q, De Boer H. Horizon scan of DNA-based methods for quality control and monitoring of herbal preparations. Front Pharmacol 2023;14.

  134. Durazzo A, Sorkin BC, Lucarini M, Gusev PA, Kuszak AJ, Crawford C et al. Analytical challenges and metrological approaches to ensuring dietary supplement quality: International Perspectives. Frontiers in Pharmacology 2022;12.

  135. Raclariu AC, Ţebrencu CE, Ichim MC, Ciupercǎ OT, Brysting AK, de Boer H. What’s in the box? Authentication of Echinacea herbal products using DNA metabarcoding and HPTLC. Phytomedicine. 2018;44:32–8.

    Article  CAS  PubMed  Google Scholar 

Download references


We are grateful to our collaborators from Agricultural Research and Development Station Secuieni, Neamt County (RO), especially to Dr. Eng. Trotuș Elena, Dr. Eng. Oana Mirzan, and Dr. Eng. Naie Margareta, as well as to Dr. Carmen Elena Tebrencu from the Research Center for Medicinal Plant and Processing Plantavorel Piatra Neamt (RO), for providing access to samples of Silybum marianum (L.) Gaertn used for scientific analysis. We are grateful for logistical support to the head engineers Audun Schrøder-Nielsen and Birgitte Lisbeth Graae Thorbek, and to the laboratory manager, Jarl Andreas Anmarkrud, from the DNA laboratory at Natural History Museum, University of Oslo (NO).


Open access funding provided by University of Oslo (incl Oslo University Hospital). This work was supported by a grant of the Romanian Ministry of Research, Innovation and Digitalization, CNCS-UEFISCDI, project number PN-III-P1-1.1-PD-2019-0522, within PNCDI III.

Author information

Authors and Affiliations



Conceptualization, ACRM; methodology, ACRM, QM, HdB, CS; software, ACRM, QM, RP, CS; formal analysis, ACRM, QM, HdB, CS; investigation, ACRM; resources, ACRM, HdB, CS; data curation, ACRM, HdB, CS; writing—original draft preparation, ACRM; writing—review and editing, ACRM. QM, RP, HdB, and CS; visualization, ACRM, and CS; supervision, ACRM; project administration, ACRM; funding acquisition, ACRM.

Corresponding author

Correspondence to Ancuța Cristina Raclariu-Manolică.

Ethics declarations

Competing interests

The authors declare no competing interests.

Ethics approval and consent to participate

All the experimental research and field studies on plants (either cultivated or wild), including the collection of plant material, were carried out in accordance with relevant institutional, national, and international guidelines and legislation.

Consent for publication

Not applicable.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Raclariu-Manolică, A.C., Mauvisseau, Q., Paranaiba, R. et al. Authentication of milk thistle commercial products using UHPLC-QTOF-ESI + MS metabolomics and DNA metabarcoding. BMC Complement Med Ther 23, 257 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: