Multivariate Discrimination of Some Grapevine Cultivars under Drought Stress in Iran

: Grapevine is one of the most important economic crops in horticulture, and drought stress is one of the most signiﬁcant threatening factors in the world. Therefore, the identiﬁcation and investigation of cultivars under drought stress are the basic steps and important goals in grapevine-breeding programs. In the present study, the 17 parameters of 14 grapevine cultivars under drought stress were ﬁrst scaled. Based on the initial information, we divided the 14 grape cultivars according to their resistance to drought stress into four groups: tolerant, semi-tolerant, semi-sensitive, and sensitive. Then, the utilization of multivariate techniques comprising principal component analysis (PCA), along with quadratic discriminant analysis (QDA), were utilized to choose the most substantial and accountable traits for the four groups’ discrimination. For the QDA, the 17 parameters were arranged into four sets. The discrimination for all parameters showed 96% correct classiﬁcation. The ﬁrst set includes shoot length (Shoot L), shoot number (Shoot N), leaf area (Leaf A), relative water content (RWC), and chlorophyll a (Chl a ) parameters that showed 71.5% correct classiﬁcation. The second set includes chlorophyll b (Chl b ), chlorophyll total, peroxidase (POX), and superoxide dismutase (Sod) parameters that had 75% correct classiﬁcation. Electrolyte leakage (EL), malondialdehyde (MDA), proline, catalase (CAT), and ascorbate peroxidase (APX) parameters were in the third set and had 87% correct discrimination. The best discrimination was obtained by the combination of the ﬁrst and third set, including the Shoot L, Shoot N, Leaf A, RWC, Chl a , EL, MDA, proline, CAT, and APX with 100% correct discrimination.


Introduction
Water shortage transpires when plant-water requirements cannot be effusively supplied due to insufficient rainfall and a reduced groundwater level, or water retention from soil particles. In this case, the level of transpired water exceeds the water absorbed by the roots [1]. Plants use different strategies due to water stress, and they create a series of morphoanatomical, physiological, and biochemical adjustments to reserve their hydric status [2]. Under stress conditions, the decrease in photosynthetic rate may be due to biochemical limitations triggered via water scarcity, such as reducing photosynthetic pigments, particularly chlorophylls, and a decline in the water prospects and comparative leaves' water content [3].
The major factor in plants that increases during drought stress is the production of reactive oxygen species (ROS) in the organs, such as chloroplasts, mitochondria, and ω = W w W s = W − W S W S ω = percentage of soil moisture, Ww = Initial weight, Ws = dry weight γ d = γ t 1 + ω γd = unit weight of dry soil volume, γt = unit weight of wet soil volume Then, the amount of water required for each treatment was calculated with the following formula: γ d × (FC 100, 75, 50 and 25% − ω) During the drought stress treatment (for two months), the pots were weighed every two days, and the water required for each pot was provided.

Growth Parameters
A leaf area meter (AM 300 Bio Scientific Ltd., Hoddesdon, UK) was used to measure the leaf area. A ruler was used to measure the length of the main branch, and the leaf number was counted.

Relative Water Content
Relative leaf-water content was measured according to the method of Yamasaki and Dillenburg [13]. First, the fresh leaf sample of 0.5 g (Wf ) was weighed and then placed in distilled H 2 O. The sample was weighed for a second time at the end of this period (Ws). They were placed in bags in the oven to dry completely (Wd) and at the end, the content of the comparative leaf water was formulated from the subsequent equation.

Photosynthetic Pigments
Photosynthetic pigments were measured according to Lichtenthaler's [14] method. In accordance with this method, the fresh leaf was ground with a mortar and pestle, and then added to 4 mL of acetone 80% and centrifuged. Chlorophyll a/b (µg/mL) = (Chlorophyll a/Chlorophyll b)

Electrolyte Leakage
To determine electrolyte leakage, fresh discs were sampled from each pot. The samples were washed three times with deionized water and placed in laboratory tubes, and deionized water was added to them. The tubes were placed at a constant temperature of 25 • C. After 12 h, the initial electrical conductivity of the environment (EC1) was measured using an EC meter (Milwaukee Mi 306 model). Then, the samples were boiled at 100 • C for 20 min to kill the tissues and release all electrolytes completely [15]. (1)

