Quantitative comparison and evaluation between aerial and underground parts of Gentiana straminea through simultaneous determination of five major compounds by RP-HPLC

TL;DR

This study compares the medicinal value of above and below ground parts of Gentiana straminea to promote sustainable harvesting.

Contribution

It introduces a method to evaluate and substitute aerial parts for underground roots in herbal medicine.

Findings

Aerial parts of G. straminea have higher composite scores than underground parts in different regions.

Aerial parts can partially replace underground roots for medicinal use without harming the plant.

Ideal growing regions for aerial parts were identified on the Qinghai-Tibetan Plateau.

Abstract

Herbal Gentiana straminea Maxim. (Family Gentianaceae), “Ma Hua Jiao” in Chinese, is a commonly used Chinese medicine. Secoiridoids and flavonoids have been identified as the major active components of herbal medicines used in the treatment of hepatitis, rheumatism and many other diseases. It is the overharvesting of the roots of this plant for medicinal purposes that has led to a drastic decline in its population. In the present study, the above and below ground parts of Gentian Bitter Glycine were quantitatively compared and evaluated for the determination of the major active constituents. Five major compounds, loganic acid, swertiamarin, gentiopicroside, sweorside and isoorientin, were extracted by solvent extraction technique and analyzed by Reversed-phase High Performance Liquid Chromatography (RP-HPLC). By analysing the principal components and calculating the composite scores, the results show that the aboveground component in different areas ranked higher compared to the underground component, with the former being able to substitute to some extent for the latter's underground component. Finally, based on hierarchical cluster analysis, we identified the ideal natural growing region for aerial parts of G. straminea distributed on the Qinghai-Tibetan Plateau. The significance of this work is that we can balance the demand for herbs with environmental preservation by selectively picking the aerial parts, which can regrow next year, instead of removing the whole plant. It protects the fragile ecological environment of the Tibetan Plateau and is important for sustainable development.

Figures

Structures of investigated five compounds and its UV chromatograms.Fig. 1

HPLC chromatograms of reference and samples. (A) Mixed standards; (B) Aerial parts of *G. straminea*; (C) Underground parts of *G. straminea*. 1, Loganic acid; 2, Swertiamarin; 3, Gentiopicroside; 4, Sweorside; 5, Isoorientin.Fig. 2

Correlation matrix.Fig. 3Table 6Variance explained and extraction sums of squared loadings.Table 6Component of Aerial partsComponent of Underground partsComponent of Aerial parts plus Underground parts1 2 31 2 31 2 3Eigenvalues2.9710.8210.7122.9441.0620.7812.4541.4710.722% of Variance59.43016.42414.23758.87221.24215.63049.08329.41714.445Cumulative %59.43075.85490.09058.87280.11495.74349.08378.50092.945Extraction Method: Principal Component Analysis.Fig. 4Component matrix.Fig. 4Table 7Component equation.Table 7SamplesComponentComponent equationAerial parts*F*_*1*_*F*_*1*_ = *0.487x*_*1*_*+0.524

Stacked histogram of five major pharmacological constituents of the above ground portion of *G. straminea*.Fig. 6

Variations in the mean content of the five active constituents (a) Loganic acid, (b) Swertiamarin, (c) Gentiopicroside, (d) Swerosideand and (e) Isoorientin, over the altitude gradient.Fig. 7

Full text

Introduction

1

