Quality evaluation of Polygonatum sibiricum slices from different regions based on appearance traits and multi-index metabolites combined with TOPSIS and gray relation analysis

Abstract

Introduction:

Traditional Chinese medicine quality control faces challenges, lacking multidimensionality and reliable quantitative evidence. Comprehensive evaluation models based on external characteristics and multiple indicator metabolites are the future research direction. This study focuses on Polygonatum sibiricum slices, aiming to establish a method for its quality evaluation.

Methods:

The appearance traits of P. sibiricum slices were quantified, and the contents of six functional metabolites were determined. With eight traits and six metabolite contents as variables, principal component analysis (PCA) and orthogonal partial least-squares discrimination analysis (OPLS-DA) were performed. A weighted TOPSIS-GRA fusion model was established by combining the technique for order preference by similarity to ideal solution (TOPSIS) and gray relation analysis (GRA).

Results:

The six metabolites showed good linear relationships (R² > 0.9992) within their respective ranges, with an average recovery rate of 98.54% - 103.07% (relative standard deviation less than 1.64%). Precision, stability, and repeatability met the relevant standards. There were significant differences in traits and metabolite contents among slices from different habitats. OPLS-DA identified differential quality-affecting markers. PCA showed that the first three principal components contributed over 80% of the cumulative variance, and 16 batches of slices were clustered into three categories by origin. The weighted TOPSIS-GRA fusion model indicated significant quality differences among slices from different regions, consistent with PCA and OPLS-DA clustering results.

Discussion:

The established multi-index content determination method is accurate and reliable for detecting metabolites in P. sibiricum slices. The PCA, OPLS-DA, and weighted TOPSIS-GRA fusion models are scientifically reliable. The correlation between appearance traits and product quality can be used to evaluate P. sibiricum slices from different regions, which is of great significance for quality control and standardization of traditional Chinese medicine.

1 Introduction

Traditional Chinese medicine (TCM) is the foundation for innovative developments in the medical industry in China and abroad, and forms the material basis for ensuring the clinical safety of associated treatments. The quality standards of TCM are key to the high-quality development of the industry. However, limitations in research ideas and methods mean that current research on TCM quality control lacks multidimensionality, chain continuity, and integration. TCM evaluation methods are based on the concept of “identifying patterns and discussing quality” (Kang et al., 2024), which makes it difficult to assess the internal metabolites and their efficacy, and results in a lack of reliable quantitative evidence. Modern physical and chemical evaluation methods are based on quantifying the content of a single or group of metabolites. Although accurate quantification can be achieved, the results are difficult to relate to clinical efficacy and safety. Moreover, the inherent characteristics of TCM mean that quality cannot be judged solely by metabolites. Although the external characteristics of medicinal materials are closely related to their internal metabolites and efficacy, the correlation and internal mechanisms acting between numerous characteristics and complex metabolites present several research challenges. Therefore, comprehensive evaluation models based on external characteristics and multiple indicator metabolites have become the mainstream research direction for future TCM quality control.

The technique for order preference by similarity to ideal solution (TOPSIS) assigns weights to multiple indicators, using the degree of proximity to the idealized target as a benchmark for the comprehensive evaluation of samples. This eliminates subjective factors brought about by “identifying patterns and discussing quality”, and ensures objectivity, accuracy, and a scientific basis (Xu et al., 2024). Gray relation analysis (GRA) has been used to identify complex relationships between multiple indicators, and can analyze and compare the factors that have the greatest impact on samples based on the degree of interrelation among multiple factors. In this way, GRA intuitively reflects the comprehensive value of each indicator (Chen et al., 2018).

Polygonatum sibiricum, a historically revered tonic botanical drug, is traditionally used to nourish vital energy (Qi), dispel rheumatism, and harmonize visceral functions. Modern pharmacological studies reveal its rich metabolites of polysaccharides, saponins, flavonoids, alkaloids, and phenylpropanoids, which collectively contribute to hypoglycemic, lipid-lowering, antioxidant, anti-aging, and antitumor activities (Yang et al., 2021; Zhang et al., 2024a). Despite polysaccharides being recognized as primary active metabolites in P. sibiricum (Chinese Pharmacopoeia, 2020), emerging evidence underscores the synergistic roles of secondary metabolites in its therapeutic efficacy. To address regional quality disparities and align with the “Quality Marker (Q-Marker)” concept (Yuan and Wang, 2024), this study integrates morphological traits (long diameter, short diameter, single weight, thickness, chromaticity values) with quantification of five bioactive metabolites: baicalein, liquiritigenin, neoliquiritin, 3′-methoxydaidzin, and diosgenin. The selection of baicalein, liquiritigenin, neoliquiritin, 3′-methoxydaidzin, and diosgenin was driven by a hierarchical rationale rooted in pharmacological relevance, alignment with TCM synergy principles, and methodological rigor. First, the pharmacological significance of these metabolites directly mirrors P. sibiricum’s traditional and modern therapeutic applications. For instance, baicalein, a flavonoid, was prioritized for its NF-κB-mediated anti-inflammatory and antioxidant properties, which align with the botanical drug’s historical use in rheumatism management (Gao et al., 2024; Mabalirajan et al., 2013). Moreover, liquiritigenin and its glycoside neoliquiritin were included due to their immunomodulatory and neuroprotective effects, supporting the botanical drug’s role in Qi replenishment and immune enhancement (Yao et al., 2022; Xu et al., 2024). Furthermore, 3′-methoxydaidzin, an isoflavone with estrogenic activity, reflects P. sibiricum’s yin-nourishing applications in alleviating menopausal symptoms (Hou et al., 2024), while diosgenin, a steroidal saponin, was selected for its antitumor and lipid-lowering mechanisms via apoptosis induction and metabolic regulation (Xu et al., 2022; Zhang et al., 2021). Critically, these metabolites collectively embody the TCM paradigm of multi-metabolite synergy, wherein therapeutic efficacy arises from interconnected interactions rather than isolated constituents (Xiao et al., 2025). Finally, methodological coherence was ensured by integrating prior Q-Marker predictions (Jiang et al., 2017; Wang et al., 2024) with TOPSIS-GRA frameworks, enabling robust correlations between morphological traits and metabolite levels to harmonize traditional and modern quality standards.