Lipid Peroxidation and Malondialdehyde Concentration
Heath and Packer's method was used to measure MDA. For this purpose, leaf tissue was ground with liquid nitrogen, then added to a phosphate buffer. One ml of 20% trichloroacetic acid (TCA) solution, containing 0.5% thiobarbituric acid (TBA), was added to 1 mL of the supernatant solution, and the mixture was exposed to a temperature of 95 • C at 30 min. Then, the mixture was immediately cooled. The supernatant absorbance was read at 532 and 600 nm [16].

Proline Content
Leaf tissue was ground with liquid nitrogen, and 10% sulfosalicylic acid was added to it. The mixture was filtered, then the filtered solution was added to the ninhydrin reagent and acetic acid and incubated for 1 h at 100 • C. Then, toluene was added to the mixture and vortexed. The light absorption of the supernatant red solution was read at 520 nm with a Genway 6705 spectrophotometer (Genway Biotech, San Diego, CA, USA) [17].

Soluble Sugar Total Content
The antron reagent method was used to measure the amount of soluble sugar. 0.2 g of fresh leaf tissue was ground and about 4 mL of 80% (v/v) ethanol was added to it and centrifuged. The supernatant was mixed with the antron reagent and placed at 100 • C for 10 min. The reaction was completed by rapid cooling in ice, and the adsorption was read at 620 nm by a spectrophotometer. Glucose was used as the standard, and a standard curve with 11 standard concentrations (0-100 g) of glucose was prepared [18].

Antioxidant Enzyme Activity Assay
About 0.5 g of fresh leaves were completely turned into powder with liquid nitrogen and were vortexed by adding 1.5 mL of extraction buffer, before being centrifuged. The supernatant was centrifuged again. The supernatant was transferred to the −80 • C freezer until use [19]. The extract was used to assay the peroxidase, superoxide dismutase, ascorbate peroxidase, catalase enzymes, and protein by the decomposition of hydrogen peroxide into water and oxygen by a spectrophotometer at a wavelength of 240 nm for one minute by Chance and Mahle's method [20], with slight changes. The enzymatic reaction rate was recorded as changes in OD/min [21]. The Giannopolitis and Ries technique measured the enzyme activity, which was based on the inhibition of superoxide dismutase enzyme reduction in nitro-blue tetrazolium (NBT). The absorbance was read at 560 nm [22]. Ascorbate enzyme activity was assayed by Nakano and Asada's method, with slight changes, by measuring the oxidation of ascorbate by a spectrophotometer at 290 nm for 3 min [23].
A total 2.5 mL of Bradford 1X solution mixed with 100 µL of leaf extract. Then, the absorbance was read at 595 nm with a spectrophotometer. In the leaf, tissues used the standard bovine serum albumin (BSA) protein to measure the concentration of soluble protein [24].

Statistical Analysis
For the statistical analysis, based on the initial information and according to Table 1, we placed 14 cultivars of grapevine genotypes into four groups and performed the discrim- ination based on these four groups. PCA was utilized to smear a discrimination process monitored by quadratic discriminant analysis (QDA) for distinguished genotypes. PCA was conducted to decrease the original data dimensions to a lesser figure of component groups by observing the liaison amid the parameters' measurements, using 95% confidence. QDA was conducted to create the discrimination model associated with the grapevine genotypes' biochemical and physiological reactions to salinity stress. The multivariate statistical analysis, comprising PCA and QDA, was conducted with JMP ® Pro software version 14.3.0 (SAS Institute, Cary, NC, USA).

Response of Grapevine Cultivars to Drought Stress
Physiological and biochemical factors from the study of 14 grape cultivars under four levels of drought stress are shown in Figure 1. That, the trend in shoot N and shoot L, leaf A, and RWC factors decreases by increasing drought stress. Additionally, the increasing trend in biochemical data is clearly shown in Figure 1. The total carbohydrate and proline content in 10 grape cultivars meaningfully increased during drought stress, and three cultivars stood at the top of the trend. Due to the importance of the evaluation of antioxidant enzymes in studies of drought stress, the increasing trend in catalase, peroxidase, ascorbate, superoxide, and peroxidase dismutase enzymes are well demonstrated in Figure 1.

Principal Components Analysis
PCA is a multivariate procedure that offers an extracting structure technique from a variance-covariance or association amongst dependent variables. It substitutes a new variable called a factor for the group of innovator qualities that are related to each other.
The determination of the number of principal components by data group can be ex-

