Oligosaccharide-based quality evaluation of Atractylodis rhizome and a strategy for simplifying its quality control

Background Atractylodis rhizoma, is the dried rhizomes of Atractylodes lancea (Thunb.) DC. or A. chinensis (DC.) Koidz. Both of two are pharmacologically and economically important, while with differences in efficacy. Therefore, an authentication system is vital for evaluation the quality and discrimination adulteration of Atractylodis rhizoma. Fructooligosaccharides (FOS), which are regarded as functional ingredients in Atractylodis rhizoma, have not been used for quality control of Atractylodis rhizoma for shortage of reference compounds. Results A HPLC-ELSD method was developed for the quantification of FOS in Atractylodis rhizoma. And chemometrics analysis showed that 2 markers including content of degree of polymerization (DP) 12 and total content of DP 3-15 could be used as the main distinctive elements for quality evaluation of Atractylodis rhizome. Actually, the separation and purification of high DP FOS, such as DP 12, is still a challenge because of high polarity. Then DP 5-based qualification evaluation was investigated for quality control of Atractylodis rhizoma. The results showed that A. lancea and A. chinensis could be clearly separated. Conclusions DP 5-based quantification method was credible and effectively adopted for solving the shortage of reference compounds and improving the quality control of Atractylodis rhizoma


Introduction
Herbal medicine has gained popularity in many countries throughout the years. With the increased trend of usage worldwide, the assessment of safety, quality, and efficacy of these medicines has been an important concern for the public [1]. Atractylodis rhizome is a traditional perennial herb, which has been known and used as medication for a very long time. It is widely distributed in eastern Asia [2]. Korean and Japanese pharmacopoeias have recorded this herb as traditional diuretic and gastric prescriptions [3,4]. Atractylodis rhizoma origin from dried rhizomes of A. lancea (Thunb.) DC. and A. Chinensis (DC.) Koidz. documented in Chinese Pharmacopoeia 2015 [5]. It is mainly distributed in Heilongjiang, Hebei, Jiangsu, Zhejiang, Anhui and Hubei province. It has been firstly reported in Shen Nong Ben Cao Jing, the earliest masterpiece about Chinese medicine [6], reputed to invigorate the spleen and eliminate dampness. And it is commonly used for the treatment of rheumatic diseases, influenza, digestive disorders and other diseases in clinical [7]. The multiple components contained in the Atractylodis rhizome contributed to its beneficial effects including a series of sesquiterpenoids, oligosaccharides, polysaccharides, monoterpenes, polyacetylenes, phenolic acids, and steroids [7]. A. lancea produced in the Maoshan area of Jiangsu Province is the geo-authentic herb, which is famous for the higher oxygenous sesquiterpenoid content and diversity. However, for ecological destruction and overexploitation, the geo-authentic herb has been endangered in the last few years [8]. Therefor many cultivated varieties were also used as Atractylodis rhizome in the market. There have been studies showed that there are obvious differences of volatile components and oxygenous sesquiterpenoid in A. lancea and cultivated varieties [3,6].
Carbohydrates is one of important medicinal compounds of Atractylodis rhizoma, which possesses various and high biological activities, for instance, enhancing immunity, releasing stress, regulating the enzyme activity, promoting probiotics growth and relieving liver injury etc. [9]. In our previous study, a series of oligosaccharides DP 3-15 were separated from A. lancea and identified as FOS, which consisting of linear β-(2,1)linked fructofuranosyl units mostly carry a terminal single α-(1-2)-linked α-glucopyranosyl unit [10][11][12]. It is a well-recognized class of prebiotics [13][14][15], which are dietary ingredients that cannot be digested by humanproduced digestive enzymes, yet they provide a health benefit to the host mediated by selectively stimulating the growth and/or activity of one or a limited number of host gut microbiota [16]. It is becoming increasingly clear that FOS can reduce or prevent gastroenteritis, inflammatory bowel disease (IBS), reduce the risk of colon cancer and reduce potentially pathogenic gastro-intestinal bacteria [17]. These activities are in line with the quality of Atractylodis rhizome for both natural and cultured products. Therefore, the level of FOS is an essential marker for quality evaluation of Atractylodis rhizome.
However, the reference compounds of FOS, which are the standards for qualitative and quantitative analysis, are still difficult to separate and purify for high polarity. Up to date, only FOS with DP 3-5 purified by some laboratories have been commercialized [18,19]. In addition, the chromatographic analysis of FOS is also challenging, since (i) the significantly different solubility, (ii) the absence of high DP standards (DP > 5), and (iii) the absence of chromophores required for UV/Vis detectors. Evaporative light scattering detector (ELSD) is a universal detector, which has been widely applied in detecting analytes without UV absorption. Moreover, ELSD is more sensitive due to gradient elution is available for this detector.
In this study, a HILIC-HPLC method couple with evaporative light-scattering detector (ELSD) was developed for the qualitative and quantitative analysis of 13 FOS in 21 samples including 11 batches of A. lancea and 10 batches of A. chinensis. Based on the quantification information, chemometrics analysis including principal component analysis (PCA), factor analysis (FA) and hierarchical clustering analysis (HCA) were also applied to discriminate different kinds of Atractylodis rhizome and found out the distinctive components to simplify the quality evaluation procedure. Also, the DP 5-based quantification results were compared with the individual quantification results to evaluate the feasibility and accuracy and provide an alternative quantification method for the quality control of oligosaccharides.

