Screening of Morphophysiological, Anatomical, and Ultrastructural Traits to Improve the Elite Genotype Selection in Sugarcane ( Saccharum ofﬁcinarum L.)

: Sugarcane is a highly productive crop with high water requirements. In addition, its growth and sugar content are limited by the water deﬁcit, a major problem affecting agriculture due to climate change. However, monitoring strategies are being developed worldwide that seek to increase productivity in the same area and with less water consumption. The sugarcane plants are produced sugar table, ethanol, and, from the hydrolysis of biomass, produce second-generation bioethanol, in addition to generating energy in thermoelectric plants. This research described the strategies and mechanisms used by sugarcane to tolerate water deﬁcit. For this, a series of physiological, biochemical, enzymatic, morphological, anatomical, and ultrastructural analyses were developed. We used four commercial varieties of sugarcane, two tolerant and two sensitives, which were tested in all phases of the study. It was concluded that the variety RB92579 is drought tolerant as well as RB867515 and RB72454 genotype is sensitive to drought stress. Therefore, we proposed that variety RB855536 be assigned as an intermediary due to tolerance and sensitivity to water deﬁcit.


Introduction
Food and water are essential elements for human existence and, along with energy, are important for economic growth, poverty reduction, and social development. However, areas for food cultivation and energy-generating species are being depleted, with an increased risk due to global climate change, provoking abiotic stressors [1]. Water stress is very limiting to crop production, including sugarcane culture [2]. Moreover, it has been estimated that soil and water salinity is pervasive throughout the globe, reaching 20% of the total irrigated land area in the world [3].
Sugarcane (Saccharum spp.) is a widely known grass that stands out in world agriculture, mainly due to its participation in the production of table sugar (80% worldwide) and bioethanol (60% worldwide) [4]. In a recent paper, Zan et al. [4] studying 236 global accessions of sugarcane, described 39 accessions from Australia, 13 from Brazil, 81 from France, 30 from the Philippines, two from Cuba, one from Malaysia, 16 from the USA, three from Bangladesh, 14 from Mexico, five from the Republic of Sudan, two from India, and 30 from China. Sugarcane is cultivated in tropical and sub-tropical areas in over 100 countries covering~27.1 million ha with a total harvest of 1.9 billion metric tons/year [5]. However, Cursi et al. [6] described that Brazil is the world's largest producer of sugarcane, followed by India, China, Thailand, Pakistan, and Mexico, reducing the main production countries to six. In the last 40 years, the cultivation of sugarcane has made significant progress, thereby making Brazil a world leader in technologies for the production and processing of table sugar [6,7]. In accordance with Cursi et al. [5] revision, Brazil has 404 sugar, flex-ethanol, historical data and precipitation for 2020 were obtained from the Meteorological Station of Campo Alegre-AL, available in Proclima [15] (more details in Supplementary File S1), while the temperature, relative humidity, and global radiation data were obtained from the Meteorological Station of Arapiraca (~40 km away), available from INMET [16]. In accordance with the Köppen classification [17], the climate is considered to have As. In this region, the De-Martone aridity index [ P T+10 , where P denote annual accumulated rainfall and T annual mean temperature] is 37.4, denoting a winter cloudy and rainy and dry summer.
As shown in Table 1, the tolerant-like (RB92579 and RB867515) and sensitive-like (RB855536 and RB72454) variables to water deficit were planted on the Santa Isabel Farm, in a commercial cultivation area of about 38 m × 30 m where eight sample plots were installed, each presenting 8 m long and 5 m wide. All culture techniques applied to the crop were the standard techniques required for sugarcane grown in the field adapted to Brazilian climatic conditions. Mature (at least 9 months) ratoon sugarcane plants (S. officinarum) i.e., sugarcane from the second year of production (2nd cycle) of the same individual, were randomly (as a bingo ball) selected for this study.
Pearson's correlation coefficient was used to examine the relationships among variables. The similarity between groups of data was calculated using an adjusted r coefficient through Sigmaplot 14.0 (Systat Software, Inc., San Jose, CA, USA). The diagram with true and false colors was constructed in Excel 360 for Windows ® as a scale color.

Leaf Water Potential Measurements (Ψ w )
Nine-month sugarcane plants of each variety were selected for this study. The leaf water potential (Ψ w ) measurements were performed using a Scholander-type pressure chamber (PMS-615D Pressure Chamber, PMS Instrument Company, Albany, OR, USA). The measurements were carried out at predawn (4:00 a.m.) and at midday (12:00 p.m.) using the +3 attached fully expanded leaf from the apex of each variety.

Leaf Gas Exchange Parameters and Chlorophyll a Fluorescence
The leaf gas exchange and chlorophyll a fluorescence was determined on the +3 attached fully expanded leaf from the apex, using a portable open-flow infrared gas analyzer (LI-6400XT; LI-COR Inc., Lincoln, NE, USA) integrated fluorescence chamber heads (LI-157 6400-40; LI-COR Inc.). The net photosynthesis (P N , µmol CO 2 m −2 s −1 ), stomatal conductance (g s , mmol H 2 O m −2 s −1 ), transpiration (E, mmol H 2 O m −2 s −1 ), leaf temperature (T leaf , • C), vapor pressure deficit (VPD, kPa), and internal CO 2 concentration (C i , µmol CO 2 mol −1 ) were obtained in rapid growth phase without flowering S. officinarum plants on July 2014 and December 2014. Two solar time (08:00-10:00 h; early morning; EM and 11:00-13:00 h; midday; MD) were evaluated. All analysis was under a clear sky and leaf irradiance of saturation of 1000 µmol photons m −2 s −1 , fixed CO 2 concentration in 390 µmol mol −1 and airflow of 400 µmol s −1 [18]. The fluorescence analysis was carried out as described by Pompelli et al. [19]. With these parameters the intrinsic water use efficiency (P N /g s ; WUE i ), electron transport rate/photosynthesis ratio (ETR/P N ) were calculated. The leaf hydraulic conductance (KL, mmol H 2 O m −2 s −1 MPa −1 ) was calculated using Ψ wpredawn to approximate soil water potential, and g s and VPD were measured at the same time as Ψ wmidday according to Dos Santos et al. [20]. Tall size, but with low growth rate. Higher responsive to irrigation and very efficient in water use, with suitable recovery after drought periods. High agroindustrial productivity, reaching more than 35% of the sugarcane area in the northeast between 2010 and 2015. It improved sugarcane productivity in the northeast region, with yields increasing up to 60% above previously cultivated varieties. Low growth rate, medium-late maturation, long period of industrialization. High sucrose level and medium fiber content. Very resistance to water deficit.

Sap Flow Measurements
The driving force, which drives sap ascent in the xylem (SAP flow velocity) was measured with heat dissipation sensors (Heat Pulse Velocity, Dual Method Approach, Implexx Sap Flow Sensor, Moorabbin, Australia) connected to an internal datalogger (Implexx Sap Flow Datalogging, Moorabbin, Australia). The measurements were carried out for 24 h in open field cultivated sugarcane plants (from 5:30 p.m. to 4:30 p.m.; +1) in on 20 July 2020, and 17 December 2020, to characterize the pattern of sap flow over 24 h in two extreme rainfall regimes.

