Next Article in Journal
Investigating the Potential of Ghee Precursor-Derived Carbon Nano Onions for Enhancing Interfacial Bonding in Thermoplastic Composites
Previous Article in Journal
Design of Molecularly Imprinted Polymers Using Supercritical Carbon Dioxide Technology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Environmental Factor Effects on the Polyphenol and Flavonoid Content in the Leaves of Chrysanthemum indicum L. and Its Habitat Suitability Prediction Mapping

1
Department of Medicinal Plant, Graduate School of Pharmaceutical Sciences, Kumamoto University, No. 5-1, Oe-Honmachi, Chuo-ku, Kumamoto 862-0973, Japan
2
BVRIT HYDERABAD College of Engineering for Women, Nizampet Rd, Hyderabad 500090, Telangana, India
3
Global Center for Natural Resources Sciences, Kumamoto University, No. 5-1, Oe Honmachi, Chuo-ku, Kumamoto 862-0973, Japan
*
Authors to whom correspondence should be addressed.
Molecules 2024, 29(5), 927; https://doi.org/10.3390/molecules29050927
Submission received: 28 December 2023 / Revised: 24 January 2024 / Accepted: 5 February 2024 / Published: 20 February 2024

Abstract

:
The leaves of Chrysanthemum indicum L. are known to have various bioactive compounds; however, industrial use is extremely limited. To overcome this situation by producing high-quality leaves with high bioactive content, this study examined the environmental factors affecting the phytochemical content and antioxidant activity using C. indicum leaves collected from 22 sites in Kochi Prefecture, Japan. Total phenolic and flavonoid content in the dry leaves ranged between 15.0 and 64.1 (mg gallic acid g−1) and 2.3 and 11.4 (mg quercetin g−1), while the antioxidant activity (EC50) of the 50% ethanol extracts ranged between 28.0 and 123.2 (µg mL−1) in 1,1-Diphenyl-2-picrylhydrazyl radical scavenging assay. Among the identified compounds, chlorogenic acid and 1,5-dicaffeoylquinic acid were the main constituents in C. indicum leaves. The antioxidant activity demonstrated a positive correlation with 1,5-dicaffeoylquinic acid (R2 = 0.62) and 3,5-dicaffeoylquinic acid (R2 = 0.77). The content of chlorogenic acid and dicaffeoylquinic acid isomers varied significantly according to the effects of exchangeable magnesium, cation exchange capacity, annual temperature, and precipitation, based on analysis of variance. The habitat suitability map using the geographical information system and the MaxEnt model predicted very high and high regions, comprising 3.2% and 10.1% of the total area, respectively. These findings could be used in future cultivation to produce high-quality leaves of C. indicum.

Graphical Abstract

1. Introduction

Chrysanthemum species have been widely used for a long time as food, beverages, cosmetics, and medicines in Asian countries [1]. The flowers are a part of Chinese, Korean, and Japanese traditional medicines used to treat various conditions, such as inflammation, infection, fever, headache, and eye fatigue [1,2]. Among the genus Chrysanthemum, only two species, C. indicum L. and C. morifolium Ramat. are listed in the Japanese Pharmacopoeia [3]. Many studies have investigated the therapeutic potential of these flowers as antioxidant [4], anti-inflammatory [5], hepatoprotective [6], antidiabetic [7], neuroprotective [8], tyrosinase inhibitory [9], and anti-allergic agents [10]. These activities were attributed to the abundance of polyphenols and flavonoids [11,12].
Polyphenols, such as chlorogenic acid and dicaffeoylquinic acid isomers and flavonoids, such as the glycosides of luteolin, apigenin and acacetin are of great dietary importance in Chrysanthemum flowers. Sesquiterpenoids, with their various pharmacological activity, are also considered to be important constituents in Chrysanthemum flowers [13]. Chlorogenic acid has various activities, such as anti-inflammatory [14], antioxidant [4], and antidiabetic activities [15]. Similarly, dicaffeoylquinic acid isomers exhibit antioxidant [4] and anti-inflammatory [16]. Luteolin glucoside exhibits cardioprotective [17] and anti-inflammatory activities [18]. Likewise, apigenin glucoside has antifungal [19] and anti-inflammatory activities [20]. Acacetin derivatives exhibit cardioprotective [21] and anticancer effects [22].
Regarding the leaves of Chrysanthemum species, previous studies have reported that Chrysanthemum leaves exhibit histamine release inhibition [10], antimicrobial, anti-inflammatory, and antioxidant activities [23,24,25]. Recently, C. morifolium leaves were reported to contain chlorogenic acid, 3,5-dicaffeoylquinic acid, luteoloside, and quercetin [26]. These studies confirmed that Chrysanthemum leaves have medicinal value.
Despite the pharmacological interest in the leaves, industrial use is extremely limited. Producing high-quality Chrysanthemum leaves with high phytochemical content would be an effective approach for the wider use of the leaves in various fields such as food, beverages, cosmetics, and medicines. This study, therefore, examined environmental factors affecting the phytochemical content and antioxidant activity of C. indicum leaves.
C. indicum is widely distributed in western Japan (Kyushu, Shikoku, and western Honshu) and grows under various environmental conditions that may significantly affect the accumulation of secondary metabolites [27]. For instance, previous studies have shown that water stress enhances the phenolic content of C. indicum [28]. Furthermore, the antioxidant capacity of C. indicum was reported to be affected by light quality [29]. Statistical analysis and geographic information systems (GIS) are useful for examining various environmental factors [27].
GIS provides efficient statistical and spatial modeling approaches for preparing habitat suitability maps (HSMs) by integrating several environmental datasets. HSMs define the relationships between species and their environments by understanding their geographical distributions. These spatial models involve statistical analyses to simplify the probability of species existence based on identified correlations between various environmental conditions and the presence or absence of species [30]. The MaxEnt model, utilizing presence occurrence data and employing stochastic generation for associated points, proficiently explains the maximum entropy distribution [31]. Its efficacy extends notably to presence-only datasets featuring constrained sample sizes, demonstrating robust performance even under conditions of data incompleteness [32].
In this study, C. indicum leaves were collected from 22 sites in Kochi Prefecture, one of the natural habitats of C. indicum. Total phenolic and flavonoid content, antioxidant activity, and the chemical composition were analyzed. Topographical, soil, environmental, and climatic factors were recorded and used for statistical analyses and GIS studies. The HSM of the plant in Kochi was created to predict locations with a high probability of plant growth.

2. Results and Discussion

2.1. Soil Analysis

Soil samples collected from the study area showed differences in macro- and micronutrient contents and other soil properties (Table S1). The soil in the study area was a mixture of gravel, coarse sand, and rocks. The soil was slightly acidic, with a pH range of 5.2–6.8. It was rich in ammonium nitrogen (1.7–7.8 mg (100 g)−1), nitrate nitrogen (0.1–4.6 mg (100 g)−1), phosphorus (5–46 mg (100 g)−1), exchangeable potassium (11–54 mg (100 g)−1), calcium (39–1230 mg (100 g)−1), and magnesium (10–514 mg (100 g)−1). The pH (5.2–6.8), soil moisture sensor SM150T output (0.11–0.51 V), and soil bearing capacity (3.6–16.8 t sf−1) were measured during the field investigation. The soil moisture (8.0–36.2% vol) was calculated from the SM150T output.