Chemicals and materials
The materials of A. lancea and A. chinensis were collected from different cultivation locations in China, which are listed in Table 1. Species identification was performed by Associate Prof. Guang-ping Lv. Voucher specimens of the samples were deposited at the School of Pharmaceutical Sciences, Nanjing Tech University, Nanjing, China. The herbs were dried below 50 °C and then ground into fine powder before use.
Inulin-type FOS with DP 3-15 (all purities determined by HPLC-ELSD were more than 90%) were separated and purified in our laboratory (Fig. 1). The structures were confirmed by comparing their methylation analysis, MS, and NMR data with the literature [20][21][22]. Acetonitrile for HPLC was purchased from Merck (Darmstadt, Germany). Deionized water for HPLC was prepared by an Aquaplore 2S system (Aquapro, USA). All other chemicals and reagents were of analytical grade.

Preparation of standard solutions
Mixed standard stock solution containing inulin-type FOS (DP 3-15) was prepared in 60% ethanol. The concentrations of DP 3-15 were about 3 mg/mL. The standard stock solution was stored in a refrigerator at 4 °C before use. Working standard solutions were prepared from the stock solution by dilution with the appropriate volume of 60% ethanol.

Sample preparation
Atractylodis rhizoma powders (1 g), passed through the 40-mesh sieve, was extracted by refluxing under the optimized conditions. In brief, the extraction was operated with 20 mL extract solvent of ethanol-water (60:40, v/v) for 2.5 h at 80 °C. All extractions were performed with magnetic stirring. After extraction, appropriate volume of 60% ethanol was added to the extract in order to complement the weight. Then, the extracted slurry was centrifuged at 3500 r/min for 20 min. And 1 mL of supernatant was filtered through a 0.45 μm filter before injection into the HPLC system for analysis.

HPLC-ELSD analysis
Chromatographic analyses were performed on a Thermo UltiMate 3000 liquid chromatography system (Thermo, USA), equipped with a quaternary pump, an on-line degasser, an auto-sampler, a column temperature controller and UV detectors, coupled with an Alltech ELSD 6000 (Alltech, Deerfield, IL, USA). Data processing was carried out with Dionex ™ Chromeleon ™ 7.0 software. A gas gen-

Calibration curves, LOD and LOQ
Standard stock solutions containing reference compounds (DP 3-15) were prepared and diluted to appropriate concentrations for the establishment of calibration curves. At least six concentrations of 13 analyte solutions were injected in triplicate, and then the calibration curves were constructed and their linear ranges were determined. Since ELSD response was nonlinear and the calibration of ELSD could be constructed by a double logarithmic plot as in previous reports, the calibration curves were established by plotting the logarithm peak area versus the logarithm concentration (mg/mL) of each analyte detected by ELSD. The stock solutions containing reference compounds were diluted with 60% ethanol to appropriate concentrations, aliquots of the diluted solutions were injected into HPLC for determining the limit of detection (LOD) and limit of quantification (LOQ). The LOD and LOQ for each analyte under the present chromatographic conditions were determined at a signal-to-noise ratio (S/N) of about 3 and 10, respectively.

Precision, repeatability, accuracy, and stability
Precision of the HPLC-ELSD method was determined by evaluating the repeatability of intra-day and inter-day measurements at low, middle and high level. Three known concentration solutions of 13 standards were prepared. For intra-day variability test, the mixed standards solutions were analyzed for six replicates within 1 day, while for inter-day variability test, the solutions were examined in duplicates for three successive days. Variations were expressed by the relative standard deviations (RSD) for intra-and inter-day. Repeatability of the developed method was confirmed with preparation and analysis of six parallel samples. In order to check the accuracy of the developed method, recovery of the investigated components was carried out. Known amounts of the mixed standards were added into the samples (0.25 g) and then the mixture was extracted, processed, and analyzed in accordance with methods mentioned above. Three replicates were performed for the test. Average recoveries were determined by the equation recovery (percent) = 100 × (observed amount − original amount)/ spiked amount, and relative standard deviation (RSD, percent) = 100 × standard deviation (SD)/mean. Stability was tested and analyzed at 0, 2, 4, 8, 16, 24, and 48 h.

