Abstract
Epimedii Folium (EF) is a commonly used traditional Chinese drug that includes many species. In Chinese Pharmacopoeia (2020 edition), EF includes Epimedium brevicornu Maxim, E. Sagittatum (Sieb. et Zucc) Maxim, E. pubescens Maxim, E. koreanum Nakai, and Epimedium Wushanense T.S. Ying. It has been reported that the active constituents in EF are different, resulting in the uneven quality of commercial medicinal materials. To explore the specific differences and make a comprehensive quality evaluation of EF, we established an analytical method to simultaneously detect 45 constituents including 23 flavonoids, 4 phenolic acids, 12 amino acids, 5 nucleosides, and 1 alkaloid in 5 species of EF, based on ultrafast performance liquid chromatography coupled with triple quadrupole-linear ion trap mass spectrometry (UFLC-QTRAP-MS/MS). In addition, orthogonal partial least squares discriminant analysis (OPLS-DA), analysis of variance (ANOVA), and grey correlation analysis (GRA) were used to discriminate and evaluate different species of samples. The results showed that EF could be divided into three categories and EBM was superior to the other four species. This study provides data basis for the comprehensive evaluation and a new perspective on the quality control of EF from different species.
1. Introduction
Epimedii Folium (EF), also known as Xianlingpi and Gangqian, is first recorded in Commentary of Shennong’s Herbal Classic of Materia Medica [1]. The medicinal part of EF is leaf. EF belongs to the family of Epimedium of Berberidaceae, which has the traditional functions of invigorating the kidney and strengthening the yang, strengthening bones and muscles, dispelling wind, and removing dampness [2]. In Pharmacopoeia (2020 edition), there are 5 species of EF, including Epimedium brevicornu Maxim (EBM), E. Sagittatum (Sieb. et Zucc) Maxim (ESM), E. pubescens Maxim (EPM), E. koreanum Nakia (EKN), and E. Wushanense T.S. Ying(EWY) [3]. Among them, EWY, as a plant of the same genus, is listed separately. About 68 species of Epimedium were found in the world, and there were 58 species in China [4]. EF is widely distributed in the regions of central, southwest, northwest, and northeast China [5], and 5 species are included in Pharmacopoeia (2020 version).. In recent years, new species have also been reported, such as E. pseudowushanense B. L. Guo, E. qingchengshanense G. Y. Zhong and Guo et al. [6], and E. myrianthum Stearn [7].
Over the past few decades, the majority of the literature on EF has primarily focused on pharmacological actions identifying [8], chemical constituents [9], and quality control [10]. The research showed that EF contains a variety of bioactive constituents, but not limited to flavonoids [9, 11], phenolic acids [12], volatile oil [13], polysaccharides [14, 15], nucleic acid, amino acids, alkaloids, and lignans [16, 17]. These bioactive constituents play a significant role in the exceptional clinical function, such as antiosteoporosis [18], bone repair [19, 20], hormone-like effects [21], antitumor [22], anti-inflammatory [23], antiaging [24], immunomodulation [15, 25], and treatment of cardiovascular and cerebrovascular diseases [26, 27]. Flavonoid glycoside derivatives with 8-isopentenyl flavonol as the basic mother nucleus, represented by icaritin and icariin, have been proved to be the main substances that exert the above pharmacological effects. Therefore, it is essential to ensure the quality of EF for its remarkable clinical effect.
Nowadays, there are many reports on the analysis and quality assessment of chemical components in EF, such as the use of reversed-phase high-performance liquid chromatography (RP-HPLC) [28], gas chromatography-mass spectrometry (GC-MS) [13], ultraperformance liquid chromatography (UPLC) [9], and rapid resolution liquid chromatography (RRLC) [29]. However, due to the various species and wide origin of this medicinal material, the quality of EF is uneven and the chemical components content is differential [30]. Previous studies have primarily concentrated on the quantitative and qualitative analysis of flavonoids with significant pharmacological effects in a single species. However, no study has used LC-MS technology to perform a comprehensive quantitative analysis of phenolic acids, amino acids, nucleic acids, and alkaloids in EF for quality assessment. As we all know, the components system of traditional Chinese medicine (TCM) is complex, and the efficacy is the result of the integration and regulation of multiple components [31]. Assessing the quality of herbal medicines based on the determination of a single or a few bioactive constituents is biased and may not fully reflect its intrinsic quality. According to the characteristics of multiple active constituents in EF, establishing a method of simultaneous determination of multiple bioactive constituents was necessary, which was scientific and practical for exploring the comprehensive evaluation system of multiple index constituents.
Hence, the aim of this study was to identify and evaluate the quality of EF from different species, including EBM, ESM, EPM, EKN, and EWY. Simultaneous determination of 45 components in 5 species of EF was conducted using ultrafast liquid chromatography-triple quadrupole-linear ion trap tandem mass spectrometry (UFLC-QTRAP-MS/MS). In addition, the extraction conditions are optimized by single-factor analysis and the Box–Behnken design (BBD) along with response surface methodology (RSM). Orthogonal partial least squares discriminant analysis (OPLS-DA) and analysis of variance (ANOVA) were applied to classify, distinguish, and reveal the differential constituents between different species of EF. Then, the quality of EF was assessed using grey association analysis (GRA). This study not only established a reliable and accurate method based on UFLC-QTRAP-MS/MS for the simultaneous determination of multiple bioactive constituents in EF but also provide a basis for the comprehensive evaluation and intrinsic quality control of EF from different species in Pharmacopoeia. This study not only developed a dependable and precise UFLC-QTRAP-MS/MS approach for the concurrent determination of multiple bioactive constituents in EF but also established a foundation for the comprehensive evaluation and intrinsic quality control of various species of EF in Pharmacopoeia.
2. Materials and Methods
2.1. Plant Materials
Plant material was collected from nine provinces of China and North Korea, the domestic region including Gansu, Shanxi, Shaanxi, Sichuan, Guizhou, Chongqing, Hubei, Liaoning, and Jilin. Table 1 details the geographical habitats of each sample. Figure 1 shows different species of EF. All samples were authenticated by Professor Xunhong Liu (Nanjing University of Chinese Medicine, Nanjing, PR China). The Laboratory of Chinese medicine identification, Nanjing University of Chinese Medicine, deposited the plant material.