Gentiana straminea Maxim. [1,2], a very famous Chinese folk herb of the family Gentianaceae, commonly known as “Ma Hua Jiao,” is quite frequently used to treat diverse diseases, such as fungal and bacterial infections, jaundice, hepatitis, constipation, pain and rheumatism [3,4]. This medicine contains several bitter compounds, mainly secoiridoids and flavonoids, which are regarded as the major active components [[5], [6], [7], [8]]. Flavonoids have been reported to have antioxidant activity and antimicrobial potential [[9], [10], [11]], while its ethanol extract produces an inhibitory effect on anthrax bacillus, staphylococci, typhoid bacillus, pneumococci, bacillus dysenteriae, and vibrio [12,13]. Historically, Chinese physicians prepared this medicine using the root of G. straminea, while Tibetan doctors typically used the entire plant. Following our previous phytochemical investigation of G. straminea [[14], [15], [16], [17], [18]], we determined that its major phytochemical constituents are loganic acid, swertiamarin, gentiopicroside, sweorside, and isoorientin. These active constituents, especially gentiopicroside, a major bitter secoiridoid glycoside widespread in the Gentianaceae family, were found to be capable of protecting against hepatitis by inhibiting the production of tumor necrosis factor [19,20], and by promoting the secretion of gastric juice [21]. Another important biologically active compound is swertiamarin, which was observed to improve skin function, promote hair growth, and function as an antipsychotic [21,22]. Loganic acid was observed to exhibit analgesic, anti-inflammatory, and febrifuge pharmacological effects [23,24]. Both sweorside and isoorientin exhibit the ability to suppress chemically and immunologically induced hepatic injuries [[25], [26], [27]].

This species is a perennial herb distributed in the high mountains and alpine environment of the Qinghai-Tibet Plateau, at altitudes from 2, 000 to 5, 000 m [28]. The alpine plants growing in these regions are often subject to serious stresses to which they have had to adapt, including low temperatures and strong solar radiation [29]. In recent years, the natural resources of G. straminea have declined dramatically due to overharvesting of the flowering plants, without leaving sufficient seeds to maintain its population. This species has now been listed as endangered by the local governments, and further harvesting has been prohibited in some areas. In order to devise adequate conservation and management strategies for this species, it is important to seek new technology that will allow us to balance the continued utilization of this resource with the need to protect the fragile ecological environment of Qinghai-Tibet Plateau.

The medicinal importance and lack of resources of the plant prompted our interest in finding alternatives to one-time destructive harvesting; specifically, the selective harvesting of aerial parts of G. straminea while retaining its underground parts for continued perennial growth and maximum resource utilization. A number of published methods have focused on the qualitative and quantitative determination of cyclic enol ether terpenoids and flavonoids in G. straminea using Liquid Chromatography (LC) or Liquid Chromatography-Mass Spectrometry (LC-MS) [30,31]. However, few of these methods aim at quantitatively comparing and evaluating the above-ground and below-ground parts of the plant by determining the major components. The purpose of this paper was to statistically compare the concentrations of the five major active constituents in the aerial and underground parts of G. straminea using materials collected from fifteen different high-altitude populations on the Qinghai-Tibetan Plateau. Principal component analysis in SPSS software has been previously used as a method to analyze and explain data prior to HPLC analysis. Results of this investigation indicated that the aerial parts yielded considerable concentrations of the five active principles that are also found in the underground parts, which are currently harvested for the treatment of diverse ailments. Furthermore, hierarchical clustering analysis [32] also helped to identify the most suitable regions for producing the aerial parts of G. straminea.

Methods

2

Chemicals, reagents and materials

2.1

Some findings have shown that polar organic solvents are the most efficient in the extraction of flavonoid and phenolic components [33]. In our laboratory, methanol was used for the extraction of the main components and the isolation and purification of the marker components, loganic acid, swertiamarin, gentiopicroside, sweorside, and isoorientin, were carried out. The purity of each compound was determined by HPLC to be greater than 98 %. Compound structures were confirmed by comparing their UV, MS, ^1^HNMR and ^13^CNMR data with data from the literature [[34], [35], [36], [37], [38]]. Loganic acid was generously provided by the Northwest Institute for Plateau Biology, Chinese Academy of Sciences (Xining, China). The chemical structures of these reference compounds are summarized in Fig. 1.Fig. 1. Structures of investigated five compounds and its UV chromatograms.Fig. 1

HPLC-grade methanol was purchased from Yuwang Co., Ltd (Shandong, China). The other chemicals and reagents used in the experiment were of Chromatographic or analytical grade quality.