Principal Components Analysis
PCA is a multivariate procedure that offers an extracting structure technique from a variance-covariance or association amongst dependent variables. It substitutes a new variable called a factor for the group of innovator qualities that are related to each other.
The determination of the number of principal components by data group can be expected from the principal components (PCs) analysis.
The first three PCs generated from this analysis that showed eigenvalues > 1 [25] accounted for 94.4% of the total variance in the data group. The first three PCs' eigenvalues were 6.91, 4.26, and 1.20 ( Figure 2a). Table 2 indicates that the factor loads represent the connections between the PC and the trait measurements. Loads with an absolute value greater than 0.435 have a sturdy influence ( Table 2). PCA could recognize the features that generate the primary dissimilarity amongst samples. The primary key module vindicated 40.7% and the subsequent 25.1% of the total variance share. The grapevine cultivar's projections under four levels of water stress with the two principal components (PC) were imaged ( Figure 2). An apparent separation of drought treatments was shown by PC1, so by increasing drought stress from normal conditions to severe drought stress, the score of the grapevine cultivar changed from negative to positive values. Water-availability reduction changes biochemical properties, such as EL, MDA, proline, CAT, and APX, and has a significant influence on the discrimination of grapevine cultivars under drought stress. EL, MDA, proline, CAT, and APX are variables in the first component with positive coefficients, and Shoot L, Shoot N, Leaf A, RWC, and Chl a, with positive coefficients had a significant outcome on group variables. In the second component, Chl b, Chl total, protein, POX, and Sod with positive coefficients had a significant effect on variable grouping. Sugar with positive coefficients had the greatest effect in the variables of the third component group.

Discriminant Analysis
The QDA converted our traits group into several canonical components, and our default tests for components make our cultivars into four groups: tolerant, semi-tolerant, semi-sensitive, and sensitive. The likelihood and entropy R square factor show every canonical power for the cultivar's discrimination. The possible factor (−2log likelihood) was utilized to designate each canonical statistical significance. The closer the likelihood was to zero and the entropy R square value was to 100 was indicative of a better discriminant. Figure 3 shows the result of quadratic discriminant analysis for all measured traits (all traits are in the training set). Another studied trait and the type of response to drought stress was prioritized. Based on drought response, grapevine cultivars were divided into four sets. Discrimination analysis was conducted based on each set.

Discriminant Analysis
The QDA converted our traits group into several canonical components, and our default tests for components make our cultivars into four groups: tolerant, semi-tolerant, semi-sensitive, and sensitive. The likelihood and entropy R square factor show every canonical power for the cultivar's discrimination. The possible factor (−2log likelihood) was utilized to designate each canonical statistical significance. The closer the likelihood was to zero and the entropy R square value was to 100 was indicative of a better discriminant. Figure 3 shows the result of quadratic discriminant analysis for all measured traits (all traits are in the training set). Another studied trait and the type of response to drought stress was prioritized. Based on drought response, grapevine cultivars were divided into four sets. Discrimination analysis was conducted based on each set.

First Set
The discrimination of the four groups of grapevine cultivars was conducted based on their Shoot L, Shoot N, Leaf A, RWC, and Chl a reactions to drought stress. Figure 4 shows the three-dimensional (3D) common diagram of the first three normal values acquired with QDA. The first canonical variance correspondence was 59%, the second canonical variance was 27%, and the third canonical variance was 13%: the initial three canonical variances explained 99% of the total variance.
The first and second canonical values showed a substantial influence on discrimination. The possible first, second, and third canonicals ratios were 0.63, 0.82, and 0.93, respectively, which indicated that the initial canonical value had a greater obligation to the cultivars' classification than the others. These reactions, as utilized PCs for discrimination, had a 71% accurate categorization through QDA using a −2Log possible factor equivalent to 64. The four groups discriminated by Shoot L, Shoot N, Leaf A, RWC, and Chl a had 17 misclassifications (16 misclassifications) ( Table 3).

Second Set
The discrimination of the four groups of grapevine cultivars was established on their Chl b, Chl total, POX, and Sod reactions to drought stress. Figure 5 demonstrates a threedimensional (3D) canonical diagram of the first three common values acquired with QDA.
The variance corresponding to canonical values 1, 2, and 3 was 85%, 13%, and 0.4%, respectively. The first three canonical values clarified 99% of the overall variance of the three canonical values; only the first canonical value had a substantial influence on discrimination. The possible three canonical ratios were equal to 0.48, 0.87, and 0.99, respectively, showing that the first canonical value had the maximum impact on grape cultivar discrimination (Table 3). These reactions, as utilized PCs for discrimination, had a 75%

