Abstract

Licorice is a commonly used traditional Chinese medicine and natural sweetening agent, rich in numerous bioactive compounds. Moreover, it is one of the oldest and most frequently employed folk medicines in both eastern and western countries. It is prescribed for the treatment of asthma, fever, and cough. However, with the increasing demand of licorice, its quality and safety become the important issue. The content in licorice varies significantly in materials from different geographical origins. In this study, a reasonable and feasible evaluation method for the quality assessment of licorice was developed based on the analysis of high-performance liquid chromatography (HPLC) fingerprint, combined with the quantitative analysis of multicomponents by single marker (QAMS) method. Glycyrrhizic acid was selected as the internal reference substance, and ten components were simultaneously determined based on relative correction factors. The contents of eleven components in 21 batches of licorice were determined by the QAMS and the ESM (external standard method); there was no significant difference by comparison of the quantitative results between the QAMS and the ESM method; the cosine value (Cir > 0.9999) confirmed the consistency of the two methods. According to the outcomes of 21 batches of licorice samples, the contents of the eleven components were used for further chemometric analysis. All of the samples of licorice from various geographical origins were divided into five categories based on hierarchical cluster analysis, which indicated the crucial influence of geographical origins on licorice. This study showed that QAMS combined with HPLC fingerprint and chemometrics methods could effectively control the quality of licorice. Hence, QAMS is a feasible and promising method for promoting the quality control standardization process of herbal medicines.

1. Introduction

Licorice is the dried root and rhizome of Glycyrrhiza uralensis Fisch., Glycyrrhiza inflata Bat., and Glycyrrhiza glabra L., used as an herbal medicine and natural sweetener in ancient Egyptian, Greece, and China [1, 2]. It is a “necessary herbal medicine” in traditional Chinese medicine (TCM) and almost constitutes about 60% of all TCM prescriptions [3]. This botanical drug has been extensively used to treat the allergic-inflammatory disease, gastrointestinal problems, cardiovascular disease, cancer, and bladder and kidney ailments [4]. Modern pharmacological studies verified that it has the efficacy of the anti-inflammatory, anticancer, antibacterial, antispasmodic, antiallergic, antioxidative, anti-AIDS, and so on [5].

At present, more than 400 compounds have been isolated and identified from licorice. The flavones and saponins are the main representative components in licorice. With the continuous expansion of the licorice market, the quality control of licorice is of increasing importance. Only glycyrrhizic acid is used for the indicators to assess the quality of licorice in European Pharmacopoeia (EP 9.0) [6], United States Pharmacopeia (USP 42) [7], and Japanese Pharmacopoeia (JP 17) [8]. Other important ingredients have been ignored, resulting in insufficient information on the integral chemical composition of licorice. In addition, Chinese Pharmacopeia (2020 version) records that liquiritin (content should not be less than 0.50%) and glycyrrhizic acid (content should not be less than 2.00%) were selected as quality control indicators [9]. However, each Chinese herbal medicine is a comprehensive complex with multiplicity of active components, contributing to pharmacological effects through synergy and mutual effect based on the theory of TCM [10]. Therefore, it is inadequate to use a few components to control the quality of complex plant products and traditional Chinese medicine. According to phytochemical studies, liquiritin, isoliquiritin apioside, isoliquiritin, neoisoliquiritin, isoliquiritigenin, formononetin, licochalcone C, licochalcone A, glabridin, and 18β-glycyrrhetinic acid widely exist in licorice [11]. And their pharmacological effects are significant, which are closely related to the efficacy of licorice. For example, liquiritin and glycyrrhizic acid could treat the inflammation-related disorders [12]; neoisoliquiritin and formononetin inhibit inflammatory cytokines generation and improve antioxidant activity [13, 14]; isoliquiritigenin, licochalcone C, licochalcone A, and glabridin have the efficacy of anticancer [15, 16]; isoliquiritin has therapeutic efficacy on neuropathic hypersensitivity [2]. Therefore, selecting these 11 components as new quality control indicators can comprehensively characterize the chemical properties, medicinal functions, and intrinsic quality of licorice.