Natural G. straminea was obtained from different counties of Qinghai province of China. Cultured G. straminea was obtained from Huangnan prefecture of Qinghai province of China (Table 1). The identities of these G. straminea samples were confirmed by a member of our research group, Professor Pengcheng Lin. All corresponding voucher specimens were deposited at the experimental center of the college of Chemistry and Life Science, Qinghai University for Nationalities, Xining, China. Dried samples were ground to 0.30–0.45 mm.Table 1. Summary for the tested samples of G. straminea.Table 1NO.CodeSamplesSoursesAltitudeCollection Date1MYG. stramineaMenyuan county, qinghai, China3100 m2008.82ZK-1G. stramineaZeku county1, qinghai, China3350 m2008.83HZG. stramineaHuzhu county, qinghai,China2200 m2008.84NQG. stramineaNangqian county, qinghai,China3700 m2008.85ZK-2G. stramineaZeku county2, qinghai, China3350 m2008.86TDG. stramineaTongde county, qinghai, China3400 m2008.87HN-1G. stramineaHenan county1, qinghai, China3200 m2008.88YS-1G. stramineaYushu county1, qinghai, China3500 m2008.89HN-2G. stramineaHenan county2, qinghai, China3200 m2008.810QL-1G. stramineaQilian county1, qinghai, China3300 m2008.811TR-1G. stramineaTongren county1, qinghai, China3200 m2008.812HNPaG. stramineaaHuangnanPrefecture, qinghai, China3150 m2008.813QL-2G. stramineaQilian county2, qinghai, China3300 m2008.814TR-2G. stramineaTongren county2, qinghai, China3200 m2008.815YS-2G. stramineaYushu county2, qinghai, China3500 m2008.8aSamples of G. straminea were cultured, while other species were wild.

Sample preparation

2.2

The Analytes were first dissolved in methanol solution at approximately 5 mg/mL as a stock solution. Aliquots of the stock solution were transferred to a 2-mL volumetric flask and methanol added to obtain the desired concentration.

An accurately weighed sample (0.50 g) of ground powder (0.30–0.45 mm) was transferred into a 25 mL round-bottomed flask with stopper, and 10 mL methanol was added. The flask was immersed in an ultrasonic water-bath (Kunshan Instrumental, Kunshan, Zhejiang Province, P.R.China) and sonicated for 30 min at 30 °C. The sample was extracted twice and the extracts combined. Finally, the solution was filtered into a 25 mL volumetric flask, followed by addition of methanol to a specified volume, and filtration through a 0.45 μm filter (Agilent Technologies, USA), prior to HPLC analysis.

HPLC analysis

2.3

Analysis was performed on an Agilent Series 1100 liquid chromatograph equipped with a vacuum degasser, a quaternary pump, and a diode-array detector (DAD), connected to an Agilent ChemStation running ChemStation software. A Hypersil ODS column (250 mm × 4.6 mm, 5 μm) was used. The standards and samples were separated using a gradient mobile phase consisting of water with 3 mM phosphate (A) and methanol (B). The elution conditions applied were as follows: 0–10 min, isocratic 20 % B; 10–30 min, linear gradient 20–50 % B; and 30–35 min, 50 % B isocratic. The column was then reconditioned with 20 % B isocratic for 15 min after washing for 10 min with 100 % B. The flow rate was 1 mL/min and the injection volume was 10 μL. The column temperature was maintained at 35 °C. The analytes were monitored at 254 nm.

Data analysis

2.4

Principal Components Analysis (PCA) is a useful statistical technique and a way to identify patterns in data by expressing the data in such a way as to highlight their similarities and differences. Since patterns in data can be difficult to identify in data of high dimensions, where the luxury of graphical representation is not available, PCA is a powerful alternative tool for data analysis. The other main advantage of PCA is that once you have identified patterns in the data, data can be compressed by reducing the number of dimensions, without significant loss of information. For our analyses, PCA was performed on sets of data using SPSS 15 software for Windows (SPSS Inc., USA).

