Abstract

A method, quantitative analysis of multicomponents by single marker (QAMS), was established and fully verified based on high-performance liquid chromatography (HPLC) for simultaneous determination of six chromone indicators of Saposhnikoviae Radix (SR). In the present study, cimifugin (C), 5-O-methylvisamminol (V), hamaudol (H), and their corresponding glycosides, prim-O-glucosylcimifugin (GC), 4′-O-β-D-glucosyl-5-O-methylvisamminol (GV), and sec-O-glucosylhamaudol (GH), were selected as bioactive constituents and indicators for the quality evaluation of SR. GV was chosen as the unique reference standard, and relative correction factors (RCF) between GV and the other five chromones were calculated. The feasibility of QAMS for the analysis of chromones was investigated by comparing with the traditional external standard method (ESM). Furthermore, the method was proven to have accuracy (96.98%–102.50%), repeatability (RSD <3%), stability (RSD <3%), precision (RSD <3%), and desirable linearity (R2 ≧0.9999). Subsequently, 55 batches of commercial SR from different regions were determined by QAMS, and their contents were analyzed by principal component analysis (PCA), correlation analysis, and hierarchical cluster analysis (HCA), respectively. Based on the results, a more refined quality standard of commercial SR was proposed: SR was qualified when the total contents of six chromones were greater than 3 mg·g−1. Furthermore, SR could initially be regarded as a superior medicine when it satisfied both conditions at the same time: the total content of GC, C, GV, V, GH, and H was greater than 8 mg·g−1, and the proportion of the total content of C, V, and H was greater than 10%. This study demonstrated that the quality of SR could be successfully evaluated by the developed QAMS method; meanwhile, valuable information was provided for improving the quality standard of SR.

1. Introduction

Saposhnikoviae Radix (SR), known as “FangFeng” in China, is derived from the dried root of the plant Saposhnikovia divaricate (Turcz.) Schischk, belonging to the family Umbelliferae. This herb is pungent in flavour, sweet, and lukewarm, and enters the bladder, liver, and spleen meridian. According to the theory of Traditional Chinese Medicine (TCM), SR has significant effects on dispelling wind to relieve exogenous syndrome, removing dampness to kill pain, and stopping spasms [1]. As an important ingredient in many traditional Chinese prescriptions such as Yu-Ping-Feng-San, Fang-Feng-Tong-Sheng pill, and Tong-Xie-Yao-Fang [2, 3], SR has great application value. Furthermore, many pharmacological studies have indicated that a number of curative effects, including antipyretic, analgesic, anti-inflammatory, antibacterial, antitumor, antiallergic, and antioxidation, are existed in SR [48].

SR has several components: chromones, polysaccharides, coumarins, volatile oils, and other components [912]. It is worth mentioning that chromones are the most representative components in SR. On the one hand, there are extremely abundant content [13]. On the other hand, it was closely related to the pharmacological efficacy of SR, for example, anti-inflammatory, analgesic, and antioxidation [1416]. With the development of SR market, increasing importance has been given to the quality control of SR. According to the record of Chinese Pharmacopoeia (2015 edition), prim-O-glucosylcimifugin (GC) and 4′-O-β-D-glucosyl-5-O-methylvisamminol (GV) were selected as quality-control indicators. However, each Chinese herbal medicine is an integrated complex with diversity of active components, contributing to particular efficacy through synergy and mutual effect based on the theory of TCM [17]. Thus, using a few components to control the quality is insufficient for the complex botanical products and traditional Chinese medicines. In line with the importance status of the chromones mentioned above, multiple representative chromone components could be selected as indicators to evaluate the quality of SR. According to phytochemical studies, GV, GC, sec-O-glucosylhamaudol (GH), cimifugin (C), 5-O-methylvisamminol (V), and hamaudol (H) are widely existing in SR [18, 19], and their pharmacological effects are significant and most related to the efficacy of SR [20, 21]. Therefore, when G, GC, V, GV, H, and GH are selected as new quality-control indicators, the chemical characterization, medicinal function, and inner quality of SR could be represented comprehensively. Some studies have shown that the methods of refluxing and ultrasonic were always used for the extraction of chromones [22, 23]. Considering the polarity and solubility of the six chromone indicators, the heating refluxing method could be applied for sample preparation, and the methanol could be used as the extraction solvent based on an optimization of sample preparation process in this work.

