Influence of Environmental Factors on the Active Substance Production and Antioxidant Activity in Potentilla fruticosa L. and Its Quality Assessment

Environmental factors may influence types and contents of active substances. This study investigated the influence of environmental factors on the active substance contents and antioxidant activity of Potentilla fruticosa L. from different regions of China. Also, HPLC fingerprint similarity analysis (SA) coupled with hierarchical cluster analysis (HCA) and discriminant analysis (DA) were further introduced for the accurate classification and quality assessment of P. fruticosa. The results showed that altitude was significantly and negatively correlated to the content of tannin (P < 0.05). Annual sunshine duration and altitude were significantly and positively correlated to the flavonoids content, rutin content and antioxidant activity (P < 0.05). Annual mean temperature was significantly and negatively correlated to the content of total phenolics, while altitude was significantly and positively correlated to the content of total phenolics (P < 0.05). Eight samples were unambiguously separated into three groups. Two types of discriminant functions with a 100% discrimination ratio were constructed. All data consistently supported the conclusion that P. fruticosa produced from Kangding, Sichuan Province had high quality among all samples, therefore, Kangding in Sichuan Province with favorable environmental conditions is recommended as a preferable production location.

The quality of traditional medicine, which plays a very important role in the health system of China, is determined by its active substances produced by the plants. Potentilla species have been used for a long time in traditional medicine for its curative properties. In Chinese traditional medicine Potentilla extracts have been used to treat diarrhoea, hepatitis, rheuma and scabies and as a remedy for detoxification 1,2 . Potentilla fruticosa L.is a species of hardy deciduous flowering shrub in the Potentilla genus of the family Rosaceae, native to the cool temperate and subarctic regions of the northern hemisphere, often growing at high altitudes in mountains 3 . Apart from common application as a garden plant, it also has numerous medicinal virtues 4 . Extracts of P. fruticosa have been shown to possess relatively high concentrations of phenolic acids and flavonoids and powerful radical scavenging capacity [5][6][7] . The activity of some extracts was higher than that of the synthetic antioxidant BHT (butylated hydroxytoluene) and of extracts isolated from Salvia officinalis L., which contains powerful antioxidants 7 . Moreover the leaves of P. fruticosa have applications as food additives and an ingredient in cosmetic products 8 . Similarly the same or other (local) Potentilla species have been used in traditional medicine of different cultures in Asia, Europe and Northern America. Modern scientific researches have confirmed that the material foundations of traditional Chinese herbal drugs are the different chemical constituents (most of them are secondary metabolites) that are contained in different raw plant materials 9 .
In China, P. fruticosa is commonly known as the "Jinlaomei medicine" and "Gesanghua" and is widely distributed in Qinghai, Gansu, Sichuan, Yunnan, Tibet and Heilongjiang, its altitude ranges from 400 to 5000 m 10,11 . For