Hierarchical clustering analysis (HCA) was also performed using SPSS 15 for windows (SPSS Inc., USA). HCA is comprised of a number of “procedures” – graphical, statistical, reporting, processing and tabulating procedures – that enable simple and rapid data evaluation. Between-groups Linkage, a very efficient method for the analysis of variance between clusters, was applied, and the Squared Euclidean distance was selected as a measurement for hierarchical clustering analysis.

Results

3

Calibration curves and recovery

3.1

For every calibration curve, the calibration concentrations were back-calculated from the peak area of the analytes (Table 2). Deviation from the nominal concentration defined the accuracy (Table 3).Table 2. Linear regression data, LOD and LOQ of investigated compounds from G. straminea.Table 2. AnalytesLinear regression dataLOD (ng)LOQ (ng)Regressive equationLinear range(mg/mL)R^2^Loganic acidy = 173.4x-1220.364–2.180.99991.203.0Swertiamariny = 195.6x-54.00.275–1.650.99990.581.8Gentiopicrosidey = 1002x+34.40.614–3.680.99990.501.3Sweorsidey = 5744x-1950.0656–0.3940.99990.842.1Isoorientiny = 1730x-1930.0899–0.5390.99990.481.5R^2^= Squares of correlation coefficients for the standard curves.Table 3. Recoveries for the assay of the five compounds in G. straminea.Table 3. AnalytesOriginal(mg)Spiked(mg)Found(mg)Recovery (%)RSD (%)Loganic acid0.8880.9101.820102.40.50Swertiamarin0.6270.6881.320100.71.10Gentiopicroside1.5601.5353.100100.70.55Sweorside0.1450.1640.310100.60.33Isoorientin0.2160.2250.44099.60.49a, Data calculated from the average of five determinations.Recovery (%) = 100*(amount found-original amount)/amount spiked.RSD (%) = 100*SD/mean.

A recovery test was used to evaluate the accuracy of the developed method. Known amount of the five analytes were added to approximately 0.50 g of G. straminea powder, and then extracted and analyzed as described above. Each sample was analyzed five times. The recovery was calculated based on the total amount of the individual, investigated components, and the average recoveries were determined by the following formula: recovery (%) = first amount found/total amount of individual investigated components × 100 %.

Limit of detection and quantification

3.2

The LOD was defined as the concentration of analyte giving S/N = 3, and the LOQ as that giving S/N = 10. The LOD and LOQ values of the investigated compounds are summarized in Table 2.

Precision, repeatability and stability

3.3

The variability test was chosen to determine the precision of the developed assay. The intra-day precision was examined on the mixed standards solutions for six replicates within a single day, and the relative standard deviation (RSD) taken as a measure of precision was <1 % of the expected concentration. The repeatability of the developed method for all investigated compounds was determined based on six trials. The repeatability was presented as RSD (n = 6). Stability of sample solution was tested, by analysis of the sample solution every 2 h for 12 h. Variation was expressed as RSD (n = 3). These results are summarized in Table 4.Table 4. Precision, repeatability and stability.Table 4. AnalytesPrecision (RSD,%,n = 6)Repeatability (RSD,%,n = 6)Stability (RSD,%,n = 3)Loganic acid0.251.220.38Swertiamarin0.891.630.88Gentiopicroside0.311.580.45Sweorside0.681.880.82Isoorientin0.551.450.47

Quantitation of the five components of G. straminea from 15 different regions

3.4

