Water Stress Alters Physiological, Spectral and Agronomic Indexes of Wheat Genotypes

1. Introduction

Wheat (Triticum aestivum L.) is the world's second most widely cultivated crop, with an estimated cultivated area of about 200 million hectares [1]. The consumption of wheat grains accounts for 19% of calories in the global human diet since they are rich in carbohydrates and have a higher protein content than other cereals such as rice, corn and rye [1, 2]. Brazil has an annual demand of about 12 million tons of wheat grains but produces less than half of its needs and imports the rest [3]. Therefore, looking for alternatives to increase crop productivity and include new areas to reduce foreign dependence and expenditure on wheat imports is important.

The Cerrado, a vast savanna biome in the Brazilian highlands, has been an alternative for rainfed or irrigated wheat cultivation [4], as it has favorable climatic and soil conditions. During the winter in the Cerrado region, plants have high yields under irrigation conditions and should have efficiency in water use. Differently, cultivation during the summer at the end of the rainy season results in wheat sowing in the off-season, and the main limitation is the dry periods, which require drought-tolerant plants [4]. The cultivation of summer wheat (early crop) depends, among other factors, on the previous crop cycle, water availability, and drought-tolerant cultivars, and diseases such as brusone (Pyricularia oryzae) are the main limiting factors for wheat cultivation in the off-season [5].

However, in rainfed (off-season) cultivation, irregular rainfall and prolonged dry periods, known as 'veranicos' [6], are common in this region. Among the environmental changes that affect crop development, the ones caused by the water deficit, which intensified in recent years due to global warming [7, 8] stand out. Drought is a phenomenon that occurs in large extensions of cultivated areas, causing considerable losses in crop yield [9].

The reduction or irregularities in precipitation in the critical stages of crop development have been reported, especially in the Midwest region of Brazil, one of the largest grain producers in the country [10], with drastic reductions in crop productivity [11]. In this scenario, understanding the mechanisms used by plants to tolerate water stress is crucial to identifying soybean and wheat cultivars with tolerance and/or efficiency in water use and, consequently, decreasing the productivity and grain quality losses of these crops.

Drought can be minimized through proper management [12] and genetic tolerance. To unravel genetic tolerance, it is necessary to identify genotypes that can maintain or minimally reduce their yield potential and grain quality under conditions of low water availability [3]. Therefore, increasing drought tolerance in wheat genotypes is a sustainable and economical solution to increase production in environments and seasons with water limitations [13] without using natural water resources.

The survival of plants under drought depends on several mechanisms acting simultaneously: reduction of the growing cycle, stomatal opening and transpiration rate, and the development of antioxidants that maintain osmotic adjustments at the tissue level [14]. Consequently, morphophysiological and crop yield changes may occur [15].

Because of the water deficit effects on plants, genotypes tolerant to water stress can be identified by combining several variables, namely drought tolerance indices [4] (Soares et al., 2021), water use efficiency [16], productivity components [17], physiological index [14, 18], and vegetation (spectral) indices [19, 20], and these variables may be correlated with grain yield and quality [21].

New technologies and methodologies that measure these variables should offer advantages in relation to yield levels, applicability in field conditions (in various climatic conditions), and speed in obtaining this information with existing and consolidated. The use of tools that provide a rapid assessment and are non-invasive and destructive at any stage of crop development and with high precision on the characteristics of plants in relation to abiotic and biotic factors has been used in the studies of the interaction plant environment [22].

Among these tools, the use of sensors in agriculture can provide valuable information on plant physiology under water stress by detecting changes related to plant structure, pigments and photosynthetic efficiency [23].

Using sensors coupled on land platforms [24] or unmanned aerial vehicles [25] has the advantage of a rapid, non-destructive evaluation with high precision. However, there is a need to validate the efficiency of these tools. Among these sensors, there are multispectral cameras coupled in unmanned vehicles, which can be used to select cultivars adapted to different conditions [26], an indispensable tool in breeding programs.