Results
Validation of the HPLC procedure. The precision and repeatability of the method were assessed using seven injections of sample solutions and six replicates of dry solid plant samples (re-extracting), respectively 23 . For the common peaks ( Fig. 1, peaks 1-10) from duplicate injections, the relative standard deviations (RSD) of relative retention times (RRT) and relative peak areas (RPA) were found to be in the range of 0.02-0.06% and 0.22-2.91%, respectively (n = 7), which were calculated, respectively, to be 0.03-0.13% and 1.29-2.74% for six replicates of the solid samples (n = 6). To evaluate the accuracy of the method, we conducted a recovery experiment in which quantified analytes were mixed with specific amounts of standard components. The average percent recoveries for common peaks were in the range of 97.24 ± 0.02% to 104.31 ± 0.04%. The RSDs shifted from 1.19% to 2.31% (n = 6). The limit of detection (LOD) (signal/noise = 3) and the limit of quantification (LOQ) (signal/noise = 10) of the ten compounds (ten common peaks) varied within the range 1.95-3.12 ng/mL and 7.82-10.36 ng/mL, respectively. The stability of common peaks in the sample solutions that were maintained for 0-24 h was evaluated by determining their RPAs. The RSDs of RRTs and RPAs were found to be less than 3%. These results demonstrated that the conditions for the fingerprint analysis were optimal. Differences in active ingredient contents. The active ingredient contents of P. fruticosa leaves differed significantly because of their various growing locations (P < 0.05, Fig. 2). The highest tannin content (10.68%), Scientific RepoRts | 6:28591 | DOI: 10.1038/srep28591 flavonoids content (4.37%), and rutin content (0.79%) were found in leaves from the Kangding, Sichuan (S8), whereas the lowest tannin content (7.64%) and rutin content (0.19%) were observed in leaves from Shangri-la, Yunnan (S6). The lowest flavonoids content (3.29%) was found in leaves from Jingyuan, Ningxia (S4). Moreover, the leaves from Shangri-la, Yunnan (S6) had higher phenolic contents than any other growing location, with total phenolic contents at 9.40%. By contrast, leaves from the Jingyuan, Ningxia (S4) had the lowest total phenolic contents (5.29%). These rich differences may be due to ecological factors, genetic factors and the status of secondary metabolism in different growing locations. DPPH radical scavenging activity. For evaluation of antioxidant activity of P. fruticosa, DPPH radical scavenging activity 24 of different samples were compared and showed in Fig. 3. The scavenging effects of different locations increased with concentrations between 1 and 100 μ g/mL. Among all samples, the highest radical scavenging activity was obtained for P. fruticosa collected from S8 (Kangding, Sichuan) with the lowest average IC 50 value of 9.24 ± 0.04 μ g/ml, followed by S2 (Diebu, Gansu, 10.01 ± 0.06 μ g/ml). For P. fruticosa samples from eight regions, the IC 50 values changed with the order: Rutin (5.85 ± 0.02 μ g/ml) < S8 (9.24 ± 0.04 μ g/ml) < S2 (10.01 ± 0.06 μ g/ml) < S5 (14.25 ± 0.11 μ g/ml) < S1 (14.76 ± 0.21 μ g/ml) < S3 (17.63 ± 0.37 μ g/ml) < S7 (18.94 ± 1.04 μ g/ml) < S6 (23.68 ± 1.28 μ g/ml) < S4 (32.35 ± 2.65 μ g/ml). Compared with rutin standard, there was no significant difference (P < 0.05) on the antioxidant activities of S8, implying that S8 could have same scavenging effect with rutin standard at adequate concentration. These data indicated that P. fruticosa had high antioxidant activity, and the antioxidant activity of the same P. fruticosa species of different locations varied immensely from region to region.   For the same variable, bars with no letters in common are significantly different (P < 0.05). Y T (%), tannin content; Y F (%), total flavonoids content; Y R (%), rutin content; Y P (%), total phenolics content. S1, Mei county, Shaanxi; S2, Diebu, Gansu; S3, Huzhu, Qinghai; S4, Jingyuan, Ningxia; S5, Yongdeng, Gansu; S6, Shangri-la, Yunnan; S7, Ningchi, Tibet; S8, Kangding, Sichuan. from eight samples were significantly different (p < 0.05), which presented an inverted parabola (Fig. 4). In details, as to P. fruticosa derived from eight regions, the FRAP value ranged from 195.58 ± 6.86 to 416.58 ± 8.62 μ mol equiv. Trolox/g. P. fruticosa sampled from S8 (Kangding, Sichuan) possessed the highest antioxidant capacity with a FRAP value value of 416.58 ± 8.62 μ mol equiv. Trolox/g, followed by samples from S2 (Diebu, Gansu) and S1 (Mei county, Shaanxi) with values of 374.43 ± 6.47 and 337.91 ± 3.41 μ mol equiv. Trolox/g, respectively.
The FRAP values showed the same order of activity observed in the DPPH method: S8 > S2 > S1 > S3 > S5 > S7 > S6 > S4. Based on these results, we concluded that P. fruticosa samples presented not only free radical scavenge capacity but also reducing capacity. Also, P. fruticosa samples from Kangding, Sichuan (S8) was strongest in antioxidant capacity among all samples, which was close to the rutin and trolox standard. Furthermore, the antioxidant capacity of P. fruticosa samples exhibited significant differences in the same species from different regions (p < 0.05).