2.2. Chemicals and Reagents
The standards of serine, glutamic acid, threonine, proline, adenosine, valine, leucine, phenylalanine, caffeic acid, hyperoside, epimedin C, and icariin were purchased from the Institute of Food and Drug Administration of China (Beijing, China). Histidine, glycine, alanine, uracil, 2′-deoxyadenosine, inosine, tyrosine, 2′-deoxyinosine, quercitrin, quercetin, and luteolin were purchased from Yuanye Biotechnology Co., Ltd. (Shanghai, China). Isoquercitrin was purchased from Jiangsu Yongjian Pharmaceutical Technology Co., Ltd. (Jiangsu, China). Astragalin was purchased from Chengdu Desite Biotechnology Co., Ltd. (Chengdu, China). Neochlorogenic acid and cryptochlorogenic acid were purchased from Chengdu Purifa Technology Development Co., Ltd. (Chengdu, China). Magnoflorine, trifolin, afzelin, epimedoside A, epimedin A1, epimedin A, epimedin B, sagittatoside A, epimedin-I, ikarisoside A, icariside-I, icaritin, sagittatoside B, 2′-O-rhamnosyl icariside-II, baohuoside-I, and anhydroicaritin were purchased from Chengdu Alfa Biotechnology Co., Ltd. (Chengdu, China). Lysine and chlorogenic acid were purchased from Baoji Chenguang Biotechnology Co., Ltd. (Baoji, China).
2.3. Preparation of Standard Solutions
A stock solution of reference standards consisting of 45 compounds was prepared using methanol, and their concentrations were as follows: 1.008, 1.00, 1.180, 1.160, 5.295, 1.100, 1.160, 4.010, 2.478, 0.800, 2.928, 1.700, 2.583, 2.185, 2.615, 1.102, 3.346, 1.055, 1.290, 1.830, 1.052, 1.620, 1.094, 0.936, 1.004, 3.470, 1.032, 1.840, 1.035, 1.545, 1.128, 1.360, 2.420, 2.460, 2.715, 1.190, 1.565, 3.130, 1.770, 2.125, 1.505, 0.735, 1.080, 0.765, and 0.153 mg/ml; the stock solution was then diluted with ethanol to generate calibration curves at various concentrations. All solutions were kept at 4°C until LC-MS analysis.
2.4. Preparation of Sample Solutions
Approximately 0.5 g of sample powder from various species of EF were weighed accurately and sonicated with 50% ethanol (5 mL) for 60 minutes, followed by cooling to room temperature. And we supplemented the weight loss with 50% ethanol, shaken well, and filtered. Then, the filtrate was centrifuged at 12000 r/min for 10 min, and the supernatant was filtered through a 0.22 μm microporous membrane and stored at 4°C in a refrigerator.
2.5. Chromatographic and Mass Spectrometric Conditions
The SIL-20A XR system (Shimadzu, Kyoto, Japan) was employed for chromatographic analysis of EF. The chromatographic column used was an XBridge® C18 column (4.6 mm × 100 mm, 3.5 μm) maintained at 30°C, with an injection volume of 2 µL. The mobile phase consisted of a 0.4% formic acid water solution (A) and acetonitrile (B) at a flow rate of 0.8 mL/min, with the following gradient elution: 0∼1 min, 5∼5% B; 1∼14 min, 5∼31% B; 14∼18 min, 31∼41% B; 18∼22 min, 41∼56% B; 22∼26 min, 56∼75% B; 26∼29 min, 75∼83% B; and 29∼32 min, 83∼5% B.
Detection was carried out using an API5500 triple quadrupole-linear ion trap tandem mass spectrometer (AB SCIEX, Framingham, MA, USA) equipped with an electrospray ionization (ESI) source. The following operating parameters were used: ion source temperature, 550°C; nebulizer gas (GS1) flow, 55 psi; auxiliary gas (GS2) flow, 55 psi; curtain gas (CUR) flow, 40 psi; and spray voltage (IS), 4500 V in positive mode and −4500 V in negative mode. Analyte detection was performed using multiple-reaction mode (MRM).
2.6. Validation of the Method
The validity of the approach was confirmed through assessments of its linearity, intraday and interday precision, repeatability, stability, recovery, and matrix effect [32]. The standard curve was established by serial dilution of mixed standards, and the equation for linear regression, coefficient of correlation, and range for linearity was computed. The limit of detection (LOD) and limit of quantification (LOQ) for 45 components were determined at signal-to-noise ratios of 3 and 10. Intraday precision was evaluated by injecting the mixed standard solution six times per day, while interday precision was assessed by injecting triplicates of the solution for three consecutive days. To verify the repeatability, six EF samples were precisely weighed, prepared individually according to the aforementioned optimal conditions, and subsequently analyzed. The stability of the sample was assessed by extracting and analyzing the same solution at 0, 2, 4, 8, 12, and 24-hour intervals, using the previously mentioned chromatographic conditions. The recovery experiment was used to evaluate the accuracy of the method; three standards with different concentration levels were added to the samples with known content, including low concentration (80%), medium concentration (100%), and high concentration (120%). The analysis of spiked samples using UFLC-QTRAP-MS/MS to determine the recovery was conducted in triplicate for each experiment. Matrix effect refers to the interference in the matrix or the enhancement or inhibition of chromatographic signals by blended constituents. The slope comparison method was used to evaluate it. Thus, the matrix effect was computed as the ratio of the slope in the matrix-matched calibration curve to the slope in the solvent standard curve. A slope close to 1.0 indicates a weaker matrix effect.
2.7. Multivariate Statistical Analysis
Following data preprocessing, OPLS-DA was utilized to observe the overall clustering trend of each group, and SIMCA-P 13.0 software (Umetrics AB, Umea, Sweden) was employed to visualize their distribution. The quantitative detection outcomes for the 45 active components were subjected to variance analysis using SPSS Statistics 22.0 software. Additionally, OriginPro 2021b (OriginLab, Northampton, MA, USA) was used to provide a more intuitive visualization of the classification and content changes of the samples. Excel for Mac 2019 (Microsoft Corporation, Seattle, WA, USA) was used to evaluate the quality of various EF species based on their content of 45 active components, using GRA. OriginPro 2021b (OriginLab, Northampton, MA, USA) was employed to generate histograms for all data.
3. Results and Discussion
3.1. Optimization of Extraction Conditions
3.1.1. Single-Factor Experiment
The selection of extraction solvent, extraction time, and solid-liquid ratio (g/ml) has a considerable impact on the constituents extracted from EF. According to the Pharmacopoeia (2020 edition), the content of icariin in leaves should not be less than 0.5%. Therefore, the extraction rate of icariin was selected as the index of optimal extraction conditions in this study. In order to obtain suitable extraction efficiency of icariin, the single-factor tests were performed for ethanol volume fraction (20% methanol/ethanol, 30% methanol/ethanol, 40% methanol/ethanol, 50% methanol/ethanol, 60% methanol/ethanol, 70% methanol/ethanol, and 80% methanol/ethanol, respectively), extraction time (30, 40, 50, 60, 70 min, respectively), and liquid-solid ratio (6 : 1, 8 : 1, 10 : 1, 12 : 1, 14 : 1, 20 : 1, 30 : 1, 1 : 40 : 1, and 50 : 1, respectively). Figure S1 (in the supplementary material) illustrates the impact of ethanol volume fraction, liquid-solid ratio, and extraction time on the extraction yields of icariin.
As shown in Figure S1, the extraction rate of icariin was affected by the concentration of ethanol volume fraction. When the volume fraction of ethanol was 50%, the maximum extraction rate was 4895.64 μg/g. With the increase of ultrasonic extraction time, the extraction of icariin first decreased and then increased. The extraction rate reached 5052.77 μg/g when the extraction time reached 60 min. The liquid-solid ratio had an irregular effect on the extraction rate of icariin, while at the ratio of 1 : 10, the extraction rate reached a peak of 4812.28 μg/g.
3.1.2. Box–Behnken Design along with Response Surface Methodology
After analyzing the outcomes of the aforementioned single-factor experiment, the Box–Behnken design (BBD) was employed in conjunction with response surface methodology (RSM) to design a three-factor, three-level experiment. The extraction rate of icariin was utilized as the response value, with the aim of optimizing the icariin extraction process. The experimental design factors and levels of the response surface method are shown in Table S1, the experimental design and results are shown in Table S2, and the analysis of variance and significance test are shown in Table S3. The results showed that the model can be used to determine the optimal extraction conditions of icariin. Then, the effects of ethanol volume fraction, extraction time, and liquid-solid ratio on the peak area of icariin were analyzed using Design-Expert 8.0.6 software, and results were shown in Figure 2.

