Metabolic Analysis of Four Cultivars of Liriope platyphylla

Liriope platyphylla (Liliaceae), a medical plant distributed mainly in China, Taiwan, and Korea, has been used traditionally for the treatment of cough, sputum, asthma, and neurodegenerative diseases. The present study involved the metabolic profiling of this plant and reports spicatoside A accumulation in four different varieties of L. platyphylla (Cheongyangjaerae, Seongsoo, Cheongsim, and Liriope Tuber No. 1) using HPLC and gas chromatography time-of-flight mass spectrometry (GC–TOFMS). A total of 47 metabolites were detected in the different cultivars using GC–TOFMS-based metabolic profiling. The resulting data were subjected to principal component analysis (PCA) for determining the whole experimental variation, and the different cultivars were separated by score plots. Furthermore, hierarchical clustering, Pearson’s correlation, and partial least-squares discriminant analyses (PLS-DA) were subsequently performed to determine significant differences in the various metabolites of the cultivars. The HPLC data revealed that the presence of spicatoside A was detected in all four cultivars, with the amount of spicatoside A varying among them. Among the cultivars, Liriope Tuber No. 1 contained the highest amount of spicatoside A (1.83 ± 0.13 mg/g dry weight), followed by Cheongyangjaerae (1.25 ± 0.01 mg/g dry weight), Cheongsim (1.09 ± 0.04 mg/g dry weight), and Seongsoo (1.01 ± 0.02 mg/g dry weight). The identification of spicatoside A was confirmed by comparing the retention time of the sample with the retention time of the standard. Moreover, the Cheongsim cultivar contained higher levels of phenolic compounds—including vanillic acid, quinic acid, gallic acid, chlorogenic acid, caffeic acid, and benzoic acid—than those of the other two cultivars. On the other hand, the levels of amino acids were higher in the Seongsoo cultivar. Therefore, this study may help breeders produce new varieties with improved nutraceutical and nutritional qualities.


Analysis of Spicatoside A
Using metabolic profiling, the present study investigated spicatoside A accumulation in four different varieties of L. platyphylla (Cheongyangjaerae, Cheongsim, Liriope Tuber No. 1, and Seongsoo). The HPLC data revealed the presence of spicatoside A in all four cultivars, with the amount of spicatoside A varying among them. Among the cultivars, Liriope Tuber No. 1 contained the highest amount of spicatoside A (1.83 ± 0.13 mg/g dry weight), followed by Cheongyangjaerae, Cheongsim, and Seongsoo (1.25 ± 0.01, 1.09 ± 0.04, and 1.01 ± 0.02 mg/g dry weight, respectively) ( Table 1). The identification of spicatoside A was confirmed by comparing the retention time of the sample with the retention time of the standard.

Analysis of Phenolic Compounds
A total of five phenolic compounds, including catechin, gallic acid, chlorogenic acid, caffeic acid, and benzoic acid, were identified and quantified in four different varieties of L. platyphylla by HPLC ( Table 2). Cheongsim contained higher concentrations of total phenolics (461.11 ± 4.96 µg/g dry weight), which were 1.36 times higher than the concentrations found in Cheongyangjaerae (340.03 ± 4.29 µg/g dry weight). Moreover, a comparison of individual phenolics showed that Cheongsim contained higher amounts of gallic acid, caffeic acid, and chlorogenic acid than those of the other varieties, respectively. On the other hand, Liriope Tuber No. 1 contained a higher level of benzoic acid than the other varieties.