Biochemical and Enzymes Activity Measurements
For biochemical and enzyme measurements, the +3 leaf fragments were collected around 10:00 a.m., immediately immersed in liquid nitrogen, and stored in the freezer (−20 • C) until used. The extraction of all biochemical compounds and the enzyme activity from the tissues was thoroughly ground with a cold mortar and pestle in an ice bath, plus 0.1 g of the sand river to improve the extraction until no fibrous residue could be seen. Chlorophyll a + b and total carotenoids were extracted with 80% (v/v) aqueous acetone and quantified spectrophotometrically [24]. Quantification of total free amino acids and protein contents of crude extracts were determined using the Moore and Stein [25] and Bradford [26] methods, respectively. To measure the superoxide dismutase (SOD; EC 1.1.5.1.1), catalase (CAT; EC 1.11.1.6), ascorbate peroxidase (APX; EC 1.11.1.11) and malondialdehyde (MDA) analysis was carried out in accordance with Pompelli et al. [27]. The proline concentration was determined via the method of Bates et al. [28] with fine adjustment described in Mendes et al. [29].
Plant height was measured from the soil surface to the highest point of insertion of the leaves using a metric tape attached to a broomstick. The stalk diameter was measured at 50 cm above the plant base through a pachymeter (Electronic Digital Vernier Caliper, Stainless Steel Caliper 150 mm/0-6 inch). The leaf area was measured using the allometric equations proposed by triticale as a gramineous plant-like [30].
For optical anatomical measurements, +3 leaf fragments were collected at around 10:00 a.m. and fixed in FAA 50 for 48 h and then kept in 70% (v/v) ethanol [31] until use. All steps from FAA fixed to optical photomicrograph as described in detail in Pompelli et al. [32]. For scanning electron microscopy (SEM) measurements, the +3 leaf fragments were collected around 10:00 a.m. and immediately fixed in Karnovsky solution [33], prepared in 0.1 M cacodylate buffer (sodium cacodylate trihydrate, Sigma Aldrich, St. Louis, MO, USA), pH 7.4 and 2.5% glutaraldehyde (part number G5882, Sigma Aldrich) for 60 h at 4 • C. Thereupon the leaf fragments were dehydrated in an ethanol series (10-100%), subjected to a critical point with carbon dioxide (CPD 030, Critical Point Dryer, Bal-tec, Balzers, Liechtenstein, German), mounted on metallic supports (stubs) and metallized with gold, for analysis in the scanning electron microscope (Zeiss Leo 1530, ETH Zurich, Zurich, Switzerland). Stomatal density (SD) was measured by using Salisbury equations [34] with modifications proposed by Pompelli et al. [19]. However, stomatal area and stomatal pore aperture measurements were performed using SEM photomicrograph. For stomatal, counting was taken into an SEM photomicrograph with at least 20 replicates per treatment and leaf surface (adaxial and abaxial), using Image-Pro Plus ® software 4.0 (Media Cybernetics, Silver Spring, MD, USA). The stomata per micrograph were converted to stomata mm −2 . To measure the stomata aperture, we measured the ostiole in equatorial aperture using Image-Pro Plus software.

Plasticity Index
Differences in the plasticity index associated with physiological, biochemical, enzymatical, anatomical, and ultrastructure ranging from 0 (no plasticity) to 1 (maximal plasticity), were calculated as described in Dos Santos et al. [35].

Statistical Procedures
The experiments were conducted in a completely randomized block design with four sugarcane varieties (RB92579, RB867515, RB855536, and RB72454), two rainfall regimes (2020, July-rainy season and 2020, December-dry season), and two solar times (08:00-10:00 h and 11:00-13:00 h), with 10 replicates. First, all the data were evaluated for normality using the Shapiro-Wilk test, and homogeneity of variance was tested with a Brown-Forsyth test. All data were analyzed by two-way ANOVA, and means were compared using an SNK test (p < 0.05) by Sigmaplot 14.0.

Climatological Data
According to Figure 1, the accumulation rainfall in 2020 (Jan to December) was 1257 mm. The rainfall tendency was followed by the temperature variation, where the rainy months recorded the lowest temperatures, while the driest months recorded the highest temperature. , maximum (red), and mean (green) temperature (right axes), and relative humidity (left axes). (B) Accumulated rainfall (left axes) and global radiation (right axes). All data were automatically registered by a climatologic tower daily between 1 June 2020 and 31 December 2020. Source INMET [36].

Gas Exchange Parameters
The gas exchange parameters perfectly reflect the climatological data, where PN, gs, E, WUEi, show the highest values in the rainy season. The PN values registered in the early morning (08:00-10:00 h) ranged from 1.17-to 4.75-fold higher in the rainy season compared to dry season plants. In the MD (11:00-13:00 h) the range of PN was greater than EM, ranging from 0.41-to 5.95-fold higher in rainy season than compared to dry season plants. Even though the fold measured in MD (mean 2.83-fold) is higher than EM (mean 2.68-fold) this range is not significant (p-value = 0.609).
In rainy season plants, the lesser PN was verified in RB72454 (12.38 μmol CO2 m −2 s −1 ) and higher in RB855536 (25.64 μmol CO2 m −2 s −1 ). In other ways, in the dry season plants, the lesser PN was verified in RB867515 (2.37 μmol CO2 m −2 s −1 ) and higher was verified in RB72454 (7.44 μmol CO2 m −2 s −1 ). This pattern is remarkably interesting because the same variety (RB72454) showed a lesser value in the rainy season and the highest value in the dry season Figure 2). In this sense, PN values registered in variety RB72454 show the lesser