Natural and cultured G. straminea were analyzed using the calibration curve of each investigated compound (Table 5). This analysis identified different distributions of the five active components when comparing the aerial parts (Fig. 2B), underground parts (Fig. 2C) and mixed standards (Fig. 2A).Table 5. Content (%) of the five components of G. straminea from 15 different regions.Table 5. SamplesLoganic acidSwertiamarinGentiopicrosideSweorsideIsoorientinSumTatolMYAerial parts0.375a0.5093.0630.0440.0894.0809.182Underground parts___b0.7154.2210.0730.0935.102ZK-1Aerial parts0.6332.05312.0480.2150.34715.29622.225Underground parts___b0.9615.7540.0800.1346.929HZAerial parts0.5501.5018.6000.1500.27011.07121.376Underground parts0.3891.3738.2600.1300.15310.305NQAerial parts3.5631.6898.5800.1440.60314.57925.185Underground parts2.9411.2226.1550.0850.20310.606ZK-2Aerial parts1.9561.3205.1350.2440.4629.11715.794Underground parts2.3610.8823.2550.0400.1396.677TDAerial parts3.5232.2488.0020.1220.52714.42222.074Underground parts1.8231.1134.5240.0820.1107.652HN-1Aerial parts3.8973.0236.3700.0620.76614.11825.525Underground parts4.6911.7064.8230.0530.13411.407YS-1Aerial parts2.3712.4634.7950.0770.59210.29829.406Underground parts5.8631.75111.3500.0290.11519.108HN-2Aerial parts2.3323.1217.5990.0830.90314.03821.715Underground parts1.2291.4624.7750.0710.1407.677QL-1Aerial parts2.6613.6168.5920.0710.50415.44424.516Underground parts2.3541.5664.9390.0420.1719.072TR-1Aerial parts4.4093.7698.7810.0750.95217.62630.541Underground parts3.9921.7566.9820.0400.14612.916HNPAerial parts4.2992.95914.4090.0370.55022.25451.174Underground parts10.1962.75015.7830.0450.14728.921QL-2Aerial parts4.3712.2458.4130.1690.36615.56431.398Underground parts5.1541.5159.0300.0420.09315.834TR-2Aerial parts3.7343.3287.7070.0590.80715.63524.716Underground parts2.6161.5434.5850.0580.2799.0805YS-2Aerial parts3.3783.3529.3920.1360.67116.92932.942Underground parts4.2851.9449.6400.0280.11616.013aThe date was present as average of three replicates (RSD<2 %).bUndetected.Fig. 2HPLC chromatograms of reference and samples. (A) Mixed standards; (B) Aerial parts of G. straminea; (C) Underground parts of G. straminea. 1, Loganic acid; 2, Swertiamarin; 3, Gentiopicroside; 4, Sweorside; 5, Isoorientin.Fig. 2

Comparison of aerial parts and underground parts of G. straminea

3.5