In many studies, the contents of multicomponents are usually determined by the external standard method (ESM). In this method, the reference standards are necessary that we need to spend more and more time and economic cost to separate and purify [24]. As an alternative method, quantitative analysis of multicomponents by single marker (QAMS) is only requiring a single reference standard to simultaneously determine the contents of multicomponents, which is more effective and appropriate for the quality control (QC) [25]. When some reference standards are unstable, low in abundance, or hard to extract from the plant, QAMS could not only reduce the cost but also reduce the difficulty in preparation [26, 27]. Besides, this method could improve the practicability of QC and expand the application for herbal or botanical products effectively. In QAMS, the content of internal standard could be obtained directly by HPLC and the other components could be calculated by using multiple conversion factors. Hence, the relative correction factor (RCF) is a critical parameter in the content computational formula about analytes. As the result of molar absorptivity of different analytes are often different, RFC plays a role of calibration when a single reference is used to determine multicomponents [28, 29]. In ESM, the concentration of analyte (Ck) can be calculated by the ratio of the peak areas of analytes in sample solution (Ak) to the peak area of its corresponding standard solution in a unit concentration (As/Cs), as shown in the following equation:

In QAMS, Ck should be calibrated by RCF of each analyte (fk) based on the calculation of ESM. The formula is as follows:

Importantly, the value of fk is calculated by the ratio of peak areas in a unit concentration of standard substance (As/Cs) to analyte (Ak/Ck) as follows:

It is worth mentioning that the final value of RCF is usually the average value of multiple RCFs by a series of determining under different concentration levels of internal referring substance [30]. Since six control components, i.e., C, GC, V, GV, H, and GH, are used, ESM will cause high costs and complicated operations. Therefore, the QAMS method is used to compute the contents of six chromones, and GV is used as the internal standard for its strong representative nature, high stability, high content, and significant pharmacological activities [31, 32].

In the present study, a new substitute method named QAMS was applied for simultaneous determination of six chromones in SR. As the reference substance, the content of GV was determined by HPLC, and the contents of C, GC, V, H, and GH were calculated with RCF based on the intrinsic function and the proportional relationship between GV and these five chromones. The feasibility could be verified by comparing the results with the ESM, and this method was validated in terms of linearity, accuracy, precision (instruction precision and intermediate precision), and stability, referring to some reliable references [33]. Subsequently, 55 batches of commercial SR were determined and a more comprehensive and reliable quality evaluation standard of SR was preliminary inferred by principle component analysis (PCA), correlation analysis, and hierarchical cluster analysis (HCA), respectively.

2. Materials and Methods

2.1. Apparatus and Chromatographic Analysis

Analyses were primarily performed by using a Waters HPLC System (Waters Crop, Milford, MA, USA) equipped with a 1525 binary pump solvent management system, 2998 PDA detector, 2707 automatic sampling device, and Breeze 2 workstation. Two additional different HPLC instruments were used: One was high performance liquid chromatography (Waters Crop, Milford, MA, USA) equipped with 1525 binary pump solvent management system, 2489 UV detector, and Breeze 2 workstation. Another was a Waters Alliance e2695-2998 HPLC system (Empower workstation, Waters Crop, Milford, MA, USA). HPLC separation was carried out on a CAPCELL PAK C18 column (4.6 mm × 150 mm, 5 μm). Column temperature was set at 25°C, and inject volume was 10 μL. The mobile phase consisted of methanol (A) and 0.3% formic acid aqueous solution (B). The gradient elution was programmed at a flow rate of 1.0 mL·min−1 as follows: 0–12 min, 32% A; 12–40 min, 32%–50% A; 40–50 min, 50%–70% A; 50–52 min, 70% A. The detection wavelength was set at 254 nm.

2.2. Chemicals and Materials

Fifty-five batches of SR were collected from different regions in China, which were identified by Professor Zhang Yuan from the Beijing University of Traditional Chinese Medicine and proved to be the dried root of Saposhnikovia divaricate (Turcz.) Schischk following the method described in Chinese Pharmacopoeia (2015 edition) [1]. GC, C, GH, and GV were obtained from Chengdu Mansite Biotechnology Co., Ltd. (Chengdu, China). The purity (≧98%) of these reference standards was assumed as provided by the suppliers. The other two compounds, V and H, were separated and purified in our lab and the purity was identified to be of not less than 98% (determined by HPLC). The structures were determined on the basis of UV, MS, and NMR data and confirmed by comparison with data from the literature. The chemical structures of all standards are shown in Figure 1.