2.2. The Phytochemical Content and Antioxidant Activity

Total phenolic and flavonoid content and 1,1-Diphenyl-2-picrylhydrazyl (DPPH) radical scavenging activity are shown in Table 1. Although these measurements demonstrated large variation, more than half of samples resulted in very high antioxidant activity (EC50 < 50 µg mL−1). These results indicated that C. indicum leaves have potent antioxidant activity. The relative strength of DPPH radical scavenging (1/EC50) correlated with both total phenolic content (R2 = 0.83) and total flavonoid content (R2 = 0.42), as shown in Figure 1. This result is consistent with that of a previous study that showed a correlation between antioxidant activity and phytochemical content [4].
Liquid chromatography-mass spectrometry (LC-MS) and high-performance liquid chromatography (HPLC) were used to analyze the samples [33,34]. The seven major peaks identified were chlorogenic acid, 1,5-dicaffeoylquinic acid, 3,5-dicaffeoylquinic acid, 4,5-dicaffeoylquinic acid, acacetin 7-O-glucoside, acacetin 7-O-malonylglucoside, and acacetin (Figure 2 and Table 2). The identified compounds were quantified by HPLC, and the results are summarized in Figure 3 and Figure 4. Chlorogenic acid and 1,5-dicaffeoylquinic acid contents varied from 0.3 to 3.3 and 0.6 to 7.5 (% dry weight), and were considered the main phenolic compounds in the C. indicum leaf extract. C. indicum leaves had little or no flavonoid compounds except for acacetin and the derivatives, whereas C. indicum flower contains various flavonoids including luteolin, apigenin, and acacetin, mainly in their glycoside forms [9]. These results suggest that the composition of C. indicum leaves chemically differs from its flowers.
Total phenolic and flavonoid content in a flower sample purchased from the market were 12.3 ± 1.0 (mg gallic acid g−1 DW) and 9.6 ± 0.0 (mg quercetin g−1 DW), respectively, while the antioxidant activity of the 50% ethanol extracts of the flower sample was low (EC50 > 125 µg mL−1). All samples of C. indicum leaves showed higher total phenolic content and antioxidant activity than the flower sample. Therefore, C. indicum leaves have great potential as a material with high polyphenol content and high antioxidant activities.
Among the identified compounds, chlorogenic acid, 1,5-dicaffeoylquinic acid, 3,5-dicaffeoylquinic acid, and 4,5-dicaffeoylquinic acid have been reported to exhibit high antioxidant activity owing to the presence of a catechol structure [4]. Catechols have two hydroxyl groups at the ortho-position and can easily donate hydrogen atoms to radical substances to form o-quinones [35].
Statistical analysis using linear correlation showed that the antioxidant activity of the sample was correlated with 1,5-dicaffeoylquinic acid (R2 = 0.62) and 3,5-dicaffeoylquinic acid (R2 = 0.77), as shown in Figure 5, which is in agreement with a previous study that found that dicaffeoylquinic acid isomers in C. morifolium greatly contribute to antioxidant activity [4].
Additionally, the variations in the content of the identified antioxidant compounds were explained by the variations in exchangeable magnesium, cation exchangeable capacity, annual precipitation (BIO12), pH, annual mean temperature, precipitation, annual mean temperature (BIO01), clay, bulk density, and topographic wetness index (TWI) based on analysis of variance (ANOVA), as shown in Table 3. Magnesium deficiency in plants can cause a decrease in total phenolic content [36], and exchangeable magnesium is crucial for generating the main phenolic compounds responsible for the antioxidant activity of C. indicum. From this result, cation exchange capacity, which shows the ability of soil to hold and supply nutrients such as magnesium and calcium, is also an important parameter for plant growth, as mentioned in previous studies [37,38]. Moreover, this result suggests the significance of climatic factors, such as temperature and precipitation, on the phenolic content of C. indicum, as a previous study reported that temperature and precipitation influenced the secondary metabolite content [39].

2.3. Habitat Suitability Map (HSM) of MaxEnt Model

The construction of the Habitat Suitability Model (HSM) utilized the MaxEnt model, integrating 32 environmental layers as conditioning factors. These factors included a comprehensive range of topographical, soil, environmental, and climatic variables. The resulting model output was classified into four suitability classes—namely, low, moderate, high, and very high—employing the natural break classification technique within the ArcGIS platform, as shown in Figure 6.
The HSM employs statistical approaches to discriminate optimal plant zones by establishing correlations between conditioning factors and the occurrence of species within the study area. The assessment of the relative importance of variables influencing these conditioning factors is conducted through the rigorous jackknife variance estimation process, specifically applied to quantify the area under the curve (AUC). The outcomes of this analysis are depicted in Figure S1.
In the study area, conditioning factors such as soil sand content, TWI, hillshade, aspect, slope, and precipitation in the driest month (BIO14) had little impact on habitat suitability. In contrast, temperature, precipitation in the wettest month (BIO13) and BIO12, pH in H2O, organic content density (OCD), and digital elevation model (DEM) were the most important factors in predicting suitable areas for C. indicum growth in the study area. The samples were collected during the water-rich summer season; therefore, precipitation did not substantially affect habitat estimation, although temperature had a significant effect on suitability estimation.
The study area covered an elevation range of 56 to 325 m above sea level, where the DEM and elevation were identified as key factors influencing habitat suitability for C. indicum. Particularly, despite considering DEM-derived topographical factors, their impact on predicting habitat suitability for C. indicum was relatively modest or inconspicuous. This underscores the intricate relationship between topographical variables and the habitat preferences of C. indicum within the studied altitude range. Among the soil factors, the OCD and pH in H2O had important effects on the prediction of suitability for C. indicum. EC, soil bearing capacity, and pH influenced the suitability estimation, followed by soil bulk density, CEC, SM150T output (soil moisture sensor), and organic carbon content (OCC).
In accordance with the HSM, areas classified as very highly and highly suitable covered 3.2% and 10.1% of the total study area, respectively. The moderately suitable zone comprised a substantial proportion, accounting for 44.2% of the total area, while the low suitability zone encompassed 42.4% of the entire study area, as demonstrated in Figure 7.
From the HSM, all collected samples belonged to high and very high suitability zones, indicating that the study area had good growth conditions for C. indicum. This suitability model assists in plant cultivation and preserves habitats with existing plants in the study area.

2.4. Validation of MaxEnt Model

The validation of the HSM for C. indicum involved a comprehensive analysis using a receiver operating characteristic (ROC) curve within the MaxEnt modeling framework. Model performance was assessed through the application of the area under the curve (AUC) metric, as shown in Figure 8. The validation results showed a significant AUC value of 96.6% for the MaxEnt model within the study area. Models with AUC > 0.5 were indicative of superior performance, and in this instance, the MaxEnt model exhibited both logical and acceptable AUC values. The model demonstrated commendable precision in predicting the habitat suitability for C. indicum within the selected study area.

3. Materials and Methods

3.1. Study Area

In this study, the southeastern part of Kochi Prefecture was selected to investigate the natural habitats of C. indicum, as shown in Figure 9. Mountains and plains cover this area, with long hours of sunshine and a mild climate throughout the year. The region has a warm and temperate humid subtropical climate with an annual rainfall of approximately 2316 mm. Rainfall was lowest in December (approximately 57 mm) and highest in June (approximately 316 mm). The annual average temperature is 16.4 °C. The highest temperature of 27.2 °C was in August, and the lowest of 5.9 °C was in January, according to the Meteorological Agency of Japan (2020) [40].
The study area, covering 702.4 km2, is geographically situated between 133°35′0″ and 134°5′0″ E longitude and 33°30′0″ and 33°50′0″ N latitude within the Kochi Prefecture. As illustrated in Figure 9, the geographical delineation is represented by the ALOS PALSAR DEM satellite image and was processed using ArcGIS version 10.5.1 (ESRI Japan Corporation, Tokyo, Japan).

3.2. Plant Materials and Data Collection

During field investigations in the study area, 114 locations of C. indicum were investigated from 2019 to 2020. Among the investigated locations, 22 were selected for plant collection based on plant maturity, abundance, and status (Figure S2). In June, ten specimens of C. indicum, each ranging from 60 to 90 cm in height and exhibiting no signs of diseases, were systematically collected during the fieldwork. Geospatial data, including latitude, longitude, and altitude, were meticulously recorded at the collection sites using a global positioning system (GPS) device (Garmin eTrex 30x from Olathe, KS, USA). The flower of C. indicum was purchased from Tochimoto-Tenkaido (Osaka, Japan). A comparison was made between the flower sample and the collected leaf samples in terms of their chemical content and antioxidative activity. The leaf samples were dried in an air-circulating oven at 50 °C for 3 days. After that, the flower and leaf samples were ground into powder separately for extract preparation.
Geographical and climatic information on collected locations of C. indicum are shown in Table 4.

3.3. Chemicals

Imtakt CD-C18 was bought from Imtakt (Kyoto, Japan). Folin–Ciocalteu phenol reagent and DPPH were obtained from Sigma-Aldrich (St. Louis, MO, USA). 3,5-dicaffeoylquinic acid was purchased from MedChem Express (Monmouth Junction, NJ, USA). 2-Morpholinoethanesulfonic acid monohydrate (MES) was obtained from Dojin Chemical Research Institute (Kumamoto, Japan). Gallic acid was bought from the Tokyo Chemical Industry (TCI) (Tokyo, Japan). Anhydrous sodium carbonate and aluminum chloride (III) were purchased from Kishida Chemical Co., Ltd. (Osaka, Japan). Quercetin was bought from Funakoshi Co., Ltd. (Tokyo, Japan). Trolox, chlorogenic acid, acetonitrile, and other reagents were purchased from FUJIFILM Wako Pure Chemical Industries (Osaka, Japan).

3.4. Instrumentation

The EYELA multi-shaker MMS-type and centrifugal evaporators were obtained from Tokyo Rika Kikai (Tokyo, Japan) and Sakuma Seisakusho Co., Ltd. (Tokyo, Japan). A SpectraMax iD5 microplate reader (Molecular Devices, San Jose, CA, USA) was used for the in vitro assays. A Nexera X2 HPLC/UHPLC system (Shimadzu, Kyoto, Japan) was used for quantitative chromatographic analysis. A Bruker amaZon Speed Ion-Trap Mass Spectrometer (Billerica, MA, USA) was used for chemical identification. For the in situ soil analysis, a soil moisture sensor kit SM150T (Delta-T Devices, Cambridge, UK) and a Yamanaka-type soil hardness tester (Fujiwara Seisakusho, Tokyo, Japan) were used. Soil macronutrients were analyzed using an EW-THA1J soil analyzer (Air Water Biodesign, Osaka, Japan).

3.5. Soil Sample Collection and Analysis

Soil samples were acquired for subsequent chemical analyses. The soil samples were subjected to laboratory drying and subsequent chemical analysis using a soil analyzer (EW-THA1J, Air Water Biodesign Co., Ltd., Osaka, Japan).
The pH and soil bearing capacity values of each plant were collected at depths of 10, 25, 50, 75, and 100 cm, and the mean values were calculated for further analysis. An SM150T sensor was used to measure the dielectric properties of moist soil. The output of the SM150T was recorded in Volts (V). The soil refractive index (√ε) was determined by applying the polynomial Equation (1). The relationship between soil moisture (θ) and √ε was found to be linear (2). In the case of mineral soils, θ was calculated based on the output of the SM150T device using Equation (3) as per manufacturer’s recommendation [41].
√ε = 1.0 + 14.4396V − 31.2587V2 + 49.0575V3 − 36.5575V4 + 10.7117V5
√ε = a0 + a1 × θ
where a0 is the intercept, and a1 is the slope (mineral soil: a0 = 1.6, a0 = 8.4).
θmineral = −0.0714 + 1.7190V − 3.7213V2 + 5.8402V3 − 4.3521V4 + 1.275V5

3.6. Preparation of Plant Extract and Samples

About 1 g of dry plant powder was mixed with 40 mL of 50% ethanol and sonicated at 50 °C for 1 h. Subsequently, it was extracted at 25 °C for 1 d with shaking at 100 rpm. The extract was filtered through filter paper no. 2 (Advantec, Tokyo, Japan) and concentrated using a centrifuge evaporator. Finally, it was dried overnight under high vacuum (<10 Pa).
The dried extracts were dissolved in dimethyl sulfoxide (DMSO) at 0.5% and used to evaluate DPPH radical scavenging and polyphenolic and flavonoid contents.

3.7. Evaluation of DPPH Radical Scavenging Activity

The DPPH radical scavenging activity of the extracts was determined using the method described by Oki et al. with modifications [42]. Briefly, 100 µL of the dilution series of samples and positive control (Trolox) in 50% ethanol were transferred to the 96-well assay plate with 50 µL of MES buffer (0.2 M, pH 6.0) and 50 µL of DPPH (800 µM). For the negative control experiment, 100 µL of 50% ethanol was used instead of a sample solution. The mixture was allowed to stand in the dark at 25 °C for 20 min, and the absorbance was measured at 520 nm. Measurements were performed in triplicate, and the DPPH radical scavenging activity was calculated using Equation (4):
DPPH radical scavenging activity% = [(A0 − A1)/A0] × 100
where A0 is the absorbance of the negative control and A1 refers to the sample absorbance.
The EC50 for DPPH radical scavenging was calculated as previously described [43].

3.8. Total Phenolic Content

The total phenolic content of the extracts was determined using the Folin–Ciocalteu method described by Julkunen-Titto et al. [44]. In a 96-well plate, 20 μL of plant extract and 100 μL of 10% Folin–Ciocalteu’s phenol reagent were added, and the plate was allowed to stand for 5 min. Next, 80 μL of 2.5% sodium bicarbonate solution was added and allowed to stand for 1 h at 25 °C. Absorbance was measured at 755 nm using a microplate reader. The total phenolic content was determined as gallic acid equivalents (mg g−1 extract) based on a standard curve constructed with gallic acid.

3.9. Total Flavonoid Content

The total flavonoid content in the extracts was determined using the reported method with some changes [45]. In total, 25 μL of plant extract and 75 μL of ethanol were added to a 96-well plate. After that, 5 μL of 10% aluminum chloride (III) and 145 μL of pure water were filled. The 96-well plate was stored in the dark at 25 °C for 30 min. Absorbance was measured at 420 nm using a microplate reader. The total flavonoid content was determined as quercetin equivalent (mg g−1 extract) according to the standard curve made with quercetin.

3.10. Chromatographic Quantification and Analysis

A 50% ethanol dry extract of C. indicum was dissolved in 50% methanol to prepare a 0.5% solution. After centrifugation at 16,000 rpm for 10 min, the supernatant was collected with HPLC using an Imtakt CD-C18 column (2 × 150 mm, 3 µm). The mobile phase comprised 10% acetonitrile containing 0.1% acetic acid (A) and acetonitrile (B). The mixing ratio (A:B) was linearly changed from 100:0 to 53.3:46.7 for 35 min after the start, and then the column was washed with 100% B solution for 9 min. The injection volume was 1.0 µL, flow rate was 0.3 mL min−1, column temperature was 40 °C, and PDA detection wavelength was 200−360 nm.
Chlorogenic acid, 3,5-dicaffeoylquinic acid, and 4,5-dicaffeoylquinic acid were dissolved in 50% ethanol at 10, 30, 100, 300, and 1000 µg mL−1 to obtain the corresponding HPLC standard curves (R2 > 0.99). 1,5-dicaffeoylquinic acid was similarly dissolved in 50% ethanol up to 2000 µg mL−1. Acacetin was dissolved in DMSO–water (8:2) to obtain the corresponding HPLC standard curve (R2 > 0.99). Acacetin 7-O-glucoside and acacetin 7-O-malonylglucoside contents were also determined as acacetin equivalents, based on a standard curve of acacetin.
For the LC-MS analysis, 0.1% extract samples were analyzed under the same elution conditions as described above.

3.11. MaxEnt Model

MaxEnt software version 3.4.1 was downloaded from https://biodiversityinformatics.amnh.org/open_source/maxent/ (accessed on 21 December 2021) and used to predict habitat suitability for C. indicum.
In this investigation, 70% of the georeferenced dataset pertaining to C. indicum species was allocated for training, while the remaining 30% served for model validation. The collected species presence-only data, along with environmental layers in both continuous and categorical formats, were input into the MaxEnt model. The model ran with ten replicates to predict the distribution of C. indicum, and the relative importance of conditioning factors was systematically assessed using the jackknife test. The resulting output from the MaxEnt model delineates the potential distribution of C. indicum within the study area. Subsequently, the ASCII file (.ASC) was imported into ArcGIS 10.1 software (licensed) to generate a spatial distribution map of C. indicum. Higher values on the map signify an elevated likelihood of species presence, while lower values indicate a diminished degree of species adaptation.

3.11.1. Dataset Preparation for Habitat Suitability Modeling