GC-TOFMS-Based Metabolic Profiling and Multivariate Analysis
Gas chromatography-mass spectrometry (GC-MS) is a useful tool with high sensitive reproducibility and handling capacity for analysis via fast spectra accumulation times [19]. Through the GC-TOFMS technique, the levels of primary metabolites including carbohydrates, amino acids, and organic acids, and secondary metabolites including phenolic compounds, were determined. A total of 47 metabolites were identified from four L. platyphylla). The metabolites identified in L. platyphylla were quantitated using selected ions and normalization, based on the signal intensity of IS (Tables  S1 and S2). The quantitation data of the 47 hydrophilic metabolites were then subjected to principal component analysis (PCA), to explore the differences in metabolite profiles among the four L. platyphylla varieties. Metabolic profile analysis using the GC-MS technique and multivariate analysis have been used to explore metabolic differentiation among various genotypes. In particular, PCA has been successfully utilized as a preliminary stage in multivariate analysis to identify the patterns in complex experimental data [20]. The score plot from PCA results presented an overview of the differences among the four L. platyphylla cultivars, and the loading plot enabled correlation examination among the 47 metabolites ( Figure 1). All the samples representing dissimilar colors were differently clustered, which were not differentiated by the two highest-ranking components. The components accounted for 64.1% of the total variance within the data set. Notably, the first component, accounting for 42.8% of the total variation, resolved the measured metabolite profiles of Cheongsim and Liriope Tuber No. 1 from the other two cultivars, Seongsoo and Cheongyanghaerae. This separation was mainly attributable to organic acids and amino acids; the corresponding loading vector showed that most of the organic acids, with the exception of three phenolic acids (p-hydroxybenzoic, vanillic, and sinapic acids), were predominant in Cheongism and Liriope cultivars, whereas the amino acids, with the exceptions of arginine and glutamine, differentiated the Seongsoo and Cheongyanghaerae cultivars from the other two. The PCA results showed that most amino acids were higher in Cheongyanghaerae compared with the other three cultivars, while Cheongsim had higher levels of organic acids than the other cultivars.
Metabolites 2019, 9, x FOR PEER REVIEW 4 of 12 hydroxybenzoic, vanillic, and sinapic acids), were predominant in Cheongism and Liriope cultivars, whereas the amino acids, with the exceptions of arginine and glutamine, differentiated the Seongsoo and Cheongyanghaerae cultivars from the other two. The PCA results showed that most amino acids were higher in Cheongyanghaerae compared with the other three cultivars, while Cheongsim had higher levels of organic acids than the other cultivars. To maximize the separation between the cultivars, we conducted partial least-squares discriminant analysis (PLS-DA) ( Figure 2). PLS-DA is a well-established chemometric approach that separates groups of observations by rotating the PCA such that a maximum separation among classes, here L. platyphylla cultivars, is obtained. Identification of the first component aided in resolving the measured metabolite profiles of the four cultivars. The contribution of variables in the projection could be explained using variables important in the projection (VIP). Shikimic acid, one of the important building blocks employed in the biosynthesis of phenylpropanoids, was the most important for creating a prediction for L. platyphylla classification. The levels of shikimic acid and phenolic acids were higher in Cheongsim than other cultivars. For a cross validation, the total data should be divided into a training set and a test set. The 47 metabolites from the 12 samples were To maximize the separation between the cultivars, we conducted partial least-squares discriminant analysis (PLS-DA) ( Figure 2). PLS-DA is a well-established chemometric approach that separates groups of observations by rotating the PCA such that a maximum separation among classes, here L. platyphylla cultivars, is obtained. Identification of the first component aided in resolving the measured metabolite profiles of the four cultivars. The contribution of variables in the projection could be explained using variables important in the projection (VIP). Shikimic acid, one of the important building blocks employed in the biosynthesis of phenylpropanoids, was the most important for creating a prediction for L. platyphylla classification. The levels of shikimic acid and phenolic acids were higher in Cheongsim than other cultivars. For a cross validation, the total data should be divided into a training set and a test set. The 47 metabolites from the 12 samples were divided into 8 training set samples and 4 test set samples. The model shown in Figure 2C has Q 2 = 0.949. In the partial least-squares (PLS) prediction model, a cross-validated correlation coefficient (Q 2 ) > 0.9 indicates an excellent model.  Figure 2C has Q2 = 0.949. In the partial least-squares (PLS) prediction model, a cross-validated correlation coefficient (Q2) > 0.9 indicates an excellent model.  To gain insights into the relationships among the concentrations of the 47 hydrophilic metabolites identified in L. platyphylla, we performed a hierarchical cluster analysis (HCA) using Pearson's correlation results on the data sets. Correlation analysis is one of the effective techniques to determine the strength of a relationship among quantitative samples. Using this method, significant correlations can be found between diverse metabolites involved in closely related pathways in a biological system [21]. The HCA results showed the degree of correlation among the identified metabolites. Compounds from the same class in a biological system containing amino acids, organic acids, or carbohydrates clustered together, which are marked by a dotted box in Figure 3. One group included all the amino acids except glutamine, whereas another group contained most organic acids, including the metabolites related to the TCA-cycle metabolism. The Pearson correlation coefficients among leucine, isoleucine, and valine, which are branched amino acids, were higher than 0.7. Similarly, a significant positive relationship was detected between serine and glycine (r = 0.7492, p = 0.005), and phenylalanine and tryptophan (r = 0.7284, p = 0.007), which are biosynthetically linked metabolites. Among the organic acids, the metabolites related to the TAC cycle, such as citric, succinic, fumaric, and malic acids, showed a strong correlation (r > 0.9, p < 0.0001). In addition, most of the identified carbohydrates, such as galactose, glucose, mannose, raffinose, sucrose, and trehalose, formed a group. However, fructose and xylose clustered with a group of organic acids rather than with other sugars. These results are in agreement with the loading plot of PCA results ( Figure 1B), which described the robustness of the present experimental system and the potential of GC-TOFMS-based metabolic profiling, combined with chemometrics, as a useful tool for investigating metabolic links in biological systems.
Sugars, including xylose, fructose, mannose, galactose, glucose, sucrose, raffinose, and trehalose, were the most abundant metabolites in the four different cultivars. The total sugar level in Cheongsim was higher compared with the other cultivars; specifically, the levels of carbohydrates, including fructose, mannitol, sucrose, glucose, mannose, and galactose, were significantly higher in Cheongsim than in the other three cultivars. On the other hand, the levels of raffinose and trehalose were higher in Cheongyangjaerae. Liriope Tuber No. 1 showed the lowest level of total carbohydrates. Furthermore, a total of 18 amino acids were identified in all the cultivars, and the total amount of amino acids in Seongsoo was higher than in the other cultivars; however, the levels of only two amino acids, namely pyroglutamate and serine, were statistically higher compared with the other cultivars.
Similarly, the levels of proline, beta-alanine, and alanine were significantly higher in Cheongyangjaerae than in the other three cultivars. These results were consistent with a previous study, which suggested Cheongsim as a potential source of amino acids among the cultivars, based on HPLC-based amino acid analysis. A total of four intermediates of the tricarboxylic acid (TCA) cycle were detected in these varieties. Liriope Tuber No. 1 and Cheongsim contained larger pools of citrate, succinate, fumarate, and malate, which were related to the levels of pyruvate. Likewise, the levels of phenolic compounds, including vanillic acid, quinic acid, gallic acid, chlorogenic acid, caffeic acid, and benzoic acid, were higher in Liriope Tuber No. 1 and Cheongsim, compared with the other two cultivars.