Ultimately, using Design-Expert 8.0.6 to solve the regression equation for the data in Table S2, the optimal extraction conditions were obtained when the extraction time was 60.1 min and the liquid-solid ratio was 10 : 1 for 49.41% ethanol, with the predicted value was 3271.00 μg/g. In order to simplify the experimental operation, the above conditions were adjusted as 50% ethanol volume fraction, 10 : 1 liquid-solid ratio, and 60 min extraction time. Under this condition, the average extraction yield was 3841.22 μg/g, which was slightly higher than the theoretical value of 3271.00 μg/g. The experimental results supported the theoretical value, indicating that the optimized conditions were feasible and could be used as extraction conditions for EF.
3.2. Optimization of UFLC Parameters
The UFLC chromatographic conditions, including the column, mobile phase, and column temperature, were optimized to achieve the most favorable separation efficiency for the target constituents in EF and to improve the peak shape. Comparing the separation effects of XBridge® C18 column (4.6 mm × 100 mm, 3.5 μm) and Agilent ZORBAXSB-C18 column (4.6 mm × 250 mm, 5 μm) on 45 index constituents, it was found that the XBridge® C18 column (4.6 mm × 100 mm, 3.5 μm) had relatively high resolution and sensitivity. In addition, four mobile phase systems (water-methanol, water-acetonitrile, water-methanol:acetonitrile (1 : 1), and 0.1%, 0.4%, 0.8% formic acid water solution-acetonitrile, respectively), flow rates (0.3, 0.4, 0.5, 0.7, 0.8, 0.9, and 1.0 ml/min, respectively), and column temperatures (25, 30, 35, and 40°C, respectively) were examined and compared. The final chromatographic conditions were established using gradient elution with 0.4% formic acid as eluent A and acetonitrile as eluent B, at a flow rate of 0.8 ml/min and a column temperature of 30°C.
3.3. Optimization of Mass Spectrometry (MS) Conditions
Both positive and negative scanning modes are used by the mass spectrometry to detect all standard solutions (about 100 ng/ml) for a sensitive and accurate quantitative method. Following multiple experimental tests, it was determined that amino acids, nucleosides, and magnoflorine exhibit strong sensitivity in positive ion mode, while flavonoids and phenolic acids demonstrate significantly greater reactivity in negative ion mode. Constituents of magnoflorine, epimedoside A, and icariin are correspondingly better in negative ion mode than those in positive ion mode. Thus, this study uses both ESI+ and ESI− modes. Although some constituents have similar retention times, they can be accurately quantified based on different precursor and product ions. Figures S2 and S3 show the spectra of the negative and positive ion modes from the mixed standards, respectively. The optimal parameters for the 45 constituents, including retention time (tR), precursor and product ions, declustering potential (DP), and collision energy (CE), are provided in Table 2. The CAS numbers of 45 target constituents are shown in Table S4. Figure 3 shows the multiple-reaction monitoring (MRM) of the 45 constituents.

