Preprint
Article

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

Altmetrics

Downloads

113

Views

69

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

15 September 2023

Posted:

19 September 2023

You are already at the latest version

Alerts
Abstract
Selecting drought-tolerant and more water-efficient wheat genotypes is a research priority, spe-cifically in regions with irregular rainfall or areas where climate change is expected to result in reduced water availability. The objective of this work was to use high-throughput measurements with morphophysiological traits to characterize wheat genotypes in relation to water stress. Field experiments were conducted from May to September 2018 and 2019, using a sprinkler bar irriga-tion system to control water availability to 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). The water regimes varied from 22 to 100% of the crop evapotranspiration replacement. Water stress negatively affected gas exchange, vegetation indices and grain yield. High throughput variables TCARI, NDVI, OSAVI, SAVI, PRI, NDRE and GNDVI had higher correlation with yield and morphophysiological measurements. The drought re-sistance index indicated that genotypes Aliança BRS 254, BRS 404, CPAC 01019, PF 020062, and PF 080492 were more drought tolerant.
Keywords: 
Subject: Biology and Life Sciences  -   Agricultural Science and Agronomy

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, gs, 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 (gs) 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 CO2 concentration, i.e., under water deficit conditions, plants close stomata to prevent transpiration, which reduces photosynthesis by lowering CO2 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 CO2, water use efficiency, and NDVI (Figure 2). The groups of genotypes showed exponential responses of grain yield and net assimilation of CO2, 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 CO2 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 CO2 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 (CaCl2) 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 cm3 cm-3, saturated water content (θs) 0.5500 cm3 cm-3, and parameters α (1.892 kPa-1) and n (1.2390). The field capacity moisture was 0.3423 cm3 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, P2O5, K2O), 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 CO2 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 cm3 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 CO2 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.

Author Contributions

Conceptualization: Ribeiro Júnior, W.Q.; Ramos, M.L.G. Metodology: Ribeiro Júnior, W.Q.; Tavares, C.J.; Pereira, L.F.; Casari, R.A.C.; Sousa, C.A.F. Formal analysis: Tavares, C.J.; Ribeiro Júnior, W.Q.; Ramos, M.L.G.; Silva, A.R.; Muller, O.; Pereira, L.F.; Casari, R.A.C. Investigation: Ribeiro Júnior, W.Q.; Ramos, M.L.G.; Sousa, C.A.F.; Pereira, L.F. Writing and editing: Ribeiro Júnior, W.Q.; Ramos, M.L.G.; Tavares, C.J.; Silva, A.R; Muller, O.; Silva, A.R.

Funding

This research received no external funding.

Acknowledgments