In many studies, the ESM has been used to analyze the concentrations of multiple components. However, this method has higher requirements of reference standard which requires more time and high experimental cost to isolation and purification [17]. As an alternative method, quantitative analysis of multicomponents by single marker (QAMS) only needs a single reference standard to simultaneously analyze the contents of multicomponents, which is more available and efficient for the quality control (QC) [18]. This method could reduce not only the cost but also the detection time of experiment [19]. Additionally, QAMS has been widely used in QC to effectively ensure the safety of herbal or botanical products effectively. And it has gradually been recognized by researchers and successfully employed for the analysis of medicinal herb [20]. Since this method was firstly proposed by Dr. Wang in China, it has received more attention and acceptance in the world. And in American Herbal Pharmacopoeia, it is even applied for the quantification study of 20 variants of plants [21].

In recent years, the chromatographic fingerprint analysis has also quite been recommended by academia community as a useful method for the quality control of medicinal plants and herbal products. Furthermore, this method has been accepted as a strategy for the quality assessment of herbal medicines and preparations by the US Food and Drug Administration, State Food and Drug Administration of China, and the European Medicines Agency [22]. Fingerprint technology emphasizes overall characterization, in order to distinguish some substitutes and adulterants from the real original medicine by chromatographic fingerprinting based on the presence or absence of a limited number of peaks. Therefore, we combine characteristic fingerprints with QAMS, which can more effectively and comprehensively improve the quality evaluation of licorice.

In this study, we established a comprehensive and effective quality evaluation method for licorice, which integrated HPLC fingerprint characterization technology, combined with the QAMS method, and chemometric analysis. This strategy not only provides comprehensive chemical information but also can quickly and reliably determine various components and evaluate quality consistency. This innovative method could be beneficial for simplifying the processes in the quality control of licorice and promoting the quality control standardization process of herbal medicines.

2. Materials and Methods

2.1. Plant Materials

21 batches of dried licorice (Table 1 and Figure 1) were collected from different areas (including Ningxia, Gansu, Inner Mongolia, Xinjiang, and Heilongjiang) of China and identified by Prof. Zhang Xinhui, Ningxia Medical University, Yinchuan, China.

2.2. Reagents and Chemicals

Acetonitrile (purity 100%) of HPLC grade was from Thermo Fisher Scientific Company. Glacial acetic acid (purity 100%, HPLC grade) was acquired from Kaixin Chemical Industry Co. Ltd. (Tianjin, China). Methanol and ethanol were analytical grade. These reference substances of liquiritin (L1), isoliquiritin apioside (L2), isoliquiritin (L3), neoisoliquiritin (L4), isoliquiritigenin (L5), formononetin (L6), licochalcone C (L7), licochalcone A (L8), glabridin (L9), glycyrrhizic acid (L10), and 18β-glycyrrhetinic acid (L11) were purchased from Yuanye Biotechnology Co., Ltd. (Shanghai, China). The purity of each compound was verified to be greater than 98% by high-performance liquid chromatography (HPLC) using the peak area normalization method. The chemical structures of these compounds are shown in Figure 2. All other chemicals used were of analytical grade.

2.3. Preparation of the Sample Solution

Powder sample of licorice was accurately weighed (0.2 g) and added in 25 mL of 70% ethanol in 50 mL triangular flask, sonicated for 30 min, and allowed to cool to room temperature. Additional 70% ethanol was complemented for weightlessness. All of the sample solutions were filtered through a 0.22 µm Millipore filter prior to HPLC analysis.

2.4. Preparation of the Standard Solution