3.4. Method Validation
All quantitative method validation was performed using the above method. The results of the verification of each method are shown in Table 3. Each standard calibration curve was constructed by plotting the peak area (Y) versus the corresponding concentration (X). All analytes showed good linearity with an appropriate coefficient of determination (r > 0.9990). The ranges of LOD and LOQ range from 0.085∼8.792 ng/ml to 0.283∼29.306 ng/ml, respectively. The relative standard deviations (RSDs) of intraday and interday precision are 0.32%∼5% and 0.62%∼4.99%, respectively. In turn, the RSDs for the repeatability and stability tests were less than 4.94% and 4.81%. The total recovery rate varied from 94.51% to 105.66%, with RSD <4.94%. The recovery of different species of EF and the RSD of recovery for each species are shown in Tables S5 and S6. The slope ratio indicated that the effect of the matrix on analyte ionization was not obvious under optimal conditions. The matrix effect of different species of EF is shown in Table S7.
3.5. Quantitative Analysis of Samples
Information regarding the samples is provided in Table 1. Using the UFLC-QTRAP-MS/MS verification analytical method, 45 constituents, which includes 23 flavonoids, 4 phenolic acids, 12 amino acids, 5 nucleosides, and 1 alkaloid, were simultaneously determined in EF. Table S8 shows quantitative results of 45 constituents. EF samples were all rich in flavonoids, with the content of flavonoid ranging from 5520.32 to 21299.16 μg/g, comprising over 80% of the total analytes analyzed in this study. Icariin, as the quality index component of EF in Pharmacopoeia, ranges from 802.48 to 6419.17 μg/g, with a percentage of 0.08%∼0.64%. In particular, the flavonol glycoside derived from 8-isoprenyl flavonol as the basic parent nucleus unique to EF accounted for 99.23% of the total flavonoid content. The flavonoid level varied from sample to sample, with the total flavonoid in the EBM being the highest of all samples. The general content of phenolic acid was 11.50∼3039.92 μg/g, of which the substance of chlorogenic acid was the highest, accounting for more than 40.19% of the total phenolic acid. There were obvious differences in the phenolic acid content in different species, and the phenolic acid content in EKN was significantly higher than that in others. The total content of amino acids was 202.26∼2201.79 μg/g, of which proline, glutamic acid, and phenylalanine were relatively high. Besides, the level of amino acid in EBM and ESM was significantly higher than that of the other three species. The total content of nucleosides was 35.28∼283.65 μg/g, accounting for 0.21%∼1.56% of the total content of all constituents. The level of nucleic acid in ESM and EBM was higher than that of other species. The only alkaloid detected was magnoflorine, with a total content of 470.45∼2165.07 μg/g; its content in ESM was significantly higher than that in other species. The content of flavonoids, phenolic acids, amino acids, nucleosides, and magnoflorine in EF is shown in Figure 4.