Procedures of principal component analysis involve the following: 1) Standardization of raw data; 2) Calculation of covariance and correlation coefficient matrix (Fig. 3); 3) Calculation of the eigenvalues and eigenvectors of the covariance matrix (Table 6); and 4) Determination of the principal components combined with professional knowledge to explain the information contained and given by the principal component analysis (Fig. 4, Table 7, Table 8, Table 9) (see Table 10).Fig. 3. Correlation matrix.Fig. 3. Table 6Variance explained and extraction sums of squared loadings.Table 6. Component of Aerial partsComponent of Underground partsComponent of Aerial parts plus Underground parts1 2 31 2 31 2 3Eigenvalues2.9710.8210.7122.9441.0620.7812.4541.4710.722% of Variance59.43016.42414.23758.87221.24215.63049.08329.41714.445Cumulative %59.43075.85490.09058.87280.11495.74349.08378.50092.945Extraction Method: Principal Component Analysis.Fig. 4. Component matrix.Fig. 4. Table 7Component equation.Table 7. SamplesComponentComponent equationAerial partsF1F1 = 0.487x1+0.524 x2+0.403 x3+0.309x4+0.481 x5F2F2 = 0.281x1-0.032 x2+0.430 x3+0.854x4-0.064 x5F3F3 = -0.038x1+0.236 x2-0.660 x3+0.398x4+0.590 x5Underground partsF1F1 = 0562x1+0.539 x2+0.507 x3-0.364x4-0.061 x5F2F2 = 0.045x1+0.238 x2-0.038 x3+0.210x4+0.947 x5F3F3 = -0.043x1+0.131x2+0.488 x3+0.840x4-0.197 x5Aerial plus Underground partsF1F1 = 0.475x1+0.574 x2+0.519 x3-0.093x4+0.408 x5F2F2 = -0.427x1+0.256 x2-0.222 x3+0.623x4+0.561x5F3F3 = 0.232x1-0.244 x2+0.465x3+0.740x4-0.352 x5Table 8Component score and composite score ranking of Aerial parts and Underground parts.Table 8NO.CodeAerial partsUnderground partsF_1_F_2_F_3_FarankingF_1_F_2_F_3_Fbranking1MY−3.338−2.061−0.753−2.42815−2.021−1.3380.078−1.462152ZK-1−2.2240.7180.551−1.12513−1.671−0.3540.408−0.996143HZ−1.9790.301−0.8715−1.25014−1.4580.5912.334−0.369104NQ−0.1090.734−0.1140.03910−0.8661.2510.351−0.19095ZK-2−2.2231.2181.044−0.97212−1.109−0.531−1.267−0.963136TD0.0500.321−0.1500.0619−1.303−0.7070.403−0.855127HN-11.122−0.7880.6850.63450.343−0.116−0.4390.10968YS-1−0.085−0.7550.137−0.154111.966−0.718−0.1580.98129HN-20.997−0.6470.8970.6146−0.905−0.0310.075−0.5281110QL-10.803−0.480−0.3200.3537−0.1850.441−0.916−0.158811TR-12.296−0.2350.7661.43410.7330.014−0.5520.349512HNP2.2270.415−2.1891.08024.3720.5541.0362.854113QL-2−0.1441.256−0.3820.06681.067−1.162−0.0520.374414TR-21.550−0.6510.4910.8843−0.5612.737−0.9370.103715YS-21.0580.6550.2070.76541.598−0.631−0.3630.7513aF = 0.594 F1 + 0.164 F2+0.142 F3.bF = 0.589 F1 + 0.212 F2+0.156F3.Table 9. Component score and composite score ranking of aerial parts plus underground parts.Table 9NO.CodeF_1_F_2_F_3_Faranking1Aerial partsMY−2.623−0.589−0.803−1.577302ZK-1−0.9512.0601.2640.322133HZ−0.8780.8251.027−0.040154NQ0.4350.8890.7910.589105ZK-2−1.1572.2211.7080.332126TD0.6040.6850.3370.547117HN-11.3580.768−1.2080.71888YS-10.4410.698−0.7280.317149HN-21.4541.546−1.1071.009510QL-11.4380.580−0.6910.777711TR-12.5541.302−1.0781.481112HNP2.544−0.6740.0691.061413QL-20.5430.6891.3510.664914TR-21.8190.838−1.1920.967615YS-21.7151.3880.2511.287216Underground partsMY−2.415−0.203−0.328−1.2932917ZK-1−1.937−0.084−0.111−0.9912718HZ−1.1990.3760.871−0.3521819NQ−0.937−0.4380.185−0.5622120ZK-2−1.819−0.844−0.773−1.2532821TD−1.684−0.345−0.087−0.9412622HN-1−0.517−1.064−0.308−0.6122223YS-10.943−2.0970.542−0.0761624HN-2−1.487−0.208−0.396−0.8482525QL-1−1.041−0.686−0.707−0.8152426TR-1−0.222−1.195−0.249−0.4962027HNP3.368−2.7801.6271.070328QL-20.166−1.7420.365−0.3781929TR-2−0.924−0.314−0.651−0.6402330YS-20.409−1.6020.030−0.26617aF = 0.491 F1 + 0.294F2+0.144 F3.Table 10. Matrix correlation of altitude and analyzed constituents.Table 1012345. Altituder0.4140.188−0.0370.0820.295p0.1250.5020.8970.7710.2861 Loganic acid, 2 Swertiamarin, 3 Gentiopicroside, 4 Sweroside, 5 Isoorientin.

Discussion

4