Second Set
The discrimination of the four groups of grapevine cultivars was established on their Chl b, Chl total, POX, and Sod reactions to drought stress. Figure 5 demonstrates a threedimensional (3D) canonical diagram of the first three common values acquired with QDA.
The variance corresponding to canonical values 1, 2, and 3 was 85%, 13%, and 0.4%, respectively. The first three canonical values clarified 99% of the overall variance of the three canonical values; only the first canonical value had a substantial influence on discrimination. The possible three canonical ratios were equal to 0.48, 0.87, and 0.99, respectively, showing that the first canonical value had the maximum impact on grape cultivar discrimination (Table 3). These reactions, as utilized PCs for discrimination, had a 75% accurate categorization through QDA using a −2Log probable feature equivalent to 106.4. The four groups discriminated by Chl b, Chl total, POX, and Sod had 25% or 14 misclassifications (Table 3).

Third Set
In this set, EL, MDA, proline, CAT, and APX responses to drought stress were used to discriminate the four groups of grape cultivars. Figure 6 displays a three-dimensional (3D) canonical diagram of the first three canonical values acquired with QDA. The first three cones accounted for 100% of the total variance, and the variances corresponding to the first, second, and third canonical values were 79%, 14%, and 5.8%, respectively. The first canonical value, with a likelihood ratio of 0.4, had a substantial effect on classification and had the highest responsibility; two other canonical values had an insignificant effect on the classification of grape cultivars. The accurate categorization was 87%, and the −2Log probable feature was 54.67. The QDA by EL, MDA, proline, CAT, and APX responses for the discrimination of four groups of grape cultivars had 18 misclassifications (23.21%) ( Table 3).

Third Set
In this set, EL, MDA, proline, CAT, and APX responses to drought stress were used to discriminate the four groups of grape cultivars. Figure 6 displays a three-dimensional (3D) canonical diagram of the first three canonical values acquired with QDA. The first three cones accounted for 100% of the total variance, and the variances corresponding to the first, second, and third canonical values were 79%, 14%, and 5.8%, respectively. The first canonical value, with a likelihood ratio of 0.4, had a substantial effect on classification and had the highest responsibility; two other canonical values had an insignificant effect on the classification of grape cultivars. The accurate categorization was 87%, and the −2Log probable feature was 54.67. The QDA by EL, MDA, proline, CAT, and APX responses for the discrimination of four groups of grape cultivars had 18 misclassifications (23.21%) ( Table 3).

Fourth Set
The fourth set was obtained from a combination of traits in sets one and three, including, the Shoot L, Shoot N, Leaf A, RWC, Chl a, EL, MDA, proline, CAT, and APX responses of the four groups of grape cultivars under different water-supply regimes, and utilized for discrimination by QDA. Figure 7 indicates the 3D canonical schematic of the first three canonicals.
The three canonical values clarified a total variance of 100%: the variance conforming to canonical values one, two, and three were 64%, 21%, and 13%, respectively. The first canonical value had a substantial influence on the cultivars' classification. Conversely, the first, second, and third canonical value had a probable ratio equivalent to 0.24, 0.55, and 0.79, respectively. The most significant impact on discrimination was made by the first canonical value. These retorts, as utilized PCs for discrimination, had a 100% accurate categorization through QDA using a −2Log probable factor identical to 2.73. The Shoot L, Shoot N, Leaf A, RWC, Chl a, EL, MDA, proline, CAT, and APX responses used for discrimination had no misclassifications (Table 3).

Fourth Set
The fourth set was obtained from a combination of traits in sets one and three, including, the Shoot L, Shoot N, Leaf A, RWC, Chl a, EL, MDA, proline, CAT, and APX responses of the four groups of grape cultivars under different water-supply regimes, and utilized for discrimination by QDA. Figure 7 indicates the 3D canonical schematic of the first three canonicals.

Discussion
The preliminary results of this research showed that the drought stress decreased the relative water content and photosynthetic pigments and increased the electrolyte leakage of the cell membrane due to the damage to it, causing more production of proline, malondialdehyde, protein, and antioxidant enzymes. Therefore, under drought stress, The three canonical values clarified a total variance of 100%: the variance conforming to canonical values one, two, and three were 64%, 21%, and 13%, respectively. The first canonical value had a substantial influence on the cultivars' classification. Conversely, the first, second, and third canonical value had a probable ratio equivalent to 0.24, 0.55, and 0.79, respectively. The most significant impact on discrimination was made by the first canonical value. These retorts, as utilized PCs for discrimination, had a 100% accurate categorization through QDA using a −2Log probable factor identical to 2.73. The Shoot L, Shoot N, Leaf A, RWC, Chl a, EL, MDA, proline, CAT, and APX responses used for discrimination had no misclassifications (Table 3).

