Garlic Ecotypes Utilise Different Morphological, Physiological and Biochemical Mechanisms to Cope with Drought Stress

Drought negatively affects plants by altering morphological, physiological and metabolic processes and ultimately reducing yields. Garlic (Allium sativum L.), an important member of the Alliaceae family, is also sensitive to drought and maximizing the yield of garlic bulbs is largely dependent on water availability. The objective of this study was to determine the effects of drought stress on morphological and physiological characteristics, as well as on phenolic, sugar, inulin and free amino acid content and antioxidant activity in two Croatian garlic ecotypes, ‘Istarski crveni’ (IC) and Istarski bijeli (IB). Drought was induced by using polyethylene glycol 8000 (PEG) solution (−0.6 MPa) starting 21 days after clove planting and lasted for 20 days. Drought reduced plant height, number of leaves and plant weight, but increased root length in both ecotypes compared to the control treatment. Among the physiological parameters, significant differences were observed between the two ecotypes studied in the spectral characteristics of the leaves, namely reflection in red, green and blue, VAL, values of the vegetation indices related to the chlorophyll content (CHI, GI), and the anthocyanin content (ARI). Ecotype IC showed higher antioxidant activity in the control treatment due to higher total phenolic content (TPC), but under drought conditions higher DPPH radical scavenging activity was determined in ecotype IB and higher values of FRAP in IC. Sucrose and glucose generally decreased under drought, while inulin increased in IB but decreased in IC. Total free amino acid content increased under drought in both ecotypes. In conclusion, drought tolerance of IB might be associated with increased accumulation of inulin and higher levels of amino acids, especially those shown to contribute to drought resistance. In IC, drought tolerance is associated with an increase in some amino acid compounds and better root growth in depth, probably due to a more efficient translocation of sucrose to the underground part of the plant.


Introduction
With the advent of global climate change, plants are increasingly exposed to various abiotic stresses, including extreme temperatures and prolonged drought [1]. Drought is one of the most important limiting factors for agricultural production, leading to a significant reduction in crop yields worldwide [2,3]. Drought affects morphological, physiological, These results are supported by the research of Chaudhry et al. [13], who reached a similar conclusion in onions (Allium cepa L.). In their study, all onion cultivars tested showed increased root length development under drought stress. They suggest that the increase in root length is due to the ability to mitigate osmotic stress by maintaining osmotic potential. Plant weight was higher in ecotype IC than in IB in both treatments ( Figure 1D). In our study, the interaction between ecotype and treatment was significant only for the number of leaves (Table 1). The reduction in the number of leaves during drought is the result of the plants primary drought defence mechanisms. By reducing leaf area and decreasing the number of leaves, the plant conserves water by reducing the area for transpiration [29]. Lack of water leads to a decrease in turgor, cell growth and division are inhibited, resulting in weaker leaf growth and faster senescence and shedding [30]. A decrease in plant height with increasing drought intensity in garlic was als mined by Marostica et al. [31]. The same result was obtained by Ajayi et al. [32] ranth and Hanci and Cebeci [33] in onion. In both studies, drought stress was f significantly affect plant height and other morphological characteristics.

Chlorophyll Fluorescence and Physiological Characteristics of Garlic
Significant differences in leaf spectral characteristics were found between studied garlic ecotypes (Table 1).
Namely, higher reflection in red, green and blue and higher VAL were obs IB compared to IC, while higher values of vegetation indices related to chlorophyll (CHI, GI) and anthocyanin content (ARI) and VAL were observed in the ecotype plementary Table S1). In addition, the ecotype IC exhibited higher minimum (F0) a imum fluorescence (Fm) compared to IB. However, this increase was proportio therefore did not affect the final Fv/Fm values, which were not significantly diffe tween these two ecotypes.
Sperdouli et al. [34], according to Moustakas et al. [35], found in Arabidopsis that the intensity of stress had a significant effect on the maximum efficiency of P tochemistry (Fv/Fm). In this study, the maximum efficiency of PSII photochemistry decreased 24 h after the beginning of drought stress, while it further decreased sixth day after the beginning of stress. However, ten days after the beginning of th the maximum efficiency of PSII photochemistry recovered. Further drought treatm severe for Arbidopsis and resulted in significantly decreased PSII [35]. PSII. Most did not find significant decreases in Fv/Fm under moderate drought stress [36][37][38], ing that ETR is unaltered under drought stress [39]. These discrepancies in Fv/ mainly due to the differences among plant species, the stage of plant growth an opment when drought occurs, the duration of the stressful period and their susce to drought stress [40]. Here we hypothesize that the more pronounced effect of treatment on examined multispectral and chlorophyll fluorescence characteris A decrease in plant height with increasing drought intensity in garlic was also determined by Marostica et al. [31]. The same result was obtained by Ajayi et al. [32] in amaranth and Hanci and Cebeci [33] in onion. In both studies, drought stress was found to significantly affect plant height and other morphological characteristics.