The reference solution of licorice was prepared by accurately dissolving weighed samples of each compound in 70% methanol, making a mixture of 43.85 μg/mL of liquiritin, 40.77 μg/mL of isoliquiritin, 36.92 μg/mL of neoisoliquiritin, 40 μg/mL of isoliquiritin apioside, 41.54 μg/mL of isoliquiritigenin, 30.77 μg/mL of formononetin, 30.07 μg/mL of glabridin, 33.07 μg/mL of licochalcone A, 30.77 μg/mL of licochalcone C, 185 μg/mL of 18β-glycyrrhetinic acid, and 240 μg/mL of glycyrrhizic acid, mixed evenly. All the standard solutions were stored in a refrigerator at 4°C, before use.

2.5. Chromatographic Procedures

Chromatographic analysis was performed by using an Agilent 1260 series system (Agilent Technologies, MA, USA) equipped with an online degassing system, an autosampler, and a thermostated column compartment. Agilent ZORBAX SB-C18 Column (4.6 × 250 mm, 5 μm) (Agilent Technologies, USA) was used for chromatographic separation of target compounds. The mobile phase was made up of 0.2% glacial acetic acid (A) and acetonitrile (B) with the following optimized gradient elution process: 0∼15 min, 15%B; 15∼25 min, 15∼20%B; 25∼70 min, 20∼50%B; 70∼90 min, 50∼70%B; 90∼95 min, 70∼15%B. The detection wavelength was 254 nm. The column temperature was set at 30°C. The flow velocity was 1.0 mL/min with a sample injection volume of 20 µL.

2.6. Method Validation

Eleven reference compounds were accurately weighted and separately dissolved in ethanol to obtain appropriate concentrations. The stock solutions were serially diluted with ethanol to a certain concentration range to obtain a calibration curve. The linear calibration curves were constructed by six different concentrations of mixed standard solutions. The limits of detection (LODs) and limits of quantification (LOQs) were determined at signal-to-noise ratio of 3 and 10, respectively. The evaluation of precision was examined from six-time analysis of the licorice sample placed in the autosampler. And the stability test was conducted by a six-time analysis of the same licorice sample at different time points (0, 3, 6, 9, 12, 18, and 24 h) after solution preparation at room temperature (25 ± 5°C). The repeatability was evaluated according to the assay of six duplicates. Relative standard deviations (RSDs) of relative peak areas (RPAs) of all the analytes were used for evaluation. As to the accuracy investigation, a recovery test was conducted with conventional standard addition method at different concentration following equation (1):

2.7. QAMS Quality Evaluation Model for Licorice
2.7.1. RCF Calculation and Its Tolerance Investigation

The calculation methods of relative correction factor (RCF) have been previously reported. Firstly, glycyrrhizic acid is used as the internal standard, and a multipoint method (equation (2)) was used to calculate the relative correction factor of liquiritin, isoliquiritin, neoisoliquiritin, isoliquiritin apioside, isoliquiritigenin, formononetin, glabridin, licochalcone A, licochalcone C, and 18β-glycyrrhetinic acid. The reproducibility of the correction factors was investigated at different column, using the different instruments. The relative correction factor is calculated aswhere is the peak area of the internal reference substance, is the peak area of sample i, is the concentration of the internal reference substance, and is the concentration of sample i.

2.7.2. Assessment of the ESM and the QAMS Method

In order to verify the feasibility of QAMS, the validated traditional ESM and QAMS methods were employed to evaluate the similarity of the 21 batches of licorice samples from different producing regions. For ESM, the content of the eleven components was determined by eleven reference standards (L1–L11), whereas for QAMS, the results were based on the properties of the calculation of RCF, the content of the selected internal substance (L10) was determined to be similar to ESM, and then the contents of the other ten components were calculated based on the RCF between analytes and the internal substance [23]. In this study, the contents of liquiritin, isoliquiritin, neoisoliquiritin, isoliquiritin apioside, isoliquiritigenin, formononetin, glabridin, licochalcone A, licochalcone C, and 18β-glycyrrhetinic acid were calculated using the external standard method and QAMS, and glycyrrhizic acid as the internal standard. The content of the multimarker components measured by QAMS was compared with results from ESM, using the correlation coefficient method (r) (equation (3)) and the cosine ratio value (Cir) (Equation (4)) method to verify the feasibility of QAMS [24, 25]. The angle cosine value (Cir) and the correlation coefficient (r) are commonly used to evaluate the similarity of TCM. A larger Cir value indicates a higher similarity of the target samples [26]. Meanwhile, the correlation coefficient from 0 (no correlation) to 1 (full correlation) reflects the closeness between the two variables [27]. Cir and r were calculated by using Microsoft Excel 2019 software in this study.where is the value of compound k in sample i in the ESM and is the value of compound k of the same sample in the QAMS method.

