Different leaf water dynamics post- anthesis affects the final kernel weight in wheat Diferentes dinámicas de agua en hojas pos-antesis afecta el peso final del grano de trigo Diferentes dinâmicas da água em folhas pós-antese afeta o peso final do grão de trigo

Drought stress during wheat (Triticum aestivum L.) grain filling can affect the grain yield and its quality; effects that in some regions can be increased by the event La Niña, the cold phase of the climatic phenomenon El Niño Southern Oscillation. The aim was to evaluate if the evolution of relative water content (RWC) and Editor

stomatal conductance (gs) can be related to grain yield and its quality. Four spring wheat cultivars grown in the southern cone of South America were subjected to two levels of irrigation applied after anthesis, a wellwatered treatment (Control), and another with 50 % of the Control (Stress). RWC and gs were determined weekly during the stress period and agronomic traits at harvest. With a mean grain yield decrease of 14.8 % (P<0.05), two types of response to water deficit were identified for RWC and gs. Two cultivars (Biointa 1001 and LE2333) had slower RWC and gs decreases, which were related to a lower ratio of kernels to aerial biomass (P<0.05), and to a non-significant effect in kernel weight. While, in the other cultivars (LE2249 and LE2331), the water deficit caused a more rapid RWC and gs decrease related to a higher ratio of kernels to aerial biomass, and a significant kernel weight decrease (10.4 and 20.7 %, respectively). The latter cultivars, that had more rapid leaf dehydration under water stress, may also have been source-limited during the grainfilling period.

Introduction
Wheat (Triticum aestivum L.) farm yield potential is increasing at a relative rate of 1.0 % per year (1) ; an insufficient increase to feed a population that annually grows at the same rate (2) . One solution to meet this demand for food is to expand the agricultural frontier, but the inclusion of soils with growth constraints can lead to higher risks of drought stress. Other risks for food security are the effects of 'climate change' adjudicated mainly to global warming (3) . From the 1970s in the tropics and subtropics, periods of drought have been more frequent, intense and in broader geographical extents (4) . Particularly, in South America when La Niña (the cold phase of the ENSO climate phenomenon) event occurs, the risk of drought stress increases in a greater wheat crop area during grain filling.
Water deficit during grain filling affects the kernel weight (5)(6)(7) (8) and the grain industrial quality, largely due to that nitrogen and carbohydrates deposition rates are affected differently (9) . Severe stressors during this period or just before anthesis, may also affect the number of fertile tillers and the number of kernels (10) . These yield components are affected differently by water deficit, depending on the relationship between the plant water consumption and the initial level of available water, the time in which the stress event is initiated, and according to the duration of that stress period.
Severe stress events affect the rate of grain filling and the length of the grain-filling period (11) ; however, mild stress events may lead to increase the rate of grain filling due to a greater sink activity, which would enhance the remobilization of carbon reserves from vegetative tissues to the kernels (12) . 'Mild water stress' refers to the level of water deficit to which the plant metabolism can be recovered after a certain period of stress. One way to identify mild water stress is the approach of Flexas and others (13) , who sustained that the level of photosynthetic activity can be recovered one day after re-watering if the previous values of gs were above 150 mmol H2O m -2 s -1 .
There is agreement in the literature that kernel size in small cereals, in the absence of water stress and other adversities (biotic and abiotic), is a yield component limited by the sink capacity of the kernel (14) (15) . Normally, under drought stress conditions the sink strength in anthers and ovaries is irreversibly reduced; although, there is germplasm that can maintain the number of kernels, but affecting the kernel size (16) . Similarly, when germplasm of high yield is grown in more restrictive environments not only the sink strength would be limiting, but also the source capacity (17) .
Cultivars of slightly early anthesis and delayed but rapid leaf senescence, promote a higher rate of grain filling but of shorter duration, a higher rate of grain water absorption, and a maximum grain water content (18) . Senescence is a process genetically controlled but also affected by the sinksource relationship and environmental stresses (19) . Hence, cultivars of delayed but rapid leaf dehydration under mild water stress could be also little affected in their kernel weight. Even though leaf senescence is more related to the green leaf area index (LAI), and wilting is more related to the interception of radiation per unit of LAI (20) , droughtinduced wilting may execute prematurely the onset of senescence (19) .
The genetics of drought stress resistance is complex and related to other abiotic stresses as salinity and frost temperatures (21) (22) . Largely due to this reason, several physiological traits have been studied; among them, the evolution of relative water content (RWC) as an indicator of wilting and of the plant water status (23) (24) . The normal physiology of the plant begins to be affected when the RWC decreases below 80 % (25) . Even though RWC may underestimate the water content when is measured in plants of high osmotic adjustment, it is a useful indicator to carry out easy measurements of the plant water status (24) . Among other studies, RWC has been used as a physiological trait to study adaptive strategies of cultivars under stress factors affecting the yield formation (10) (26) .
As an indicator of stomatal adjustment, another trait for studying adaptive strategies of cultivars to drought stress is the stomatal conductance (gs). Even though the gs is largely influenced by environmental conditions, as radiation; air temperature; vapour pressure deficit (VPD); and leaf water potential (27) (28) (29) (30) , it has been proposed as an important trait because cultivars of scarce regulation of gs show lower values of water use efficiency (WUE) (31) . Cultivars of lower WUE have also been identified in durum wheat, where landraces of high dry matter accumulation but moderate final grain yield, did higher use of water which was related to canopies of high LAI and gs. While, modern cultivars with canopies of similar gs but of lower LAI, did a lower use of water in early stages, which led to they finally had higher WUE (32) . Similarly, mediated by stomatal adjustment, conservative transpiration rates in early stages of wheat may outweigh the negative effect of decreased photosynthesis under drought stress conditions (33) (34) .
Water deficit during grain filling would promote different evolutions of RWC and gs that would cause different cultivar responses. Wheat cultivars that could maintain high values of RWC mediated by stomatal adjustment, would not necessarily affect the allocation of assimilates to the kernel. The aim was to study if the evolution of RWC and gs of four spring wheat cultivars under two levels of irrigation applied after anthesis can be related to kernel weight and the ratio of kernels to aerial biomass. The resulting information could be useful to identify wheat cultivars that under terminal drought stress could maintain the kernel weight, as a determinant component of grain yield and industrial quality.

Plant material
Four spring bread wheat cultivars, Biointa 1001, LE2249, LE2331, and LE2333 were evaluated under two levels of irrigation. These cultivars, that have been cultivated in the southern cone of South America, were chosen because in a previous study they showed different resource allocation to reproductive structures. In that study, LE2333 had a high number of kernels per spike, but a low number of spikes per plant and a low kernel weight together with LE2249. While, the cultivar LE2331 had a high number of spikes per plant, and together with Biointa 1001 a high kernel weight (35) .

Experimental design
The experiment was located at the Experimental Station INIA La Estanzuela, Colonia, Uruguay (34º20'15" S; 57º41'29" W), in a greenhouse equipped with fan and pad evaporative cooling system. The experiment was arranged in a randomized complete block design with three replications. A two-factorial arrangement of four cultivars and two levels of irrigation were evaluated. The experimental unit consisted of 24 plants planted in six pots (four plants per pot) placed on a plastic tray of 0.60×0.40×0.10 m. The pots were PVC tubes of 0.16 m of inner diameter and 0.30 m of height. A metal grid of fine mesh was placed in the base of each tube to contain the substrate. The substrate used was a mixture (1:1:1) of a silty clay loam soil, sand, and organic substrate (BioFer Almácigos, Riverfilco S.A., Montevideo, Uruguay).
The experiment was sown (27 June) with eight seeds per pot and thinned to four seedlings per pot after implantation. The plants received 8 to 6 h of artificial illumination with sodium vapour lamps of high pressure (SON-T 400 W, Philips, Belgium).
Before the anthesis stage, the experiment was watered once or twice a week according to the average water consumption of the four cultivars; this is that each cultivar was irrigated with the same amount of water ( Figure 1). When 50 % of plants by cultivar and block reached the anthesis stage, two levels of irrigation were applied: a treatment well-watered (Control) and other watered with 50 % of the Control (Stress). Both irrigation treatments were applied until the experimental units in the Control treatment reached maturity.
The initial amount of water to irrigate the Control treatment was defined according to the maximum water holding capacity of a sample of an ovendried substrate (12.8 % of initial moisture content). The dried sample substrate was put in an extra pot and was slowly watered until saturation. After 24 h, once the substrate in the pot stopped draining, it was weighted to calculate its maximum water holding capacity (1 L per pot).
The water consumption of the Control treatment was estimated following the water consumption of two check trays that were sown with the four cultivars (one cultivar per pot and tray, respectively). The check trays were placed on each side of the experiment (blocks 1 and 3, respectively). The water consumption of the check trays was calculated by the difference between the pot weights at maximum water holding capacity, and the pot weights immediately before each irrigation treatment.
At tillering stage GS 2.2 (36) , each pot was fertilized with 0.5 mL of N-P-K (12-8-5, foliar fertilizer NPK micronutrients, Industria Sulfúrica S.A., San José, Uruguay) diluted in 0.42 L of water. At the end of stem elongation (GS 3.7 to 3.9), the same dose of fertilizer was repeated, but applied on the trays and diluted in 5 L per tray.

Measurements
The RWC = [Fresh weight -Dry weight] × 100 / [Turgid weight -Dry weight] was determined according to Barrs and Weatherley (23) once a week from anthesis to physiological maturity. The gs (mmol H2O m -2 s -1 ) was measured with a steadystate diffusion porometer (SC-1, Decagon Devices, Inc., Pullman, WA) during the same period and frequency than the RWC. The measurements of both traits were always carried out in the morning and before irrigation.
The RWC was determined in sections (2 to 3 cm) of three leaf blades per experimental unit. The blades were sampled from leaves below the flag leaf (FL-1) of the main stem of three plants of different pots. The fresh, turgid and dry weights of the three blade sections were determined. The gs was measured in the middle of the adaxial side of the flag leaf (FL) of the main stem in three plants of different pots.
After harvest, sun-dried aerial biomass per plant (biomass per plant) (g), spikes per plant, kernels per spike, kernels per plant, single kernel weight (mg), grain yield per plant (g), the ratio of kernels to aerial biomass (kernels/biomass ratio) (g -1 ), harvest index (grain yield/aerial biomass), irrigation WUE on yield basis [IWUEYI (g L -1 ) = grain yield/irrigation water], IWUE on aerial biomass basis [IWUEBI (g L -1 ) = aerial biomass/irrigation water], and the kernel protein concentration (%) were determined. The average number of kernels per plant was determined by counting the kernels in each spike of all plants of each experimental unit (24 plants). The kernel weight was determined weighting all kernels of each experimental unit divided by the total number of kernels. The irrigation water, the denominator of the IWUEs, was not discriminated in transpired water and evaporated water. The kernel protein concentration (moisture 13.5 %) was determined by the Laboratory of grain quality of the Experimental Station INIA La Estanzuela, using a spectrophotometer Serie 6500 (FOSS NIRSystem Inc., Silver Spring, MD, USA) calibrated by the technique of Kjeldahl (approved method 46-12, AACC 2000).

Data analysis
For all traits, except for RWC and gs, the statistical model (M1) used was: Yijk = µ + αi + γj + (αγ)ij + βk + eijk, where µ, α, γ, β, and e, are the mean, irrigation treatments, cultivars, blocks, and the errors, respectively. To identify if there were treatment and interaction effects, a significance level of P<0.05 was considered for all traits. To identify mean differences among cultivars, the Tukey's test and the Tukey-Kramer test were performed for balanced and unbalanced data, respectively. Also, orthogonal contrast tests were performed to discriminate interactions between cultivars and irrigation treatments.
For the repeated measurements in time of RWC and gs, the statistical model (M2) used was: where the first four terms are the same that were mentioned above (M1), s is the between-subjects random effect, τ is the effect of days of initiated irrigation treatments (DIIT), and e(a) and e(b) are the errors a and b, respectively. Additionally, in the statistical model for gs, the VPD was included as a covariate. The VPD was calculated with the empirical exponential model of Prenger and Ling (37) , employing hourly records of temperature and air relative humidity at the time that the gs measurements were performed ( Figure 1).
For traits RWC and gs, the covariance structure with different models was analysed. For the estimation of means, considering the likelihood ratio test (LRT) and the Akaike and Bayesian information criteria (AIC and BIC, respectively), it was chosen the spatial power law [SP(POW)] model with the Kenward-Roger adjustment. This model is recommended for repeated measurements taken at unequal times, correlations that decline over time, unbalanced data and multiple random effects (38) . Also for both traits, the adjusted means of the repeated measurements that were carried out on the same day were paired and correlated (according to Pearson's r coefficient). Finally, all the adjusted means of RWC and gs were correlated with those traits affected by Stress treatment.
Moreover, the repeated measurements in time of RWC were fitted iteratively by least squares to the exponential model: bo + A [1 -exp(-k DIIT)] (39) . The model indicates how the RWC decreases during grain filling and how it accelerates towards the senescence. The first parameter of the model (bo) indicates the RWC at anthesis; the second (A), the linear rate of decrease of the RWC; and the third parameter (k) indicates the RWC exponential decrease. The greater the water stress, the larger are the parameters A and k, whereas greater is the influence of A respect to k in the RWC decrease. To identify differences between treatments, the parameters estimated by the fitted model were subjected to the t-test of Welch-Satterthwaite.
All data were analysed with the PROC MIXED statistical procedure of the software SAS ® version 9.2 (SAS Institute, Cary, NC, United States) (40) . In turn, the data of RWC were fitted to the model of Orskov and McDonald (39) using the software InfoStat ® version 9.0 (InfoStat Group, FCA, National University of Cordoba, Argentina) (41) .

Transparency of data
Available data: The entire data set that supports the results of this study was published in the article itself.

Results
The Control treatment was irrigated with a cumulative mean of 2 L per plant, whereas the Stress treatment was irrigated with 1.5 L per plant, from sowing to maturity date (Figure 1). The cultivars LE2331, Biointa 1001, LE2333 and LE2249, were irrigated with 1.79, 1.72, 1.70 and 1.66 L per plant, respectively, but only the means of LE2331 and LE2249 differed significantly (P<0.05) related to differences in anthesis time, length of the grain filling period, and to the different VPD values during that time (Figure 1). Interaction between cultivars and irrigation treatments was not identified (P=0.5674, data not shown); this, due to after anthesis all cultivars in the Stress treatment were always irrigated at 50 % of the Control.  Table 1. Least square means of two irrigation treatments (Control and Stress) and four spring wheat cultivars (Biointa 1001, LE2249, LE2331, and LE2333). The traits observed are spikes per plant (SP PL -1 ); kernels per spike (KE SP -1 ); kernels per plant (KE PL -1 ); grain yield per plant (YI PL -1 ); dry matter of aerial biomass per plant (BI PL -1 ); the ratio of kernels to aerial biomass (KE BI -1 ); and irrigation water use efficiencies, on yield (IWUEYI) and aerial biomass (IWUEBI), respectively.

Treatments
Least square means On the contrary, grain yield per plant, biomass per plant, and the kernels/biomass ratio were significantly affected by Stress treatment, but while grain yield and biomass per plant were reduced, 14.8 and 10.9 % respectively, the kernels/biomass ratio was increased 5.9 % (P<0.05; Table 1). For the above traits, statistically significant differences among cultivars were also observed (P<0.01). The cultivar LE2331 had higher values of grain yield and biomass per plant, although they did not differ significantly from the values obtained by Biointa 1001 (P<0.05). Concerning to the kernels/biomass ratio, the cultivars LE2249 and LE2331 had higher values than Biointa 1001 and LE2333 (P<0.05). No significant interactions were identified for the three mentioned traits (data not shown).
No interactions were also identified for IWUEYI and IWUEBI (data not shown), but both traits had a significant increase due to the Stress treatment, 12.9 and 17.7 %, respectively (P<0.05; Table 1). Significant differences among cultivars (P<0.01) showed that LE2331 had higher values in both IWUEs, although they did not differ significantly from the values obtained by Biointa 1001, and by LE2249 in the case of IWUEYI (P<0.05).
A significant interaction between cultivars and irrigation treatments was observed for kernel weight (P=0.0125), harvest index (P=0.0014), and kernel protein concentration (P=0.0286) ( Table 2). Only the cultivars LE2249 and LE2331 had a significant decrease in their kernel weight because of the Stress treatment, 10.4 % (P=0.0389) and 20.7 % (P=0.0002), respectively. With respect to the harvest index, only Biointa 1001 and LE2331 had a significant decrease, 4.0 % (P=0.0335) and 10.1 % (P<0.0001), respectively. On the other hand, except for LE2333 that had a decrease of 4.6 % (P=0.0520), the kernel protein concentration had a non-significant increase in the other three cultivars due to the Stress treatment.
At anthesis (DIIT 0), there were no significant differences in RWC and gs among cultivars and between both irrigation treatments (data not shown). The cultivar means of RWC ranged from 87.0 to 91.5 % (Figure 2), while the means of gs ranged from 136 to 207 mmol H2O m -2 s -1 in relation to the daily evolution of the VPD in the greenhouse (Figures 1 and 3). For this reason, the inclusion of the VPD as a covariate into the model had a significant effect (P=0.0024, data not shown) in the adjustment of the gs means.  The correlation analysis for all treatment means indicated a significant positive correlation between the repeated measurements of RWC and gs (r=0.32; P=0.0227; n=50, data not shown). In turn, until 27 DIIT significant and positive correlations were also identified between the RWC measurements and kernel weight, whereas significant but negative correlations were identified for the kernels/biomass ratio (Table 4). Less clear correlations were identified through the time between the gs measurements and kernel weight, but they were significant and positive on days 5, 14 and 20 after initiated the irrigation treatments. Finally, between gs and the kernels/biomass ratio significant but negative correlations were also identified on the same days (Table 5).

Discussion
In this work, the number of kernels per plant was not affected by the Stress treatment (Table 1), a yield component that can be affected by severe stresses and/or by stresses that begin before the anthesis stage (5) (10) . However, the progressive water stress applied after anthesis caused different cultivar responses for harvest index, kernel weight, and kernel protein concentration ( Figure 1; Table  2). Those responses were related to cultivar differences in the formation of the grain yield and of the biomass allocation (Table 1), and to cultivar differences in the RWC and gs evolutions in response to both irrigation treatments (Figures 2  and 3).
In this regard, it is observed that cultivars with a greater number of spikes per plant and kernels/biomass ratio (LE2331 and LE2249), were the only ones that had a significant decrease of the kernel weight (Tables 1 and 2). According to Dreccer and others (42) , the decrease of the kernel weight in genotypes of high number of tillers could be related to a lower level of stem reserves of water-soluble carbohydrates. In addition, only the harvest index of cultivars LE2331 and Biointa 1001 decreased significantly due to the Stress treatment (10.1 and 4.0 %, respectively). So, it could be interpreted that the three mentioned cultivars were affected because of a co-limitation sink-source.
Positive correlations identified between RWC measurements and kernel weight indicate that Stress treatment caused leaf dehydration (Table  4), which ultimately affected the allocation of  Table 4. Correlations between repeated measurements in time of relative water content (RWC) with kernel weight, yield per plant, biomass per plant, the ratio of kernels to aerial biomass (kernels/biomass), harvest index, and irrigation water use efficiencies, on yield (IWUEYI) and aerial biomass (IWUEBI). The data correspond to four spring wheat cultivars (Biointa 1001, LE2249, LE2331, and LE2333) subjected to two irrigation treatments (Control and Stress) after anthesis. *; ** Significant at the 0.05 and 0.01 probability level, respectively. Table 5. Correlations between repeated measurements in time of stomatal conductance (gs) with kernel weight, yield per plant, biomass per plant, the ratio of kernels to aerial biomass (kernels/biomass), harvest index, and irrigation water use efficiencies, on yield (IWUEYI) and aerial biomass (IWUEBI).
The data correspond to four spring wheat cultivars (Biointa 1001, LE2249, LE2331, and LE2333) subjected to two irrigation treatments (Control and Stress) after anthesis.
biomass into the kernels of cultivars LE2249 and mainly LE2331 (Figure 2C-D; Table 2). Final kernel size depends on the number of endosperm cells -defined for wheat around 15 to 20 days after anthesis-, and of the final endosperm cell size (43) (44) . Values of RWC below 80 % promote a general metabolic inhibition; this is, lower activity of the sucrose phosphate synthase and protein synthesis related to a lower transpiration efficiency (25) (45) ; while, values below 70 % cause non-stomatal oxidative stress in the photosynthetic machinery (25) (46) . In that sense, in the Stress treatment the RWC of LE2249 decreased below 80 % after 11 DIIT, whereas for LE2331 the same occurred, but from the beginning of the period (Figure 2 C-D).
No significant correlations between kernel protein concentration and the measurements of RWC and gs were identified (data not shown). However, there was a tendency that as the kernel weight and harvest index decreased, the kernel protein concentration increased (Table 2). This would be related to the fact that the rate of nitrogen deposition in the kernel is reduced more slowly than the rate of carbohydrates deposition (9) .
Two types of contrasting responses in the evolution in time of RWC were identified: 1) LE2333 and Biointa 1001 that had a similar evolution of RWC in both irrigation treatments (Figure 2A-B), and 2) LE2249 and LE2331 that had an early decrease of RWC in the Stress treatment ( Figure 2C-D).
In one response, the cultivars of lower allocation of biomass to the kernels, LE2333 and Biointa 1001, maintained the RWC for a longer period of time under the stress conditions (Table 1; Figure 2A-B). However, Biointa 1001 would have sensed the stress earlier because it did not decrease the kernel weight, but decreased the harvest index (Table 2). Conversely, the cultivar LE2333 was not affected in both mentioned traits. Cultivars as LE2333, of delayed but rapid leaf dehydration at the end of the cycle, would promote higher rates of grain filling but of shorter duration (18) .
In the other response, the cultivars of higher allocation of biomass to the kernels (LE2249 and LE2331) (Table 1), would have sensed the deficit at the beginning of the stress period as a mild stress. The rapid loss of water would have been caused by scarce control of the stomatal adjustment ( Figures 2C-D and 3C-D). Mild stresses during grain filling would affect the kernel weight, but not other yield components because the effects would be delayed beyond the cell division phase (7) (47) . In that sense, a rapid decrease of the water status caused a decrease in the kernel weight of both cultivars. Thus, the applied progressive water stress during the grain filling would have promoted in these cultivars that kernel weight was finally limited not only by its sink capacity, but also by the source capacity of the plant (12) .
Despite the contrasting type of response in the evolution of RWC of Biointa 1001 and LE2331, both cultivars did not differ significantly in their grain yield per plant (Table 1; Figure 2B, 2D). A low control of the closure of stomata of Biointa 1001 caused that this cultivar had a rapid wilting of the basal leaves to maintain the RWC in the upper leaves ( Figures 2B and 3B). However, a lower volume of green LAI would explain the decrease of the harvest index (Table 2) through a reduced interception of radiation per unit of LAI (20) , and probably due to a reduced sink size, since under drought stress conditions the sink strength in anthers and ovaries can be irreversibly reduced (16) .
Stomatal responses shown by the four cultivars are consistent with the two types of response identified in the RWC evolution (Figures 2 and 3). Cultivars that had a rapid decrease of the RWC in the Stress treatment (LE2249 and LE2331), had a greater stomatal adjustment, as can be observed by the gs values located below the 1:1 relationship lines ( Figures 2C-D and 3C-D). The positive and significant correlation between all the repeated measurements of RWC and gs (r=0.32; P=0.0227; n=50, data not shown) was indicative that maintaining high values of gs is only possible if high water status is also maintained (33) .
During the grain-filling period, the gs of the cultivar LE2333 was scarcely affected by the Stress treatment, as was also the RWC evolution. A possible explanation for this is that the osmotic potential in the sink may not have declined below the optimum. So, an enhanced sucrose uptake and starch synthesis related to increased carbon remobilization from vegetative tissues may have occurred (7)(12)(18)(42) (48) . In relation to this, the grain yield per plant, biomass per plant, and kernels/biomass ratio of LE2333, did not differ significantly from the cultivar means of LE2249 (Table 1). However, the Stress treatment caused in LE2249 a significant decrease of the RWC, and consequently, a 10.4 % decrease in the kernel weight (Tables 1, 2 and 3; Figure 2).
Similar to WUE, the ratio between the CO2 assimilation rate and the transpiration rate (49) , the increase of the IWUEYI and IWUEBI caused by the Stress treatment means that per each unit of grain yield and aerial biomass, more efficient was the use of irrigation water made by the plants (Table 1; Figure 3). Cultivars of high WUE per se not necessarily should be associated with drought stress resistance. This, because they may have higher yield potential with moderate use of water, or may have moderate yield potential with reduced use of water (50) . In this regard, even though LE2331 and Biointa 1001 showed higher WUEs (Table 1), they would have used much more water if it had been provided. This was observed after each irrigation event when the available water in the trays decreased faster than the water in the trays of cultivars LE2249 and LE2333.

Conclusions
A progressive water stress applied after anthesis to four spring wheat cultivars had different effects in relation to the time the stress was sensed by the plants. Using an exponential model not reported previously for the adjustment of RWC measurements, two types of cultivar responses in the RWC evolution were identified and correlated with the gs evolution, the kernel/biomass ratio, and the kernel weight. In the most notorious response, at the beginning of the stress period, the cultivars LE2249 and LE2331 would have sensed the deficit as a mild water stress. This, because an early and rapid leaf wilting related to a concomitant stomatal adjustment, finally caused grain yield losses due to the decrease of the kernel weight. Cultivars with that type of response could be limited not only by the sink strength, but also by the plant source capacity during the grain-filling period.