Yield and Quality Performance of Traditional and Improved Bread and Durum Wheat Varieties under Two Conservation Tillage Systems

: The increasing spread of conservation agriculture demands that the next generation of wheat varieties includes cultivars capable of maintaining satisfactory yields with lower inputs and under uncertain climate scenarios. On the basis of the genetic gains achieved during decades of selection oriented to yield improvements under conventional crop management, it is important that novel breeding targets are deﬁned and addressed. Grain yield, yield-related traits, and phenological and morphological characteristics, as well as functional quality parameters have been analyzed for six varieties each of bread and durum wheat, under minimum tillage and no-tillage. During the three-year experiment, the climatic conditions at the ﬁeld trial site were characterized by low rainfall, although di ﬀ erent degrees of aridity—from moderate to severe—were experienced. Di ﬀ erences were found between these two soil management practices in regard to the varieties’ yield stability. A positive inﬂuence of no-tillage on traits related to grain and biomass yield was also evidenced, and some traits among the examined seemed involved in varietal adaptation to a particular non-conventional tillage system. The study also conﬁrmed some breeding targets for improved performance of wheat genotypes in conservation agroecosystems. These traits were represented in the small set of traditional varieties analysed.


Introduction
Cropland soils have low organic matter content in many dry and semi-arid areas, like the Mediterranean zone, due to climatic limitations that affect net primary productivity. Water deficit and low soil fertility make these agroecosystems vulnerable, both to cyclic drought events and the current process of global climate change. With degraded soils and unpredictable weather patterns, return on investment is uncertain, which increases the economic and social risks for farmers. Under such conditions, one approach to improve crop production, while minimizing soil and other natural resources degradation, is to adopt crop-nutrient-water-land management practices that are collectively included in the term 'conservation agriculture' [1]. The practice is not only effective in enhancing