Data analysis
HCA, PCA and FA were performed by SPSS 16.0 for windows (SPSS Inc., Chicago IL, USA), which comprises a number of procedures-graphical, statistical, reporting, processing and tabulating procedures that can be used for fast data analysis. A method named as Ward's method was applied and Squared Euclidean distance, a pattern similarity measure, was finally selected as measurement for HCA.

Optimization of sample preparation
The extraction solvents, extraction solvent volume, extraction time, which directly influence the extraction yields of FOS, were optimized. Refluxing extraction procedure was optimized using sample A. lancea (A9). The amount of 18 investigated compounds (DP 3-20) was used as the marker for evaluation of extraction efficiency. Different extraction solvents (40, 60, 80, and 100% ethanol) were investigated using univariate approach. Figure 2 showed that the extraction efficiency of FOS by 40, and 60% ethanol were obviously higher than that of 80, and 100% ethanol. In order to reduce the extraction of watersoluble impurities such as polysaccharides and other macromolecules, the 60% ethanol was adopted for sample preparation [18,23]. Different ratio of solid to liquid, 1:10, 1:15, 1:20 and 1:25, were compared for the extraction of FOS in A. lancea. The results suggested that the extraction efficiency of FOS was obviously higher when solvent volume was 20-fold to sample material (v/w). We also have investigated the effect of extraction time (1.5, 2, 2.5, and 3 h). It turned out that the optimal extraction time of FOS was 2.5 h, whereas longer extraction time could not increase extraction yields of FOS. Therefore, the optimized refluxing extraction method was solvent, 60% ethanol; solvent volume, 20-fold to sample material (v/w); extraction time, 2.5 h.

Method validation
The method of HPLC-ELSD was validated for the quantification. Linearity, regression, and linear ranges of the 13 analytes were determined by the developed HPLC method. The calibration curves with the R 2 , linear range and regression equation, LOD and LOQ of 13 targeted analytes are listed in Table 2. The results indicated the calibration curves of the 13 analytes had good linearity (R 2 > 0.9993) between the analytes concentrations and their peak areas within the test ranges, and their LODs and LOQs were in the ranges of 0.01-0.04 mg/mL and 0.04-0.14 mg/mL, respectively (Table 2). Moreover, for intra-day precision at low, middle and high concentration, RSDs of the peak area were less than 4.43% and RSDs of the retention time were less than 0.50%; for inter-day precision, the corresponding RSDs were less than 10.91% and 0.48% (Table 3). The recoveries of 13 analytes were between 83.16 and 99.66%, and the analytes were stable in 60% ethanol solution during the tested period (Table 4). Repeatability with RSD less than 2.66% suggested that the developed method had a good repeatability for the quantitative evaluation of the analytes. These results showed that the developed HPLC-ELSD method was sensitive, precise and accurate for quantitative analysis of 13 inulin-type FOS (DP 3-15) in Atractylodis rhizome.

Quantitative determination of FOS in plants
The established HPLC-ELSD method was applied for analysis of inulin-type FOS in different samples of A. lancea and A. chinensis. Typical HPLC-ELSD chromatograms of mixed standards and crude extracts were shown in Fig. 3. The identification of the investigated compounds was carried out by comparison of their retention time with those obtained by injecting standards under the same conditions. Besides the DP 3-15, FOS with higher DP were also found in A. lancea and A. chinensis (Fig. 3b, c). DP 16-18 were identified by chemical standards separated in our lab, which the amount was limited and therefore did not validate for the quantification method. DP 3-15 were validated for the quantification analysis. The 13 investigated oligosaccharides in A. lancea and A. chinensis were well separated using the developed HPLC method. The content of these 13 investigated oligosaccharides were calculated based on individual calibration curve (Method 1), and the content were summarized in Table 5. As shown in Table 5 and Fig. 4, the content of FOS in A. chinensis were relatively higher than that of A. lancea. The content of the 13 analytes differed among these two species, which might lead to variances in the pharmacologic actions, even their therapeutic effects. Thus, determination of multiple components is   essential for the quality evaluation of Atractylodis rhizome for both natural and cultured products.

Chemical characteristics of A. lancea and A. chinensis
In order to evaluate the differences of FOS in A. lancea and A. chinensis and identify the distinctive ingredients, PCA was performed. The PCA was carried out by using the content of 13 analytes derived from quantification analysis. The scatter points (Fig. 5a) showed that A. lancea and A. chinensis were well separated. The results indicated that significant difference existed in A. lancea and A. chinensis, which probably induce the different pharmacological effects. FA was employed to identify the distinctive ingredients and reduce the markers for quality control of A. lancea and A. chinensis. Eventually, 2 markers including content of DP 12 and total content of FOS 3-15 were identified as the main distinctive elements, which were available for distinguishing and quality control of A. lancea and A. chinensis (Fig. 5b).
Hierarchical cluster analysis of 21 selected samples of Atractylodis rhizome was performed using content of DP 3-15 (Fig. 6a) and 2 selected markers based on FA (Fig. 6b), respectively. Their results were very similar, A. lancea and A. chinensis were well separated, which were also in accordance with that of PCA. Therefore, the content of DP 12 and total content of FOS 3-15 could be used as markers for quality control of Atractylodis rhizome. Actually, the separation of high DP FOS such as DP 12 is difficult and time-consuming, therefore an alternative method for accurate determination of total content of FOS 3-15 is an ideal means to solve the shortage of reference compound.

Quantification method assessment
Take the advantage of similar structure of FOS and homogeneous response on ELSD [24,25], we have calculated the FOS 3-15 based on the calibration curve of DP 3, 4 and 5, respectively (Method 2). In order to evaluate the method feasibility, percent difference (PD), and cosine similarity between two vectors Cos(θ), were employed [26,27]. The calculation of PD is as follows: where x1 and x2 are the content produced by Methods 1 and 2.
The calculation of Cos(θ) is as the following equation: where X and Y are the content produced by Methods 1 and 2, and n is the number of data sets [26,27]. Based on the calibration curve of DP 3 and 4, the average PDs of FOS 3-15 were up to 26.00% and 25.84%, then the corresponding Cos(θ) were as low as 0.99797 and 0.99785, respectively. Therefore, we calculated the FOS 3-15 based on the calibration curve of DP 5. The results showed that the average PDs of 13 analytes were all less than 13.76% and the Cos(θ) were 0.99986, 1.00000, 1.00000, 1.00000, 1.00000, 1.00000, 1.00000, 1.00000, 0.99998, 0.99995, 0.99993, 0.99996       (Table 5), which demonstrated that the similarities of pairwise arrays between Methods 1 and 2 were very high. In addition, the HCA analysis based on the total content of FOS 3-15 calculated by DP 5-based quantification (Fig. 6c) also showed high degree of consistency with the results based on individual calibration curve (Fig. 6a). The results indicated that DP 5-based quantitative method could be adopted to simplify the analytical procedure for quantification of all the investigated components.
Actually, the variation of analytes' content among different samples derived from natural material was usually high, which significantly indicated that the quality control of FOS in Atractylodis rhizome was crucial. Also, the variation of content was much higher than the quantification errors using DP 5 for calculation and the quantification errors should be acceptable for quality control of Atractylodis rhizome.

Conclusions
In this study, a HPLC-ELSD quantification method based on the individual calibration curve (Method 1) has been developed. Coupled with chemometrics analysis, we demonstrated that 2 markers including content of DP12 and total content of DP 3-15 could be used as the main distinctive elements, which were available for distinguishing and quality control of Atractylodis rhizome. Actually, the separation and purification of high DP FOS, such as DP 12, is still a challenge because of high polarity. However, the quantitative evaluation of total content of DP 3-15 based on low DP standards  Table 1 is relatively easy to implement. DP 5 could be purified on Bio-Gel P-2 column eluted with water at the flow rate of 0.3 mL/min according to our previous work and at present DP 5 has been commercialized. Thus, DP 5-based qualification of multiple analytes (Method 2) is especially suitable for the quality control of Atractylodis rhizome. In addition, using relative content of 13 oligosaccharides based on Method 2 as input data matrix, A. lancea and A. chinensis were well separated, which were also in accordance with that of HCA based on Method 1. Therefore, DP 5-based qualification and quantitative evaluation of FOS in Atractylodis rhizome was demonstrated to be credible for solving the shortage of reference compounds and improving the quality control of FOS in Atractylodis rhizome as well as other herbs containing FOS.  Table 1 Fig. 6 Dendrograms of hierarchical cluster analysis for 21 tested samples of A. lancea and A. chinensis based on a 13 investigated FOS calculated by individual calibration curve (Method 1), b 2 markers (content of DP 12, and total content of FOS 3-15) derived from their FA and c 13 investigated FOS calculated by calibration curves of DP 5 (Method 2). The sample codes were the same as in Table 1