Thus, multispectral sensors that capture wavelengths in the visible and near-infrared range of the electromagnetic spectrum allow the calculation of various vegetation indices [22] and can be used for rapid and non-destructive selection of drought-tolerant wheat genotypes [27] with high precision [28]. Multispectral traits derived from NIR, red and green bands showed a strong relationship with wheat biomass, water use efficiency, photosynthesis and grain yield. The study of morphophysiological, spectral and productivity changes due to water stress in wheat genotypes is important to select cultivars more tolerants to water stress [29].

Therefore, the characterization of genotypes tolerant to water stress in rainfed cultivation and greater water use efficiency in irrigated cultivation under Cerrado [4] conditions is important to increase grain yield and quality. This work hypothesizes that non-destructive physiological responses of wheat, as well as hightroughput vegetative indices, are related to soil water availability and can be used to discriminate tolerant genotypes to water stress. Thus, this work aimed to validate the use of sensors as a tool for selecting wheat genotypes for drought tolerance through morphophysiological and agronomic assessments under field conditions.

2. Results

2.1. Variable contributions in the multivariate response

Data were submitted to joint multivariate analysis of variance, and significant differences were found for the sources of variation genotypes (p < 0.01), water regime (p < 0.01), year of cultivation (p < 0.01), and the interaction genotypes x water regime (p < 0.01) (Table 1). However, there was no significant (p > 0.05) effect of the interaction genotype x water regime x cropping year (p = 1). The significance of genotype x water regime interaction shows that genotypes respond differently to water availability.

After a joint multivariate statistical analysis, the wheat genotypes were grouped in Group 1: Aliança; BRS 254; BRS 404; CPAC 01019; PF 020062; PF 080492; group 2: BR 18_Terena; MGS Brilhante; PF 020037; PF 120337; group 3: BRS 264; BRS 394; CPAC 01047; CPAC 07258; CPAC 8318; CPAC 9110; PF 100368; TBIO. Average values of the characteristics analyzed in the wheat crop according to the groups of genotypes and water regimes in 2018 and 2019 are presented in Table 2 and Table 3. In general, changes in water availability affected the morphophysiological, spectral characteristics and crop yield, with different levels in the genotypes studied.

In Figure 1, three groups of genotypes were identified. Group 1 with six genotypes (33.33%) (Aliança; BRS 254; BRS 404; CPAC 01019; PF 020062; PF 080492); group 2 has four genotypes (22.22%) (BR 18_Terena; MGS Brilhante; PF 020037; PF 120337) and group 3 has eight genotypes (44.45%) (BRS 264; BRS 394; CPAC 01047; CPAC 07258; CPAC 8318; CPAC 9110; PF 100368; TBIO Sintonia).

The genotypes in group 1 had a higher drought resistance index (DRI= 1.01) than the genotypes in group 3 (DRI = 0.93), with intermediate values in group 2 (DRI = 0.95) (Table 4). Vegetative indices accounted for 97% of the total divergence observed among genotypes under water stress (Table 3). Indices based on the near-infrared band (OSAVI (33%), SAVI (33%), NDRE (11%), DVI (11%), GNDVI (3%) and NDVI (1%)) were most strongly associated with the differentiation of genotype groups. In addition, they were among the variables most affected by water stress.

The singular value decomposition of means of combinations of water regime levels and genotype groups retained 86% of the total variability in the first principal coordinate (Figure 2). However, there is a strong contrast in responses among genotype groups that received WR1 (22% of CET replacement) compared to the same groups with WR4 (100% of CET replacement). This means that the selection process must be carried out exactly under the conditions at the field, either under irrigation or under rainfed conditions.

The variables gas exchange (A, g_s, E) and vegetative indices (NDVI, GNDVI, GRVI, DVI, NDRE, SAVI, PRI, OSAVI, TCARI) are correlated (Figure 2 and Figure 3), and the higher correlations were between A and NDGI 90.78), gs and RVI (0.73) and E and NDGI (0.70).