Chlorophyll Fluorescence and Physiological Characteristics of Garlic
Significant differences in leaf spectral characteristics were found between the two studied garlic ecotypes (Table 1).
Namely, higher reflection in red, green and blue and higher VAL were observed in IB compared to IC, while higher values of vegetation indices related to chlorophyll content (CHI, GI) and anthocyanin content (ARI) and VAL were observed in the ecotype IC (Supplementary Table S1). In addition, the ecotype IC exhibited higher minimum (F 0 ) and maximum fluorescence (F m ) compared to IB. However, this increase was proportional and therefore did not affect the final F v /F m values, which were not significantly different between these two ecotypes.
Sperdouli et al. [34], according to Moustakas et al. [35], found in Arabidopsis thaliana that the intensity of stress had a significant effect on the maximum efficiency of PSII photochemistry (F v /F m ). In this study, the maximum efficiency of PSII photochemistry (F v /F m ) decreased 24 h after the beginning of drought stress, while it further decreased on the sixth day after the beginning of stress. However, ten days after the beginning of the stress, the maximum efficiency of PSII photochemistry recovered. Further drought treatment was severe for Arbidopsis and resulted in significantly decreased PSII [35]. PSII. Most authors did not find significant decreases in F v /F m under moderate drought stress [36][37][38], suggesting that ETR is unaltered under drought stress [39]. These discrepancies in F v /F m were mainly due to the differences among plant species, the stage of plant growth and development when drought occurs, the duration of the stressful period and their susceptibility to drought stress [40]. Here we hypothesize that the more pronounced effect of drought treatment on examined multispectral and chlorophyll fluorescence characteristics are absent due to the morphological and biochemical changes which probably had a protective role regarding physiological functions. The response of plants to drought is complex and is reflected on the morphological, physiological and biochemical levels. Osmotic adjustment (OA), because of the accumulation of cellular compatible solutes, is an important adaptive mechanism to drought stress and is widely recognized to have a role in plant adaptation to dehydration mainly through turgor maintenance and the protection of specific cellular functions by defined solutes [41]. The well-known organic osmolytes include sugars, fructans (inulin) and free amino acids, especially proline. The substances mentioned are not only osmolytes but, together with phenolic compounds, act as ROS scavengers that reduce the oxidative stress that occurs as a result of primary abiotic stress. Higher levels of phenolic compounds are therefore associated with high antioxidant capacity in plants, including garlic [42][43][44]. To estimate the antioxidant capacity of plant extracts, several assays have been used, including DPPH (2,2-dyphenyl-1-picrylhydrazyl) and FRAP (ferric reducing antioxidant potential) [45]. Drought has different effects on the antioxidant activity of plant genotypes. One group of authors found an increase in antioxidant activity under the influence of drought [46,47], while another group found a decrease in antioxidant activity [48,49].
This study determined significant differences in total phenolic content (TPC) and antioxidant activity between the two studied garlic ecotypes. The interaction between ecotype and treatment was also significant ( Table 2). In the control treatment, the ecotype IC was found to have a higher total phenolic content as well as a higher antioxidant capacity as measured by DPPH radical scavenging activity and FRAP. Under drought conditions, higher DPPH was determined in IB and higher values of FRAP in IC ecotype (  Some authors have reported that accumulation of phenolic compounds in plants under abiotic stress leads to an increase in antioxidant activity [50,51]. In contrast, Astaneh et al. [52] showed that garlic loses phenolic compounds under drought induced by salinity, as evidenced by phenol leakage and lower total phenolic content of the stressed plants compared to the control. In a subsequent study, the same authors reported a decrease in antioxidant capacity, measured as DPPH, of garlic plants grown under drought induced by salinity [52]. Phenolic compounds play a key role in radical scavenging activity. However, since the content of phenolic compounds under drought is the same for both ecotypes in this study, the differences in antioxidant capacity measured by DPPH and FRAP under drought treatment are probably caused by differences in other antioxidant compounds such as sugars and/or inulin. Regardless, a very strong, significant correlation (average of both treatments) was determined between antioxidant capacity, measured by Some authors have reported that accumulation of phenolic compounds in plants under abiotic stress leads to an increase in antioxidant activity [50,51]. In contrast, Astaneh et al. [52] showed that garlic loses phenolic compounds under drought induced by salinity, as evidenced by phenol leakage and lower total phenolic content of the stressed plants compared to the control. In a subsequent study, the same authors reported a decrease in antioxidant capacity, measured as DPPH, of garlic plants grown under drought induced by salinity [52]. Phenolic compounds play a key role in radical scavenging activity. However, since the content of phenolic compounds under drought is the same for both ecotypes in this study, the differences in antioxidant capacity measured by DPPH and FRAP under drought treatment are probably caused by differences in other antioxidant compounds such as sugars and/or inulin. Regardless, a very strong, significant correlation (average of both treatments) was determined between antioxidant capacity, measured by DPPH and FRAP, and TPC in the garlic ecotypes. The positive correlation between TPC and antioxidant activity was observed in cereals [53], soybean [54] and lettuce seedlings [55], suggesting that the degradation of DPPH and FRAP largely depends on the presence of large amounts of phenolic compounds [56]. A significant difference in sugar content was determined between ecotypes and treatments. The interaction between ecotypes and treatment was also significant ( Table 2). Higher sucrose content was determined in the control treatment in ecotype IC ( Figure 3A), while ecotype IB had higher sucrose content in drought and higher glucose content in both treatments ( Figure 3A,B). Soluble sugars trigger the proliferation of organs and produce larger and thicker leaves [61]. In this study, glucose content showed a significant positive correlation with the number of leaves and plant height (Figure 4). This especially applies to ecotype IB in  The reason for the lower accumulation of sucrose in the drought treatment in the ecotype IC compared to IB could be a better translocation of sucrose into the root during drought and thus a stronger root growth, and IC indeed had much longer roots in drought than in the control ( Figure 1C). The deep root system is adapted for absorbing deep soil water under extreme drought conditions [57]. Since the major vegetative organs are underground, the plant root system relies on translocation of photoassimilates to sustain growth. Sucrose is the most common form of transported photoassimilate that is transported by phloem from the source tissue (e.g., shoots) to the sink tissue (e.g., roots) [58]. As Durand et al. [59] and Chen et al. [60] found in Arabidopsis, key players in phloem loading and sucrose transport from the shoot to the root are the sucrose transporters SUT and SWEET. The SWEET11 and 12 transporters in Arabidopsis are rapidly phosphorylated under drought and abscisic acid treatment. This phosphorylation enhances the oligomerization and sucrose transport activity of SWEETs, resulting in increased sucrose content in roots, improved root growth under drought and drought resistance [60]. Importantly, this enhanced root growth did not compromise shoot growth, as can be seen from plant height, number of leaves and plant weight of IC under drought stress ( Figure 1A,B,D). Drought reduces the intensity of photosynthesis and reduces CO 2 uptake, but does not significantly reduce respiration, so metabolism requires more soluble sugars.
Soluble sugars trigger the proliferation of organs and produce larger and thicker leaves [61]. In this study, glucose content showed a significant positive correlation with the number of leaves and plant height ( Figure 4). This especially applies to ecotype IB in the control treatment ( Figure 1B). At the same time, there was no correlation with plant weight. In this study, the drought treatment reduced the glucose content in the leaves in both ecotypes ( Figure 3B). This is in agreement with findings of Hlahla et al. [62] in vegetable-type soybean. The content of the monosaccharide glucose was reduced in favour of the disaccharide sucrose, which is confirmed by the significant negative correlation between glucose and sucrose ( Figure 4). Soluble sugars are required for osmotic regulation to maintain phloem turgor and to maintain phloem transport under drought stress [63]. However, a high concentration of soluble sugars reverses physiological processes in a concentration-dependent manner [64] and inhibits photosynthesis as a result of sugar accumulation [65].
Increased accumulation of inulin could protect cells from dehydration during drought [66]. The protective action of fructans can be attributed to their ability to insert themselves between the lipids of the membrane [66], and thus stabilize it during drought [67]. Vereyken et al. [68] demonstrated that inulin-type fructans protect membranes. Besides this, fructans have also been reported to have antioxidant properties [69]. Fructans, as well as other sugars, have been shown to be better scavengers of •OH radicals compared to O 2 − [70]. Since plants lack enzymatic •OH-scavenging mechanisms, high concentrations of non-enzymatic antioxidant mechanisms are used to neutralize ROS [71]. In this study, the ecotype IC had a significantly higher inulin content in the control, but in drought conditions IC reduced the inulin content. In contrast, IB showed a tendency to increase inulin content under drought conditions ( Figure 3C). The antioxidant activity of inulin in garlic is obviously very high, as indicated by the very high correlation coefficients of inulin with DPPH and FRAP (Figure 4). Based on the results obtained, the antioxidant properties of the ecotype IB in drought treatment are not only the result of TPC but also of inulin content.

Plant Material and Growth Conditions
In this study, garlic cloves (Allium sativum L.) of 'Istarski crveni' (IC) and Istarski bijeli (IB) ecotypes provided by the Institute for Agriculture and Tourism Poreč were used. The cultivation of plants was carried out in the growth chambers at the Faculty of Agriculture, University of Zagreb. Cloves were treated with POLYRAM ®® DF (Chromos-agro, Zagreb, Croatia) fungicide before planting. Individual cloves were planted in a substrate of a mixture of black and brown peat Kekkila TSM 3 (Kekkila Proffesional, Vantaa, Finland) and placed in the growth chamber. The temperature in the chamber was 22 °C, the relative humidity of 50%, and the photoperiod was 16/8 h (light/darkness) and light intensity 250 µmol m −2 s −1 provided by PhenoLight3 (PhenoVation, Wageningen, The Netherlands). A total of 50 IC cloves and 50 IB cloves were planted. After the initial growth, before inducing drought stress, plants of the same habitus (equal height, number of leaves, good general condition) were selected. A total of 54 plants were used in the experiment, of which 30 were IC and 24 plants were IB. In each treatment, 15 IC plants and 12 IB plants were tested.
The drought stress induction started 21 days after clove planting when plants had three developed leaves. The experiment was carried out in two treatments: (1) Control (C)plants were watered twice a week with tap water; (2) Drought stress (D)-plants were watered with a 15% solution of polyethylene glycol 8000 (PEG) (−0.6 MPa) twice a week. Drought stress was induced for 20 days.

Amino acid Compounds
Free amino acids act as osmoprotectants, and can function as ROS scavengers (e.g., Pro) [24], but also have an ameliorating role in drought-induced nitrogen uptake, through its reassimilation and maintenance of protein homeostasis [72]. In situations of insufficient carbohydrate supply due to a decrease in photosynthesis rates, which usually occur during stress conditions, plants can use amino acids as alternative substrates for mitochondrial respiration [24].
In this study, ecotype, treatment and ecotype x treatment interaction significantly affected amino acid content (except for proline, histidine, threonine and alanine) ( Table 2). In both ecotypes, IC and IB, an increase in almost all amino acids is observed in plants exposed to drought compared to the control, which is a strong signal that the plants are suffering from drought (Supplementary Table S2). Numerous studies have established a significant increase in amino acid content under drought [23,[73][74][75][76]. Namely, drought increases proteolysis in plants and consequently increases the content of free amino acids and their metabolites [75].
Proline (Pro) synthesis in plants is a common physiological response to drought. Pro plays a role as an osmolyte, antioxidant, chaperone, stabilizer of membranes and enzymes [76]. However, an extreme increase in proline is not good for plants because it comes at the cost of other amino acids [75]. Schafleitne et al. [77] found higher proline accumulation in drought-susceptible than in drought-resistant potato genotypes. In contrast, Ghodke [78] found that a higher increase in proline levels reflects the adaptation mechanism to drought present in the tolerant onion genotype compared to the droughtsensitive genotype. In both ecotypes studied here, Pro content in drought treatment was about twice as high as in plants exposed to control treatment (Table 3), indicating a similar stress response in these two ecotypes. The upregulation of genes related to tyrosine, glycine, serine and threonine metabolism in garlic under drought was found by Zhou et al. [25]. Our results are in partial agreement with those of Zhou et al. [25]. Significantly higher levels of serine and threonine were found under drought in both ecotypes, but the levels of tyrosine and glycine did not change significantly under drought in IC (Table 3). Apart from these two amino acids, drought had no effect on the content of methionine, phenylalanine, isoleucine and leucine in IC. In the ecotype IB, on the other hand, the content of isoleucine and leucine was significantly higher under drought compared to the control treatment. These two amino acids belong to the branched-chain amino acids that have been shown to be associated with drought-tolerant genotypes of potato [73] and tomato [72]. A significant positive correlation was found between total free amino acid content and root length (Figure 4).

Plant Material and Growth Conditions
In this study, garlic cloves (Allium sativum L.) of 'Istarski crveni' (IC) and Istarski bijeli (IB) ecotypes provided by the Institute for Agriculture and Tourism Poreč were used. The cultivation of plants was carried out in the growth chambers at the Faculty of Agriculture, University of Zagreb. Cloves were treated with POLYRAM ®® DF (Chromos-agro, Zagreb, Croatia) fungicide before planting. Individual cloves were planted in a substrate of a mixture of black and brown peat Kekkila TSM 3 (Kekkila Proffesional, Vantaa, Finland) and placed in the growth chamber. The temperature in the chamber was 22 • C, the relative humidity of 50%, and the photoperiod was 16/8 h (light/darkness) and light intensity 250 µmol m −2 s −1 provided by PhenoLight3 (PhenoVation, Wageningen, The Netherlands). A total of 50 IC cloves and 50 IB cloves were planted. After the initial growth, before inducing drought stress, plants of the same habitus (equal height, number of leaves, good general condition) were selected. A total of 54 plants were used in the experiment, of which 30 were IC and 24 plants were IB. In each treatment, 15 IC plants and 12 IB plants were tested.
The drought stress induction started 21 days after clove planting when plants had three developed leaves. The experiment was carried out in two treatments: (1) Control (C)-plants were watered twice a week with tap water; (2) Drought stress (D)-plants were watered with a 15% solution of polyethylene glycol 8000 (PEG) (−0.6 MPa) twice a week. Drought stress was induced for 20 days.

Morphological Analysis
Morphological measurements were performed 20 days after the induction of drought stress. The following morphological characteristics were analysed: plant height (the length of the bulb together with above ground part, cm), the number of leaves, root length (the length from the end of the bulb to the tips of the roots using a ruler, cm) and plant fresh weight (g).

Chlorophyll Fluorescence and Multispectral Analysis
Chlorophyll fluorescence and multispectral analysis were performed at the end of stress induction period (41 days after planting). All plants were recorded using a CropReporter TM (PhenoVation B.V., Wageningen, The Netherlands) device. All images are captured with the same lens (10 Mp lens, 200 Lp mm −1 resolution, 400-1000 nm spectral range) and CCD-camera (1.3 Mp, 1296 × 966 pixels), with real 14-bit signal resolution. The output is 16-bit RAW format, and automatic analysis of chlorophyll fluorescence, colour, and multispectral images was performed by DA TM software (PhenoVation B.V., Wageningen, The Netherlands). For chlorophyll fluorescence analysis, plants were imaged using an optimized slow fluorescence induction protocol [79], which includes dark adaptation, measuring the induction curve of dark-adapted plants followed by turning on the light for light adaptation and measuring the induction curve of light-adapted plants. For chlorophyll fluorescence measurements of dark-adapted plants (plants were dark adapted for 30 min), a saturating pulse of light of 4000 µmol m −2 s −1 for 800 ms was used. The minimum chlorophyll fluorescence (F 0 ) was measured after 10 µs, and the maximum chlorophyll fluorescence (F m ) was measured after saturation. After measuring the dark-adapted plants, the plants were left in the dark for 15 s, and then the actinic lights (300 µmol m −2 s −1 ) were turned on to allow the plants to light-adapt for 5 min. Steady-state fluorescence yield (F s ) was measured at the beginning of the saturation pulse, and maximum chlorophyll fluorescence (F m ) of light-adapted plants was measured at saturation, using saturation pulse intensity (4000 µmol m −2 s −1 ).
After the measurement, actinic light was turned off, and in the presence of far-red light, the minimal fluorescence yield of the illuminated plant (F 0 ) was estimated. These measured parameters were used for the calculation of different chlorophyll fluorescence parameters, which are shown in Table 4. Table 4. Chlorophyll fluorescence parameters, the equation for calculation, and the reference.
A list of all measured multispectral parameters and calculated vegetation indices are given in Table 5.

Biochemical Analysis
Plant leaves for analysis of biochemical parameters were collected 20 days after induced drought stress from all plants included in each treatment. The lyophilized plant tissue was milled to a particle size of 0.2 mm on a centrifugal mill (Retsch ZM-200, Haan, Germany). Samples (50 mg) were homogenized with 2.4 mm glass beads (Omni kit 19-670, Kennesaw, GA, USA) for 1 min at 5 ms −1 in 1.5 mL of aqueous methanol (80:20, methanol: water, v/v) using a bead mill (Omni Bead Ruptor Elite, Kennesaw, GA, USA). The samples were left to macerate for 1 h on a rotator (Biosan RS-60, Riga, Latvia) and subsequently centrifuged for 5 min at 16,000× g. The supernatants were filtered through a 0.22 µm nylon filter prior to analysis.

Total Phenolic Content
The total phenolic content (TPC) was determined according to Singleton and Rossi [86]. Briefly, 20 µL of the sample was mixed with 140 µL of 0.2 M Folin-Ciocalteu reagent, and, after 1 min, 140 µL of 6% sodium carbonate was added. The reaction mixture was incubated at 25 • C for 1 h and the absorbance was read at 750 nm (Tecan Infinite 200 Pro M Nano+, Männedorf, Switzerland). The TPC was standardised against the gallic acid and expressed as the mg of gallic acid equivalents per g sample in DW. The results were calculated against a standard curve of gallic acid (serial dilutions of gallic acid: 12.5, 25, 50, 75, 100, 150 and 250 mg L −1 ) and expressed as mg GAEQ g −1 DW. Determined R 2 was 0.9999.

Total Antioxidant Activity
The total antioxidant activity was evaluated using the FRAP assay (ferric-reducing antioxidant capacity) [87] and the DPPH radical scavenging activity (2,2-diphenyl-1picrylhydrazyl) assay [88]. Briefly, 100 µL of the sample was mixed with 200 µL of either freshly prepared FRAP reagent or 0.02 M DPPH radical for the FRAP or DPPH assays, respectively. The antioxidant capacity using the FRAP assay was evaluated by reading the absorbance at 593 nm (Tecan Infinite 200 Pro M Nano+, Männedorf, Switzerland) after 10 min of reaction time at 25 • C, while the DPPH radical scavenging capacity was evaluated by reading the absorbance at 517 nm (Tecan Infinite 200 Pro M Nano+, Männedorf, Switzerland) after 30 min of reaction time at 25 • C. Both the FRAP and DPPH values were calculated against a standard Trolox calibration curve (serial dilutions of Trolox-2, 5, 10, 25, 50, 75 and 100 µM) and expressed as µmol Trolox g −1 DW. Determined R 2 were 0.9998 and 0.9997, respectively.

Statistical Analysis
All experiments were set up in a completely randomized design. The normality of the residuals for each model were assessed by visual inspection of the diagnostic plots as suggested by Kozak and Piepho [91]. The data from all measurements were subjected to ANOVA. The Bonferroni post-hoc test at p ≤ 0.05 was used for means comparison. Correlation analysis was performed using RStudio (ver. 5.2; RStudio Desktop, Boston, MA, USA) at p < 0.05. The other statistical analysis of the data was carried out using the SAS/STAT ®® [92] program package.

Conclusions
Based on the results obtained, we determined that drought had a significant effect on the morphological characteristics of both ecotypes of garlic. While plant height, number of leaves and plant weight decreased, root length increased in drought compared to the control treatment. Although significantly higher plant weight was observed at IC, this is influenced by the genotype rather than the environment, since this ecotype also had higher weight in the control treatment. However, no significant influence of drought was observed in the physiological characteristic, except for CHI, F v /F m and F 0 . Apparently, the ecotypes were able to overcome the moderate drought stress without major damage to the photosystem due to their excellent antioxidant potential and the increase in the content of osmoprotectants. The garlic ecotypes responded differently to drought in terms of biochemical parameters. The drought tolerance of IB could be related to increased accumulation of inulin and higher levels of amino acids, especially those shown to contribute to drought resistance. In IC, drought tolerance is associated with an increase in some amino acid compounds and a decrease in the leaves' sucrose content due to translocation to the underground part of the plant to intensify root growth at depth.

Supplementary Materials:
The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/plants12091824/s1, Table S1. Average values of morphological and physiological characteristics, and chlorophyll fluorescence parameters from main factors, Table  S2. Average values of biochemical characteristics from main factors. Funding: This research was partially funded by the Operational Program Competitiveness and Cohesion 2014-2020 and the European regional fund under the specific Scheme "Scheme to strengthening applied research in proposing actions for climate change adaptation", grant number KK.05.1.1.02.0005.