Analysis of environmental factors influencing the active ingredients and antioxidant activity. Principal component analysis (PCA) of environmental factors.
The PCA was conducted to identify the principal variables from a lot of independent variables (environmental factors). Contribution rate reflects the quantity of original information contained within each factor. The accumulated contribution rate of the first three eigenvalues reached 95.867% (Table 1), which indicates that the first three main components nearly covered total original information of the eighteen environmental factors. Thus, these components can be extracted to obtain the loading level of each environmental factor based on SPSS 19.0 software (Fig. 5). Figure 5 indicated that important environmental factors that effected the first principal component (F 1 ) were X N (0.872), X AAP (0.833), X ALTI (0.753), X pH (0.651), X AAT (0.434), X AMT (− 0.132), and X P (− 0.472). For the second principal component (F 2 ), important environmental factors were the X 7 (0.853), X 1 (0.815), X ASD (0.757), and X TP (− 0.123). The third principal component (F 3 ) accounted for a larger proportion in the X OM (0.761), X TN (0.723), X ALT (0.423), X AHT (0.411), X FFP (0.146), X TK (− 0.387) and X K (− 0.614) than that in other factors. However, its contribution rate was only 1.111%, and the F 3 was not be considered. Two principal components of environmental factors, F 1 (X N , X AAP , X ALTI , X pH , X AAT , X AMT , and X P ) and F 2 (X 7 , X 1 , X ASD , and X TP ) were thus screened for further analysis.
Gray correlation analysis (GCA). A GCA was conducted between the contents of active ingredients, antioxidant activity, and the principal components of environmental factors (F 1 , F 2 ) to deduce primary environmental factors 25,26 . The results of GCA (Table 2) demonstrated that environmental factors imposed different degrees of influences on different active ingredients and antioxidant activity: the primary environmental factor for the  accumulation of tannin was altitude (X ALTI ); for the accumulation of total flavonoids, rutin and antioxidant activity, were annual sunshine duration (X ASD ) and altitude (X ALTI ); for the accumulation of total phenolics, were annual mean temperature (X AMT ), annual average precipitation (X AAP ), and altitude (X ALTI ).
factors were further explored by PA, which makes the multiple statistical analyses more rational because the path coefficient possesses directional features 27 . The results of PA (Table 3) in this study demonstrated the effect of primary factors was significant on active ingredient contents and antioxidant activity. X ALTI (altitude) was key environmental factor for the content of tannin, X ASD (annual sunshine duration) and X ALTI (altitude) were key ones for flavonoids, rutin and antioxidant activity, while X AMT (annual mean temperature) and X ALTI (altitude) were key ones for total phenolics. The content of tannin was significantly and negatively correlated to altitude (P < 0.05). The flavonoids content, rutin content and antioxidant activity were significantly and positively correlated to annual sunshine duration and altitude (P < 0.05). The content of total phenolics was significantly and negatively correlated to annual mean temperature, while the content of total phenolics was significantly and positively correlated to altitude (P < 0.05). No significant correlation was observed between the soil factors and the contents of each active substance or antioxidant activity. The statistical significance demonstrated that P. fruticosa samples from Kangding, Sichuan Province (S8) was selected as the most valuable materials due to the higher phytochemical contents and significant antioxidant activity. Furthermore, Kangding, Sichuan Province was proved to be the most suitable region for producing P. fruticosa through the comprehensive investigation and analysis of the effect of environmental factors.  Table 3. Path analysis between the active ingredient contents, antioxidant activity and primary factors. Primary factors including X AMT , X AAP , X ASD , X ALTI were performed in the  (Fig. 1). These herbs showed different chromatographic profiles. After carefully analyzing the fingerprint profiles of these samples, ten common peaks with acceptable heights and good resolution were selected as characteristic peaks for the identification and classification of the P. fruticosa samples from different production regions 16,19,21,22 . The relative retention time and relative areas of these constituents (characteristic peaks) were calculated with respect to the reference peak eluting at a retention time of 18.0 min (peak 4). The similarities among the 8 chromatograms were evaluated based on a matrix comprising the relative areas of the 10 constituents of samples 1-8 by the CASE software. The results showed that the similarity of chromatograms of eight samples had a wide range and varied within the range of 0.21-0.89 (Table 4), which indicated that the chromatographic patterns of the samples varied greatly in different production regions.
Hierarchical clustering analysis (HCA). Prior to HCA, we compared the fingerprints visually and simply divided the samples into three distinct groups, namely A, B, and C ( Fig. 6). Although it was possible to differentiate chromatograms on the basis of visual comparison, this process was subjective and non-quantitative. HCA can provide a quantitative and objective analysis of fingerprints 20 . HCA was performed according to the relative peak areas of ten constituents and a dendrogram was acquired ( Fig. 7). At the rescaled distance of 12.5, HCA can divide the eight samples from different locations into three groups consistent with the results of visual comparisons. S1 and S2 were merged into a group G1; S3, S4, S5, S6, and S7 were grouped in G2; and S8 was found in G3. The correlation coefficients of each chromatogram within groups G1, G2 and G3 corresponding to software-generated group simulated mean chromatograms, and the correlation coefficients between these simulated mean chromatograms were shown in Table 5. The chromatograms within a particular group were generally consistent. Correlation coefficients for each chromatogram which were classified into a particular group to the corresponding simulative mean were higher than 0.90. However, the chromatograms within a particular group were markedly different from the chromatograms in different groups, and the differences in the similarity values for the three groups (Table 5) reflected the HCA-generated data. Meanwhile, the correlation coefficient between the software-generated group simulated mean chromatograms was lower than 0.65, and ANOVA investigations indicated the difference between groups was significant (P = 0.041). These results demonstrate that HCA can distinguish P. fruticosa species with the same origin from different districts.
Discrimination analysis (DA). Ten common peaks were selected from the fingerprints, thereby creating 10 variables. However, not all these variables are of value to establish discriminant function. The DA will generate discriminant functions only by use of valuable variables. Two types of discriminant functions [equations (1) and (2)] were obtained using the SPSS software (SPSS for Windows 19.0, SPSS Inc., USA).
Canonical discriminant function: Discrimination standard: and Y Y G 0 : 0 : 0 : Fisher's discrimination function:    Discriminant standard: each sample has three functional values and is assigned to the group corresponding to the highest of these function values. G1 [equation (3)], G2 [equation (4)] and G3 [equation (5)] denote the samples from groups G1, G2 and G3 respectively, and X denotes the variable.
Only three variables, corresponding to the areas of peaks 2, 5 and 8, were used to generate the discriminant functions. In order to place an unknown sample, the values of the three variables are inserted into the equations, and the unknown sample was grouped according to the discriminant standard value obtained. Use of the three most discriminating variables enable test samples belonging to groups G1, G2 and G3 to be classified with 100% accuracy. The well-resolved DA plots for the three groups are shown in Fig. 8. DA approach supported a division of the 8 samples into 3 major groups for which HPLC SA confirmed the considerable variation in the active substances and quality of P. fruticosa from different growing locations. Furthermore, the samples from Kangding, Sichuan Province (S8) was separately assigned into a special group Group 3 (Fig. 8), and also far away from Group 1 and 2. There was a huge overlap between the center of the samples contained by Group3 (S8) and group centroid, indicating that S8 was in high quality.