To develop an efficient and suitable species habitat model, various climatic, environmental, topographical, and soil layers in the study area are required [46]. Therefore, 32 factors that influenced the habitat suitability for C. indicum were identified in the study area. The selected conditioning factors were classified as topographical, soil, environmental, or climatic. Topographical factors were extracted from an Advanced Land Observing Satellite-Digital Elevation Model (ALOS-DEM) satellite image with a 12.5 × 12.5 spatial resolution. This image was downloaded from the ALOS PALSAR (The Phased Array Type L-band Synthetic Aperture Radar) satellite website (https://search.asf.alaska.edu/#/ (accessed on 17 July 2022)). In addition, soil factors except for SM150T output and soil bearing capacity were downloaded from SoilGrids250m 2.0 (http://soilgrids.org (accessed on 11 May 2022)) at a depth of 0.15 m. Climatic factors, including annual mean rainfall and temperature were obtained from the Japan Meteorological Agency (https://www.jma.go.jp/ (accessed on 5 May 2022)). Additionally, 6 bioclimatic variables out of 19 [47], including BIO01, BIO05, BIO06, BIO12, BIO13, and BIO14, for 2010–2018 were downloaded from the WorldClim website (https://www.worldclim.org/ (accessed on 3 May 2022)). The environmental factors were prepared using topographical maps at a scale of 1:250,000 (Table S2).
Thematic maps of soil factors, climatic data, and environmental layers were prepared using the interpolation method of the inverse distance-weighted (IDW) spatial analysis technique in ArcGIS 10.1 version. All thematic raster layers were converted into the ASCII format with the same spatial resolution. All layers were then imported into the MaxEnt Version 3.4.1 software tool as environmental layers to predict habitat suitability for C. indicum species. The thematic maps and various conditioning factors associated with each plant location, used in predicting the habitat suitability of the study area, are illustrated in Figure S3 and detailed in Table S2. Table S3 presents a comprehensive list of conditioning factors, including their categories and data scales.

3.11.2. Validation of HSM

The validation process for the HSM involved employing the ROC curve to assess the predicted HSMs generated by the MaxEnt model within the study area. The AUC served as a metric to gauge the prediction accuracy of the MaxEnt model. Importantly, AUC values ranging from 0.5 to 1 were observed, with higher values indicative of more precise results. Within the ROC method, the x-axis portrays the cumulative percentage of suitability classes, while the y-axis represents the cumulative percentage of the training set within those classes [48]. It is important to highlight that the ROC method is essential for accurately modeling the estimated distribution of plant species [49].

3.12. Statistical Tests

Mean values and standard deviations were calculated for all measurements. The linear correlation between the examined variables, plant activity, and secondary metabolite content was calculated using the Data Analysis function of the Excel Add-in and XLSTAT statistical and data analysis solutions (Addinsoft, 2020, New York, NY, USA). Correlations with a p < 0.05 were considered significant.

4. Conclusions

In the present study, C. indicum leaves were revealed to have a potent antioxidant activity, which was possibly explained by the high content of chlorogenic acid and dicaffeoylquinic acid isomers. The leaves had little or no flavonoid compounds except for acacetin and the derivatives. This is in contrast with the flavonoid composition of the flower part, which was reported to contain various flavonoids including luteolin, apigenin, and acacetin, mainly in their glycosides.
The content of chlorogenic acid and dicaffeoylquinic acid isomers was significantly affected by exchangeable Mg and cation exchange capacity as well as BIO12, pH, annual mean temperature, precipitation, BIO01, clay, bulk density, and TWI. The MaxEnt model provided accurate HSM results for assessing the habitat suitability for C. indicum. Overall, climatic factors such as temperature, BIO13, and BIO12 had the strongest effect on habitat suitability. Soil factors including pH in H2O and OCD also significantly affected habitat predictions. These findings could be used for the future cultivation and investigation of C. indicum to obtain high phenolic content and bioactivity to meet industrial needs.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules29050927/s1. Table S1: Soil data of the plant collection sites. Table S2: Various conditioning factors were used to predict habitat suitability maps using the MaxEnt model. Table S3: Conditioning factors used to predict habitat suitability for C. indicum. Figure S1: Analysis of the relative importance of effective environmental variables. Figure S2: Photographs of Chrysanthemum indium L. in situ (at the collection sites). Figure S3: Thematic maps for predicting habitat suitability using the MaxEnt model for the study area.

Author Contributions

Conceptualization, R.U., R.A., D.H.M.A. and T.W.; formal analysis, R.U. and R.A.; investigation, R.U., R.A., D.W., H.O., K.S. and T.W.; methodology, R.U., R.A., D.H.M.A., K.K. and T.W.; software, R.A. and D.H.M.A.; validation, R.A., D.H.M.A. and K.K.; writing—original draft preparation, R.U., R.A. and D.H.M.A.; writing—review and editing, R.U., R.A., D.H.M.A. and K.K. All authors have read and agreed to the published version of the manuscript.

Funding

This study received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article and supplementary materials.

Acknowledgments

The authors would like to acknowledge the support of the School of Pharmacy and the Global Center for Natural Resources Sciences, Kumamoto University, Japan, for providing laboratory facilities to conduct the experiments.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Hadizadeh, H.; Samiei, L.; Shakeri, A. Chrysanthemum, an Ornamental Genus with Considerable Medicinal Value: A Comprehensive Review. S. Afr. J. Bot. 2022, 144, 23–43. [Google Scholar] [CrossRef]
  2. Yang, W.S.; Kim, D.; Yi, Y.-S.; Kim, J.H.; Jeong, H.Y.; Hwang, K.; Kim, J.-H.; Park, J.; Cho, J.Y. AKT-Targeted Anti-Inflammatory Activity of the Methanol Extract of Chrysanthemum indicum Var. Albescens. J. Ethnopharmacol. 2017, 201, 82–90. [Google Scholar] [CrossRef]
  3. Society of Japanese Pharmacopoeia. The Japanese Pharmacopoeia, 18th ed.; Pharmaceuticals and Medical Devices Agency: Tokyo, Japan, 2021. [Google Scholar]
  4. Gong, J.; Chu, B.; Gong, L.; Fang, Z.; Zhang, X.; Qiu, S.; Wang, J.; Xiang, Y.; Xiao, G.; Yuan, H.; et al. Comparison of Phenolic Compounds and the Antioxidant Activities of Fifteen Chrysanthemum morifolium Ramat Cv. ‘Hangbaiju’ in China. Antioxidants 2019, 8, 325. [Google Scholar] [CrossRef]
  5. Lee, D.Y.; Choi, G.; Yoon, T.; Cheon, M.S.; Choo, B.K.; Kim, H.K. Anti-Inflammatory Activity of Chrysanthemum indicum Extract in Acute and Chronic Cutaneous Inflammation. J. Ethnopharmacol. 2009, 123, 149–154. [Google Scholar] [CrossRef]
  6. Jeong, S.C.; Kim, S.M.; Jeong, Y.T.; Song, C.H. Hepatoprotective Effect of Water Extract from Chrysanthemum indicum L. Flower. Chin. Med. 2013, 8, 7. [Google Scholar] [CrossRef]
  7. Chen, M.; Wang, K.; Zhang, Y.; Zhang, M.; Ma, Y.; Sun, H.; Jin, Z.; Zheng, H.; Jiang, H.; Yu, P.; et al. New Insights into the Biological Activities of Chrysanthemum morifolium: Natural Flavonoids Alleviate Diabetes by Targeting α-Glucosidase and the PTP-1B Signaling Pathway. Eur. J. Med. Chem. 2019, 178, 108–115. [Google Scholar] [CrossRef]
  8. Ma, D.; Wako, Y. Evaluation of Phenolic Compounds and Neurotrophic/Neuroprotective Activity of Cultivar Extracts Derived from Chrysanthemum morifolium Flowers. Food Sci. Technol. Res. 2017, 23, 457–467. [Google Scholar] [CrossRef]
  9. Choi, K.T.; Kim, J.H.; Cho, H.T.; Lim, S.S.; Kwak, S.S.; Kim, Y.J. Dermatologic Evaluation of Cosmetic Formulations Containing Chrysanthemum indicum Extract. J. Cosmet. Dermatol. 2016, 15, 162–168. [Google Scholar] [CrossRef] [PubMed]
  10. Yamada, K.; Sawa, T.; Murayama, K.; Yamaguchi, S.; Miyazaki, Y.; Tachibana, H. Histamine Release-Suppressive Effect of Water Extracts Prepared from Flower and Leaf of Chrysanthemum, Shiranui Himekiku. Nippon. Shokuhin Kagaku Kogaku Kaishi 2012, 59, 394–400. [Google Scholar] [CrossRef]
  11. Shim, S.-Y.; Kang, H.-S.; Sun, H.-J.; Lee, Y.-J.; Park, J.-R.; Chun, S.-S.; Song, Y.-H.; Byun, D.-S. Isolation and Identification of Flavonoids from Gujeolcho (Chrysanthemum zawadskii Var. latilobum) as Inhibitor of Histamine Release. Food Sci. Biotechnol. 2012, 21, 613–617. [Google Scholar] [CrossRef]
  12. Lai, J.-P.; Lim, Y.H.; Su, J.; Shen, H.-M.; Ong, C.N. Identification and Characterization of Major Flavonoids and Caffeoylquinic Acids in Three Compositae Plants by LC/DAD-APCI/MS. J. Chromatogr. B 2007, 848, 215–225. [Google Scholar] [CrossRef]
  13. Jiang, S.; Wang, M.; Jiang, Z.; Zafar, S.; Xie, Q.; Yang, Y.; Liu, Y.; Yuan, H.; Jian, Y.; Wang, W. Chemistry and Pharmacological Activity of Sesquiterpenoids from the Chrysanthemum Genus. Molecules 2021, 26, 3038. [Google Scholar] [CrossRef]
  14. Liu, C.-C.; Zhang, Y.; Dai, B.-L.; Ma, Y.-J.; Zhang, Q.; Wang, Y.; Yang, H. Chlorogenic Acid Prevents Inflammatory Responses in IL-1β-Stimulated Human SW-1353 Chondrocytes, a Model for Osteoarthritis. Mol. Med. Rep. 2017, 16, 1369–1375. [Google Scholar] [CrossRef]
  15. Hunyadi, A.; Martins, A.; Hsieh, T.-J.; Seres, A.; Zupkó, I. Chlorogenic Acid and Rutin Play a Major Role in the In Vivo Anti-Diabetic Activity of Morus alba Leaf Extract on Type II Diabetic Rats. PLoS ONE 2012, 7, e50619. [Google Scholar] [CrossRef]
  16. Hong, S.; Joo, T.; Jhoo, J.-W. Antioxidant and Anti-Inflammatory Activities of 3,5-Dicaffeoylquinic Acid Isolated from Ligularia fischeri Leaves. Food Sci. Biotechnol. 2015, 24, 257–263. [Google Scholar] [CrossRef]
  17. Yao, H.; Shang, Z.; Wang, P.; Li, S.; Zhang, Q.; Tian, H.; Ren, D.; Han, X. Protection of Luteolin-7-O-Glucoside Against Doxorubicin-Induced Injury Through PTEN/Akt and ERK Pathway in H9c2 Cells. Cardiovasc. Toxicol. 2016, 16, 101–110. [Google Scholar] [CrossRef] [PubMed]
  18. Guo, X.; Cao, W.; Yao, J.; Yuan, Y.; Hong, Y.; Wang, X.; Xing, J. Cardioprotective Effects of Tilianin in Rat Myocardial Ischemia-Reperfusion Injury. Mol. Med. Rep. 2015, 11, 2227–2233. [Google Scholar] [CrossRef] [PubMed]
  19. De Stefano, A.; Caporali, S.; Di Daniele, N.; Rovella, V.; Cardillo, C.; Schinzari, F.; Minieri, M.; Pieri, M.; Candi, E.; Bernardini, S.; et al. Anti-Inflammatory and Proliferative Properties of Luteolin-7-O-Glucoside. Int. J. Mol. Sci. 2021, 22, 1321. [Google Scholar] [CrossRef] [PubMed]
  20. Fuchs, J.; Milbradt, R. Skin Anti-Inflammatory Activity of Apigenin-7-Glucoside in Rats. Arzneimittelforschung 1993, 43, 370–372. [Google Scholar] [PubMed]
  21. Smiljkovic, M.; Stanisavljevic, D.; Stojkovic, D.; Petrovic, I.; Marjanovic Vicentic, J.; Popovic, J.; Golic Grdadolnik, S.; Markovic, D.; Sankovic-Babice, S.; Glamoclija, J.; et al. Apigenin-7-O-Glucoside versus Apigenin: Insight into the Modes of Anticandidal and Cytotoxic Actions. EXCLI J. 2017, 16, 795–807. [Google Scholar] [CrossRef]
  22. Singh, R.P. Acacetin Inhibits Cell Growth and Cell Cycle Progression, and Induces Apoptosis in Human Prostate Cancer Cells: Structure-Activity Relationship with Linarin and Linarin Acetate. Carcinogenesis 2005, 26, 845–854. [Google Scholar] [CrossRef]
  23. Rajalakshmi, G.; Komathi, S.; Raviganesh, B.; Poongodi, N.; Sasikala, T. In-Vitro Micropropagation and Antimicrobial Activity of Chrysanthemum indicum. Sch. Acad. J. Pharm. 2013, 2, 285–288. [Google Scholar]
  24. Li, Y.; Liu, X.-J.; Su, S.-L.; Yan, H.; Guo, S.; Qian, D.-W.; Duan, J.-A. Evaluation of Anti-Inflammatory and Antioxidant Effects of Chrysanthemum Stem and Leaf Extract on Zebrafish Inflammatory Bowel Disease Model. Molecules 2022, 27, 2114. [Google Scholar] [CrossRef]
  25. Kim, Y.-S.; Hwang, J.-W.; Park, P.-J.; Jeong, J.-H. Antioxidant Activity and Protective Effects of Extracts from Chrysanthemum boreale on t-BHP Induced Oxidative Stress in Chang Cells. J. Korean Soc. Food Sci. Nutr. 2014, 43, 60–66. [Google Scholar] [CrossRef]
  26. Wang, T.; Shen, X.-G.; Guo, Q.-S.; Zhou, J.-S.; Mao, P.-F.; Shen, Z.-G. Comparison of Major Bioactive Components from Leaves of Chrysanthemum morifolium. China J. Chin. Mater. Medica 2015, 40, 1670–1675. [Google Scholar] [CrossRef]
  27. Alsaadi, D.H.M.; Raju, A.; Kusakari, K.; Karahan, F.; Sekeroglu, N.; Watanabe, T. Phytochemical Analysis and Habitat Suitability Mapping of Glycyrrhiza glabra L. Collected in the Hatay Region of Turkey. Molecules 2020, 25, 5529. [Google Scholar] [CrossRef]
  28. Hodaei, M.; Rahimmalek, M.; Arzani, A.; Talebi, M. The Effect of Water Stress on Phytochemical Accumulation, Bioactive Compounds and Expression of Key Genes Involved in Flavonoid Biosynthesis in Chrysanthemum morifolium L. Ind. Crops Prod. 2018, 120, 295–304. [Google Scholar] [CrossRef]
  29. Zheng, L.; Van Labeke, M.-C. Chrysanthemum Morphology, Photosynthetic Efficiency and Antioxidant Capacity Are Differentially Modified by Light Quality. J. Plant Physiol. 2017, 213, 66–74. [Google Scholar] [CrossRef] [PubMed]
  30. Hirzel, A.; Guisan, A. Which Is the Optimal Sampling Strategy for Habitat Suitability Modelling. Ecol. Modell. 2002, 157, 331–341. [Google Scholar] [CrossRef]
  31. Phillips, S.J.; Anderson, R.P.; Schapire, R.E. Maximum Entropy Modeling of Species Geographic Distributions. Ecol. Model. 2006, 190, 231–259. [Google Scholar] [CrossRef]
  32. Convertino, M.; Muñoz-Carpena, R.; Chu-Agor, M.L.; Kiker, G.A.; Linkov, I. Untangling Drivers of Species Distributions: Global Sensitivity and Uncertainty Analyses of MaxEnt. Environ. Model. Softw. 2014, 51, 296–309. [Google Scholar] [CrossRef]
  33. Han, A.-R.; Nam, B.; Kim, B.-R.; Lee, K.-C.; Song, B.-S.; Kim, S.; Kim, J.-B.; Jin, C. Phytochemical Composition and Antioxidant Activities of Two Different Color Chrysanthemum Flower Teas. Molecules 2019, 24, 329. [Google Scholar] [CrossRef] [PubMed]
  34. Sugawara, T.; Igarashi, K. Identification of Major Flavonoids in Petals of Edible Chrysanthemum Flowers and Their Suppressive Effect on Carbon Tetrachloride-Induced Liver Injury in Mice. Food Sci. Technol. Res. 2009, 15, 499–506. [Google Scholar] [CrossRef]
  35. Dey, S.K.; Mukherjee, A. Catechol Oxidase and Phenoxazinone Synthase: Biomimetic Functional Models and Mechanistic Studies. Coord. Chem. Rev. 2016, 310, 80–115. [Google Scholar] [CrossRef]
  36. Li, J.; Li, Q.-H.; Zhang, X.-Y.; Zhang, L.-Y.; Zhao, P.-L.; Wen, T.; Zhang, J.-Q.; Xu, W.-L.; Guo, F.; Zhao, H.; et al. Exploring the Effects of Magnesium Deficiency on the Quality Constituents of Hydroponic-Cultivated Tea (Camellia sinensis L.) Leaves. J. Agric. Food Chem. 2021, 69, 14278–14286. [Google Scholar] [CrossRef] [PubMed]
  37. Manrique, L.A.; Jones, C.A.; Dyke, P.T. Predicting Cation-Exchange Capacity from Soil Physical and Chemical Properties. Soil. Sci. Soc. Am. J. 1991, 55, 787–794. [Google Scholar] [CrossRef]
  38. Cui, X.; Mao, P.; Sun, S.; Huang, R.; Fan, Y.; Li, Y.; Li, Y.; Zhuang, P.; Li, Z. Phytoremediation of Cadmium Contaminated Soils by Amaranthus hypochondriacus L.: The Effects of Soil Properties Highlighting Cation Exchange Capacity. Chemosphere 2021, 283, 131067. [Google Scholar] [CrossRef]
  39. Sampaio, B.L.; Edrada-Ebel, R.; Da Costa, F.B. Effect of the Environment on the Secondary Metabolic Profile of Tithonia diversifolia: A Model for Environmental Metabolomics of Plants. Sci. Rep. 2016, 6, 29265. [Google Scholar] [CrossRef]
  40. Japan Meteorological Agency. General Information on Climate of Japan. Available online: https://www.jma.go.jp/jma/indexe.html (accessed on 10 April 2021).
  41. Delta-T Devices Ltd. User Manual for the SM150T-UM-0.f. Soil Moisture Sensor; Delta-T Devices Ltd.: Cambridge, UK, 2016. [Google Scholar]
  42. Oki, T.; Masuda, M.; Furuta, S.; NishibaI, Y.; Suda, I. Radical Scavenging Activity of Fried Chips Made from Purple-Fleshed Sweet Potato. Nippon Shokuhin Kagaku Kogaku Kaishi 2001, 48, 926–932. [Google Scholar] [CrossRef]
  43. Zar Wynn Myint, K.; Kido, T.; Kusakari, K.; Prasad Devkota, H.; Kawahara, T.; Watanabe, T. Rhusflavanone and Mesuaferrone B: Tyrosinase and Elastase Inhibitory Biflavonoids Extracted from the Stamens of Mesua ferrea L. Nat. Prod. Res. 2021, 35, 1024–1028. [Google Scholar] [CrossRef]
  44. Julkunen-Tiitto, R. Phenolic Constituents in the Leaves of Northern Willows: Methods for the Analysis of Certain Phenolics. J. Agric. Food Chem. 1985, 33, 213–217. [Google Scholar] [CrossRef]
  45. Pękal, A.; Pyrzynska, K. Evaluation of Aluminium Complexation Reaction for Flavonoid Content Assay. Food Anal. Methods 2014, 7, 1776–1782. [Google Scholar] [CrossRef]
  46. Zhang, K.; Zhang, Y.; Zhou, C.; Meng, J.; Sun, J.; Zhou, T.; Tao, J. Impact of Climate Factors on Future Distributions of Paeonia ostii across China Estimated by MaxEnt. Ecol. Inform. 2019, 50, 62–67. [Google Scholar] [CrossRef]
  47. Fick, S.E.; Hijmans, R.J. WorldClim 2: New 1-km Spatial Resolution Climate Surfaces for Global Land Areas. Int. J. Climatol. 2017, 37, 4302–4315. [Google Scholar] [CrossRef]
  48. Blahut, J.; van Westen, C.J.; Sterlacchini, S. Analysis of Landslide Inventories for Accurate Prediction of Debris-Flow Source Areas. Geomorphology 2010, 119, 36–51. [Google Scholar] [CrossRef]
  49. Zhang, X.; Li, X.; Feng, Y.; Liu, Z. The Use of ROC and AUC in the Validation of Objective Image Fusion Evaluation Metrics. Signal Process. 2015, 115, 38–48. [Google Scholar] [CrossRef]
Figure 1. Linear relationship between relative strength of DPPH radical scavenging (1/EC50) and the total phenolic and flavonoid content.
Figure 1. Linear relationship between relative strength of DPPH radical scavenging (1/EC50) and the total phenolic and flavonoid content.
Molecules 29 00927 g001
Figure 2. Typical HPLC chromatogram of 50% ethanol extract of C. indicum leaves showing peaks of chlorogenic acid (1), 1,5-dicaffeoylquinic acid (2), 3,5-dicaffeoylquinic acid (3), 4,5-dicaffeoylquinic acid (4), acacetin 7-O-glucoside (5), acacetin 7-O-malonylglucoside (6), and acacetin (7).
Figure 2. Typical HPLC chromatogram of 50% ethanol extract of C. indicum leaves showing peaks of chlorogenic acid (1), 1,5-dicaffeoylquinic acid (2), 3,5-dicaffeoylquinic acid (3), 4,5-dicaffeoylquinic acid (4), acacetin 7-O-glucoside (5), acacetin 7-O-malonylglucoside (6), and acacetin (7).
Molecules 29 00927 g002
Figure 3. Phenolic content of chlorogenic acid, 1,5-dicaffeoylquinic acid, 3,5-dicaffeoylquinic acid, and 4,5-dicaffeoylquinic acid in the collected leaves.
Figure 3. Phenolic content of chlorogenic acid, 1,5-dicaffeoylquinic acid, 3,5-dicaffeoylquinic acid, and 4,5-dicaffeoylquinic acid in the collected leaves.
Molecules 29 00927 g003
Figure 4. Flavonoid content of acacetin and its glucoside in the collected leaves.
Figure 4. Flavonoid content of acacetin and its glucoside in the collected leaves.
Molecules 29 00927 g004
Figure 5. Relationship between relative strength of DPPH radical scavenging (1/EC50) and dicaffeoylquinic acid isomers (1,5-dicaffeoylquinic acid and 3,5-dicaffeoylquinic acid).
Figure 5. Relationship between relative strength of DPPH radical scavenging (1/EC50) and dicaffeoylquinic acid isomers (1,5-dicaffeoylquinic acid and 3,5-dicaffeoylquinic acid).
Molecules 29 00927 g005
Figure 6. Habitat suitability map of C. indicum using MaxEnt model.
Figure 6. Habitat suitability map of C. indicum using MaxEnt model.
Molecules 29 00927 g006
Figure 7. Percentages of habitat suitability classes.
Figure 7. Percentages of habitat suitability classes.
Molecules 29 00927 g007
Figure 8. Receiver operator characteristic (ROC) curve with AUC value of C. indicum habitat for MaxEnt model.
Figure 8. Receiver operator characteristic (ROC) curve with AUC value of C. indicum habitat for MaxEnt model.
Molecules 29 00927 g008
Figure 9. Location map of the study area with the 22 collected plant samples. Kochi Prefecture of Japan shown on ALOS PALSAR 90 M DEM satellite image (https://search.asf.alaska.edu/#/ (accessed on 17 July 2022)).
Figure 9. Location map of the study area with the 22 collected plant samples. Kochi Prefecture of Japan shown on ALOS PALSAR 90 M DEM satellite image (https://search.asf.alaska.edu/#/ (accessed on 17 July 2022)).
Molecules 29 00927 g009
Table 1. Total phenolic and flavonoid content and DPPH radical scavenging activity of 50% ethanol extracts of C. indicum leaves *.
Table 1. Total phenolic and flavonoid content and DPPH radical scavenging activity of 50% ethanol extracts of C. indicum leaves *.
ID aTotal Phenolic Content
(mg Gallic Acid g−1 DW) b
Total Flavonoid Content
(mg Quercetin g−1 DW) b
EC50 for DPPH Radical Scavenging (µg mL−1) c
Ko-127.2 ± 0.16.2 ± 0.169.8 ± 1.5
Ko-217.5 ± 0.16.3 ± 0.692.2 ± 1.3
Ko-315.0 ± 0.72.3 ± 0.0123.2 ± 3.5
Ko-428.3 ± 0.15.5 ± 0.049.0 ± 3.9
Ko-518.9 ± 0.24.1 ± 0.161.3 ± 10.0
Ko-622.7 ± 0.28.0 ± 0.074.2 ± 2.7
Ko-723.6 ± 0.06.9 ± 0.060.8 ± 6.6
Ko-864.1 ± 0.26.4 ± 0.129.4 ± 4.7
Ko-938.6 ± 0.78.9 ± 0.139.8 ± 4.5
Ko-1015.6 ± 0.56.9 ± 0.265.9 ± 1.8
Ko-1148.6 ± 1.16.7 ± 0.031.9 ± 0.2
Ko-1261.9 ± 0.310.6 ± 0.132.7 ± 0.8
Ko-1321.9 ± 0.12.5 ± 0.257.7 ± 1.9
Ko-1455.8 ± 0.810.5 ± 0.034.1 ± 0.8
Ko-1540.1 ± 1.67.4 ± 0.138.0 ± 0.6
Ko-1653.0 ± 0.611.4 ± 0.132.2 ± 3.0
Ko-1741.5 ± 0.26.9 ± 0.128.0 ± 1.6
Ko-1856.4 ± 0.310.3 ± 0.130.7 ± 4.3
Ko-1953.1 ± 0.69.1 ± 0.036.2 ± 7.6
Ko-2030.0 ± 0.15.9 ± 0.138.3 ± 1.2
Ko-2117.5 ± 0.04.2 ± 0.076.0 ± 1.3
Ko-2215.2 ± 0.23.2 ± 0.076.5 ± 2.2
* Data represent the average ± standard deviation from triplicate measurements of a mixed sample collected in the same location. a ID represents C. indicum collection sites, where Ko refers to Kochi Prefecture. b DW represents dry weight of the leaves sample of C. indicum. c EC50 represents effective concentration for 50% scavenging of DPPH radicals.
Table 2. Identified compounds in LC-MS.
Table 2. Identified compounds in LC-MS.
Peak No.Identified CompoundObserved Ion (Ion–Trap MS)Retention Time (min.)
1Chlorogenic acidm/z [M + H]+: 355.0, m/z [M − H]: 353.14.1
21,5-dicaffeoylquinic acidm/z [M + H]+: 517.1, m/z [M − H]: 515.210.7
33,5-dicaffeoylquinic acidm/z [M − H]: 515.111.1
44,5-dicaffeoylquinic acidm/z [M + H]+: 517.1, m/z [M − H]: 515.112.1
5Acacetin 7-O-glucosidem/z [M + H]+: 447.0, m/z [M + H − 162]+: 285.018.4
6Acacetin 7-O-malonylglucosidem/z [M + H]+: 533.1, m/z [M + H − 248]+: 285.020.0
7Acacetinm/z [M + H]+: 285.0, m/z [M − H]: 283.029.2
Table 3. ANOVA analysis of variables effects on chlorogenic acid, 1,5-dicaffeoylquinic acid, 3,5-dicaffeoylquinic acid, 4,5-dicaffeoylquinic acid, and acacetin 7-O-malonylglucoside.
Table 3. ANOVA analysis of variables effects on chlorogenic acid, 1,5-dicaffeoylquinic acid, 3,5-dicaffeoylquinic acid, 4,5-dicaffeoylquinic acid, and acacetin 7-O-malonylglucoside.
VariablesCGA
(%)
1,5-DCQA
(%)
3,5-DCQA
(%)
4,5-DCQA
(%)
Acacetin
7-O-Malonylglucoside (%)
F ValueR2F ValueR2F ValueR2F ValueR2F ValueR2
Exchangeable
magnesium
2.250.854.12 *0.9117.87 ***0.9890.23 ***0.991.210.75
Cation exchangeable
capacity
3.61 *0.695.01 *0.762.98 *0.6516.61 ***0.912.050.65
BIO121.260.491.310.502.660.6713.21 ***0.911.090.45
pH0.700.481.790.711.48 ***0.942.070.734.73 *0.86
Annual mean temperature and precipitation1.220.573.11 *0.775.34 *0.869.73 **0.921.470.62
BIO0110.532.090.704.74 *0.849.65 **0.910.880.49
Clay0.660.171.130.262.910.489.09 ***0.741.130.26
Bulk density2.210.884.250.938.70 *0.974.030.930.870.74
TWI2.150.904.200.957.34 *0.973.040.934.540.95
* p < 0.05, ** p < 0.01, and *** p < 0.001.
Table 4. Geographical and climatic information of collection sites of C. indicum.
Table 4. Geographical and climatic information of collection sites of C. indicum.
ID aLongitude (E)Latitude (N)Elevation (m)Mean Annual
Temperature (°C) b
Mean Annual Rainfall (mm) b
Ko-1133°43′50.4″33°35′33.12″32615.92059
Ko-2133°43′35.5″33°35′20.82″29616.32054
Ko-3133°43′14.59″33°34′48.65″19716.52055
Ko-4133°43′14.34″33°34′48.1″19916.52055
Ko-5133°42′58.3″33°34′38.06″22616.52055
Ko-6133°42′59.08″33°34′35.9″22016.52055
Ko-7133°42′57.57″33°34′38.28″23216.52055
Ko-8133°44′6.93″33°38′28.55″14115.72017
Ko-9133°43′55.26″33°38′36.82″19115.32042
Ko-10133°43′59.79″33°38′33.49″21415.22045
Ko-11133°44′21.54″33°38′46.85″30115.22045
Ko-12133°44′11.05″33°38′44.76″32215.22045
Ko-13133°43′49.03″33°38′19.7″8616.11995
Ko-14133°43′6.83″33°37′18.78″17616.62000
Ko-15133°52′20.23″33°40′50.78″27514.32049
Ko-16133°52′23.51″33°40′50.72″26314.32049
Ko-17133°52′18.18″33°40′5″27314.32049
Ko-18133°52′36.35″33°40′30.73″29213.92073
Ko-19133°52′36.33″33°40′30.54″29313.92073
Ko-20133°52′59.64″33°40′5.54″29614.02078
Ko-21133°44′24.09″33°38′32.58″9715.22045
Ko-22133°43′8.01″33°38′35.53″7515.62022
a ID represents the location of the C. indicum collection, where Ko refers to the Kochi Prefecture. b Acquired from Japan Meteorological Agency, General Information on Climate of Japan (https://www.jma.go.jp/ (accessed on 5 May 2022)).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Uranishi, R.; Aedla, R.; Alsaadi, D.H.M.; Wang, D.; Kusakari, K.; Osaki, H.; Sugimura, K.; Watanabe, T. Evaluation of Environmental Factor Effects on the Polyphenol and Flavonoid Content in the Leaves of Chrysanthemum indicum L. and Its Habitat Suitability Prediction Mapping. Molecules 2024, 29, 927. https://doi.org/10.3390/molecules29050927

AMA Style

Uranishi R, Aedla R, Alsaadi DHM, Wang D, Kusakari K, Osaki H, Sugimura K, Watanabe T. Evaluation of Environmental Factor Effects on the Polyphenol and Flavonoid Content in the Leaves of Chrysanthemum indicum L. and Its Habitat Suitability Prediction Mapping. Molecules. 2024; 29(5):927. https://doi.org/10.3390/molecules29050927

Chicago/Turabian Style

Uranishi, Rei, Raju Aedla, Doaa H. M. Alsaadi, Dongxing Wang, Ken Kusakari, Hirotaka Osaki, Koji Sugimura, and Takashi Watanabe. 2024. "Evaluation of Environmental Factor Effects on the Polyphenol and Flavonoid Content in the Leaves of Chrysanthemum indicum L. and Its Habitat Suitability Prediction Mapping" Molecules 29, no. 5: 927. https://doi.org/10.3390/molecules29050927

APA Style

Uranishi, R., Aedla, R., Alsaadi, D. H. M., Wang, D., Kusakari, K., Osaki, H., Sugimura, K., & Watanabe, T. (2024). Evaluation of Environmental Factor Effects on the Polyphenol and Flavonoid Content in the Leaves of Chrysanthemum indicum L. and Its Habitat Suitability Prediction Mapping. Molecules, 29(5), 927. https://doi.org/10.3390/molecules29050927

Article Metrics

Back to TopTop