Methanol, acetonitrile, and formic acid of HPLC grade were purchased from Thermo Fisher Scientific Inc. (Waltham, MA, USA), and other reagents (Beijing Chemical Industry Factory) were of analytical grade. HPLC grade water was prepared using a Pall Cascada IX system (Pall, USA). All other reagents were of analytical grade.

2.3. Preparation of Mixed Standard Solutions

Substances of C, GC, H, GH, V, and GV were weighed precisely and dissolved into methanol to prepare the mixed stock solutions of reference standards with the concentrations of 0.1220 mg·ml−1, 0.2096 mg·ml−1, 0.0369 mg·ml−1, 0.1636 mg·ml−1, 0.0808 mg·ml−1, and 0.3340 mg·ml−1, respectively. Then, a series of concentrations of calibration standard solutions were produced by diluting the mixed stock solution (dilution factor = 1, 2, 5, 10, 20, 50, and 100) with the same methanol. The solutions were stored at 4°C in a refrigerator and filtered through a 0.45 μm membrane filter before injection. All samples being injected into HPLC system were prepared right before analysis.

2.4. Preparation of Sample Solutions

An appropriate amount of the samples to be tested was crushed and passed through an 80-mesh screen. Next, 0.25 g of the sample powder was precisely weighed and placed in a stuffed flask with accurate addition of 10 ml of methanol, subjected to heating reflux for 120 min. It is worth noting that the stuffed flask was weighed before and after refluxing, then added the solvent to keep the weight equivalent at room temperature if required. By filtering through the 0.45 μm filter and discarding the first 2 ml, the remaining filtrate was used as the test solution of samples.

2.5. Method Validation
2.5.1. Specificity

The mixed reference standard solution and test solution were separately injected into HPLC under the optimized chromatographic conditions (Section 2.1).

2.5.2. Linearity

The prepared stock standard solutions mentioned above (Section 2.3) with a series of appropriate concentration levels were used for HPLC based on the chromatographic conditions (Section 2.1), respectively. The limits of detection (LOD) and quantification (LOQ) were measured based on a signal-to-noise ratio (S/N) at about 3 and 10, respectively.

2.5.3. Precision

To ensure the validity of this newly developed method, the tests of instrument precision and intermediate precision were performed. For instrument precision, test solutions prepared (Section 2.4) were examined by HPLC for six replicates within one day, and on the purpose of detecting intermediate precision, the prepared test solutions were injected into HPLC by different operators with different instruments on different dates.

2.5.4. Stability

The same tested solutions, which were prepared (Section 2.4) and placed at room temperature, were injected into HPLC at different time points (0, 2, 4, 8, 12, and 24 h).

2.5.5. Repeatability

Six parallel sample solutions with the same batch were prepared following the method (Section 2.4) individually and determined by HPLC according to the chromatographic conditions (Section 2.1).

2.5.6. Accuracy

Six copies of the same batch of SR powder (0.125 g) with known content were weighed, respectively. Then, a certain amount of control standard solution was added into the samples according to the proportion of sample content to reference substance, about 1 : 1. Preparation and determination of six sample solutions were conducted in parallel referring to the method as described in Sections 2.1 and 2.4. Recoveries were computed.

2.6. QAMS Method
2.6.1. Calculation of RCF

Seven concentration levels of mixed standard solution were prepared (Section 2.3) and injected into HPLC under the chromatographic conditions (Section 2.1), respectively. Besides, the chromatographic peak areas of each component were recorded.

2.6.2. Durability Test of RCF

Three different instruments (as listed in Section 2.1), three kinds of chromatographic columns (Capcell Pak C18, Water SunFire C18, and Water Symmetry C18) (4.6 mm × 150, 5 μm), different flow rates (0.9, 1.0, and 1.1 ml·min−1), and different column temperatures (25, 30, and 35°C) were used to investigated the influence of different conditions on RCF.

2.6.3. Location of the Chromatographic Peak of Measured Component