Site Description and Experimental Procedure
The study was performed at the INIA experimental farm La Canaleja (Alcalá de Henares, Madrid, Spain: 40 • 32 N and 3 • 20 W; 600 m). The soil is a sandy-loam Calcic Haploxeralf, pH 8 with low inherent fertility. The climate of the area is semi-arid continental-Mediterranean with dry summers. Precipitations are irregularly distributed over the year and from one year to another. The average rainfall over the last 25 years is 365 mm per year, the annual precipitation being below the average 11 years during that period . Daily meteorological data were recorded over the period of study (autumn 2015 to summer 2018) at the weather station located in the growing area. The average monthly precipitation, and minimum and maximum temperatures from October to June are shown in Figure S1 (Supplementary Material). For environment characterization, some agroclimatic variables were estimated for two growing periods: from sowing to anthesis (SA) (Zadoks stage 65; [24] and from anthesis to physiological maturity (AM) (Zadoks stage 87). The mean dates of anthesis and physiological maturity of all plots were used to define the SA and AM periods in each environment (i.e., agronomic year) [25]. The variables calculated for these two growing periods were: length of growing season expressed as thermal time (LGS, GDD growing degree-days), average daily mean, maximum and minimum temperatures (Tm, Tmax and Tmin, respectively, • C), mean daily relative air humidity (RH, %), average daily rainfall (Rainfall, mm) and potential evapotranspiration (ETo, mm) estimated by the Pemnan-Monteith method [26]. Thermal-time was calculated by summing the GDD accumulated during the period of interest. For a given day, GDD were calculated as the difference between the average air temperature and the base temperature [27], considering 0 • C as the base temperature. When calculating the average air temperature, the upper limit for Tmax was established at 25 • C [28].
The study, conducted over three seasons (2015-2016, 2016-2017 and 2017-2018), was designed on the basis of a split-plot randomized complete block, with three blocks divided into two plots of 1008 m 2 (67.2 m width and 15 m length). Each plot was divided into two equal parts (504 m 2 ), where fallow-wheat rotation was alternated during the study. Experimental treatments included two tillage practices (minimum tillage and no-tillage), which were assigned to the main plots, and 12 wheat genotypes in subplots of 42 m 2 (16 rows of 15 m length and 17.5 cm row spacing). Minimum tillage (MT) consisted of chisel ploughing at a depth of 15 cm, while no-tillage (NT) consisted of direct seeding. Fallow and tillage practices were established two years before the study.
Wheat was sown in autumn and harvested in June. Fertilizer was applied before sowing (200 kg ha −1 of NPK 8-24-8) and in spring (200 kg ha −1 of 26% ammonium nitrosulphate). Weeds were controlled according to the tillage type: by a cultivator pass before sowing and post-emergence herbicides in MT, and by herbicide in NT. During the fallow seasons, weeds were controlled with chisel and herbicides in MT and NT plots, respectively. After harvesting, all crop residues were chopped and left on the soil surface, regardless of the tillage system. No fungicides or insecticides were applied during the experiment. Figure S2 (Supplementary Material) shows images of the field trials during the three-year experiment.

Soil Properties
Soil samples were collected at the end of March in the second season (2016)(2017), between the wheat stages of stem elongation and ear development (Zadoks stages 39 to 47). Soil samples were taken from the six crop plots, covering the two tillage practices and the three blocks. For each plot, soil samples were taken from three locations in a zig-zag pattern at depths of 0-15 cm and 15-30 cm. Composite soil samples were air-dried and ground to pass a 2-mm mesh. A portion of the dried ground soil was then milled for soil organic carbon (SOC) and nitrogen analysis. SOC concentration was determined using the Walkley-Black wet oxidation method. Total soil N (SN) was determined by the Kjeldahl method [29]. Soil ammonium (N-NH 4 + ) and nitrate (N-NO 3 − ) content was first extracted with KCl 2 M and then determined with an autoanalyser (FIAStar 5000). Additional soil variables were also determined, including C/N ratio, and P (Olsen method), assimilable K, Mg and Ca concentration [30]. After grain harvest in the third year, the saturation time (s) and rate of water infiltration in the soil surface (L s −1 ) were measured in wheat and fallow plots using a metal ring 53 cm in diameter [31]. Saturation time was the time required for visible soil surface saturation-when the water poured into the ring starts to run off outside the ring. Water infiltration rate was determined as the quotient between the volume of water infiltrated (water poured into the ring) and the saturation time. At least three measures following a zig-zag pattern were taken from each tillage system and block to calculate average values.

Yield Components and Agronomic Traits
The plants contained in a one-meter-long row were cut at the ground level. The number of tillers and spikes in each sample were counted, and the aerial portion was weighted after being oven-dried at 70 • C for 48 h. The following yield components were calculated for each replicate: Grain weight expressed on a 12% moisture basis (GW, g m −2 ), tiller number (TN, m −2 ), spike number (SN, m −2 ), grain number (GN, m −2 ), number of grains per spike (GS), biomass weight (BW, g m −2 ), harvest index (HI, GW: BW ratio), thousand kernel weight (TKW, g) and test weight (TW, kg hL −1 ).
Days to heading (DH) and maturity (DM) were recorded when more than 50% of the main spikes within a plot had reached Zadoks stage 55 and 87, respectively. The grain filling period (GF, days) was estimated as the number of days between anthesis and maturity. Grain filling rate (GFR, mg GDD −1 ) was then obtained as the quotient between average dry weight per grain and thermal time from anthesis to maturity. At harvest maturity, five plants chosen at random from the centre of the plots were used for determining plant height (PH, cm) and peduncle length (PL, cm). PH was measured from the ground level to the top of the main spike, excluding awns, and PL from the uppermost (last) node on the stem to the spike collar.

Technological Quality Evaluation
Grain protein content (GPC, % at 14% moisture basis) was measured in wholemeal flour samples by near-infrared reflectance analysis using the Technichon Infranalyzer 300. Gluten strength was evaluated by the sodium dodecyl sulfate sedimentation volume (SDSS, mm) test [32]. GPC was analyzed individually in the three replicated plots, while one single flour sample was analyzed for SDSS per genotype and tillage system. For both parameters, technical duplicates were prepared.
To analyse the glutenin composition of the wheat genotypes, endosperm proteins were extracted from wholemeal flours [33]. Electrophoresis of reduced and alkylated glutenins (high molecular weight subunits, HMW-GS, and low molecular weight subunits, LMW-GS) was performed on sodium dodecyl sulphate polyacrylamide gel electrophoresis SDS-PAGE (12% polyacrylamide gels). The alleles at the Glu-1 loci, controlling the HMW-GS, and at Glu-3 and Glu-B2 loci, controlling LMW-GS, were identified [34].

Data Analyses
Analysis of variance of plant-evaluated data (i.e., agronomic and quality traits) was performed separately for bread and durum wheat, using the generalized linear mixed model with tillage treatments, genotype and year (i.e., environment) considered as fixed effects and the block as a random variable. Genotype means were compared within the same species by Fisher's protected least significant difference. Similar statistical approaches were followed to determine the effects of depth and tillage treatment, and their interactions, on soil properties and to compare soil composition data.
Multivariate analysis with principal component analyses (PCA) was applied to analyse the relationships between agronomic variables and wheat genotypes under minimum tillage and no tillage in different environments. The sampling adequacy of individual and set variables were verified by the Kaiser-Meyer-Olkin measure (>0.50) and Bartlett's test of sphericity (<0.05). Selected main components had eigenvalues over 1.0. All variables had communality values > 0.5, except for TW in durum wheat under MT in 2017-2018 (TW was removed from that analysis). Cophenetic correlation coefficients were higher than 0.9 in all the PCAs.
Relationships between variables were examined using Pearson correlation coefficients (p < 0.05). Yield stability of the varieties, based on grain weight per square meter (GW) data, was analyzed by Wricke's ecovalence (W i ; [35]). This statistic estimates the sums of squares contribution of each genotype to the overall G × E term. The lower the value of W i , the smaller the fluctuations from the predictable yield response in different environments: For grain yield, the genotype × year interaction was partitioned according to the AMMI model (additive main effects and multiplicative interaction model) [36], and the percentage of the sum of squares explained by each interaction principal component axis (IPCA) was calculated for each species. Factorial regression analyses were performed to identify the agroclimatic variables with a major effect on the grain yield of wheat genotypes in different years. For these analyses, the agroclimatic data corresponding to the sowing to anthesis and anthesis to maturity periods were considered.
Data were analyzed using Infostat version 6.12 software [37]. This was also the software employed for graphical representations of the results.

Environment Conditions and Soil Properties
Rainfall during the crop cycle (October to June) was 292 mm, 203 mm and 374 mm for the seasons 2015-2016, 2016-2017 and 2017-2018, respectively. For reasonable wheat growing, 300 mm is enough, provided that most water input is during spring [38]. Accordingly, water input was adequate in the first and third years, but in the second year, characterised by being especially dry and hot, wheat plants likely suffered severe drought stress ( Figure S1, Supplementary Material). These climatic differences affected the duration of the crop cycle, the first and third years having longer pre-anthesis and post-anthesis periods measured in thermal time (Table S1, Supplementary Material).
The soil analyses carried out in the second year indicated that SN, N-NO 3 − , SOC, P and K decreased significantly with an increase in soil depth. Higher N-NO 3 − and P concentration at surface soil depth (0-15 cm) was observed under no-tillage ( Table 2). Significantly higher values of saturation time were observed in the crop plots under NT when water infiltration properties were determined on soils after the crop harvest of the third year. The rate of water infiltration from the surface also increased under NT, although the differences compared to MT were not statistically significant ( Figure  S3, Supplementary Material).   Table 3 shows the mean values of all the variables analysed for each year. The mean values obtained for each variety of bread and durum wheat in the three years of experiment are provided as Supplementary Material ( Figure S4).

Grain Yield and Agronomic Traits
These results evidenced the remarkable influence of the environmental conditions that occurred during the study on yield and all the other agronomic traits evaluated, which was further confirmed by the analyses of variance performed (Table S2, Supplementary Material). Variability of climatic conditions among the years significantly affected grain weight (GW) and yield components. Lower values were recorded in the driest and hottest second year (2016-2017) for all the traits in both species ( Table 3). The highest values for HI, TKW and TW were obtained in the rainiest third year (2017-2018), while no significant differences were found for GW, GN and BW between the first and third years. The morphology and phenology of bread and durum wheat plants were also influenced by extreme water stress in the second year, which reduced all the measured values.
All the yield-related traits were significantly affected by the genotype (except TN, SN and BW in durum wheat), and some of them showed a significant genotype × year interaction (Table S2 36 and Senatore Capelli showed generally low tiller and spike numbers but relatively better values in the driest than in the two wetter years. It can be added that less between-genotype differences for GN, GS, TN and SN seemed to exist among the durum wheat varieties analysed than within the set of bread wheat varieties.  The bread wheat landrace Aragon-03 produced high biomass in the three years, whereas BW of Marius and Chamorro was relatively increased in the wetter years. In durum wheat, differences among varieties were again less evident, although the landrace Jerez-36 showed higher BW values across the three years and Vitron was the variety that produced less biomass. The spring bread wheat varieties (Califa Sur and Yecora) had significantly high HI values in the three years, which was also observed for the intermediate variety Vitron among durum wheat genotypes. In both wheat crops, the landraces showed the lowest HI values, especially bread wheat Chamorro and durum wheat Senatore Capelli. In general, TKW values decreased in modern cultivars and increased in the landraces and intermediate varieties. It is remarkable the high TKW that found for the durum wheat landrace Jerez-36. Two groups of bread wheat varieties were separated according to TW values: the spring genotypes and the landrace Aragon-03 (higher values), and the winter genotypes and the landrace Chamorro (lower values). In durum wheat, the landraces' kernels had lower TW in the two former years, but TW of Jerez-36 was among the best in the third year.
Differences among varieties were quite consistent for morphology and phenology traits between years ( Figure S4, Supplementary Material). In bread or durum wheat, the landraces showed higher PH and PL values than the intermediate and modern breeding varieties. The bread and durum wheat landraces were also invariably later than the breeding cultivars. Nevertheless, significant influence of the growth habit of bread wheat breeding lines on DH and DM was found. Thus, the winter varieties Marius and Berdun were always later that the spring genotypes Yecora and Califa Sur. Additionally, Aragon-03 headed significantly earlier than Chamorro. In general, the landraces of both species had lower GF values in the three-year study. Between-varieties differences for grain filling rates were similar in the first and third year. In bread wheat, the highest, medium and lowest values corresponded to the landraces, and intermediate and modern varieties, respectively. In durum wheat, the landraces and the intermediate variety Vitron showed higher GFR, whereas lower values were obtained for the modern varieties and the intermediate Cocorit. Such trends were not observed in the second year, but the landraces Aragon-03 (bread wheat) and, especially, Jerez-36 (durum wheat) maintained quite high GFR values.
The genotype × year interaction was not statistically significant for GW. However, some genotypes exhibited remarkable differences in grain yield between years ( Figure S4, Supplementary Material). To identify the agroclimatic variables with a major effect on yield performance of the genotypes, an AMMI analysis was performed. Since no differences were evident from biplots performed separately for bread and durum wheat varieties, the combined AMMI biplot of the two species was represented ( Figure 1). The agroclimatic variables contributing the most to the yield performance of the genotypes in the three years (largest projection over IPCA1) were related to the length of the growing periods and the humidity (TT, Rainfall and RH) in the positive direction, and to heat (Tmax and Tm) in the negative direction. The mean minimum temperatures (Tmin, SA and AM) had the lowest influence on the yield performance across the three years. In agreement with their placement in the AMMI biplot, the second year (2016-2017) was characterised by having high temperatures and low humidity during both growing periods, whereas the third year (2017-2018) had opposite agroclimatic conditions (see Table S1, Supplementary Material). The first year (2015-2016) had an intermediate position, more similar to that of the third year, but with higher Tmin (SA) and lower Tmin (AM).
The location of the bread wheat cultivars in the biplot confirmed that the intermediate variety Marius (winter genotype) performed the best in the third year, whereas the spring genotypes (Yecora and Califa Sur) had the best yield performance in the hottest and driest second year. The landraces Aragon-03 and Chamorro had an intermediate position between these two most contrasting years. Among durum wheat genotypes, Cocorit and Avispa were the best adapted to the wetter environments whereas, the remaining varieties performed relatively better in the second-driest year.
According to the results of the ANOVA tests, the conservation tillage system employed did not seem to affect yield nor most agronomic traits evaluated, only TKW in bread wheat and DM in both species being affected (Table S2, Supplementary Material). However, mean comparisons analyses revealed that, in bread wheat, no-tillage practice significantly increased TKW and GFR, and TW in the third year (Table 4). In durum wheat, higher DM values were observed under no-tillage, whereas minimum tillage increased GFR in the second year and PH in the third year. The location of the bread wheat cultivars in the biplot confirmed that the intermediate variety Marius (winter genotype) performed the best in the third year, whereas the spring genotypes (Yecora and Califa Sur) had the best yield performance in the hottest and driest second year. The landraces Aragon-03 and Chamorro had an intermediate position between these two most contrasting years. Among durum wheat genotypes, Cocorit and Avispa were the best adapted to the wetter environments whereas, the remaining varieties performed relatively better in the second-driest year.
According to the results of the ANOVA tests, the conservation tillage system employed did not seem to affect yield nor most agronomic traits evaluated, only TKW in bread wheat and DM in both species being affected (Table S2, Supplementary Material). However, mean comparisons analyses revealed that, in bread wheat, no-tillage practice significantly increased TKW and GFR, and TW in the third year (Table 4). In durum wheat, higher DM values were observed under no-tillage, whereas minimum tillage increased GFR in the second year and PH in the third year.
For those traits showing significant genotype × tillage interaction (Table S2, Supplementary Material), further mean comparison analyses considering variety and tillage system were done ( Figure 2). Overall, GFR was the trait for which more varieties showed a significant influence of the tillage system, the effect seeming to be dependent on the environment conditions. Higher GFR values were mostly obtained under no-tillage in wetter years, whereas minimum tillage seemed to be more beneficial during the extremely driest, second year. Nevertheless, this trend could not be generalized for all the varieties examined (i.e., see Yecora in 2016-2017). No-tillage had also a significant positive effect on SN, GN, GS, PH and/or PL in some bread wheat genotypes during one of the wetter years. In the case of durum wheat varieties, and limited to the first year, significant increasing values of TN   For those traits showing significant genotype × tillage interaction (Table S2, Supplementary Material), further mean comparison analyses considering variety and tillage system were done ( Figure 2). Overall, GFR was the trait for which more varieties showed a significant influence of the tillage system, the effect seeming to be dependent on the environment conditions. Higher GFR values were mostly obtained under no-tillage in wetter years, whereas minimum tillage seemed to be more beneficial during the extremely driest, second year. Nevertheless, this trend could not be generalized for all the varieties examined (i.e., see Yecora in 2016-2017). No-tillage had also a significant positive effect on SN, GN, GS, PH and/or PL in some bread wheat genotypes during one of the wetter years. In the case of durum wheat varieties, and limited to the first year, significant increasing values of TN under no-tillage were also detected for Senatore Capelli, while SN of the other landrace, Jerez-36, was positively influenced by minimum tillage. It is remarkable that GFR was the only trait for which significant differences depending on the tillage treatments were detected in the driest year. It agreed with the finding that the effect of genotype × tillage interaction on data recorded in 2016-2017 was only significant for this trait (Table S2, Supplementary Material).
Given that the environment strongly affected the agronomic variables examined, multivariate PCA analysis was used to determine the relations among these traits in the two most contrasting environments, the driest second year (2016-2017) and the wettest third year (2017-2018), under each of the conservation tillage systems tested. The first two axes of the PCAs accounted for 73% to 84% of the total variance, indicating that most of the information contained in those two seasons' data could be summarized by projecting the points on the plane that these axes determined ( Figure 3).
Overall, GN was the trait most closely related with crop yield (GW) independent of the year considered and under any of the two conservation tillage systems.
In bread wheat, the relation between yield and the remaining variables under minimum tillage was quite similar in the two analysed years, GW being positively correlated with yield components GS and HI, and negatively correlated with DH and DM. Only for BW it was found that the positive relation with GW was restricted to the driest season (2016-2017). Relevant differences between years were observed, however, under no-tillage. The same trends as noted for minimum tillage were observed in 2016-2017. However, in the wettest year, GW appeared closely associated with SN and TN, whereas it was unrelated to HI and the phenology traits. The position of the bread wheat varieties in the biplot spaces obtained for the driest year showed that the earliest spring varieties, Califa Sur and Yecora, had high PC1 values (related to lower DH and DM, and higher GW), regardless of the tillage system. Among the remaining varieties, Aragon-03 and Marius performed better under no-tillage, while the opposite was observed for Chamorro and Berdun. In the wetter year, grain yield of the landraces and the intermediate variety Marius improved under no-tillage, whereas the modern varieties and Yecora were placed closer to GW under minimum tillage.
In durum wheat, the relation between yield and the remaining traits remained quite similar regardless of the tillage system and the year under study. In general terms, GW was positively correlated with HI and GS, while negatively correlated with DH and DM. Only under no-tillage, some relevant differences were detected between years for TN and SN, which occupied opposite positions along the PC2 axe in 2016-2017 and 2017-2018, and for TKW (inversely related with GW in the wetter year). The traditional varieties were always placed apart from the improved varieties. These were separated by the PC2, which was mostly associated with tillering ability and biomass development.
The yield stability of the genotypes was assessed by Wricke's ecovalence method ( Table 5). The analysis was separately conducted for GW values recorded across the three-year study in no-tillage and minimum tillage plots. For all the varieties tested, lower ecovalence values, and therefore higher stability, were generally obtained under minimum tillage. Among bread wheat genotypes, landraces had usually lower ecovalence values, especially Chamorro, although no clear effect of breeding intensity was detected. The winter variety Marius and the spring variety Califa Sur were invariably the least stable, evidencing that yield stability was not influenced by growth habit. The relative stability of some durum wheat genotypes seemed to be dependent on the tillage system. Thus, landrace Jerez-36 was the most stable under minimum tillage but quite unstable under no-tillage, while the opposite trend was found for Vitron. The high-yielding variety Avispa and the traditional Senatore Capelli were unstable with all of the tillage treatments. Given that the environment strongly affected the agronomic variables examined, multivariate PCA analysis was used to determine the relations among these traits in the two most contrasting environments, the driest second year (2016-2017) and the wettest third year (2017-2018), under each of the conservation tillage systems tested. The first two axes of the PCAs accounted for 73% to 84% of the total variance, indicating that most of the information contained in those two seasons' data could be summarized by projecting the points on the plane that these axes determined (Figure 3).    An alternative approach was followed to compare the yield stability among genotypes. The grain yield reduction of each variety with respect to the GW value of the best yielding genotype was calculated for the two more contrasting years, 2016-2017 and 2017-2018, and represented in scatter plots (Figure 4). The closer to the diagonal a genotype is placed in the scatter plot, the less its relative yield is affected by environmental changes.  An alternative approach was followed to compare the yield stability among genotypes. The grain yield reduction of each variety with respect to the GW value of the best yielding genotype was calculated for the two more contrasting years, 2016-2017 and 2017-2018, and represented in scatter plots (Figure 4). The closer to the diagonal a genotype is placed in the scatter plot, the less its relative yield is affected by environmental changes.    Overall, the bread wheat varieties were less stable under no-tillage. The spring varieties Califa Sur and Yecora were the most stable under minimum tillage, but very unstable under no-tillage. The landrace Aragon-03 was among the most stable genotypes whatever tillage system used. In general, the relative stability of durum wheat varieties was less dependent on the tillage system. However, Avispa and Vitrón were more stable under minimum tillage and no-tillage, respectively. These results, based on GW data from the wettest and driest years, were broadly in agreement with data for the whole set of observations obtained in the three-year experiment by the ecovalence method (Table 5).

Technological Quality
The ANOVA test indicated that the quality traits were significantly affected by year, genotype and genotype × year interaction in the two species, while the tillage system and interactions of tillage with year and with genotype did not have any relevant effect (Table S2).
The highest GPC values were obtained in 2016-2017, whereas the lowest contents were found in 2017-2018, either in bread or durum wheat genotypes (Table 3), which supported an inverse influence of the environment on GPC and GW. Correlation analyses performed separately for varieties within each breeding period detected a significant negative correlation between GPC and GW, except for the landraces of durum wheat. However, when the analyses were performed for each of the years of study, the correlation was only significant (and negative) for bread wheat in 2015-2016 (Table S3). The landraces of both species had higher protein content than modern and intermediate genotypes ( Figure S4, Supplementary Material), which could also be related to their generalized low yield.
The average SDSS values were significantly lower in the rainiest third year than in the first and second years (Table 3). Among the bread wheat varieties, Chamorro and Yecora had the best quality performance across the study ( Figure S4, Supplementary Material). The modern varieties Califa Sur and Berdun suffered greater penalties in the first and third wetter year, while Marius and Aragon-03 invariably showed the lowest SDSS values. In durum wheat, the old genotype Jerez-36 showed the best SDSS values in the three years, the modern variety Don Pedro performing the worst. A significant positive correlation between GPC and SDSS was detected for varieties within each breeding period, except for the intermediate varieties of bread wheat (p = 0.056). Separate analyses for each year revealed significant positive correlations in the first year for both species and in the third year for durum wheat varieties, whereas no correlation between GPC and SDSS was detected in the driest season's analyses (Table S3).
The electrophoretic analyses indicated that the bread wheat varieties had different allelic composition for the HMW-GS loci. For the Glu-D1 locus, main genetic determinant of bread wheat rheological properties, Califa Sur and Yecora expressed subunits 5 + 10, Marius and Chamorro expressed subunits 4 + 12, and Berdun and Aragon-03 expressed subunits 2 + 12. In durum wheat, the analyses of the LMW-GS revealed that all the genotypes possessed the alleles Glu-A3j (subunit 6), Glu-B3r (subunits 2 + 4 + 15 + 19) and Glu-B2a (subunit 12). Higher variability was observed for the HMW-GS, although all the varieties except one possessed the null allele at Glu-A1, and only three alleles were detected at the highly variable Glu-B1 locus ( Figure S5).

Discussion
The next generation of crop cultivars will have to satisfactorily respond to the main concerns of current agroecosystems: increasing food demand, soil degradation, overuse of inputs and climate uncertainty. In the case of wheat, functional quality requirements must also be kept in mind. Our study evaluated the yield and quality performance of a set of bread and durum wheat varieties under two sustainable soil tillage systems. The panel of genotypes included representatives of traditional and improved varieties to identify particular traits for which old materials can aid the genetic improvement of higher-yielding cultivars developed for intensive, high-input agriculture systems.

Overall Varieties' Performance under Conservation Soil Management
The climatic conditions during the field trial experiment were quite typical of dry regions of Spain, characterised by erratic amounts and distribution of rainfall [39]. Even in the rainiest season (2017-2018), rainfall was within the range of a semi-arid climate (e.g., [40,41]). It surely implies that the plant-growing periods were always under some degree of drought stress, from moderate to severe, which might be the future climate scenario in many template zones that have traditionally received adequate water input [41].
All the traits examined were significantly influenced by the environment. The most important agroclimatic variables affecting genotype adaptation to the environment, in terms of grain yield (i.e., GW), were temperature (daily mean and maximum), rainfall, and length of the cycle expressed as thermal time (Figure 1). No significant differences in GW between the first (2015-2016) and third (2017-2018) years were detected in any of the two species, probably due to adequate water input during the reproductive period and grain filling stage (from March to May) [38]. In contrast, the extremely drought environment in the second year (2016-2017) greatly penalised grain yield and yield components ( Table 3). The water deficit also reduced the growth cycle (DH, DM and GF) and decreased the daily rates of translocation of carbohydrate reserves from vegetative organs to the grain (GFR), which is in agreement with lower TKW values, as reported here and in other studies [42,43]. These results confirmed the already contrasted role of the amount and distribution of rainfall during the growth period in wheat yield under Mediterranean rainfed conditions [43][44][45].
With the exception of Aragon-03, landraces were always the genotypes with the lowest grain yields. This reflects the improved partitioning of assimilates to the grain that provoked the introduction of dwarfing alleles in wheat cultivars [46]. However, the relative yield performance of each commercial, breeding-derived cultivar varied depending on the year. Among the modern and intermediate varieties of bread wheat, grain yield in response to the different environmental conditions was associated with growth habit, which is in agreement with Sanchez-Garcia et al. [45]. Under severe water stress, the earliest spring varieties performed relatively better than the winter varieties, whereas the latter had the best grain yields in more humid environments. The landrace Aragon-03 behaved acceptably well across the environments, with yield decrements of 20-25% respect to the best performing variety at any of the trial seasons (see Figure 1 and Figure S4). Chamorro behaved also as very stable based on the ecovalence values, but this landrace showed the lowest yield among all the varieties tested. As confirmed in our study, Aragon-03 is characterized by high tillering and spike number, which may compensate for the low spike fertility inherent to landraces. It can be noted that this old genotype, selected from a Spanish landrace before 1940, was the leading variety in Spain during the period 1960-1976 due to its wide adaptation to drought and cool environments [47].
In durum wheat, the tested varieties showed less overall differences in yield stability than in bread wheat. The high stability of the intermediate genotypes and the finding that the landraces performed relatively better during the driest and hottest year concurred with findings reported by Subira and co-workers [25]. These authors suggested that Senatore Capelli has a good adaptation to drought stresses, which agrees with its lowest score for IPCA1, strongly related to rainfall. However, the Spanish landrace Jerez-36 showed better yield performance and stability, with yield loss rates clearly lower than those of Senatore Capelli (17-37% for Jerez-36; 36-46% for Senatore Capelli). Its higher GFR and TKW values, especially in the driest environment, support the fact that Jerez-36 has a greater water use efficiency during grain filling, thus resulting in heavy grain formation.
The presence of dwarfing alleles in varieties released from the Green Revolution times has been associated with physiological changes leading to superior number of grains rather than with changes in individual kernel weight [48]. TKW has indeed been included among the traits for which landraces can provide genetic diversity that may contribute to improved yields in rainfed agroecosystems, where heavier kernels may compensate for lower spike fertility [49,50]. The direct physiological relation between grain filling rate and weight provides another target trait for adaptation to lower input wheat cropping (e.g., [51,52]). The landraces examined here have generally showed the greatest GFR values, especially under moderate water scarcity. Furthermore, even in the driest season, GFR of bread wheat Aragon-03 and durum wheat Jerez-36 have been significantly higher that that found in modern cultivars.
Crop biomass is another interesting trait in the context of non-conventional management agro-systems. Decreased biomass and the subsequent reduced ability to compete against weeds represents a problematic issue in conservation agriculture, where weed pressure can adversely affect crop yield, mainly in the long term (e.g., [53]). On the other hand, unharvested crop residues contribute to generate above-ground and below-ground biomass, whose decomposition and recirculation substantially enhances soil quality. It can be noted that soil organic matter constitutes a store of nutrients and water, and controls the common erosive processes, its low level in dryland representing an important vulnerability factor for the sustainability of Mediterranean agriculture [54,55]. Additionally, increased aerial and root biomass contribute to improve wheat plant's ability to capture resources and can lead to a greater ground cover, reducing soil temperature and water lost by evaporation. All these beneficial effects of plant biomass under conservation agriculture systems are, however, constrained when the crop cultivar shows a genotype-dependent limited vegetative growth. The higher biomass production that characterizes landraces (this study; [1,56,57]) not only helps crop plants compete better with weeds, thus reducing herbicide applications, but also enhances the quality of degraded rainfed soils, clearly being advantageous for sustainable agroecosystems [1,58]. Among the landraces examined in the study, Aragon-03 and Jerez-36 can again be highlighted as good genotypes under either moderate or severe dry conditions. In the case of Aragon-03, it might be related to the remarkable tillering ability that characterizes this old bread wheat genotype.
As already noted for yield and agronomic traits, technological quality parameters were also greatly affected by climatic conditions. Some yield diluting effect is surely on the basis of both, the straight inverse relation between GPC and annual rainfall [59], and the higher GPC values obtained in the usually lower yielding landraces. Aragon-03 is indeed the only genotype among all examined that presented high GPC and acceptable grain yield. The SDSS values decreased significantly in the rainiest year. The remarkable reduction observed in durum wheat flour (Table 3) was probably influenced by the extremely low protein content, known to have a deleterious effect on wheat rheological properties [60]. No relation with the period of release was detected for gluten strength, as estimated by the SDSS test. It is well known that this quality trait has a remarkable genotype influence. In the case of bread wheat, the HMW-GS encoded by the Glu-D1 locus are the main genetic determinants of gluten quality [61,62]. This explains the consistently higher SDSS values of Yecora, intermediate of Berdun, and lower of Aragon-03 and Marius, associated with the presence of the subunits 5 + 10 (good quality), 2 + 12 (medium quality) and 4 + 12 (bad quality), respectively [61][62][63]. The good values of Chamorro are inconsistent with its glutenin composition (4 + 12 Glu-D1 subunits) but were not unexpected. This traditional variety, largely appreciated by its bread-making quality, is currently being used in the elaboration of artisanal baking products [64]. In durum wheat, all the varieties had the same allelic combination for LMW-GS, the most relevant prolamins for gluten quality. The worse values of Senatore Capelli, Don Pedro and Cocorit can be due to the negative effect on gluten strength of HMW-GS 20x + 20y encoded by the Glu-B1 locus [65].

Minimum Tillage versus No-Tillage
The results discussed above refer to the average behaviour of the varieties under any of the two non-conventional tillage systems tested. However, our study has revealed some differences between soil conservation management practices regarding varieties' performance and stability as well as on the relevance of some particular traits (i.e., TKW and GFR) for improved grain yields. Grain yield itself was not significantly affected by tillage management, although some positive effect of no-tillage was evidenced on GW and most yield-related traits, especially in bread wheat (Table 4). This could be related to the overall beneficial effect of no-tillage over reduced tillage on soil components (especially, for N-NO 3 − and P), saturation time and water infiltration. These observations are in agreement with other studies that have reported significant increases in soil N or water content under no-tillage [4,7,9,10]. Some modifications of soil properties might be associated not merely with the reduced breakdown of macro-aggregates but also with the slower decomposition of crop residues when no-tillage is practiced [6,9,66]. Earlier results on the differential effects of conservation tillage practices on wheat yield are quite controversial. Some studies have reported a positive effect of no-tillage [7,12,67], whereas others have found negative effects [9,68] or no significant differences between both systems [69,70]. It must be noted that only the studies of Khorami et al. [9] and Carr et al. [70] were based on the analysis of more than one genotype (4 and 5, respectively) in different environments. Comparison of the mean GW values obtained for each variety under minimum tillage and no-tillage reveals different trends among the genotypes tested in the present study ( Figure 5). Furthermore, significant genotype × tillage treatment interactions were detected for some relevant yield-related traits in less severe drought conditions (Table S2, Supplementary Material). Therefore, not only the specific varieties analysed but also the different climatic conditions that affected the field trials conducted might well explain such disparate conclusions. reduced breakdown of macro-aggregates but also with the slower decomposition of crop residues when no-tillage is practiced [6,9,66]. Earlier results on the differential effects of conservation tillage practices on wheat yield are quite controversial. Some studies have reported a positive effect of no-tillage [7,12,67], whereas others have found negative effects [9,68] or no significant differences between both systems [69,70]. It must be noted that only the studies of Khorami et al. [9] and Carr et al. [70] were based on the analysis of more than one genotype (4 and 5, respectively) in different environments. Comparison of the mean GW values obtained for each variety under minimum tillage and no-tillage reveals different trends among the genotypes tested in the present study ( Figure 5). Furthermore, significant genotype × tillage treatment interactions were detected for some relevant yield-related traits in less severe drought conditions (Table S2, Supplementary Material). Therefore, not only the specific varieties analysed but also the different climatic conditions that affected the field trials conducted might well explain such disparate conclusions. Yield stability of bread wheat genotypes was markedly higher when reduced tillage was practised (Table 5; Figure 4), which was in agreement with the conserved position of the varieties on the planes determined by the PCA axes (Figure 3a, minimum tillage). This is likely linked with the finding that the relation between GW and relevant yield-related and adaptive traits, like HI and precocity, may differ under no-tillage depending on the climate conditions, whereas remains quite constant under minimum tillage across the environments (Figure 3a). A similar though less marked trend of higher stability under minimum tillage has been found in durum wheat, where tillering ability and TKW were the traits whose relation with yield differed between years under no-tillage.
Honsdorf and co-workers determined the effects of genotype, tillage, and genotype × tillage interactions on yield performance of a panel of bread and durum wheat genotypes, all of them developed by CIMMYT and released after 1964 [71]. The materials were cultivated on raised beds under different tillage systems (conventional soil management and no-tillage) and water regimes (full and reduced irrigation). The results of that study cannot be directly compared with the results found here because of the different experimental conditions and materials assayed. Nevertheless, these authors also found that genotype × tillage interactions were more frequent under wetter conditions, and more relevant in bread wheat than in durum wheat. One reason adduced for this latter finding was because the group of bread wheat varieties analyzed might have comprised a larger genetic Yield stability of bread wheat genotypes was markedly higher when reduced tillage was practised (Table 5; Figure 4), which was in agreement with the conserved position of the varieties on the planes determined by the PCA axes (Figure 3a, minimum tillage). This is likely linked with the finding that the relation between GW and relevant yield-related and adaptive traits, like HI and precocity, may differ under no-tillage depending on the climate conditions, whereas remains quite constant under minimum tillage across the environments (Figure 3a). A similar though less marked trend of higher stability under minimum tillage has been found in durum wheat, where tillering ability and TKW were the traits whose relation with yield differed between years under no-tillage.
Honsdorf and co-workers determined the effects of genotype, tillage, and genotype × tillage interactions on yield performance of a panel of bread and durum wheat genotypes, all of them developed by CIMMYT and released after 1964 [71]. The materials were cultivated on raised beds under different tillage systems (conventional soil management and no-tillage) and water regimes (full and reduced irrigation). The results of that study cannot be directly compared with the results found here because of the different experimental conditions and materials assayed. Nevertheless, these authors also found that genotype × tillage interactions were more frequent under wetter conditions, and more relevant in bread wheat than in durum wheat. One reason adduced for this latter finding was because the group of bread wheat varieties analyzed might have comprised a larger genetic variation than the durum wheat genotypes. This might also hold for the differences found in the present study between bread and durum wheat regarding not only yield stability but also the interactions of genotypes with either tillage systems and environment conditions.

Conclusions
The present study has provided some relevant clues that may be of help for wheat varietal choice when conservation tillage practices are implemented in drought-prone zones. All the old varieties examined, except bread wheat cultivar Aragon-03, have shown too low yields to be competitive in terms of crop productivity. However, their advantages over commercial varieties have been described for some yield-related traits that have been identified as potential targets to achieve better grain yields under soil conservation management such as TKW, GFR, biomass production and tillering ability [11]. The high variability for these traits reported in Mediterranean landraces (e.g., [1,15,16]) supports the potential of landraces to improve variety adaptation to conservation agriculture. Further research is then strongly recommended to evaluate more diversified landrace in order to select the most adapted to sustainable agriculture practices under different environments. No significant influence of tillage management, tillage × environment, or genotype × tillage management has been detected for GPC or SDSS values. This makes it possible to assume that wheat quality of any variety eventually selected will be unaffected whatever the conservation tillage practice utilized.
Supplementary Materials: The following are available online at http://www.mdpi.com/2071-1050/11/17/4522/s1, Table S1: Values of agroclimatic variables for SA and AM growing periods, Table S2: Results of ANOVAs tests for all variables examined, Table S3: Correlation coefficients between GPC and GW, and between GPC and SDSS, Figure S1: Monthly rainfall, mean Tmax and mean Tmin during the three years of study, Figure S2: Images of MT and NT plots, Figure S3: Water infiltration rate and saturation time in the crop and fallow plots, Figure S4: Mean comparison of all traits between genotypes per year of study, Figure