The variables iWUE, WUE, and Rm/Sm ratio positively correlate. In contrast, WUE and net CO2 assimilation (A) negatively correlated (Figure 2 and Figure 3). There are correlations between vegetative and gas exchange indices as they formed sharp angles between these variables, showing a correlation with grain yield, which should be the main criteria for the selection because yield is the most important variable.

All IVs except TCARI/OSAVI (TO) have a strong positive correlation (>0.6, p<0.05) with grain yield, photosynthetic rate, stomatal conductance, and transpiration (Figure 3). Correlations between 0.2 and 0.6 with chlorophyll a and b, shoot mass and thousand-grain mass and correlations between 0.05 and 0.19 with root mass. Negative correlations above 0.5 occur between these variables with WUE, iWUE, TO and ratio chlorophyll a to chlorophyll b. Positive correlations are observed between the WUE and iWUE variables (0.3, p<0.05) (Figure 3).

Photosynthesis (A) and stomatal conductance (g_s) showed a high correlation with productivity (0.62 and 0.63, respectively), but these evaluations can not be done on a large scale. On the other hand, the vegetative indices TCARI, NDVI, OSAVI, SAVI, PRI, NDRE and GNDVI showed the highest correlation with productivity (>0.7, p<0.05). They can be useful in breeding programs (Figure 3) as they can be evaluated on a large scale. The vegetation indices were generally correlated, and not all of them need to be assessed for selecting wheat genotypes.

Variables related to vegetative indices, gas exchange, chlorophyll, MTG, GY, WUE and iWUE contents were the most important to differentiate genotype groups and water regime levels (Figure 2), as they have weights (length of arrows (≅ 1.0)) greater, indicating that these variables can be used to select wheat genotypes more productive under water stress and are also important components in the process of selecting more productive materials. Root mass (0-20 cm) had a low weight (0.25) for differentiating water regimes and genotype groups.

Without water stress, genotype groups with WR4 (100% of CET replacement) are mostly on the right side of the biplot, meaning that they present higher vegetative indices and photosynthetic activity, promoting higher grain yields (Figure 2).

The groups of genotypes showed linear positive responses of the latent variable (principal coordinate 1) (Table 5) as a function of the water regime. And we noticed that group 3 has a response rate (slope = 0.0326) 18% higher on average. In contrast, the other two groups have similar slopes, probably because the genotypes grouped in 3 are classified as biotypes of irrigated crops, except for PF 100368 and TBIO Sintonia. The changes in water availability lead to significant changes in all traits according to the weights (horizontal arrow length) in the principal coordinate 1 (Figure 3).

3. Discussion

This paper studied the effect of four water regimes in twenty-seven variables and eighteen wheat genotypes. Through the joint multivariate analysis, we could separate the wheat genotypes in three groups (Table 2 and Table 3) through genetic distance studies. This indicates that genetic variability in wheat genotypes among the studied variables is essential for conducting genetic distance studies (Figure 1). In addition, different responses to water regimes also occur among irrigated and rainfed wheat biotypes and other species, such as triticale and common beans [30].

The most similar genotypes were Aliança and BRS 254, and those that differed the most from one another were PF 020037 and PF 120337. Aliança is a rainfed genotype and, in a previous study, was considered the most drought-tolerant genotype in severe water stress [4]. In addition, Tavares et al. [31] studied six soybean genotypes in the Brazilian Cerrado region during the dry season and also showed genotypes adapted to different water regimes.

The study carried out by Soares et al. [4] confirms the clustering results from this work, in which rainfed genotypes (Aliança, BRS 404 and MGS Brilhante) had higher DRI than irrigated genotypes (BRS 394, BRS 254 and BRS 264) after two years of cultivation. This was expected since the genotypes classified in group 1, except for CPAC 01019 and BRS 254, are rainfed biotypes. In contrast, the members of group 3, except for genotypes PF 100368 and TBIO Sintonia, are biotypes for irrigated cultivation, which generally have a low DRI since they were developed for environments without water restriction. The wheat genotypes with higher DRI can be selected as a reference in breeding programs as they were more tolerant to water stress.