2.8. Data Analysis

The data were analyzed and assessed by the specialized software, named Similarity Evaluation System for Chromatographic Fingerprint of Traditional Chinese Medicine (TCM) (Version 2012.130723, Chinese Pharmacopoeia Committee). The similarity among different chromatograms was calculated by the correlative coefficient. HCA is a chemical pattern recognition and classification evaluation method, which is used to classify objects (samples) into clusters, such that each object is similar to the others within a cluster, but different from those in other clusters with respect to a predetermined selection criterion [28]. Using the peak areas of the eleven compounds, quality variations of the 21 batches of licorice samples were evaluated by the HCA with the software of SPSS (International Business Machines Corporation, New York, NY, USA) 22.0 to further investigate the difference among the licorice samples.

3. Results

3.1. The Optimization of the Chromatographic Separation Method

It is important to select a suitable extraction method and HPLC parameters for effective chromatograms with good separation. Extraction solvent (70% ethanol, 70% methanol, and water), extraction method (ultrasonic extraction, heated reflux extraction, and decocting method), extraction time (30, 60, and 90 min), and extraction times (once, twice, and three times) were investigated for the higher extraction efficiency. The specific steps were eventually demonstrated to be the most suitable as the section of preparation of sample preparation described.

Several mobile phase systems including methanol-water, acetonitrile-water, methanol with 0.2% glacial acetic acid, and acetonitrile with 0.2% glacial acetic acid were tested to optimize the chromatographic conditions. Flow velocity (0.8, 1, and 1.2 mL·min−1) was examined and compared. An optimized HPLC condition was determined with a thorough comparison on the resolution, base line, elution time, and number of characteristic peaks.

3.2. HPLC Method Validation
3.2.1. Linearity, LOD, and LOQ

The calibration curves of the eleven standard substances were plotted with a series of mixed standard solutions with gradient concentration. The linearity of each component was evaluated with different concentrations and the corresponding peak areas by plotting its calibration curve. Then, standard solution of the individual component was diluted gradually; the limit of quantification (LOQ) and limit of detection (LOD) values were determined by using signal-to-noise ratios of 3 : 1 and 10 : 1, respectively. As shown in Table 2, the high correlation coefficient values (R2 > 0.9990) show a good linear relationship over a wide concentration range, and the calibration curve could be used for quantitative analysis within a given concentration range.

3.2.2. Stability, Precision, Repeatability, and Accuracy

The RSD values for the precision are summarized in Table 3. From the results, the RSD results of eleven components for the instrument precision were in the range of 0.27%–2.69%, which were all lower than 3%. This indicates that this method has a good precision.

The peak areas of L1-L11 were recorded, and the values of RSD were 0.29%, 0.22%, 0.28%, 1.68%, 1.53%, 1.37%, 0.44%, 1.55%, 0.61%, 0.74%, and 2.52%, respectively. Therefore, the sample solution was stable at room temperature within 24 h (Table 3).

There was good accuracy because the average recovery rates of eleven compounds varied in the range of 86.5%–102.7%. Meanwhile, RSD values of recovery rates for each compound were lower than 3% totally.

3.3. Chromatographic Fingerprint Analysis and Similarity

