Traits Explaining Durum Wheat (Triticum turgidum L. spp. Durum) Yield in Dry Chilean Mediterranean Environments

Yield under water stress (YS) is used as the main criterion in the selection of wheat varieties for dry Mediterranean environments. It has been proposed that selection of genotypes using YS assisted by morphological and physiological traits associated with YS is more efficient in selecting high yielding genotypes for dry environments. A study was carried out at the Antumapu Experiment Station of the University of Chile, located in Santiago, Chile (33° 40′S and 70° 38′ W). The objective was to evaluate the extent to which morpho physiological traits could explain YS. For this purpose, grain yield and yield components of 185 durum wheat genotypes from ICARDA (International Center for Agricultural Research in the Dry Areas) and INIA (Chilean National Institute for Agricultural Research) were evaluated along with seed size and weight, days to heading (DH), glaucousness (GLAU), plant height (PH) and ¹³C discrimination (Δ). The design was an α-lattice with two replications, the genotypes were grown in two different water conditions (high and low irrigation) during two seasons (2011-2012/2012-2013). Grain weight (GW) was the only yield component with high H associated with YS, but it was not associated with yield under high irrigation (YI). The combination of YI with DH+GLAU+PH+Δ+GW obtained in LI environments explained a greater fraction of YS (38%) across years; these traits had lower genotype x environment interaction than YS, they also explained a higher proportion of yield under drought than YI. None of the traits studied could replace YS in selections for grain yield. It is concluded that these traits could aid in the selection of durum wheat subject to water stress, particularly in early generations.

Introduction

Obtaining high yielding genotypes of Durum wheat (Triticum turgidum L. ssp. durum) for different environments, especially Mediterranean rainfed areas (rainy, cold winters and dry, hot summers), is considered difficult due to high genotype x environment interaction (Annicchiarico, 2002; Acevedo and Silva, 2007). For this reason, several authors have proposed that genotype selection under Mediterranean rainfed conditions may be improved by selecting for traits associated with yield under water stress (YS) having higher heritability than yield (Donald, 1968; Acevedo, 1992; Quarrie et al., 1999; Reynolds, 2006; McIntyre et al., 2010). These traits increase transpiration, transpiration efficiency and harvest index, to increase YS (Passioura, 1977, 2006).

Yield components are agronomic traits that provide an opportunity to improve yield under water stress, along with morphological and physiological traits. The yield components which together cover all the period of crop development are the number of grains per m² (GM²) and grain weight (GW). Most studies agree that GM² best explains yield (Fischer, 1973; Mouret et al., 1991; Slafer et al., 1996, 2003, 2014; Ribaut et al., 1997), and although some studies did not find an association of GW with yield under irrigation (Craufurd et al., 1991; Slafer et al., 1996, 2014), GW was found to be associated with yield in rainfed areas (Craufurd et al., 1991; Nachit and Ketata, 1991; Ribaut et al., 1997; Araus et al., 1998; Rharrabti et al., 2001; Royo et al., 2002; Nouri-Ganbalani et al., 2009). There is a negative association between GW and GM² (Slafer et al., 1996, 2014; Saeidi and Abdoli, 2015) indicating compensation, although it has been shown that in the case of wheat the compensation is not competitive (Slafer et al., 1996) and a positive association between these two yield components has also been reported (Nouri-Ganbalani et al., 2009; Motzo et al., 2010). Several studies reported GW as the yield component with the highest heritability in bread wheat and Durum wheat (Nachit et al., 1995; Royo et al., 2002; García del Moral et al., 2003; Rharrabti et al., 2003; Wang et al., 2009; Li et al., 2015). In contrast, GM² has shown similar heritability to yield in bread wheat and Durum wheat (Nachit and Ketata, 1991; Acevedo, 1992) thus yield improvements under water stress using GM² have had limited success, since GM² is as unpredictable as yield itself (Dolferus et al., 2011).

The morphological and physiological traits reported by the literature with high heritability and used in breeding programs are, (a) early vigor, determined mostly by the size of the seed. This trait increases soil coverage and plant water use, decreasing direct evaporation of soil water (Ritchie, 1972; Acevedo and Ceccarelli, 1989; López-Castañeda et al., 1996; Richards et al., 2001, 2011; Smith et al., 2003; Reynolds et al., 2004; Reynolds, 2006; Dodig et al., 2008; Blum, 2011; Elazab et al., 2015); (b) early flowering, which provides the ability to escape to terminal drought stress (Acevedo and Ceccarelli, 1989; Quarrie et al., 1999; Richards et al., 2001; Reynolds, 2006; Reynolds et al., 2007; Araus et al., 2008; Motzo et al., 2010; Blum, 2011; Crossa et al., 2016); (c) leaf glaucousness, produced by a white wax as a result of the β-diketone that appears on the plant cuticle and reduces water loss by increasing resistance to the movement of water to the leaf surface and altering the leaf energy balance favorably due to a greater reflection of short wave radiation, decreasing leaf temperature (Quarrie et al., 1999; Richards et al., 2001; Reynolds et al., 2004, 2007; Reynolds, 2006; Saint Pierre et al., 2010; Blum, 2011; Yeats and Rose, 2013; Zhang et al., 2013; Gosney et al., 2016); (d) plant height, modern genotypes are semi-dwarfs since they carry one of the Rht genes (Mellado, 2007; Sip et al., 2009), and this has been positive for grain yield, however it has been found that greater plant height has a positive effect on yield in dry Mediterranean rainfed conditions; for this reason it is desirable that the genotypes maintain their height in arid rainfed areas (Zapata et al., 2004; Royo et al., 2014) compared to irrigated areas; e) ¹³C isotopic discrimination (Δ), C₃ plants discriminate ¹³C, preferring ¹²C during CO₂ fixation (Sarievaa et al., 2010). In wheat, Δ has a positive relation to CO₂ intercellular concentration and a negative relation to transpiration efficiency (TE) (Acevedo et al., 1997; Araus et al., 2001, 2013; Condon et al., 2002; Rebetzke et al., 2002; Monneveux et al., 2005; Khazaei et al., 2009; Barbour et al., 2010; Sarievaa et al., 2010; Blum, 2011; Elazab et al., 2015). Although a negative association between Δ and yield has been found under very dry Mediterranean rainfed conditions, a positive association is more common because genotypes capable of sustaining greater stomatal conductance and water consumption, are more productive and better adapted (Araus et al., 2013; Elazab et al., 2015).