Discussion
The vague boundary between primary and secondary metabolism, comprising the intermediate components biosynthesized by primary metabolism, which are also involved in secondary metabolism pathways, suggests interactions between both metabolisms. In this study, the connection between primary and secondary metabolites was investigated in four cultivars of L. platyphylla using GC-TOFMS and HPLC.
Carbohydrates play an important role in plant development, growth, and morphogenesis, and function as energy sources and metabolic precursors [22]. The lowest level of total carbohydrates in Liriope Tuber No. 1 reflected the requirement of carbons and energy to produce higher levels of spicatoside A. This finding is consistent with previous studies explaining the relationship between carbohydrates and secondary metabolites: Yeo et al. (2018) reported that the lower levels of carbohydrates in the sprouts of Vigna unguiculata exposed to blue light-emitting diodes reflected the higher production of phenolics and carotenoids [23]. Similarly, the total amount of carbohydrates in red flowers was lower than that in violet and white flowers of Rhododendron schlippenbachii due to

Discussion
The vague boundary between primary and secondary metabolism, comprising the intermediate components biosynthesized by primary metabolism, which are also involved in secondary metabolism pathways, suggests interactions between both metabolisms. In this study, the connection between primary and secondary metabolites was investigated in four cultivars of L. platyphylla using GC-TOFMS and HPLC.
Carbohydrates play an important role in plant development, growth, and morphogenesis, and function as energy sources and metabolic precursors [22]. The lowest level of total carbohydrates in Liriope Tuber No. 1 reflected the requirement of carbons and energy to produce higher levels of spicatoside A. This finding is consistent with previous studies explaining the relationship between carbohydrates and secondary metabolites: Yeo et al. (2018) reported that the lower levels of carbohydrates in the sprouts of Vigna unguiculata exposed to blue light-emitting diodes reflected the higher production of phenolics and carotenoids [23]. Similarly, the total amount of carbohydrates in red flowers was lower than that in violet and white flowers of Rhododendron schlippenbachii due to their anthocyanin production [24].  described that red lettuce contains lower levels of sugars and higher levels of phenolics and carotenoids than green lettuce [25], and that there is a negative correlation between carbohydrates and phenolic compounds in red radishes [6]. Furthermore, the enhancement of alkaloid production in cell cultures of the opium poppy, derived from treatment with a fungal elicitor, showed a more rapid depletion of carbohydrate pools [26]. Similarly, TCA intermediates play a major role in cellular catabolism and flavonol synthesis [27]. Therefore, the larger pools of TCA intermediates in Liriope Tuber No. 1 and Cheongsim suggest a greater accumulation of phenolic compounds. This finding is consistent with a previous study [20].
Amino acids are considered protein building blocks and are required for the biosynthesis of secondary metabolites, including phenolics and glucosinolates [28]. In this study, the Seongsoo cultivar contained the highest total amount of amino acids, which is consistent with a previous study reporting that the total amount of amino acids in Seongsoo was 1.76, 1.77, and 2.18 times higher than in Cheongyangjaerae, Cheongsim, and Liriope Tuber No. 1 cultivars, respectively [29].
Spicatoside A is a representative chemical found in L. platyphylla. In this study, spicatoside A was detected and quantified in four cultivars of L. platyphylla. Among them, Liriope Tuber No. 1 contained the highest amount of spicatoside A. Similarly, previous studies reported the identification of spicatoside A in L. platyphylla cultivated in South Korea; Cho et al. (1998) identified spicatoside A as well as spicatoside B using NMR [30]. Similarly, Shin (2002) detected and quantified both chemicals in Liriope Tuber No. 1 grown in South Korea using HPLC [31]. Furthermore, Lee et al. (2009) reported the changes in spicatoside A content in L. platyphylla tubers, with varying drying processes [32]. Choi et al. (2015) established the high performance counter-current chromatography (HPCCC) coupled with evaporative light scattering detection (ELSD) method for the separation of spicatoside A and spicatoside D from L. platyphylla [33].
In this study, the Liriope Tuber No. 1 cultivar had the highest level of spicatoside A, and the Cheongsim cultivar had the highest levels of phenolic compounds. On the other hand, the Seongsoo cultivar had the highest total amount of amino acids. Therefore, this study suggests that the Liriope Tuber No. 1 and Cheongsim cultivars could be a potential source of spicatoside A and phenolic compounds, respectively. As well, the Seongsoo cultivar could be considered a good source of amino acid for human intake. Furthermore, this study may help breeders establish a breeding strategy for new varieties with improved nutraceutical and nutritional qualities.