In the Pharmacopoeia, G. straminea is regarded as one of the "Ma Hua Jiao", and the minimum content of gentianoside in its dried roots must be 2.0 % [39]. According to the results of our study, the minimum content of gentianosides in G. straminea populations was not less than 2.0 % (Table 5). This finding supports the use of all these populations as an original medicinal resource. This analysis of the constituents of G. straminea also suggested a plausible explanation as to why Tibetan doctors utilize the whole plant while Chinese doctors typically only the root of G. straminea. The principal component analysis was performed on data from the aerial and underground parts of G. straminea, respectively. After three components were extracted, the cumulative variance contributions from the two parts exceeded 85 % and satisfied the minimum requirements of analysis (Table 6). This indicated that the extracted components were sufficiently representative of the total information related to the five pharmacological components. Combining the results calculated by principal component analysis of the aerial and underground parts of G. straminea, respectively (Table 7, Table 8), confirmed that the concentrations of the five active components did not differ significantly between the aerial and underground parts from different regions. In brief, the principal component analysis is a single composite index representing all of the determined components from G. straminea; it is a collection of components and playing a leading role.

Unfortunately, calculations of data comparing the main components of the aerial and underground parts were insufficient to explain how the aerial parts of G. straminea were able to entirely replace the underground parts. Therefore, another principal component analysis was performed to comprise and unify all measured data from both the aerial parts and underground parts under a standard (Table 9). The results indicated that component F1 represented loganic acid, swertiamarin and gentiopicroside; component F2 represented sweorside; and component F3 represented isoorientin. Because of the pharmacological activity associated with loganic acid, swertiamarin and gentiopicroside, component F1 represented the major constituents responsible for G. straminea activity related to analgesic, anti-inflammatory, and febrifuge properties, as well as its protection against hepatitis. Sweorside can function as an anticonvulsant, so component F2 represented pluralities of constituents possessing the properties of an anticonvulsant in G. straminea. The pharmacological activity of isoorientin includes anticoagulation, anti-bacterial and anti-oxidation, so component F3 represented several compounds exhibiting the same pharmacological activity in G. straminea. Analysis of these three combined principal components and calculation of composite scores indicated higher ranking of the aerial parts from different regions, compared to the underground parts, thereby indicating that the aerial parts can replace the underground parts.

Comparison of different regions of G. straminea

4.1

Hierarchical clustering analysis was performed based on the data determined by HPLC of aerial parts of G. straminea from 15 different regions. Between-groups Linkage, a very efficient method for the analysis of variance between clusters, was applied, and Squared Euclidean distance was selected as a measurement. Principal component stacked contrast analysis and correlation analysis of the concentrations of the five components in relation to elevation were also carried out.

The dendrogram of clustering analysis showed that the G. straminea populations could be divided into four main clusters quality production regions according to the content of four marker compounds during the harvest season(Fig. 5): HN-1(7), TR-2(14), TR-1(11), HN-2(9), QL-1(10), YS-2(15), YS-1(8), NQ(4), TD(6), and QL-2(13) are cluster one; ZK-1(2), ZK-2(5), and HZ(3) are cluster two; MY(1) is cluster three; and HNP(12) is cluster four. Generally, the results of cluster analysis are more closely realistic. Interestingly, of the 15 populations, four populations contained superior constituent content to the others and these populations (Fig. 5). This suggests that there are several “geo-authentic” production regions of “Ma Hua Jiao”. The regions of HN-1(7), TR-2(14), TR-1(11), and HN-2(9) in cluster one is indeed good growth areas for G. straminea, and the determined contents of each component in G. straminea from those regions are higher. A sub-group consisting of NQ (4), TD (6), QL-2 (13), QL-1 (10), YS-2 (15), and YS-1 (8) from cluster one more closely approached the results suggested by the composite score of principal component analysis from aerial parts of G. straminea. HNP (12) represents a cluster for cultivation. In addition, using data from the component score of three components extracted from aerial parts, results of the hierarchical cluster analysis of the tested 15 samples were similar to those derived from the five active components (Fig. 6). Furthermore, stacked analyses of the five main components of G. straminea from the 15 regions were performed (Fig. 6), and it is clear and intuitive that HNP (12) corresponds to cluster four in Fig. 5, which is a well-grown region although it is an artificially cultivated sample, and that MY (1) corresponds to cluster three in Fig. 5, which is the region with the least amount of G. straminea components. It was also clear to see that gentianoside was the highest (see Fig. 6).Figs. 5 and 6Dendrograms of hierarchical cluster analysis for the 15 tested samples of aerial parts of G. straminea. The hierarchical clustering was done using SPSS software. A method called average linkage between groups was applied, and Squared Euclidean distance was selected as the measurement. In figure 5 (left), the dendrogram represents five active components, derived from data determined by HPLC of the tested 15 G. straminea samples. In figure 6 (right), the dendrogram represents data from a score of three components extracted from aerial parts of the tested 15 G. straminea.Figs. 5 and 6Fig. 6Stacked histogram of five major pharmacological constituents of the above ground portion of G. straminea.Fig. 6