Using vegetative indices through cameras coupled in VANTS are quick and non-invasive complementary methods or can be replaced to the traditional selection of genotypes more tolerant to water stress [32, 33, 34].

Vegetative indices based on near-infra-red bands could separate the groups of wheat genotypes and were more affected by water stress, indicating the possibility of using vegetative indices as a quick and non-destructive way to differentiate and select wheat genotypes. Overall, these results are consistent with data from other authors who have found that NIR-based spectral indices have high accuracy and efficiency in selecting more productive wheat genotypes under water stress [28, 35]. In contrast, indices based on the visible range are more strongly correlated with plant growth [20, 36].

In addition, spectral indices collected by cameras coupled to aircraft consistently show a stronger association with grain yield than indices obtained by proximal methods, indicating greater precision, speed, and scale gains [28].

Positive correlations of the NDVI, GNDVI, TCARI and OSAVI vegetation indices and negative correlations between the TCARI/OSAVI index with shoot biomass and grain yield in wheat are reported by Frels et al. [37]. TCARI/OSAVI index is very sensitive to chlorophyll variation and resistant to the effect of soil reflectance and non-photosynthetic matter, distorting the response in relation to the other indices [38]. Non-destructive, simple and fast measurement characteristics of PRI and NDRE constitute an important advantage over the physiological parameter (contents of photosynthetic pigments, chlorophyll fluorescence, and gas exchange rate) and will make it useful in stress detection, especially under severe level or late-stage of heat and water stress [39].

Genotypes under stressed conditions (WR1: 22% CET replacement) have higher values of WUE, iWUE, root mass and TO (TCARI/OSAVI) (Table 2 and Table 3) because there was less water lost, and TO, for instance, is a spectral predictor that is sensitive to chlorophyll variation [40].

The roots contribute only 10 to 20% of the total plant weight, but a well-developed root system is essential for the nutrition of nutrients and water and, therefore, for the growth and final yield of the plants [2]. Under conditions of limited water availability during vegetative growth, such as plants seeking to increase root volume to improve water uptake and avoid yield losses.

Water availability promoted significant alterations in all traits, according to the weights (horizontal arrows length) in the principal coordinate 1 (Figure 2). These observations confirm the findings of other authors who reported that under water deficit conditions, water relations and plant metabolism are impaired [18, 41]. Photosynthesis can be affected by Rubisco activity and limited by chloroplast CO₂ concentration, i.e., under water deficit conditions, plants close stomata to prevent transpiration, which reduces photosynthesis by lowering CO₂ influx [14].

In addition to direct limitation, water stress can cause morphological changes in plants, such as changes in the organization of chloroplasts and the number of pigments, especially chlorophylls and xanthophylls, which consume energy and lead to lower grain yields [9, 15]. Besides, water stress reduces the number and length of the ears of wheat [4]. Abscisic acid (ABA) generally regulates and stimulates reactions such as stomata closure and maintaining water balance and stimulates transcription and activities of antioxidant enzymes under water deficit conditions [42].

In addition to drought tolerance, escape may occur because of some mechanism that prevents drought, including reduced leaf area, chlorophyll content, number of tillering, plant height, and stomatal conductance to avoid water loss; increase root length and root-to-shoot ratio, probably to increase water uptake capacity; increase proline content in leaves and reduce cycling [13]. Besides, water stress reduced the plant cycle in our experiments, similar to Tavares et al. [31] working with soybean.

Other studies show that water stress alters the morphology of the root system of wheat plants. For instance, changes in root angle, primary root length, number of lateral roots [43], mean root diameter [44] and root mass [27] have been reported in the literature.

In the present study, with WR1 (22% CET replacement), a 31% reduction in shoot mass was observed in relation to the non-stressed water regime (WR4 - 100% CET replacement) (Table 3). A decrease in the shoot mass ratio and root mass suggests that wheat plants increase the root system to absorb more water. Furthermore, under water deficit, wheat develops nodal roots with a narrower angle, which tend to grow deeper than those with a larger root angle and increase root density and decrease diameter, as a consequence, roots achieve deeper layers of soil [45].