By synthesizing morphological and multi-metabolite data, this work establishes a scientifically robust framework for evaluating P. sibiricum quality, offering insights for authenticity verification, origin tracing, and standardized industrial practices. The findings advance TCM quality standardization and underscore the necessity of integrating traditional wisdom with contemporary analytical methodologies.

2 Materials and methods

2.1 Test drugs

The reference substances of baicalein (batch No. A1119L021, purity 98.5%), liquiritigenin (batch No. A105J021, purity 98.5%), neoliquiritin (batch No. 2230816001, purity 98.5%), 3′-methoxydaidzin (batch No. Y14F10w79687, purity 98.5%), and diosgenin (batch No. A922E024, purity 98.5%) were procured from Beijing Solarbio Science & Technology Co., Ltd. High-performance liquid chromatography-grade acetonitrile, methanol, ethanol, and glacial acetic acid were supplied by Tianjin Fuyu Fine Chemical Co., Ltd. The P. sibiricum slices were authenticated as the dried rhizomes of Polygonatum by Professor Cairang Nanjia from the School of Pharmacy at Qinghai Minzu University. The P. sibiricum specimens used in this study are preserved in the Specimen Laboratory of the School of Pharmacy, Qinghai Minzu University. Table 1 details the information regarding the sample collection sites. As shown in Figure 1, the slices of P. sibiricum are presented, which can provide a visual reference for the morphological characteristics of this TCM.

FIGURE 1

2.2 Instruments

To prepare, examine, and characterize the samples, we used a Shimadzu LC-20A High Performance Liquid Chromatograph (Shimadzu, Japan); Analytical High-Performance Liquid Chromatography Column (Hypersil ODS-C18: Elite Suzhou Analytical Instruments Co., Ltd., China); KQ-50DA CNC Ultrasonic Cleaner (Kunshan Hechuang Ultrasonic Instrument Co., Ltd., China); XS105DU Electronic Balance (Shanghai Mettler Toledo Co., Ltd., China); SHE-D Vacuum Pump (Gongyi Yuhua Instrument Co., Ltd., China); UPH-IV Ultra Pure Water Machine (Chengdu Ultra Pure Science & Technology Co., Ltd., China); and 3nh Color Difference Meter (Guangdong Sannshi Intelligent Technology Co., Ltd., China).

2.3 Appearance characteristic measurement

Slices of P. sibiricum from diverse origins were thoroughly homogenized and arranged to form a square configuration. Subsequently, a partition board was employed to divide the square into four congruent sections along the diagonal. Two diagonally opposed sections were selectively retained, while the remaining two sections were discarded. The two selected portions were then re-combined, and this procedure was iteratively repeated until the number of samples per batch attained 20. A precision balance and a vernier caliper were utilized to measure the mass, thickness, major diameter, and minor diameter of the samples. Subsequently, the average values of these measured parameters were calculated.

2.4 Decoction pieces chromaticity measurement

A colorimeter was employed to conduct measurements under the D65 light source, with a measurement wavelength range spanning from 180 to 740 nm and a measurement field of 10°. Prior to sample measurement, black and white calibration was performed on the colorimeter. Each batch, consisting of 20 samples, was placed within the test chamber, ensuring that the colorimeter’s measuring spot was precisely aligned at the center of the test chamber. The colorimetric values L^*, a^*, and b^* were then recorded. Here, L^* denotes the brightness of the sample color, a^* represents the red - green gradient axis, and b^* represents the yellow - blue gradient axis, where larger absolute values correspond to more intense colors (Yang et al., 2020). Each measurement was replicated three times, and the average values were calculated. Subsequently, the comprehensive color difference value ∆E for each batch of samples was computed according to Equation 1.

2.5 Determination of indicator metabolite contents

2.5.1 Preparation of standard solution

The requisite amounts of the reference substances, namely, baicalein, liquiritigenin, neoliquiritin, diosgenin, and 3′-methoxydaidzin, were accurately weighed. These substances were then dissolved in chromatographic-grade methanol and transferred to a 25 mL volumetric flask. The volume was adjusted to the mark with the same solvent, thereby preparing a mixed reference substance solution. The resulting mass concentrations of baicalein, liquiritigenin, neoliquiritin, diosgenin, and 3′-methoxydaidzin in the solution were 0.0816, 0.094, 0.0828, 0.098, and 0.0796 mg·mL⁻¹, respectively.