Yield under water stress (YS) has low heritability, and an important unexplained residual remains when classical selection criteria including yield potential, YS and flowering date are used (Bidinger et al., 1987) to select for dry environments. We need to know the proportion of YS explained by each trait, i.e., the ability of a morphological or physiological trait or combination to explain yield variation (Ceccarelli et al., 1991; Acevedo, 1992; Falconer and Mackay, 1996; Reynolds et al., 2007).

This study tests the hypothesis that Durum wheat yield growing under water stress can be improved when grain yield along with non-competitive yield components and morphological and physiological traits are used in selection. The overall objective is to determine how much YS of Durum wheat improves under Mediterranean arid conditions if grain weight along with the morpho-physiological traits seed size, flowering date, leaf glaucousness, plant height and ¹³C isotopic discrimination are present.

Materials and Methods

Experimental Design

The study was carried out in a dry environment, at the Antumapu Experiment Station of the Faculty of Agronomic Sciences of the University of Chile (33° S, 70° W, 625 m elevation). The soil is a sandy clay, alluvial Mollisol with a slope lower than 1%, good drainage and an effective depth between 30 and 73 cm. 185 unrelated high-yield Durum wheat genotypes (Annex I in Supplementary Material) of different origins (ICARDA and INIA) adapted to Mediterranean conditions were sown in a 14 × 14 α-lattice design with two replications for two seasons (2011-2012/2012-2013), under 2 different water conditions (high irrigation and low irrigation). The experimental unit (EU) was a plot of 2 m² (1 m wide × 2 m long) having five 2 m long and five 20 cm apart rows.

Agronomic Management

Planting was done on May 31 of each year; about 382 seeds per m² were sown manually at a depth of 3 cm. Weeds were controlled manually during the two seasons; there were no crop pests or diseases to be controlled in either season. Fertilization each season followed soil analysis and yield expectations for the site (yield target 5,000 kg/ha and 12% protein). To determine the timing and amount of irrigation, soil water content was evaluated twice per week with FDR probes (Frequency Domain Reflectometry Diviner 2000 probe), inserting two 60 cm depth access tubes in 2 plots per replication of each water condition. Climatic data were obtained from the nearby (500 m) weather station at La Platina, INIA (33° S, 70° W, and 625 m elevation) next to the test area.

Accumulated rainfall for the first season was 173 mm, which occurred between June and September; there was no rain in the sowing month (May), thus 15 mm of water were applied to facilitate sowing. In September rainfall was very low (5 mm), so irrigation (sprinkler) started at that time. In the second season rainfall was 230 mm, accumulated between May and November; the months of July, September and November had very little rain, thus irrigation began in August. The low irrigation trial was partially irrigated to avoid crop failure. The total water received by the low irrigation and high irrigation plots in the first year was 263 and 413 mm, and in the second year 275 and 440 mm, respectively (Table 1).

TABLE 1

Table 1. Rainfall and Water applied (mm) monthly for the two irrigation treatments (environments).

Evaluations

Yield, its components and morphological and physiological traits were evaluated in each genotype and each environment using the following criteria.

Yield and Its Components

Grain yield (kg ha⁻¹) was determined by the weight of grains obtained from mechanical threshing spikes harvested manually in 3 rows of 1.6 m length (0.96 m²) in the central area of each experimental unit (EU) to avoid edge effect. The yield components determined were grain weight (GW), biomass (BIO, above ground biomass), harvest index (HI), grains per spike (GS), spikes per m² (SM²), grains per m² (GM²), grain production rate (GPR, grain yield divided by the number of days from heading to physiological maturity) and biomass production rate (BPR, total above ground biomass divided by the number of days from sowing to physiological maturity). It was assumed that the small area harvested would increase grain yield similarly in all genotypes.

Morphological and Physiological Traits

Seed size was evaluated by measuring the weight (SW) and length of 10 seeds of each genotype (López-Castañeda et al., 1996). Grain weight (GW) was measured at harvest.

Flowering date

Anthesis in Durum wheat is difficult to recognize, since the anthers of some genotypes are not visible; for this reason we measured the number of days from sowing to heading, determining this state when the 75% of the EU had the spikes completely emerged (Zadoks 59) (Zadoks et al., 1974). The number of days between sowing and heading corresponds to days to heading (DH).

Physiological maturity

The date of physiological maturity was determined when 75% of the EU showed a yellow color of the flag leaf, the grains had a dry appearance and pinching them did not leave a mark (Zadoks 92), the number of days between sowing date and physiological maturity date (DM) corresponds to days to physiological maturity (Zadoks et al., 1974).

Glaucousness (GLAU) was recorded visually on a 1 to 5 scale (notes 1 and 5 for the least and most glaucous respectively) at Zadoks 59 (Zadoks et al., 1974), when it was at its maximum expression in leaves, stems, and spikes (Richards et al., 2001).

Plant height (PH) was measured in 5 plants of each EU at the end of the growing period from the ground to the apex of the main stem, without considering the awns (Zapata et al., 2004; Castañeda-Saucedo et al., 2009).

¹³C discrimination

Carbon discrimination analysis was performed on samples of 52 randomly selected genotypes; 10 g of grain at harvest were separated and dried in a stove at 70°C for 48 h and then encapsulated in tin capsules for ¹³C composition isotope analysis. The samples were sent to the Stable Isotope Facility of the University of California at Davis, U.S.A. to be analyzed using a PDZ Europe ANCA-GSL elemental analyzer in interface with an isotope ratio mass spectrometer (IRMS) PDZ Europe 20-20 (Sercon Ltd., Cheshire, UK). Values of δ¹³C (product) obtained were used to calculate ¹³C isotopic discrimination (Δ) as follows:

$\begin{array}{ll}\Delta =\frac{{\delta }^{13}C(source)-{\delta }^{13}C(product)}{1+{\delta }^{13}C(source)÷1000}& (1)\end{array}$

Atmospheric δ¹³C is usually between 6.5 and 8.0‰, according to geographical location and altitude (Acevedo et al., 1997; Monneveux et al., 2006). A 7‰ was used in this case.

The proportion of grain yield explained by a given trait in the presence of all possible combinations of the other traits was calculated using the following formula (Ceccarelli et al., 1991; Acevedo, 1992; Reynolds et al., 2007).

$\begin{array}{l}\%Yieldexplainedbytrait=\\ \frac{\left[Yieldwithtrait(+)-Yieldwithtrait(-)\right]}{(10\%Higheryield-10\%Loweryield)}& (2)\end{array}$