G. straminea was characterised by a combination of dichogamy and herkogamy [40] and had high population genetic diversity [41]. And correlation analyses were carried out to determine the relationship between the concentrations of the five components and altitude. Although the concentrations of the loganic acid (Fig. 7a), swertiamarin (Fig. 7b), gentiopicroside (Fig. 7c), sweorside (Fig. 7d), and isoorientin (Fig. 7e) fluctuated between altitudes, Pearson correlation showed no significant relationship between the concentrations of the five components and altitude. The correlation coefficients were 0.414, 0.188, −0.037, 0.082 and 0.295, respectively (Table 9). Our results suggest that the populations of this species also have significant chemical variation among different populations on the Qinghai-Tibetan Plateau. Chinese medicinal material cultured in different localities is believed to differ in therapeutic potency [42]. Although genetic differences have been suggested to guide research and conservation of G. straminea [43], genetic proof is very limited and therefore concentrations of bioactive components should also be considered. Due to the high vulnerability of the region to climate change and human activities, over-harvesting of plants may lead to population scarcity. It is important to pay attention to the concentration of bioactive components when implementing policies for harvesting and conserving wild species, especially in areas where they are more suitable for growth. Therefore, the findings of this paper provide a strong reference value for sustainable species development and conservation.Fig. 7. Variations in the mean content of the five active constituents (a) Loganic acid, (b) Swertiamarin, (c) Gentiopicroside, (d) Swerosideand and (e) Isoorientin, over the altitude gradient.Fig. 7

Conclusions

5

This work first developed a simple and reliable HPLC method for simultaneous, quantitative determination of five active compounds in the aerial parts of G. straminea. Active constituents’ analysis of G. straminea showed that the aerial parts are slightly better than the underground parts, and the former can replace the latter in traditional medicines. Gentianolide was the most abundant constituent, with all populations having levels well above national standard levels. Finally, Hierarchical clustering analysis suggests that suitable growth areas for the above-ground parts of broad-leaved gentian were identified based on the content of labelled compounds. Although there was no significant relationship between the concentrations of the five components and altitude, our findings provide definitive phytochemical evidence for the conservation of resources and more careful utilization of this traditional Tibetan medicinal herb. Based on the analyses of the five constituents, it also seems reasonable to establish quality standards within the Tibetan Plateau Taoist medicinal culture region.

Limitations

6

This study has potential limitations. It should be noted that this paper only represents investigation of five active compounds from G. straminea, and it is slightly one-sided on behalf of all the ingredients. However, this study represents a preliminary attempt to apply statistical approaches towards solving practical issues, and further study will be needed in the future.

CRediT authorship contribution statement

Junlei Hao: Writing – review & editing, Methodology, Investigation, Data curation. Jiang Zhou: Data curation. Pengcheng Lin: Writing – review & editing. Jiang Wu: Methodology, Conceptualization.

Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.