2.5.2 Preparation of test Sample solution

Appropriate quantities of P. sibiricum samples were meticulously pulverized and sieved. Subsequently, precisely 2.00 g of the resultant powder was weighed and transferred into a conical flask. To this, 40 mL of 70% ethanol was added, and the mixture was subjected to sonication for a duration of 30 min. After sonication, the mixture was cooled and then filtered. A portion of the filtrate was collected, and the above-mentioned procedures were reiterated for the filter residue. The filtrates were then pooled together and concentrated under reduced pressure to yield a paste. This paste was re-dissolved in methanol to a final volume of 10 mL. Subsequently, 100 µL of a 10 g·L⁻¹ chitosan-glacial acetic acid solution was added dropwise, and the resulting solution was left to stand overnight. Finally, the solution was filtered through a 0.22 - µm filter membrane to obtain the test solution.

2.5.3 High-performance liquid chromatography conditions

High-performance liquid chromatography (HPLC) analysis was performed using a Hypersil ODS-C18 column (250 mm × 4.6 mm, 5 µm particle size). The detection wavelength was set at 230 nm. The mobile phase, consisting of acetonitrile (A) and water (B), was applied for gradient elution with the following program: 0–20 min, 5%–30% A; 20–30 min, 30%–40% A; 30–40 min, 40%–60% A; 40–80 min, 60%–100% A. The injection volume was 10 μL, the column temperature was maintained at 30°C, and the flow rate was 1.0 mL·min⁻¹. Chromatograms of the reference and test samples of P. sibiricum slices under these conditions are shown in Figure 2.

FIGURE 2

2.5.4 Investigation of the linear correlation

In accordance with the chromatographic conditions described in Section 2.5.3, injection volumes of 4, 8, 12, 16, 20, and 24 µL were prepared. The corresponding chromatographic peak areas were recorded. Linear regression equations were established, with the injection volume as the independent variable (X, abscissa) and the peak area as the dependent variable (Y, ordinate).

2.5.5 Precision test

An appropriate volume of the mixed standard solution was subjected to six consecutive injections (10 μL per injection) under the chromatographic conditions specified in Section 2.5.3. The chromatographic profiles were recorded and the relative standard deviation (RSD) values of the peak areas for five target analytes - baicalein, liquiritigenin, neoliquiritin, diosgenin, and 3′-methoxydaidzein - were systematically calculated to evaluate the precision of the analytical method.

2.5.6 Stability test

A sample (S1) solution of P. sibiricum slices was extracted and detection was performed every 4 h over a 24-h period. The baicalein, liquiritigenin, neoliquiritin, diosgenin, and 3′-methoxydaidzin peak areas were recorded.

2.5.7 Repeatability test

Six samples of P. sibiricum powder (S1) were prepared in parallel according to the method described in Section 2.5.2. The samples were injected according to the chromatographic conditions detailed in Section 2.5.3, and the peak area was determined.

2.5.8 Sample recovery test

Nine aliquots of S1 powder containing known amounts of the five target analytes (baicalein, liquiritigenin, neoliquiritin, diosgenin, and 3′-methoxydaidzin) were prepared. Three aliquots each of low-, medium-, and high-mass-fraction control solutions (three replicates per concentration level) were then spiked. The test solution was prepared according to the method in Section 2.5.2, and chromatographic analysis was performed under the conditions described in Section 2.5.3. Finally, the average recoveries and RSD were calculated by comparing the measured amounts with the spiked amounts.

2.5.9 Content determination

Sixteen batches of samples were pulverized and sieved, and then test solutions were prepared according to the method detailed in Section 2.5.2, after which the prepared test solutions were diluted with methanol at a ratio of 1:10, vortex-mixed for 30 s, and allowed to stand for 5 min.The prepared samples were injected under the chromatographic conditions specified in Section 2.5.3, and the corresponding peak areas were recorded. The contents of baicalein, liquiritigenin, neoliquiritin, diosgenin, and 3′-methoxydaidzin were calculated by applying the pre-established linear regression equations.

2.6 Determination of polysaccharide content in P. sibiricum slices

2.6.1 Drawing of standard curve

The anhydrous glucose was dried to constant weight at 105°C, and precisely 33 mg was taken and added to a 100-mL volumetric flask. Then, 0, 0.2, 0.5, 1, 1.5, and 2 mL of glucose solution was added to six 10-mL test tubes, with distilled water added up to 2 mL where necessary, followed by the addition of 0.2% anthrone–sulphuric acid solution in an ice-water bath up to 10 mL. The test tubes were left at room temperature before being subjected to a boiled water bath for 10 min, followed by a room temperature water bath for 10 min. The absorbance was then measured at 620 nm. The standard curve was plotted, with the concentration c of glucose standard solution as the horizontal coordinate and the absorbance A as the vertical coordinate. The results showed that there was a good linear relationship between absorbance and concentration, with the regression equation A = 34.661c+0.8079 achieving R² = 0.9998.

2.6.2 Determination of polysaccharide content