Discussion
Variation in active substance contents and antioxidant activity. Altitude, temperature, illumination, and moisture are important factors to influence the metabolism and accumulation of secondary metabolites 28 . Environmental differences (such as altitude, temperature, illumination, precipitation, humidity, soils) in different production locations contribute to the differences in active ingredient contents and antioxidant activity of medicinal plants 9,16,19 . In this study, significant differences were observed in the contents of the chemical compositions and antioxidant activity of P. fruticosa leaves obtained from different growing locations (Figs 1~4). The contents of tannin in the leaves of P. fruticosa planted in S8 (Kangding, Sichuan) were higher than those planted in other locations, followed by S1 (Mei county, Shaanxi). The contents of total flavonoids were high in those planted in S8 (Kangding, Sichuan), followed by S7 (Ningchi, Tibet). Li et al. also discovered that the contents of total flavonoids in P. fruticosa leaves from different environments displayed great differences 29 , which agrees with 0.999 ± 0.000 a Table 5. Correlation coefficients between individual chromatograms within a group and the group simulative mean chromatogram, and between the group simulative mean chromatograms. a Correlation coefficient of individual chromatograms to the simulative mean chromatogram of the corresponding group. b Correlation coefficient between simulative mean chromatograms. our findings. The contents of rutin were also high in S8 (Kangding, Sichuan), followed by S7 (Ningchi, Tibet). The contents of total phenolics were high in S6 (Shangri-la, Yunnan), followed by S8 (Kangding, Sichuan). The antioxidant activity was high in S8 (Kangding, Sichuan), followed by S2 (Diebu, Gansu). The results indicate that the contents and properties of active substances are closely related to the growing locations, and are determined by environmental factors of the growing locations. Integrating the active substance contents, antioxidant activity ability and the effect of environmental factors, P. fruticosa from Kangding, Sichuan Province (S8) was selected as the most valuable material, and Kangding, Sichuan Province was considered as the most suitable region for producing P. fruticosa. The rich variation in chemical compositions and antioxidant activity may lead to significant differences in effectiveness as medicines, functional foods and nutritional supplements. HPLC fingerprint SA coupled with HCA and DA were further introduced to establish an effective method for the accurate classification and quality assessment of P. fruticosa. This method performed well for the quality evaluation of P. fruticosa samples and also confirmed that S8 had high quality.  Table 6. Similarly, in this study, the altitude in Kangding is higher than other locations, the contents of total flavonoids were the highest, and significantly different from other locations (Fig. 2). GCA and PA showed that maximum degree of correlation was found between altitude and the contents of total flavonoids, and positive relationship existed (Tables 2 and 3). Flavonoids contain ortho-dihydroxylated structure, which exhibits the characteristics of absorbing ultraviolet, and scavenging free radicals. Significant positive correlation between the contents of total flavonoids and altitude are responses of the P. fruticosa species distributed at higher altitude to the strong ultraviolet radiation. Also, UV-absorbing compounds such as flavonoids and phenolic acids have been proven to be important for plants to protect their vegetative organs from harmful UV radiation 13 . By contrast, the total amount of flavonoids of some medicinal plants, such as Buxus sempervirens L., decreased with increased altitude 31 . The reduction of flavonoids in B. sempervirens leaves with altitude might be a consequence of an increased synthesis of wax in response to the decrease in temperature 32 since both, cuticular waxes and flavonoids, are derived from acetyl-CoA 33 .

Altitude and active substance contents.
Temperature and active substance contents. The production of secondary metabolites by temperature stress is the expression of self defense mechanism of medicinal plants. Recently a significant increase of the biosynthesis of phenolics was demonstrated for plants growing in a low temperature regime 34 . Similarly, in this study, annual mean temperature was significantly and negavely correlated to the contents of total phenolics. This result agreed with the previous researches for total phenolics contents involved in other medicinal plants. For example, in walnut (Juglans regia L.), a negative correlation was found between temperature and total phenolic contents 30 . Albert et al. reported that enhanced UV-B radiation is probably not the key factor triggering shifts in the phenolic composition in Arnica montana L. cv. ARBO grown at higher altitudes but rather temperature, which displays negative correlation with the phenolics 35 . Reyes et al. studied total phenolics during development of Solanum tuberosum L. in cultivars grown in Texas and Colorado and found cooler temperatures in Colorado favored about a 1.4-times higher total phenolics content than in Texas-grown tubers 36 . The present study also showed phenolic compounds increased with altitude, similar the trends as the results for Buxus sempervirens obtained in a previous study 31 . Usually a decrease of the mean temperature of 0.55 °C per 100 m of altitude is observed 37 . Hence, the contents of phenolic compounds increased with the decrease of temperature. However, temperature is positively correlated to the contents of active ingredients in some plants. The contents of hypericin in Hypericum perforatum L. increased with temperature 38 . The content of anthocyanin in sugarcane that was under high temperature obviously increased 39 . Plants in alpine habitats are exposed to many environmental stresses, in particular temperature and radiation extremes. Another laboratory studies confirm that the biosynthesis of certain phenolic compounds is increased under increased UV radiation, whereas shading leads to a decrease in the biosynthesis of these compounds 40 . Therefore, plant could have a chemical adaption to the alpine environment, indicating that the impacts of ecological factors on the secondary metabolites are related to their chemical types, structures and characteristics.

Sunshine duration and active substance contents.
Illumination also affects the synthesis and accumulation of secondary metabolites in medicinal plants. The increase of illumination time can increase the contents of secondary metabolites. For example, the amount of flavonoids in Arabidopsis increased after long time illumination 41 . The contents of ginsenogsides in Panax quinquefolius L. were positively correlated to the annual sunshine duration 42 . Dong et al. investigated the influence of environmental factors on the contents of secondary metabolites in the leaves of Eucommia ulmoides and found that annual sunshine duration was significantly and positively correlated to the contents of geniposidic acid (P < 0.05) 9 . In this study, PA showed that the contents of active substances (tannin, total flavonoids, rutin, and total phenolics) were highly associated to annual sunshine duration, positive correlation between them was observed ( Table 3). The results were also verified in the study on production of active substances of other medicinal plant such as Sinopodophyllum hexandrum 16 . These data indicated that annual sunshine duration is a key environmental factor to the synthesis and accumulation of active substances. P. fruticosa is a heliophilous plant, locations with long time sunshine would be favorable for its secondary metabolism, resulting in adequate substrate to catalyze the synthesis of the active substances and thereby increasing their contents.
Precipitation and active substance contents. In this study, the contents of tannin, rutin and total phe- Zhang could be as high as 1% under dry conditions, which is 2.5-times greater than that obtained under humid conditions (0.4%) 16 . In the meantime, the study showed that the contents of total flavonoids were positively correlated to annual average precipitation. Although the contents of tannin, rutin and total phenolics were negatively correlated to the annual average precipitation, and the contents of total flavonoids were positively correlated to annual average precipitation, the impacts of annual average precipitation on each secondary metabolite were not significant (Table 3). Therefore, precipitation is not a limited factor to the synthesis and accumulation of secondary metabolites for P. fruticosa. Additionally, 500 mm annual precipitation can meet the demand of the growth of P. fruticosa according to the biological features and the demands for environmental conditions. In the field investigation, the annual average precipitations in 8 locations are all over 500 mm, they all can meet the moisture demand from the growth of P. fruticosa. This could explain why the influence of annual average precipitation was not significant on active substances of P. fruticosa.
Soil and active substance contents. Some researchers have stated that the variations in the active ingredient contents of the plants are related to soil fertility. For example, the leaf polyphenol content of young greenhouse tomato plants was increased considerably in response to low N availability 43 . The saponin contents of Panax quinquefolius were increased by 27.86% by the repeated applications of high-quality organic fertilizer and the functionality of potassic fertilizer was the same as that of a total nutrient admixture in terms of saponin contents 44   Similarly, altitude showed significant and positive correlation to the antioxidant activity in this study. These data suggest potential influences of growing conditions on the antioxidant properties of P. fruticosa and the possibility of producing a selected wild species with a desired antioxidant property by optimizing growing conditions. The significant effects of environment must, therefore, be considered. However, few specific data are available in the literature on the relationship between environmental factors and antioxidant activity of P. fruticosa. More research is required to further investigate the cause of environmental factor effects using a greater sample of both field locations and varieties. Going forward, we would conduct further studies in cells and living bodies to fully reflect the antioxidant properties of P. fruticosa, which could not be achieved by the DPPH scavenging capacity assay and FRAP assay carried out above.
Quality assessment of P. fruticosa. Traditional Chinese medicine (TCM), which has a 5000-years history of application, is still in wide demand 48 . TCMs are being closely examined for the development of novel pharmaceuticals 49 . However, quality control of TCM has always been a bottleneck for their adoption worldwide because of their complexity, the presence of unknown components and lack of quality control. Chromatographic fingerprint analysis has been introduced as a rational strategy for the assessment of complex TCMs. Unlike the traditional practice that one or more compounds were selected as active markers for identification and quality assessment, fingerprint relies on the inherent relationship of multiple compounds and displays the chemical pattern of TCM. Among available quality control methods, chromatographic fingerprint has gained more and more attention recently. In this study, HPLC fingerprint coped with chemometric methods were applied to deal with the classification for the quality assessment of P. fruticosa from different regions of China. The results showed that the eight samples were divided into 3 clusters. Two types of discriminant functions were generated using three selected predictor variables and the ratio of discrimination was 100%. The results of HPLC fingerprint SA, HCA, and DA in this study were in agreement, and the techniques can identify and classify the same species from different regions and performed well for the quality assessment of P. fruticosa. Also, it can be used to compare and control other natural products prepared from them. Many similar studies on other medical plants have been published for the quality evaluation [20][21][22][23] . The quality of herbal drugs is affected by many composite factors. Climatoecological type, symbiosis ecotype, and soil ecotype in growing locations are biologically essential to produce high-quality geoauthentic herbal drugs 9 . For the widespread plant species, because the large span in geological zones, differences in chemical constituents within the same plant species would occur. P. fruticosa belongs to one of widespread plant species, distributed in the Asia, Europe and North America. Currently, the influence of environmental factors in other locations on active substances of P. fruticosa, and its quality assessment are less studied. Further research include more samples from other regions of China and other countries should be carried out. The present study provides meaningful information for the collection and application of wild P. fruticosa herbs in both healthcare and the food industry. Moreover, although this present study relates only to P. fruticosa, the methodology has wider relevance for quality assessment of other edible and medicinal plants.

Methods
Plant materials and soil samples. P. fruticosa were collected from eight representative locations located in seven provinces of China during 10-25th July in 2014 (Fig. 9). Specifically, a total of eighty healthy specimens from four populations (each population was separated geographically by at least 30 km, and 5 m for adjacent individuals) with similar growth stature were collected in every test region, as shown in Table 7. Matured leaves were picked up respectively from four directions (north, south, east, and west) of three positions (up, middle and low part) of the plant and then mixed as one sample for each test location, thereby obtaining eight test samples. All samples were dried under the vacuum at 40 °C. After grinding into powders, the samples were stored in the dark at − 20 °C before further use. Simultaneously, soil rhizosphere were collected and treated for each study site, thereby obtained eight soil samples from eight study sites for the measurement of soil factors. Voucher specimens from all populations were identified by professor Jianjun Liu of Northwest A&F University and were deposited at the Herbarium of Northwest A&F University (WUK0781792-0781823).
Related data of environmental factors. Soil samples were used to determine the key soil parameters including rapidly available nitrogen (X N ), rapidly available phosphorus (X P ), rapidly available potassium (X K ), organic matter (X OM ), total nitrogen (X TN ), total phosphorus (X TP ), total potassium (X TK ), and pH (X pH ) using the modified soil chemistry analysis methods described by Sparks et al. 50 . Related data of climate factors including annual mean temperature (X AMT ), january average temperature (X 1 ), july average temperature (X 7 ), annual accumulated temperature (≥ 10 °C) (X AAT ), annual highest temperature (X AHT ), annual lowest temperature (X ALT ), annual average precipitation (X AAP ), annual sunshine duration (X ASD ), and frost free period (X FFP ) in the 50 years (1964~2013) was collected from local meteorological bureaus (stations) for the eight study sites. The environmental factors of the eight study sites were summarized in Table 6.
Preparation of the extracts. Considering the impact of various factors, the extraction process was optimized using a response surface method 51

Measurement of tannin.
Tannin contents were determined using the Folin-Denis method 52 . 1.0 mL of the diluted samples (2 mg/mL) was transferred to a 25 mL calibration flask to which were added 1.0 mL of Folin-Denis reagent and 5 mL of sodium carbonate (1 mol/L). The solution was diluted to 25 mL by the addition of methanol. After incubating for 30 min at room temperature, the absorbance at 720 nm was measured against a blank. Tannin acid (1 mg/L to 10 mg/L) was used for the standard curve calibration. All measurements were performed in triplicate.
Measurement of total flavonoids. Total flavonoids contents were determined by sodium nitrite-aluminum nitrate colorimetric method 53 . Sample solution (1.0 mL, 2 mg/mL) was transferred into 25 mL volumetric flasks, then 0.3 mL NaNO 2 (5%) was added and held for 6 min. Next, 0.3 mL Al(NO 3 ) 3 (10%) was added and held for another 6 min. Finally, 4 mL NaOH (1 mol/L) was added and the solution was diluted to 25 mL with 70% ethanol solution. After 30 min of incubation at room temperature, the absorbance at 510 nm was measured against a blank. Rutin (4 mg/L to 40 mg/L) was used for the standard curve calibration. All measurements were performed in triplicate.
Measurement of total phenolics. The total phenolic contents were determined using a modified Folin-Ciocalteau colourimetric method 54 . Sample solution (1 mL, 2 mg/mL) was transferred to a 25 mL calibration flask, 0.5 mL of Folin-Ciocateu reagent and 2.5 mL of sodium carbonate (1 mol/L) was added and final volume was marked with methanol. After 30 min of incubation at 30 °C, the solution was centrifuged at 4,000 rpm (10 min) and the absorbance at 760 nm was measured against a reagent blank. Gallic acid (0 mg/L to 6 mg/L) was used for the standard curve calibration. Three replicates were made for each sample.

RP-HPLC analysis.
Sample solution (1 mg/mL) was filtered through a 0.22 μ m microporous filtering film and were then separated by RP-HPLC at ambient temperature to obtain chromatograms 55 . The amounts of rutin were also quantified by RP-HPLC. The contents of rutin were identified and calculated by a comparison of the relative retention time (RRT) and relative peak areas (RPAs) with those of the standards. HPLC analyses were conducted on an Agilent Series 1260 liquid chromatograph, equipped with a quaternary gradient pump, a variable wavelength detector system and a reversed-phase SB-C18 column (5 μ m, 4.6 × 250 mm, Agilent, USA). The mobile phase consisted of water with 0.5% acetic acid (mobile phase A) and methanol with 0.5% acetic acid (mobile phase B). The flow rate was 1 mL/min. The gradient program was as follows: 0-15 min, 25% to 35% B; 15-30 min, 35% to 60% B; 30-40 min, 60% to 100% B; 40-45 min, hold B at 100%. The injection volume and detect wavelength were 20 μ L and 360 nm, respectively 56,57 . Analyses were performed in triplicate.
DPPH radical scavenging activity. DPPH assay has been widely used for the determination of antioxidant activity of pure antioxidant compounds as well as of different plant extracts, a lower IC 50 representing stronger antioxidant capacity 24 . IC 50 obtained by interpolation from linear regression analysis is the effective sample concentration at which DPPH radicals is scavenged by 50%. A slightly modified DPPH method was used to determine the radical scavenging property of P. fruticosa leaves 58 . Amounts of 2 mL of the tested samples (1-100 μ g/mL) and the positive controls (rutin, 1-70 μ g/mL) were mixed with 2.0 mL of 0.1 mol/L DPPH in methanol. The mixture was vortexed thoroughly and allowed to stand in the dark for 30 min. The absorbance was measured at 517 nm spectrophotometrically against a blank. All measurements were performed in triplicate. DPPH free radical scavenging activity (SA) can be expressed with the following equation (6): where A i : absorbance of tested samples (2 mL) mixed with DPPH (2 mL); A j : absorbance of tested samples (2 mL) mixed with methanol (2 mL); A 0 : absorbance of methanol (2 mL) mixed with DPPH (2 mL).

Ferric reducing antioxidant power (FRAP) assay.
In this assay, the colour changed from yellow to green depending on the reducing power of the sample. Fe 3+ /ferricyanide complex was reduced to the ferrous form by reductants in the solution. Thus, Fe 2+ can be monitored as the reducing power index by measuring the absorbance. The reducing capacity of the extract was determined by the method of Oyaizu 59 . Varying concentrations of the extracts in the methanol (2.5 mL, 50-500 μ g/mL) were mixed with phosphate buffer (2.5 mL, 0.2 mol/L, pH 6.6) and potassium ferricyanide (2.5 mL, 1%), and incubated at 50 °C for 20 min. Aliquots (2.5 mL) of 10% trichloroacetic acid were added to the reaction mixture, which was then centrifuged at 4,000 rmp (10 min). The upper layer of the solution (2.5 mL) was mixed with distilled water (2.5 mL) and ferric chloride (0.5 mL, 0.1%). The absorbance was measured at 593 nm. The FRAP results were expressed in terms of micromoles trolox equivalent per gram dry weight of samples (μ mol equiv. Trolox/g). All measurements were performed in triplicate.

Data analyses.
Three methodologies were performed step-by-step to analysis systematically the influence of environmental factors on the active substance contents and antioxidant activity of P. fruticosa leaves, which also represented a quantitative comparative analysis of the dynamic process of plant and the environment. Firstly, principal component analysis (PCA) was carried out using SPSS software (SPSS for Windows 19.0, SPSS Inc., USA). PCA has the strong advantage of significantly reducing the dimension of the complex data while preserving most of the variance within by using dependencies among large numbers of variables without requiring knowledge of the data set in order to visualize high dimensional data and identify the most important variables 60 . PCA in the present study was used to find out principal components of environmental factors. Secondly, gray correlation analysis (GCA) between the contents of active substances, antioxidant activity and the principal components of environmental factors obtained by PCA were conducted by DPS7.5 software (Date Processing System, Science Press, China) to screen primary factors from principal components of environmental factors. The result of GCA reflects the close degree of the relationship between principal behavior factors and other factors, to deduce primary factors and secondary factors based on the gray correlation degrees of different factors 9 . High correlation degree indicates high influence of the factor to the accumulation of active substances and antioxidant activity 25,26 . GCA is usually applied in the indefinite system, including small sampling and poor data information system, in which limited information is available, others are not known 25,26 . Therefore, it is suitable to be used for the analysis of the data obtained in this study. Thirdly, path analysis (PA) deals with the quantitative relationship between dependent and independent variables to explain the relative significance of each factor to the dependent variables 9,27 . Figure 10 showed the path network of four independent variables included in PA. PA can determine whether the effect of X i on Y is significant or not, and can identify not only the direct effect of X i on Y (b i , X i -Y) but also indirect effect of X i on Y through X j (r ij b j , X i -X j -Y, i≠ j). Thereby, the correlation coefficient (r iy ) contains the direct path coefficient (b i , X i -Y) and the indirect path coefficient (r ij b j , X i -X j -Y, i≠ j) = + ∑ ≠ r b b r ( ) iy i j j ij 1 . PA was conducted by statistical software SAS 9.1 (SAS Institute, Cary, NC, USA) to evaluate correlation between the active ingredients, antioxidant activity and primary environmental factors Figure 10. Network of path analysis. X 1 , X 2 , X 3 and X 4 are independent variables; Y is dependent variable. r ab represents the correlation coefficient between independent variables a and b.
Scientific RepoRts | 6:28591 | DOI: 10.1038/srep28591 obtained by GCA. Environmental factors were used as independent variables, and active ingredients and antioxidant activity were used as dependent variables in each test. Furthermore, in order to further evaluate the quality of P. fruticosa from different growing locations, another three methods viz. HPLC fingerprint analysis coupled with chemometric methods were applied. Firstly, HPLC fingerprint similarity analysis (SA) was performed using Computer Aided Similarity Evaluation software (CASE 2004, Zhejiang University, Hangzhou, China) as recommended by the Chinese Pharmacopoeia Committee. The software is used for evaluating similarities between different chromatograms based on the cosine values of vectorial angel 61 . The cosine values of the two chromatograms approaching 1 means they are highly similar. This software was also used to compute the mean chromatogram as a representative standard chromatogram for a group of chromatograms. The standard HPLC fingerprint is set up with the median of all chromatograms [62][63][64] . Subsequently, hierarchical clustering analysis (HCA) and discrimination analysis (DA) were performed using SPSS 19.0 19,20 . The 'average linkage between groups' method was applied and the Pearson correlation was selected as a measurement 65 . DA can be used to build a predictive model of the group membership based on observed characteristics in each case. This procedure can generate a discrimination function or a set of discriminant functions from the samples with known membership based on linear combinations of the predictor variables that provide the best discrimination among the groups. The functions can be applied to discriminate and classify new cases with measurements for the predictor variables but with unknown group membership 20 . A detailed description/theory of statistical methods/models was provided in supplementary Text S1. Normalization of data was performed by the Z score transformation normalization method if necessary 66 .
The results were presented as the mean value ± SD (standard deviation). The data was analyzed by one-way ANOVA (Analysis of Variance) followed by Duncan multiple comparison (p < 0.05) based on SPSS software (SPSS for Windows 19.0, SPSS Inc., USA).