Gas Exchange Parameters
The gas exchange parameters perfectly reflect the climatological data, where P N , g s , E, WUE i , show the highest values in the rainy season. The P N values registered in the early morning (08:00-10:00 h) ranged from 1.17-to 4.75-fold higher in the rainy season compared to dry season plants. In the MD (11:00-13:00 h) the range of P N was greater than EM, ranging from 0.41-to 5.95-fold higher in rainy season than compared to dry season plants. Even though the fold measured in MD (mean 2.83-fold) is higher than EM (mean 2.68-fold) this range is not significant (p -value = 0.609).
In rainy season plants, the lesser P N was verified in RB72454 (12.38 µmol CO 2 m −2 s −1 ) and higher in RB855536 (25.64 µmol CO 2 m −2 s −1 ). In other ways, in the dry season plants, the lesser P N was verified in RB867515 (2.37 µmol CO 2 m −2 s −1 ) and higher was verified in RB72454 (7.44 µmol CO 2 m −2 s −1 ). This pattern is remarkably interesting because the same variety (RB72454) showed a lesser value in the rainy season and the highest value in the dry season Figure 2). In this sense, P N values registered in variety RB72454 show the lesser fold (1.72-fold higher) from 15.83 µmol CO 2 m −2 s −1 ) in the rainy season to 9.21 µmol CO 2 m −2 s −1 in the dry season, both in EM ( Figure 2). In another way, the dataset registered in the EM, variety RB92579 shows the highest fold change when comparing the rainy (19.51 ± 0.53 µmol CO 2 m −2 s −1 ) to dry (5.63 ± 0.19 µmol CO 2 m −2 s −1 ) or 3.46-fold higher in the rainy season when compared to dry season plants. This last variety also showed the highest amplitude in the P N values (6.32-fold higher), registered in rainy (20.74 ± 0.35 µmol CO 2 m −2 s −1 ) and dry season (3.28 ± 0.18 µmol CO 2 m −2 s −1 ) on the MD dataset. In this scenario, we highlighted in the first plan one of the most tolerant varieties (RB92579) and in the second plan to more sensitive (RB72454) variety to water stress. When the comparison was taken P N values between EM to MD values in the same variety and the same season, we showed that the RB867515 varieties show a ratio of 1.04-fold higher in the EM when compared to MD. For clarity, only RB867515 variety shows higher P N values in the EM, while the other three variety RB855536, RB92579, and RB72454, showed 0.73-, 0.94-, and 0.99-fold higher in the EM values when compared to MD values. The less than 1-fold values denote that EM P N values were lesser than those registered in MD (Figure 2A,B).
The same tendency to show largest amplitudes in the MD dataset to g s and E. In this scenario we highlighted the RB867515 variety that showed 6.2-and 3.9-fold higher, respectively to g s and E when compared to the rainy than dry season dataset. In another way the variety RB72454 showed a lesser amplitude between values registered in the rainy and dry seasons, when rainy season g s and E dataset were 1.6-and 0.8-fold higher than dry season ones. It is worth noting that this 0.8-fold higher or 1.25 ± 0.08 mmol H 2 O m −2 s −1 (rainy season) is not statistically different to 1.28 ± 0.06 mmol H 2 O m −2 s −1 (dry season) (p-value = 0.33). In the EM were RB92579 and RB855536 varieties that showed, respectively the highest (3.1-fold) and lesser (1.1-fold) amplitude to g s , values registered, respectively in the rainy and dry seasons. The RB92579 also registered a moderated amplitude on E (2.1-fold), when comparing the rainy (1.83 ± 0.07 mmol H 2 O m −2 s −1 ) to dry (0.87 ± 0.03 mmol H 2 O m −2 s −1 ) season plants. In this variable, the RB72454 shows a lesser amplitude (0.76-fold), in the rainy season (1.25 ± 0.08 mmol H 2 O m −2 s −1 ) than those registered in the dry season (1.64 ± 0.02 mmol H 2 O m −2 s −1 ).
The values for WUE i , measured in EM were statistically different from those registered in the dry season, while in the MD dataset only RB92579 and RB855536 show higher values in the rainy season compared to from those registered in the dry season. When we compared the varieties using only rainy season plants on an EM, we verified that the highest and lowest value were registered in the RB867515 (159.73 ± 2.4 µmol CO 2 mol −1 H 2 O) and RB855536 (90.30 ± 6.6 mmol CO 2 mol −1 H 2 O). In another way, when we compared the dry season, it was verified that except for the RB855536 variety all other means that do not were statistically different. All sugarcane plants showed a high value of WUE i , measured in the rainy season when compared to the dry season, except RB867515 variety where the dry season plants showed an increase of 50.3% in the dry season compared to rainy season plants. Figure 3 showed that all the varieties measured at MD reinforce that higher P N provokes lower C i :C a ratio, with depletion of C i in the substomatal space except in RB867515, RB855536. This hypothesis is corroborated by a higher C i :C a ratio, which was higher than 0.5, due to lower P N in these varieties, at least on MD. Another fact that reinforces this hypothesis is the lower C i :C a ratio in the rainy season, in the same moments when the P N is higher. This fact also is improved by the strong negative correlation between P N and C i :C a ratio (r = −0.731; p = 5.71 × 10 −28 ). The values of C i :C a ratio in the dry season is partially due to lower g s , which shows a moderate and negative correlation between C i :C a ratio and g s (r = −0.485; p = 8.02 × 10 −11 ). A similar pattern was verified with T leaf , where higher values are provoked to lower E (r = −0.242; p = 2.09 × 10 −3 ). Another feature that confirms our hypothesis that higher T leaf is provokes by lower g s , is verified with the highest correlationship between T leaf and C i :C a ratio (r = 0.509; p = 6.43 × 10 −12 ). Therefore, we can infer that there are three forces governing gas exchange, (i) the increase in g s and E, (ii) the higher VDP, (iii) the low rates of hydraulic conductivity of the stems, which, in turn, is strongly decreased by the very negative water potencies (see more details in Discussion). The strong influence of VPD on lower g s can be confirmed by both strong negative correlation between these value (r = −0.638; p = 1.15 × 10 −19 ) ( Figure 4). The values for WUEi, measured in EM were statistically different from t registered in the dry season, while in the MD dataset only RB92579 and RB855536 s higher values in the rainy season compared to from those registered in the dry sea When we compared the varieties using only rainy season plants on an EM, we ver that the highest and lowest value were registered in the RB867515 (159.73 ± 2.4 μmol mol −1 H2O) and RB855536 (90.30 ± 6.6 mmol CO2 mol −1 H2O). In another way, when compared the dry season, it was verified that except for the RB855536 variety all o     Figure 5 showed us that between 5:30 p.m. and 4:00 a.m. the sap not significantly different with others. After predawn, the varieties sta patterns. While RB86715 rose moderately, accelerating later and from sharp drop, returning the initial values, the RB92579 showed a continuous increase until 4:30 p.m. with no significant reduction in s varieties are considered to be more drought-tolerant than others (Figu RB855536 and RB72454 showed an increase in the sap flow from 5 showed a more abrupt rise than RB72454. However, while RB855536 flow around 11:30 a.m. and showed an abruptly decrease to null sap f the RB72454 remain increase the sap flow to 12:00 p.m. as it slowly dro In the dry season, all plants do not show a significant increase in ( Figure 5B). The rate of sap flow in the rainy season was 2-to 5-fold h in the dry season, independent of the variable.  Figure 5 showed us that between 5:30 p.m. and 4:00 a.m. the sap flow in all variety is not significantly different with others. After predawn, the varieties start to show different patterns. While RB86715 rose moderately, accelerating later and from 2:30 p.m. showed a sharp drop, returning the initial values, the RB92579 showed a more delayed and continuous increase until 4:30 p.m. with no significant reduction in sap flow. These two varieties are considered to be more drought-tolerant than others ( Figure 5A). The varieties RB855536 and RB72454 showed an increase in the sap flow from 5 a.m., but RB855536 showed a more abrupt rise than RB72454. However, while RB855536 start a drop in sap flow around 11:30 a.m. and showed an abruptly decrease to null sap flow from 1:00 p.m., the RB72454 remain increase the sap flow to 12:00 p.m. as it slowly dropped to initial value. In the dry season, all plants do not show a significant increase in its sap flow values ( Figure 5B). The rate of sap flow in the rainy season was 2-to 5-fold higher than sap flow in the dry season, independent of the variable.

Water Potential (Ψw)
The Ψw closes the water flux in the soil-plant-atmosphere system. In the rain we described a slight difference to Ψwpredown. However, in Ψwmidday the varieties R (Ψwmidday −1.39 MPa) and RB72454 (Ψwmidday −1.38 MPa) varieties benefited du negative shown Ψwmidday in comparison to RB92579 and RB867515, values that we higher than those shown to RB92579 and RB867515. The Ψwpredown was 4-to 5-fol than those values registered in the midday, independently of the season (rainy The rainy season Ψwpredown ranged from −0.

Water Potential (Ψ w )
The Ψ w closes the water flux in the soil-plant-atmosphere system. In the rainy season, we described a slight difference to Ψ wpredown . However, in Ψ wmidday the varieties RB855536 (Ψ wmidday −1.39 MPa) and RB72454 (Ψ wmidday −1.38 MPa) varieties benefited due to less negative shown Ψ wmidday in comparison to RB92579 and RB867515, values that were 10.8% higher than those shown to RB92579 and RB867515. The Ψ wpredown was 4-to 5-fold higher than those values registered in the midday, independently of the season (rainy or dry). The rainy season Ψ wpredown ranged from −0.

Chlorophyll a Fluorescence Measurements
Without exception, all rainy season sugarcane plants showed higher v ΦPSII, and ETR, when compared to those registered in the dry season, indepe sample time (Figure 7). The Fv/Fm measured EM in rainy season plants is n Contrasting, the dry season plants showed a 2-level pattern: one joining the RB867515 and the other joining the RB855536 and RB72454, with the last showing an increase of 20.9% of the Fv/Fm measured in the first varieties. measurements in the rainy season tend to be not physiologically differen values ranged from 0.72 (RB72454) to 0.78 (RB92579). In another way, t sugarcane plants, the Fv/Fm values registered in MD plants were higher in sharp decrease in other varieties RB867515 (4.9%), RB855536 (9.4%), and RB The pattern verified in ΦPSII and ETR tended to be same as the pattern desc measurements. So, we presented more integration of these parameters in t However, in suitable agreement, the reductions in AN may be associated w availability of NADPH and ATP, since both ΦPSII and ETR decreased due to (Figure 7).

Chlorophyll a Fluorescence Measurements
Without exception, all rainy season sugarcane plants showed higher values of F v /F m , Φ PSII , and ETR, when compared to those registered in the dry season, independently of the sample time ( Figure 7). The F v /F m measured EM in rainy season plants is nonsignificant. Contrasting, the dry season plants showed a 2-level pattern: one joining the RB92579 and RB867515 and the other joining the RB855536 and RB72454, with the last two varieties showing an increase of 20.9% of the F v /F m measured in the first varieties. The MD F v /F m measurements in the rainy season tend to be not physiologically different because the values ranged from 0.72 (RB72454) to 0.78 (RB92579). In another way, the dry season sugarcane plants, the F v /F m values registered in MD plants were higher in RB92579 with sharp decrease in other varieties RB867515 (4.9%), RB855536 (9.4%), and RB72454 (10.4%). The pattern verified in Φ PSII and ETR tended to be same as the pattern described to F v /F m measurements. So, we presented more integration of these parameters in the discussion. However, in suitable agreement, the reductions in A N may be associated with the lower availability of NADPH and ATP, since both Φ PSII and ETR decreased due to drought stress ( Figure 7).

Figure 7.
Variable-to-initial Chl fluorescence (Fv/Fm; (A,B)), the actual PSII efficiency during the exposure to sunlight (ΦPSII; (C,D)), and electron transport rate (ETR; (E,F)) measured in four varieties (RB92579, RB867515, RB855536, and RB72454) of sugarcane. The evaluations were made in dry season (gray bars) and rainy season (green bars) at 8:00 to 10:00 h (A,C,E) or at 12:00 to 13:00 h (B,D,F). Different lowercase letters denote significance within season for each variety, and different capital letters denote significance within varieties for the same season. An asterisk (*) denotes significance within 8:00 to 10:00 h and 12:00 to 13:00 h measurements in the same treatment.

Chlorophyll and Carotenoids
Independently of variety, the chlorophyll (Chl) a showed a decrease of ~34.1% on dry season plants when compared to rainy season ones. The highest and lesser drop was shown in RB72454 (41.3%) and RB867515 (28.4%). Both within rainy and dry season plants, variety RB92579 stood out for showing greater concentrations than the other varieties that did not show significant differences between them ( Table 2).

Figure 7.
Variable-to-initial Chl fluorescence (F v /F m ; (A,B)), the actual PSII efficiency during the exposure to sunlight (Φ PSII ; (C,D)), and electron transport rate (ETR; (E,F)) measured in four varieties (RB92579, RB867515, RB855536, and RB72454) of sugarcane. The evaluations were made in dry season (gray bars) and rainy season (green bars) at 8:00 to 10:00 h (A,C,E) or at 12:00 to 13:00 h (B,D,F). Different lowercase letters denote significance within season for each variety, and different capital letters denote significance within varieties for the same season. An asterisk (*) denotes significance within 8:00 to 10:00 h and 12:00 to 13:00 h measurements in the same treatment.

Chlorophyll and Carotenoids
Independently of variety, the chlorophyll (Chl) a showed a decrease of~34.1% on dry season plants when compared to rainy season ones. The highest and lesser drop was shown in RB72454 (41.3%) and RB867515 (28.4%). Both within rainy and dry season plants, variety RB92579 stood out for showing greater concentrations than the other varieties that did not show significant differences between them ( Table 2).
A different situation was shown in the Chl b concentration, where the two droughttolerant-like varieties had an average of 0.40 g kg −1 DW, while the variety RB855536 had 0.42 ± 0.03 g kg −1 DW and variety RB72454 had 0.36 ± 0.02 g kg −1 DW, only when compared between rainy season plants. In this sense, there was a strong drop in Chl b in the tolerant varieties (RB92579 and RB867515) of 49.2% and 28.4%, while in RB855536 there was a nonsignificant decline (5.5%) and an increase of~14% in the variety RB72454 (Table 2). Table 2. Biochemical and antioxidative features measured of four varieties (RB92579, RB86751, RB855536, and RB72454) of sugarcane. The evaluations were made in rainy season (2020, July) and dry season (2020, December) at 8:00 to 10:00 h. Different lowercase letters denote significance within season for each variety, and different capital letters denote significance within varieties for the same season. All value denotes median ± SE (n = 10). With differences in Chl a and Chl b patterns, the Chl total was modulated more by reductions in Chl a than in Chl b. However, the rainy season plants showed a level of Chl total of 1.6-, 1.4-, 1.3-, and 1.4-fold higher in favor of the rainy season plants, respectively to RB92579, RB867515, RB855536, and RB72454 (Table 2).

Characteristics
While Chl levels, in general, decreased in dry season plants, total carotenoids were increased in dry season ones by 37% (RB92579) and 56.2% (RB72454), while the varieties RB867515 and RB855536 presented a median and very low, but not significant fall of 15.5% and 2.4% in dry season plants when compared to rainy season plants ( Table 2).
With the significant decreases in total chlorophylls and a slight increase in carotenoid levels in dry season plants, we showed an increase in the total carotenoids/chlorophylls (Car/Chl) ratio. In this sense, it was shown that the varieties RB92579 and RB72454 presented a strong increase in the Car/Chl ratio in the order of 120.9% and 120.8%, while in the varieties RB867515 and RB855536 the increase was slight in the variety RB867515 (17%) and moderate in the variety RB855536 (32%).

Free Amino Acids and Total Proteins
The concentration of free amino acids [aa] showed a quite different pattern between tolerant and sensitive sugarcane plants to water deficit. While varieties RB92579 and RB867515 showed an increase in concentrations in dry season plants (21.6% and 9.9%, respectively), the sensitive plants showed a moderate (RB855536; 29.2%) and negligible (RB72454; 1.6%) decrease in free [aa]. However, when comparing the varieties within the same season, in the rainy season sugarcane plants, the variety RB855536 showed an [aa] of 16.73 ± 0.25 mmol kg −1 DW, which was 1.1%, 5.6%, and 12.6% higher than those presented in the varieties RB72454, RB867515, and RB92579, respectively. A completely different pattern was verified in the dry season sugarcane plants, where the variety RB92579 presented the highest [aa] than the other varieties, values that were 2.4%, 8.4%, and 33.4% higher than the other varieties (RB867515, RB72454, and RB855536, respectively). However, the difference between [aa] in RB92579 and RB867515 or RB72454 is not significant.
Regardless of variety, dry season sugarcane plants showed a moderately lower protein concentration than rainy season sugarcane plants. The varieties RB867515, RB92579, RB72454 and RB855536 presented a decrease of 13.6%, 19.9%, 20.6%, and 23.1% in dry season, respectively. Moreover, when the comparison is performed by the variety within rainy season plants, the variety RB92579 presented a protein concentration of 112.36 ± 3.93 g kg −1 DW, a value that was 13.8%, 18.4%, and 26.2% higher than the concentration of total proteins in the varieties RB855536, RB867515, and RB72454, respectively. However, when comparing the dry season sugarcane plants, the RB92579 variety showed the highest protein concentration (89.97 ± 1.16 g kg −1 DW), a value that was 11.9%, 17.2%, and 26.8% higher than the total protein concentration of the varieties RB867515, RB855536 and RB72454, respectively ( Table 2).

Antioxidant System
The antioxidant system was strongly modulated by water deficit, especially in the varieties RB92579 and RB855536. SOD activity was respectively increased by 52.4% and 74% in these two varieties. However, the other two varieties also showed a very expressive increase in SOD activity, with 60% (RB867515) and 107.8% (RB72454). Moreover, it was the RB855536 variety that showed the highest SOD activity in both rainy and dry season plants. When comparing SOD activity within rainy season plants, the variety RB855536 showed an activity of 23.7%, 28.5%, and 39.1% higher when compared to RB867515, RB92579, and RB72454, respectively. When SOD activity was compared only in dry season plants, the RB855536 was increased in the greatest proportion, being 16.5%, 34.5%, and 46.8% higher than the SOD activity of the varieties RB72454, RB867515, RB92579, respectively ( Table 2).
Like SOD activity, the dry season plants showed an increase in CAT activity of 15.8%, 34.3%, 70.8%, and 85.3%, respectively in RB72454, RB867515, RB855536, and RB92579 varieties. When comparing the CAT activity within the rainy season plants, it was verified that the variety RB92579 presented a 21.2% lower CAT activity than the other varieties that did not differ statistically between them. However, in dry season plants, the variety RB855536 showed the highest activity among all other varieties (36.26 ± 0.53), or 13.2%, 27.5%, and 48.3% higher than the activity shown, respectively, in varieties RB92579, RB867515, RB72454 ( Table 2).
APX activity was distinctly modulated in the varieties. In varieties RB92579 and RB855536 the APX activity was mild, but statistically increased by 5.1% and 11.4%, while in the varieties RB867515 and RB72454, an APX activity in dry season plants was 7.3% and 28.8% lower than those activities showed in rainy season plants. When comparing the varieties in dry season plants, the variety RB855536 presented a 10.8%, 20.1%, and 62.6% higher APX activity than the varieties RB92579, RB867515, and RB72454, respectively ( Table 2).
As a result of water stress, the dry season sugarcane shows a higher concentration proline when compared to the rainy season plants, with expressive increases in the RB92579 (53.1%) and RB855536 (71.8%). In the varieties RB867515 and RB72454, the proline concentration in the dry season plants was respectively increased by 16.1% and 15.9%. Within the rainy season plants variety RB72454 had the highest proline concentration (51.10 ± 0.28 mmol kg −1 DW), which was 1.1-, 1.3-, and 1.3-fold higher proline than accumulated in RB867515, RB92579, and RB855536, respectively. While the dry season plants, RB855536 plants had a higher proline concentration (66.81 ± 0.32 mmol kg −1 DW), which was 8.8%, 12.8%, and 21.4% higher than those concentrations in RB92579, RB72454, and RB867515, respectively ( Table 2).
The MDA was strongly (91.3% and 79.2%) increased in RB72454 and RB867515 and slightly increased (10.5%) in the RB92579 variety, while in the RB855536 variety the value of rainy and dry season plants was not statistically different, and the dry season plants accumulated 0.61% less MDA than rainy season plants. In relative values, the MDA accumulation in the rainy season was 24.4% (RB92579), 34.9% (RB867515), and 39.8% (RB855536) lower than that accumulated in RB72454 (25.23 ± 0.61 mmol kg −1 DW). However, when only the dry season plants are compared, the variety RB72454 presented a concentration of 48.27 ± 0.15 mmol kg −1 DW or of 44 g kg −1 DW %, 115.3%, and 169.2% more MDA, respectively in relation to the varieties RB867515, RB92579, and RB855536 (Table 2).
It can be considered that the RB855536 genotype, even with a lower F v /F m ratio, did not have their photosynthetic and hydraulic apparatus negatively affected, since the predawn potential was 62% higher in the dry season in comparison with the rainy season, which in part can be explained by the anticipation of the decrease in sap flow ( Figure 5).
Thus, during the hottest moments of the day, these plants reduce their sap flow, causing a stomatal closing that is more affected in dry season plants ( Figure 8). The mechanisms that command this modulation could have the higher concentration of Car/Chl ratio, or the higher Φ PSII , since higher the Φ PSII reduce the need to dissipate heat, promoting a security thermal dissipation. These hypothesis gives more strength by positive correlation with F v /F m ratio and Chl a (r = 0.522; p = 7.03 × 10 −7 ) or with total chlorophyll (r = 0.428; p = 7.47 × 10 −5 ). The antioxidative enzymes also corroborates this hypothesis, where F v /F m ratio and SOD (r = −0.758; p = 3.71 × 10 −16 ), CAT (r = −0.546; p = 1.37 × 10 −7 ) and proline (r = −0.694; p = 9.53 × 10 −13 ) had strong negative correlations.

Plant Anatomy and Ultrastructure
The total leaf thickness measured in dry season plants (Table 3) showe (RB855536), 22% (RB72454), 24.1% (RB867515), and 25.8% (RB92579) thinner leaf t than the rainy season plants. Within the rainy season plants the variety RB72454 thickness being 14.7%, 27.8%, and 67.9% thicker than RB855536, RB867515, and R respectively. Within the dry season plants, the variety RB72454 presented a 7.1% and 69.5% thicker than RB855536, RB867515, and RB92579, respectively (Table 3) Morphologically, the transversally cut leaf ( Figure 9) demonstrates that thickness is quite irregular, and, for this reason, we had to increase the sampling thickness measurement to 30 true repetitions instead of 10 as shown in Table 3. Re of that, were performed 20 measurements (false repetition) along with the leaf f to give more fidelity to the data. The leaf morphology has a characteristic patte species of the Poaceae family, without a distinction between palisade and parenchyma. It features very prominent bundle sheath cells arranged with 2where the second part of photosynthesis takes place. Mesophyll cells are rarely bundle sheath cells. In the epidermis, the presence of bulliform cells is quite c which are large and favor thermal dissipation (more information, see Discussion
Morphologically, the transversally cut leaf (Figure 9) demonstrates that the leaf thickness is quite irregular, and, for this reason, we had to increase the sampling for leaf thickness measurement to 30 true repetitions instead of 10 as shown in Table 3. Regardless of that, were performed 20 measurements (false repetition) along with the leaf fragment to give more fidelity to the data. The leaf morphology has a characteristic pattern of the species of the Poaceae family, without a distinction between palisade and spongy parenchyma. It features very prominent bundle sheath cells arranged with 2-3 layers, where the second part of photosynthesis takes place. Mesophyll cells are rarely far from bundle sheath cells. In the epidermis, the presence of bulliform cells is quite common, which are large and favor thermal dissipation (more information, see Discussion). Table 3. Morphological and anatomical features of four varieties (RB92579, RB86751, RB855536, and RB72454) of sugarcane. The evaluations were made in rainy season (2020, July) and dry season (2020 December) at 8:00 to 10:00 h. Different lowercase letters denote significance within season for each variety, and different capital letters denote significance within varieties for the same season. All value denotes median ± SE (n = 10).   Table 3.  Table 3.

Characteristics
Except for the variety RB72454, which presented an adaxial stomatic density (SD adaxial ) of 7.85 ± 0.17 stomata mm −2 (rainy) and 6.62 ± 0.12 stomata mm −2 (dry) with a reduction of 15.7% in dry season than compared to rainy season, all other varieties presented a higher SD adaxial in dry season plants in this order: RB92579 (11.1%), RB867515 (17.1%), and RB855536 (25.6%). When comparing only between rainy season plants, the variety RB72454 showed an SD adaxial of 7.85 ± 0.17 stomata mm −2 , a value that was 14.6%, 15.4%, and 20.7% higher than RB855536, RB867515, and RB92579 varieties, respectively. In another way in the dry season plants, the varieties RB867515 and RB855536 showed an average SD adaxial of 8.3 stomata mm −2 , the varieties RB92579 and RB72454 showed an average SD adaxial of 6.9 stomata mm −2 . A similar pattern was verified in abaxial stomatal density (SD abaxial ), where the variety RB72454 showed no significant difference (reduction of 2.4%) between the rainy and dry season plants. In another way, the varieties RB855536, RB92579, and RB867515 had a moderate reduction in SD abaxial in the order of 10.6%, 11.5%, and 14.1%. When the comparison is made only between the rainy season plants, the RB855536 variety, showed an SD abaxial of 16.50 ± 0.22 stomata mm −2 , a value that was 7.2%, 8.7%, and 14.7% higher than the varieties RB72454, RB867515, and RB92579. However, when the analyses are carried out between the dry season plants, there was a strong significant variation between the means, establishing the following sequence: RB855536 > RB867515 > RB92579 > RB72454.
The ostiole opening (OO) showed a similar pattern on both sides of the sugarcane epidermis. Thus, in the tolerant-like varieties, (RB92579 and RB867515) there was a decrease in 34.8% and 34.4% to OO adaxial and 21.2% and 9.1% to OO abaxial in dry season plants in comparison with rainy season plants. On the other hand, in sensitive-like varieties (RB855536 and RB72454), dry season plants showed an increase in ostiole opening in the order of 35% and 7.8% (OO adaxial ) and 8.4% and 62.2% (OO abaxial ). However, when comparing the OO adaxial only among the rainy season plants, the varieties RB867515 (5.36 ± 0.10 µm) and RB72454 (5.48 ± 0.18 µm) showed the highest OO adaxial . In another way the varieties RB92579 and RB855536 presented an OO adaxial 10.1% and 39.7% higher in relation to RB72454. When the comparison is only between the dry season plants, the RB72454 OO adaxial (5.91 ± 0.07 µm) was 11.6%, 68%, and 82.1%, respectively for RB855536, RB867515, RB92579. To comparison in the OO abaxial of rainy season plants, the variety RB855536 presented an OO abaxial of 5.93 ± 0.19 µm, a value that was 56.4%, 83.7%, and 149.8% higher than that registered in the varieties RB92579, RB867515, and RB72454, respectively. However, when the comparison is made only between the dry season plants, the variety RB855536 presented an OO abaxial of 6.43 ± 0.05 mm, a value that was 67%, 115.3%, and 119% higher than that registered in the varieties RB72454, RB92579, and RB867515, respectively.
Ultrastructural, the epidermis is quite rough with an expressive deposition of epicuticular wax, which is more prominent in the abaxial epidermis. As seen in Figures 10 and 11, the epicuticular wax partially and sometimes even completely covers the stomata. Especially on the abaxial face, the leaf showed a strong increase in fiber concentration.  In (e,l,m) the arrow points to the epidermal hair. s. denotes stomata. For more information, see Table 3. Figure 11. SEM micrograph of abaxial epidermis in four varieties (RB92579, (a-d); RB867515, (e-RB855536, (i-k); and RB72454, (i-n)) of sugarcane. The leaf fragments were collected in rainy sea (a,b,e,f,i,j,l,m) and dry season (c,d,g,h,k,n) at 8:00 to 10:00 h. Scale = 10 μm. In (g,i,k), the arr points to the epidermal hair. s. denote stomata. For more information, see Table 1. Figure 11. SEM micrograph of abaxial epidermis in four varieties (RB92579, (a-d); RB867515, (e-h); RB855536, (i-k); and RB72454, (l-n)) of sugarcane. The leaf fragments were collected in rainy season (a,b,e,f,i,j,l,m) and dry season (c,d,g,h,k,n) at 8:00 to 10:00 h. Scale = 10 µm. In (g,i,k), the arrow points to the epidermal hair. s. denote stomata. For more information, see Table 1.

Physiological and Biochemical Features
Abiotic stress, such as drought stress, high salt, extreme temperature, and heavy metal toxicity, frequently leads to oxidative stress. In the present study, physiological, biochemical, enzymatical, leaf anatomical, and ultrastructural features were analyzed with sugarcane varieties to understand the drought-tolerance mechanism responses when they were exposed to the rainy and dry seasons. It is known that gas exchange variables and plant ecophysiology have a suitable correlation with each other and may help to identify genotypes with drought-tolerance characteristics [18,20]. Leaf Ψ wpredown is directly related to soil water potential, which in times of drought leads to a reduction in plant water potential, which can indirectly lead to a decrease in stomatal conductance [27]. When the objective is to compare the capacity of a plant to recover after a water stress event, the measurement of the Ψ wpredown is the most indicated; moreover, when the purpose is to measure the degree of impact of the abiotic stressor, it should consider the Ψ wmidday . The fact that all varieties recovered their Ψ wpredown suggests that all varieties could be considered slight or mildly drought tolerant. However, the fact that the variety RB867515, considered by many authors as resistant [20,37,38], presented a 2.61-fold higher Ψ wpredown in the dry season plants when compared to the rainy season is very interesting. To make a comparison, the variety RB855536, which, unlike the previous one, is considered as sensitive, the potential was modulated by 0.64-fold, i.e., Ψ wpredown of dry plants was higher than that of wet plants.
For photosynthesis to start, it is necessary to have light that is captured by leaf chlorophyll and carotenoid molecules [39]. However, in situations of abiotic stress, two ways can happen to protect the plant from the excess of captured photons: (1) adjustments in light capture, use, and dissipation are required to provide photoprotection to the photosynthetic apparatus, particularly under full sunlight, a process called photobleaching; (2) increase in the concentration of carotenoids, which act more intensely in the dissipation of heat and fluorescence as it was observed to affect the major protective carotenoids and antioxidants [40][41][42]. Commonly, the chlorophyll content of drought-sensitive varieties showed a more decline than those of drought-tolerant varieties [39,42]. The data presented in this study contradicts this description because, in chlorophyll a, the drop is moderate, while in chlorophyll b does not show a statistical difference between tolerant and sensitive sugarcane plants, reflected in the level of chlorophyll a + b, which does not shows the statistical difference between tolerant and sensitive sugarcane plants. It should be noted that Devi et al. [42] did not measure the concentration of chlorophyll by a standard method, instead measured the intensity of green through the SPAD, a very old technique that is currently in disuse.
Net photosynthesis is the preferred pathway for electrons under optimal conditions; however, P N is extremely sensitive and decreases more rapidly than ETR. Under stress condition and this promotes an increase in the ETR/P N ratio, which indicate that energy compromise takes place, and thus, it can be linked to ROS production [43]. ETR/P N ratio reflects the energy transfer to P N and an increase in this ratio, as noted in water-stressed sugarcane, indicates that alternative electron sinks most likely increased due to the stressful conditions [43]. When ETR/P N is correlated with g s , a lesser ETR/P N ratio denotes better use of electrons for photosynthesis, indicating high efficiency of radiation use by carbon assimilation reactions [18,43]. In another way, the higher ETR/P N ratio led to inferring that other pathways than photosynthesis using the electrons that are being transported by the photosynthetic electron transfer system (PETS) to NADP + . Figure 12 clearly shows this hypothesis, since the tolerant varieties (RB92579 and RB867515) showed higher ETR/P N values, even with low g s values, while the sensitive varieties (RB855536 and RB72454) showed intermediary-to-high ETR/P N at the expense of higher g s . Dos Santos [20] suggests that lower KL with lower sap flow may be translated to stomata closing since the plants, even if the atmosphere is dry and with very negative water potential or more, accumulated ABA. ABA not only directly mediates stomatal closure (chemical signaling) but also systemically decreases (KL) upstream of the stomata [44]. This work showed that KL has a strong influence on gas exchange, significantly and positively influencing PN (r = 0.843; p = 2.19 × 10 −44 ), gs (r = 0.885; p = 2.18 × 10 −54 ), and E (r = 0.762; p = 2.27 × 10 −31 ), while it is reduced by lowering the Ψw (r = 0.262; p = 0.03) and rainfall (r = 0.595; p = 1.85 × 10 −3 ). The lowering of KL in the same proportion is hardly discussed previously in this specie [20].
The main sap flow in 24 h tended to be more favorable to drought tolerant-like (RB92579 and RB86715), which shows 265.95 g H2O plant −1 h −1 , while the sensitive-like drought (RB855536 and RB72454) varieties showed 189.5 g H2O plant −1 h −1 , a reduction of 28.8% in sap flow around 24 h of measurement. Noteworthy that even in the rainy season, the anticipation of drop sap flow in the RB855536 should be a safe water strategy, corroborating other measurements that pointed to this variable as more intermediary than sensitive to water deficit. A wide range of literature covers that high VPD commands the water loss by leaf through stomatal conductance or transpiration [18,38,45,46] that leads to the increased sap flow. If it is true, many correlations may be translated as VPD and gs (r = −0.638; p = 1.14 × 10 −19 ) or E (r = −0.286; p = 0.019), impacting to lower KL (r = −0.532; p = 4.67 × 10 −13 ), which is provoked lower gs (r = 0.885; p = 2.18 × 10 −54 ), that along soil-plantatmosphere was provoke a lesser Ψw (r = −0.857; p = 2.77 × 10 −47 ). These correlations confirm our hypothesis that higher VPD supports higher stomata pressure (r = 0.937; p = 5.07 × 10 −96 ) that requires higher KL, which should exert strong pressure on the sap flow; but, this does not happen because the soils are dry with very negative water potential (mean Ψw = −2.14 MPa and −1.31 MPa, respectively, in dry and rainy plants evaluated in the MD). These features may lead to cavitation, which leads to inhibit the sap flow in the cavitated vessel.
The Fv/Fm ratio is often used as an indicator of the degree of photoinhibition in PSII [47]. We believed that the lower predawn Fv/Fm, particularly in dry season plants, may be partially associated with sustained xanthophyll cycle-dependent thermal energy dissipation, as previously reported in Coffea arabica [41], Alpinia oxyphylla [48], Cicer arietinum L. [49], and Spondias tuberosa [50] plants. In this study, we only measured total Dos Santos [20] suggests that lower KL with lower sap flow may be translated to stomata closing since the plants, even if the atmosphere is dry and with very negative water potential or more, accumulated ABA. ABA not only directly mediates stomatal closure (chemical signaling) but also systemically decreases (KL) upstream of the stomata [44]. This work showed that KL has a strong influence on gas exchange, significantly and positively influencing P N (r = 0.843; p = 2.19 × 10 −44 ), g s (r = 0.885; p = 2.18 × 10 −54 ), and E (r = 0.762; p = 2.27 × 10 −31 ), while it is reduced by lowering the Ψ w (r = 0.262; p = 0.03) and rainfall (r = 0.595; p = 1.85 × 10 −3 ). The lowering of KL in the same proportion is hardly discussed previously in this specie [20].
The main sap flow in 24 h tended to be more favorable to drought tolerant-like (RB92579 and RB86715), which shows 265.95 g H 2 O plant −1 h −1 , while the sensitive-like drought (RB855536 and RB72454) varieties showed 189.5 g H 2 O plant −1 h −1 , a reduction of 28.8% in sap flow around 24 h of measurement. Noteworthy that even in the rainy season, the anticipation of drop sap flow in the RB855536 should be a safe water strategy, corroborating other measurements that pointed to this variable as more intermediary than sensitive to water deficit. A wide range of literature covers that high VPD commands the water loss by leaf through stomatal conductance or transpiration [18,38,45,46] that leads to the increased sap flow. If it is true, many correlations may be translated as VPD and g s (r = −0.638; p = 1.14 × 10 −19 ) or E (r = −0.286; p = 0.019), impacting to lower KL (r = −0.532; p = 4.67 × 10 −13 ), which is provoked lower g s (r = 0.885; p = 2.18 × 10 −54 ), that along soil-plant-atmosphere was provoke a lesser Ψ w (r = −0.857; p = 2.77 × 10 −47 ). These correlations confirm our hypothesis that higher VPD supports higher stomata pressure (r = 0.937; p = 5.07 × 10 −96 ) that requires higher KL, which should exert strong pressure on the sap flow; but, this does not happen because the soils are dry with very negative water potential (mean Ψ w = −2.14 MPa and −1.31 MPa, respectively, in dry and rainy plants evaluated in the MD). These features may lead to cavitation, which leads to inhibit the sap flow in the cavitated vessel.
The F v /F m ratio is often used as an indicator of the degree of photoinhibition in PSII [47]. We believed that the lower predawn F v /F m , particularly in dry season plants, may be partially associated with sustained xanthophyll cycle-dependent thermal energy dissipation, as previously reported in Coffea arabica [41], Alpinia oxyphylla [48], Cicer arietinum L. [49], and Spondias tuberosa [50] plants. In this study, we only measured total carotenoids, where plants of the RB855536 had the lowest F v /F m ratio but linked to the highest Car/Chl ratio (0.65 ± 0.02) which was 1,6-, 2.1, and 2.1-fold higher than RB92579, RB72454, RB867515 varieties, respectively. The correlation between Car/Chl ratio and rainfall was strongly negatively correlated (r = −0.741; p = 1.36 × 10 −8 ). In the same way, the Car/Chl ratio was strongly and negatively correlated with P N (r = −0.698; p = 1.08 × 10 −13 ). Summarily, it can be considered that RB855536 genotype, even with a lower F v /F m ratio, did not have their photosynthetic and hydraulic apparatus negatively affected, since the predawn potential was 62% higher in the dry season in comparison with the rainy season, which in part can be explained by the anticipation of the decrease in sap flow ( Figure 5). Thus, during the hottest moments of the day, these plants reduce their sap flow, causing a stomatal closing affecting the high magnitude in dry season plants (Figure 8). The mechanisms that command this modulation could have the higher concentration of Car/Chl ratio, or the higher Φ PSII , since higher the Φ PSII reduce the mandatory dissipate heat, promoting a security thermal dissipation. These hypothesis gives more strength by positive correlation with F v /F m ratio and Chl a (r = 0.522; p = 7.03 × 10 −7 ). The antioxidative enzymes also corroborates this hypothesis, where F v /F m ratio and SOD (r = −0.758; p = 3.71 × 10 −16 ), CAT (r = −0.546; p = 1.37 × 10 −7 ) and proline (r = −0.694; p = 9.53 × 10 −13 ) had strong negative correlations.
Less rainfall impacts more synthesis of proline (r = −0.802; p = 9.94 × 10 −11 ), where proline synthesis is provoked by less P N (r = −0.828; p = 2.70 × 10 −21 ), g s (r = −0.592; p = 7.58 × 10 −9 ), ETR (r = −0.841; p = 1.68 × 10 −22 ), Φ PSII (r = −0.668; p = 1.25 × 10 −11 ) and had a very strong participation in the increase activities of antioxidant enzymes, like SOD (r = 0.722; p = 1.96 × 10 −37 ) and CAT (r = 0.841; p = 3.93 × 10 −61 ). In accordance with Naheed et al. [39], the degree of tolerance to drought stress is strongly correlated with the increase in the concentration of proline. In the present study, the sugarcane varieties RB92579 and RB855536 showed an increase in proline to~53% and~72%, translated to low production of MDA than the other two sugarcane varieties. Here, we reported that under water stress, some proteins can be cleaved, increasing the level of amino acids that could be important in osmoregulation or used as building blocks to generate other types of proteins, such as membrane proteins, defense proteins [18,51], or used as a carbon source to maintain the cellular respiration and ATP synthesis [18,51]. Furthermore, contrary to Dos Santos et al. [20], here, we described the most recent version of the role of proline, not as an osmoprotectant but as a stress preventer as reported by others [2,39,52] that defends the proline synthesis as a consequence of stress. In this way, the synthesis of proline comes from an alternative route of metabolizing glutamate, which is phosphated by 1 mole of ATP and reduced to proline from the consumption of 2 moles of NADPH, by activity of glutamate dehydrogenase and pyrroline-5-carboxylate reductase, which uses one mole of FADH 2 as a cofactor. Consequently, a stressed plant tends to proportionally reduce the stomatal opening, which limits the input of CO 2 . With less CO 2 in the system, the Calvin-Benson cycle is compromised, and the ATPs and NADPHs that formed and continue to form through the PETS remain in their reduced form [18]. At the end of the PETS, ATP and NADPH were reduced and remained in the reduced form because the Calvin-Benson cycle is operating at the lowest potential oxidizing less ATP and NADPH. Once these molecules are reduced, there are no oxidized forms of them to receive more electrons that come from the PETS, and these would be quickly captured by O 2 , forming the terrible reactive oxygen species (ROS) [53]. So, proline may have strong importance because it acts as stress prevention than osmoregulation.
Superoxide dismutase forms the very first line of defense against reactive oxygen species. Naheed et al. [39] describe the higher activity of SOD in drought-sensitive sugarcane plants as described in this study. Strong evidence for the role of chloroplast SOD has been shown from the interruption of PETS by methyl viologen [39]. After that, intense SOD activity was verified, both directly measuring the activity of the enzyme and indirectly by the production of H 2 O 2 . The latter can inactivate PETS, generating even more H 2 O 2 . In accordance with these authors, the restriction of electron transport from PSII to NADP + may cause the total interruption of PETS by peroxidation induced by ROS. The highest antioxidative enzymes are inversely proportional to the damage caused by oxidative stress in sugarcane as previously described [39] and reported in this study were less precipitation increase in SOD (r = −0.756; p = 4.51 × 10 −9 ), CAT (r = −0.775; p = 1.03 × 10 −9 ), and APX (r = −0.248; p = 0.011) activity. To prevent the deleterious effect of H 2 O 2 produced by SOD, CAT and APX promote the dismutation of this molecule as H 2 O + O 2 . To plant defense, the activity of APX is a better choice than CAT [54] because while CAT only dismutates H 2 O 2 into H 2 O + O 2 , APX needs the reduced ascorbate as a cofactor. In the recycling of ascorbate, the cycle will need monodehydroascorbate, which is reduced to ascorbate with the spend of 1 mole of NADPH to contribute as another line of defense as it acts to prevent ROS production. So, the APX act in two ways: promoting reoxidation of NADPH and dismutating H 2 O 2 . Thus, it is expected that the activity of APX is more efficient in plants under severe stress since the synthesis of ascorbate is a costly process for the plant. In plants with a lower degree of stress, the APX activity is expected to be overwhelmed by the CAT activity. This hypothesis is also described here since genotypes RB867515 and RB72454 showed an increase in CAT activity by 34.3% and 15.8%, while APX activity in these same two varieties was decreased by 7.3% and 28.8%.

Morphological, Anatomical, and Ultrastructural Features
Circumstantial evidence that the lower photosynthetic rates verified in the tolerant varieties may not be exclusively linked to the stomata, but the greater leaf thickness, around 30%, in the tolerant varieties compared to the related varieties. However, less transpiration in the tolerant varieties may have been caused by the presence of bulliform cells on the adaxial surface of the epidermis (Figure 9) compared to sensitive varieties. Larger bulliform cells facilitate leaf curling [55], allowing the abaxial, the more hairy, reflective, and stomatabearing surface to become sun-facing ( Figure 11) [55,56]. Our SEM micrograph is supported by the description by Kadioglu et al. [55] since the abaxial face, easily exposed in the case of tolerant varieties, with a marked presence of epidermal hair trichomes and hook-shaped trichomes (Figures 10 and 11). Furthermore, the more prominent arrangement of bulliform cells in the tolerant varieties presented in this study corroborates previous studies with the species [57,58]. These authors described that bulliform cells are often located slightly above the level of other epidermal cells. Regardless of this, it is known that the size and height of epidermal cells and their accessory structures are related to ecological aspects, presenting responses to environmental changes [55].
The cuticle that covers the leaf is nearly impermeable to CO 2 , so the main input of CO 2 into the leaf is the stomatal pore, hampered by epicuticular waxes and hydrophobic by nature. Even with the intense activity of carbonic anhydrase (EC 4.2.1.1), CO 2 has difficulty in permeability to the substomatal chamber. The lower stomatal conductance of the tolerant varieties may be due not only to the lower degree of stomatal opening in these varieties but also to the remarkable deposition of waxes around and on the stomatal pore, which is significantly higher in the tolerant varieties compared to the sensitive ones as previously reported in sensitive and tolerant wheat varieties [59]. Studies with Nicotiana glauca [60], Carex hirta [61], and sugarcane [62] also described the influence of wax deposition on the stomatal pore. The cuticular layer can vary widely in its ultrastructure and the number of cuticles deposited on leaves [63]. The cuticle consists of two types of lipophilic substances: cutin and cuticular waxes (intracuticular or extracuticular). These waxes can crystallize and leak out of the cuticle giving the surface a grainy or ornate appearance [64].
Zan et al. [4], studying 236 sugarcane accessions, including 13 Brazilian accessions, described many characteristics to search elite accessions, to be cultivated because they present characteristics of agronomic interest and from the point of view of molecular biology. However, de Paula et al. [65], studying the establishment of gene pools for systematic heterosis exploitation in sugarcane, indicated that 7 of 13 accessions studied by Zan et al. [57] are no longer cultivation in Brazil, data that is corroborated by other studies [6,[66][67][68]. While RB83-5054 cultivation in Brazil is currently inadvisable [65], both RB85-5113 and SP80-181 currently is used as mother parents in crop breeding programs, and SP81-3250 is currently used as father parents in crop breeding programs. Table 4 summarizes this information and shows the expected sugarcane Brazilian production of SP80-1816, SP80-3280, SP81-3250, and SP83-2847 genotypes in 2021 when compared to 2020 cultivation. The genotypes RB92579 and RB867515 have been detached, with 3% and 26.4% more harvested in 2021 when compared to 2020 (Table 4). Apart from genotypes described by de Paula et al. [58], Carneiro et al. [69] describe another Brazilian cultivar with high agroindustrial yield and disease resistance.  [59]. 4 = In accordance with Cursi et al. [61]. 5 = In accordance with de Paula et al. [58]. 6 = RIDESA-BRASIL [70]. 7 = In accordance with Barreto et al. [62]. ♀= Currently used as mother parent in crop breeding. ♂= Currently used as father parent in crop breeding. Table adapted to data available in Matsuoka et al. [59], de Paula et al. [58], and Cursi et al. [61].

Conclusions
As signaling of stress-tolerant plants is the very first priority of the researchers. In the current study, employing sugarcane, the physiological strategy involving gas exchange, fluorescence parameters, sap flow, water potential, and biochemical and enzymatic features linked to drought stress were evaluated. In addition, morphological, anatomical, and ultrastructural features allow us to assess how the growth and tolerance of sugarcane behave to low rainfall and how these same variables are modulated when the plant receives a higher rainfall.
So, our data support that sugarcane productivity is largely affected by photosynthesis, and the identification of the most photosynthetically efficient variables under water stress conditions is extremely important. So, sugarcane plants must be continuously improved to produce more sugar per square meter and be less affected by environmental stressors.
The data presented here support the idea that the RB92579 genotype is drought tolerant, as well as RB867515, and introduce the RB855536, as a tolerant-like, due to its high capacity to increase its antioxidant system to cope and protect the PETS from ROS, disruption, and less level of MDA, even so, many characteristics lead RB855536 to drought tolerance.
The morphological features gave us confirmation that sugarcane under water stress produced smaller leaves, with shorter plants and a large increase in the bulliform cells that help the leaf curl, which allows it to expose its abaxial face, which has a higher concentration of reflective trichomes and hook-shaped trichomes.
The anatomical features confirmed that all physiological and biochemical features could be translated to an anatomical modulation, like thinner leaves, as a strategy to promote latent heat loss. Linking this idea with gas exchange parameters permits us to infer that the lower net photosynthesis in dry season plants could be a result of thinner leaves because the gas exchange parameters are given in m −2 s −1 , so a thinner leaf, with less photosynthetic tissue per square meters, could be underestimating of net photosynthesis.
We confirmed the hypothesis that the sugarcane, as well as the connected soil-plantatmosphere plant, could be modulated by some features. To complete this idea, we measured each degree of each physiological, biochemical, morphological, anatomical, and ultrastructure feature that has an plant acclimation to stress. To our surprise, all these features contribute to plant plasticity and modulation of water deficit. This dataset can be uploaded in Supplementary File S2.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/horticulturae8111069/s1, Supplementary File S1: The Supplementary S1 shows us day-by-day rainfall registered in Campo Alegre city-between 1999, January and 2020, December; Supplementary File S2: Plasticity index (PI) of physiological, biochemical, morphological, and ultrastructural features measured in four varieties (RB92579, RB867515, RB855536, and RB72454) of sugarcane. Numbers in bold denote the mean (±SE). Supplementary Data File: Correlation matrix showing all pairwise correlations between physiological, morphological, and ultrastructure features. The color code following the scale. An asterisk (*) denotes significant correlation (p < 0.05).
Author Contributions: All authors contributed equally. All authors have read and agreed to the published version of the manuscript.