Agro-Climatic Suitability Evaluation for Saffron Production in Areas of Western Himalaya

Saffron (Crocus sativus L.) is an expensive spice crop cultivated successfully in Iran, Afghanistan, India, Greece, Morocco, Spain, and Italy. The present study was conducted during the periods 2018–2019 and 2019–2020 to evaluate the morphological, yield and quality parameters of saffron in six different regions of non-traditional areas of the western Himalayas. The two experimental factors were “year” and “location.” The experiment was conducted using a factorial randomized block design with three replications. Yield attributes, viz., number of flowers, fresh flower yield, fresh stigma yield and dry stigma yield were significantly higher in location L₃ compared to other geographical locations. Dry stigma yield in location L₃ was higher by 50.0, 41.2, 33.3, 14.3, and 9.1% compared to locations L₆, L₅, L₁, L₄, and L₂, respectively. These were characterized by the appropriate climatic conditions, viz., high altitude, sandy-loam texture of the soil, optimum temperature, lesser relative humidity and total rainfall, demonstrating that it is possible to cultivate this spice even in non-traditional areas of the western Himalaya. Positive correlations were established for stigma yield with increased altitude and lesser rainfall. Secondary metabolites viz., crocin and picrocrocin increased significantly with the increase in altitude; however, a reverse trend was recorded for safranal content. Total phenolics and flavonoids were significantly higher in the geographical location of Kinnaur, H.P. (L₁) and Bharmour, H.P. (L₄). In conclusion, the assessment of different geographical locations and soil types is particularly necessary to encourage saffron production and its qualitative traits. Based on current findings, saffron can be grown successfully in some non-traditional locations of the western Himalayan regions.

Introduction

Saffron (Crocus sativus L.) “(family Iridaceae),” commonly known as red gold, is one of the most expensive and valuable spice crops in the world market (Gomez-Gomez et al., 2012). It is a triploid plant that is propagated through corms (Bayat et al., 2016). More than 418 tons/annum of saffron are produced worldwide from an area of 108,000 ha in Iran, 7,557 ha in Afghanistan, 3,674 ha in India, 1,000 ha in Greece, 850 ha in Morocco, 150 ha in Spain, 70 ha in Italy and 37 ha in France (Cardone et al., 2020a). Saffron flowers in the range of 75–100 provide 225–300 stigma threads, which produce only 0.5 g of dry stigma (Gohari et al., 2013). The crop price is very high and varies from 1,500 to 2,200 Euro/kg (Mykhailenko et al., 2020). Nowadays, saffron cultivation is gaining interest due to its major uses in industries, viz., the textile, dye, drug and culinary adjunct, food additive, coloring, and flavoring industries; furthermore, it has also gained interest for it’s pharmacological properties, e.g., antioxidant, antitumor activity, anticancer, and antimutagenic activity (Moradi Rikabad et al., 2019; Mykhailenko et al., 2020).

Altitude, soil characteristics, temperature, photoperiod, and topographical locations are the critical environmental parameters that affect saffron production (Siracusa et al., 2010; Rahimi et al., 2017; Cardone et al., 2020b). This spice crop grows well in friable, loose, low density, well-irrigated and well-drained clay calcareous soils with an optimum pH range between 6.8–7.8 and electrical conductivity (E.C.) below 2 dS m^–1 (Zarghani et al., 2016). It is cultivated mainly in Iran, Afghanistan, Morocco, India, Spain, Greece, and Italy (Cardone et al., 2019). Previous studies suggest that 600 mm seasonal rainfall is almost sufficient for saffron cultivation under rainfed conditions, which can vary depending on soil characteristics and fertilization practices (Fallahi and Mahmoodi, 2018a,b). Gresta et al. (2009) reported the highest flower number and dry stigma yield when corms were planted with high density in sandy soil, while the highest stigma weight was obtained on sowing corms in clay soil with high density.

The quality of the saffron entirely depends on the environmental conditions, content and composition of secondary metabolites. The major marker compounds crocins, picrocrocin and safranal are responsible for the coloring, bittering and aromatic powers, respectively (Siracusa et al., 2010; Garcia-Rodriguez et al., 2017; Giorgi et al., 2017; Mykhailenko et al., 2020). Besides the significant amount of information that exists on aspects of physiological, morphological, genotype, chemical, and genetic diversity, there remains a lack of information on different environmental conditions that alter the yield and secondary metabolite profile of the saffron crop. In India, it is only cultivated in traditional areas, viz., Pulwama and Kishtwar districts in Jammu and Kashmir, in silty clay loam textures with an electrical conductivity ranging from 0.09 to 0.30 dS m^–1, pH between 6.3 and 8.3, calcium carbonate content of 4.61% and average organic carbon of 0.35% (Ganaie and Singh, 2019). Although many studies have been performed, the adaptability of this crop under diverse environmental conditions for studying different vegetative, qualitative and quantitative traits needs to be further investigated under the western Himalayan conditions (Cardone et al., 2020b).

To meet the increased demand for saffron, an attempt was made in India in Himachal Pradesh and Uttarakhand with the following objectives: (i) to identify the suitable locations for the cultivation of saffron in different regions of western Himalayas and (ii) to test the quality of the harvested spice using ultra pressure liquid chromatography (UPLC).