For better authentication as well as convenience to quality control of the commercial SR, the chromatographic peak positions of GC, C, GV, V, GH, and H were investigated using different instruments and different columns.

2.6.4. Comparison of the Results between QAMS and ESM

In order to assess and validate QAMS feasibility of multicompounds in SR, the contents of C, GC, V, GV, H, and GV were determined by ESM and QAMS in 15 batches, respectively. For ESM, the determination of the six chromones was carried out with six reference standards (GC, C, GV, V, GH, and H), whereas for QAMS, the results were based on the nature of the calculation of fx, the intrinsic function, and the proportional relation between the selected reference analyte and other analytes. The content of the selected internal substance (GV) was determined like ESM, and then the contents of the other five chromones were calculated in accordance with relative conversion factors between analytes and the internal substance [34].

2.7. Application and Data Analysis

The developed QAMS method was applied for the quantitative assessment of 55 batches of commercial SR from different regions. The contents of GC, C, GV, V, GH, and H were determined and then analyzed by PCA, correlation analysis and HCA, respectively. Meanwhile, the further analysis of three chromone glycosides (C, H, and V) was carried out to make clarification of their importance for the overall quality of SR. The figures presented were developed by exploration of the analysis function using SPSS 22.0 software package.

3. Results and Discussion

3.1. Method Validation
3.1.1. Specificity

As shown in Figure 2, the analytes had good separation because it has no interference in the corresponding position of the six components, and the target peaks of the test solution corresponded to the peaks of reference standard solution according to retention time in the chromatogram. So, it indicated that this method had specificity.

3.1.2. Linearity

The standard curves of six reference substances were established by using the chromatographic peak area (y) as the vertical axis and the concentration of the reference solution (x) as abscissa, respectively. There was a good linearity as the result of all the correlation coefficients (R) which were not less than 0.999 over the concentration range. Furthermore, LOD and LOQ of six substances were also calculated as shown in Table 1.

3.1.3. Precision

The areas of chromatographic peak and RSD of six compounds were recorded and calculated, respectively. In Table 2, the RSD results of six compounds for the instrument precision were in the range of 0.42%–0.75%. Apart from this, the RSDs which were calculated by intermediate precision were all lower than 3%. It indicated that this method has a good precision.

3.1.4. Stability

The peak areas of GC, C, GV, V, GH, and H were recorded, and the values of RSD were 0.84%, 0.56%, 054%, 0.57%, 0.34%, and 0.69%, respectively. Hence, the sample solution was stable at room temperature within 24 h (Table 2).

3.1.5. Repeatability

The results showed that average mass fractions of GC, C, GV, V, GH, and H were 1.478, 0.808, 2.510, 0.173, 1.053, and 0.151 mg·g−1, respectively. The RSD of corresponding average mass fraction was 1.03%, 0.73%, 0.49%, 0.80%, 0.96%, and 1.35% which proved that this quantitative method was of good repeatability.

3.1.6. Accuracy

There was favorable accuracy because the average recovery rates of six marker compounds were varied in the range of 96.98%–102.5%. Meanwhile, RSD values of recovery rates for each compound were lower than 3% totally.

3.2. QAMS Method
3.2.1. Calculation of RCF

GV was selected as the internal standard, and the values of RCF (fx) for other five indicators were computed in different concentrations according to the equation (3) mentioned above. The average RCF of each compounds was shown as follows: fGC = 1.047 (RSD% = 2.32), fC = 0.6489 (RSD% = 0.24), fv = 0.5909 (RSD% = 0.30), fGH = 0.7223 (RSD% = 0.68), and fH = 0.7223 (RSD% = 0.28).

3.2.2. Durability Test of RCF

The RCF of five chromones in different conditions (instruments, chromatographic columns, flow rates, and column temperatures) were obtained, and the RSDs were all less than 5%, which could clearly demonstrate that the RCF calculated by the proposed method has good durability and system suitability for routine testing.

3.2.3. Location of the Chromatographic Peaks of Measured Components

The chromatographic peak position was identified by the relative retention time (RRT), which was calculated according to the following equation:where tRk is the retention time of measured components, tRS is the retention time of internal reference, and △tRks is the difference of retention time in both. Among them, the chromatographic peak position of GV was explicitly identified and designated as the reference peak in the SR samples. RRT was measured between GV and the other five components by different instruments and columns at the same time, and their RSDs were all less than 5%. It indicated that the calculation of RRT was stable and could be used for identifying the chromatographic peaks location of measured components.

