Original Research ARTICLE
Comprehensive Quality Assessment Based Specific Chemical Profiles for Geographic and Tissue Variation in Gentiana rigescens Using HPLC and FTIR Method Combined with Principal Component Analysis
- 1Institute of Medicinal Plants, Yunnan Academy of Agricultural Sciences, Kunming, China
- 2College of Traditional Chinese Medicine, Yunnan University of Traditional Chinese Medicine, Kunming, China
Roots, stems, leaves, and flowers of Longdan (Gentiana rigescens Franch. ex Hemsl) were collected from six geographic origins of Yunnan Province (n = 240) to implement the quality assessment based on contents of gentiopicroside, loganic acid, sweroside and swertiamarin and chemical profile using HPLC-DAD and FTIR method combined with principal component analysis (PCA). The content of gentiopicroside (major iridoid glycoside) was the highest in G. rigescens, regardless of tissue and geographic origin. The level of swertiamarin was the lowest, even unable to be detected in samples from Kunming and Qujing. Significant correlations (p < 0.05) between gentiopicroside, loganic acid, sweroside, and swertiamarin were found at inter- or intra-tissues, which were highly depended on geographic origins, indicating the influence of environmental conditions on the conversion and transport of secondary metabolites in G. rigescens. Furthermore, samples were reasonably classified as three clusters along large producing areas where have similar climate conditions, characterized by carbohydrates, phenols, benzoates, terpenoids, aliphatic alcohols, aromatic hydrocarbons, and so forth. The present work provided global information on the chemical profile and contents of major iridoid glycosides in G. rigescens originated from six different origins, which is helpful for controlling quality of herbal medicines systematically.
For centuries, Gentianae Radix et Rhizoma (Longdan in Chinese) has been prominent in treating liver diseases in the traditional Chinese medicine system, one of origins of which is the roots or rhizomes of Gentiana rigescens Franch. ex Hemsl (State Pharmacopoeia Commission, 2015). G. rigescens (GR) belonging to Gentianaceae family, with the Chinese name of Dianlongdan or Jianlongdan, is a perennial herb native to Yunnan Province, China (Editorial Committee of Flora China Nica of the Academy of Sciences of China, 1988). In addition to the document in Chinese Pharmacopoeia, a wealth of clinical experience from Yi ethnic minority group verified its treatment uses of hepatitis and cholecystis (Suyama et al., 2017). The major compounds in GRs include iridoid glycosides and xanthones (Pan et al., 2016), active ingredients responsible to various pharmacological effects such as anti-inflammation, antioxidant, anti-cancer, antiviral, and so forth (Ma et al., 2001; Wang et al., 2013; Guo et al., 2014; Fabiani, 2016; Balkan et al., 2017). Iridoid glycosides belong to the monoterpene family and are divided into iridoid glycosides and seco-iridoid glycosides, depending on the bond between C-7 and C-8 of cyclopentene skeleton. GRs typically contain three seco-iridoid glycosides and one iridoid glycoside: gentiopicroside, sweroside, swertiamarin, and loganic acid (Pan et al., 2015). A number of studies have revealed that these four ingredients, and especially their derivate compounds, function as agents against hepatitis, inflammation, bacteria, free radical, etc. (Kumarasamy et al., 2003; Suyama et al., 2013; He et al., 2015). For example, gentiopicroside act the hepatoprotective effect by changing bile acids metabolism to correct the dyshomeostasis induced by pathogen (Tang et al., 2016), while the swertiamarin and sweroside instead by enhancing the activities of superoxide dismutase and catalase and the level of glutathione (Mihailović et al., 2014). Compounds with bioactive effects have been increasingly important to the quality regulation of herbal medicines in recent years, due to properties of inherence and abundance in favor of detection (Liu et al., 2016; Ding et al., 2017).
It has been widely recognized that environmental factors appear to play important roles in the accumulation of secondary metabolites in plant tissues, a consequence of that adjusting the growth environment may change the transcription and expression of genes related to the biosynthesis (Zoratti et al., 2014; Vuorinen et al., 2015). Zanatta et al. (2017) investigated the phytochemicals in soybeans cultured in two microregions, the result of which indicated that contents of protein, fatty acid, total carotenoid and γ-tocopherol were higher in microregion II where the main features are the relatively low altitude and precipitation. Melito et al. (2016) suggested that chemical profiles of Helichrysum italicum grew in seaside and mountains were different, particularly when plants lived in the farther climatic locations was diverse. The Yunnan area is characterized by the three-dimensional climate on low latitude plateau, where a changing climate and complex terrain have been created naturally as a place for displaying the plant diversity. The strong adaptive capability to different environmental conditions makes GR become a widespread species in Yunnan area, with three kinds of climatic cones i.e., tropical, subtropical, and temperate zones, large scope of altitude, annual average temperature and annual precipitation, and various soil types (Zhang et al., 2012). Maintaining the homogeneous quality of GRs, of which original plants live in complex growth condition, is a challenge. Thus, the quality control and assessment is paramount importance. It is not clear that how the geographical factors influence the content of secondary metabolites in GR, not just gentiopicroside, despite Wu et al. (2017) establishing an evaluation strategy for discriminating GRs collected from different geographical origins by FT-IR coupled with PLS-DA and SVM-GS.
Methods that can reflect the chemical features and accurately quantify target compounds in samples is critical to the quality control and assessment of herbal medicines. Providing the convenience for both qualitative and quantitative analysis, high performance liquid chromatography (HPLC) in combination with different spectrometric detectors has traditionally been the workhorse of herbal medicines' studies (Song et al., 2013). Zhang et al. (2016) simultaneously determined nine active ingredients isolated from Salvia miltiorrhiza and its variety using HPLC-DAD method, followed by revealing latent difference in components via principal components analysis (PCA) and partial least squares discriminant analysis (PLS-DA). One of the most widely utilized techniques to obtain information of chemical structures is Fourier transform infrared spectra (FTIR), especially when combined with pattern recognition methods (Hirri et al., 2016). The most useful chemical constituents for diagnostic monitoring of storage effects on quality control of traded saffron was achieved using FT-MIR coupled to PCA, according to the study by Ordoudi et al. (2014). Moreover, other tools that have found considerable success in herbal medicines uses are gas chromatography coupled with mass spectrometry (GC-MS), ultra-performance liquid chromatography-quadrupole time of flight mass spectrometry (UPLC-QTOFMS), high performance thin-layer chromatography (HPTLC) and so on (Nguyen et al., 2016; Guzelmeric et al., 2017; Liu J. et al., 2017).
In the present work, the quality assessment of GR samples collected from six geographic origins of Yunnan area was carried out using HPLC combined with IR rather than other expensive techniques, meanwhile, in conjunction with the PCA. Gentiopicroside, sweroside, swertiamarin, and loganic acid were selected to be standard compounds for the quantitative analysis of chromatography. This study may provide some potential supports on the reasonable application and exploitation of GR.
Materials and Methods
Chemicals and Plant Materials
The extraction solvent (methanol) and potassium bromide (KBr) were purchased from Tianjin Fengchuan Chemical Reagent Technolodies Co., Lid. (Tianjin, China). Acetonitrile and formic acid (HPLC grade) were obtained from Thermo Fisher Scientific (Fair Lawn, NJ, USA) and Dikmapure (Lake Forest, CA, USA), respectively. Purified water was provided by Hangzhou Wahaha Group (Huangzhou, China). All other chemicals and reagents were analytical grade. The standard compounds (1, loganic acid; 2, swertiamarin; 3, gentiopricroside; 4, sweroside) were purchased from the National Institutes for Food and Drug Control (Beijing, China). Purities of standard compounds were all > 98%.
One hundred and seventy eight specimens of G. rigescens (Gentianaceae) were collected from 15 different sites characterized by varied environmental conditions in Yunan and Guizhou Provinces, China (Table 1). The root, stem, leaf, and flower of each individual were separated. All of samples were authenticated by Prof. Hang Jin (Institute of Medicinal Plants, Yunnan Academy of Agricultural Sciences, Kunming, China), and voucher specimens (DLDYS20121112 and DLDYS20131011) were deposited at the specimen room of Institute of Medicinal Plant, Yunnan Academy of Agricultural Sciences.
Standard solutions were dissolved in HPLC grade methanol individually to achieve a stock with concentration of 3.9 mg/mL (gentiopicroside), 1.1 mg/mL (loganic acid), 0.7 mg/mL (sweroside), and 1.0 mg/mL (swertiamarin). The calibration standards were prepared at 11 levels from 0.003 to 1.9 mg/mL (gentiopicroside), 0.001–0.5 mg/mL (loganic acid), 0.001–0.3 mg/mL (sweroside) and 0.0009–0.5 mg/mL (swertiamarin) for establishment of external standard calibration curves. The calibration curve was obtained by plotting the chromatographic peak area of standard compounds (Y) vs. the corresponding concentrations (X). All standard compounds were stored at −20°C when not in use.
Root, stem, leaf and flower samples were dried at 50°C in a 101A electric thermostatic drying oven (Experimental Instrument Factory, Shanghai, China). Dried samples were milled and sifted through a 60 mesh sieve, respectively. An accurately weighted sample powder (0.025 g) was extracted by 1.5 mL 80% methanol for 45 min. The extract solution was filtered through 0.22 μm membrane filters. Then, the filters collected in auto sampler were analyzed directly by LC systems. 1.2 mg sample powder was mixed evenly 100.0 mg KBr crystal. The mixture was ground and pressed into a tablet.
High Performance Liquid Chromatography (HPLC) Analysis
Analyses were performed on an Agilent 1260 HPLC system (Agilent Technologies, Santa Clara, California, USA) composed of a G1315D diode-array detector (UV-vis, DAD, 190-400 nm), a G1311C VL quaternary gradient pumps equipped with a vacuum degasser, a G1316A thermostatted column compartment and a G1329B ALS auto sampler. The modified version of our published HPLC method (Chu et al., 2016) was used as separation conditions in this study. The separation was carried out on Agilent Intersil-C18 column (150 × 4.6 mm, 5 μm). The mobile phase consisted of 0.1% formic acid in water (A) and acetonitrile (B). The samples were eluted with the following gradient: 0.0–0.4 min, 0–7% B; 0.4–2.5 min, 7–10% B; 2.5–20.0 min, 10–26% B; 20.0–29.0 min, 26–58.3% B; 29.0–30.0 min, 58.3–90% B; 30.0–34.0 min, 90% B. The flow rate was 1.00 mL/min. The injection volume was 10 μL. The detective wavelength was set at 246 nm. The temperature-controlled column oven was set at 30°C. Methanol and water supplemented 0.1% formic acid were degassed by ultrasonication for 30 min to avoid bubbles in solutions prior to analysis.
Fourier Transform Infrared (FTIR) Spectroscopy Analysis
Infrared absorption spectra of samples were recorded using a FTIR spectrometer (Perkin Elmer, Foster City, CA, USA) equipped with a deuterted triglycine sulfate (DTGS) detector. Typically, the accumulation spectra of 16 scans per samples was collected and averaged. Absorption spectra in the area between 4,000 and 400 cm−1, at a resolution of one data point every 4 cm−1, were obtained. Interferences of CO2 and H2O in air were eliminated automatically.
For reference standards, linearity, the limit of detection (LOD), the limit of quantification (LOQ), precision and accuracy were experimentally verified. The LOD and LOQ, signal-to-noise ratios (S/N) of 3 to 10, were determined by serial dilution of each standard solution using the described conditions.
Precision was evaluated by intra- and inter-day variation determined by analyzing mixed standard solutions with known concentration six times within a day and on three consecutive days in triplicate. Accuracy was validated by recovery test that was performed by accuracy adding three different amounts (low, medium, and high spike) of the reference standards to the crude samples. The recovery rate was calculated as follow:
%R = [(measured amount – original amount) / amount added] × 100%.
PCA were performed for assessing the difference of chemical constituents in samples from various geographic origins using FTIR data sets, which were second derivative-transformed before analysis. ANOVA followed by Tukey's test at p < 0.05 were employed to evaluate the statistical significance of differences in the variables among different samples. Correlations between gentiopicroside, loganic acid, sweroside, and swertiamarin were achieved by Pearson's correlation analysis. All assays were unfolded using R 3.4.0 program (R Core Team, 2017).
To assess possible correlations of chemical profiles with environment conditions, annual mean temperature, annual mean relative humidity, and annual mean precipitation data from 1981 to 2010 for the sampling sites of G. rigescens were downloaded from Climatic Data Center, National Meteorological Information Center, China Meteorological Administration.
Result and Discussion
Quantitative Analysis of Four Target Compounds in GR
Figure 1 shows the plot of chromatographic peaks of different tissue samples (roots, flowers, leaves, and stems) from six geographic origins. An appreciable amount of compounds were isolated under the gradient elution by HPLC-DAD, according to the chromatogram. Among them, four target compounds (peak 1–4) were identified unambiguously as loganic acid, swertiamarin, gentiopicroside, and sweroside, respectively, by comparing with the retention times of reference standards. The content of each compound was calculated by their calibration curve, and the result was shown in Figure 2 and Table S1.
Figure 1. HPLC chromatograms plot of different tissue samples of G. rigescens from six geographic origins.
Figure 2. The contents of gentiopicroside, loganic acid, sweroside and swertiamarin in roots (A–D), stems (E–H), leaves (I–L) and flowers (M–P) of G. rigescens from six geographic origins (mg/g, n = 10).
From the results, gentiopicroside (14.55 ± 0.75–53.40 ± 2.34 mg/g) was the predominant ingredient in GR regardless of tissue types and geographic origins, of which the content was significantly more than other constituents. Of four tissues, the root maximized accumulations of gentiopicroside (39.41 ± 2.24–53.40 ± 2.34 mg/g) and loganic acid (5.24 ± 0.71–1.13 ± 1.23 mg/g). Remarkable variations (p < 0.05) were observed in gentiopicroside and loganic acid in tissues from six geographic origins, except for the root which has historically been regarded as the medicinal part for treating various ailments. Contents of gentiopicroside in roots from Dali and Diqing (Figure 2A) and in leaves from Diqing (Figure 2I) were significantly higher than that from other origins, while the low level of that was found in stems from Qujing (Figure 2E). The root still has been regarded as the main site to accumulate gentiopicroside despite the geographical factors changing the contents more or less, combined with the study performed by Qi et al. (2017a,b), which focused on the chemical diversity of GR on different cultivation years. Still, leaves from Diqing contained higher contents of loganic acid compared with others (Figure 2J). Leaves from Dali and Qujing showed significant low contents of this compound. The fairly low contents were found on the rest of two compounds in all samples, particularly the swertiamarin (< 1.78 ± 0.15 mg/g) which even could not be detected in some samples such as stems from Kunmig and Qujing (Figure 2H). Six origins in present work are geographically divided into three large producing areas: Dali is located in the western Yunnan characterized by the high precipitation (Figure S1A); Yuxi, Kunming, and Qujing lie on the central Yunnan featured as the high relative humidity (Figure S1B); Lijiang and Diqing are situated in the northwestern Yunnan famous with the low temperature (Figure S1C). From the tissue point of view, there were high contents of gentiopicroside, loganic acid, and sweroside in leaves from the cold northwestern Yunnan. Then, high levels of gentiopicroside existed in roots and flowers from western Yunnan with abundant precipitation. Combining content data described above, we speculated that the temperature and precipitation are two important factors influencing the accumulation of secondary metabolites. The temperature was the major factor influenced the accumulation of secondary metabolites in the leaf, whereas that in the root, stem, and flower was mainly depended on the precipitation. These results also indicated that secondary metabolites express a tight link to environmental factors, in which the high humidity showed an adverse effect on the accumulation of (seco-) iridoid glycosides.
Gentiopicroside is the major compound in G. rigescens as well as the unique indicator of quality assessment recommended by Chinese Pharmacopoeia. Our results showed that plants from Dali have the potential powerful pharmaceutical effects not only due to high content of gentiopicroside but also loganic acid, swertiamarin, and sweroside.
The total contents of four target compounds were calculated for eliminating the interference of tissues (Figure S2 and Table S1). Contents of gentiopicroside in whole plants from central Yunnan were significantly lower than other two big producing areas (Figure S2A), whereas the relatively high level of loganic acid was observed in samples from northwestern Yunnan with low temperature (Figure S2B). Enhanced water treatment decreased iridoid glycoside concentrations, according to the precipitation patterns established by Jamieson et al. (2013). Martz et al. (2009) indicated that low temperature favored higher contents of iridoids concomitantly reduced the content of flavonols in Menyanthes trifoliate. Also, Liang et al. (2014) reported the clear-cut inhibitory effect of high temperature on the synthesis of iridoid glycosides in Scrophularia ningpoensis. With strong support from these works, which could explain the relatively low content gentiopicroside from central Yunnan, the combination of multiple environmental factors should be considered as a vital role to influence the accumulation of iridoid glycosides in GR. Nevertheless, impacts of three climate factors on contents of sweroside and swertiamarin were not evident, although the low temperature of northwestern Yunnan caused slightly high and low levels of them, respectively. More studies are needed to unfold for compounds just like sweroside and swertiamarin using more sensitive techniques, because accumulative rules may be hidden by their small amounts.
Chemicals-Geographic Origins Association Analysis
Quality definition of GR is generally depended on levels of getiopicroside whose accumulation was apparently distinct in samples from different geographical origins. Association analysis of target ingredients across specimens of inter- and intra-tissues provides the approach to further understand the influence of environmental factors on the quality of this herbal medicine. To address this issue, Pearson's correlation analysis was performed between contents of four major ingredients (i.e., gentiopicroside, loganic acid, sweroside, and swertiamarin) belonging to the biosynthesis pathway of (seco-) iridoid glycosides in Gentianaceae species (Li et al., 2017). Pearson's correlation analysis is based on the Pearson's correlation coefficient (PCC), a range of values from 1 (a perfect positive correlation) to −1 (a perfect negative correlation) with 0 representing a random distribution (Barlow et al., 2010), which is designed for describing the linear relationship of intensity distributions between two channels in pattern recognition (Zinchuk et al., 2005).
From the Table 2, the focus of association analysis in samples of Dali (western Yunnan) was the bond between target ingredients in intra-tissues, particularly in stems. There was a significantly positive correlation (p < 0.01) between loganic acid and sweroside in both leaves and roots from Dali. Biosyntheses of sweroside, gentiopicroside, and swertiamarin were found to rely on loganic acid whose creation mainly happens in leaves, as the precursor, according to most reports about (seco-) iridoid pathway in recent years (Miettinen et al., 2014; Tham et al., 2016) However, discussing in the study associated with seco-iridoid biosynthesis in Swertia mussotii, Liu Y. et al. (2017) indicated that gene expressions associated with synthesis were extremely low in roots except for one named #SM7DLGT responsible to the coding of 7-deoxyloganetic acid glucosyl transferase who catalyzes the reaction of producing 7-deoxyloganetic acid (loganic acid's precursor). Hua et al. (2014) also hypothesized that seco-iridoids were transformed to the roots after finishing the synthesis in other tissues. Even yet, relatively few is known about which organ the biosynthesis of intermediate products take place in Gentiana species. Therefore, we speculated that the correlation between loganic acid and sweroside tends to represent the conversion relationship in intra-tissue. Even leaf's advance for producing loganic acid is likely to attribute to its synthesizing factory role known as “source,” with our data showing that correlation between loganic acid and sweroside were significant in leaves, relationships between loganic acid and other compounds could be evident in stems. Loganic acid in stems as well as gentiopicroside, revealed highly positive correlations (p < 0.05) with both swertiamarin and sweroside which simultaneously played a positive correlation (p < 0.05) with its counterpart in roots. Suggested by Inouye et al. (1970); Inouye (1971), gentiopicroside is generated from sweroside via swertiamarin in G. triflora. As a result of successively chemical conversion in one pathway, it was acceptable that the enhancement of loganic aicd's content would increase the amount of swertiamarin and sweroside, whose levels induced the gentiopicroside's content to elevate. Gentiopicroside generated obviously positive correlations (p < 0.01) with sweroside in leaves and swertiamarin in roots. The flower produced adverse results that sweroside showed a remarkably inverse correlation with swertiamarin in this tissue (p < 0.01) and gentiopicroside in roots (p < 0.01). What's more, the positive correlation (p < 0.05) was found between sweroside in stems and roots. However, while a large amount of gene related to the pathway express in leaves, there was not striking connection between gene expression profile and accumulation of (seco-) iridoid glycosides in Gentianaceae plants (Courdavault et al., 2014; Luca et al., 2014; Liu Y. et al., 2017). Thus, information obtained from different organs in present work might reflect the transporting relationship in inter-tissues.
Table 2. Pearson's correlation coefficients for contents of gentiopicroside, loganic acid, sweroside, and swertiamarin in samples collected from Dali.
The results of association analysis for samples from other two big producing areas were listed in Tables S2–S6. From the Tables S2, S3, compounds' bonds in inter-tissues accounted for around 50% of total relationships in samples from northwestern Yunnan, a low-temperature area. Connections between gentiopicroside, loganic acid, sweroside, and swertiamarin in roots and leaves were dominant in samples from Lijiang (Table S2) and Diqing (Table S3), followed by that between stems and leaves and between stems and roots. Of the broad positive correlations of target compounds in roots and leaves, the connection between loganic acid and sweroside was the most significant (p < 0.01) in both two areas' samples along with that between loganic acid and swertiamarin in samples from Diqing. It was true that ties between compounds in stems and leaves from Lijiang (p < 0.01) were solidified compared with Diqing (p < 0.05). The result of that sweroside in roots separately showed apparent links to gentiopicroside and sweroside (p < 0.01) in stems from Lijiang was agree reasonably well with that of Diqing. Differently, the core relationship of central Yunnan's samples (Tables S4–S6) was the correlation between componuds in inter-tissues, which was completely on the opposing side of western Yunnan's materials. Tissues where the compounds' bonds happened in were slightly different despite the three origins belonging to central Yunnan, a high-humidity area, finding them on the same inter-tissues relationship. The significantly positive correlation (p < 0.01) between loganic acids in stems and roots was found in samples from Qujing (Table S4), whereas that between swertiamarins in stems and leaves was observed in samples from Yuxi (Table S5). For samples from Kunming (Table S6), loganic acids in flowers correlated positively (p < 0.01) to sweroside and swertiamarin in roots, respectively. These results revealed that chemical variation in G. rigescens showed strongly detectable associated with environmental conditions impacting the accumulation of iridoid glycosides by adjusting the transformation of intermediates but also the transport of them between different tissues.
Principal Component Analysis
For a long time, GR materials collected from different places by cultivators or farmers have been sold to customers and patients in markets, whose different pharmacological activities have been proven by quiet a number of clinical experiences for centuries. The important goal of quality control of GR is to gain insight into the comprehensive information of secondary metabolites in GR, which can be accomplished by the application of FTIR that provides global messages into the overall secondary metabolites through an easy and rapid process (Surewicz et al., 1993). In this work, IR data set subjected to second derivative were loaded into R language for PCA, a non-supervision algorithm concerned with finding classes of similar objects by reducing dimensionality of the data set while retaining most of variation (Wold et al., 1987; Ringnér, 2008). In our previous paper (Wu et al., 2017), classification prediction models were established using PLS-DA, a method required the allocation of observation to sets of a priori defined classes, which generated a confusion that the classification result is due to the chemistry properties or the priori definition (Härdle et al., 2017). The artificial factor could be eliminated through the PCA method, which makes the classification more objective compared to the PLS-DA. Samples second derivative spectra of IR from each geographic origin are given in Figure 3.
Figure 3. Second derivative spectral details of root (A), stem (B), leaf (C) and flower (D) samples of G. rigescens from six geographic origins.
The classification data from the 240 samples of different tissues derived from six geographic origins are displayed as scores plots of two dominant dimensions given in Figure 4. Parameters including percentages of two dominant dimensions in total variability, cumulative percentage of variance, eigenvalues were showed in Figure 5. The Figure 4A reveals root samples from northwestern Yunnan can be separated well with the Dim 1. The Dim 2, in turn, can separate samples from central Yunnan well; but there was poor aggregation between Kunming samples and Qujing and Yuxi samples which completely overlapped. Stem samples from central Yunan were recognized very well with the Dim 2 except for Yuxi ones, according to Figure 4B. Instead, there exists the poor separation among Dali (western Yunnan), Diqing (northwestern Yunnan) and Yuxi (central Yunnan) samples. Despite long distance between leaf samples from Lijiang and Diqing, the good discriminant ability for the three large producing areas was clearly observed (Figure 4C). For flowers, samples from central Yunnan were separated well with the Dim 1, except those from Yuxi (Figure 4D). The Dim 2 allowed for good separation for Diqing and Lijiang samples and Kunming and Qujing samples, reflecting the inter-class variability in the Dim 2.
Figure 4. Scores plot for the first and second dimensions for root (A), stem (B), leaf (C) and flower (D) samples of G. rigescens from three large producing areas i.e. Central Yunnan, Western Yunnan and Northwestern Yunnan.
Figure 5. Eigenvalue, percentage of variance and cumulative percentage of variance of PCA for root (A–C), stem (D–F), leaf (G–I) and flower (J–L) samples collected from three large producing areas including six geographic origins.
From the analytical data set of IR, 52 important variables i.e., wavelength numbers were filtered and assigned by the cos2 method (Figure S3) and previous literatures, respectively. Cos2 values for these wavelength numbers are high, indicating large discriminatory ability in the classification. Unabridged lists of main peak assignments are given in Table S7. The characteristics in root samples from northwestern Yunnan were regions at 620–630, 890, 920–1,160, 1,200–1,210, 2,970, and 3,490–3,520 cm−1. Many of these were consistent with carbohydrates, phenols, benzoates, terpenoids, and xanthones (Kizil et al., 2002; Ivanova and Singh, 2003; Xu et al., 2007; Gao et al., 2010; Zhao et al., 2015; Miskolczi et al., 2003). Reporting in the paper that analyzed chemical compounds of three Gentianaceae species, Aberham et al. (2011) suggested that absorption bands at 3,413, 2,935, 1,208, 1,147, 923, and 627 cm−1 indicated the existence of xanthones. The root samples from central Yunan specially contained absorption regions at 620–630, 920–990, 2,970, and 3,490–3,520 cm−1, demonstrating existences of carbohydrates, benzoates and triterpenoids. Stem samples from central Yunnan had unique bands at 890, 920–930, 970–1,070, and 1,160 cm−1, which attributed to saccharides and aliphatic alcohols (Cael et al., 1975; Jouraiphy et al., 2008; Shen et al., 2016). The frequencies at 560–620, 670–710, 1,160–1,190, and 1,367 cm−1 were classified as important variables for the leaf samples from central Yunnan. They were indicative of phosphate, aromatic hydrocarbons, polyols, and celluloses (Hren et al., 2000; Garside and Wyeth, 2003; Chandran et al., 2006; Gerçel et al., 2007). The same important variables were observed for leaf samples from northwestern Yunnan, although no overlapping was found in samples from Lijiang and Diqing. Flower samples from central Yunnan could be distinguished by aliphatic alcohols, saccharides and terpenes (Jouraiphy et al., 2008; Qi et al., 2017b), namely, the bands at 1,010–1,110 and 1,200–1,210 cm−1. Characteristic compounds of flower samples from Diqing and Lijiang and Kunming and Qujing were large amount of saccharides, silicates, sulfates and terpenes (Zhbankov et al., 1997; Tatzber et al., 2007; Wang et al., 2011), i.e., the bands at 622, 1,016, 1,180–1,190, and 1,210 cm−1. The absorption band at 2,000–2,550 cm−1 attributed to the CO2 originated from air was removed. Elucidating these facts in conjunction with the habitat information (Table 1 and Figure S1), the observed classification results can be easily rationalized in terms of environmental conditions. The difference on chemical compounds between Diqing and Lijiang was likely influenced by geographic factors, especially the elevation.
Particularly problematic was the reorganization by scores values of samples from Kunming and Qujing, which grew with pine trees and weeds, respectively. In order to understand what makes poor aggregation between these samples, the PCA was performed using second derivatives of IR data set of samples lived in fir woods, pine woods and weeds. Figure S4 of the Supplemental files gives parameters to demonstrate the robustness of PCA. The classification data from the 472 samples stemmed from different habitats including fir woods, pine woods, and weeds were showed as scores plot of the two dominant dimensions given in Figure S5. That samples were divided into three groups along the habitats was not found even overlapped completely, regardless of tissues. Therefore, the likelihood that associated plants influenced chemical constituents in GR was extremely low, which was adverse with the result of study by Chu et al. (2016). We emphasize here that this result has been tentatively discovered through samples collected from natural setting. Single factor i.e., different associated plants experiments combined with agronomy, phytophysiology and molecular biology is needed. With strong supports from these facts, geographic factors could be credited as playing a big part in chemical constituents. On the other hand, roots, stems, leaves and flowers of GR displayed different responses to the environment.
The external standard method was employed to calculate the contents of four compounds in G. rigescens. The calibration curves of four compounds showed good linearity (R2 ≥ 0.99) in all cases. The LODs and LOQs were defined in the range 2.38–48.69 μg/mL and 8.06–158.70 μg/mL, respectively. The relative standard deviation (RSD%) for the peak area (Pa) and retention time (Rt) were below 2.90% for the intraday experiment and 3.80% for the interday experiment. The analytes demonstrated acceptable recovery efficiency (96–106%). The result was listed in Table 3.
The HPLC-DAD combined with FTIR analysis of different tissues of G. rigescens collected from six geographic origins further provided information on content of four iridoid glycosides and function group of other compounds. Gentiopicroside was the major compound in GR. The content of iridoid glycoside and chemical profile were highly depended on the climate and geography of the growth sites, resulting in the fluctuation of medicinal quality. The present work provided global information on the chemical profile and content of major iridoid glycosides in G. rigescens originated from six different origins. Temperature, water and altitude factors are worth valuing in the cultivation, a useful tool for controlling quality of herbal medicines systematically.
H-YH and Y-ZW designed the project, processed plant collection. JL processed literature review, experimental work, data collection and analysis and manuscript drafting. JZ and Y-LZ helped in statistical analysis and manuscript revision. All authors read and approved the final manuscript for publication.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors are grateful to financial support provided by the National Natural Science Foundation of China (Grant No. 81760684 and Grant No. 81660638) and the Key Project of Yunnan Province Natural Science Foundation (Grant No. 2017FA049).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2017.00125/full#supplementary-material
Figure S1. Annual mean precipitation (A), annual mean relative humidity (B) and annual mean temperature (C) data of six geographic origins of G. rigescens downloading from Climatic Data Center, National Meteorological Information Center, China Meteorological Administration.
Figure S2. The total contents of gentiopicroside (A), loganic acid (B), sweroside (C) and sweritiamarin (D) in whole plant samples from three large producing areas composed of six geographic origins (mg/g, n = 10).
Figure S3. Variable plots of PCA for the FTIR wavelength number of samples collected from three big producing areas including six geographic origins.
Figure S4. Scores plot for the first and second dimensions for the samples living in fir woods, pine woods and weeds.
Figure S5. Eigenvalue, percentage of variance and cumulative percentage of variance of PCA for samples living in fir woods, pine woods and weeds.
Table S1. The contents of chemical compositions (mg/g raw materials ± SE) in four iridoid glycosdies in samples from six geographic origins.
Table S2. Pearson's correlation coefficients for contents of gentiopicroside, loganic acid, sweroside, and swertiamarin in samples collected from Lijiang.
Table S3. Pearson's correlation coefficients for contents of gentiopicroside, loganic acid, sweroside, and swertiamarin in samples collected from Diqing.
Table S4. Pearson's correlation coefficients for contents of gentiopicroside, loganic acid, sweroside, and swertiamarin in samples collected from Qujing.
Table S5. Pearson's correlation coefficients for contents of gentiopicroside, loganic acid, sweroside, and swertiamarin in samples collected from Yuxi.
Table S6. Pearson's correlation coefficients for contents of gentiopicroside, loganic acid, sweroside, and swertiamarin in samples collected from Kuming.
Table S7. Fourier transformed infrared (FTIR) wavelength numbers, in cm-1, and their respective tentative assignments.
Aberham, A., Pieri, V., Croom, E. M. Jr., Ellmerer, E., and Stuppner, H. (2011). Analysis of iridoids, secoiridoids and xanthones in Centaurium erythraea, Frasera caroliniensis and Gentiana lutea using LC-MA and RP-HPLC. J. Pharmaceut. Biomed. 54, 517–525. doi: 10.1016/j.jpba.2010.09.030
Balkan, I. A., Akülke, A. Z., Bagatur, Y., Telci, D., Gören, A. C., Kirmizibekmez, H., et al. (2017). Sambulin A and B, non-glycosidic iridoids from Sambucus ebulus, exert significant in vitro anti-inflammatory activity in LPS-induced RAW 264.7 macrophages via inhibition of MAPKs's phosphorylation. J. Ethnopharmacol. 206, 347–352. doi: 10.1016/j.jep.2017.06.002
Barlow, A. L., MacLeod, A., Noppen, S., Sanderson, J., and Guérin, C. J. (2010). Colocalization analysis in fluorescence micrographs: verification of a more accurate calculation of pearson's correlation coefficient. Microsc. Microanal. 16, 710–724. doi: 10.1017/S143192761009389X
Cael, J. J., Koenig, J. L., and Blackwell, J. (1975). Infrared and raman spectroscopy of carbohydrates. Part VI: normal coordinate analysis of V-amylose. Biopolymers 14, 1885–1903. doi: 10.1002/bip.1975.360140909
Chandran, S. P., Chaudhary, M., Pasricha, R., Ahmad, A., and Sastry, M. (2006). Synthesis of gold nanotriangles and silver nanoparticles using Aloe vera plant extract. Biotechnol. Prog. 22, 577–583. doi: 10.1021/bp0501423
Chu, B. W., Zhang, J., Li, Z. M., Zhao, Y., Zuo, Z., Wang, Y., et al. (2016). Evaluation and quantitative analysis of different growth periods of herb-arbor intercropping systems using HPLC and UV-vis methods coupled with chemometrics. J. Nat. Med. 70, 803–810. doi: 10.1007/s11418-016-1009-x
Courdavault, V., Papon, N., Clastre, M., Giglioli-Guivarc'h, N., St-Pierre, B., and Burlat, V. (2014). A look inside an alkaloid multisite plant: the Catharanthus logistics. Curr. Opin. Plant. Biol. 19, 43–50. doi: 10.1016/j.pbi.2014.03.010
Ding, G., Wang, Y., Liu, A., Hou, Y., Zhang, T., Bai, G., et al. (2017). From chemical markers to quality markers: an integrated approach of UPLC/Q-TOF, NIRS, and chemometrics for the quality assessment of honeysuckle buds. RSC. Adv. 7, 22034. doi: 10.1039/C6RA28152D
Gao, L., Li, J., and Qi, J. (2010). Gentisides A and B, two new neuritogenic compounds from the traditional Chinese medicines Gentiana rigescens Franch. Bioorg. Med. Chem. 18, 2131–2134. doi: 10.1016/j.bmc.2010.02.004
Gerçel, Ö., Özcan, A., Özcan, A. S., and Gerçel, H. F. (2007). Preparation of activated carbon form a renewable bio-plant of Euphorbia rigida by H2SO4 activation and its adsorption behavior in aqueous solutions. Appl. Surf. Sci. 253, 4843–4852. doi: 10.1016/j.apsusc.2006.10.053
Guo, Z. Y., Li, P., Huang, W., Wang, J. J., Liu, Y., Liu, B., et al. (2014). Antioxidant and anti-inflammatory caffeoyl phenylpropanoid and secoiridoid glycosides from Jasminum nervosum stems, a Chinese folk medicine. Phytochemistry 106, 124–133. doi: 10.1016/j.phytochem.2014.07.011
Guzelmeric, E., Ristivojević, P., Vovk, I., Milojković-Opsenica, D., and Yesilada, E. (2017). Quality assessment of marketed chamomile tea products by a validated HPTLC method combined with multivariate analysis. J. Pharmaceut. Biome. 132, 35–45. doi: 10.1016/j.jpba.2016.09.030
He, Y., Zhu, S., Ge, Y., Kazuma, K., Zou, K., Cai, S., et al. (2015). The anti-inflammatory secoiridoid glycosides from Gentianae Scabrae Radix: the root and rhizome of Gentiana scabra. J. Nat. Med. 69, 303–312. doi: 10.1007/s11418-015-0894-8
Hirri, A., Bassbasi, M., Platikanov, S., Tauler, R., and Oussama, A. (2016). FTIR spectroscopy and PLS-DA classification and prediction of four commercial grade virgin olive oils from Morocco. Food. Anal. Methods. 9, 974–981. doi: 10.1007/s12161-015-0255-y
Hua, W., Zheng, P., He, Y., Cui, L., Kong, W., and Wang, Z. (2014). An insight into the genes involved in secoiridoid biosynthesis in Gentiana macrophylla by RNA-seq. Mol. Biol. Rep. 7, 4817–4825. doi: 10.1007/s11033-014-3352-x
Jamieson, M. A., Quintero, C., and Blumenthal, D. M. (2013). Interactive effects of simulated nitrogen deposition and altered precipitation patterns on plant allelochemical concentrations. J. Chem. Ecol. 39, 1204–1208. doi: 10.1007/s10886-013-0340-x
Jouraiphy, A., Amir, S., Winterton, P., El Gharous, M., Revel, J. C., and Hafidi, M. (2008). Structural study of the fulvic fraction during composting of activated sludge-plant matter: element analysis, FTIR and 13C NMR. Bioresour. Technol. 99, 1066–1072. doi: 10.1016/j.biortech.2007.02.031
Kumarasamy, Y., Nahar, L., Cox, P. J., Haspars, M., and Sarker, S. D. (2003). Bioactivity of secoiridoid glycosides from Centaurium erythraea. Phytomedicine 10, 344–347. doi: 10.1078/094471103322004857
Li, J., Zhao, Y., Huang, H., and Wang, Y. (2017). Phytochemistry and pharmacological activities of the genus Swertia (Gentianaceae): a review. Am. J. Chin. Med. 45: 667–736. doi: 10.1142/S0192415X17500380
Liang, Z., An, Y., and Liu, H. (2014). High temperature stress decreases root iridoid glycosides biosynthesis of Scrophularia ningpoensis during florescence. J. Med. Plants Res. 8, 392–394. doi: 10.5897/JMPR12.336
Liu, C., Chen, S., Xiao, X., Zhang, T., Hou, W., and Liao, M. (2016). A new concept on quality marker of Chinese material medica: quality control for Chinese medicine products. Chin. Tradint. Herbal Drugs 47, 1443–1457. doi: 10.1016/S1674-6384(17)60070-4
Liu, J., Liu, Y., Wang, Y., Abozeid, A., Zu, Y., and Tang, Z. (2017). The integration of GC-MS and LC-MS to assay the metabolomics profiling in Panax ginseng and Panax quinquefolius reveals a tissue- and species-specific connectivity of primary metabolites and ginsenosides accumulation. J. Pharmaceut. Biomed. 135, 176–185. doi: 10.1016/j.jpba.2016.12.026
Liu, Y., Wang, Y., Guo, F., Zhan, L., Mohr, T., Cheng, P., et al. (2017). Deep sequencing and transcriptome analyses to identity genes involved in secoiridoid biosynthesis in the Tibetan medicinal plant Swertia mussotii. Sci. Rep. 7:43108. doi: 10.1038/srep43108
Luca, V., Salim, V., Thamm, A., Masada, S., and Yu, F. (2014). Making iridoids/secoiridoids and monoterpenoid indole alkaloids: progress on pathway elucidation. Curr. Opin. Plant. Biol. 19, 35–42. doi: 10.1016/j.pbi.2014.03.006
Ma, S., He, Z., Deng, X., But, P. P., Ooi, V. E., Xu, H., et al. (2001). In vitro evaluation of secoiridoid glucosides from the fruits of Ligustrum lucidum as antiviral agents. Chem. Pharm. Bull. 49, 1471–1473. doi: 10.1248/cpb.49.1471
Martz, F., Turunen, M., Julkunen-Tiitto, R., Lakkala, K., and Sutinen, M. (2009). Effect of the temperature and the exclusion of UVB radiation on the phenolics and iridoids in Menyanthes trifoliate L. leaves in the subarctic. Environ. Pollut. 157, 3471–3478. doi: 10.1016/j.envpol.2009.06.022
Melito, S., Petretto, G. L., Podani, J., Foddai, M., Maldini, M., Chessa, M., et al. (2016). Altitude and climate influence Helichrysum italicum subsp. microphyllum essential oils composition. Ind. Crop. Prod. 80, 242–250. doi: 10.1016/j.indcrop.2015.11.014
Mihailović, V., Katanić, J., Mišić, D., Stanković, V., Mihailović, M., Uskoković, A., et al. (2014). Hepatoprotective effects of secoiridoid-rich extracts from Gentiana cruciata L. against carbon tetrachloride induced liver damage in rats. Food. Funct. 2, 1795–1803. doi: 10.1039/C4FO00088A
Nguyen, H. T., Lee, D., Lee, W. J., Lee, G., Yoon, S. J., Shin, B., et al. (2016). UPLC-QTOFMS based metabolomics followed by stepwise partial least square-discriminant analysis (PLS-DA) explore the possible relation between the variations in secondary metabolites and the phylogenetic divergences of the genus Panax. J. Chromatogr. B 1012, 61–68. doi: 10.1016/j.jchromb.2016.01.002
Ordoudi, S. A., Pascual, M. M., and Tsimidou, M. Z. (2014). On the quality control of traded saffron by means of transmission Fourier-transform mid-infrared (FT-MIR) spectroscopy and chemometrics. Food. Chem. 150, 414–421. doi: 10.1016/j.foodchem.2013.11.014
Pan, Y., Shen, T., Zhang, J., Zhao, Y., Wang, Y., and Li, W. (2015). Simultaneous determination of six index constituents and comparative analysis of four ethnomedicines from genus Gentiana using a UPLC-UV-MS method. Biomed. Chromatogr. 29, 87–96. doi: 10.1002/bmc.3243
Pan, Y., Zhao, Y., Zhang, J., Li, W., and Wang, Y. (2016). Phytochemistry and pharmacological activities of the genus Gentiana (Gentianaceae). Chem. Biodivers. 13, 107–150. doi: 10.1002/cbdv.201500333
Qi, L., Zhang, J., Zhao, Y., Zuo, Z., Jin, H., and Wang, Y. (2017a). Quantitative and qualitative characterization of Gentiana rigescens Franch (Gentianaceae) on different parts and cultivation years by HPLC anf FTIR spectroscopy. J. Anal. Methods Chem. 2017. doi: 10.1155/2017/3194146
Qi, L., Zhang, J., Zuo, Z., Zhao, Y., Wang, Y., and Hang, J. (2017b). Determination of iridoids in Gentiana rigescens by infrared spectroscopy and multivariate analysis. Anal. Lett. 50, 389–401. doi: 10.1080/00032719.2016.1178757
Shen, Y. X., Zhao, Y. L., Zhang, J., Zuo, Z. T., Wang, Y. Z., and Zhang, Q. Z. (2016). Study on different parts of wild and cultivated Gentiana rigescens with fourier transform infrared spectroscopy. Spectrosc. Spect. Anal. 36, 667–671.
Song, Z., Hashi, Y., Sun, H., Liang, Y., Lan, Y., Wang, H., et al. (2013). Simultaneous determination of 19 flavonoids in commercial trollflowers by using high-performance liquid chromatography and classification of samples by hierarchical clustering analysis. Fitoterapia 91, 272–279. doi: 10.1016/j.fitote.2013.09.006
Surewicz, W. K., Mantsch, H. H., and Chapman, D. (1993). Determination of protein secondary structure by fourier transform infrared spectroscopy: a critical assessment. Biochemistry 32, 389–394. doi: 10.1021/bi00053a001
Suyama, Y., Kurimoto, S., Kawazoe, K., Murakami, K., Sun, H., Li, S., et al. (2013). Rigenolide A, a new secoiridoid glucoside with a cyclobutane skeleton, and three new acylated secoiridoid glucosides from Gentiana rigescens Franch. Fitoterapia 91, 166–172. doi: 10.1016/j.fitote.2013.08.006
Suyama, Y., Tanaka, N., Kawazoe, K., Murakami, K., Li, S., Sun, H., et al. (2017). Rigenolides B and C, conjugates of norsecoiridoid and secoiridoid glucoside from Gentiana rigescens Franch. Tetrahedron. Lett. 58, 1459–1461. doi: 10.1016/j.tetlet.2017.02.075
Tang, X., Yang, Q., Gong, J., Han, H., Yang, L., and Wang, Z. (2016). Target profiling analysis of bile acids in the evaluation of hepatoprotective effect of gentiopicroside on ANIT-induced cholestatic liver injury in mice. J. Ethnopharmacol. 194, 63–71. doi: 10.1016/j.jep.2016.08.049
Tatzber, M., Stemmer, M., Spiegel, H., Katzlberger, C., Haberhauer, G., Mentler, A., et al. (2007). FTIR-spectroscopic characterization of humic acids and humin fractions obtained by advanced NaOH, Na4P2O7 and Na2CO3 extraction procedures. J. Plant. Nutr. Soil. Sci. 170, 522–529. doi: 10.1002/jpln.200622082
Tham, A. M. K., Qu, Y., and Luca, V. (2016). Discovery and metabolic engineering of iridoid/secoiridoid and monoterpenoid indole alkaloid biosynthesis. Phytochem. Rev. 15, 339–361. doi: 10.1007/s11101-016-9468-y
Vuorinen, A. L., Markkinen, N., Kalpio, M., Linderborg, K. M., Yang, B., and Kallio, H. P. (2015). Effect of growth environment on the gene expression and lipid related to triacylglycerol biosynthesis in sea buckthorn (Hippophaë rhamnoides) berries. Food. Res. Int. 77, 608–619. doi: 10.1016/j.foodres.2015.08.023
Wang, S., Xu, Y., Jiang, W., and Zhang, Y. (2013). Isolation and identification of constituents with activity of inhibiting nitric oxide production in raw 264.7 macrophages from Gentiana triflora. Planta. Med. 79, 680–686. doi: 10.1055/s-0032-1328460
Wang, Y., Li, W., Wu, D., Wang, J., Wu, C., Liao, J., et al. (2011). In vitro activity of 2-methoxy-1,4-naphthoquinone and stimasta-7,22-diene-3β-ol from Impatiens balsamina L. against multiple antibiotic-resistant Helicobacter pylori. Evid. Based. Compl. Alt. 2011. doi: 10.1093/ecam/nep147
Wu, Z., Zhao, Y., Zhang, J., and Wang, Y. (2017). Quality assessment of Gentiana rigescens from different geographical origins using FT-IR spectroscopy combined with HPLC. Molecules 22:1238. doi: 10.3390/molecules22071238
Zanatta, T. S. C., Manica-Berto, R., Ferreira, C. D., Cardozo, M. M. C., Rombaldi, C. V., Zambiazi, R. C., et al. (2017). Phosphate fertilizer and growing environment changes the phytochemicals, oil quality, and nutritional composition of roundup ready genetically modified and conventional soybean. J. Agric. Food. Chem. 65, 2661–2669. doi: 10.1021/acs.jafc.6b05499
Zhang, J., Shen, T., Yang, W., Zhang, J., and Wang, Y. (2012). Investigation and evaluation of Gentiana rigescens resource in Yunnan Province. J. Plant. Genet. Resour. 13, 890–895. doi: 10.13430/j.cnki.jpgr.2012.05.027
Zhang, L., Liu, Y., Liu, Z., Wang, C., Song, Z., Liu, Y., et al. (2016). Comparison of the roots of Salvia miltiorrhiza Bunge (Danshen) and its variety S. miltiorrhiza Bgef. Alba (Baihua Danshen) based on multiwavelength HPLC-fingerprinting and contents of nine active components. Anal. Methods 8, 3171–3182. doi: 10.1039/C6AY00185H
Zhao, Y., Zhang, J., Jin, H., and Zhan, J. (2015). Discrimination of Gentiana rigescens from different origins by fourier transform infrared spectroscopy combined with chemometric methods. J. AOAC. Int. 98, 22–26. doi: 10.5740/jaoacint.13-395
Zhbankov, R. G., Andrianov, V. M., and Marchewka, M. K. (1997). Fourier transform IR and raman spectroscopy and structure of carbohydrates. J. Mol. Struct. 436, 637–654. doi: 10.1016/S0022-2860(97)00141-5
Zinchuk, V., Zinchuk, O., and Okada, T. (2005). Experimental LPS-induced cholestasis alters subcellular distribution and affects colocalization of Mrp2 and Besp protein: a quantitative colocalization study. Microsc. Res. Tech. 67, 65–70. doi: 10.1002/jemt.20184
Keywords: Gentiana rigescens, HPLC, FTIR, principal component analysis, quality assessment
Citation: Li J, Zhang J, Zhao Y-L, Huang H-Y and Wang Y-Z (2017) Comprehensive Quality Assessment Based Specific Chemical Profiles for Geographic and Tissue Variation in Gentiana rigescens Using HPLC and FTIR Method Combined with Principal Component Analysis. Front. Chem. 5:125. doi: 10.3389/fchem.2017.00125
Received: 21 October 2017; Accepted: 12 December 2017;
Published: 22 December 2017.
Edited by:Nicole J. Jaffrezic-Renault, Claude Bernard University Lyon 1, France
Reviewed by:Xia Guan, Louisiana State University, United States
Pawel Pohl, Wrocław University of Technology, Poland
Copyright © 2017 Li, Zhang, Zhao, Huang and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.