The water regimes most affected the grain yield, net assimilation rate of CO₂, water use efficiency, and NDVI (Figure 2). The groups of genotypes showed exponential responses of grain yield and net assimilation of CO₂, and a logarithm response of NDVI, water use efficiency and Intrinsic water use efficiency (Table 5).

Increasing water availability increased NDVI and photosynthetic activity, promoting higher grain yield. On a descending scale, we observed that genotypes of groups II and III respond faster to water availability than genotypes of group I for the variables grain yield, net CO₂ assimilation and NDVI.

There is greater water use efficiency under lower water availability. However, this efficiency is reduced more quickly in the genotypes of groups III and II than in the genotypes of group I.

Soares et al. [4] also obtained differences between wheat genotypes for drought tolerance. The number of ears, mass of a thousand grains, grain yield, net CO₂ assimilation and hectoliter weight are important characteristics for selecting more productive genotypes. In addition, the variables DRI and efficiency in water use are important to classify genotypes for drought tolerance.

Water use efficiency is a variable related to the ability of a plant to produce grains with a lower amount of water. Water use efficiency by wheat plants increases under conditions of moderate water stress and decreases under stress-free conditions [15, 16]. Using genotypes with greater water use efficiency [46], especially for species with high water demand cultivated in rainfed agriculture [47], reduces crop losses in areas with water shortages.

4. Materials and Methods

4.1. Experimental design and conducting the experiment

The study was carried out in the wheat cultivation areas in Planaltina, DF, Brazil (15º35'30''S and 47º42'30''W, an altitude of 1006 m) between May and September 2018 and 2019. The climate is classified as Aw (Koeppen-Geiger) - Tropical, with rainfall concentrated in summer (October to April) and the dry period during the winter (May to September). The annual rainfall is between 1200 and 1500 mm. Precipitation and air temperature data from the experimental area in 2018 and 2019 are presented in Figure 4 and were obtained from a meteorological station close to the experiment. The experimental area had been cultivated for the last four years with wheat under different water regimes in winter and fallow in summer.

The soil is classified as Oxisol according to [48]. Before the installation of the experiment, soil samples were collected at a depth of 0-20 cm, with the following physicochemical properties: pH (CaCl₂) of 5.7, 11 mg dm^-3 P, 186 mg dm^-3 K, 5.77 cmolc dm^-3 Ca, 1.83 cmolc dm^-3 Mg, 0.02 cmolc dm^-3 Al, 26.7 g kg^-1 organic matter. The soil texture comprised 46, 10 and 44% of clay, silt, and sand, respectively.

The model fitted the soil water retention curve [49]. The following estimates were obtained: Residual water content (θs) 0.0839 cm³ cm^-3, saturated water content (θs) 0.5500 cm³ cm^-3, and parameters α (1.892 kPa^-1) and n (1.2390). The field capacity moisture was 0.3423 cm³ cm^-3.

The experimental design consisted of randomized blocks in a split-plot scheme with three replicates. The plots consisted of eighteen wheat genotypes: BRS 254; BRS 264; CPAC 01019; CPAC 01047; CPAC 07258; CPAC 08318; CPAC 9110; BRS 394 (irrigated biotypes) and Aliança; BR 18_Terena; BRS 404; MGS Brilhante; PF 020037; PF 020062; PF 120337; PF 100368; PF 080492; TBIO Sintonia (rainfed biotypes). Subplots were composed of four water regimes (WR). In 2018, 123.68 mm, 241.74 mm, 455.38 mm, and 562.2 mm were applied during the cycle, corresponding to WR1, WR2, WR3, and WR4, respectively. In 2019, 119.1 mm, 232.15 mm, 438.55 mm, and 541.43 mm were applied, corresponding to WR1, WR2, WR3, and WR4, respectively. In 2018 and 2019, the irrigation regimes used corresponded to 22%, 43%, 81% and 100% of crop evapotranspiration (CET) replacement, respectively.