3.6. OPLS-DA of Samples
To begin with, the quality differentiation among samples of varied EF species was achieved through the application of principal component analysis (PCA). However, the PCA results showed that samples of different species could not be completely distinguished, and it was deemed unsuitable for establishing a foundation to differentiate and assess the quality of EF. Therefore, orthogonal partial least squares discriminant analysis (OPLS-DA) was used. This analytical method is a supervised discriminant analysis statistical method, which can maximize the variation of between-groups and minimize the separation of within-groups. Using OPLS-DA, a model of the relationship between the content of the target constituents and the sample class was successfully established to predict all samples. Figure 5 shows the OPLS-DA scoring plot. The samples of different species of EF were obviously divided into three categories. EBM was the category I, three species, EPM, EKN, and EWY, were grouped into category II, and ESM was the category III. This outcome indicated the presence of disparities in the bioactive constituents of EBM, ESM, and other three species, which was also consistent with the results of HCA. Furthermore, the statistical parameters of OPLS-DA R2X (cum), R2Y (cum), and Q2 (cum) were 0.824, 0.934, and 0.794, respectively, indicating that the model had good reproducibility and predictability. When a variable’s VIP is greater than 1, it signifies that the variable plays a significant role in the classification of the samples. As depicted in Figure 6, twenty constituents were identified to have a significant impact on the cluster based on their VIP values, including anhydroicaritin, epimedin-I, cryptochlorogenic acid, epimedin C, 2′-O-rhamnosyl icariside-II, isoquercitrin, alanine, leucine, epimedin A, epimedin B, luteolin, quercetin, icaritin, sagittatoside A, epimedoside A, tyrosine, icariside-I, ikarisoside A, lysine, and caffeic acid.