The sample powder was dried at 60°C to constant weight. Subsequently, 0.25 g of the dried powder was accurately weighed and transferred into a round-bottomed flask. Then, 150 mL of 80% ethanol was added. Ultrasonic extraction was performed for 15 min, followed by reflux extraction in a boiling water bath for 1 h. The hot solution was filtered, and the residue was washed three times with 10 mL of 80% ethanol each. The residue and filter paper were transferred to another flask, and 150 mL of distilled water was added. Reflux extraction in a boiling water bath was repeated for 1 h, followed by hot filtration. The filtrate and washings were combined and further washed four times with 10 mL of distilled water each (filtered while hot). The residue and flask were also washed four times with 10 mL of hot distilled water. All filtrates and wash solutions were combined, cooled to room temperature, and transferred quantitatively to a 250 mL volumetric flask. The volume was adjusted to the mark with distilled water, and the flask was shaken thoroughly to obtain the test solution.

Exactly 1 mL of the test solution was pipetted into a 10 mL stoppered test tube, and 2 mL of distilled water was added. After thorough mixing, 0.2% anthrone-sulfuric acid solution was slowly added to the marked volume in an ice-water bath. The mixture was shaken gently, cooled, and then placed in a boiling water bath for 10 min. After cooling in a room-temperature water bath for 10 min, the solution was equilibrated at ambient conditions. Using distilled water as the blank, absorbance was measured at 620 nm. The absorbance value was substituted into the pre-established linear regression equation to calculate the polysaccharide content.

2.7 Statistical analysis

In the SIMCA 14.1 software, the orthogonal partial least-squares discrimination analysis (OPLS-DA) function was used to analyze the variable importance in the projection (VIP) of the indicator metabolites and polysaccharides. Duncan’s multiple range test (P < 0.05) was applied to determine the level of statistical significance after principal component analysis (PCA) and one-way analysis of variance in IBM’s SPSS Statistics 27.0 software. All data (mean ± standard deviation) were determined independently in triplicate.

2.8 Weighted TOPSIS-GRA fusion modeling

2.8.1 Weighted TOPSIS data normalization

The analysis described in the previous subsections suggests that larger values of the mass, thickness, long diameter, short diameter, baicalein, liquiritigenin, neoliquiritin, diosgenin, 3′-methoxydaidzin, and polysaccharides of P. sibiricum slices are preferable. These variables were normalized using Equation 2. In contrast, smaller values of L^*, a^*, b^*, and ∆E are desirable. These variables were normalized using Equation 3 (Liao et al., 2024).

Note: is the normalization result and is the value of each indicator measurement.

2.8.2 Weighted TOPSIS Calculation for each evaluation indicator

The variable importance in projection (VIP) values of the 14 indicators derived from orthogonal partial least-squares discriminant analysis (OPLS-DA) were assigned as the weighting coefficients (Q_j). The weighting matrix values (Z_ij) for each indicator were then computed using Equation 4. Based on this matrix, the positive and negative ideal samples (Z⁺) and (Z⁻) were defined. Finally, the distances from each evaluation indicator to the ideal samples (d_i⁺) and (d_i⁻) were calculated via Equations 5, 6 (Zheng et al., 2024).

2.8.3 GRA standardization

Because of the significant differences in the results of the 14 indicators measured in the 16 batches of P. sibiricum slices, data standardization was required prior to GRA (Li et al., 2024). The data were standardized using Equation 7.

Note: is the result of the standardization process for each indicator and is the average of the results for each indicator.

2.8.4 Calculation of GRA correlation coefficients

The GRA method evaluates the interrelationships among research objects by quantifying their degrees of relevance. Higher GRA correlation values indicate stronger associations and superior quality of the objects (Fu et al., 2019). Among the 16 batches of P. sibiricum slices, parameters including mass, thickness, long diameter, short diameter, and the contents of baicalein, liquiritigenin, neoliquiritin, diosgenin, 3′-methoxydaidzin, and polysaccharides were classified as positive indicators. For these positive indicators, the maximum values observed across all batches were selected as the optimal reference sequence, while the minimum values were designated as the worst reference sequence during GRA standardization. Conversely, the color parameters L^*, a^*, b^*, and ∆E were considered negative indicators. For these, the maximum values from GRA standardization were assigned as the worst reference sequence, and the minimum values as the optimal reference sequence. The correlation coefficients between each evaluation index and both reference sequences (optimal and worst) were calculated using Equations 8, 9.

Note: are the correlation coefficients of the best and worst reference sequences; are the best and worst reference sequences; , , , ; α is the resolution coefficient, which takes a value of 0.5 (Zhang et al., 2020).

2.8.5 Calculation of relative correlation degree of GRA

Using the process explained in Section 2.8.4, the relative correlation of each batch was calculated. The relative correlation is the mean value of the correlation coefficients of all evaluation indicators of the research object. Thus, the correlation of P. sibiricum slices with respect to the optimal and the worst reference sequences was calculated using Equations 10, 11 (Qu et al., 2019), respectively.

2.8.6 Weighted TOPSIS-GRA fusion Modeling

Using Equations 12–15, the Euclidean distance and relative correlation of the reference sequence for 16 batches of P. sibiricum slices were subjected to dimensionless processing. Subsequently, post-processing by the weighted TOPSIS-GRA fusion model was applied using Equations 16, 17. D_i⁺ and Q_i⁻ denote the extent to which samples deviate from the ideal sample; smaller values indicate closer proximity to the ideal sample. Conversely, D_i⁻ and Q_i⁺ signify the extent to which samples approach the ideal sample; larger values suggest a closer resemblance to the ideal sample. G_i⁺ offers a comprehensive assessment of how closely all samples approximate the ideal sample; higher values equate to superior sample quality. G_i⁻ provides a comprehensive evaluation of the distance of all samples from the ideal sample; higher values correspond to lower sample quality. The relative closeness κ was determined using Equation 18; higher values signify higher comprehensive sample quality.

Note: α and β are correlation coefficients, both taking a value of 0.5 (Yan et al., 2024).

3 Results

3.1 Appearance characteristics of P. sibiricum slices

The measurements of appearance-characteristic indices are presented in Table 2. As shown in the table, the 16 batches of samples exhibit significant differences among different origins based on eight evaluated traits.

3.2 Methodological review

3.2.1 Linear range Inspection

The standard curve of five quality metabolites was plotted, and the linear regression equation, correlation coefficient, and linear range of each reference substance was calculated. The results are presented in Table 3.

3.2.2 Precision test results

The RSD values for baicalein, liquiritigenin, neoliquiritin, diosgenin, and 3′-methoxydaidzin were determined to be 1.75%, 0.65%, 0.31%, 1.56%, and 1.06%, respectively. These findings indicate that the precision of the analytical instrument meets the requirements of the test.

3.2.3 Stability test results

The RSD values of baicalein, liquiritigenin, neoliquiritin, diosgenin, and 3′-methoxydaidzin were found to be 1.78%, 1.71%, 1.96%, 1.10%, and 1.98%, respectively. These results show that the solutions of P. sibiricum slices attain good stability within 24 h.

3.2.4 Repeatability test results

The RSD values of baicalein, liquiritigenin, neoliquiritin, diosgenin, and 3′-methoxydaidzin were found to be 1.90%, 1.86%, 1.96%, 1.66%, and 0.92%, respectively. These results indicate that the method has good repeatability.

3.2.5 Sample recovery test results

The baicalein, liquiritigenin, neoliquiritin, diosgenin, and 3′-methoxydaidzin average recoveries were 99.17%, 100.54%, 98.54%, 103.07%, and 99.08%; the RSD values were 1.51%, 0.87%, 1.64%, 1.09%, and 1.62%, respectively. These results demonstrate that the method satisfies the requirements of the spiked sample recovery test, indicating a high degree of accuracy.

3.3 Determination of indicator metabolites and polysaccharide content

The contents of the five indicator metabolites in P. sibiricum slices were quantified using the method detailed in Section 2.5. The polysaccharide content was assayed by the technique described in Section 2.6. The resultant data are presented in Table 4.

The results in Table 4 indicate that P. sibiricum slices of different origins exhibit significant differences in all indicator metabolites, except for diosgenin. The highest contents of baicalein, liquiritigenin, neoliquiritin, 3′-methoxydaidzin, and polysaccharides originated from samples S16, S5, S16, S6, and S8, respectively.

3.4 PCA and OPLS-DA

3.4.1 PCA

PCA was applied to the 16 batches of P. sibiricum slices, with the appearance index data (mass, thickness, long diameter, short diameter, and chromaticity values L^*, a^*, b^*, ∆E) and the contents of baicalein, liquiritigenin, neoliquiritin, diosgenin, 3′-methoxydaidzin, and polysaccharides as the variables. The analysis was performed using SPSS27.0 software. The results are presented in Tables 5, 6.

As shown in Table 5, the first three principal components exhibit eigenvalues greater than 1 (8.011, 2.324, and 1.100) (Zhang et al., 2024b), contributing variances of 52.219%, 16.601%, and 7.858%, respectively. Collectively, these components account for 81.678% of the total variation in P. sibiricum slices. Table 6 demonstrates that mass, thickness, long diameter, short diameter, L^*, b^*, ∆E, baicalein, liquiritigenin, neoliquiritin, 3′-methoxydaidzin, and polysaccharides exhibit higher loadings on the first principal component, while a^* and diosgenin are more strongly associated with the second principal component. A PCA model was constructed using SIMCA 14.1 software to explore the clustering patterns of P. sibiricum slices. The results are presented in Figure 3.

FIGURE 3

As shown in Figure 3, the 16 batches of P. sibiricum slices can be divided into three categories with significant clustering and dispersion, indicating that there are large differences in quality between the categories. Samples S4 and S11–S13 are clustered into one category; S2, S8, S14, and S16 are clustered into another category; and the remaining samples are clustered into a third category.

3.4.2 OPLS-DA

To further investigate the factors influencing the quality differences between categories, the data matrix of 14 indicators in the 16 batches of P. sibiricum slices was imported into SIMCA 14.1 and OPLS-DA was applied. The results are presented in Figures 4, 5. The model parameters R²_X = 0.986, R²_Y = 0.971, and Q² = 0.853 indicate that the model has good stability, reliability, and predictive ability.

FIGURE 4

FIGURE 5

Figure 4 reveals that the 16 batches of P. sibiricum slices are categorized into three distinct groups, consistent with the PCA results. Figure 5 presents the VIP values. Applying a VIP threshold >1 as the screening criterion (Gan et al., 2021), polysaccharides, long diameter, baicalein, b^*, and 3′-methoxydaidzin were identified as key markers for quality classification of P. sibiricum slices. These parameters exhibit substantial contributions to the overall quality evaluation and can effectively discern quality variations. The OPLS-DA model was validated through 200 permutation tests, with results displayed in Figure 6.

FIGURE 6

Figure 6 shows that the resultant R² fitted curve intercepts the Y-axis below Y = 0.3 and the Q² fitted curve intercepts the Y-axis below Y = 0. The results of cross validation analysis of variance (CV-ANOVA) indicated that the F value was 135.76 (P < 0.01), suggesting that the model was significantly effective. Therefore, the OPLS-DA model can be applied to evaluate the quality of P. sibiricum slices from different regions.

3.5 Matrix calculation of weighted TOPSIS

Table 7 presents the TOPSIS weighting matrix for each index, and Table 8 lists the distances between each evaluation index and the positive/negative ideal samples (d_i⁺, d_i⁻).

The data presented in Tables 7,8 demonstrate that the distances between the 16 batches of samples and the positive ideal samples span from 0.5138 to 3.4690, while the distances to the negative ideal samples range from 0.6597 to 3.4573. Samples S2, S8, and S14 exhibit relatively higher quality, whereas S4, S11, and S13 show markedly lower quality. These findings indicate that P. sibiricum slices sourced from diverse geographical origins exhibit a heterogeneous quality profile, with statistically significant disparities.

3.6 Relative correlation degree of GRA

As shown in Table 9, the correlation coefficients of the 16 batches of P. sibiricum slices range from 0.3895 to 0.8094 for the optimal reference sequence and 0.3811 to 0.8085 for the worst reference sequence. Notably, samples S2, S8, S14, and S16 exhibit significantly stronger correlations with the optimal sequence (coefficients >0.7000), indicating superior quality. In contrast, samples S4, S11, S12, and S13 demonstrate higher correlations with the worst sequence (coefficients >0.7000), suggesting markedly inferior quality. These findings are consistent with the results of the weighted TOPSIS model.

3.7 Calculation of relative closeness κ

The results of the weighted TOPSIS-GRA fusion model demonstrate that the average relative closeness values of the 16 batches of P. sibiricum slices range from 0.2526 to 0.7635 (Table 10), indicating substantial quality variation across different geographical origins. The top four samples (S14, S8, S16, and S2) exhibit relative closeness values exceeding 0.7000, while the remaining samples fall below 0.6000. Notably, P. sibiricum slices from Zunyi (Guizhou Province), Yibin (Sichuan Province), Jinggu (Yunnan Province), and Shizhu (Chongqing Municipality) show significantly higher quality compared to other regions within the same provinces. These findings align with historical records of P. sibiricum’s traditional Dao-di producing areas (Ma et al., 2019), thereby validating the reliability of the TOPSIS-GRA fusion model.

4 Discussion

4.1 Relationship between morphological traits and quality attributes in P. sibiricum slices

Morphological traits are critical yet undervalued in P. sibiricum quality assessment. Single-marker methods are insufficient; integrating morphological and metabolic perspectives is essential (Chinese Pharmacopoeia, 2020; Yuan and Wang, 2024). Processed P. sibiricum quality links to traits like slice long diameter and chromatic parameter b^*, reflecting regional environments and bioactive metabolite accumulation (Liang et al., 2025; Shao et al., 2025). Longer diameters indicate better-developed tissues (potentially higher active ingredients), while b^* correlates with specific metabolites, signaling environmental impacts on biosynthesis. These traits are indispensable for quality evaluation, offering insights into metabolic processes and therapeutic efficacy.

4.2 Regional differences in polysaccharide and metabolite profiles of P. sibiricum slices

Geographical differences significantly shape P. sibiricum polysaccharide and metabolite profiles. OPLS-DA revealed distinct patterns: Sichuan and Guizhou samples exhibited elevated polysaccharides and region-specific metabolites (e.g., flavonoids, saponins) (Chen et al., 2024). Polysaccharides, a key pharmacopeial metabolite, underpin pharmacological effects, while metabolites like 3′-methoxydaidzin and baicalein—regionally variable—exhibit anti-inflammatory activities. Aligned with the “Dao-di botanical drugs” concept, these differences facilitate standardized sourcing and quality control in pharmaceutical applications.

4.3 Multi-model integrated approach for quality evaluation of P. sibiricum slices

This study introduces a TOPSIS-GRA methodology weighted by OPLS-DA-derived VIP values, providing a robust framework for P. sibiricum quality assessment. By incorporating VIP values—quantifying variable importance in sample differentiation—the model mitigates subjectivity, enabling objective, comprehensive evaluations. Industrially, it supports quality prioritization, raw material selection, and origin traceability—critical for pharmaceutical compliance. Consistent with empirical evidence, this approach ensures safety and quality in P. sibiricum-based products by integrating multi-dimensional factors.

In summary, this research establishes a morphology-metabolism linkage paradigm, advancing TCM quality control. Future studies should dissect biosynthetic pathways underlying regional traits to optimize cultivation, ensuring clinical consistency and enhancing TCM product quality.

5 Conclusion