3.2.4. Comparison of the Results between QAMS and ESM

The contents of 15 batches of commercial SR were determined by the QAMS method and ESM. The results are shown in Table 3. Relative error (RE) was built between the two component variables to examine the deviations between QAMS and ESM. By comparing two sets of contents of five components between QAMS and ESM, respectively, their content variations were found to be within the range of 5%. It met the requirement of Chinese Pharmacopoeia. To assess the consistency of the results, correlation analysis was used to evaluate the similarity between the QAMS method and ESM. Correlation coefficient value is a commonly used parameter in the similarity evaluation. The larger the values, the higher the similarity of the target sample will be. When they are equal to 1, the targets are identical. In this work, the data, as shown in Table 4, were above 0.900, which indicated that there were no significant differences between the QAMS and ESM, and the identified RCF and parameters of chromatographic peak location for measured chromones of SR were reliable. In conclusion, QAMS can be applied in the determination of six chromones.

3.3. Application and Data Analysis
3.3.1. Sample Analysis and Characteristics of Six Chromone Compounds in 55 Batches of Commercial SR

In 55 batches of commercial SR from different regions, the contents of six chromones, GC, C, GV, V, GH, and H, were determined by QAMS. Results are listed in Table 4. It was observed that the maximum total content of six compounds was 11.22 mg·g−1 in no. S53, while the minimum total content of those was 3.047 mg·g−1 in no. S18. Such a wide concentration variance of these 55 batches of commercial SR may be attributed to a variety of factors, including plant sources, genetic variation, and geography differences. To further verify the relationships among the samples and evaluate the variation of six compounds, PCA, coefficient analysis, and HCA were performed using the SPSS 22.0 software (IBM, USA).

3.3.2. PCA and Correlation Analysis

The contents of GC (X1), C (X2), GV (X3), V (X4), GH (X5), and H (X6) of 55 batches of commercial samples were subjected to PCA. The results are shown in Figure 3. A two-component PCA model was established accounting for the accumulated variation of 80.260%, where the first principal component (Z1) was 62.614% and the second (Z2) was 17.646%. According to the component score coefficient matrix (Table 5; Figure 4), every coefficient between Z1 and six indicators was significant, which could indicate that Z1 represented the total content of six components. That is to say, each component was indispensable for the quality evaluation of SR. Furthermore, Z2 was mainly reflecting X6 as the result of their coefficient which was the largest (0.733). The comprehensive score (Z) for each batch sample could be obtained as follows:

The details of 55 batche samples are shown in Table 6. Through correlation analysis of results from every batch sample, a good correlation between Z and the total content of six components could be included in accordance with the correlation coefficient which was 0.875. Thus, it was feasible to evaluate the quality of SR comprehensively if the total content of six components was used as indicator. Meanwhile, we preliminary set that the total content of six chromones should not be less than 3 mg·g−1 for qualified SR based on the above determined content range (3.047–11.22 mg·g−1).

3.3.3. HCA

The HCA was applied to analyze the concentrations of GC, G, GV, V, GH, and H. The result indicated that 54 batches of commercial SRs (the 51st batch of SR was self-contained, not considered) were divided into two categories as shown in Figure 5. On the one hand, Cluster 1 consisted of 45 batch samples (34, 55, 14, 54, 1, 11, 4, 15, 36, 39, 10, 40, 41, 42, 35, 18, 47, 7, 28, 43, 46, 23, 27, 50, 16, 29, 2, 3, 5, 21, 9, 12, 8, 38, 19, 44, 13, 37, 24, 48, 26, 20, 32, 33, 30), in which their total content of six components was not more than 8 mg·g−1. On the other hand, the remaining 9 batch samples (6, 31, 25, 22, 45, 53, 17, 49, 52) were classified for second category, which were greater than 8 mg·g−1 of the indicator. Contacting with the results of Section 3.1.1, it could strongly prove that taking the total content of GC, G, GV, V, GH, and H as an indicator can reasonably, comprehensively, and objectively evaluate the quality of SR.

For the quality evaluation of SR, the chromone aglycone should not be ignored because this plays an important role in the pharmacological activity of SR. In addition, as shown in Table 5, the correlation coefficients of C, V, and H were significant (0.866, 0.900, and 0.0599), respectively. Using the ratio of the total content of C, V, and H to the total content of six chromones as an indicator, 53 batches of commercial SRs were divided into two categories by HCA (the 13th and 37th batches were self-contained, not considered). The details are shown in Table 7 and Figure 6. Cluster 1 included 41 batches (32, 54, 14, 19, 20, 4, 33, 15, 48, 25, 30, 51, 49, 52, 26, 1, 6, 31, 45, 24, 41, 11, 42, 27, 10, 17, 50, 34, 44, 53, 55, 22, 2, 18, 3, 7, 39, 29, 40, 36, 35), in which the ratio was greater than 10%. Meanwhile, the remaining 12 batch samples (12, 23, 5, 38, 9, 8, 47, 28, 46, 21, 16, 43) were attributed to Cluster 2, where the ratio was less than 10%. It was indicated that these three chromone aglycones were so important that they could be used as quality indicators for the quality evaluation of SR.

Last but not least, based on the results of the twice cluster analysis, samples of the Cluster 2 (6, 31, 25, 22, 45, 53, 17, 49, 52) in the first HCA were all contained to the Cluster 1 (32, 54, 14, 19, 20, 4, 33, 15, 48, 25, 30, 51, 49, 52, 26, 1, 6, 31, 45, 24, 41, 11, 42, 27, 10, 17, 50, 34, 44, 53, 55, 22, 2, 18, 3, 7, 39, 29, 40, 36, 35) in the second HCA. Thus, we could initially infer a conclusion that the SR could be regarded as a superior medicine when the total concentration of GC, C, GV, V, GH, and H was greater than 8 mg·g−1. Meanwhile, the proportion of aglycons (C, V, and H) was greater than 10%. In this study, 9 batch samples (6, 31, 25, 22, 45, 53, 17, 49, 52) were the superior medicine with the six chromone content which was greater than 8 mg·g−1, and meanwhile, the total content of C, V, and H was greater than 0.8 mg·g−1.

3.4. Optimization of Sample Preparation

Optimization of extraction methods, solvents, solvent volume, and extraction time were investigated by single-factor test to obtain the best extraction efficiency. The results revealed that extraction efficiency of refluxing was more efficient than ultrasonic extraction for analytes, so the remaining experiment was carried out by refluxing, and methanol was chosen as the best solvent by comparing with various solvents including methanol, 70% methanol, and 50% methanol. In addition, the extraction volume (5, 10, and 15 ml) and the extraction times (60, 120, and 240 min) were tested subsequently. Consequently, the optimal sample preparation parameter was refluxing with 10 mL methanol for 120 min as for 0.25 g powder of SR.

4. Conclusion

A method named QAMS was established to evaluate the quality of SR based on routine HPLC apparatus. In this method, GV was chosen as the internal standard to determine the RCF between GV and other five chromones (GC, C, V, GH, and H) of SR. QAMS was accurate and feasible for the quality evaluation according to the results of method validation, and no significant difference existed in the content results obtained by QAMS and ESM. Using this developed method, 55 batches of commercial SR from different regions were determined, and the results were analyzed by PCA, correlation analysis, and HCA, respectively. It could include that chromones played an important role for the quality of SR; meanwhile, the total content of GC, C, GV, V, GH, and H was used as the evaluation indicator that was comprehensive, objective, and reliable. Other than this, the SR could be regarded as qualified medicine if the total content of six chromones was not less than 3 mg·g−1. Moreover, the importance of chromone aglycones for the quality evaluation of SR was further demonstrated. It could initially infer that it was a superior SR medicine if the total content of GC, C, GV, V, GH and H was greater than 8 mg·g−1; meanwhile, the proportion of C, V, and H was greater than 10%.

All in all, these results provided useful information for the development of commercial SR, and the quality of SR from different origins or different purchase locations was confusing and unstable. Therefore, a compulsive processing standard for SR should be established and standardized. Last but not least, the abovementioned series of analyses will play a positive role in the improvement of the quality evaluation system of SR.

Data Availability

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

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research was supported by the Beijing Municipal Education Commission Common Construction Project (no. 2016022).