The fingerprint chromatograms of 21 batches of licorice samples were recorded as shown in Figure 3(a). A chromatogram of licorice sample (Figure 3(c)) was developed with the mean of all the chromatograms by the Similarity Evaluation System for Chromatographic Fingerprint of Traditional Chinese Medicine (version 2012). Based on comparison of the chromatograms of the mixed reference standard solution (Figure 3(b)), eleven common peaks (L1–L11) of fingerprint of licorice were identified, respectively, being liquiritin, isoliquiritin apioside, isoliquiritin, neoisoliquiritin, isoliquiritigenin, formononetin, licochalcone C, licochalcone A, glabridin, glycyrrhizic acid, and 18β-glycyrrhetinic acid. The peak value of L10 was selected as the reference peak with relative retention time and the ratio of peak area as shown in Figure 3(c). As a result, their similarity values calculated based on similarity system theory of 21 samples were >0.617, which meant that the common peaks were in good correlation. Therefore, glycyrrhizic acid was chosen as the internal standard to determine the contents of liquiritin, isoliquiritin, neoisoliquiritin, isoliquiritin apioside, isoliquiritigenin, formononetin, glabridin, licochalcone A, licochalcone C, and 18β-glycyrrhetinic acid.

3.4. QAMS Method
3.4.1. RCF Calculation and Its Robustness Investigation

The RCF for component I (internal standard) is the ratio of the absolute correction factor for component I to standard substance S (reference compound). It is of great significance to select appropriate internal reference standards for the accurate determination of various components in TCM. Glycyrrhizic acid was used as an internal reference because it is inexpensive and easily obtained. And the values of RCF for the other ten compounds were computed in different concentrations according to equation (2). The average RCF of each compound is shown in Table 4.

In previous studies of QAMS, researchers have investigated the robustness of different columns, instruments, and flow rates [29]. In this study, the experiments on different instruments were carried out in different labs, using Agilent 1100, Agilent 1260, and Waters Alliance E2695 instruments. DiKMA-C18 (4.6  mm × 250  mm, 5 µm), Agilent-C18 (4.6  mm × 250  mm, 5 μm), and Thermo-C18 (4.6 mm × 250 mm, 5 µm) columns were used to study the effects of different columns. The results (Table 5) show that the RCF from different columns and instruments had good reproducibility with RSD < 3%.

3.4.2. Comparative Study of the QAMS and ESM Methods

In order to verify the availability of the QAMS method, the contents of eleven compounds in 21 batches of licorice were evaluated using the QAMS and ESM methods, respectively. The results revealed that Cir≥0.9998, r ≥ 0.9936 in Table 6, suggesting that no significant differences existed between the ESM and the QAMS method. Hence, the present QAMS method was reliable for the simultaneous quantification of the eleven components of the licorice.

3.5. Hierarchical Cluster Analysis

In order to quantitatively and objectively analyze licorice, the contents of eleven components in 21 batches of licorice samples were selected as a clustering variable to obtain a dendrogram. Euclidean distance is usually used to distinguish the similarity of samples; a larger Euclidean distance represents a greater difference and the similarity is relatively low [30]. From Figure 4, it can be seen that the 21 samples from different origins were separated into five clusters based on the dendrogram cut at a distance of 4.5. Group 1 is S5 from Ningxia in China; group 2 is S7 and S11 from Gansu in China; group 3 is S9 from Gansu in China; group 4 is S3, S4, S13, S16, S17, and S19; and group 5 contains the rest of the samples. Particularly, Ningxia has been considered as the Daodi production area (area which produces authentic and superior original medicine) of licorice in China. However, samples 1 to 6 from Ningxia showed different levels and ratios of chemical components. The possible explanation is that the licorice had been grown in different environments or in a different manner, resulting in the contents of the eleven components being different. The factors included soil type, climate, humidity, and temperature, harvesting time, planting patterns, dying methods, and other factors [31]. Therefore, the key of quality control of licorice is strictly controlling the different geographical sources to prevent adulteration.

4. Discussion