We demonstrated that high-quality P. sibiricum slices are typically distinguished by their elongated, yellowish-white appearance, and contain elevated levels of polysaccharides, baicalein, and 3′-methoxydaidzin. Manufacturers should concentrate on managing the five quality-determining factors identified above, while sellers and consumers can perform preliminary selection based on the length and color attributes of the slices. Our group intends to further investigate the correlation between the regional characteristics of P. sibiricum slices and their appearance, chemical composition, and pharmacological activities. This will lay the groundwork for the development of a quality control system for botanical drugs that embody the three excellences of shape, quality, and superior effect.

References

• 1

  ChenJ. L.TanM. X.ZouL. S.LiuX. H.ChenS. Y.ShiJ. J.et al (2018). Stimultaneous determination of multiple bioactive constituents in Panacis Japonici Rhizoma processed by different methods and grey relational analysis. China J. Chin. Mater. Med.43, 4274–4282. 10.19540/j.cnki.cjcmm.20180726.001

  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  ChenZ. L.WangJ. K.QinC. Y.LiW.DuH. Z. (2024). Preliminary study on the grade quality of Polygonatum Cyrlonema. Chin. J. Ethnomed. Ethnopharm.33, 33–40. 10.3969/j.issn.1007-8517.2024.02.zgmzmjyyzz202402008

  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  FuS. P.ShenH. W.WangQ. B.LiJ. P.WangC.GuoS. L.et al (2019). Study on fingerprint and gray correlation analysis of Paeoniae Radix Rubra from different habitats. Chin. Tradit. Herb. Drugs50, 5865–5871. 10.7501/j.issn.0253-2670.2019.23.030

  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  GanL.YanY. J.DingQ.HongW. M.JiY. H.YangX. D.et al (2021). Research on quality difference of standard decoction of raw and fried Paeoniae Radix Alba based on fingerprint and multicomponent determination. China J. Chin. Mater. Med.46, 1410–1416. 10.19540/j.cnki.cjcmm.20201130.301

  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  GaoJ. T.LuC. Y.ShiX. C.HeX. (2024). Progress on pharmacological effects and solubilization technology of baicalein. Prog. Vet. Med.45, 83–87. 10.3969/j.issn.1007-5038.2024.11.015

  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  HouT. H.ZhaoY. M.WangF.ZhangB. W.ZhangY. H.LiuX. F. (2024). The physiological and Pharmacological activities of isoflavones. Shandong Chem. Ind.53, 70–71. 10.3969/j.issn.1008-021X.2024.21.022

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  JiangC. X.ZhangT. J.ChenC. Q.LiX. K.LiuX. C. (2017). Research progress in Polygonati Rhizoma and predictive analysis on Q-marker. Chin. Tradit. Herb. Drugs48, 1–16. 10.7501/j.issn.0253-2670.2017.01.001

  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  KangT. G.ZhangH.LiF.WangT. M.XuL.ZhangJ. K.et al (2024). Discussion on Chinese materia medica characterology. Chin. Tradit. Herb. Drugs55, 3579–3587. 10.7501/j.issn.0253-2670.2024.11.001

  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  LiH. Y.WangH. R.NaL. S.JinC. H.WangJ. L. (2024). Quality evaluation of Panacis Majoris Rhizoma from different producing areas by PCA, OPLS-DA and weighted TOPSIS-GRA fusion model. Chin. Tradit. Herb. Drugs55, 3116–3126. 10.7501/j.issn.0253-2670.2024.09.025

  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  LiangM. F.LiaoN.ZhuS. S.WangZ. J.WanX. F. (2025). Fusion model of multi-index component quantitative combined with Chemometrics,Weighted TOPSIS and grey relational analysis for quality evaluation of Euphorbia hirta from different regions. Tradit. Chin. Drug Res. Clin. Pharmacol.36, 125–133. 10.19378/j.issn.1003-9783

  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  LiaoC. L.WuY. H.WuY. H.WangZ. G. (2024). Application of QAMS, chemometrics and weighted TOPSIS in multi-index quality control and difference evaluation of Humulus scandens from different producing areas. Chin. Tradit. Herb. Drugs55, 4907–4916. 10.7501/j.issn.0253-2670.2024.14.026

  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  MaY. L.ZhangY. S.ZhangF.ZhangY. Y.ThakurK.ZhangJ. G.et al (2019). Methyl protodioscin from Polygonatum sibiricum inhibits cervical cancer through cell cycle arrest and apoptosis induction. Food Chem. Toxicol.132, 1110655. 10.1016/j.fct.2019.110655

  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  MabalirajanU.AhmadT.RehmanR.LeishangthemG. D.DindaA. K.AgrawalA.et al (2013). Baicalein reduces airway injury in allergen and IL-13 induced airway inflammation. PLos One.8, e62916. 10.1371/journal.pone.0062916

  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  QuL. H.CaoG. S.TuJ. Y.ShiK.WangJ.XuY. Y.et al (2019). Spectrum-effect relationship of dryness of Atractylodis Rhizoma volatile oil based on grey relational grade and orthogonal partial least squares analysis. Chin. Tradit. Herb. Drugs50, 150–156. 10.7501/j.issn.0253-2670.2019.01.023

  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  ShaoL.QiaoL.KangB. W.LiuD. J.ZhuL. C. (2025). Quality evaluation of Sambuci Williamsii Ramulus from different origins based on chemometrics combined with weighted TOPSIS and GRA fusion model. Nat. Prod. Res. Dev.37, 438–449. 10.16333/j.1001-6880.2025.3.006

  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  WangC.YeJ.HeX. G.TangY. Y. (2024). Research progress on chemical components and pharmacological effects of Polygonati Rhizoma and prediction analysis of quality marker. Nat. Prod. Res. Dev.36, 881–899+855. 10.16333/j.1001-6880.2024.5.017

  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  XiaoW. T.ZengQ. L.ZhangH.HanY.ZhouH. R.XiaoX.et al (2025). Research progress on new drugs derived from effective components of single herb of traditional Chinese medicine. Chin. New Drugs J.34, 264–269. 10.20251/j.cnki.1003-3734.2025.03.007

  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  XuJ. H.ZhouH.LiZ. M.LuoC. (2024). Quality evaluation of Pteris multifida from different producing areas based on QAMS combined with chemical pattern recognition method and weighted TOPSIS and grey correlation degree fusion model. Chin. Tradit. Herb. Drugs55, 6738–6748. 10.7501/j.issn.0253-2670.2024.19.026

  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  XuS. L.MaZ.XingL. H.ChengW. Q. (2024). Polygonatum sibiricum component liquiritigenin restrains breast cancer cell invasion and migration by inhibiting HSP90 and chaperone-mediated autophagy. Korean J. Physiol. Pharmacol.4, 379–387. 10.4196/kjpp.2024.28.4.379

  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  XuY. S.GaoM. J.ZhaoN. D.HuangZ. Q.TianY. Q.HuZ. Y.et al (2022). Establishment of quality standard for freeze-dried tablets of Polygonatum sibiricum Delar. ex Redoute and study on the anti-tumor activity of diosgenin. Hubei Agric. Sci.61, 133–137+141. 10.14088/j.cnki.issn0439-8114.2022.13.025

  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  YanQ. Q.WangJ. C.ZhangW. Y. (2024). Evaluation of quality of Valeriana officinalis from different origins based on multi indicator component quantitative combined chemometrics, weighted TOPSIS, and grey relational degree fusion model. Drug Eval. Res.47, 1305–1315. 10.7501/j.issn.1674-6376.2024.06.015

  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  YangB. F.XuF.LiS. L.LiL. F.LiL. B.LiZ. F.et al (2021). Research advances on chemical composition,physiological function and industrial development of Polygonatum sibiricum. J. Anhui Agr. Sci.49, 8–12. 10.3969/j.issn.0517-6611.2021.11.003

  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  YangL.GongY. T.XuM. S.YangL.ChenJ. B.DongL. (2020). Research on relationship between color and composition changes in processing of rhubarb charcoal based on “external and internal correlation”. Chin. Tradit. Herb. Drugs51, 5705–5713. 10.7501/j.issn.0253-2670.2020.22.006

  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  YaoD.ShiB.WangS. C.BaoL. X.TanM.ShenH.et al (2022). Isoliquiritigenin ameliorates ischemia-induced myocardial injury via modulating the Nrf2/HO-1 pathway in mice. Drug Des. devel. Ther.16, 1273–1287. 10.2147/DDDT.S362754

  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  YuanY. J.WangQ. H. (2024). Research progress on the chemical components and pharmacological effects of Polygonati rhizoma and prediction analysis of its quality markers. Inf. Tradit. Chin. Med.41, 72–80+86. 10.19656/j.cnki.1002-2406.20240213

  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  ZhangN.WengW. F.WeiK. S.HuangC.WangQ. (2020). Quality evaluation of Curcumae Radix from different areas based on grey correlation degree and TOPSIS fusion model. Chin. J. Exp. Tradit. Mde. Form.26, 137–145. 10.13422/j.cnki.syfjx.20200215

  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  ZhangS. Q.SongY. X.ZhangW. X.ChenM. J.ManS. L. (2021). Research on anti-tumor natural product diosgenin. China J. Chin. Mater. Med.46, 4360–4366. 10.19540/j.cnki.cjcmm.20210610.701

  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  ZhangX. Y.HuangZ. Y.HeJ. W.JiangH. M.LiuY.DiJ. X. (2024a). Quality evaluation of Sanguisorba officinalis from different origins by fingerprint combined with multi-component quantification and chemical pattern recognition. Chin. Tradit. Herb. Drugs55, 5649–5655. 10.7501/j.issn.0253-2670.2024.16.026

  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  ZhangY. H.WuY. W.XuY. Y.ZhaoD.XuT. H. (2024b). Pharmacological effects and patent development progress of Polygonatum sibiricum. Mod. Chin. Med.26, 181–188. 10.13313/j.issn.1673-4890.20230407002

  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  ZhengW. L.ZhangQ.FengF.WuW. F.ZhangW.XieH.et al (2024). Optimization of extraction process for classical prescription Qianghuo Shengshi Decoction based on standard relation and entropy weight-TOPSIS method. Chin. Tradit. Herb. Drugs55, 3310–3320. 10.7501/j.issn.0253-2670.2024.10.010

  • CrossRef
  • Google Scholar
  • View reference in article