3.7. ANOVA of Samples
Analysis of variance (ANOVA) is a technique for analyzing the way in which the mean of a variable is affected by different types and combinations of factors, which can be used to compare any number of groups or treatments [33]. In this study, the content of bioactive constituents was analyzed using ANOVA followed by the least significant difference (LSD) test (assuming equal variance) or Tamhane’s test (not assuming equal variance) to evaluate the changes of 45 constituents in 5 EF species. The values with values less than 0.05 were deemed as having significant differences. As shown in Table 4, among the 45 constituents, there were 5 constituents, including proline, astragalin, trifolin, epimedin A1, and sagittatoside B, showed no significant differences among different varieties (), which indicates that the content of the above five constituents was consistent among different species without significant differences. The other 40 constituents showed significant differences among different varieties (). Figure 7 shows the box plot of 45 constituents’ content in 5 species of EF. It was worth noting that there were 19 constituents, including 10 flavonoids, 8 amino acids, and uracil, which were significantly higher in EBM than other species. The highest content of 14 constituents was found in ESM, including 7 flavonoids, 4 amino acids, 2 nucleosides, and magnoflorine. The content of five constituents in EPM was significantly higher than that in other species, including four flavonoids and proline. The content of 4 phenolic acids, epimedin-I, and 2-deoxyinosine was the highest in EKN, while the content of 2′-O-rhamnosyl icariside-II in EWY was the highest.

(a)

(b)