The research group of Embrapa Cerrados, the University of Brasília, the Institute for Bio-and Geosciences (Julich), the Federal Institute Goiano for the partial release and support for doctoral-level training of the first author, and the National Council for Scientific Development and Technology – CNPq for for the scientific productivity fellowships granted to the third and eighth authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zörb, C.; Ludewig, U.; Hawkesford, M.J. Perspective on wheat yield and quality with reduced nitrogen supply. Trends Plant Sci. 2018, 23, 1029–1037. [Google Scholar] [CrossRef]
  2. Liu, W.; Wang, J.; Wang, C.; Ma, G.; Wei, Q.; Lu, H.; Xie, Y.; Ma, D.; Kang, G. Root growth, water and nitrogen use efficiencies in winter wheat under different irrigation and nitrogen regimes in North China Plain. Front. Plant Sci. 2018, 9. [Google Scholar] [CrossRef]
  3. Soares, D.C.; Bonfim-Silva, E.M.; Silva, T.J.A.; Anicésio, E.C.A.; Duarte, T.F.; Oliveira, J.F. Growth and production of wheat cultivars under water tensions in Cerrado soil. Rev. Bras. Eng. Agríc. Ambiental. 2023, 27, 279–286. [Google Scholar] [CrossRef]
  4. Soares, G.F.; Ribeiro Júnior, W.Q.; Pereira, L.F.; Lima, C.A.; Soares, D.S.; Muller, O.; Rascher, U.; Ramos, M.L.G. Characterization of wheat genotypes for drought tolerance and water use efficiency. Scient. Agricola. 2021, 78, e20190304. [Google Scholar] [CrossRef]
  5. Pereira, J.F.; Cunha, G.R.; Moresco, E.R. Improved drought tolerance in wheat is required to unlock the production potential of the Brazilian Cerrado. Crop. Breed. Appl. Biotechnol. 2019, 19, 217–225. [Google Scholar] [CrossRef]
  6. Cattelan, A.J.; Dall’agnol, A. The rapid soybean growth in Brazil. OCL. 2018, 25, 1–12. [Google Scholar] [CrossRef]
  7. Schauberger, B.; Rolinski, S.; Schaphoff, S.; Muller, C. Global historical soybean and wheat yield loss estimates from ozone pollution considering water and temperature as modifying effects. Agric. For. Meteorol. 2019, 265, 1–15. [Google Scholar] [CrossRef]
  8. Peña-Gallardo, M.; Vicente-Serrano, S.M.; Quiringb, S.; Svoboda, M.; Hannafordd, M.; Hannafordd, J.; Tomas-Burguera, M.; Martín-Hernándeza, N.; Domínguez-Castro, F.; Kenawy, E.A. Response of crop yield to different timescales of drought in the United States: Spatio-temporal patterns and climatic and environmental driver. Agric. For. Meteorol. 2019, 264, 40–55. [Google Scholar] [CrossRef]
  9. Ahmed, H.G.M.-D.; Sajjad, M.; Li, M.; Azmat, M.A.; Rizwan, M.; Maqsood, R.H.; Khan, S.H. Selection Criteria for Drought-Tolerant Bread Wheat Genotypes at Seedling Stage. Sustainability 2019, 11, 2584. [Google Scholar] [CrossRef]
  10. Conab. (2023). Companhia Nacional de Abastecimento. Acompanhamento da safra brasileira de grãos. Safra 2022/2023, n. 1 – Boletim de monitoramento agrícola – cultivos de verão e inverno, p. 1-20, 2023. Available online: https://www.conab.gov.br/info-agro/safras/graos/monitoramento-agricola (accessed on 20 June 2023).
  11. Hatfield, J.L. Increased temperatures have dramatic effects on growth and grain yield of three maize hybrids. Agric. Environ. Lett. 2016, 1. [Google Scholar] [CrossRef]
  12. Wang, L.; Wang, S.; Chen, W.; Li, H.; Deng, X. Physiological mechanisms contributing to increased water-use efficiency in winter wheat under organic fertilization. PLoS ONE. 2017, 12, 1–21. [Google Scholar] [CrossRef]
  13. Li, P.; Ma, B.; Palta, J.A.; Ding, T.; Cheng, Z.; Lv, G.; Xiong, X. Wheat breeding highlights drought tolerance while ignores the advantages of drought avoidance: A meta-analysis. Eur. J. Agron. 2021, 122, 126196. [Google Scholar] [CrossRef]
  14. Ahmed, H.G.M.-D.; Zeng, Y.; Yang, X.; Anwaar, H.A.; Mansha, M.Z.; Hanif, C.M.S.; Ikram, K.; Ullah, A.; Alghanem, S.M.S. Conferring drought-tolerant wheat genotypes through morpho-physiological and chlorophyll indices at seedling stage. Saudi J. Biol. Sci. 2020, 27, 2116–2123. [Google Scholar] [CrossRef]
  15. Pour-Aboughadareh, A.; Mohammadi, R.; Etminan, A.; Shooshtari, L.; Maleki-Tabrizi, N.; Poczai, P. Effects of drought stress on some agronomic and morpho-physiological traits in Durum wheat genotypes. Sustainability 2020, 12, 5610. [Google Scholar] [CrossRef]
  16. Zhao, W.; Liu, L.; Shen, Q.; Yang, J.; Han, X.; Tian, F.; Wu, J. Effects of water stress on photosynthesis, yield, and water use efficiency in winter wheat. Water. 2020, 12, 1–19. [Google Scholar] [CrossRef]
  17. Anwaar, H.A.; Perveen, R.; Mansha, M.Z.; Abid, M.; Sarwar, M.; Aatif, H.M.; Umar, U.U.D.; Sajid, M.; Aslam, H.M.U.; Alam, M.M.; Rizwan, M.; Ikram, R.M.; Alghanem, S.M.S.; Rashid, A.; Khan, K.A. Assessment of grain yield indices in response to drought stress in wheat (Triticum aestivum L.). Saudi J. Biol. Sci. 2020, 27, 1818–1823. [Google Scholar] [CrossRef] [PubMed]
  18. Ouyang, W.; Struik, P.C.; Yin, X.; Yang, J. Stomatal conductance, mesophyll conductance, and transpiration efficiency in relation to leaf anatomy in rice and wheat genotypes under drought. J. Exp. Bot. 2017, 68, 5191–5205. [Google Scholar] [CrossRef]
  19. Anderegg, J.; Yu, K.; Aasen, H.; Walter, A.; Liebisch, F.; Hund, A. Spectral vegetation indices to track senescence dynamics in diverse wheat germplasm. Front. Plant Sci. 2020, 10, 1749. [Google Scholar] [CrossRef] [PubMed]
  20. Yang, M.; Hassan, M.A.; Xu, K.; Zheng, C.; Rasheed, A.; Zhang, Y.; Jin, X.; Xia, X.; Xiao, Y.; He, Z. Assessment of water and nitrogen use efficiencies through uav-based multispectral phenotyping in winter wheat. Front. Plant Sci. 2020, 11, 927. [Google Scholar] [CrossRef] [PubMed]
  21. El-Hendawy, S.E.; Al-Suhaibani, N.A.; Elsayed, S.; Hassan, W.M.; Dewir, Y.H.; Refay, Y.; Abdella, K.A. Potential of the existing and novel spectral reflectance indices for estimating the leaf water status and grain yield of spring wheat exposed to different irrigation rates. Agri. Water Manag. 2019, 217, 356–37. [Google Scholar] [CrossRef]
  22. Qiu, R.C.; Wei, S.; Zhang, M.; Sun, H.; Li, H.; Liu, G.; Li, M. Sensors for measuring plant phenotyping: A review. Int. J. Agric. Biol. Eng. 2018, 11, 1–17. [Google Scholar] [CrossRef]
  23. Sobejano-Paz, V.; Mikkelsen, T.N.; Baum, A.; Mo, X.; Liu, S.; Köppl, C.J.; Johnson, M.S.; Gulyas, L.; García, M. Hyperspectral and Thermal Sensing of Stomatal Conductance, Transpiration, and Photosynthesis for Soybean and Maize under Drought. Remote Sens. 2020, 12, 3182. [Google Scholar] [CrossRef]
  24. Andrade-Sanchez, P.; Gore, M.A.; Heun, J.T.; Thorp, K.R.; Carmo-Silva, A.E.; French, A.N.; SalvuccI, M.E.; White, J.W. Development and evaluation of a field-based high-throughput phenotyping platform. Funct. Plant Biol. 2014, 41, 68–79. [Google Scholar] [CrossRef] [PubMed]
  25. Haghighattalab, A.; Gonzalez, L.; Mondal, S.; Singh, D.; Schinstock, D.; Rutkoski, J.; Ortiz-Monasterio, I.; Singh, R.; Goodin, D.; Poland, J. Application of unmanned aerial systems for high throughput phenotyping of large wheat breeding nurseries. Plant methods. 2016, 35. [Google Scholar] [CrossRef] [PubMed]
  26. Keller, B.; Vass, I.; Matsubara, S.; Paul, K.; Jedmowski, C.; Pieruschka, R.; Nedbal, L.; Rascher, U.; Muller, O. Maximum fluorescence and electron transport kinetics determined by light-induced fluorescence transients (LIFT) for photosynthesis phenotyping. Photosynth. Res. 2018, 140, 221–233. [Google Scholar] [CrossRef] [PubMed]
  27. Khadka, K.; Earl, H.J.; Raizada, M.N.; Navabi, A. A Physio-Morphological Trait-Based Approach for Breeding Drought Tolerant Wheat. Front. Plant Sci. 2020, 11, 715. [Google Scholar] [CrossRef] [PubMed]
  28. Hu, Y.; Knapp, S.; Schmidhalter, U. Advancing High-Throughput Phenotyping of Wheat in Early Selection Cycles. Remote Sens. 2020, 12, 574. [Google Scholar] [CrossRef]
  29. Yang, M.; Hassan, M.A.; Xu, K.; Zheng, C.; Rasheed, A.; Zhang, Y.; Jin, X.; Xia, X.; Xiao, Y.; He, Z. Assessment of water and nitrogen use efficiencies through uav-based multispectral phenotyping in winter wheat. Front. Plant Sci. 2020, 11, 927. [Google Scholar] [CrossRef] [PubMed]
  30. Silva, A.D.N.; Ramos, M.L.G.; Ribeiro Junior, W.Q.; Alencar, E.R.; Silva, P.C.; Lima, C.A.; Vinson, C.C.; Silva, M.A.V. Water stress alters physical and Chemical quality ingrains of common bean, triticale and wheat. Agric. Water Manag. 2020, 231, 16023. [Google Scholar] [CrossRef]
  31. Tavares, C.J.; Ribeiro Junior, W.Q.; Ramos, M.L.G.; Pereira, L.F.; Casari, R.A.C.N.; Pereira, A.F.; de Sousa, C.A.F.; Silva, A.R.; Silva Neto, S.P.; Mertz-Henning, L.M. Water Stress Alters Morphophysiological, Grain Quality and Vegetation Indices of Soybean Cultivars. Plants. 2022, 11, 559. [Google Scholar] [CrossRef]
  32. Casari, R.A.C.N.; Paiva, D.S.; Silva, V.N.B.; Ferreira, T.M.M.; Souza, M.T.; Oliveira, N.G.; Kobayashi, A.; Molinari, H.B.C.; Santos, T.T.; Gomide, R.L.; Magalhães, P.C.; Sousa, C.A.F. Using thermography to confirm genotypic variation for drought response in maize. Int. J. Mol. Sci. 2019, 20, 273. [Google Scholar] [CrossRef]
  33. Vitek, P.; Veselá, V.; Klem, K. Spatial and temporal variability of plant leaf responses cascade after PSII inhibition: raman, chlorophyll fluorescence and infrared thermal imaging. Sensors. 2020, 20, 1015. [Google Scholar] [CrossRef] [PubMed]
  34. Silva, A.D.N.; Ramos, M.L.G.; Ribeiro Junior, W.Q.; da Silva, P.C.; Soares, G.F.; Casari, R.A.D.C.N.; Vinson, C.C. Use of Thermography to Evaluate Alternative Crops for Off-Season in the Cerrado Region. Plants. 2023, 12, 2081. [Google Scholar] [CrossRef] [PubMed]
  35. Bowman, B.; Chen, J.; Zhang, J.; Wheeler, J.; Wang, Y.; Zhao, W.; Nayak, S.; Heslot, N.; Bockelman, H.; Bonman, J. Evaluating grain yield in spring wheat with canopy spectral reflectance. Crop Sci. 2015, 55, 1881–1890. [Google Scholar] [CrossRef]
  36. Prey, L.; Schmidhalter, U. Simulation of satellite reflectance data using high-frequency ground based hyperspectral canopy measurements for in-season estimation of grain yield and grain nitrogen status in winter wheat. ISPRS J. Photogramm. Remot. 2019, 149, 176–187. [Google Scholar] [CrossRef]
  37. Frels, K.; Guttieri, M.; Joyce, B.; Leavitt, B.; Baenziger, P.S. Evaluating canopy spectral reflectance vegetation indices to estimate nitrogen use traits in hard winter wheat. Field Crops Res. 2018, 217, 82–92. [Google Scholar] [CrossRef]
  38. Damm, A.; Paul-Limoges, E.; Haghighi, E.; Simmer, C.; Morsdorf, F.; Schneider, F.D.; Van Der Tol, C.; Migliavacca, M.; Rascher, U. Remote sensing of plant-water relations, An overview and future perspectives. J. Plant Physiol. 2018, 227, 3–19. [Google Scholar] [CrossRef]
  39. Cao, Z.; Yao, X.; Liu, H.; Liu, B.; Cheng, T.; Tian, Y.; Cao, W.; Zhu, Y. Comparison of the abilities of vegetation indices and photosynthetic parameters to detect heat stress in wheat. Agric. For. Meteorol. 2019, 265, 121–136. [Google Scholar] [CrossRef]
  40. Ballester, C.; Zarco-Tejada, P.J.; Nicolás, E.; Alarcón, J.J.; Fereres, E.; Intrigliolo, D.S.; Gonzalez-Dugo, V. Evaluating the performance of xanthophyll, chlorophyll and structure-sensitive spectral indices to detect water stress in five fruit tree species. Precis. Agric. 2018, 19, 178–193. [Google Scholar] [CrossRef]
  41. Marček, T.; Hamow, K.; Végh, B.; Janda, T.; Darko, E. Metabolic response to drought in six winter wheat genotypes. PLoS ONE. 2019, 14, 1–23. [Google Scholar] [CrossRef]
  42. Wang, X.; Li, Q.; Xie, J.; Huang, M.; Cai, J.; Zhou, Q.; Dai, T.; Jiang, D. Abscisic acid and jasmonic acid are involved in drought priming-induced tolerance to drought in wheat. Crop J. 2021, 9, 120–132. [Google Scholar] [CrossRef]
  43. Zhan, A.; Lynch, J. P. Reduced frequency of lateral root branching improves N capture from low-N soils in maize. J. Exp. Bot. 2015, 66, 2055–2065. [Google Scholar] [CrossRef] [PubMed]
  44. Haling, R.E.; Brown, L.K.; Bengough, A.G.; Young, I.M.; Hallett, P.D.; White, P.J. Root hairs improve root penetration, root – soil contact, and phosphorus acquisition in soils of different strength. J. Exp. Bot. 2013, 64, 3711–3721. [Google Scholar] [CrossRef] [PubMed]
  45. Chen, X.; Ding, Q.; Blłaszkiewicz, Z.; Sun, J.; Sun, Q.; He, R.; Li, Y. Phenotyping for the dynamics of field wheat root system architecture. Sci. Rep. 2017, 7, 37649. [Google Scholar] [CrossRef]
  46. Yasi, H.; Attia, H.; Alamer, K.; Hassan, F.; Ali, E.; Elshazly, S.; Siddique, K.H.M.; Hessini, K. Impact of drought on growth, photosynthesis, osmotic adjustment, and cell wall elasticity in Damask rose. Plant Physiol. Biochem. 2020, 150, 133–139. [Google Scholar] [CrossRef]
  47. Tavera, V.C.; Mancilla, C.L.A.; Gallegos, J.A.C.; Pimentel, J.G.R.; Arriaga, A.I.M.; Ruiz Nietto, J.E. Mechanisms for water use efficiency between bean cultivars tolerant to drought are different. Acta Sci. Agron. 2018, 40, e39378. [Google Scholar] [CrossRef]
  48. Soil Survey Staff. 2014. ‘Keys to Soil Taxonomy’. In: Burt, R., Soil Survey Staff (Eds.), Soil Survey Field and Laboratory Methods Manual. Soil Survey Investigations Report No. 51, Version 2.0. U.S. Department of Agriculture, Natural Resources Conservation Service. Washington, DC.
  49. Van Genuchten, M.T. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. 1980, 44, 892–898. [Google Scholar] [CrossRef]
  50. Embrapa. 2011. Programa de Monitoramento da Irrigação. EMBRAPA, Brasília.
  51. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop evapotranspiration: guidelines for computing crop water requirements. FAO Irrigation and Drainage Paper, 1998, 56, Rome, p. 300.
  52. Jayme-Oliveira, A.; Ribeiro Junior, W.Q.; Ramos, M.L.G.; Ziviani, A.C.; Jakelaitis, A. Amaranth, quinoa, and millet growth and development under different water regimes in the Brazilian Cerrado. Pesqui. Agropecu. Bras. 2017, 52, 561–571. [Google Scholar] [CrossRef]
  53. Genty, B.; Briantais, J.M.; Baker, N.R. The relationship between the quantum yield of photosynthetic electron transport and quenching of chlorophyll fluorescence. Biochim Biophys Acta. 1989, 990, 87–92. [Google Scholar] [CrossRef]
  54. Rigon, J.P.G.; Capuani, S.; Beltrão, N.M.; Brito Neto, J.F.; Sofiatti, V.; França, V.F. Non-destructive determination of photosynthetic pigments in the leaves of castor oil plants. Acta Sci. Agron. 2012, 34, 325–329. [Google Scholar] [CrossRef]
  55. Ratke, R.F.; Pereira, H.S.; Santos Junior, D.G.; Frazão, J.J.; Barbosa, J.M.; Dias, B.O. Root growth, nutrition and yield of maize with applied different limestone particle size in the cerrado soil. Am. J. Plant Sci. 2014, 5, 463–472. [Google Scholar] [CrossRef]
  56. Hatfield, J.L.; Sauer, T.J.; Prueger, J.H. Managing soils to achieve greater water use efficiency: a review. Agronomy. 2001, 93, 271–280. [Google Scholar] [CrossRef]
  57. Fischer, R.A.; Maurer, R. Drought resistance in spring wheat cultivars. I. Grain yield responses. Aust. J. Agric. Res. 1978, 29, 897–912. [Google Scholar] [CrossRef]
  58. Mojena, R. Hierarchical grouping methods and stopping rules: an evaluation. J. Comput. 1977, 20, 359–363. [Google Scholar] [CrossRef]
  59. Singh, D. The relative importance of characters affecting genetic divergence. Indian J. Genet. Plant Breed. 1981, 41, 237–245. [Google Scholar]
  60. Gabriel, K.R. The biplot graphical display of matrices with application to principal component analysis. Biometrika. 1971, 58, 453–467. [Google Scholar] [CrossRef]