Materials and Methods

Experimental Site Description and Planting Material

The experiment was conducted during the years 2018–2019 and 2019–2020 at six geographical locations of different altitudinal zones of Himachal Pradesh and Uttarakhand, characterized by different environmental conditions. Potential locations for the experiment were identified with the maximum entropy (MAXENT) model on a map and depicted in Figure 1. The study sites were Moorang, Kinnaur, Himachal Pradesh (L₁), Kapkote, Bageshwar, Uttarakhand (L₂), Suppa, Bharmour, Himachal Pradesh (L₃), Sathli, Bhamour, Himachal Pradesh (L₄), Langha, Palampur, Himachal Pradesh (L₅) and the Council of Scientific and Industrial Research-Institute of Himalayan Bioresource Technology, Palampur, Himachal Pradesh (L₆). Geographic coordinates (latitude and longitude), altitude and an average of 2 years’ standard meteorological data (average minimum and maximum temperatures, relative humidity and total rainfall) during cultivation seasons from each site were extrapolated from meteorological stations nearest each sample site. Saffron corms were procured from Kishtwar, Jammu and Kashmir, India. The average size of each saffron corm was 10–15 g.

FIGURE 1

Figure 1. Identification of different locations through the MAXENT model.

Agronomic Management

At each location, the land was prepared by manual labor. Raised beds 8 m² in size were prepared. Farmyard manure (25 t/ha) was added before plantation at each site and corms were treated with a fungicide solution to reduce fungal diseases. During September, corms were planted at different geographical locations; the plots were managed, and crop yield data were collected for 2 years. Saffron corms were planted at 12 cm depth with the spacing of 20 cm × 10 cm (row to row X plant to plant). Two factorial approaches were considered with 2 years and six planting sites in a factorial randomized block design (RBD). The experiment was executed in three replications, which were repeated for 2 years. Similar cultural practices designed by the Council of Scientific and Industrial Research-Institute of Himalayan Bioresource Technology, Palampur, Himachal Pradesh, India were implemented at different sites. Therefore, the differences between their performances are primarily due to different environmental conditions. The data on the number of flowers/m², fresh flower yield (kg/ha), fresh stigma yield (kg/ha) and dry stigma yield (kg/ha) were recorded. Growth parameters, viz., the number of leaves and leaf lengths, were recorded from January to February at all locations during its life cycle. During the flowering season, flowers were plucked daily at early hours in the morning before perianth openings to minimize the loss of volatile compounds of economic importance. After plucking, fresh stigmas were removed from the remaining flower immediately and were shade dried at room temperature.

Quality Analysis

Saffron dry stigmas were stored in the dark inside closed glass jars at room temperature (18–22°C) until further analysis (Gresta et al., 2008). Dried saffron stigmas (50 mg) were crushed in a mortar pestle to form a powder. Saffron powder was extracted using 70% methanol, followed by 5 min of vortex mixing and centrifugation at 8,000 rpm at 4°C for 10 min. The final volume of the extract was made up to 5 mL, with 70% methanol, and the samples were kept at 4°C for further analysis (Joshi et al., 2015).

Total Flavonoid Content

The flavonoid content of the saffron extract was estimated by using Joshi et al. (2015) method with minor alteration. The saffron methanol extract (1 mL) was diluted with distilled water, and then 5% of sodium nitrate was added to the solution. After 5 min., 3 mL of 10% aluminum chloride was added. After 6 min of incubation 2 mL of 1 M sodium hydroxide was added, and the sample was kept for 1-h incubation. Absorbance was taken at 510 nm wavelength against blank with a Shimadzu 8000 UV-Vis spectrophotometer. Total flavonoid content was expressed as mg/g equivalent of quercetin. Readings were taken in triplicate.

Total Phenolic Content

Total phenolic content was calculated using the method of Kumar et al. (2020). In 1.5 mL Folin-Ciocalteu Reagent (1 F.C.: 9 Water), 300 μL of methanol extract was added and incubated at room temperature for 10 min in dark conditions (Kumar et al., 2020). After that 1.2 mL of 2.5% Na₂CO₃ solution was added and the mixture was kept idle for 1 h. Absorbance was taken at 765 nm using a Shimadzu 8000 UV-Vis spectrophotometer. Gallic acid was expressed in the total phenolic content as mg/g equivalent. Readings were taken in triplicates.

Ultra Pressure Liquid Chromatography Analysis

Acquity UPLC-H Class System was used to analyze marker compounds, viz., crocin, picrocrocin and safranal, with an auto sampler, binary solvent manager and a column heater equipped with a PDA detector (Joshi et al., 2015). For separation, the BEH C18 column (2.1 mm × 50 mm, 1.7 μm) was used. The column-heater was set at a temperature of 35°C with flow rate of 0.2 mL per min). Eluents A and B were 0.1% formic acid in water and 0.1% formic acid in acetonitrile, respectively. The non-linear separation gradient was 0–0.5 min (90% A), 4 min (50% A), 9 min (20% A), 10 min (10% A), 11 min (10% A), 12 min (90% A), and 15 min (90% A). One μL of the sample was injected for analysis. Set the eλ PDA detector at 354, 254, and 340 nm. The mean area of three replicate injections was considered for quantification against calibration curves for each standard (Joshi et al., 2015).

Soil Analysis

Soil samples collected from each location were dried at room temperature and passed through a 30 mesh panel. The method of Mehlich No. 3 (Mehlich, 1984) was used to estimate available P₂O₅ and K₂O. The Walkley and Black method was used to determine the organic carbon, and the macro Kjeldahl method was used for the analysis of available nitrogen content (Black, 1965). The hydrometer method and pH meter (Black, 1965) were used to determine the soil texture and pH of the soil.

Statistical Analysis

Data collected on different yield and growth parameters of saffron was subjected to analysis of variance (ANOVA) in factorial randomized block design (RBD). Treatment variance was measured by the values of the least significant difference (LSD) at P = 0.05 and P = 0.01 by multiplying standard error of mean (SEM) values. For correlation studies, software XLSTAT 2017 was used to explore the relationship between growth and yield parameters. Past three software was used to analyze the principal component analysis (PCA). Performance measuring and ranking by main component analysis (PSR-PCA) provides a useful tool for evaluating performance under various factors (Coussement et al., 2016).

Results and Discussion

Climate and Soil Characteristics at Experimental Sites

Weather data of six different altitudinal locations is presented in Table 1. The highest altitude was observed for location L₁ (2,591 m) followed by L₂ (2,400 m), L₃ (2,195 m), L₄ (2,019 m), L₅ (1,565 m) and L₆ (1,472 m). This study identified and validated favorable environmental conditions under the western Himalayas for saffron cultivation through the MAXENT model. The altitude of selected locations varied from 1,472 to 2,591 m amsl. It is understood that saffron’s productivity and efficiency have a major impact on altitudinal variability (Mykhailenko et al., 2020). Throughout India, only the Kashmir valley represents one of the most significant saffron growing areas situated at an altitude of 1,585–1,677 m amsl under temperate climatic conditions (Ganaie and Singh, 2019).

TABLE 1

Table 1. Location, geographical and climatic characteristics of the study area of saffron crop.

The mean value of air temperature (maximum and minimum) was lowest in locations L₁ (20.9°C) and L₂ (23.9°C). In the case of maximum temperature, location L₃ showed the highest temperature (28.2°C) followed by location L₆ (27.8°C), L₄ (27.5°C), and L₅ (26.2°C). Location L₆ showed the highest average range of minimum temperature (17.7°C) compared with other locations. During the crop period, total rainfall received was 998.9 mm in L₁, 742.7 mm in L₃, 1,223.8 mm in L₂, 598.1 mm in L₆, 1,195.3 mm in L₄ and 654.8 mm in L₅, respectively (Table 1). Reports of previous studies explain the significant influence of climatic factors, viz., rainfall, and temperature on soil organic matter and other nutrients, ultimately affecting saffron cultivation (Kamyabi et al., 2014; Mykhailenko et al., 2020). Relative humidity was observed maximum in location L₂ (60.6%) compared with other locations; however, it was lowest in location L₃ (38.3%).

The soil pH was lowest in location L₁ (5.2) followed by L₅ (5.3), L₆ (5.6), L₂ (6.0), L₃ (6.2), and L₄ (6.4) (Table 2). Electrical conductivity was observed highest in L₁; however, lowest in L₂ and L₆ locations. According to Ganaie and Singh (2019), slightly alkaline soil with pH in the range of 6.3 to 8.3 and electrical conductivity in the range of 0.09 to 0.30 ds/m is most suitable for increasing saffron productivity. The soil textures of the selected study sites were sandy loam to clay loam. In Iran’s climatic conditions, sandy soil has been reported as a useful soil for raising the size of corms and flowers along with pH levels 6.8 and 7.8 (Mykhailenko et al., 2020). The percentage of soil organic carbon was significantly higher in location L₁ (1.2%) followed by L₂ (1.0%), L₃ (0.9%), and L₄ (0.8%). The soil of location L₅ was medium in organic carbon (0.7%) but low at location L₆ (0.4%). This increase in soil organic matter at high altitudes might be due to the higher input of organic matter and limited decomposition rate by lesser temperature and higher water retention capacity. The available nitrogen was medium in locations L₁ (269.4 kg/ha), L₃ (448.4 kg/ha), L₂ (313.8 kg/ha), and L₄ (283.5 kg/ha) yet low in locations L₆ (150.6 kg/ha) and L₅ (167.7 kg/ha). Across six different locations, the value of available phosphorus was medium in the range at all the locations except for location L₃. The content of available potassium was very high for location L₁ (887.5 kg/ha) followed by L₃ (856.3 kg/ha) and L₂ (640.8 kg/ha). Available potassium for location L₄ (439.3 kg/ha) was recorded higher in range, while, L₆ (229.4 kg/ha) and L₅ (246.5 kg/ha) locations were medium in potassium availability. This might be due to the variability of temperature and total rainfall in different geographical locations. The soils of most of the selected locations were sandy loam in texture except for locations L₆, which were sandy clay loam (Table 2).

TABLE 2

Table 2. Variation in soil physicochemical properties at different altitudinal locations.

Morphological and Productive Traits

The number of flowers/m² and fresh flower yield (kg/ha) was significantly higher during the 2019–2020 compared with the 2018–2019 season (Table 3). The increase in the number of flowers and fresh flower yield during the 2019–2020 was 9.9 and 14.7% compared with the 2018–2019 crop season. The rest of the parameters were not significantly affected by two different cropping years. The better quality corms were produced in the subsequent year, responsible for the higher number of flowers (Bayat et al., 2016; Mykhailenko et al., 2020). Significant differences in weather parameters in given years ensure a direct relationship between studied traits and the environment and weather conditions during both years.

TABLE 3

Table 3. Different altitudinal locations affect growth, yield and yield attributes of saffron.

The different altitudinal locations showed that the results of growth and yield traits were location-specific (Table 3). The number of flowers/m² was significantly higher in L₃ (32.8); however, it was lowest in L₆ (21.5). The interaction effect of years and locations highlighted the best performance of saffron when cultivated in location L₃ in terms of the number of flowers during the second year of crop season (Table 3). Fresh flower yield showed a significant high value in geographical location L₃ (117.5 kg/ha) followed by L₂ (111.6 kg/ha), L₄ (100.9 kg/ha), L₁ (85.3 kg/ha), and L₅ (82.0 kg/ha).

The number of flowers and fresh flower yield was significantly higher in a location situated at an altitude of 2,195 m (Suppa, Bharmour, H.P.) as compared with other places (Table 3). It might be because of the favorable environmental conditions, viz., average temperature (28.2°C), total rainfall (742.7 mm) and relative humidity (38.3%). It means that saffron required lower relative humidity to produce a higher yield of spice. Rainfall and temperature (23.0–27.0°C) are the crucial climatic factors controlling the growth and flowering in Crocus species (Aghhavani Shajari et al., 2015), and in this study, we have recorded temperature from 20.9 to 27.8°C, which is suitable for the cultivation of a saffron crop in selected locations.

Fresh stigma yield was significantly higher in location L₃ (12.5 kg/ha) and lowest in L₆ (6.4 kg/ha). The percentage increase in the L₃ location was 95.3% as compared with the L₆ location (Table 3). Dry stigma yield was significantly affected by different altitudinal variations. In location L₃, significantly higher dry stigma yield was observed compared with other different geographical locations. The percentage increase in L₃ was 50.0, 41.2, 33.3, 14.3, and 9.1% compared with L₆, L₅, L₁, L₄, and L₂ locations, respectively (Table 3). In our studies the average dry stigma yield varied from 1.6 to 2.4 kg/ha, however, in Greece the average yield ranges from 4.0–7.0 kg/ha, Italy from 3.4 to 10.0 kg/ha, in Morocco 2.0–2.5 kg/ha; in Spain 2.5–6.0 kg/ha and in Iran 3.0–5.0 kg/ha (McGimpsey et al., 1997; Gresta et al., 2008). The geographical location of Suppa, Bharmour, H.P. (L₃) significantly recorded higher dry stigma yield (2.4 kg/ha) as compared with other sites; however, the dry yield of stigma behaved statistically at par with Kapkote, Bageshwar, Uttarakhand (2.2 kg/ha) and Sathli, Bharmour, H.P. (2.1 kg/ha). The high input of organic matter and low decomposition processes is carried out at high altitudes, as temperature and temporal water saturation contributed to increased yield (Davidson and Janssens, 2006).

The lowest dry stigma yield was recorded in Palampur, Himachal Pradesh, as the altitude was relatively low and rainfall was higher as compared with other locations. It might also be related to low pH, organic carbon, available nitrogen, phosphorous, potassium, soil texture (Table 2) and other external environmental factors. Thus, in our studies saffron yield is close to the yield in Morocco and Spain, which showed high environmental adaptability to India’s latitude and produced higher cultivation rates. Specifically, the climatic conditions of these locations were desirable for the production of the saffron corm and other vegetative characteristics. Previous studies reported a 70% increase in saffron flowering when planted in soils with lighter consistency (Aghhavani Shajari et al., 2015); however, an 18% increase was reported in saffron stigma yield when grown in sandy soil as compared with the heavy soil (Khorramdel et al., 2015).

Growth parameters viz., the number of leaves/plant and leaf lengths were significantly affected by different geographical locations (Table 3). A significantly higher number of leaves/plant and leaf length were recorded in L₃ followed by L₅, L₂, L₄, L₁, and L₆. The interaction effect of year and geographical locations was significant on the number of leaves/plant and leaf length. A significantly higher number of leaves and leaf length were observed in location L₃ during the 2019–2020 crop season. Sand-particles usually cause higher pores, improving the soil’s permeability (Banihabib and Iranpour, 2012). This increases the root growth, improves the production and development of leaves (Fallahi et al., 2017). Panwar et al. (1995) in their study reported efficient cultivation of saffron at the highest altitude, which is most desirable for corm production, suggesting that elevation can play a crucial role. The development of corms is directly dependent on shoots or leaves at appropriate environmental conditions, including temperature, altitude and soil texture for saffron growth (Mykhailenko et al., 2020). Thus, to grow saffron as a valued commodity, one must take measures to achieve ideal conditions for the cultivation and processing of saffron corms in this area.

Regression and Correlation Analysis

Regression equations have been developed between independent variables, viz., altitude and total rainfall, and dependent variables, viz., fresh stigma yield and dry stigma yield (Figure 2). In the current study, fresh stigma yield and dry stigma yield increased with the corresponding increase in altitude (Figure 2A) and was found highest in location L₃ at altitude 2,195 m. Thus, a strong relationship was formed for different altitude with fresh stigma yield and dry stigma yield with an equation of y = −45.99 + 0.055x-1E-05x² (R² = 0.818; P ≤ 0.01) and y = −5.594 + 0.007x-2E-06x² (R² = 0.793; P ≤ 0.01), respectively (Figure 2). Total rainfall, fresh stigma yield and dry stigma yield enhanced with a decrease in altitude (Figure 2B). Therefore, a strong relationship was established by total rainfall with fresh stigma yield; y = 11.68 + 0.001x-1E-06x² (R² = 0.510; P ≤ 0.01) and dry stigma yield; y = 2.712-0.000x + 8E-06x² (R² = 0.566; P ≤ 0.01). The current study revealed that with the increase in altitude, fresh and dry stigma yield increases; however, with the increase in rainfall, yield attributes decrease (Figure 2).

FIGURE 2

Figure 2. Regression equation between independent variable, (A) altitude (m), (B) total rainfall and dependent variables, i.e., fresh stigma yield (kg/ha) and dry stigma yield (kg/ha). Altitude and rainfall are represented in the primary X-axis. Fresh stigma yield and dry stigma yield are presented in the primary Y-axis and secondary Y-axis, respectively.

Yield attributes, genotype and phenotypic correlation were measured to study the relation between saffron yields (Figure 3). Marker compound safranal showed a negative correlation with all the dependent and independent variables. Picrocrocin reported a higher positive correlation to altitude (r = 0.62) at a 5% significance level; however, rainfall (r = −0.58) and fresh flower yield (r = −0.24) showed a negative correlation. The compound crocin reported higher positive correlation with leaf length (r = 0.75), dry stigma yield (r = 0.75), fresh stigma yield (r = 0.89) and number of flowers (r = 0.81) at a 1% significance level and showed positive correlation with other variables. Saffron leaf length and leaf number have positive correlations for all characteristics studied except for rainfall. Dry stigma yield had the highest correlation with the number of flowers (r = 0.99) and fresh stigma yield (r = 0.93) at a 1% significance level. Dry stigma yield and fresh stigma yield have created a strong correlation with the increase in altitude and fresh flower yield, while they are negatively correlated with total rainfall. Fresh stigma yield recorded a substantially greater positive correlation with the number of flowers (r = 0.96) at a 1% significance level and altitude (r = 0.58) at a 5% significance level. Fresh flower yield reported positive correlation with number of flowers (r = 0.50) and altitude (r = 0.23), while, showed negative correlation with temperature (r = −0.16). A positive correlation was reported for the number of flowers with an increase in altitude. Total rainfall produced a negative correlation with yield attributes, yield and marker compounds of the saffron crop. The results indicated that selecting these traits is useful to improve the dry stigma yield of saffron (Figure 3). Genetic associations are more significant in correlation studies than phenotypic correlations due to removing environmental effects in calculating genetic correlation coefficients (Singh et al., 2005). Our study indicates that saffron yield attributes and the yield itself correlate greatly with altitude. Hence, optimum environmental conditions improve the yield attributes to produce a higher saffron yield (Figure 3).

FIGURE 3

Figure 3. Correlation matrix between dependent and independent variables. ^∗ and ^∗∗indicate that the corresponding values are significant at P = 0.05 and P = 0.01, respectively.

Quality Analysis

Principal Component Analysis

In this trial, we used PCA to analyze and acknowledge the relationship between independent variables (altitude and rainfall) with secondary metabolites (marker compounds, total phenolics and flavonoids) and how they vary according to the different altitudinal locations. The composition of secondary metabolites was taken as dependent variables. Three marker compounds, total phenolics and flavonoids of dry stigma samples from six different geographical areas were subjected to PCA for analyzing compositional variation. The chromatogram of marker compounds of different geographical locations has been presented in Figure 4. The correlation coefficient of secondary metabolites in PCA is determined by the cosine of the angle between their vectors (Dehghani et al., 2006). PCA plots jointly explained the variance of 74.5% (Figure 5). PC-1 explained 54.0% of the total variation, accounting for the positive contribution of altitude, crocin, picrocrocin, total phenolics and flavonoids; however, the negative contribution of total rainfall and safranal compound (Figure 5A). PC-2 explained 20.5% of the total variance, clearly distinguishing altitude, total precipitation, crocin and total phenolics in negative contribution with picrocrocin, safranal and total flavonoids in positive contribution (Figure 5B). Figures 5C,D showed the distribution of variables and treatments, respectively.

FIGURE 4

Figure 4. Chromatogram of marker compounds of saffron in different environmental locations.

FIGURE 5

Figure 5. The multivariate analysis of the mean value of independent variables, viz., altitude and total rainfall and secondary metabolites were conducted through principal component analysis. PC-1 and PC-2 jointly explained the variations of 74.5% (A–E). (A) represent variables of loading plot PC-1, (B) variables of loading plot PC-2, (C) variables presented as vectors in the space of the PCA, (D) loading scores of the treatments with PC-1 and PC-2, and (E) represent eigenvalue.

In the present study, three distinct clusters were observed in the score plot (Figures 5C,D). Cluster I from Langha, Palampur, Himachal Pradesh (L₅), explained a higher concentration of crocin (48.8 mg/g), picrocrocin (22.8 mg/g), safranal (3.8 mg/g), and total phenolics (5.9 mg/g); however, the total flavonoid count (4.9 mg/g) was higher in CSIR-IHBT, Palampur, Himachal Pradesh. In cluster II location of Sathli, Bharmour, Himachal Pradesh (L₄) explained the higher concentration of crocin (52.5 mg/g) and total flavonoid (mg/g); however, geographical location L₁ produced a higher concentration of picrocrocin (36.5 mg/g) and total phenolics (7.1 mg/g). Cluster III included two geographical locations (L₂ and L₃) showing a higher concentration of all the secondary metabolites in location L₃ (57.8 mg/g), except for crocin (57.8 mg/g) which was higher in the L₂ location (Table 4). This study also showed that the first three PCs with eigenvalues 3.8, 1.4, and 1.2 were most informative accounting for approximately 92.2% of the overall variance for all variables (Figure 5E). The saffron’s quality is a difficult parameter because it depends on several factors, in particular environmental conditions, including altitude, temperature, precipitations that affect the anatomy and quantity viz., crocin, picrocrocin and safranal (Cardone et al., 2019). Climatic factors, viz., altitude might also influence the quality of saffron spice, especially crocin (Lage and Cantrell, 2009). Fertile soils and a favorable climate distinguish these regions.

TABLE 4

Table 4. Cluster variability in secondary metabolites of saffron affected by altitudinal locations.

Factors affecting the accretion of compounds in plants used as medicinal and pharmaceutical raw materials and food additives are dependent on altitude, temperature, soil type, irrigation cycles, plant quality, and harvest times of saffron crop (Mykhailenko et al., 2020). The content of crocin was significantly higher in location Sathli, Bharmour, H.P. (L₄), and picrocrocin was markedly higher in Kapkote, Bageshwar, Uttarakhand (L₂); however, safranal was higher in the sample from CSIR-IHBT, Palampur, H.P. (Figures 4, 5). Increased crocin and picrocrocin at high altitudes might be due to total rainfall, air temperature, solar radiation, and soil characteristics, significantly affecting the accumulation of marker compounds in plants (Mykhailenko et al., 2020). The drying process is also a vital course ensuring a product’s quality (Maghsoodi et al., 2012). Previous studies indicated that with the increase in altitude, the crocin content increased. Thus, agronomical and climatic factors affect the quality of saffron (Lage and Cantrell, 2009). Total phenolics and flavonoids were significantly higher in high altitude locations as compared with lower altitude. The annual precipitation in comparison with high altitude regions was higher at low altitudes. Therefore, it was clear that different environmental conditions influenced the alteration of secondary metabolites of saffron.

Conclusion

This study revealed significant differences in different altitudinal locations for all studied traits for the saffron crop. Thus, the accurate selections of the geographical locations are considered important and essential factors. The yield attributes, viz., the flower number/m², fresh flower yield, and dry stigma yield were significantly higher in the second year at the geographical location of Suppa, Bharmour, H.P. (L₃), which confirmed that these traits are significantly influenced by environmental factors. Results of independent variables, genotypic and phenotypic coefficients of variation revealed a positive correlation for saffron yield with most characteristics. Thus, it was concluded that saffron could successfully be cultivated in the hilly regions of Himachal Pradesh. The studies suggest a need to have more research work on multiple sites in the same geographical locations with other factors that influence the quality of saffron production with altitude. Moreover, for each geographical location, there is a need to study the effect of saffron quality with meteorological data.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author Contributions

DK: experiment execution, data collection, soil and quality analysis, literature search, and manuscript writing. MT: data collection, statistical analysis and data presentation, literature search, and manuscript writing. RJ: quality analysis. AK: identification of sites through MAXENT model. RK: develop the idea, overall guidance, and manuscript editing. All authors contributed to the article and approved the submitted version.

Funding

The authors acknowledge financial support to the present study by the Council of Scientific and Industrial Research (CSIR), Government of India, under CSIR funding project MLP-0127. This is IHBT publication number 4699.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We would like to thank the Director, CSIR-IHBT, Palampur, (H.P.), India, for providing necessary facilities during study. We would also like to thank Kuldip Singh and Ramdeen Prasad for their technical support.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.657819/full#supplementary-material

References

Aghhavani Shajari, M., Rezvani Moghaddam, P., Koocheki, A., Fallahi, H. R., and Taherpour Kalantari, R. (2015). Evaluation of the effects of soil texture on yield and growth of saffron (Crocus sativus L.). Saffron Agron. Technol. 2, 311–322. doi: 10.22048/jsat.2015.8624

CrossRef Full Text | Google Scholar

Banihabib, M. E., and Iranpour, M. (2012). Effect of some soil physical and chemical properties on permeability in field conditions. Irrig. Water Eng. 2, 14–25.

Google Scholar

Bayat, M., Amirnia, R., Tajbakhsh, M., and Ramezani, M. (2016). Evaluation of saffron ecotypes for stigma yield and yield components using different maternal corm weights. J. Plant Physiol. Breed. 6, 53–64.

Google Scholar

Black, C. A. (1965). Method of soil analysis part 2. Chem. Microbiol. Prop. 9, 1387–1388.

Google Scholar

Cardone, L., Castronuovo, D., Perniola, M., Cicco, N., and Candido, V. (2019). Evaluation of corm origin and climatic conditions on saffron (Crocus sativus L.) yield and quality. J. Sci. Food Agric. 99, 5858–5869. doi: 10.1002/jsfa.9860

PubMed Abstract | CrossRef Full Text | Google Scholar

Cardone, L., Castronuovo, D., Perniola, M., Cicco, N., and Candido, V. (2020a). Saffron (Crocus sativus L.), the king of spices: an overview. Sci. Hortic. 272:109560. doi: 10.1016/j.scienta.2020.109560

CrossRef Full Text | Google Scholar

Cardone, L., Castronuovo, D., Perniola, M., Scrano, L., Cicco, N., and Candido, V. (2020b). The influence of soil physical and chemical properties on saffron (Crocus sativus L.) growth, yield and quality. Agronomy 10:1154. doi: 10.3390/agronomy10081154

CrossRef Full Text | Google Scholar

Coussement, A., Isaac, B. J., Gicquel, O., and Parente, A. (2016). Assessment of different chemistry reduction methods based on principal component analysis: comparison of the MG-PCA and score-PCA approaches. Combust. Flame 168, 83–97. doi: 10.1016/j.combustflame.2016.03.021

CrossRef Full Text | Google Scholar

Davidson, E. A., and Janssens, I. A. (2006). Temperature sensitivity of soil carbon decomposition and feedbacks to climate change. Nature 440, 165–173. doi: 10.1038/nature04514

PubMed Abstract | CrossRef Full Text | Google Scholar

Dehghani, K., Jenkins, S., Fischer, D. J., and Skalinski, M. (2006). “Application of integrated reservoir studies and probabilistic techniques to estimate oil volumes and recovery, Tengiz Field, Republic of Kazakhstan,” in Proceedings of the Society of Petroleum Engineers-1st International Oil Conference and Exhibition in Mexico Mexico, 81–98.

Google Scholar

Fallahi, H. R., Aghhavani-Shajari, M., Sahabi, H., and Feizi, H. (2017). Mother corm weight and soil amendment improves the vegetative and reproductive growth of saffron (Crocus sativus L.). J. Med. Spice Plants 22, 110–114.

Google Scholar

Fallahi, H. R., and Mahmoodi, S. (2018a). Influence of organic and chemical fertilisation on growth and flowering of saffron under two irrigation regimes. Saffron Agron. Technol. 6, 147–166. doi: 10.22048/jsat.2017.71511.1207

CrossRef Full Text | Google Scholar

Fallahi, H. R., and Mahmoodi, S. (2018b). Impact of water availability and fertilisation management on saffron (Crocus sativus L.) biomass allocation. J. Horti. Postharvest Res. 1, 131–146. doi: 10.22077/jhpr.2018.1487.1017

CrossRef Full Text | Google Scholar

Ganaie, D. B., and Singh, Y. (2019). Saffron in Jammu & Kashmir. Int. J. Res. Geogr. 5, 1–12. doi: 10.20431/2454-8685.0502001

CrossRef Full Text | Google Scholar

Garcia-Rodriguez, M. V., Lopez-Corcoles, H., Alonso, G. L., Pappas, C. S., Polissiou, M. G., and Tarantilis, P. A. (2017). Comparative evaluation of an ISO 3632 method and an HPLC-DAD method for safranal quantity determination in saffron. Food Chem. 221, 838–843. doi: 10.1016/j.foodchem.2016.11.089

PubMed Abstract | CrossRef Full Text | Google Scholar

Giorgi, A., Pentimalli, D., Giupponi, L., and Panseri, S. (2017). Quality traits of saffron (Crocus sativus L.) produced in the Italian Alps. Open Agric. 2, 52–57. doi: 10.1515/opag-2017-0005

CrossRef Full Text | Google Scholar

Gohari, A. R., Saeidnia, S., and Mahmoodabadi, M. K. (2013). An overview on saffron, phytochemicals, and medicinal properties. Pharmacogn. Rev. 7, 61–66. doi: 10.4103/0973-7847.112850

PubMed Abstract | CrossRef Full Text | Google Scholar

Gomez-Gomez, L., Trapero-Mozos, A., Gomez, M. D., Rubio-Moraga, A., and Ahrazem, O. (2012). Identification and possible role of a MYB transcription factor from saffron (Crocus sativus). J. Plant Physiol. 169, 509–515. doi: 10.1016/j.jplph.2011.11.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Gresta, F., Lombardo, G. M., and Avola, G. (2009). “Saffron stigmas production as affected by soil texture,” in Proceedings of the III International Symposium on Saffron: Forthcoming Challenges in Cultivation, Research and Economics, Vol. 850, (Leuven: Acta Horticulture), 149–152. doi: 10.17660/actahortic.2010.850.23

CrossRef Full Text | Google Scholar