The rainfed biotypes were MGS Brilhante, Aliança, TBIO Sintonia and BR18_Terena (all traditional material), and BRS 404 (registered in 2014). The rainfed lines developed for the Cerrado region were PF020037 [(strong wax formation on leaves and stems] and PF020062 (but without wax formation). In addition, the rainfed biotype PF080492 was adapted to the Brazilian Cerrado region [4]. The irrigated cultivars were BRS 254, BRS 264 (traditional materials) and BRS 394 (registered in 2014) and in the experimental phase for commercialization, the irrigated biotypes CPAC 01019; CPAC 01047; CPAC 07258; CPAC 08318 and CPAC 9110.

The highest-level irrigation was carried out according to the Cerrado Irrigation Monitoring Program [50], that is, by replacing evapotranspiration using regional agrometeorological indicators, soil texture, and date of total plant emergence date of the plants and estimating the reference evapotranspiration with based on the equation proposed by Penman-Monteith [51]. Irrigation was carried out approximately every five days, according to weather conditions and the phenological phase of the plants. The amount of water applied in each irrigation was estimated by collectors placed parallel to the irrigation bar.

Twenty days before planting, glyphosate was applied (1440 g.e.a ha^-1). Seeding was done mechanically on May 22, 2018, and May 23, 2019, in a no-till system with 90 seeds per meter at a depth of 3 cm. Fertilization was carried in the furrows at 400 kg ha-1 of 04-30-16 (N, P₂O5, K₂O), and at 25 days after wheat emergence (DAE), nitrogen was applied at a dose of 90 kg ha^-1 N, as urea.

Furthermore, at 36 and 38 DAE, in 2018 and 2019, the growth reducer trinexapac-ethyl at a dose of 125 g ha^-1 was applied in the first visible node and with the second palpable node. In addition, the herbicide metsulfuron-methyl was applied at 4 g ha^-1 in both years for weed control on 15 DAE.

In both experiments, a homogeneous water layer was applied during the first 35 DAE at the tillering stage. An average of 150 mm of water was applied to obtain a homogeneous plant stand. After this period, the “line source” method was applied modified by Jayme-Oliveira et al. [52]. The water regime (WR) was obtained by an irrigation bar (IrrigaBrasil model 36/42) 20 m wide on each side of the bar. The bar was connected to a self-propelled TurboMaq 75/GB) with controlled speed.

Each experimental unit consisted of a genotype with 18.0 m formed by eight cultivation rows spaced 0.17 m apart. Each water regime consisted of an experimental subunit 2.0 m in length formed by eight rows 0.17 m apart. The usable area consisted of the six central rows, omitting the margins and 2 m on each side.

4.2. Analyzed variables

At the flowering stage (75 DAE,) gas exchange and electron transport rate were evaluated from 8:30 am to 12:30 pm at an irradiance of 1200 μmol photons m^-2s^-1 and an external CO₂ concentration (Ca) of 400 μmol mol^-1. For each subplot, three fully expanded flag leaves were used to evaluate gas exchange using a portable open-flow gas exchange system (IRGA, model LI-6400xt LI-COR Inc., Lincoln, NE). In addition, chlorophyll fluorescence and maximum quantum yield of photosystem II (Fv/Fm) were evaluated with a modulated portable fluorometer coupled to IRGA. Evaluations were performed on dark-adapted leaves for at least three hours after 10:30 pm. During this period, initial fluorescence (F0), maximum fluorescence (Fm), and potential quantum yield of photosystem II [Fv/Fm = (F0-Fm)/Fm)] were estimated [53].

Chlorophyll a and b content were measured with the ClorofiLOG portable chlorophyll meter (CFL-1030, Falker), which provides relative measurements (0 to 100) of total chlorophyll but is linearly correlated with total chlorophyll content [54].