10% Higher yield refers to the highest 10% yield of the nursery and 10% Lower yield to the lowest 10% yield. Yield with trait (+) is the mean yield of genotypes having a positive expression of a trait or combination of traits. Yield with trait (−) is the mean yield of the nursery having a negative expression of a trait or combination of traits.

The data used in the calculation of the proportion of grain yield explained by a given trait were obtained through selections made by a selection index software of CIMMYT (1999) for the different combinations of the traits. To make the selection it is necessary to specify the target and intensity for each trait included in the combination, where the target is the objective which is pursued in the selection, expressed in units of standard deviation of the mean (from −3.0 to 3.0); the intensity reflects the relative importance of the different variables to be used in selection expressed in values from 0 to 10; the trait that you want to progress more must receive higher intensity (CIMMYT, 1999).

Statistical Analysis

For statistical analysis we used the InfoStat statistical software (Di Rienzo et al., 2008), which has a friendly interface of the R platform to estimate the General Linear Mixed Model, which allows differentiating random effects from fixed effects (Di Rienzo et al., 2011). Genotypes, years, water conditions and the interaction between them were introduced in the model as fixed or random effects depending on the objective of the analysis. When the objective was the comparison of the tested material for selection, the fixed model was used. To control the experimental error from random effects the cross effects of row and column of the α-lattice design were introduced in the model as random effects nested within the blocks, which in turn were nested within the year-water condition effect. These models best fit the characteristics of this experiment. The presence or absence of significant variance was determined (p ≤ 0.05) in terms of traits, genotypes, environments, and interactions. For ¹³C isotope discrimination, we used 52 out of the 185 genotypes what allowed the original lattice design to be maintained. Therefore the analysis of variance for ¹³C isotope discrimination was carried out similarly to the other traits but with a lower number of genotypes (n = 52). The percentage of the total of the sum of squares corresponding to each source of variation in each variable was calculated. The variances needed to calculate broad sense heritability (H) of each trait were calculated using the random model and the following relation:

$\begin{array}{ll}H=\frac{{\sigma }_g^2}{{\sigma }_g^2+\frac{{\sigma }_{ge}^2}{e}+\frac{{\sigma }^2}{re}}& (3)\end{array}$

where r = number of replications, e = number of environments, σ² = error variance, ${\sigma }_g^2=$ genetic variance, ${\sigma }_{ge}^2$ = variance of genotype X environment interaction (Falconer and Mackay, 1996).

Pearson correlations were calculated using the statistical software InfoStat (Di Rienzo et al., 2008) to determine the existence and magnitude of associations (p ≤ 0.05) between yield and traits; the software was also used to perform multiple linear regressions to estimate yield, using the traits and hence to determine how much the residuals changed with the presence or absence of a trait.

The worth of a trait for plant improvement was tested with the following relation (Falconer and Mackay, 1996).

$\begin{array}{ll}Cry/Dry=r_{\text{G}}sqr{H}_{\text{c}}/H_{\text{s}}& (4)\end{array}$

Where r_G is the genotypic correlation between the indirect trait (Cry) and the direct trait (Dry), H_c represents H of the indirect trait and H_s represents H of the direct trait.

Results

Associations between Traits and Yield

Table 2 shows the Pearson correlations with grain yield of the variables of concern for all environments. All variables were associated with yield in at least one year under the low irrigation treatment, but not all the variables were associated with yield under high irrigation. In all environments yield components had higher association with yield than the morphological and physiological traits, except for grain weight, which although associated with yield in both years under low irrigation, was not associated with yield in either year when the crop was fully irrigated. Δ was the only physiological trait associated with yield under high irrigation.

TABLE 2

Table 2. Mean of yield, its components and morphological and physiological traits in each environment, Pearson correlation (r) of the morpho-physiological traits with yield in each environment; percent of genotypic variability (G); genotype x environment interaction (GxE) of the total sum of squares of the analysis of variance for each variable calculated with the fixed linear mixed model.

We analyzed the correlation between years of yield and the variables under low irrigation. The results are shown in Table 3, where it can be seen that traits were correlated between years under low irrigation except for yield (p ≤ 0.10).

TABLE 3

Table 3. Pearson correlation coefficients among years of yield and of high heritability variables under low irrigation.

Seed length was measured in the first year only and was correlated (r = 0.60) with seed weight (SW) in that year, while SW was correlated with grain weight (GW) in both environments with a Pearson correlation coefficient above 0.33; in addition, a negative association was found between GM² and GW in all environments (Table 4).

TABLE 4

Table 4. Pearson correlation coefficients among grain variables in each environment.

Trait Contributions to Yield under Water Stress

To assess how much each trait contributed to yield under low irrigation, multiple linear regressions of the traits including YI, with yield under low irrigation as the dependent variable were performed for each year. The residuals of the multiple regressions corresponded to the yield variation not explained by the regression; the best regression been the one with the smallest residual. The analysis was performed with the 52 genotypes in which Δ was measured, considering the traits listed in Table 5.

TABLE 5

Table 5. Percentages of the residuals of the sum of squares of the multiple linear regression along with the coefficient of determination.

Using classical selection criteria (YS = DH+YI) the minimum residual was 81%, significant only for one of the 2 years under low irrigation. The residuals were reduced adding traits independently, but the reduction was significant only in one of the years, except for SW and Δ, whose reduction was significant in both years under low irrigation. The combination of two or more traits reduced the residuals to less than 74%, reaching 51% in the first year when all traits were combined (Table 5). The traits along with DH and YI that most decreased the residuals were GLAU, PH, Δ and GW, with no differences between years (Table 5).

Percentage of YS Explained by Each Trait or Combination

Table 6 shows the proportion of YS explained by each trait separately in the year of greater water stress, the target of selection was + 3 and −3, and the intensity was 10 for each of these traits. The traits with the lowest explanation of YS were SW and YI (0.06 and 0.01, respectively), and those with the higher explanation of YS were Δ and GW, followed by PH and GLAU.

TABLE 6

Table 6. Explanation of YS by trait selecting 10% of the total genotypes having maximum and minimum expression of each character, and the 10% of the genotypes of maximum and minimum yield in the year of greater water stress (Equation 2).

Table 7 shows the proportion of YS explained by combinations of the traits in the year with greater water stress. All traits were considered, including Δ, so the corresponding population was 52 genotypes, target was + 3 and −3, and intensity 10. Almost all the combinations shown in this table explain a higher proportion of YS than any trait individually (0.37 in Table 6); except for those in which YI is added to the combination (0.35) and Δ combined with DH (0.31). Of the 3 greatest explanations of YS, the combination Δ+DH+GLAU+PH has the greatest number of traits.

TABLE 7

Table 7. Explanation of YS by combinations of traits selecting 10% of 52 genotypes having maximum and minimum expression of each combinations obtained under LI environments, and the 10% of the genotypes of maximum and minimum yield in the year of greater water stress.

The genotypes of Durum wheat showing the highest expression of each of the traits of the combination Δ+DH+GLAU+PH were in the environment of greater water stress. Note that lower values of DH are considered greater expression. These genotypes showed yields in the Mediterranean low irrigation conditions above the mean (2,502 kg/ha) of the original population (185 genotypes). The number of each genotype is given in Table 8 and its description in the Annex I; these genotypes are ordered within columns from greater to lower trait expression. Genotypes marked with * are among the 5 genotypes of greatest YS out of a total of 185.

TABLE 8

Table 8. Durum wheat genotypes selected by the combination of the maximum expression of the traits ¹³C discrimination, glaucousness, plant height and days to heading, in the environment of greater water stress (n = 185, except for the genotypes with ¹³C discrimination determination, n = 52).

Implications for Plant Breeding

Traits useful for breeding plants for water stressed environments are required to have a high genetic correlation with YS along with high heritability. Formula 4 from Falconer and Mackay (1996) puts these two concepts together. Table 9 shows the ratio of indirect (selection for a trait other than grain yield) vs direct selection (selection for yield under water stress) and shows that selection for SW, DH, GLAU, PH and Δ, all make a positive contribution to selection for grain YS (recall that lower days to heading is the desirable trait for DH in Mediterranean environments). The positive contribution of these traits occurs under water stress only. The last column of Table 9 shows that the contribution drops to practically zero or is even negative under high irrigation. The heritabilities in this table were calculated with the random linear mixed model.

TABLE 9

Table 9. Genetic correlation (r_g) of each trait with yield in low irrigation environments.

None of the indirect traits by itself will lead to higher yield (Table 9) and they may be only used in early stages of selection and only when trials are conducted under water stress conditions.

Discussion

A considerable decrease in durum wheat yield occurs with lower than 300 mm annual rainfall which is associated with greater genotype x environment interaction (Annicchiarico, 2002). When breeding for high yielding genotypes in these environments, the use of indirect selection traits, associated with greater grain yield under water stress, having lower genotype × environment interaction would make results more reliable and repeatable in most of the environments (Falconer and Mackay, 1996).

The criteria for choosing a trait of indirect selection included genetic variability, easy to measure, association with the direct trait (grain yield), high heritability and low genotype x environment interaction (Falconer and Mackay, 1996; Reynolds et al., 2007). Traits which met these criteria were grain weight, days to heading, glaucousness, plant height and Δ (Tables 2, 9).

All agronomic traits, except SM² had high merit for indirect selection under the dry environment (Table 9). GW is not competitive with other traits, making wheat improvement using other agronomic traits relatively useless. The highest heritabilities (Table 9) were obtained for the morpho physiological traits, except for Δ which had a very low H (0.19).

Although the number of grains (GM²) has been traditionally used to improve wheat yield under irrigation (Slafer et al., 1996, 2014), and shows a high correlated response with grain yield under stress (Crx/Cry, YS Table 9), it is extremely vulnerable to stress, which causes abortion and floret sterility, affecting the number of grains (Dolferus et al., 2011). Because of this, yield improvement by selection using GM² has had limited success in stress environments (Dolferus et al., 2011).

Another important point to consider in the choice of traits to combine guiding the selection process is their association with YS (Table 9). All yield components had a significant genotypic correlation with yield under low irrigation but they lose their value due to the compensatory interactions with other yield components. There was no association of GW with yield under high irrigation (data not shown), in accordance with other studies (Slafer et al., 1996, 2014), instead grain weight was associated with yield under low irrigation (Tables 2, 9) which has also been reported in other studies (Reynolds, 2006). Days to heading was negatively correlated with yield, indicating that the most precocious genotypes would be desirable, in accordance with other reports for Mediterranean environments (Acevedo and Ceccarelli, 1989; Quarrie et al., 1999; Richards et al., 2001; Reynolds, 2006).

According to Table 2 some traits were associated with yield under low irrigation in one year but not in the other, what could result from the inter-annual variation of yield and not of the traits. Table 3 shows that yield under low irrigation was different between years; by contrast all indirect traits correlated with LI yield in one year were also correlated between years under LI (Table 3), demonstrating that this effect was mainly due to the inter-annual variation of YS and no to the variation of the traits. Furthermore, the merit of the indirect selection traits holds under water stress only (last two columns, Table 9).

As reported by other authors (López-Castañeda et al., 1996; Smith et al., 2003) we found that seed weight was associated with seed length (Table 4). The correlation found between seed weight (SW) and grain weight (GW) indicates that sowing large seeds will result in large grains at harvest (Richards et al., 2001; Reynolds et al., 2004; Reynolds, 2006; Sadras, 2007; Table 4). If genotypes of greater GW are chosen in the first generation, in the second generation there will be greater GW under water stress, since GW has high heritability (Table 9). GW would be a good estimator of grain weight at sowing (SW) and YS (Table 2).

There was a negative relation between GW and GM² in all environments of this study (Table 4), which indicates compensation between GW and GM², thus trying to improve YS using GM² would be counterproductive because individual grain weight will be reduced. It has been argued that this compensation is not competitive (Slafer et al., 1996), i.e., there is no competition between grains for available assimilates, since these components are determined at different growth stages of the crop; GM² is determined before GW begins and there is no feedback, thus theoretically variation in GM² should have a minor effect on GW and grain size could be increased while maintaining the number of grains under stress conditions (Motzo et al., 2010). This was not the case in this study.

Indirect traits selected that complied with the criteria to be used in plant breeding for YS included days to heading (DH), glaucousness (GLAU), plant height (PH), ¹³C Discrimination (Δ), grain weight (GW) and seed weight (SW).

In plant breeding programs for dry Mediterranean environments the trait most used is flowering date in conjunction with yield potential, but while it is true that flowering date has high merit for selection and is negatively associated with yield in Mediterranean environments (DH, Table 9) when improving using these criteria much unexplained residual remained, thus it was important to know what would happen with the residual if one or all traits were added to a multiple linear regression (Table 5). The indirect selection traits for YS chosen in this study, GLAU, PH, Δ, GW, and SW, together with DH determined under drought stress and yield under full irrigation (YI) decreased the residuals (Table 5). When these traits were combined, the greatest decrease in residuals was the one that contained all the traits except seed weight (SW) (Table 5). When GW was not introduced in the combination, the residuals increased, as well as the residuals difference between years, showing the importance of grain weight in the combination. Thus the traits ¹³C discrimination (Δ), plant height (PH), glaucousness (GLAU), and grain weight (GW) combined with yield under irrigation (YI) and flowering date (DH) were those which explained a greater fraction of yield under low irrigation (Table 5).

The measure of the percentage of the overall yield explained by a trait is an index that gives a measure of how good the trait is to explain yield indirectly (Ceccarelli et al., 1991). The objective of using traits to select yield under stress indirectly is to achieve an increase in yield heritability and a decrease in the genotype x environment interaction, thus understanding what trait or combination of traits explains a higher percentage of YS allows choosing those which most contribute to yield and will lead to better selection with greater contribution to heritability (Table 9).

By using selection criteria for each single trait we achieved explanation of YS greater than 0.23 in almost all the traits except for seed weight (0.06) and yield under HI (0.01) (Table 6). When YI was included in trait combinations the percent of explanation of YS decreased (Table 7) despite being correlated with yield under low irrigation, indicating that it is not a good indirect selection trait for genotypes under water stress. While the morpho-physiological traits together without considering YI (Table 7) always showed a higher explanation of YS than separately (Table 6), the use of more traits did not significantly improve the explanation of YS (Table 7). This is not to say that more traits should not be used, because the contribution of the traits in terms of heritability (Table 9) should also be considered. The greatest explanation of YS (0.59) was obtained by combining Δ and glaucousness (GLAU), and it was followed by the combination that also contained DH (0.56) (Table 7). Although these appear to be the best combinations, analyzing trait heritabilities (Table 9), it is clear that the gain in explanation of YS does not imply better heritability, therefore it would be preferable to use the combination that does not significantly decrease explanation of YS and employs the largest number of traits with high heritability, as in the case of the combination Δ + DH + GLAU + PH (Table 7). While it is true that using a large number of traits increases economic cost, to obtain the best estimate one must choose between the more effective method (combined selection) and the less expensive method (single selection) (Falconer and Mackay, 1996).

In crop improvement by backcrossing the final product is a line containing the gene of interest from the donor parent, with the genotype of the recurrent parent in the rest of the genome (Hospital, 2005). Therefore, a good donor genotype for a trait is one that provides its best expression, because the recurrent parent provides the rest of the genome. This study allowed us to identify genotypes of Durum wheat as donors of each trait of the best proposed combination (Δ+DH+GLAU+PH) for higher YS (Table 8). Some of these genotypes are also the genotypes with greatest yield in the low irrigation conditions (genotypes 40, 53, and 137). To understand in detail the features of each genotype selected, the Annex I provides a description of each of the 185 genotypes used in this study.

Conclusions

The selection for yield under drought of Durum wheat under Mediterranean arid conditions, using grain weight along with the morpho-physiological traits seed size, flowering date, leaf glaucousness, plant height and ¹³C isotopic discrimination can assist in selection for grain yield particularly in early generations but the nursery must be grown under drought stress. Possible donor parents of each trait of the best proposed combination (Δ+DH+GLAU+PH) for higher YS were identified.

None of the indirect selection traits tested here will lead per se to higher grain yield. They may assist in the selection process, particularly in early generations.

The hypothesis that selection of Durum wheat YS improves when grain yield along with non-competitive yield components and morphological and physiological traits are used is sustained.

Author Contributions

GG and EA conceived the experiment, GG, PS, and MO conducted the Experiment, GG analyzed the results. GG and EA wrote the manuscript. All authors reviewed the manuscript.

Funding

During this study first author was funded by The National Commission for Scientific and Technological Research of Chile (CONICYT) scholarship N° 63100017 and The Soil-Plant-Water Relations Laboratory of the University of Chile. The genotypes of durum wheat used were provided by Dr. Ivan Matus from The National Agricultural Research Institute of Chile (INIA); trials carried out in field were funded by the Laboratory of Soil-Water-Plant Relationships (SAP) of the University of Chile.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors thank Mr. Marco Garrido for critical reading of the manuscript and useful suggestions.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2017.01781/full#supplementary-material

References

Acevedo, E. (1992). Increasing the Yield Potential of Irrigated Bread Wheat: Basis for Physiological Research at CIMMYT (No. 12). México, D.F: CIMMYT.

Acevedo, E., Baginsky, C., Solar, B., and Ceccarelli, S. (1997). Discriminación Isotópica de ¹³C y su Relación con el Rendimiento y la Eficiencia de Transpiración de Genotipos Locales y Mejorados de Cebada Bajo Diferentes Condiciones Hídricas. Investigación Agrícola Chile 17, 41–54.

Google Scholar

Acevedo, E., and Ceccarelli, S. (1989). “Role the physiologist-breeder in a breeding program for drought resistance conditions,” in Drought Resitance in Cereal, ed F. W. G. Baker (Wallingford, UK: C.A.B International), 117–139.

Google Scholar

Acevedo, E., and Silva, P. (Eds.). (2007). Trigo Candeal: Calidad, Mercado y Zonas de Cultivo. Universidad de Chile, Facultad de Ciencias Agronómicas.

Annicchiarico, P. (2002). Genotype x Environment Interactions: Challenges and Opportunities for Plant Breeding and Cultivar Recommendations (No. 174). Rome: Food & Agriculture Org.

Google Scholar

Araus, J. L., Amaro, T., Casadesús, J., Asbati, A., and Nachit, M. M. (1998). Relationships between ash content, carbon isotope discrimination and yield in durum wheat. Aust. J. Plant Physiol. 25, 835–842. doi: 10.1071/PP98071

CrossRef Full Text | Google Scholar

Araus, J. L., Cabrera-Bosquet, L., Serret, M. D., Bort, J., and Nieto-Taladriz, M. T. (2013). Comparative performance of δ13C, δ18O and δ15N for phenotyping durum wheat adaptation to a dryland environment. Funct. Plant Biol. 40, 595–608. doi: 10.1071/FP12254

CrossRef Full Text | Google Scholar

Araus, J. L., Casadesus, J., and Bort, J. (2001). “Recent tools for the screening of physiological traits determining yield,” in Application of Physiology in Wheat Breeding, eds M. Reynolds, J. Ortiz-Monasterio and A. McNab (México, D.F: CIMMYT), 59–77.

Google Scholar

Araus, J. L., Slafer, G. A., Royo, C., and Serret, M. D. (2008). Breeding for yield potential and stress adaptation in cereals. CRC Crit. Rev. Plant Sci. 27, 377–412. doi: 10.1080/07352680802467736

CrossRef Full Text | Google Scholar

Barbour, M. M., Warren, C. R., Farquhar, G. D., Forrester, G. U. Y., and Brown, H. (2010). Variability in mesophyll conductance between barley genotypes, and effects on transpiration efficiency and carbon isotope discrimination. Plant Cell Environ. 33, 1176–1185. doi: 10.1111/j.1365-3040.2010.02138.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Bidinger, F. R., Mahalakshmi, V., and Rao, G. D. P. (1987). Assessment of drought resistance in pearl millet [Pennisetum americanum (L.) Leeke]. I. Factors affecting yields under stress. Aust. J. Agric. Res. 38, 37–48. doi: 10.1071/AR9870037

CrossRef Full Text | Google Scholar

Blum, A. (2011). Plant Breeding for Water-Limited Environments, 254. New York, NY: Springer Science+Business Media, LLC.

Google Scholar

Castañeda-Saucedo, M. C., López-Casta-eda, C., Colinas-De León, M. T. B., Molina, J. C., and Hernández-Livera, A. (2009). Rendimiento y calidad de la semilla de cebada y trigo en campo e invernadero. Interciencia 24, 286–292.

Ceccarelli, S., Acevedo, E., and Grando, S. (1991). Breeding for yield stability in unpredictable environments: single traits, interaction between traits, and architecture of genotypes. Euphytica 56, 169–185. doi: 10.1007/BF00042061

CrossRef Full Text | Google Scholar

CIMMYT (1999). Manual de Usuario Para Fieldbook 5.117.1 y Alfa. México, D.F: CIMMYT.

Condon, A. G., Richards, R. A., Rebetzke, G. J., and Farquhar, G. D. (2002). Improving intrinsic water-use efficiency and crop yield. Crop Sci. 42, 122–131. doi: 10.2135/cropsci2002.1220

PubMed Abstract | CrossRef Full Text | Google Scholar

Craufurd, P. O., Clipson, N. J., Austin, R. B., and Acevedo, E. (1991). “Defining an ideotype for barley in low rainfall Mediterranean environments,” in Improvement and Management of Winter Cereals Under Temperature, Drought and Salinity Stresses: Proceedings of the International Symposium, Vol. Monografías INIA, ed E. Acevedo (Córdba: Instituto Nacional de Investigaciones Agrarias), 309–325.

Google Scholar

Crossa, J., de los Campos, G., Maccaferri, M., Tuberosa, R., Burgueño, J., and Pérez-Rodríguez, P. (2016). Extending the marker × environment interaction model for genomic-enabled prediction and genome-wide association analysis in durum wheat. Crop Sci. 56, 2193–2209. doi: 10.2135/cropsci2015.04.0260

CrossRef Full Text | Google Scholar

Di Rienzo, J. A., Casanoves, F., Balzarini, M. G., Gonzalez, L., Tablada, M., and Robledo, C. W. (2008). InfoStat version 2008. Córdoba: InfoStat Group, College of Agricultural Sciences, National University.

Di Rienzo, J. A., Macchiavelli, R. E., and Casanoves, F. (2011). Modelos Lineales Mixtos: Aplicaciones en InfoStat, Córdoba: Grupo Infostat.

Dodig, D., Zoric, M., Knezevic, D., King, S. R., and Surlan-Momirovic, G. (2008). Genotype environment interaction for wheat yield in different drought stress conditions and agronomic traits suitable for selection. Aust. J. Agric. Res. 59, 536–545. doi: 10.1071/AR07281

CrossRef Full Text | Google Scholar

Dolferus, R., Ji, X., and Richards, R. A. (2011). Abiotic stress and control of grain number in cereals. Plant Sci. 181, 331–341. doi: 10.1016/j.plantsci.2011.05.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Donald, C. M. (1968). The breeding of crop ideotypes. Euphytica 17, 385–403. doi: 10.1007/BF00056241

CrossRef Full Text | Google Scholar

Elazab, A., Bort, J., Zhou, B., Serret, M. D., Nieto-Taladriz, M. T., and Araus, J. L. (2015). The combined use of vegetation indices and stable isotopes to predict durum wheat grain yield under contrasting water conditions. Agric. Water Manag. 158, 196–208. doi: 10.1016/j.agwat.2015.05.003

CrossRef Full Text | Google Scholar

Falconer, D. S., and Mackay, T. F. C. (1996). Introduction to Quantitative Genetics. Harlow: Longman.

Google Scholar

Fischer, R. A. (1973). “The effect of water stress at various stages of development on yield processes in wheat,” in Plant Response to Climatic Factors: Proceedings of the Uppsala Symposium, ed R. O. Slatyer (Paris: Ecology and Conservation 5, UNESCO), 233–241.

Google Scholar

García del Moral, L. F., Rharrabti, Y., Villegas, D., and Royo, C. (2003). Evaluation of grain yield and its components in durum wheat under mediterranean conditions: an ontogenic approach. Agron. J. 95, 266–274. doi: 10.2134/agronj2003.0266

CrossRef Full Text | Google Scholar

Gosney, B. J., Potts, B. M., O'Reilly-Wapstra, J. M., Vaillancourt, R. E., Fitzgerald, H., Davies, N. W., et al. (2016). Genetic control of cuticular wax compounds in Eucalyptus globulus. New Phytol. 209, 202–215. doi: 10.1111/nph.13600

PubMed Abstract | CrossRef Full Text | Google Scholar

Hospital, F. (2005). Selection in backcross programmes. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 360, 1503–1511. doi: 10.1098/rstb.2005.1670

PubMed Abstract | CrossRef Full Text | Google Scholar

Khazaei, H., Mohammady, S.-D., Zaharieva, M., and Monneveux, P. (2009). Carbon isotope discrimination and water use efficiency in Iranian diploid, tetraploid and hexaploid wheats grown under well-watered conditions. Genet. Resour. Crop Evol. 56, 105–114. doi: 10.1007/s10722-008-9349-3

CrossRef Full Text | Google Scholar

Li, Q., Zhang, Y., Liu, T., Wang, F., Liu, K., Chen, J., et al. (2015). Genetic analysis of kernel weight and kernel size in wheat (Triticum aestivum L.) using unconditional and conditional QTL mapping. Mol. Breed. 35, 194. doi: 10.1007/s11032-015-0384-4

CrossRef Full Text | Google Scholar

López-Castañeda, C., Richards, R. A., Farquhar, G. D., and Williamson, R. E. (1996). Seed and seedling characteristics contributing to variation in early vigor among temperate cereals. Crop Sci. 36, 1257–1266. doi: 10.2135/cropsci1996.0011183X003600050031x

CrossRef Full Text | Google Scholar

McIntyre, C. L., Mathews, K. L., Rattey, A., Chapman, S. C., Drenth, J., Ghaderi, M., et al. (2010). Molecular detection of genomic regions associated with grain yield and yield-related components in an elite bread wheat cross evaluated under irrigated and rainfed conditions. Theor. Appl. Genet. 120, 527–541. doi: 10.1007/s00122-009-1173-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Mellado, M. (2007). El Trigo en Chile. Cultura, Ciencia y Tecnología. (Colección Libros INIA; 21), 684. Chillán: Centro Regional de Investigación Quilamapu.

Monneveux, P., Rekika, D., Acevedo, E., and Merah, O. (2006). Effect of drought on leaf gas exchange, carbon isotope discrimination, transpiration efficiency and productivity in field grown durum wheat genotypes. Plant Sci. 170, 867–872. doi: 10.1016/j.plantsci.2005.12.008

CrossRef Full Text | Google Scholar

Monneveux, P., Reynolds, M. P., Trethowan, R., González-Santoyo, H., Pe-a, R. J., and Zapata, F. (2005). Relationship between grain yield and carbon isotope discrimination in bread wheat under four water regimes. Eur. J. Agron. 22, 231–242. doi: 10.1016/j.eja.2004.03.001

CrossRef Full Text | Google Scholar

Motzo, R., Giunta, F., and Pruneddu, G. (2010). The response of rate and duration of grain filling to long-term selection for yield in Italian durum wheats. Crop Pasture Sci. 61, 162–169. doi: 10.1071/CP09191

CrossRef Full Text | Google Scholar

Mouret, J. C., Conesa, A. P., Bouchier, A., Gaid, M., and Sennour, M. (1991). “Identification des facteurs de variabilité du rendement du blé dur en conditions hydriques limitantes dans la région de Sidi -bel-Abbés (Algérie),” in Physiology Breeding of Winter Cereals for stressed Mediterranean Environments, eds E. Acevedo, A.P. Conesa, P. Monneveux, and P. Srivastava (Montpellier: INRA), 35–47.

Nachit, M. M., Baum, M., Impiglia, A., and Ketata, H. (1995). “Studies on some grain quality traits in durum wheat grown in mediterranean environments,” in Durum wheat quality in the Mediterranean region. Seminar on Durum Wheat Quality in the Mediterranean Region, eds N. Di Fonzo, F. Kaan, and M. Nachit (Zaragoza: CIHEAM, Options Méditerranéennes: Série A. Séminaires Méditerranéens; n. 22),181–187.

Google Scholar

Nachit, M. M., and Ketata, H. (1991). “Selection of morpho-physiological traits for multiple abiotic stresses resistance in durum wheat (Triticum turgidum L. var. Durum),” in Physiology Breeding of Winter Cereals for stressed Mediterranean Environments, eds E. Acevedo, A.P. Conesa, P. Monneveux, and Srivastava (Montpellier: INRA), 35–47.

Nouri-Ganbalani, A., Nouri-Ganbalani, G., and Hassanpanah, D. (2009). Effects of drought stress condition on the yield and yield components of advanced wheat genotypes in Ardabil, Iran. J. Food Agric. Environ. 7, 228–234.

Google Scholar

Passioura, J. B. (1977). Grain yield, harvest index and water use of wheat. J. Aust. Inst. Agric. Sci. 43, 117–120.

Google Scholar

Passioura, J. B. (2006). Increasing crop productivity when water is scarce—from breeding to field management. Agric. Water Manag. 80, 176–196. doi: 10.1016/j.agwat.2005.07.012

CrossRef Full Text | Google Scholar

Quarrie, S. A., Stojanović, J., and Pekić, S. (1999). Improving drought resistance in small-grained cereals: a case study, progress and prospects. Plant Growth Regul. 29, 1–21. doi: 10.1023/A:1006210722659

CrossRef Full Text | Google Scholar

Rebetzke, G. J., Condon, A. G., Richards, R. A., and Farquhar, G. D. (2002). Selection for reduced carbon isotope discrimination increases aerial biomass and grain yield of rainfed bread wheat. Crop Sci. 42, 739–745. doi: 10.2135/cropsci2002.0739

CrossRef Full Text

Reynolds, M. (2006). “Trigo que utiliza el agua de manera eficiente y produce grano de buena calidad,” in Informe Annual del CIMMYT 2006-2007 (Mexico, D.F: CIMMYT), 14–15.

Google Scholar

Reynolds, M. P., Pierre, C. S., Saad, A. S. I., Vargas, M., and Condon, A. G. (2007). Evaluating potential genetic gains in wheat associated with stress-adaptive trait expression in elite genetic resources under drought and heat stress. Crop Sci. 47, S-172–S-189. doi: 10.2135/cropsci2007.10.0022IPBS

CrossRef Full Text | Google Scholar

Reynolds, M. P., Rubeena, S., and Trethowan, R. (2004). “Using “smart” physiological-trait based crossing strategies to accumulate drought-adaptive genes,” in Resilient Crops for Water Limited Environments: Proceedings of a Workshop Held at Cuernavaca, eds P. Poland, M. Sawkins, J. Ribaut, and D. Hoisington (Mexico D.F: CIMMYT), 185–187.

Rharrabti, Y., Garcia del Moral, L. F., Villegas, D., and Royo, C. (2003). Durum wheat quality in mediterranean environments: III. Stability and comparative methods in analysing G × E interaction. Field Crops Res. 80, 141–146. doi: 10.1016/S0378-4290(02)00178-8

CrossRef Full Text | Google Scholar

Rharrabti, Y., Villegas, D., García del Moral, L. F., Aparicio, N., Elhani, S., and Royo, C. (2001). Environmental and genetic determination of protein content and grain yield in durum wheat under Mediterranean conditions. Plant Breed. 120, 381–388. doi: 10.1046/j.1439-0523.2001.00628.x