(c)
3.8. GRA of Samples
Grey Correlation Degree Analysis (GRA) is an influence measure in grey system theory that analyzes the uncertain relationship between one major factor in a given system and all other factors. Therefore, according to the content of 45 bioactive constituents, a comprehensive GRA evaluation of different species of EF was performed. Table 4 presents the grey composite evaluation values (ri) and quality rankings obtained from the results. The relative correlation degree reflects the association between the component content and the samples, with those having higher relative correlation considered to be of superior quality.
Table 5 shows the grey composite evaluation values (ri) and quality rankings. Overall, the quality of EBM ranked first among the 5 species, which was followed by ESM. While the interspecific quality difference of EPM, EKN, and EWY was not significant, and the overall quality of the three species was inferior to that of EBM and ESM. In addition, there were also differences in the quality of samples of EF of the same species, which may be related to differences in the growing environment, climate, and harvest time. For example, in EBM, samples from Longnan, Gansu Province, were better than those from other sources. It can also be seen that the ri value of each sample varies greatly, with a maximum value of 22.41%, which effectively differentiates the quality among the 5 species of EF. In summary, GRA can successfully assess the quality of different species of EF based on the content of multiple constituents.
3.9. Discussion
As mentioned above, EF is widely used in clinical practice with significant pharmacological effects [18, 19, 21, 23]. Nonetheless, there are a limited number of comprehensive studies on the quality of EF both domestically and internationally, particularly articles that provide a comprehensive analysis of the bioactive constituents of various species. The majority of studies have concentrated on the quantitative and qualitative analysis of flavonoids with significant pharmacological effects and high content in a single species. However, it is widely recognized that the efficacy of TCMs is the outcome of the integrated adjustment of multiple components [31, 32], emphasizing the need for ensuring the quality of TCMs to achieve their remarkable clinical effects. Yin et al. [34] established a method based on UFLC-QTRAP-MS/MS technology for simultaneous determination of multiple active components in Abelmoschi Corolla, which provided an objective method reference for comprehensive evaluation and comprehensive control of Abelmoschi Corolla internal quality. Wei et al. [35] explored the distribution patterns of multiple active constituents in the pericarp, stalk, and seed of Forsythiae Fructus, which offers a scientific foundation for further exploration of the quality formation mechanism of Forsythiae Fructus. Overall quality control of TCMs is more in line with the characteristics of complex system of TCMs. Hence, given the characteristics of multiple active constituents in EF, the method for simultaneous determination of multiple active constituents was established, which was of practical and scientific significance for exploring the comprehensive evaluation system of multiple marker components.
This study established an effective and dependable UFLC-QTRAP-MS/MS method, which can simultaneously determine 23 flavonoids, 12 amino acids, 4 phenolic acids, 5 nucleosides, and 1 alkaloid in 40 different samples of EF. Single-factor analysis and the Box–Behnken design (BBD) along with response surface methodology (RSM) were used to optimize the extraction conditions (Figure 2). OPLS-DA results showed that EBM, ESM, EPM, EWY, and EKN can be clearly divided into 3 categories (Figure 5). The VIP values indicated that there were significant differences in 20 bioactive constituents, which can be thought of as chemical markers that distinguish different species of EF (Figure 6). The results of ANOVA showed that the content of 40 of the 45 target constituents was significantly different among different species (). In addition, based on the ri value of the samples in the GRA analysis, the overall quality of the EBM was optimal, followed by ESM. While the interspecific quality difference of EPM, EKN, and EWY was not significant, and the overall quality of the three species was inferior to that of EBM and ESM. This suggests that the species is indeed an important factor influencing the accumulation of bioactive constituents. However, the GRA analysis found that there were also significant differences between samples of the same EF species, such as the content of total phenolic acids in different samples of EWY was 11.50–3068.36 μg/g, which may be related to external factors such as the growing environment and harvest time. Some articles have also reported on this [12, 36], but there are still shortcomings of researching single species, which cannot form a comprehensive evaluation system.
In Pharmacopoeia (2020 edition), the content of icariin is used as the index component for quality identification of EF, and it is required that icariin should not be less than 0.50%. However, in our research, the content determination results showed that only 5% of the samples met the Pharmacopoeia standard, and the content of icariin was significantly different among different species. The same problem was found in the previous articles [37, 38]. Therefore, in order to further establish a more comprehensive quality evaluation system for EF, multiple investigations should be carried out. Besides the detection of icariin, the constituent’s determination may be added as new indexes for the quality control of EF. In our research, in addition to icariin, epimedin A, epimedin B, epimedin C, sagittatoside A, and magnoflorine also had high content, with the maximum content of 0.24%, 3%, 0.4%, 0.43%, and 0.2%, respectively. In some samples, the content of epimedin C and sagittatoside A was higher than that of icariin. Also, the above constituents have been proved to have obvious pharmacological effects [39–41]. Therefore, in the future research and application, the composition determination may be considered as new quality control indexes of EF.
Compared with previous studies, our research not only analyzed the flavonoids with obvious pharmacological effects but also carried out quantitative analysis of magnoflorine, phenolic acids, amino acids, and nucleosides, which filled the blank of quantitative analysis of constituents with nonobvious pharmacological effects in EF [42]. In addition, the UFLC-QTRAP-MS/MS method was faster than the previous HPLC method, and mass spectrometer detector adopted the positive and negative ion mode for simultaneous determination, which was more sensitive than the conventional detector and more accurate in determining the molecular weight of the target components [43].
It has been reported that China is the center of modern differentiation and diversity of Epimedium. A total of 68 species have been published and 58 from China, among which 57 are unique to China. In addition to the 5 species included in Pharmacopoeia as EF, more than 10 species are commonly used in the folk [44]. Therefore, in the future research, for the problem of uneven quality of commercial medicinal materials on the market, the established method can be applied to the quality evaluation of more species of EF, which can establish a more perfect quality control and evaluation system of EF. In short, the combination of the established UFLC-QTRAP-MS/MS method and multivariate statistical analysis furnished fundamental insights into the quality assessment of EF from diverse species.
4. Conclusion
We established an efficient and reliable UFLC-QTRAP-MS/MS method to simultaneously determinate 23 flavonoids, 12 amino acids, 4 phenolic acids, 5 nucleosides, and 1 alkaloid in 40 different samples of EF. Single-factor analysis and the Box–Behnken design (BBD) along with response surface methodology (RSM) were used to optimize the extraction conditions. Multivariate statistical analysis made use of OPLS-DA, ANOVA, and GRA to classify and distinguish the differential constituents between different species of EF. The experimental results showed that the overall quality of the EBM was optimal, followed by ESM while the interspecific quality difference of EPM, EKN, and EWY was not significant and the overall quality of the three species was inferior to that of EBM and ESM. This study not only established a reliable and accurate method for simultaneous determination of multiple bioactive components in EF by UFLC-QTRAP-MS/MS but also provided a basis for the comprehensive evaluation and intrinsic quality control of EF from different species.
Data Availability
The data that support the findings of this study are included in the supplementary material of this article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Authors’ Contributions
J.X., C.C., and X.L. conceived and designed the experiments; J.X., Y.Z., H.C., J.Y., S.Y., and Y.S. performed the experiments; J.X., J.Y., W.Y., N.W., and Y.S. analyzed the data; J.X. and H.C. wrote the manuscript; Z.C., C.C., X.L., and J.C. revised the manuscript; X.L. and L.Z. provided the samples of Epimedii Folium. All authors read and approved the final manuscript.
Acknowledgments
This research was supported by the General Project of Natural Science Research in Universities of Jiangsu Province (20KJD360001).
Supplementary Materials
Supplementary information is available online. Figure S1: effects of ethanol concentration, solid-liquid ratio, and extraction time on extraction yields of icariin (extraction yield (%) = weight of analyte (mg)/weight of dried sample (g) × 100). Figure S2: the spectra of the negative ion modes from the mixed standards. The numbers denoted are the same as those in Table 2. Figure S3: the spectra of the positive ion modes from the mixed standards. The numbers denoted are the same as those in Table 2. Table S1: experimental design factors and levels of response surface experiment. Table S2: experimental design and results of Box–Benhnken. Table S3: variance analysis and significance test of response surface experiment. Table S4: the CAS numbers of 45 target constituents. Table S5: the recovery of different species. Table S6: the RSDs of recovery for each specie. Table S7: the matrix effect of different species of EF. Table S8: the content of 45 constituents in samples (μg/g, n = 3). (Supplementary Materials)