Using a multispectral camera, model Micasense RedEdge, which captures images in five different spectral ranges, that is, range: 465-485 nm; range: 550-570 nm; range: 663-673 nm; range: 712-722 nm; range: 820-860 nm, coupled to an unmanned aerial vehicle (UAV) with a rotary-wing, images were taken to estimate the vegetation indices in the flowering phase of the crop. This camera features an optical resolution of 1280 x 960 pixels, with images recorded in RAW12 bits. The flight took place at a height of 45 m at 10 am. The reflectance maps were calculated by generating mosaics in the Pix4D Mapper software (v5.4.6, Pix4D, Lausanne, Switzerland).

The maps were processed in R software using the raster package, and the following vegetative indices were extracted as in Tavares et al. [31]: Normalized Difference Vegetation Index - NDVI; Green Normalized Difference Vegetation Index - GNDVI; Green-Red Vegetation Index - GRVI; Difference Vegetation Index - DVI; Normalized Difference Red Edge - NDRE; Soil Adjusted Vegetation Index - SAVI; Photochemical Physiological Reflectance Index - PRI; Optimized Soil Adjusted Vegetation Index - OSAVI; Chlorophyll Absorption and Reflectance Index - TCARI; TCARI/OSAVI - TO ratio.

Soil samples were collected to quantify the root mass of wheat genotypes during flowering. An auger with sharp edges, an internal diameter of 9.8 cm, and a length of 20 cm was used for root sampling. For each experimental subunit, 1508 cm³ of soil with roots was sampled, and the auger was inserted over the wheat crop line. The aerial part of the wheat plants corresponded to the area of the auger where the soil samples with roots (9.8 linear cm) were removed, cut close to the soil with a knife, and placed in a paper bag. The soil auger is made of stainless steel and has an iron shaft to facilitate rotation and application of force. This sampling procedure was adopted as in [55].

The root samples with soil were kept in plastic bags, and the shoot in paper bags. The soil mass with roots was determined, and a subsample of 100 g was used to determine soil moisture content. Samples were stored in cold rooms at a temperature of -5°C possible until the roots were separated from the soil to avoid loss of mass or drying the roots. Roots were separated from soil by washing with water three times in a 500 µm sieve and later separated from other organic materials. The roots and shoots were kept in a convection oven at 60ºC for 72 hours to determine the dry matter. The results were expressed as root dry mass to soil mass ratio and root to shoot ratio.

At harvest, grain yield (GY) on the usable area of each experimental subunit and mass of thousand grains (MTG) was measured and standardized to grain moisture content of 13% on a wet basis. Water use efficiency agronomic (WUE) was calculated using the ratio of grain yield to crop water requirement [56]. The intrinsic water-use efficiency (iWUE) was calculated by the ratio between the net assimilation of CO₂ and stomatal conductance (A/gs). The drought resistance index was calculated according to [57] Fischer and Maurer (1978).

4.3. Statistical Analysis

The data were subjected to joint multivariate analysis of variance by harvest based on singular value decomposition (SVD). Genotypes were grouped based on Mahalanobis distance using Ward's method. The Mojena criterion [58] was used to define the cutoff point in the dendrogram, and the relative importance (proportion) of variables in the distance between genotypes was determined by the Singh [59] criterion. In both years, treatments (combinations of genotypes and water regime levels) were analyzed graphically in a biplot [60]. This allows visualizing the relationship between genotypes and treatments. Pearson correlation analysis (t-test, p<0.05) was performed with the residuals. The statistical analyses were performed using the R v3.6.1. software.

5. Conclusions

High-throughput and non-destructive vegetative indices based on the near-infrared band can be used to detect drought tolerance of wheat genotypes, as they are correlated with physiological and agronomic traits. The wheat genotypes can be clustered according to their responses to water stress. The higher drought resistance index genotypes are Aliança; BRS 254; BRS 404; CPAC 01019; PF 020062; and PF 080492. On the other hand, the genotypes BRS 264; BRS 394; CPAC 01047; CPAC 07258; CPAC 8318; CPAC 9110; PF 100368; and TBIO Sintonia are the most affected by water availability. There is less variability in wheat genotypes' physiological, agronomic and spectral responses when providing water based on 100% evapotranspiration.