CrossRef Full Text | Google Scholar

Ribaut, J.-M., Jiang, C., Gonzalez-de-Leon, D., Edmeades, G. O., and Hoisington, D. A. (1997). Identification of quantitative trait loci under drought conditions in tropical maize. 2. Yield components and marker-assisted selection strategies. Theor. Appl. Genet. 94, 887–896. doi: 10.1007/s001220050492

CrossRef Full Text

Richards, R. A., Condon, A. G., and Rebetzke, G. J. (2001). Traits to Improve Yield in Dry Environments in Application of Physiology in Wheat Breeding, eds M. P. Reynolds, J. I. Ortiz-Monasterio and A. McNab (Mexico, D.F: CIMMYT), 88–100.

Richards, R. A., Rebetzke, G. J., Condon, A. G., and Watt, M. (2011). “Breeding to improve grain yield in water limited environments: the CSIRO experience with wheat,” in Crop Stress Management and Global Climate Change, eds J. L. Araus and G. A. Slafer (Wallington, UK: CABI), 105–121.

Google Scholar

Ritchie, J. T. (1972). Model for predicting evaporation from a row crop with incomplete cover. Water Resour. Res. 8, 1204–1213. doi: 10.1029/WR008i005p01204

CrossRef Full Text | Google Scholar

Royo, C., Nazco, R., and Villegas, D. (2014). The climate of the zone of origin of Mediterranean durum wheat (Triticum durum Desf.) landraces affects their agronomic performance. Genet. Resour. Crop Evol. 61, 1345–1358. doi: 10.1007/s10722-014-0116-3

CrossRef Full Text | Google Scholar

Royo, C., Villegas, D., Del Moral, L. G., Elhani, S., Aparicio, N., Rharrabti, Y., et al. (2002). Comparative performance of carbon isotope discrimination and canopy temperature depression as predictors of genotype differences in durum wheat yield in Spain. Aust. J. Agric. Res. 53, 561–569. doi: 10.1071/AR01016

CrossRef Full Text | Google Scholar

Sadras, V. O. (2007). Evolutionary aspects of the trade-off between seed size and number in crops. Field Crops Res. 100, 125–138. doi: 10.1016/j.fcr.2006.07.004

CrossRef Full Text | Google Scholar

Saeidi, M., and Abdoli, M. (2015). Effect of Drought Stress during grain filling on yield and its components, gas exchange variables, and some physiological traits of wheat cultivars. J. Agric. Sci. Technol. 17, 885–898.

Google Scholar

Saint Pierre, C., Crossa, J., Manes, Y., and Reynolds, M. P. (2010). Gene action of canopy temperature in bread wheat under diverse environments. Theor. Appl. Genet. 120, 1107–1117. doi: 10.1007/s00122-009-1238-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Sarievaa, G. E., Kenzhebaevab, S. S., and Lichtenthalerc, H. K. (2010). Adaptation potential of photosynthesis in wheat cultivars with a capability of leaf rolling under high temperature conditions. Rus. J. Plant Physiol. 57, 28–36. doi: 10.1134/S1021443710010048

CrossRef Full Text | Google Scholar

Sip, V., RuŽek, P., Chrpová, J., Vavera, R., and Kusá, H. (2009). The effect of tillage practice, input level and environment on the grain yield of winter wheat in the Czech Republic. Field Crops Res. 113, 131–137. doi: 10.1016/j.fcr.2009.04.013

CrossRef Full Text | Google Scholar

Slafer, G. A., Calderini, D. F., and Miralles, D. J. (1996). “Yield components and compensation in wheat: opportunities for further increasing yield potential,” in Increasing Yield Potential in Wheat: Breaking the Barriers: Proceedings of a Workshop Held in Ciudad Obregón, (Sonora: CIMMYT). 101–133.

Google Scholar

Slafer, G. A., Miralles, D. J., Savin, R., Whitechurch, E. M., and González, F. G. (2003). “Ciclo Ontogénico, Dinámica del Desarrollo y Generación del Rendimiento y la Calidad en Trigo,” in Producción de Granos: Bases Funcionales para su Manejo, eds H. Satorre, R. L. Benech Arnold, G. A. Slafer, E. B. de la Fuente and D. J. Miralles (Buenos Aires: Editorial Facultad de Agronomía, Universidad de Buenos Aires), 101–132.

Slafer, G. A., Savin, R., and Sadras, V. O. (2014). Coarse and fine regulation of wheat yield components in response to genotype and environment. Field Crops Res. 157, 71–83. doi: 10.1016/j.fcr.2013.12.004

CrossRef Full Text | Google Scholar

Smith, K. F., McFarlane, N. M., Croft, V. M., Trigg, P. J., and Kearney, G. A. (2003). The effects of ploidy and seed mass on the emergence and early vigour of perennial ryegrass (Lolium perenne L.) cultivars. Aust. J. Exp. Agric. 43, 481–486. doi: 10.1071/EA01130

CrossRef Full Text | Google Scholar

Wang, R. X., Hai, L., Zhang, X. Y., You, G. X., Yan, C. S., and Xiao, S. H. (2009). QTL mapping for grain filling rate and yield-related traits in RILs of the Chinese winter wheat population Heshangmai× Yu8679. Theor. Appl. Genet. 118, 313–325. doi: 10.1007/s00122-008-0901-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Yeats, T. H., and Rose, J. K. C. (2013). The Formation and function of plant cuticles. Plant Physiol. 163, 5–20. doi: 10.1104/pp.113.222737

PubMed Abstract | CrossRef Full Text | Google Scholar

Zadoks, J. C., Chang, T. T., and Konzak, C. F. (1974). A decimal code for growth stages of cereals. Weed Res. 14, 415–421. doi: 10.1111/j.1365-3180.1974.tb01084.x

CrossRef Full Text | Google Scholar

Zapata, C., Silva, P., and Acevedo, E. (2004). Comportamiento de Isolíneas de Altura en Relación con el Rendimiento y Distribución de Asimilados en Trigo. Agric. Técnica 64, 139–155. doi: 10.4067/S0365-28072004000200003

CrossRef Full Text | Google Scholar

Zhang, Z., Wang, W., and Li, W. (2013). Genetic interactions underlying the biosynthesis and inhibition of β-diketones in wheat and their impact on glaucousness and cuticle permeability. PLoS ONE 8:e54129. doi: 10.1371/journal.pone.0054129

PubMed Abstract | CrossRef Full Text | Google Scholar