Figure 1. Dendrogram of eighteen wheat genotypes, based on Ward's clustering using the generalized Mahalanobis distance. 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 Sintonia.
Figure 1. Dendrogram of eighteen wheat genotypes, based on Ward's clustering using the generalized Mahalanobis distance. 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 Sintonia.
Preprints 85251 g001
Figure 2. Biplot for scores of combined levels of water regimes (22, 43, 81 and 100% CET replacement, respectively) and wheat genotype groups (1, 2, 3), based on the scores of principal coordinates from the variables: Root mass (RM); Shoot mass (SM); Ratio RM/SM (RM.SM); Grain humidity (GH); mass of thousand grain (MTG); Grain yield (GY); Water use efficiency (WUE); Intrinsic water use efficiency (iWUE - A/gs: net assimilation of CO2/stomatal conductance); The Transformed Chlorophyll Absorption and Reflectance Index (TCARI); Green red vegetation index (GRVI); Normalized difference vegetation index (NDVI); Optimized Soil Adjusted Vegetation Index (OSAVI); Soil adjusted vegetation index (SAVI); Photochemical Physiological Reflectance Index (PRI); Green normalized difference vegetation index (GNDVI); Red Edge Normalized Difference (NDRE); Difference vegetation index (DVI); Net Assimilation of CO2 (A); Stomatal Conductance (gs); Transpiration (E); Internal CO2 concentrations (Ci); Electron transport rate (ETR); Maximum Quantum Yield of Photosystem II (Fv/Fm); Effective Quantum Yield of Photosystem (Fv’/Fm’); Chlorophyll a (Cha); Chlorophyll b (Chb); Ratio Chlorophyll a/b (Cha.b). 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 Sintonia.
Figure 2. Biplot for scores of combined levels of water regimes (22, 43, 81 and 100% CET replacement, respectively) and wheat genotype groups (1, 2, 3), based on the scores of principal coordinates from the variables: Root mass (RM); Shoot mass (SM); Ratio RM/SM (RM.SM); Grain humidity (GH); mass of thousand grain (MTG); Grain yield (GY); Water use efficiency (WUE); Intrinsic water use efficiency (iWUE - A/gs: net assimilation of CO2/stomatal conductance); The Transformed Chlorophyll Absorption and Reflectance Index (TCARI); Green red vegetation index (GRVI); Normalized difference vegetation index (NDVI); Optimized Soil Adjusted Vegetation Index (OSAVI); Soil adjusted vegetation index (SAVI); Photochemical Physiological Reflectance Index (PRI); Green normalized difference vegetation index (GNDVI); Red Edge Normalized Difference (NDRE); Difference vegetation index (DVI); Net Assimilation of CO2 (A); Stomatal Conductance (gs); Transpiration (E); Internal CO2 concentrations (Ci); Electron transport rate (ETR); Maximum Quantum Yield of Photosystem II (Fv/Fm); Effective Quantum Yield of Photosystem (Fv’/Fm’); Chlorophyll a (Cha); Chlorophyll b (Chb); Ratio Chlorophyll a/b (Cha.b). 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 Sintonia.
Preprints 85251 g002
Figure 3. Pearson correlogram between variables representing physiological, spectral and agronomic as a function of water regimes (22, 43, 81 and 100% CET replacement - WR1, WR2, WR3, and WR4, respectively) and wheat genotype groups (1, 2, 3), based on the scores of principal coordinates from the variables: Ratio RM/SM (RM.SM); Ratio Chlorophyll a/b (Cha.b); Water use efficiency (WUE); TCARI/OSAVI Ratio (TO); Intrinsic water use efficiency (iWUE - A/gs: net assimilation of CO2/stomatal conductance); Electron transport rate (ETR); ); Maximum Quantum Yield of Photosystem II (Fv/Fm); Grain humidity (GH); Chlorophyll b (Chb); Chlorophyll a (Cha); Chorophyll total (Cht); Drought Resistance Index (DRI); Root mass (RM); Internal CO2 concentrations (Ci); Mass of thousand grain (MTG); Shoot mass (SM); Effective Quantum Yield of Photosystem (Fv’/Fm’); Stomatal Conductance (gs); Transpiration (E); Grain yield (GY); Net Assimilation of CO2 (A); The Transformed Chlorophyll Absorption and Reflectance Index (TCARI); Green red vegetation index (GRVI); Normalized difference vegetation index (NDVI); Optimized Soil Adjusted Vegetation Index (OSAVI); Soil adjusted vegetation index (SAVI); Photochemical Physiological Reflectance Index (PRI); Green normalized difference vegetation index (GNDVI); Red Edge Normalized Difference (NDRE); Difference vegetation index (DVI); Normalized Differential Greenness Index (NDGI); Ratio Vegetation Index (RVI). 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 Sintonia.
Figure 3. Pearson correlogram between variables representing physiological, spectral and agronomic as a function of water regimes (22, 43, 81 and 100% CET replacement - WR1, WR2, WR3, and WR4, respectively) and wheat genotype groups (1, 2, 3), based on the scores of principal coordinates from the variables: Ratio RM/SM (RM.SM); Ratio Chlorophyll a/b (Cha.b); Water use efficiency (WUE); TCARI/OSAVI Ratio (TO); Intrinsic water use efficiency (iWUE - A/gs: net assimilation of CO2/stomatal conductance); Electron transport rate (ETR); ); Maximum Quantum Yield of Photosystem II (Fv/Fm); Grain humidity (GH); Chlorophyll b (Chb); Chlorophyll a (Cha); Chorophyll total (Cht); Drought Resistance Index (DRI); Root mass (RM); Internal CO2 concentrations (Ci); Mass of thousand grain (MTG); Shoot mass (SM); Effective Quantum Yield of Photosystem (Fv’/Fm’); Stomatal Conductance (gs); Transpiration (E); Grain yield (GY); Net Assimilation of CO2 (A); The Transformed Chlorophyll Absorption and Reflectance Index (TCARI); Green red vegetation index (GRVI); Normalized difference vegetation index (NDVI); Optimized Soil Adjusted Vegetation Index (OSAVI); Soil adjusted vegetation index (SAVI); Photochemical Physiological Reflectance Index (PRI); Green normalized difference vegetation index (GNDVI); Red Edge Normalized Difference (NDRE); Difference vegetation index (DVI); Normalized Differential Greenness Index (NDGI); Ratio Vegetation Index (RVI). 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 Sintonia.
Preprints 85251 g003
Figure 4. Precipitation and average temperature measured by an automatic weather station near the experiment in 2018 and 2019, and a historical series (1974 - 2019).
Figure 4. Precipitation and average temperature measured by an automatic weather station near the experiment in 2018 and 2019, and a historical series (1974 - 2019).
Preprints 85251 g004
Table 1. Summary of Multivariate Analysis of Variance (MANOVA) joint with all interactions.
Table 1. Summary of Multivariate Analysis of Variance (MANOVA) joint with all interactions.
Source Df Pillai Approx F Num Df Den Df Probabilty
Year 1 0.882 44.831 31 186 2.2e-16 **
Genotypes 17 6.5224 4.056 527 3434 2.2e-16 **
Water regime 3 2.6453 45.234 93 564 2.2e-16 **
Year:Block 4 2.2584 7.906 124 756 2.2e-16 **
Year:Genotypes 17 1.6662 0.708 527 3434 1ns
Year:Water regime 3 0.9002 2.600 93 564 7.026e-12 **
Genotypes:Water regime 51 9.6958 1.928 1581 564 2.2e-16 **
Year:Block:Genotypes 68 9.9413 1.500 2108 6699 2.2e-16 **
Year:Genotypes:Water regime 51 1.8329 0.266 1581 6699 1ns
Residuals 216
Table 2. Mean values of vegetative indices and yield of wheat according to the groups of genotypes and water regimes in 2018 and 2019 cropping years.
Table 2. Mean values of vegetative indices and yield of wheat according to the groups of genotypes and water regimes in 2018 and 2019 cropping years.
G1/ WR Variables
NDVI SAVI PRI DVI GRVI GNDVI NDRE TCARI OSAVI TO GH MTG GY
1 22 0.27 0.18 0.1 1.81 3.09 0.5 0.12 0.12 0.22 0.54 10.5 32.97 2337
43 0.45 0.28 0.16 3.17 3.89 0.58 0.25 0.14 0.36 0.4 10.8 35.19 3691
81 0.62 0.38 0.21 5.58 4.7 0.63 0.34 0.18 0.49 0.38 11.3 37.56 5076
100 0.65 0.39 0.22 6.27 5.06 0.65 0.35 0.18 0.51 0.37 11.4 38.74 5604
2 22 0.29 0.19 0.12 1.93 3.15 0.51 0.14 0.12 0.24 0.54 11.2 33.68 1983
43 0.47 0.3 0.17 3.38 3.97 0.58 0.26 0.15 0.38 0.41 10.9 36.2 3913
81 0.63 0.39 0.23 5.91 4.82 0.63 0.35 0.18 0.5 0.37 11.6 36.99 5095
100 0.65 0.41 0.23 6.6 5.2 0.66 0.37 0.18 0.52 0.36 11.8 38.43 5439
3 22 0.24 0.16 0.09 1.67 2.93 0.49 0.11 0.11 0.2 0.55 10.3 32.43 2366
43 0.4 0.24 0.14 2.56 3.51 0.55 0.21 0.13 0.31 0.43 10.3 36.31 4013
81 0.6 0.37 0.21 5.04 4.38 0.61 0.31 0.19 0.47 0.4 10.6 39.97 5653
100 0.64 0.39 0.22 6.04 4.86 0.64 0.33 0.19 0.51 0.39 10.8 40.43 6052
Mean 0.49 0.31 0.17 4.16 4.13 0.59 0.26 0.16 0.39 0.43 10.9 36.57 4269
SE 0.05 0.03 0.01 0.55 0.24 0.02 0.03 0.01 0.04 0.02 0.15 0.76 414.5
CV 3.28 3.1 2.97 4.62 1.98 1.03 3.67 2.02 3.19 1.66 0.46 0.72 3.36
1/G: Groups of genotypes; WR: Water regimes (22, 43, 81, and 100% CET - WR1, WR2, WR3 e WR4, respectively); NDVI: Normalized Difference Vegetation Index; SAVI: Soil Adjusted Vegetation Index; PRI: Photochemical Physiological Reflectance Index; DVI: Difference Vegetation Index; GRVI: Green-Red Vegetation Index; GNDVI: Green Normalized Difference Vegetation Index; NDRE: Normalized Difference Red Edge; TCARI: Chlorophyll Absorption and Reflectance Index; OSAVI: Optimized Soil Adjusted Vegetation Index; TO: TCARI/OSAVI; GH: grain humidity (%); MTG: mass of thousand grains (g); GY; grain yield (Kg.ha-1); DRI: Drought Resistance Index; 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 Sintonia. SE: standard error; CV: coefficient of variation (%).
Table 3. Mean values morphophysiological and drought tolerance index of wheat according to the groups of genotypes and water regimes in 2018 and 2019 cropping years.
Table 3. Mean values morphophysiological and drought tolerance index of wheat according to the groups of genotypes and water regimes in 2018 and 2019 cropping years.
G1/ WR Variables
WUE Cha Chb RM SM A gs iWUE Ci E Fv’/Fm’ ETR Fv/Fm DRI
22 22.7 30.9 7.2 4.9 10.4 8.7 0.09 121.9 188 2.5 0.46 122 0.8
1 43 19.5 35.5 10.5 5 12.2 15.8 0.13 96.8 204 3.2 0.5 144 0.82 1.01
81 14.3 39.4 14.1 5.8 15.7 22.3 0.33 67.8 248 6.3 0.57 148 0.82
100 12.5 41.0 16.0 5.7 15.6 22.7 0.48 47.2 275 8.7 0.59 151 0.83
22 20.3 32.0 8.2 5.4 11.1 9.7 0.09 125.6 193 2.5 0.47 133 0.81
2 43 19.1 35.3 11.4 9.8 12.9 15.1 0.12 108.4 175 3.2 0.51 146 0.83 0.95
81 14.1 40.2 16.1 6.1 13.6 21.6 0.39 55.5 263 7.8 0.6 150 0.83
100 12.2 41.6 16.7 5.6 15.2 22.1 0.53 41.8 273 9.9 0.6 164 0.83
22 22.0 27.4 5.4 4.8 9.9 8.1 0.09 127.1 168 2.6 0.4 99 0.8
3 43 20.3 32.6 8.2 5.1 11.6 16.0 0.1 112.0 144 2.9 0.5 116 0.81 0.93
81 16.0 38.6 13 5.2 15.7 23.1 0.41 56.4 259.9 7.68 0.6 159 0.82
100 13.6 40.2 14.2 5.5 17.1 23.9 0.61 49.3 294.9 11.55 0.6 142 0.83
Mean
SE
CV
17.2 36.1 11.7 5.7 13.4 17.4 0.28 84.1 224.2 5.75 0.5 140 0.82 0.96
1.1 1.34 1.1 0.4 0.7 1.7 0.06 4.87 14.56 0.96 0.02 5.49 0.09 0.01
2.2 1.28 3.3 2.3 1.8 3.4 7.05 1.12 2.25 5.77 1.29 1.36 0.13 0.11
1/G: Group of genotypes; WR: water regimes (22, 43, 81, and 100% CET - WR1, WR2, WR3 e WR4, respectively); WUE: water use efficiency; Cha: chlorophyll a; Chb: chlorophyll b; RM: root mass; SM: shoot mass; A: net assimilation of CO2 (µmol CO2 m-2 s-1); gs: stomatal conductance (mol H2O m-2s-1); Intrinsic water use efficiency (iWUE - A/gs: net assimilation of CO2/stomatal conductance); Ci: internal CO2 concentrations; E: transpiration (mmol H2O m-2 s-1); Fv’/Fm’: Effective Quantum Yield of Photosystem II; (ETR) electron transport rate; Fv/Fm: Maximum Quantum Yield of Photosystem II. 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 Sintonia. SE: standard error; CV: coefficient of variation (%).
Table 4. Relative importance (percentage) of variables for distance between genotypes according to Singh (1981) criterion.
Table 4. Relative importance (percentage) of variables for distance between genotypes according to Singh (1981) criterion.
Variables Percentage (%)
SAVI 33
OSAVI 33
DVI 11
NDRE 11
GNDVI 3
NDGI 2
RVI 2
NDVI 1
GRVI 1
GY 1
Other < 2
Table 5. Regression equation for wheat variables as a function of water regime (x) for groups* of genotypes.
Table 5. Regression equation for wheat variables as a function of water regime (x) for groups* of genotypes.
Variable Group of genotype Equation R squared
Main Coordinate 1 (latent variable) 1 y = 1.6339 + 0,0268 x   0.98
2 y = 1.4824 + 0,0266 x 0.98
3 y = 2.1757 + 0,0326 x 0.98
Grain yield (kg ha-1) 1 y = 5864.37 / ( 1 + exp 0.04 * x 31.3 ) 0.99
2 y = 5364.0 / ( 1 + exp 0.07 * x 29.36 ) 0.99
3 y = 6194.41 / ( 1 + exp 0.05 * x 31.32 ) 0.99
Net CO2 assimilation (A, µmol CO2 m-2 s-1) 1 y = 23.13 / ( 1 + exp 0.062 * x 30.24 ) 0.99
2 y = 23.18 / ( 1 + exp 0.047 * x 29.12 ) 0.99
3 y = 24.20 / ( 1 + exp 0.064 * x 32.53 ) 0.99
Normalized Difference Vegetation Index (NDVI) 1 y = 0.52 + 0.2562 * l o g ( x )   0.99
2 y = 0.46 + 0.2454 * l o g ( x ) 0.99
3 y = 0.62 + 0.2739 * l o g ( x ) 0.99
Water use efficiency (WUE) 1 y = 25.37 0.131 * l o g ( x ) 0.98
2 y = 23.22 0.109 * l o g ( x ) 0.99
3 y = 24.65 0.108 * l o g ( x ) 0.99
Intrinsic water use efficiency (iWUE) 1 y = 135.88 0.846 * l o g ( x ) 0.96
2 y = 133.23 0.833 * l o g ( x ) 0.98
3 y = 126.47 0.824 * l o g ( x ) 0.97
*Group (1: Aliança; BRS 254; BRS 404; CPAC 01019; PF 020062; PF 080492; 2: BR 18_Terena; MGS Brilhante; PF 020037; PF 120337; 3: BRS 264; BRS 394; CPAC 01047; CPAC 07258; CPAC 8318; CPAC 9110; PF 100368; TBIO Sintonia). Water regime (WR1, WR2, WR3 and WR4 - 22, 43, 81, 100% CET replacement, respectively).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated