Simultaneous quantification of bioactive components in Chinese herbal spirits by ultra-high performance liquid chromatography coupled to triple-quadrupole mass spectrometry (UHPLC–QQQ–MS/MS)

Background The Chinese medicinal wine made from herbal medicines became prevalent among Chinese people. The Chinese herbal spirit is composed of several herbal extracts, and has the certain health functions, such as anti-fatigue and immune regulation. The quality evaluation of Chinese herbal spirit is greatly challenged by the enormous and complex components with great structural diversity and wide range of concentration distribution. Methods An ultra-high performance liquid chromatography coupled to triple quadrupole mass spectrometry (UHPLC-QQQ-MS/MS) with multiple reaction monitoring (MRM) method was developed to simultaneously determine forty-three bioactive components in the Chinese herbal spirits produced by year 2014 and 2018. Results Quantitative results showed that 11 components, i.e.., puerarin (5), purpureaside C (7), daidzin (8), echinacoside (9), acteoside (15), epimedin B (22), epimedin C (23), icariin (24), eugenol (27), chikusetsusaponin iva (30) and Z-ligustilide (40), significantly decreased along with the increasing years of storage, while 5 compounds, i.e.., geniposidic acid (1), protocatechuic acid (2), crustecdysone (14), daidzein (18) and icariside I (35), were basically stable in all samples across the years. Concusion The established method allowing to simultaneously determined 43 components with wide structural diversity and trace amounts will facilitate the quality control research of Chinese herbal spirits. Supplementary Information The online version contains supplementary material available at 10.1186/s13020-021-00435-0.


Introduction
functions of anti-fatigue and immune enhancement [1]. It was shown to strengthen the immune system in the Shen-yang deficient rats, which was associated with the activation of hypothalamic-pituitary-adrenal axis [2]. Furthermore, a clinic study indicated that the Chinese herbal spirit could relieve both physical fatigue and mental fatigue of the patients with the fatigued subhealth status [3].
According to enterprise criterion, total saponins, total flavonoids and icariin are current chemical markers for quality control of this herbal spirit. However, its quality evaluation is greatly challenged by the enormous and complex components with great structural diversity and wide range of concentration distribution. More than 150 ingredients were identified from the herbal spirit, including flavonoids, saponins, alkaloids, phenylethanoids, coumarins, anthraquinones and volatile oil, etc. To date, gas chromatography-mass spectrometry (GC-MS) [1,4,5] and liquid chromatography-mass spectrometry (LC-MS) [6], have been established to quantify the volatile flavoring substances and active ingredients in the spirit. However, the very limited quantity of ingredients were quantified. Ultra-high performance liquid chromatography coupled to triple quadrupole mass spectrometry (UHPLC-QQQ-MS) operated in multiple reaction monitoring (MRM) mode is an effective quantification method owing to its well-known high sensitivity and specificity, which could avoid the interference from the background matrix [7][8][9]. It has been successfully utilized to quantify bioactive components in the complex systems [10][11][12][13][14][15]. In the current study, therefore, an UHPLC-QQQ-MS method was developed to simultaneously quantify 43 bioactive components in the Chinese herbal spirit samples produced by year 2014 and 2018, and their concentrations in different peoduction years were also compared.

Materials and chemicals
The Chinese herbal spirit samples produced by year 2014 (n = 20) and 2018 (n = 20) were kindly provided by Jing brand Co. Ltd (Hubei, China). The voucher specimens were deposited at room temperature and shielded from light at Institute of Chinese Medical Sciences, University of Macau, Macao. Forty-three reference compounds listed in Table 1 were purchased from Baoji herbest Bio-Tech Co. Ltd (Baoji, China). Their chemical structures were shown in Additional file 1: Figure S1, and the purities of reference standards were over 98 % as confirmed by HPLC-UV. HPLC-grade acetonitrile and methanol were purchased from Merck (Darmstadt, Germany). Formic acid was obtained from Aladdin Industrial Inc. (Shanghai, China). Deionized water was purified using a Millipore Milli-Q purification system (Bedford, MA, USA).

Sample preparation and standard solution preparation
An aliquot of 1mL of the spirit was diluted with the equivalent volume of acetonitrile, and vortexed for 1 min. The mixture was centrifuged at 14,800 rpm for 20 min. After the centrifugation, the supernatant was filtered through a 0.22 μm filter (PVDF Millex-GV, 13 mm, Millipore) prior to the quantitative analysis. All of 43 reference standards were accurately weighed and dissolved in methanol to prepare individual stock solutions at the concentrations ranging from 0.58 to 2.29 mg mL − 1 . The mixed standard solution was prepared by mixing appropriate volumes of the individual stock solutions and further diluted to a series of proper concentrations with methanol.

LC-MS/MS analysis
Forty bathes of the spirit samples were analyzed by an ACQUITY ™ UPLC system (Waters, Milford, MA, US) coupled with a Xevo TQD triple-quadrupole tandem mass spectrometry (QQQ-MS/MS, Waters Co., Manchester, UK). Chromatographic separation was implemented on an ACQUITY UPLC BEH C 18 column (150 × 2.1 mm, 1.7 μm). The mobile phases were consisted of 0.1 % (v/v) aqueous formic acid solution (phase A) and acetonitrile containing 0.1 % formic acid (phase B) at the flow-rate of 0.3 mL min −1 using a gradient elution program as follows: 5-40 % B at 0-12 min, 40-100 % B at 12-16 min, isocratic 100 % B for 2 min, and the re-equilibrated by 5 % B for 3 min. The sample injection volume was set at 2 µL. The temperatures of column and injector were set at 35 ℃ and 8 ℃, respectively.
Data acquisition was performed by a Xevo TQD QQQ-MS equipped with an electrospray ionization (ESI) using MRM mode. MS was operated in either the positive (ESI+) or negative mode (ESI-) to obtain satisfactory MS response for 43 investigated compounds due to their different properties. The MS and MS/MS spectra of each compound were acquired using the mixed standard solution. The optimized MS parameters were as follows: capillary voltage, 3.5 kV (positive ion mode) and − 3.0 kV (negative ion mode); source temperature, 140 ℃; desolvation gas flow and temperature, 650 L h −1 and 350 ℃; cone flow, 50 L h −1 . The ion transitions, cone voltage, and collision energy for each compound were optimized and shown in Table 1. All instrumentations were synchronized and controlled by Waters Masslynx software (version, 4.1).

Method validation
To evaluate sensitivity and precision of the established UHPLC-QQQ-MS/MS method, the linearity, limit of detection (LOD), limit of quantitation (LOQ), precision and recovery of 43 analytes were tested. The calibration curve of each compound was constructed by plotting the peak areas against the concentrations using the mixed standard solution at a series of concentrations. The precision was examined by calculating intra-and inter-day variations of each analyte using the mixed standards for five replicates within a day and three consecutive days. The LOD and LOQ for each analyte were estimated at signal-to-noise ratio (S/N) of about 3 and 10, respectively. The accuracy of the established method was assessed by spike recovery experiments. A known amount (equal to the content for each analyte in the sample) of the mixed standards was spiked into the random spirit sample (S-2018-08). The sample was prepared with six replicates and analyzed by the method mentioned above. The recovery (%) was calculated as the following equation:

Statistical analysis
The concentrations of 43 analytes in the spirit samples were presented as mean ± standard deviation (SD). The difference between groups was assessed by student t-test using a GraphPad Prism package (version 6.0, San Diego, CA, USA), and a p-value of less than 0.05 was considered statistically significant. Orthogonal partial least squaresdiscriminant analysis (OPLS-DA), a supervised multiple regression, was conducted to discriminate the spirit samples manufactured in different years according to the levels of investigated analytes by SIMCA-P software (version 14.1, Umetrics, Umeå, Sweden).

LC-MS/MS method development
The quality control research related to this Chinese herbal spirit mainly focused on the determination of volatile components by GC-MS [1,4,5], and very limited number of nonvolatile bioactive components were quantified by LC-MS [6]. Due to consisting of multiple herbal extracts and the faint-scented Xiaoqu liqueur, the herbal spirit is a very complex matrix containing the numerous organic and inorganic compounds with wide range of concentrations. Thus, the simultaneous quantitation of forty-three compounds with various chemical types encounters the great challenge in short running time using UHPLC. MRM is a highly specific technique for quantifying the targeted analyte, regardless of baseline chromatographic separation. The targeted analyte in chromatographic co-elution could be accurately quantified if they have different MS or MS/MS characteristics. However, the co-eluted analytes may cause the potential mutual ionization suppression in ESI, leading to the low MS response. Therefore, it is also necessary to optimize the chromatographic conditions, including column and mobile phase, to achieve the high sensitivity and fast separation in LC-MS/MS analysis. Three UHPLC columns, such as BEH C 18 column, BEH HILIC column and HSS T3 C 18 column, were examined. As results, an ACQUITY BEH C 18 column was most suitable for the separation of the targeted compounds in the sample owing to the best resolution and the most peak capacity. Furthermore, several types of mobile phases, including methanol/water and acetonitrile/water system supplemented with various modifiers, were tested. The results shown that 0.1 % aqueous formic acid solution / acetonitrile with 0.1 % formic acid was the optimum mobile phases to obtain the chromatogram with the best resolution. It is critical to design ion transition of each analyte, including precursor ions and their corresponding product ions, in MRM analysis. The full scan was used to select the precursor ions using their reference standards, and the dominated fragment ion in daughter scan was chosen as the corresponding product ion (Fig. 1). Taking sagittatoside A (32, MW = 676.24), a main component derived from Epimedii folium, as an example, the  (Fig. 1 a, b). The mass difference between parent and product ions was 308.09, which corresponded to the loss of one glucose and one rhamnose units [15]. Herein, the ion transition (m/z 677.33 → 369.24) was selected to quantify sagittatoside A in the liqueur by MRM. Likewise, the ion transition (m/z 417.36 → 255.22) was optimized to determine daidzin (8, Fig. 1 c, d). To achieve maximum signal, both positive and negative ion modes were tested and compared. The results indicated that all compounds, except Chikusetsusaponin iva (30), shown the higher sensitivities in the positive ion mode than in the negative mode. For example, MRM chromatograms of daidzein (18) and puerarin (5) in positive ion mode were remarkably higher than that in negative mode (Fig. 2a, b). However, the higher MS intensity of chikusetsusaponin iva was observed in negative mode (Fig. 2c). Furthermore, the transition of chikusetsusaponin iva in the negative mode were further optimized. We found that multiple ion monitoring (MIM), i.e. m/z 793.57 → 793.57, was more suitable for the detection of chikusetsusaponin iva, compared to other transitions, such as m/z 793.57 → 631.67 (Fig. 2d).
Additionally, cone voltage (CV) and collision energy (CE), the important factors that affect the sensitivity of UHPLC-QQQ-MS/MS analysis, were also optimized for each analyte using reference standard. CV and CE were optimized from 10 V to 50 V with a step of 10 V and from 5 V to 50 V with a step of 5 V, respectively. Taking sagittatoside A as an example, the signal intensity of product ion m/z 369.24 increased along with CV from 10 V to 30 V or CE from 5 V to 15 V, then decreased with the increasing voltages (Fig. 3a). Therefore, 30 V of CV and 15 V of CE were chosen for the quantification of sagittatoside A. Likewise, 30 V of CV and 20 V of CE were used to determine the daidzin (Fig. 3b). In a similar manner, the mass spectrometry conditions of all analytes were optimized and listed in Table 1. Under the optimized LC-MS/MS conditions, 43 analytes were well separated and detected in 18 min. The representative MRM chromatograms of

Method validation
The established LC-MS/MS method was validated by a series of experiments, including linearity, sensitivity, precision, and accuracy. As shown in Table 1, the calibration curves of all analytes exhibited good linear regression (R 2 ≥ 0.9911) within the wide dynamic range. The LODs and LOQs of analytes were less than 8.54 and 27.7 ng mL − 1 , respectively. Other than that of geniposide was 23.1 and 46.2 ng mL − 1 , respectively, which was higher than other compounds. The overall intra-day and inter-day variations were lower than 4.89 % and 5.81 %, respectively. Additionally, the developed method had the acceptable accuracy with recoveries ranging from 90.3 to 111.7 %. Taken together, the proposed LC-MS/MS method is sensitive, precise and accurate for the simultaneous determination of these 43 compounds in the Chinese herbal spirit produced in different years.

Quantification of 43 compounds in the Chinese herbal spirit
The validated method was successfully applied to quantify 43 compounds in the liqueur samples. As shown in Table 2  In order to further visualize the difference between samples, OPLS-DA, a supervised multivariate data analysis, was further constructed to characterize the differences between groups on the basis of the levels of 35 analytes determined. As illustrated in Fig. 5a, the spirit samples produced at 2014 and 2018 were unambiguously segregated into two tight clusters (R 2 X = 0.826, R 2 Y = 0.96, Q 2 = 0.887), suggesting that the great alteration in the investigated compounds were presented. The cumulative values of R 2 X, R 2 Y, and Q 2 were close to optimal value of 1.0, indicating the established models with excellent predictive capability and fitness [16]. To identify the differentiated components that contribute most to the group separation, the differentiated compounds were selected by S-plot derived from the constructed OPLS-DA. Eleven compounds with variable importance in the projection (VIP) of more than 1 (VIP > 1) and p value of less than 0.05 were highlighted in the S-plot (Fig. 5b). The group separation of the samples with the different production years could also be clearly classified in OPLS-DA score plot according to the contents of these 11 highlighted components (Fig. 5c). Compared Oleanolic acid 0.14 ± 0.07 0.09 ± 0.03** to the spirit produced in 2018, the contents of the highlighted compounds, puerarin (5), purpureaside C (7), daidzin (8), echinacoside (9), acteoside (15), epimedin B (22), epimedin C (23), icariin (24), eugenol (27), chikusetsusaponin iva (30) and Z-ligustilide (40), significantly decreased in the samples produced in year 2014 (Fig. 6a). Geniposidic acid (1), protocatechuic acid (2), crustecdysone (14), daidzein (18) and icariside I (35) around the origin of S-plot were also highlighted as the stable compounds, which shown the comparable contents with less variation between samples from two production years (Fig. 6b). Additionally, a heatmap according to the relative contents of 35 detedcted analytes was constructed to display the changes in the analyte levels between groups (Fig. 6c).