Plant Material
Four varieties of L. platyphylla (Seongsoo, Cheongsim, Cheongyangjaerae, and Liriope Tuber No. 1) were cultivated in a greenhouse at the Cheongyang Boxthorn Experiment Station, Cheongyang, South Korea. The tubers of the different varieties were harvested on 30 April 2016. The tubers were immediately placed in liquid nitrogen at −196 • C and then freeze-dried at −80 • C for at least 72 h (Ilshin Lab Co., Ltd., Daejeon, South Korea). Afterwards, mortars and pestles were used to grind the dried samples for further chemical analysis.

Spicatoside A Extraction and Quantification
Ten grams of dried sample powder was placed in a 250-mL flask and extracted with 80% MeOH (200 mL × 3) by reflux (80 • C, 30 min). After vortexing, the mixture was filtered using filter paper (Whatman No. 2) and left to stand for at least 1 h at room temperature. It was then centrifuged at 15,000 rpm at 4 • C for 5 min. Subsequently, the supernatant was evaporated in vacuo. The residue was dissolved in 1 mL of MeOH and filtered with a 0.45-µm pore-size hydrophilic polyvinylidene difluoride syringe filter (Ø 13 mm, Cat. no. 6779-1304, Whatman Int. Ltd., Clifton, NJ, USA) into an HPLC vial. HPLC analysis of spicatoside A was performed on a Waters 1525 Binary HPLC (Miami, FL, USA) equipped with a photodiode array detector. Isolation was achieved on a SunFire C-18 stainless steel column (2.1 × 50 mm, 5 µm thickness) and absorption determinations were generated at 210 nm, with an oven temperature of 30 • C and a flow rate of 1.0 mL/min. A methanol (solvent A) and acetonaitrile (solvent B) isocratic solution (A:B (30:70), 15 min) was used.