Discussion
The preliminary results of this research showed that the drought stress decreased the relative water content and photosynthetic pigments and increased the electrolyte leakage of the cell membrane due to the damage to it, causing more production of proline, malondialdehyde, protein, and antioxidant enzymes. Therefore, under drought stress, cultivar increases the accumulation of adaptive osmolytes, such as proline, antioxidant enzymes and proteins, regulates intracellular osmosis, maintains cellular turgor, and reduces damage to cell membranes, and prevents enzyme degradation (Figure 1). According to the results of Garashiligh and Tukilgan cultivars, the most drought-tolerant cultivars, as well as Aghshiligh, Kechi Amjaie and Copake Bogan cultivars, are drought-sensitive cultivars, and it is recommended that further molecular and biochemical analyses be performed on these cultivars. For molecular and breeding research on grapes in arid regions, it is recommended that these resistant cultivars are considered and used in future research.
Examining four sets of traits for the discrimination of grapevine cultivars showed that in the first set (Shoot L, Shoot N, Leaf A, RWfC, and Chl a), the second set (Chl b, Chl total, POX, and Sod), and the third set (EL, MDA, proline, CAT, APX) there are an element of grapevine growth, and that physiological and biochemical parameters cannot discriminate well. The fourth set, in the combined first and second sets, showed a clear discrimination of grapevine cultivar with 100 correct classifications. In the studies of Rezazad Bari in 21 grapevine rootstocks under salinity stress, the best discrimination was shown by the fourth set with 100% correct discrimination, including chlorophyll a, and chlorophyll b, total protein content, total carbohydrate content, proline, and glycine-betaine characteristics [26].
The results showed that with the increase in the level of drought stress in 14 studied grape cultivars, plant growth parameters decreased significantly. By reducing the leaf area, the sunlight absorption level and the plant photosynthetic level is reduced, finally leading to a reduction in dry matter and plant yield. Leaf relative water content (RWC) is a vital water status indicator in plants: drought stress troubles the equilibrium amid the H 2 O supply to the leaf tissue and transpiration rate, and the leaves' comparative water content is reduced [27]. It has been found that the relative water content of leaves shows a good correlation with drought resistance and osmotic stresses in plants [28]. Tsegay, in 2014, identified that under stress conditions in the grapevine, drought-tolerant plants have a high-water potential by various mechanisms, such as an increased water uptake and a reduced water loss and maintenance of turgor, primarily by osmotic adjustment, activation of respiration, and accumulation of osmolytes and proteins [29].
There was a positive correlation between Shoot L, Shoot N, Leaf A, RWC, and Chl a in the first set that displayed a decrease in the volume of RWC under the impact of stress conditions, which can decrease Shoot L, Shoot N, Leaf A, and Chl a.
Drought stress and water deficit decreased photosynthetic activity, leaf chlorophyll synthesis, and chlorophyll a/b proportion in soybeans [30]. Recent study outcomes showed a reduction in chlorophyll a/b content. Drought stress also caused a decrease in the photochemical efficiency of photosystem II, a change in stomatal movement, a decrease in the abundance of Calvin cycle proteins, and a disruption in the plants' water status, which led to a weakening in plant production [31]. In addition, drought stress causes damage to enzymes associated with photosynthesis and result in a loss of photosynthetic pigment content [32]. The decrease in photosynthetic pigment content due to drought stress in grapes cv. has been investigated by Romeiko [33]. The results showed that increasing drought stress to 25%, the amount of chlorophyll a, chlorophyll b, and chlorophyll a/b in the leaves of all studied cultivars was significantly reduced.
One of the most critical parameters for identifying plants with high resistance to drought stress is electrolyte leakage [34]. By increasing drought stress, ion leakage in all grape cultivars is increased significantly. Drought-stress disruption of the cell membranes and inactivation of ion pumps located in cell membranes leads to an increase in the rate of electrolyte leakage [35]. Similar findings have been obtained from the role of stress in electrolyte leakage on wine grape cultivars, [35], fig [36], kiwifruit [37], spinach, and lettuce in vitro [34]. In addition, cell membrane integrity during drought stress is a sign of the control mechanism in the tolerance of desiccation, so cultivars such as Garashiligh, and Tukilgan, which have higher membrane stability, are more resistant to dehydration stress.
Excessive ROS stages are intensified due to severe drought stress on the plant's capacity to scavenge, and due to lipid peroxidation (LP) in biological membranes [38]. In the present study, drought stress significantly affected the physiological processes of the cell and increased MDA content due to increased lipid peroxidation and cell membrane damage. This finding is consistent with the studies for oats (Avena spp. L) [39], and apple trees (Malus spp. Mill.) [40] under drought.
During drought stress, proline is activated as a cellular osmotic regulator between the cytoplasm and the vacuole, and by detoxifying ROS, it protects the membrane and establishes membrane antioxidant enzymes [41]. Proline content increases significantly during drought stress in leaves of different grape rootstocks, such as Rasha [26], 110 Richter [41], and Narinc grape leaves [35]. Ous's findings are consistent with the findings of Kamangar and Haddad, reporting that proline content in the leaves of the Vitis vinifera grapevine increased in response to drought stress [42].
With the increase in drought stress in the studied cultivars, the amount of total soluble sugar increased. In most of the studied cultivars, a significant increase was observed in total soluble sugar of up to 50%. No significant change was observed at 25% field capacity, which may be due to high stress and reduced photosynthetic activity in plant products.
Statistical analysis showed that, with the decrease in irrigation level, antioxidant enzymes, such as CAT and SOD, and glutathione-ascorbate cycle enzymes, such as POD and APX, increased to protect against the effects of ROS (such as superoxide, hydrogen peroxide, and singlet oxygen). An excessive production of reactive oxygen species in drought stress leads to the peroxidation of cell membrane lipids and the degradation of enzyme proteins and nucleic acids [5]. The primary mechanisms of drought resistance are the ROS elimination through enzymatic as well as non-enzymatic methods, comprising stress protein expression, and cell membrane stability [4]. The CAT enzyme in peroxisomes and APX enzymes scavenged the resulting hydrogen peroxide [33]. In this study, a significant increase in catalase activity was observed in the stressed leaves of the control plant. The high increase in ascorbate peroxidase indicated that this enzyme is responsible for protecting chloroplasts, enhancing electron flow under stress, and is the main product of AOS activity [43]. Wang reported that APX is an essential enzyme for H 2 O 2 elimination. Additionally, the results showed that, with the increase in stress levels, the amount of peroxidase and ascorbate peroxidase enzymes increases significantly [44]. Increasing drought stress levels, the activity of antioxidant enzymes called superoxide dismutase, ascorbate peroxidase, catalase, and guaiacol peroxidase increases in cv. Romeiko [33] and ascorbate peroxidase, catalase, and superoxide dismutase increases in cv. Tempranillo [45]. In this study, antioxidant enzymes (CAT, SOD, POD, and APX) content increased in 14 grape cultivars, increasing the level of drought stress to 25% FC; these results are consistent with, Kamangar and Haddad [43] and Mirzaee [46].
Drought stress can also modify the volume and leaf protein conformation, resulting in alterations in soluble proteins and structural proteins [47]. In our study, with increasing drought stress, more proteins were observed in resistant cultivars. Several studies have shown that, under drought stress, some proteins, such as aquaporins, dehydrins, MYB proteins, and mitogen-activated protein kinases, are transcribed in grapevines, and play a significant role in plant protection under stress, which are copied differently depending on the cultivar [48].

Conclusions
In summary, we found statistically significant associations between compatible osmolytes, antioxidant enzymes, proteins, and drought tolerance after analyzing 14 grapevine genotypes. Catalase and superoxide dismutase and glutathione-ascorbate cycle enzymes, such as peroxidase and ascorbate peroxidase, increased to protect against the effects of reactive oxygen species. We conclude that physiological traits alone could not do the separation of the grapevine cultivars well, and biochemical traits should be studied along with physiological traits. Using multivariate statistical analysis and incorporating physiological and biochemical properties of grapevine cultivars under drought stress, our results imply that discrimination of four groups of grapevine cultivars leads to the conclusion that Garashiligh and Tukilgan are the most drought-tolerant cultivars, whereas Aghshiligh, Kechi Amjaie, and Copake Bogan are drought-sensitive ones. The information obtained in our study can be helpful to farmers in selecting grapevine cultivars suitable for local environmental conditions, as well as for grapevine breeders.

Data Availability Statement:
The data that support the findings of this study are available from the corresponding author upon reasonable request.