In this study, QAMS method was used to verify that the established method has high sensitivity and good precision, which is sufficient to simultaneously determine 11 compounds of 21 batches in licorice. Combined with the similarity of HPLC fingerprints, the contents of licorice from different sources show considerable differences. Although the contents of liquiritin and glycyrrhizic acid in the 21 batches of licorice were determined and accorded with the 2020 edition of the Chinese Pharmacopoeia, the licorice content of different sources is different. Furthermore, it is hard to comprehensively control the quality of licorice by using only two detection indicators. In previous studies, the contents of 4 compounds were determined by HPLC for the quality control of licorice [32]. In this study, considering the shortcomings of the multicomponents’ determination, we have combined characteristic fingerprint with QAMS to improve the quality control and evaluation of licorice. We first used QAMS to simultaneously determine the content of 11 compounds from licorice, and the feasibility of the method was verified. In addition, the contents’ variations of 11 compounds between QAMS and ESM are within the range of 5%, which indicated that there is no significant difference between the two methods. Therefore, the HPLC-QAMS method could be used for replacement of ESM in case of lack of reference substance. And it is an innovative method to provide for the comprehensive quality control of licorice.

Hierarchical cluster analysis is a chemical pattern recognition and classification evaluation method, which is used to find relatively isomorphic cluster objects based on measurement characteristics. The most similar samples are clustered preferentially [28]. In order to obtain the HCA dendrogram, the 11 components were set as variables, and Ward’s method was selected as the cluster method and squared Euclidean distance as the interval measurement [33]. The results showed that citrus and grape were divided into five main groups, and the content of licorice is significantly different in the same area. The factors included soil type, climate, humidity, and temperature, harvesting time, planting patterns, and dying methods, and other factors influenced the content of licorice. Therefore, the quality control of licorice should be taken seriously by researchers.

5. Conclusions

The quality evaluation method of the eleven components in 21 batches of licorice was established using QAMS methods, combined with HPLC fingerprints analyses and HCA method. The contents of the eleven components in licorice samples were simultaneously determined by the QAMS method. Glycyrrhizic acid was selected as the internal reference substance. In order to compare the consistency of QAMS and EMS, the angle cosine value and the correlation coefficient were calculated, and the results indicated that no significant difference existed between the newly established QAMS method and the traditional ESM.

Licorice has traditional efficacy of strengthening the spleen and benefiting vital qi and heat-clearing and detoxifying. The eleven compounds have been reported to be correlated with the traditional efficacy of licorice. So, eleven components were determined simultaneously in this experiment by QAMS-HPLC, which reflected more information than the one component (glycyrrhizic acid) determined in pharmacopoeia. Meanwhile, only one single reference substance was used in this experiment, which greatly reduced the cost of testing and was easily used for the scientific research and supervised application. On the basis of these results, the quality of licorice could be visible and better identified comprehensively by HCA and similarity analysis. Therefore, QAMS method combined with the HPLC fingerprint might provide an overall chemical characteristic, along with similarity analysis and HCA of synthesis. And this effective method has promising wide applications in the quality control of licorice and medicinal plants.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

Tingting Li, Hao Zhou, and Jiahua Ma contributed equally to this work. All the authors developed and studied the manuscript. Xueyan Fu designed this article. Tingting Li and Hao Zhou conducted the literature survey and drafted the paper. Jiahua Ma oversaw the work and revised the manuscript.

Acknowledgments

This work was supported by the Ningxia Natural Science Foundation Project (Grant no. 2018AAC03097), the 12th Five-Year National Science and Technology Support Project (Grant no. 2013BAI05B01), the National Nature Science Foundation of China (Grant nos. 81860708 and 81560647), Ningxia Key Research and Development Program (2018BEG03016), the Undergraduate Innovation and Entrepreneurship Training Program of Ningxia Medical University in 2019(19), the Ningxia Higher Education School First-Class Discipline Construction Project (Chinese Medicine) Funded Project (Grant no. NXYLXK2017A06), and School Level Project of Ningxia Medical University (XM2020022). The authors would like to thank Dr. Dong Lin and Fang Xu from Ningxia Medical University for correcting the English grammar of this article.