Phenylpropanoid Extraction and Quantification
Phenolic compounds were extracted using previously published procedures [34]. First, 0.2 g of dried sample powder was placed in a 15-mL tube and extracted with 80% MeOH (200 mL × 3). After sonication for 1 h, the mixture was filtered using filter paper (Whatman No. 2) and the extract was then evaporated. Next, the residue was dissolved in 2 mL of MeOH and filtered with a syringe filter into a brown vial. The HPLC analysis system, gradient program, and condition were carried out as detailed in a previous study [27]. We employed an NS-4000 HPLC system (Futecs, Daejeon, Korea), equipped with a C 18 [34,35]. Quantitation was performed using the respective calibration curves.

Metabolic Analysis
Extraction of hydrophilic compounds was carried out as described previously [36]. For the polar metabolite extraction, 10 mg of the powdered samples were measured and added to 1 mL of water:chloroform:methanol (2.5:1:1). Ribitol (60 µL, 0.2 mg/mL) was used as an internal standard (IS). Methoxime derivatization was executed through the addition of methoxyamine hydrochloride in pyridine (80 µL) and vigorous shaking for 90 min at 30 • C. Subsequently, trimethylsilyl etherification was carried out by adding N-methyl-N-trimethylsilyltrifluoroacetamide (80 µL). After incubation at 37 • C for 30 min, GC-TOFMS analysis was performed as described by Kim et al. (2013) [19]. The conditions for GC-TOFMS analysis were set as follows: split ratio, 1:25; injector temperature, 230 • C; flow rate of helium through the column, 1.0 mL/min; detector voltage, 1700 V; scanned mass range, 85-600 m/z; ion-source temperature, 200 • C; transfer line, 250. Furthermore, the temperature program was set as follows: initial temperature of 80 • C for 2 min, followed by an increase to 320 • C at 15 • C /min, and a 10 min hold at 320 • C. Leco ChromaTOF software was used for the peak detection and automated deconvolution of reference mass spectra. In-house libraries for standard compounds [35,36] and the National Institute of Standards and Technology (NIST) were used to identify the metabolite. The results were filtered with retention time, signal-to-noise ratio (>5:1), and mass spectral matching (based on a match >700). As a result, a total of 47 metabolites were identified (i.e., the metabolomics standards initiative (MSI) level 1) [37] (Table S2 and Figure S1). Concentration calculations of all metabolites were based on the ratios determined from the peak area of each metabolite/the peak area of the IS (Tables S1 and S3).

Statistical Analysis
Quantification data obtained from GC-TOFMS were subjected to PCA and PLS-DA (SIMCA-P version 12.0; Umetrics, Umeå, Sweden) to assess the similarity between groups of multivariate data. The data file was scaled with unit-variance scaling (standardization) before all the variables were subjected to PCA. The PCA output is made up of score plots for visualizing the contrast between different samples and loading plots to describe the cluster separation. Furthermore, the SAS 9.4 software package (SAS Institute, Cary, NC, USA) was used for Pearson correlation analysis. Hierarchical clustering analysis and heat map visualization of the calculated correlation coefficients were executed using the software Multi-Experiment Viewer version 4.9.0 (http://www.tm4.org/mev/). The HPLC data were analyzed using SAS 9.4 software, applying an analysis of variance (ANOVA) evaluation and a Duncan's multiple range test (DMRT).

Conclusions
To our knowledge, no previous studies have reported a comprehensive description of multiple metabolites in different cultivars of L. platyphylla. This study indicated that the Liriope Tuber No. 1 cultivar contained the highest level of spicatoside A and low amounts of carbohydrates, which are indicative of the carbon requirement and energy demand. Moreover, the observed high levels of TCA intermediates involved in flavonol metabolism were in agreement with the high levels of phenolic compounds. In this study, the Cheongsim cultivar contained the highest levels of phenolic compounds, while the highest total amount of amino acids was investigated in the Seongsoo cultivar. In contrast, the Liriope Tuber No. 1 cultivar contained the higher level of spicatoside A. Therefore, this study provides valuable information to help breed new cultivars with desired specific metabolite compositions.