Gresta, F., Lombardo, G. M., Siracusa, L., and Ruberto, G. (2008). Saffron, an alternative crop for sustainable agricultural systems. a review. Agron. Sustain. Dev. 28, 95–112. doi: 10.1051/agro:2007030

CrossRef Full Text | Google Scholar

Joshi, R., Rana, A., and Gulati, A. (2015). Studies on quality of orthodox teas made from anthocyanin-rich tea clones growing in Kangra valley. India. Food Chem. 176, 357–366. doi: 10.1016/j.foodchem.2014.12.067

PubMed Abstract | CrossRef Full Text | Google Scholar

Kamyabi, S., Habibi Nokhandan, M., and Rouhi, A. (2014). Effect of climatic factors affecting saffron using analytic hierarchy process(AHP); (Case Study Roshtkhar Region. Iran). Saffron Agron. Technol. 2, 75–90. doi: 10.22048/jsat.2014.6187

CrossRef Full Text | Google Scholar

Khorramdel, S., Eskandari Nasrabadi, S., and Mahmoodi, G. (2015). Evaluation of mother corm weights and foliar fertiliser levels on saffron (Crocus sativus L.) growth and yield components. J. Appl. Res. Med. Aromat. Plants 2, 9–14. doi: 10.1016/j.jarmap.2015.01.002

CrossRef Full Text | Google Scholar

Kumar, P., Partap, M., Rana, D., Kumar, P., and Warghat, A. R. (2020). Metabolite and expression profiling of steroidal alkaloids in wild tissues compared with bulb derived in vitro cultures of Fritillaria roylei – High value critically endangered Himalayan medicinal herb. Ind. Crops Prod. 145:111945. doi: 10.1016/j.indcrop.2019.111945

CrossRef Full Text | Google Scholar

Lage, M., and Cantrell, C. L. (2009). Quantification of saffron (Crocus sativus L.) metabolites crocins, picrocrocin and safranal for quality determination of the spice grown under different environmental Moroccan conditions. Sci. Hortic. 121, 366–373. doi: 10.1016/j.scienta.2009.02.017

CrossRef Full Text | Google Scholar

Maghsoodi, V., Kazemi, A., and Akhondi, E. (2012). Effect of different drying methods on saffron (Crocus sativus L.) quality. Iran J. Chem. Chem. Eng. 31, 85–89.

Google Scholar

McGimpsey, J. A., Douglas, M. H., and Wallace, A. R. (1997). Evaluation of saffron (Crocus sativus L.) production in New Zealand. N. Z. J. Crop Hortic. Sci. 25, 159–168. doi: 10.1080/01140671.1997.9514002

CrossRef Full Text | Google Scholar

Mehlich, A. (1984). Mehlich 3 soil test extractant: a modification of mehlich 2 extractant. Comm. Soil Sci. Plant Anal. 15, 1409–1416.

Google Scholar

Moradi Rikabad, M., Pourakbar, L., Siavash Moghaddam, S., and Popovic-Djordjevic, J. (2019). Agrobiological, chemical and antioxidant properties of saffron (Crocus sativus L.) exposed to TiO₂ nanoparticles and ultraviolet-B stress. Ind. Crops Prod. 137, 137–143. doi: 10.1016/j.indcrop.2019.05.017

CrossRef Full Text | Google Scholar

Mykhailenko, O., Desenko, V., Ivanauskas, L., and Georgiyants, V. (2020). Standard operating procedure of Ukrainian saffron cultivation according with good agricultural and collection practices to assure quality and traceability. Ind. Crops Prod. 151:112376. doi: 10.1016/j.indcrop.2020.112376

CrossRef Full Text | Google Scholar

Panwar, K. S., Saroch, K., and Vashist, G. D. (1995). Potential and prospects of saffron in temperate hills of Himachal Pradesh. Agric. Situat. India 50, 13–16.

Google Scholar

Rahimi, H., Shokrpour, M., Tabrizi Raeini, L., and Esfandiari, E. (2017). A study on the effects of environmental factors on vegetative characteristics and corm yield of saffron (Crocus sativus). Iran J. Hortic. Sci. 48, 45–52.

Google Scholar

Singh, S., Kumar, M., and Mahajan, R. K. (2005). The study of indoor radon in dwellings of Bathinda district, Punjab, India and its correlation with uranium and radon exhalation rate in soil. Radiat. Meas. 39, 535–542. doi: 10.1016/j.radmeas.2004.10.008

CrossRef Full Text | Google Scholar

Siracusa, L., Gresta, F., Avola, G., Lombardo, G. M., and Ruberto, G. (2010). Influence of corm provenance and environmental condition on yield and apocarotenoid profiles in saffron (Crocus sativus L.). J. Food Compos. Anal. 23, 394–400. doi: 10.1016/j.jfca.2010.02.007

CrossRef Full Text | Google Scholar

Zarghani, F., Karimi, A., Khorasani, R., and Lakzian, A. (2016). To evaluate the effect of soil physical and chemical characteristics on the growth characteristics of saffron (Crocus sativus L.) corms in Tornbat-e Heydariyeh area. J. Agroecol. 8, 120–133. doi: 10.22067/jag.v8i1.48511

CrossRef Full Text | Google Scholar