Next Article in Journal
The Attenuated Secretion of Hyaluronan by UVA-Exposed Human Fibroblasts Is Associated with Up- and Downregulation of HYBID and HAS2 Expression via Activated and Inactivated Signaling of the p38/ATF2 and JAK2/STAT3 Cascades
Next Article in Special Issue
Genes Associated with Foliar Resistance to Septoria Nodorum Blotch of Hexaploid Wheat (Triticum aestivum L.)
Previous Article in Journal
Generation of Knockout and Transgenic Zebrafish to Characterize Abcc4 Functions in Detoxification and Efflux of Lead
Previous Article in Special Issue
Comparative Proteomic Analysis of Developmental Changes in P-Type Cytoplasmic Male Sterile and Maintainer Anthers in Wheat
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Abiotic Stress Response of Near-Isogenic Spring Durum Wheat Lines under Different Sowing Densities

1
Centre for Agricultural Research, Agricultural Institute, ELKH, 2462 Martonvásár, Hungary
2
Department of Agricultural and Food Sciences, University of Bologna, 40126 Bologna, Italy
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(4), 2053; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22042053
Submission received: 30 November 2020 / Revised: 2 February 2021 / Accepted: 7 February 2021 / Published: 19 February 2021
(This article belongs to the Special Issue Genetics and Breeding of Wheat)

Abstract

:
A detailed study was made of changes in the plant development, morphology, physiology and yield biology of near-isogenic lines of spring durum wheat sown in the field with different plant densities in two consecutive years (2013–2014). An analysis was made of the drought tolerance of isogenic lines selected for yield QTLs (QYld.idw-2B and QYld.idw-3B), and the presence of QTL effects was examined in spring sowings. Comparisons were made of the traits of the isogenic pairs QYld.idw-3B++ and QYld.idw-3B−− both within and between the pairs. Changes in the polyamine content, antioxidant enzyme activity, chlorophyll content of the flag leaf and the normalized difference vegetation index (NDVI) of the plot were monitored in response to drought stress, and the relationship between these components and the yield was analyzed. In the case of moderate stress, differences between the NIL++ and NIL−− pairs appeared in the early dough stage, indicating that the QYld.idw-3B++ QTL region was able to maintain photosynthetic activity for a longer period, resulting in greater grain number and grain weight at the end of the growing period. The chlorophyll content of the flag leaf in phenophases Z77 and Z83 was significantly correlated with the grain number and grain weight of the main spike. The grain yield was greatly influenced by the treatment, while the genotype had a significant effect on the thousand-kernel weight and on the grain number and grain weight of the main spike. When the lines were compared in the non-irrigated treatment, significantly more grains and significantly higher grain weight were observed in the main spike in NIL++ lines, confirming the theory that the higher yields of the QYld.idw-3B++ lines when sown in spring and exposed to drought stress could be attributed to the positive effect of the “Kofa” QTL on chromosome 3B.

1. Introduction

In regions with a continental climate prone to drought, adaptability, drought tolerance and yield potential are closely related traits. The maintenance and improvement of productivity therefore require the breeding of drought-tolerant cultivars with better water-use efficiency and adaptability in changing environments and the more widespread introduction of water-conserving agricultural systems [1].
The temporal distribution of rainfall and the quantity of available water exhibit periodical changes not only at a given geographical location, but also within the crop growing season. Short periods of water deficit may occur at any time during the vegetation period, exerting varying effects on the plants, and even a mild level of water deficit may lead to severe drought symptoms if it is prolonged. The unforeseeable appearance of drought exerts an effect on yield levels: genotypes adapted to drought in one phenophase may be unable to tolerate water deficit in another phenophase [2]. The complexity of drought stress is aggravated by the fact that it may be caused by many phenomena, which also influence its severity. These include water deficit, low atmospheric moisture content, intense radiation, heat shock and high salt concentration, all of which interact with the soil type [3,4]. The combined occurrence of two different types of abiotic stress, such as drought and heat, may have an additive effect, thus causing much greater damage [5,6]. Water deficit in plants induces physiological and biochemical responses which depend not only on the intensity and duration of the stress, but also on the developmental stage of the plant and the genetically determined stress tolerance [7,8]. A decline in soil moisture content is first sensed by the root hairs, which then pass this information on to the shoots via a signal transduction system (a rise in abscisic acid content). The first detectable symptom is a loss of hydrostatic pressure (turgor) in the cells, as a consequence of which the plants close their stomata, leading to a cessation of tissue tension, cell wall ductility and growth [9]. The reduced light fixation and photosynthetic activity result in smaller biomass [10,11], lower plant height [12,13], smaller leaf area [14], fewer leaves per plant, smaller leaves and shorter leaf duration [15]. Tardieau et al. (2000) reported a shortening of the elongation zone of maize leaf tissue in response to poor water supplies and reduced transpiration, leading to a drop in the relative elongation rate [16]. Xiong et al. (2006) recorded restricted cell division in the lateral roots of Arabidopsis plants, with greater primary root growth [17]. Similar results were recorded by Sharp et al. (1994) in maize [18]. Lower chlorophyll content leads to earlier senescence, and this aging process is accompanied by the gradual cessation of respiration and photosynthetic processes, especially that of the sensitive photosystem II (PSII) photochemistry [19]. In the course of drought stress the photosynthetic organelles of the plant cells (the chloroplasts) are damaged by the formation of reactive oxygen species (ROS) [20]. The collapse of the photosynthetic apparatus leads to a reduction in the capacity and efficiency of photosynthetic energy storage [21,22]. Leaf senescence is a regulated process leading to the transfer of nutrients from older into younger leaves, and finally into the flag leaf. Photosynthesis in the flag leaf produces 30–50% of the assimilates required for grain filling, so this organ makes the greatest contribution to the supply of nutrients and photo-assimilates [23]. Both precocious and greatly delayed senescence may cause disturbances in nutrient remobilization [22,24], thus influencing the final yield [25]. As the result of more rapid senescence, Yang and co-workers (2000, 2001) observed increased nutrient mobilization from the leaves in wheat and rice cultivars, resulting in a shorter grain-filling period and smaller grain weight [26,27]. Many authors have reported that early and/or efficient nutrient remobilization is correlated with higher contents of protein and useful microelements—such as iron and zinc—in the grain [28].
The oxidative stress occurring in the course of water deficit and aging processes leads to the peroxidation of membrane lipids and the degradation of the ribulose-1-5-bisphosphate-carboxylase/oxygenase enzyme (RubP/carboxylase or Rubisco) in the chloroplasts due to the substantial increase in reactive oxygen species (e.g., oxygen-centered free radicals, H2O2) [29]. Rubisco is one of the first enzymes to be degraded, with a negative influence on plant photosynthesis and on the nitrogen metabolism [30]. Toxic metabolites (13-hydroperoxylinoleic acid, malondialdehyde, 4-hydroxyalkenals) are produced in the course of peroxidation, causing not only severe yield losses but also, in some cases, the death of the plant [31].
Drought stress during the reproductive phase of plant development, especially during meiosis and flowering, may result in yield losses as high as 40–70% [32]. Intact pollen grains and ovaries, an undisturbed double fertilization process, and the optimum development of the embryo and endosperm are all required for grain yield formation. The carbohydrate metabolism is disturbed by drought stress, resulting in poor development of the inner layer of the pollen wall (intine) and the reduced accumulation of reserve nutrients in the cytoplasm of the vegetative cells in the pollen, leading to poor attachment of the pollen on the surface of stigma papilla cells and weak pollen tube growth [33]. In experiments conducted by Jäger et al. (2013), water withholding reduced anther size in the drought-sensitive winter wheat cultivar Cappelle Desprez, with varying pollen sizes within a single flower [34]. In addition, the development rate of some pollen grains was retarded compared to the control, resulting in the formation of sterile pollen grains incapable of accumulating starch. When drought was accompanied by heat stress, the development of abnormal pollen forms was also observed. Water deficit also exerts a negative effect on the elongation of tissue cells in the ovary after fertilization and on the rapid division of cell nuclei in the endosperm. The consequent reduction in the number of endosperm cells leads to lower yields [35,36].
The values of yield components such as spike number per square meter, grain number per spike, number of fertile spikelets, thousand-grain weight and harvest index have been reported to decline as the result of physiological processes occurring in response to drought stress [12,13,37].
Creating a suitable population for examining QTL effects is a complex task because differential gene expression is caused not only by the trait of interest but also by the variation present in the genetic background [38]. One solution for establishing the functional association between the level of gene expression and a given trait is the use of a set of near-isogenic lines (NILs), which are genetically similar except for a single gene, marker or trait [39]. NILs can serve many functions, ranging from breeding purposes to the genetic analysis of complex quantitative traits. The ultimate objective of the lines determines to a large degree the choice of starting material, crossing scheme and eventually the genomic composition [40].
During the experimental period (2012/2013–2013/2014), isogenic lines of eight spring durum wheat cultivars were sown to study the effect of drought stress under field conditions using a combination of four approaches: (1) the plant responses to drought stress were tested in rain-fed and well-watered plots; (2) weather data were recorded throughout the growing period; (3) changes in the soil were monitored; and (4) the effect of the QYld.idw-3B QTL region on the yield was analyzed. The aim of this plant–environment–soil–genetics analysis was to give an accurate picture of the magnitude of changes in individual parameters and of the effect of drought stress on morphological, physiological, biochemical and yield biology traits. The present study was undertaken (1) to explore the correlation between morphological, physiological and biochemical parameters and yield components under drought conditions and (2) to analyze the effect of the QYld.idw-3B QTL region on the yield by studying near-isogenic lines of spring durum wheat with different parental alleles at the QYld.idw locus responsible for yielding ability on the 3B chromosome.

2. Results

2.1. Effect of Drought Stress and Plant Sowing Densities on Morphological, Physiological, Biochemical and Yield Component Traits in Near-Isogenic Lines of Spring Durum Wheat

The year, the genotype and the sowing density all had a significant effect (p < 0.001) on the heading, flowering and ripening dates in the field experiment, while the water supplies only had a significant influence (p < 0.001) on the ripening date (Table 1). In the years tested, low sowing density caused an average two-day delay in the heading, flowering and ripening dates under both irrigated and rain-fed conditions. The significant effect of the year, genotype, irrigation and sowing density could also be detected in the yield and yield components, the only exceptions being the number of sterile basal and apical spikelets in the main spike and the number of spikelets per main spike, which were not influenced by irrigation or sowing density. Low sowing density had a significant positive effect on the grain yield (t/ha), the grain weight per ear and the seed width. The seed length was significantly influenced by the year and the genotype, while irrigation had no measurable effect and low sowing density resulted in significantly longer grains in 2013 (Supplementary Materials Tables S1 and S2). When the NIL−− and NIL++ groups were compared in both field experiments the grain number, grain weight and thousand-kernel weight were found to be significantly higher for the NIL++ lines (Table 2 and Table 3).
Year and genotype had a significant effect on the chlorophyll content recorded in the Z45, Z65, Z77 and Z83 phenophases, while the significant effect of water supplies was only observed in phenophases Z77 and Z83 (Table 1). The effect of sowing density on the chlorophyll content was significant in all the phenophases tested, with the exception of flowering (Supplementary Materials Table S1). In the Z83 phenophase the chlorophyll content of NIL++ isogenic lines was significantly higher in the rain-fed treatments than that of NIL−− lines in both years. The genotype had a significant effect on the normalized difference vegetation index (NDVI) values in all three phenophases tested (Z45, Z65, Z83), while the year, water supplies and sowing density only exerted a significant effect in the Z83 phenophase. In both the rain-fed and well-watered treatments NIL++ lines had significantly higher NDVI values than NIL−−lines in the early waxy ripeness stage (Z83) (Table 1 and Table 2).
In rain-fed plots, low sowing density resulted in significantly longer stems but had no effect on the spike length or spikelet number. Great variation was observed between the individual isogenic pairs with respect to plant height. Genotype, irrigation and sowing density also had a significant effect on the plant height measured to the flag leaf and to the base and the tip of the ear. In response to water deficit there was a reduction in spike length but not in the number of spikelets. Line NIL3++ was the tallest, while lines NIL1−− and NIL1++ had the longest spikes in both treatments (Table 2 and Table 3). The relative water content of the flag leaf at flowering was significantly influenced by irrigation in both years, with 3% higher values on average in the irrigated plots. Sowing density had no effect on these values, as no significant difference in RWC was recorded within the treatment. There was no genotype effect, no significant difference being observed between the NIL++ and NIL1−− lines. The fresh aboveground biomass recorded at flowering (WW65) was not influenced by irrigation or sowing density, but the dry weight (DW65) was significantly higher in the irrigated plots in both years. No significant difference was detected between the NIL++ and NIL1−− lines. Lines NIL1++ and NIL3++ had the highest values of WW65 and DW65 in both years in both the rain-fed and well-watered plots.
The ascorbate peroxidase (APX) activity of the flag leaf declined significantly in response to drought stress in plots with normal sowing density, while the guaiacol peroxidase (G-POD) activity and putrescine (PUT)exhibited a significant rise. The reduction in APX enzyme activity was of a similar order of magnitude for all the isogenic lines (averaging 70%). Compared to the well-watered plots, water deficit caused a substantial rise in G-POD activity in the tested lines. The greatest relative increase was detected for the lines NIL1++, NIL2−− and NIL4−−. Water deficit induced the accumulation of PUT, SPD and SPM in both the NIL++ and NIL1−− lines, with the greatest change in the quantity of putrescine per unit dry matter (Figure 1). The agmatine and cadaverine contents were below the detection limit. Lower values of PUT, SPD and SPM content were recorded for lines NIL1−−, NIL1++ and NIL3++ in the irrigated plots, but in response to drought stress higher accumulation was observed for these lines than for the other isogenic lines.

2.2. Interaction between the Morphological, Physiological, Biochemical and Yield Component Traits in Near-Isogenic Lines of Spring Durum Wheat under Drought Stress

Principal component analysis on correlations between the traits revealed six factors with eigenvalues greater than unity, which together explained 97.89% of the total variance. The first background variable (34.05% of total variance) was most closely correlated with the yield (r = 0.92 ***), and also with various yield components (SNM, SWM, TGW, SW; r = 0.70 **−0.93 ***), with morphological (FLC, BE, TE, DW65; r = 0.84 ** to −0.87 ***) and physiological parameters (SPAD45, SPAD77, SPAD83, NDVI83; r = 0.70 **–0.92 ***), and with putrescine (r = 0.86 ***). The second factor (26.19% of total variance) was correlated with the phenophases (HD, FD, MD; r = 0.88 ***−0.97 ***), morphological traits (SKNM, FTN, FLA, EL; r = 0.62 ** to −0.91 ***), physiological parameters (NCVI45, NDVI65: r = 0.76 **–0.80 **) and ascorbate peroxidase activity (r = 0.92 ***). The chlorophyll content at flowering (SPAD65; r = 0.63), the fresh biomass (WW65; r = 0.75 **), the test weight (r = 0.80 **) and the grain width (r = 0.70 **) exhibited a correlation with the third background variable. Based on the factor–variable correlation pattern, eight of the traits tested—the grain number, grain weight, seed width and thousand-kernel weight recorded for the main spike after harvest, the dry biomass at flowering, the NDVI value recorded at early waxy ripeness, the plant height to the flag leaf and the putrescine content—could be grouped with the yield (Figure 2).
The yield of the isogenic lines exhibited a significant correlation with the thousand-kernel weight (r = 0.735 *), the number of grains per main spike (r = 0.846 **), the grain weight (r = 0.844 **), the height parameters (FLC: r = 0.800 *; BE: r = 0.679 *; TE: r = 0.717 *), the vegetation activity of the plots at flowering and grain filling (r = 0.684 *, 0.863 **), the chlorophyll content of the flag leaf in phenophases Z77 and Z83 (r = 0.786 **, 0.772 **), the number of side-tillers per square meter (r = 0.689 *) and the dry biomass at flowering (r = 0.648 *). The APX enzyme activity was in significant positive correlation with the number of spikelets on the main spike (r = 0.803 ***), the guaiacol peroxidase (GPX) enzyme activity with the number of sterile basal spikelets (r = 0.847 ***), the spermine content with the grain number per main spike (r = 0.658 **) and the putrescine content with the thousand-kernel weight (r = 0.664 **), the yield (r = 0.931 ***), the grain number per main spike (r = 0.734 **) and the grain weight (r = 0.780 **) (Figure 3).
The correlations demonstrated that the yield was determined most significantly by the number of side-tillers per square meter (which explained 40.41% of phenotypic variance), the grain number per main spike (48.87%) and the grain weight in the main spike (22.79%). The plant height to the flag leaf and to the tip of the ear correlated significantly with the yield and the thousand-kernel weight. The chlorophyll content of the flag leaf in the Z65, Z77, Z83 and Z85 phenophases was in significant correlation with the thousand-kernel weight, the grain weight in the main spike and the seed width. The grain weight of the main spike was also significantly influenced by the chlorophyll content in phenophases Z77 and Z83, while the number of grains per main spike was correlated with the chlorophyll content at the milk development stage and the flag leaf area. The light reflectance of the plots at the early dough stage (NDVI83) exhibited a significant correlation with the yield, the grain weight of the main spike and the grain number per main spike. The aboveground fresh and dry biomass at flowering was significantly correlated with the yield, the thousand-kernel weight, the grain weight of the main spike and the seed width (Table 4).

3. Discussion

Near-isogenic lines are efficient tools for investigations of the function and regulation of single genes [41,42]. The reason for using NILs is that pairs of near-isogenic lines, which have the same genetic background but differ at a single locus, can be used to turn quantitative traits into Mendelian factors, thus providing a better insight into genotypic and phenotypic interactions [43]. In the present experiment the joint effect of year, genotype, sowing density, QTL region and irrigation was investigated in the QYld.idw-3B QTL region on the short arm of chromosome 3B using near-isogenic lines of spring durum wheat. Low soil moisture content and high temperatures were recorded from the start of the milk development stage in both years of the field experiment. Significant differences in NDVI values between the NIL++ and NIL1−− pairs were already visible at the early dough stage (Z83), indicating that the QYld.idw-3B++ QTL region was able to maintain photosynthetic activity for a fairly long period, resulting in higher grain numbers and grain weight at the end of the vegetation period. Naser et al. (2020) also reported that NDVI values obtained from canopy sensing could be useful for assessing the yields of individual wheat genotypes [44]. The overall agreement between NDVI values and grain yield classes for different growth stages suggested that wheat genotypes could best be differentiated and classified between anthesis and the middle of grain filling. Fu et al. (2020) constructed a multifactor model to increase accuracy and stability and demonstrated that the correlation between the vegetation index and the yield was better in later phenophases than in the early growth stage [45]. This is particularly true in the Carpathian Basin, where drought stress occurs most frequently in the grain-filling and early dough stages after the cessation of May rainfall. The chlorophyll content of the flag leaf in phenophases Z77 and Z83 was in significant correlation with the grain number and grain weight of the main spike, the thousand-kernel weight and the plant height, confirming the likelihood that correlations between vegetation indices, yield and yield components should be sought in phenophases after flowering. In the course of breeding, the physiological parameters that are in positive correlation with the yield during drought stress (the chlorophyll content of the flag leaf, remote sensing-based indices) could be suitable for the identification of genotypes with stay-green traits, which could make them more tolerant of drought [46,47].
If the soil does not contain sufficient water for plant development, growth will be inhibited and less biomass will be formed, resulting in poor grain filling. The yield will decline due to smaller grain size or, if the stress occurred in an earlier phenophase, due to a smaller number of grains per spike. Positive turgor is essential for the maintenance of growth [48], which is why higher values of plant height parameters (FLC, BE, TE) and dry biomass at flowering (DW65) were recorded for the well-watered treatment in both years [49]. The factor–variable correlation pattern showed that morphological traits related to biomass and the vegetation activity at grain filling played an important role in the development of the yield components that determine the yield. Measurements of the vegetation activity of the plots gave a good reflection of the number of side-tillers per unit area and the area of the flag leaf, as proved by the fact that these traits were put into the same group in the course of principal component analysis.
Measurements on changes in stem dry weight can be used to estimate the level of nutrient mobilization into the grains [50,51] and many authors have reported the detection of genetic variability in the storage capacity of vegetative organs and in the efficiency of remobilization. In the present work, the dry biomass exhibited a significant positive relationship with the yield and yield components (GY, SNM, SWM) and with the height parameters (FLC, BE, TE), so its genetic variability during drought stress could make it useful in breeding [52]. These results are in agreement with the findings of Korohou et al. (2020), who reported that wheat grain yield could be estimated from the total dry matter of the leaves, stem and ear. When six models developed by these authors were evaluated on a test dataset the results showed that the models were generally accurate in terms of grain yield prediction [53].
The grain yield was greatly influenced by the water supplies, while the genotype had a significant effect on the thousand-kernel weight of the main spike. When the lines were compared in the rain-fed treatment, significantly more grains and significantly higher grain weight and thousand-kernel weight were recorded in the main spike for the NIL++ lines, confirming the possibility that the higher yield of the QYld.idw-3B++ lines in spring sowings subjected to drought stress could be attributed to the positive effect of the Kofa QTL on the 3B chromosome. Bennett et al. (2012) identified two QTL regions on chromosome 3B that exhibited a close correlation with yield and plant temperature in a doubled haploid population arising from a cross between a drought-tolerant line and a sensitive wheat cultivar (RAC875/Kukri) [54]. A Meta-QTL linked to yield [55], a QTL closely associated with plant temperature, the chlorophyll content at the milk development stage and the NDVI value during heat stress [56] were also found on this chromosome. A region correlated to yield, thousand-kernel weight and early growth vigor was identified by Bonneaue et al. (2013) in 21 different environments in a doubled haploid winter wheat population [57]. The qGYWD.3B.1 and qGYWD.3B.2 genomic regions detected on chromosome 3BS exhibited significant associations with yield, grain number per square meter, number of fertile side-tillers and spike number in a winter wheat RIL mapping population [58]. In a set of substitution lines with varying levels of drought tolerance, Dudziak et al. (2019) found that chromosomes 3B, 5A, 7B and 7D had the greatest impact on the initial wheat response to short-term stress and could therefore be critical for the genetic control of drought tolerance in wheat [59].
Li et al. (2016) reported a significant reduction in various yield components as the sowing density rose [60]. At low plant densities, Whaley et al. (2000) also found an increase in the relative growth rate of the crop, including a longer period of tiller production, better radiation use efficiency, greater green area per shoot and improved shoot survival, which prevented a decline in crop dry matter production [61]. The greater area available to each plant at low sowing density clearly had a significant positive effect on the grain weight and seed width per spike and on dry biomass (DW65). This was promoted by a number of indirect traits (higher chlorophyll content and NDVI value, longer stems, larger root system). Genotypes with good tillering ability are able to perform better at low sowing density, probably due to the better utilization of water, solar radiation and nitrogen [62,63]. As an early signal of cell division, light influences the development of new leaves and tillers, thus enhancing the photosynthetically active green plant area. This was confirmed by the significantly higher biomass, vegetation index and flag leaf chlorophyll content. It should not be forgotten, however, that for genotypes with poor tillering ability and a lower number of tillers, higher sowing density is required for the development of a satisfactory yield [64]. If wheat yields are to continue to improve, it is essential to define the optimum agronomic plant density, i.e., the minimum number of plants per unit area required to maximize yield [65].
Strong correlative and analytical evidence assigns a major role in drought tolerance to the redox regulatory and antioxidant system. The strong antioxidant system found in resurrection plants underlines the importance of the antioxidant system for surviving long periods of dehydration [66]. Changes in antioxidant enzyme activity during drought stress are greatly dependent on the plant species and cultivar and on the severity and duration of the stress. The wheat species investigated by Pour-Benab et al. (2019) exhibited diverse responses to drought stress, particularly as regards the ability of APX, GPX and SOD (superoxide dismutase) to scavenge reactive oxygen species [67]. In barley plants the activity of the antioxidant enzymes catalase, guaiacol peroxidase and superoxide dismutase was significantly higher in drought-stressed treatments, while the activity of ascorbate peroxidase did not change significantly with the treatment [68]. In the present experiment the reduction in the APX activity of the flag leaf in the non-irrigated treatment was of a similar magnitude in all the isogenic lines. By contrast, higher GPX activity was recorded in the flag leaf of all the lines in the non-irrigated treatment compared to plots with good water supplies. It would seem that sensitive plants tend to activate the glutathione-dependent scavenging system rather than the ascorbate-dependent system, which is only induced to a moderate extent, or may even be down-regulated [66]. Numerous yield components proved to be in negative correlation with the activity of the ascorbate peroxidase and guaiacol peroxidase enzymes and in significant positive correlation with the polyamine content (SW, SNM, SWM), thus confirming the role of these compounds in the stress response under field conditions. It is clear that increasing the activity of antioxidant enzymes is one strategy in the drought tolerance mechanisms in plants. When investigating changes in antioxidant enzyme activity throughout the growing period in wheat genotypes grown with normal or deficient water supplies, Huseynova et al. (2010) concluded that drought caused changes in the balance between free radical production and protective enzyme reactions [69]. Polyamines (PAs) are important for plant growth and development and for their environmental stress responses, thus acting as a novel plant biostimulant [70]. In the present experiment the accumulation of polyamines rose during drought stress in the flag leaves of the isogenic lines, with the greatest increases for putrescine, spermidine and spermine. A positive correlation was found between the polyamine content and the grain number/grain weight of the main spike in the rain-fed treatment. Water deficit induced the accumulation of PUT, SPD and SPM in both NIL++ and NIL1−− lines but the greatest change induced by the treatment in the detectable polyamine content was observed for the quantity of putrescine per unit dry matter. In experiments on the use of waste water to irrigate wheat, Aldesuquy et al. (2014) demonstrated that the application of spermine or spermidine was able to mitigate the negative effects of the treatment on the yield and yield components [71]. Liu et al. (2016) agreed that the interaction between hormones rather than the action of a single hormone was involved in the regulation of wheat-grain filling under drought [72]. However, although the overexpression of certain enzymes could enhance drought tolerance, it could also cause great delays in germination and development, thus disrupting the growing season [66]. Therefore, mapping strategies should be devised for the localization of genetic loci associated with the response of cereals to drought [59].
Due to the multiplicity of factors influencing experiments carried out under field conditions, there is limited scope for determining the degree of stress tolerance using a single parameter. Also, because of the substantial interaction between the plant and the environment, genotypic differences may also occur within the same plant [73]. In the future, the use of systems biology approaches such as genomics, transcriptomics, proteomics and metabolomics could result in new techniques for the identification of key regulators, genes, proteins and metabolites, which will help breeders to develop new, stress-tolerant genotypes [74].

4. Materials and Methods

4.1. Plant Materials

Eight near-isogenic lines of spring durum wheat, derived from four different recombinant inbred lines of a cross between the spring durum wheat cultivars ”Kofa” and ”Svevo” [75], were examined under rain-fed (RF) and well-watered conditions (WW) in the nursery of the Agricultural Institute at Martonvásár, Hungary, during the 2013 and 2014 spring cropping seasons.
The NIL pairs were obtained by MAS (Marker-assisted selection), selecting pairs differing only for the 3BS QTL (the NILs were all fixed for “Kofa KK allele” on chromosome 2B). When the QYld.idw-3B region on chromosome 3B was “Kofa KK” (KK2BLKK3BS), the NILs were coded as ++ (“NIL1++”, “NIL2++”, “NIL3++”, ”NIL4++”), and when the allele was “Svevo SS” (KK2BLSS3BS), they were coded as −− (”NIL1—“, ”NIL2−−”, ”NIL3−−”, ”NIL4−−”) [38].

4.2. Site Description and Set-Up for Field Experiment

The experiments were laid out in an unbalanced incomplete alpha lattice block design with eight replications for rain-fed and eight for well-watered conditions under two sowing densities (480 seeds/m2 normal, 320 seeds/m2 low). The geo-coordinates of the field were latitude 47°18′39.81″ N and longitude 18°46′42.83″ E (Martonvásár, Hungary), and the soil was classified as chernozem with forest residues. The genotypes were sown in the middle of March each year. Individual plots consisted of eight rows, 12 cm apart, in 0.96 m x 4 m plots. No disease was observed during the growth period and weeds were controlled using chemicals. Capacitive soil moisture sensors (5TE—Decagon Devices, Inc., Pullman, WA, USA) were placed at depths of 30 and 60 cm at eight points, all in the plot containing the line ”NIL1++”. Measurements of the moisture content (v/v %), temperature (°C) and electrical conductivity (dS/m) of the soil were made every hour throughout the growing season (Figure A1a,b).
Based on the pF curve, the water stress state was associated with a value of 21.5 v/v % (pF 3.4). The change in tension was constantly monitored with an MPS-1 tensiometer. Values of −30 to −40 kPa represent satisfactory water supplies for the crop, while values lower than −40 kPa stress the plants and irrigation is required to ensure optimum development and yields. Irrigation was carried out by means of microirrigation with microjets (microsprinkler technology), which distributed a uniform quantity of water in the well-watered replications. An ECRN-100 rainfall meter with a sensitivity of 0.2 mm was used for the accurate measurement of the quantity of water applied. In April 2013 (the shooting period) the plants suffered from rainfall deficiency and despite rainfall in May the soil moisture content dropped below the critical 22.3 tf% level (at depths of 30 and 60 cm) after the start of the dough development stage (Z83). This period was very hot, with maximum temperatures averaging 29.7 °C in the second and third weeks of June. Irrigation started in the booting stage (Z45) and was continued regularly till the end of July to keep the soil tension around −35 kPa. In 2014 the soil moisture content began to fall after flowering (Z65), partly due to the high temperatures recorded from May 20th onwards. In the rain-fed plots, values below 22.3 tf% were measured during grain filling, demonstrating that little water was available to the plants. Water supplies were provided on five occasions during the 2014 growing season, from the boot stage (Z45) to the early dough development stage (Z83). Air temperature (min, max, mean), humidity, global radiation, wind speed and rainfall quantity were recorded daily throughout the growing season.

4.3. Data Collection

The phenology of each accession was evaluated on the basis of the Zadoks scale [76]. The chlorophyll content of the flag leaf was measured using a SPAD-502 chlorophyll meter (Konica Minolta Sensing, Inc., Osaka, Japan) and expressed as a relative value (SPAD value) at the boot stage (SPAD45), at flowering (SPAD65), in the late stages of milk development (SPAD77), at early dough (SPAD83) and at the end of dough development (SPAD85). Readings were taken at three positions on the flag leaf—the base, the middle and the tip of the leaf lamina—in 24 replications per line for each water regime. The average per plot was computed and used in the analysis. The NDVI was measured in three phenophases: at the boot stage (NDVI45), at flowering (NDVI65) and at early dough development (NDVI83), using a GreenSeeker-505 handheld sensor (Trimble, Inc., Sunnyvale, CA, USA). Both the mean plant height, measured to the flag leaf collar (FLC), to the base of the ear (BE) and to the tip of the ear (TE) at the end of the physiological maturity phenophase, and the ear length (EL) were determined in 12 replications. The flag leaf area (FLA) at flowering (Z65) was determined as the mean value of the flag leaf area on the main tiller of ten representative plants (taken from positions where no plants were missing) using a LI-3100C leaf area meter (LI-COR, Inc., Lincoln, NE, USA). The relative water content (RWC%) of the flag leaf was calculated at flowering (Z65) on the flag leaves on the main tiller of three representative plants. The biomass was recorded in the flowering (Z65) and physiological maturity (Z91) phenophases, based on measurements of the whole aboveground fresh (WW65; WW91) and dry (DW65; DW91) biomass of ten representative plants. The ascorbate peroxidase and guaiacol peroxidase activity and the polyamine content (PUT, SPD, SPN) were measured in the flag leaves of the main tiller in five replications on samples collected from rain-fed and well-watered plots at flowering (Z65), as described by Pál et al. (2013) [77]. The plots were harvested with a small-plot combine, and data were recorded for the yield (GY) and yield components, including the thousand-kernel weight (TGW) and the average grain width (SW) and length (SL) on the main spike, determined using the Marvin digital seed analyzer system. Measurements were made on the spikelet number per spike for 16 main spikes (SKNM) per line, the grain number and grain weight per main spike (SNM), and the number of sterile apical (ASM) and basal (BSM) spikelets per main spike.

4.4. Statistical Analysis

Analysis of variance on the phenotypic data was performed for each individual year and for all the years combined using the linear random/mixed effect model in the Genstat 18 software package. Interactions and correlations between the traits and treatments were determined using principal component analysis, and means within a treatment were designated by different letters (p = 0.05) based on Duncan’s multiple range test. A modified equation was used to calculate the area under the SPAD value decline curve and the area under the vegetation index curve, which allowed genotypes that stayed green for a longer period to be selected [78].
A U S D C = i = 1 n 1 G i + 1 + G i 2   D i + 1 D i
Gi = SPAD value recorded on the ith day, Di = ith day, n = number of days.
A U V I C = i = 1 n 1 V i + 1 + V i 2   D i + 1 D i
Vi = NDVI value recorded on the ith day, Di = ith day, n = number of days.

5. Conclusions

In the present study, a complex field phenotyping system was described which, by a joint plant–environment–soil–genetics analysis, gave an accurate picture of the magnitude of the changes occurring in individual parameters and of the changes caused by drought stress in the morphology, physiology and yield biology of the plants. Lower plant density, which resulted in a larger growing area per plant, clearly had a significant positive effect on the grain weight per spike and the grain width. This may have been promoted by a number of indirect traits (higher chlorophyll content and NDVI value, longer stems, larger root system) and shows that genotypes with good tillering ability perform better at lower plant density. This was probably also aided by the better utilization of water, solar irradiation and nitrogen. The significantly higher grain number, grain weight and thousand-kernel weight in the main spike for the NIL++ lines confirmed the possibility that the higher yield of the QYld.idw-3B++ lines in spring sowings subjected to drought stress could be attributed to the positive effect of the Kofa QTL on the 3B chromosome. The better yield stability of lines with the KK2BLKK3BS allele combination indicates that it is worth giving further attention to the role of the QYld.idw-2B and QYld.idw-3B QTL regions in yield enhancement in durum wheat.

Supplementary Materials

Supplementary Materials can be found at https://0-www-mdpi-com.brum.beds.ac.uk/1422-0067/22/4/2053/s1.

Author Contributions

M.M. (Marco Maccaferri), J.B. and L.L. contributed to the conceptualization and designed all the research. J.B. analyzed and interpreted all the obtained data, prepared graphical presentations, and wrote the manuscript. J.B., M.M. (Marianna Mayer), V.T., M.C. and G.V. helped in material collection and participated in the analysis of data. J.B., M.M. (Marianna Mayer), V.T. and M.P. participated in experiments related to the polyamines measurements. K.M. and G.V. contributed to the editing and supervising. All authors have read and agreed to the published version of the manuscript.

Funding

Funding from the EU FP7-244374 DROPS and EU BONUS 12-1-2012-0017 projects is gratefully acknowledged.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

AUSDCarea under SPAD value decline curve
AUVICarea under vegetation index curve
HDheading date
SPADSPAD unit
NDVInormalized difference vegetation index
GYgrain yield (t/ha)
FWbiomass fresh weight (kg)
DWbiomass dry weight (kg)
HLtest weight (HL)
FLCplant height up to the flag leaf collar (cm)
BEplant height up to the base of the ear (cm)
TEplant height up to the tip of the ear (cm)
FTNfertile tiller number (pcs./m2)
SNMseed number/main spike
SWMseed weight/main spike (g)
SWseed width (cm)
SLseed length (cm)
TGWthousand-grain weight (g)
RWCrelative water content (%)
APXascorbate peroxidase (nkatal g-1 DW)
GPXguaiacol peroxidase (nkatal g-1 DW)
PUTputrescine (mg g-1 DW)
SPDspermidine (mg g-1 DW)
SPNspermine (mg g-1 DW)
ZZadoks scale

Appendix A

Figure A1. (a) Trends in soil moisture content, soil temperature and electrical conductivity at depths of 30 and 60 cm in the rain-fed plots (a) 2013 growing season; (b) 2014 growing season.
Figure A1. (a) Trends in soil moisture content, soil temperature and electrical conductivity at depths of 30 and 60 cm in the rain-fed plots (a) 2013 growing season; (b) 2014 growing season.
Ijms 22 02053 g0a1

References

  1. Trethowan, R.M.; Reynolds, M.; Sayre, K.; Ortiz-Monasterio, I. Adapting wheat cultivars to resource conserving farming practices and human nutritional needs. Ann. Appl. Biol. 2005, 146, 405–413. [Google Scholar] [CrossRef]
  2. Fukai, S.; Pantuwan, G.; Jongdee, B.; Cooper, M. Screening for drought resistance in rainfed lowland rice. Field Crops Res. 1999, 64, 61–74. [Google Scholar] [CrossRef]
  3. Fleury, D.; Jefferies, S.; Kuchel, H.; Langridge, P. Genetic and genomic tools to improve drought tolerance in wheat. J. Exp. Bot. 2010, 61, 3211–3222. [Google Scholar] [CrossRef] [Green Version]
  4. Grigorova, B.; Vaseva, I.; Demirevska, K.; Feller, U. Combined drought and heat stress in wheat: Changes in some heat shock proteins. Biol. Plant. 2011, 55, 105–111. [Google Scholar] [CrossRef]
  5. Keles, Y.; Oncel, I. Response of antioxidative defence system to temperature and water stress combinations in wheat seedlings. Plant Sci. 2002, 163, 783–790. [Google Scholar] [CrossRef]
  6. Barnabás, B.; Jäger, K.; Fehér, A. The effect of drought and heat stress on reproductive processes in cereals. Plant Cell Environ. 2008, 31, 11–38. [Google Scholar] [CrossRef]
  7. Rizhsky, L.; Liang, H.; Mittler, R. The combined effect of drought stress and heat shock on gene expression in tobacco. Plant Physiol. 2002, 130, 1143–1151. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Bartels, D.; Sunkar, R. Drought and salt tolerance in plants. Plant Sci. 2005, 24, 23–58. [Google Scholar] [CrossRef]
  9. Lobell, B.D.; Gourdji, M.S. The influence of climate change on global crop productivity. Plant Physiol. 2012, 160, 1686–1697. [Google Scholar] [CrossRef] [Green Version]
  10. Alaei, M.; Farboodi, M.; Khorshidi, M.B.; Zaeifizadeh, M. Durum wheat land race screening for drought tolerance. Middle-East J. Sci. Res. 2010, 63, 289–292. [Google Scholar]
  11. Dixit, S.; Singh, A.; Kumar, A. Rice breeding for high grain yield under drought: A strategic solution to a complex problem. Int. J. Agron. 2014. [Google Scholar] [CrossRef]
  12. Kilic, H.; Yagbasanlar, T. The effect of drought stress on grain yield, yield components and some quality traits of durum wheat (Triticum turgidum ssp. durum) cultivars. Not. Bot. Horti. Agrobot. Cluj-Napoca 2010, 38, 164–170. [Google Scholar]
  13. Nawaz, H.; Hussain, N.; Yasmeen, A. Growth, yield and antioxidants status of wheat (Triticum aestivum L.) cultivars under water deficit conditions. Pak. J. Agric. Sci. 2015, 52, 953–959. [Google Scholar]
  14. Rucker, K.S.; Kevin, C.K.; Holbrook, C.C.; Hook, J.E. Identification of peanut genotypes with improved drought avoidance traits. Peanut Sci. 1995, 22, 14–18. [Google Scholar] [CrossRef] [Green Version]
  15. Shao, H.B.; Chu, L.Y.; Jaleel, C.A.; Zhao, C.X. Water-deficit stress-induced anatomical changes in higher plants. Comptes Rendus 2008, 331, 215–225. [Google Scholar] [CrossRef] [PubMed]
  16. Tardieau, F.; Reymond, M.; Hamard, P.; Granier, C.; Muller, B. Spatial distributions of expansion rate, cell division rate and cell size in maize leaves: A synthesis of the effects of soil water status, evaporative demand and temperature. J. Exp. Bot. 2000, 51, 1505–1514. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Xiong, L.; Wang, R.-G.; Mao, G.; Koczan, J.M. Identification of drought tolerance determinants by genetic analysis of root response to drought stress and abscisic acid. Plant Physiol. 2006, 142, 1065–1074. [Google Scholar] [CrossRef] [Green Version]
  18. Sharp, R.E.; Wu, Y.J.; Voetberg, G.S.; Saab, I.N.; Lenoble, M.E. Confirmation that abscisic-acid accumulation is required for maize primary root elongation at low water potentials. J. Exp. Bot. 1994, 45, 1743–1751. [Google Scholar] [CrossRef]
  19. Czyczyło-Mysza, I.; Marcinska, I.; Skrzypek, E.; Chrupek, M.; Grzesiak, S.; Hura, T.; Stojałowski, S.; Myskow, B.; Milczarski, P.; Quarrie, S. Mapping QTLs for yield components and chlorophyll a fluorescence parameters in wheat under three levels of water availability. Plant Genet. Resour. C. 2011, 9, 291–295. [Google Scholar] [CrossRef]
  20. Chachar, M.H.; Chachar, N.A.; Chachar, Q.; Mujtaba, S.M.; Chachar, S.; Chachar, Z. Physiological characterization of six wheat genotypes for drought tolerance. Int. J. Res. Granthaalayah 2016, 4, 184–196. [Google Scholar]
  21. Srivalli, S.; Khanna-Chopra, R. Delayed wheat flag leaf senescence due to removal of spikelets is associated with increased activities of leaf antioxidant enzymes, reduced glutathione/oxidized glutathione ratio and oxidative damage to mitochondrial proteins. Plant Physiol. Biochem. 2009, 47, 663–670. [Google Scholar] [CrossRef]
  22. Zhu, C.; Zhu, J.; Zeng, Q.; Liu, G.; Xie, Z.; Tang, H.; Cao, J.; Zhao, X. Elevated CO2 accelerates flag leaf senescence in wheat due to ear photosynthesis which causes greater ear nitrogen sink capacity and ear carbon sink limitation. Funct. Plant Biol. 2009, 36, 291–299. [Google Scholar] [CrossRef]
  23. Gregersen, P.L.; Holm, P.B.; Krupinska, K. Leaf senescence and nutrient remobilization in barley and wheat. Plant Biol. 2008, 10, 37–49. [Google Scholar] [CrossRef]
  24. Falqueto, A.R.; Cassol, D.; De Magalhães, J.M.A.; De Oliveira, A.C.; Bacarin, M.A. Physiological analysis of leaf senescence of two rice cultivars with different yield potential. Pesqui. Agropecu. Bras. 2009, 44, 695–700. [Google Scholar] [CrossRef]
  25. Saeedipour, S. Effect of drought at the post-anthesis stage on remobilization of carbon reserves in two wheat cultivars differing in senescence properties. Int. J. Plant Physiol. Biochem. 2011, 3, 15–24. [Google Scholar]
  26. Yang, J.; Zhang, J.; Huang, Z.; Zhu, Q.; Wang, L. Remobilization of carbon reserves is improved by controlled soil-drying during grain filling of wheat. Crop Sci. 2000, 40, 1645–1655. [Google Scholar] [CrossRef]
  27. Yang, J.; Zhang, J.; Wang, Z.; Zhu, Q.; Liu, L. Water deficit-induced senescence and its relationship to the remobilization of pre-stored carbon in wheat in grain filling. Agron. J. 2001, 93, 196–206. [Google Scholar] [CrossRef]
  28. Waters, B.M.; Uauy, C.; Dubcovsky, J.; Grusak, M.A. Wheat (Triticum aestivum) NAM proteins regulate the translocation of iron, zinc, and nitrogen compounds from vegetative tissues to grain. J. Exp. Bot. 2009, 60, 4263–4274. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Sedigheh, H.G.; Mortazavian, M.; Norouzian, D.; Atyabil, M.; Akbarzadeh, A.; Hasanpoor, K.; Ghorbani, M. Oxidative stress and leaf senescence. BMC Res. Notes 2011, 4, 477. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Ishida, H.; Makino, A.; Mae, T. Fragmentation of the large subunit of ribulose-1,5-bisphosphate carboxylase by reactive oxygen species occurs near Gly-329. J. Biol. Chem. 1999, 274, 5222–5226. [Google Scholar] [CrossRef] [Green Version]
  31. Bray, E.A. Classification of genes differentially expressed during water-deficit stress in Arabidopsis thaliana: An analysis using microarray and differential expression data. Ann. Bot. 2002, 89, 803–811. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Habash, D.Z.; Kehel, Z.; Nachit, M. Genomic approaches for designing durum wheat ready for climate change with a focus on drought. J. Exp. Bot. 2009, 60, 2805–2815. [Google Scholar] [CrossRef] [PubMed]
  33. Barnabás, B.; Fehér, A. Az embriófejlődés és a szemfeltöltődés molekuláris folyamatai a vízhiányban. In A búza Nemesbítésének Tudománya; Dudits, D., Ed.; MTA Szegedi Biológiai Központ-Winter Fair Kft: Szeged, Hungary, 2006; pp. 180–190. [Google Scholar]
  34. Jäger, K.; Fábián, A.; Barnabás, B. A tartós vízhiány és a magas hőmérséklet hatása a búza pollen életképességére. Martonvásár 2013, 55, 24–25. [Google Scholar]
  35. Gao, X.; Francis, D.; Ormrod, J.C.; Bennett, M.D. Changes in cell number and cell division activity during endosperm development in allohexaploid wheat, Triticum aestivum L. J. Exp. Bot. 1992, 43, 1603–1609. [Google Scholar] [CrossRef]
  36. Yang, J.C.; Zhang, J.H.; Wang, Z.Q.; Zhu, Q.S.; Liu, L. Involvement of abscisic acid and cytokinins in the senescence and remobilization of carbon reserves in wheat subjected to water stress during grain filling. Plant Cell Environ. 2003, 26, 1621–1631. [Google Scholar] [CrossRef]
  37. Jäger, K.; Fábián, A.; Barnabás, B. Effect of water deficit and elevate temperature on pollen development of drought sensitive and tolerant winter wheat (Triticum aestivum L.) genotypes. Acta Biol. Szeged. 2008, 52, 67–71. [Google Scholar]
  38. Bányai, J.; Maccaferri, M.; Cané, M.A.; Monostori, I.; Spitkó, T.; Kuti, C.; Mészáros, K.; Láng, L.; Pál, M.; Karsai, I. Phenotypical and physiological study of near-isogenic durum wheat lines under contrasting water regimes. S. Afr. J. Bot. 2017, 108, 248–255. [Google Scholar] [CrossRef]
  39. Varshney, R.K.; Graner, A.; Sorrells, M.E. Genomics-assisted breeding for crop improvement. Trends Plant Sci. 2005, 10, 621–630. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Kooke, R.; Wijnker, E.; Keurentjes, J.J.B. Backcross populations and near isogenic lines. In Quantitative Trait Loci (QTL): Methods and Protocols; Methods in Molecular Biology; Rifkin, S.A., Ed.; Humana Press: New York, NY, USA, 2012; Volume 871, pp. 1–14. [Google Scholar]
  41. Sevanthi, A.M.; Prakash, C.; Shanmugavadivel, P.S. Recent progress in rice varietal development for abiotic stress tolerance. In Advances in Rice Research for Abiotic Stress Tolerance; Hasanuzzaman, M., Ed.; Woodhead Publishing: Cambridge, UK, 2019; pp. 47–68. [Google Scholar]
  42. Tsujimoto, H. Production of near-Isogenic lines and marked monosomic lines in common wheat (Triticum aestivum) cv. Chinese Spring. Heredity 2001, 92, 254–259. [Google Scholar] [CrossRef] [Green Version]
  43. Mia, M.S.; Liu, H.; Wang, X.; Yan, G. Multiple near-isogenic lines targeting a QTL hotspot of drought tolerance showed contrasting performance under post-anthesis water stress. Front. Plant Sci. 2019, 10, 271. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Naser, M.A.; Khosla, R.; Longchamps, L.; Dahal, S. Using NDVI to differentiate wheat genotypes productivity under dryland and irrigated conditions. Remote Sens. 2020, 12, 824. [Google Scholar] [CrossRef] [Green Version]
  45. Fu, Z.; Jiang, J.; Gao, Y.; Krienke, B.; Wang, M.; Zhong, K.; Cao, Q.; Tian, Y.; Zhu, Y.; Cao, W.; et al. Wheat growth monitoring and yield estimation based on multi-rotor unmanned aerial vehicle. Remote Sens. 2020, 12, 508. [Google Scholar] [CrossRef] [Green Version]
  46. Goodwin, A.W.; Lindsey, L.E.; Harrison, S.K.; Paul, P.A. Estimating wheat yield with normalized difference vegetation index and fractional green canopy cover. Crop Forage Turfgrass Manag. 2018, 4, 180026. [Google Scholar] [CrossRef] [Green Version]
  47. Vannoppen, A.; Gobin, A.; Kotova, L.; Top, S.; De Cruz, L.; Viksna, A.; Aniskevich, S.; Bobylev, L.; Buntemeyer, L.; Caluwaerts, S.; et al. Wheat yield estimation from NDVI and regional climate models in Latvia. Remote Sens. 2020, 12, 2206. [Google Scholar] [CrossRef]
  48. Dodig, D.; Savic, J.; Kandic, V.; Zoric, M.; Radovic, B.V.; Popovic, A.; Quarrie, S. Responses of wheat plants under post-anthesis stress induced by defoliation: I. Contribution of agro-physiological traits to grain yield. Exp. Agric. 2015. [Google Scholar] [CrossRef] [Green Version]
  49. Cossani, C.M.; Slafer, G.A.; Savin, R. Yield and biomass in wheat and barley under a range of conditions in a Mediterranean site. Field Crops Res. 2009, 112, 205–2013. [Google Scholar] [CrossRef]
  50. Cruz-Aguado, J.A.; Rodés, R.; Pérez, J.P.; Dorado, M. Morphological characteristics and yield components associated with accumulation and loss of dry matter in internodes of wheat. Field Crops Res. 2000, 66, 129–139. [Google Scholar] [CrossRef]
  51. Ehdaie, B.; Alloush, G.A.; Madore, M.A.; Waines, J.G. Genotypic variation for stem reserves and mobilization in wheat. I. Post anthesis changes in internode dry matter. Crop Sci. 2006, 46, 735–746. [Google Scholar] [CrossRef]
  52. Majer, P.; Sass, L.; Lelley, T.; Cseuz, L.; Vass, I.; Dudits, D.; Pauk, J. Testing drought tolerance of wheat by a complex stress diagnostic system installed in greenhouse. Acta Biol. Szeged. 2008, 52, 97–100. [Google Scholar]
  53. Korohou, T.; Okinda, C.; Li, H.; Cao, Y.; Nyalala, I.; Huo, L.; Potcho, M.; Li, X.; Ding, Q. Wheat grain yield estimation based on image morphological properties and wheat biomass. J. Sens. 2020, 2020, 1571936. [Google Scholar] [CrossRef]
  54. Bennett, D.; Reynolds, M.; Mullan, D.; Izanloo, A.; Kuchel, H.; Langridge, P.; Schnurbusch, T. Detection of two major grain yield QTL in bread wheat (Triticum aestivum L.) under heat, drought and high yield potential environments. Theor. Appl. Genet. 2012, 125, 1473–1485. [Google Scholar] [CrossRef]
  55. Zhang, L.Y.; Liu, D.C.; Guo, X.L.; Yang, W.L.; Sun, J.Z.; Wang, D.W.; Zhang, A. Genomic distribution of quantitative trait loci for yield and yield-related traits in common wheat. J. Integr. Plant Biol. 2010, 52, 996–1007. [Google Scholar] [CrossRef] [PubMed]
  56. Pinto, R.S.; Reynolds, M.P.; Mathews, K.L.; Mcintyre, C.L.; Olivares–Villegas, J.J.; Chapman, S.C. Heat and drought adaptive QTL in a wheat population designed to minimize confounding agronomic effects. Theor. Appl. Genet. 2010, 121, 1001–1021. [Google Scholar] [CrossRef] [Green Version]
  57. Bonneau, J.; Taylor, J.; Parent, B.; Bennett, D.; Reynolds, M.; Feuillet, C.; Langridge, P.; Mather, D. Multi-environment analysis and improved mapping of a yield-related QTL on chromosome 3B of wheat. Theor. Appl. Genet. 2013, 126, 747–761. [Google Scholar] [CrossRef]
  58. Shukla, S.; Mishra, B.K.; Mishra, R.; Siddiqui, A.; Pandey, R.; Rastogi, A. Comparative study for stability and adaptability through different models in developed high thebaine lines of opium poppy (Papaver somniferum L.). Ind. Crops Prod. 2015, 74, 875–886. [Google Scholar] [CrossRef]
  59. Dudziak, K.; Zapalska, M.; Börner, A.; Szczerba, H.; Kowalczyk, K.; Nowak, M. Analysis of wheat gene expression related to the oxidative stress response and signal transduction under short-term osmotic stress. Sci. Rep. 2019, 9, 2743. [Google Scholar] [CrossRef] [Green Version]
  60. Li, Y.; Cui, Z.; Ni, Y.; Zheng, M.; Yang, D.; Jin, M.; Chen, J.; Wang, Z.; Yin, Y. Plant density effect on grain number and weight of two winter wheat cultivars at different spikelet and grain positions. PLoS ONE 2016. [Google Scholar] [CrossRef]
  61. Whaley, Y.M.; Sparkes, D.L.; Foulkes, M.J.; Spink, J.H.; Semere, T.; Scott, R.K. The physiological response of winter wheat to reductions in plant density. Ann. Appl. Biol. 2000, 137, 165–177. [Google Scholar] [CrossRef]
  62. Wood, G.A.; Welsh, J.P.; Godwin, R.J.; Taylor, J.C.; Earl, R.; Knight, S.M. Real-time measures of canopy size as a basis for spatially varying nitrogen applications to winter wheat sown at different seed rates. Biosyst. Eng. 2003, 84, 513–531. [Google Scholar] [CrossRef] [Green Version]
  63. Valério, I.P.; De Carvalho, F.I.F.; Benin, G.; Da Silveira, G.; Da Silva, J.A.G.; Nornberg, R.; Hagemann, T.; De Souza Luche, H.; De Oliveira, C.A. Seeding density in wheat: The more, the merrier? Sci. Agric. 2013, 40, 176–184. [Google Scholar] [CrossRef]
  64. Ozturk, A.; Caglar, O.; Bulut, S. Growth and yield response of facultative wheat to winter sowing, freezing sowing and spring sowing at different seeding rates. J. Agron. Crop Sci. 2006, 192, 10–16. [Google Scholar] [CrossRef]
  65. Bastos, L.M.; Carciochi, W.; Lollato, R.P.; Jaenisch, B.R.; Rezende, C.R.; Schwalbert, R.; Prasad, P.V.V.; Zhang, G.; Fritz, A.K.; Foster, C.; et al. Winter wheat yield response to plant density as a function of yield environment and tillering potential: A review and field studies. Front. Plant Sci. 2020, 11, 54. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Laxa, M.; Liebthal, M.; Telman, W.; Chibani, K.; Dietz, K.J. The role of the plant antioxidant system in drought tolerance. Antioxidants 2019, 8, 94. [Google Scholar] [CrossRef] [Green Version]
  67. Pour-Benab, S.M.; Fabriki-Ourang, S.; Mehrabi, A.A. Expression of dehydrin and antioxidant genes and enzymatic antioxidant defense under drought stress in wild relatives of wheat. Biotechnol. Equip. 2019, 33, 1063–1073. [Google Scholar] [CrossRef] [Green Version]
  68. Salekjalali, M.; Haddad, R.; Jafari, B. Effects of soil water shortages on the activity of antioxidant enzymes and the contents of chlorophylls and proteins in barley. Am. Eurasian J. Agric. Environ. Sci. 2012, 12, 57–63. [Google Scholar]
  69. Huseynova, I.M.; Suleymanov, S.Y.; Rustamova, S.M. Response of photosynthetic apparatus and antioxidant defense systems in Triticum aestivum L. genotypes subjected to drought stress. Proc. ANAS (Biol. Sci.) 2010, 65, 49–59. [Google Scholar]
  70. Chen, D.; Shao, Q.; Yin, L.; Younis, A.; Zheng, B. Polyamine function in plants: Metabolism, regulation on development, and roles in abiotic stress responses. Front. Plant Sci. 2019, 9, 1945. [Google Scholar] [CrossRef] [PubMed]
  71. Aldesuquy, H.; Haroun, S.; Abo-Hamed, S.; El Saied, A.-W. Involvement of spermine and spermidine in the control of productivity and biochemical aspects of yielded grains of wheat plants irrigated with waste water. Egypt. J. Basic Appl. Sci. 2014, 1, 16–28. [Google Scholar] [CrossRef] [Green Version]
  72. Liu, Y.; Liang, H.; Lv, X.; Liu, D.; Wen, X.; Liao, Y. Effect of polyamines on the grain filling of wheat under drought stress. Plant Physiol. Biochem. 2016, 100, 113–129. [Google Scholar] [CrossRef]
  73. Sahitya, U.L.; Krishna, M.S.R.; Deepthi, R.S.; Prasad, G.S.; Kasim, D.P. Seed antioxidants interplay with drought stress tolerance indices in chilli (Capsicum annuum L) seedlings. Biomed. Res. Int. 2018, 2018, 1605096. [Google Scholar] [CrossRef] [Green Version]
  74. Hasanuzzaman, M.; Bhuyan, M.H.M.B.; Zulfiqar, F.; Raza, A.; Mohsin, S.M.; Al Mahmud, J.; Fujita, M.; Fotopoulos, V. Reactive oxygen species and antioxidant defense in plants under abiotic stress: Revisiting the crucial role of a universal defense regulator. Antioxidants 2020, 9, 681. [Google Scholar] [CrossRef]
  75. Maccaferri, M.; Sanguineti, M.C.; Corneti, S.; Ortega, J.L.; Salem, M.B.; Bort, J.; DeAmbrogio, E.; del Moral, L.F.; Demontis, A.; El-Ahmed, A.; et al. Quantitative Trait Loci for grain yield and adaptation of durum wheat (Triticum durum Desf.) across a wide range of water availability. Genetics 2008, 178, 489–511. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Zadoks, J.C.; Chang, T.T.; Konzak, C.F. A decimal code for the growth stage of cereals. Weed Res. 1974, 14, 415–421. [Google Scholar] [CrossRef]
  77. Pál, M.; Kovács, V.; Vida, G.; Szalai, G.; Janda, T. Changes induced by powdery mildew in the salicylic acid and polyamine contents and the antioxidant enzyme activities of wheat lines. Eur. J. Plant Pathol. 2013, 135, 35–47. [Google Scholar] [CrossRef]
  78. Rosyara, U.R.; Pant, K.; Duveiller, E.; Sharma, R.C. Variation in chlorophyll content, anatomical traits and agronomic performance of wheat genotypes differing in spot blotch resistance under natural epiphytotic conditions. Australas. Plant Pathol. 2007, 36, 245–251. [Google Scholar] [CrossRef]
Figure 1. Changes in the (a) ascorbate peroxidase (APX) and (b) guaiacol peroxidase (G-POD) enzyme activity and (c) in the putrescine (PUT) content in the flag leaf in the well-watered (WW) and rain-fed (RF) treatments. *, ** and *** denote significant differences between the treatments at the p < 0.05, 0.01 and 0.001 levels of probability, respectively (2014).
Figure 1. Changes in the (a) ascorbate peroxidase (APX) and (b) guaiacol peroxidase (G-POD) enzyme activity and (c) in the putrescine (PUT) content in the flag leaf in the well-watered (WW) and rain-fed (RF) treatments. *, ** and *** denote significant differences between the treatments at the p < 0.05, 0.01 and 0.001 levels of probability, respectively (2014).
Ijms 22 02053 g001
Figure 2. Correlation system demonstrated by principal component analysis between plant development, morphological, physiological and biochemical traits and the yield and yield components based on the parameter matrix formed from the mean values of rain-fed treatments over two years.
Figure 2. Correlation system demonstrated by principal component analysis between plant development, morphological, physiological and biochemical traits and the yield and yield components based on the parameter matrix formed from the mean values of rain-fed treatments over two years.
Ijms 22 02053 g002
Figure 3. Pearson’s correlation coefficients for the plant development, morphological, physiological and biochemical traits and the yield components of eight near-isogenic lines of spring durum wheat in the rain-fed treatment (2013–2014). Critical r values for the correlation coefficients: 0.6319 (p < 0.05); 0.7646 (p < 0.01); 0.8721 (p < 0.001).
Figure 3. Pearson’s correlation coefficients for the plant development, morphological, physiological and biochemical traits and the yield components of eight near-isogenic lines of spring durum wheat in the rain-fed treatment (2013–2014). Critical r values for the correlation coefficients: 0.6319 (p < 0.05); 0.7646 (p < 0.01); 0.8721 (p < 0.001).
Ijms 22 02053 g003
Table 1. Significance of year, genotype, treatment and sowing density effects and interactions on plant development, morphological and physiological parameters and on yield components, based on the Chi-squared test for eight near-isogenic lines of durum wheat. In the REML analysis, F = Wald statistic/d.f., combined over the two years. *, **, *** significant at the p < 0.05, 0.01 and 0.001 probability levels, respectively.
Table 1. Significance of year, genotype, treatment and sowing density effects and interactions on plant development, morphological and physiological parameters and on yield components, based on the Chi-squared test for eight near-isogenic lines of durum wheat. In the REML analysis, F = Wald statistic/d.f., combined over the two years. *, **, *** significant at the p < 0.05, 0.01 and 0.001 probability levels, respectively.
Year (Y)Genotype (G)Treatment (T)Density (D)G × YG × TG × DY × TY × DT × DY × G × TY × G × DY × T × DG × T × DY × G × T × D
d.f.171177711177177
HD 12156.85 ***502.57 ***3.71525.87 ***27.06 ***2.026.53 ***6.77 *0.0812.71 ***3.14*5.26 ***1.821.270.48
FD 23167.93 ***332.91 ***1.88393.75 ***20.86 ***1.243.97 ***8.72 **6.44 *9.22 **2.53 *3.78 ***1.881.431.12
MD 33798.87 ***157.56 ***377.06 ***119.30 ***25.58 ***4.84 ***3.81 ***21.11 ***27.96 ***20.32 ***3.21 **4.53 ***26.16 ***1.561.03
GY 472.11 ***30.10 ***247.04 ***109.42 ***8.23 ***1.582.62 *78.62 ***32.88 ***37.09 ***2.24 *1.2433.14 ***2.12 *2.75 **
TGW 5133.71 ***88.35 ***685.52 ***56.52 ***20.26 ***9.06 ***2.87 **487.73 ***0.080.522.46 **0.831.952.81 **3.02 **
SNM 6826.46 ***8.50 ***2.020.081.270.271.090.032.840.940.790.160.951.060.79
SWM 7132.34 ***27.68 ***414.71 ***133.16 ***1.761.421.730.713.460.191.030.9524.78 ***0.510.48
SKNM 8644.12 ***53.33 ***7.9724.600.4519.62 **3.140.380.041.391.430.100.550.98
BSM 91.340.580.280.070.690.400.430.110.260.510.910.460.150.910.28
ASM 100.370.591.1700.890.770.681.660.030.070.560.830.030.850.75
SL 114.06 *13.07 ***0.6831.23 ***4.09 ***1.901.17026.32 ***15.52 ***4.66 ***3.50 *4.06 *0.761.99
SW 12600.45 ***220.04 ***34.51 ***16.23 ***44.47 ***1.777.44 ***600.45 ***30.68 ***28.85 ***3.141.708.14 ***1.380.99
FLC 13876.31 ***18.48 ***16.41 ***13.94 ***12.82 ***1.771.1627.49 ***4.27 *31.21 ***1.690.740.11.340.74
BE 14978.21 ***128.18 ***52.24 ***16.60 ***32.02 ***3.31 *1.2012.97 ***0.5841.13 ***2.88 *2.002.171.211.66
TE 15979.14 ***102.85 ***62.92 ***22.70 ***35.24 ***0.280.417.33 **1.2622.70 ***1.180.551.312.021.37
EL 16101.32 ***4.38 ***25.08 ***8.56 ***0.720.430.410.420.020.890.130.150.020.170.07
FTN 170.419.97 ***0.357.77 **8.19 ***2.27 *1.12 *0.291.680.111.760.380.151.040.26
FLA 1828.4571.7202.0117.073.212.280.460.150.062.610.700.211.090.82
SPAD45 191610.61 ***18.52 ***236.93 ***1.861.751.210.250.2011.43 ***1.061.041.230.871.27
SPAD65 20220.84 ***16.70 ***07.103.141.020.660.111.111.830.651.061.401.440.46
SPAD77 21202.70 ***26.52 ***70.18 ***28.00 *10.75 ***2.29 *0.680.7352.70 ***3.251.910.183.791.141.38
SPAD83 2214.41 ***11.60 ***336.60 ***19.00 *5.23 ***3.15 **0.351.5155.48 ***18.87 ***1.200.2300.460.49
SPAD85 232.633.7023.0719.43 ***3.66 ***1.201.0311.20 ***27.94 ***0.561.400.6210.69 **0.210.96
NDVI45 241.0516.91 ***0.020.023.82 ***0.870.150.130.200.010.670.490.060.330.64
NDVI65 251.182.13 *0.020.041.240.090.060.020.070.010.150.070.010.080.22
NDVI83 2630.83 ***7.45 ***5.54 *56.25 ***3.86 ***0.480.880.3537.11 ***0.710.280.150.160.490.49
RWC 27858.81 ***0.9245.38 ***1.471.050.480.360.560.530.080.411.473.92 *0.560.68
WW65 281.323.60 **0.180.261.740.200.240.110.020.020.600.120.010.240.36
DW65 29407.65 ***1.8449.66 ***34.36 ***1.161.260.994.041.434.601.230.490.190.870.63
WW91 30202.57 ***3.5047.54 ***26.68 ***1.380.420.7645.45 ***6.93 *7.22 *0.730.543.261.190.08
DW91 31261.27 ***3.3539.82 ***17.72 ***1.480.340.6540.35 ***2.946.07 *0.610.754.951.080.12
1 Days to heading; 2 days to flowering; 3 days to maturity; 4 grain yield; 5 thousand-grain weight; 6 seed number/main spike; 7 seed weight/main spike; 8 spikelet number/main spike; 9 basal sterile spikelet number/main spike; 10 apical sterile spikelet number/main spike; 11 seed length; 12 seed width; 13 plant height up to the flag leaf collar; 14 plant height up to the base of the ear; 15 plant height up to the tip of the ear; 16 ear length; 17 fertile tiller number; 18 flag leaf area; 19 SPAD value at Z45; 20 SPAD value at Z65; 21 SPAD value at Z77; 22 SPAD value at Z83; 23 SPAD value at Z85; 24 normalized difference vegetation index (NDVI) at Z45; 25 NDVI at Z65; 26 NDVI at Z83; 27 relative water content; 28 fresh biomass at Z65; 29 dry biomass at Z65; 30 fresh biomass at Z91; 31 dry biomass at Z91.
Table 2. Mean values of plant development, morphological and physiological traits and of yield components recorded for QYld.idw-3B−− and QYld.idw-3B++ near-isogenic lines of spring durum wheat in rain-fed and well-watered treatments. *, ** and *** denote significant differences within treatments between the QYld.idw-3B−− and QYld.idw-3B++ near-isogenic lines at the p < 0.05, 0.01 and 0.001 levels of probability, respectively (2013).
Table 2. Mean values of plant development, morphological and physiological traits and of yield components recorded for QYld.idw-3B−− and QYld.idw-3B++ near-isogenic lines of spring durum wheat in rain-fed and well-watered treatments. *, ** and *** denote significant differences within treatments between the QYld.idw-3B−− and QYld.idw-3B++ near-isogenic lines at the p < 0.05, 0.01 and 0.001 levels of probability, respectively (2013).
HD 1FD 2MD 3SPAD45 4SPAD65 5SPAD77 6SPAD83 7SPAD85 8NDVI45 9NDVI65 10NDVI83 11AUSDC 12AUVIC 13FLA 14
WNWWNWWNWWNWWNWWNWWNWWNWWNWWNWWNWWNWWNWWNW
NIL1−−150.56 ***150.83156 ***156.36194.72 **193.2742.1041.1648.3048.4837.4130.7829.3221.2110.929.800.34980.35460.55150.55020.33970.32791597.951441.5314.8615.3127.65 **26.22
NIL++149.09148.94154.44154.39194.36193.1442.9742.9148.8548.8943.8141.5734.40 ***26.95 **10.7210.730.35850.35630.55150.55320.3900 ***0.3800*1839.04 ***1734.17 **15.5715.6523.5524.48
LSD5%0.452.270.662.150.191.911.063.061.372.782.413.381.443.432.682.100.00820.02470.00440.01150.01790.042972.74175.270.451.692.523.67
LSD1%0.673.360.973.180.282.821.584.532.044.113.585.012.145.083.963.120.012230.03660.00650.01700.02650.0635107.65259.400.672.513.745.44
LSD0.1%1.045.191.504.910.434.372.447.013.156.365.537.743.307.856.134.820.01890.05660.01010.02630.04090.0982166.36400.861.033.875.788.41
FTN 15FLC 16BE 17TE 18EL 19GY 20SKNM 21SNM 22SWM 23BSM 24ASM 25TGW 26SL 27SW 28
WNWWNWWNWWNWWNWWNWWNWWNWWNWWNWWNWWNWWNWWNW
NIL1−−370.59378.3440.5940.9055.6854.1461.7959.855.935.653.022.6713.5413.5135.2534.831.881.440.75820.85860.54330.541935.4343.637.016.992.893.06
NIL++400.63 **392.30 **43.13 **43.52 *62.66 **62.05 ***69.33 ***67.82 ***6.215.813.42 ***3.15 ***13.4213.3237.52 **37.15 ***2.23 **1.64 **0.61760.68610.54330.507039.82 ***46.50 **7.017.023.05 ***3.16 *
LSD5%18.596.601.152.363.162.862.912.690.100.220.150.200.220.331.020.980.150.130.49920.33130.18890.16081.121.700.030.090.050.09
LSD1%27.529.771.713.504.674.234.313.980.160.330.230.300.320.481.511.450.230.190.73880.49020.27960.23801.662.520.050.140.070.24
LSD0.1%42.5315.112.645.417.236.546.676.150.240.520.350.470.500.752.332.240.350.301.14170.75760.43210.36782.573.900.080.220.110.37
HL 29P 30G 31Z 32WW65 33WW91 34DW65 35DW91 36RWC 37APX 38GPX 39PUT 40SPD 41SPN 42
WNWWNWWNWWNWWNWWNWWNWWNWWNWWNWWNWWNWWNWWNW
NIL1−−76.7482.5417.31 ***16.7636.18 **35.0962.0964.110.14440.14390.03400.03390.01950.01400.02400.016191.8087.62824.89576.573040.494899.19192.87355.05394.47524.73217.78296.51
NIL++78.9382.5716.1617.0734.0135.3961.3465.910.14370.14420.0389 *0.03800.01860.01610.0281*0.028091.0787.52982.15504.773438.564840.85177.70460.20*419.43572.83305.50342.78
LSD5%1.541.360.451.391.172.361.012.870.00390.00540.00450.00940.00140.00390.00340.00691.752.88282.9078.581331.833291.5786.8891.8924.13171.95113.3859.48
LSD1%2.282.010.672.061.733.501.494.250.00580.00800.00670.01390.00210.00580.00510.01022.594.26519.22144.232444.386041.20159.46168.6644.30315.60208.10109.05
LSD0.1%3.523.111.043.192.675.412.316.560.00900.01240.01040.02150.00330.00900.00780.01584.016.581148.88319.145408.6113367.16352.84373.2098.03698.31460.47241.29
1 Days to heading; 2 days to flowering; 3 days to maturity; 4 SPAD value at Z45; 5 SPAD value at Z65; 6 SPAD value at Z77; 7 SPAD value at Z83; 8 SPAD value at Z85; 9 NDVI at Z45; 10 NDVI at Z65; 11 NDVI at Z83; 12 area under SPAD value decline curve; 13 area under vegetation index curve; 14 flag leaf area; 15 fertile tiller number; 16 plant height up to the flag leaf collar; 17 plant height up to the base of the ear; 18 plant height up to the tip of the ear; 19 ear length; 20 grain yield; 21 spikelet number/main spike; 22 seed number/main spike; 23 seed weight/main spike; 24 basal sterile spikelet number/main spike; 25 apical sterile spikelet number/main spike; 26 thousand-grain weight; 27 seed length; 28 seed width; 29 test weight; 30 protein content; 31 glutenin content; 32 Zeleny index; 33 fresh biomass at Z65; 34 fresh biomass at Z91; 35 dry biomass at Z65; 36 dry biomass at Z91; 37 relative water content; 38 ascorbate peroxidase; 39 guaiacol peroxidase; 40 putrescine; 41 spermidine; 42 days to spermine.
Table 3. Mean values of plant development, morphological and physiological traits and of yield components recorded for QYld.idw-3B−− and QYld.idw-3B++ near-isogenic lines of spring durum wheat in rain-fed and well-watered treatments. *, ** and *** denote significant differences within treatments between the QYld.idw-3B−− and QYld.idw-3B++ near-isogenic lines at the p < 0.05, 0.01 and 0.001 levels of probability, respectively (2014).
Table 3. Mean values of plant development, morphological and physiological traits and of yield components recorded for QYld.idw-3B−− and QYld.idw-3B++ near-isogenic lines of spring durum wheat in rain-fed and well-watered treatments. *, ** and *** denote significant differences within treatments between the QYld.idw-3B−− and QYld.idw-3B++ near-isogenic lines at the p < 0.05, 0.01 and 0.001 levels of probability, respectively (2014).
HD 1FD 2MD 3SPAD45 4SPAD65 5SPAD77 6SPAD83 7SPAD85 8NDVI45 9NDVI65 10NDVI83 11AUSDC 12AUVIC 13
WNWWNWWNWWNWWNWWNWWNWWNWWNWWNWWNWWNWWNW
NIL1−−146.50146.09150.43150.00189.09186.4252.2152.2348.3051.87187546.0341.9229.4521.4912.139.190.35270.35220.55740.55570.38290.37002170.422020.2824.8922.14
NIL++146.05145.71150149.62189.70188.0453.6753.1848.8554.20 *47.9544.9431.4724.81 **13.3110.030.35270.35230.55570.55560.4000 **0.3900 *2229.152109.61 **26.28 *22.54
LSD5%1.231.371.121.311.151.841.141.141.371.632.012.092.381.661.250.960.00320.00240.00470.00390.00870.015698.8145.390.960.86
LSD1%1.822.031.661.951.712.721.691.682.042.422.973.103.522.471.851.420.00480.00360.00700.00590.01290.0231146.2367.181.411.28
LSD0.1%2.813.142.563.012.644.212.612.613.153.744.594.795.443.812.862.190.00740.00550.01080.00910.01990.0357225.98103.822.191.98
FLA 14FTN 15FLC 16BE 17TE 18EL 19GY 20SKNM 21SNM 22SWM 23BSM 24ASM 25TGW 26
WNWWNWWNWWNWWNWWNWWNWWNWWNWWNWWNWWNWWNW
NIL1−−26.4127.37361.31362.7851.4348.3365.5763.9773.5970.556.736.313.592.5815.2514.9844.7144.262.111.710.61250.51460.48170.833037.6738.83
NIL++28.0227.49423.81 **407.77 ***52.8450.47 *73.50 **68.67 *80.00 **76.02 **7.11 **6.64 ***4.40 **3.22 ***15.2114.8946.32 **45.75 **2.51 ***1.97 **0.59810.73260.48710.538740.82 **40.93 ***
LSD5%2.591.0832.7715.153.721.543.493.784.193.460.160.140.450.220.190.311.071.000.090.130.23080.46370.16750.61901.600.52
LSD1%3.841.6048.5022.415.502.285.175.606.215.130.240.200.670.330.280.461.591.480.140.200.34150.68630.24790.91612.370.77
LSD0.1%5.942.4974.9634.648.503.527.998.669.597.920.380.311.040.510.440.722.462.280.220.310.52781.06060.38311.41583.671.19
SL 27SW 28HL 29P 30G 31Z 32WW65 33WW91 34DW65 35DW91 36RWC 37
WNWWNWWNWWNWWNWWNWWNWWNWWNWWNWWNW
NIL1−−7.057.022.972.8778.3779.6418.1817.7538.1336.9865.24 *63.250.13990.14080.07350.04560.03580.02970.06140.037173.2068.58
NIL++7.017.022.972.8878.7280.2417.4417.8236.8937.1664.0063.670.1567 ***0.15250.06850.05030.036930.02830.05900.042873.0767.80
LSD5%0.06010.03440.05020.05322.081.021.080.87541.581.351.191.030.00620.00840.01020.00940.00520.00360.01070.00882.304.11
LSD1%0.08900.05090.07430.07883.081.511.601.29562.341.991.761.520.00920.01240.01510.01390.00770.00530.01590.01313.406.08
LSD0.1%0.13760.07860.11490.12184.762.332.482.00223.633.082.732.350.01420.01920.02340.02160.01200.00830.02460.02035.269.40
1 Days to heading; 2 days to flowering; 3 days to maturity; 4 SPAD value at Z45; 5 SPAD value at Z65; 6 SPAD value at Z77; 7 SPAD value at Z83; 8 SPAD value at Z85; 9 NDVI at Z45; 10 NDVI at Z65; 11 NDVI at Z83; 12 area under SPAD value decline curve; 13 area under vegetation index curve; 14 flag leaf area; 15 fertile tiller number; 16 plant height up to the flag leaf collar; 17 plant height up to the base of the ear; 18 plant height up to the tip of the ear; 19 ear length; 20 grain yield; 21 spikelet number/main spike; 22 seed number/main spike; 23 seed weight/main spike; 24 basal sterile spikelet number/main spike; 25 apical sterile spikelet number/main spike; 26 thousand-grain weight; 27 seed length; 28 seed width; 29 test weight; 30 protein content; 31 glutenin content; 32 Zeleny index; 33 fresh biomass at Z65; 34 fresh biomass at Z91; 35 dry biomass at Z65; 36 dry biomass at Z91; 37 relative water content.
Table 4. Correlations between morphological and physiological parameters and yield components in the rain-fed treatment. *, **, *** significant at the p < 0.05, 0.01 and 0.001 probability levels, respectively.
Table 4. Correlations between morphological and physiological parameters and yield components in the rain-fed treatment. *, **, *** significant at the p < 0.05, 0.01 and 0.001 probability levels, respectively.
TraitsSimple Regression between the Traits (R2)
GY 1TGW 2SNM 3SWM 4SW 5
Grand mean for years and treatments
NDVI45 6
NDVI65 7
NDVI83 838.16 *** 18.29 *45.78 ***
SPAD45 9
SPAD65 1013.33 *16.49 * 27.78 ***25.26 **
SPAD77 1132.73 *** 38.41 ***38.98 ***
SPAD83 1244.38 *** 21.38 ***
SPAD85 1316.51 *
FLC 1415.27 *33.11 *** 16.00 *
BE 1533.35 ***
TE 1634.04 ***12.81 *
EL 17 37.42 *** 34.36 ***
FTN 1840.41 *** 48.87 ***22.79 ***
FLA 19 14.65 *
WW65 2017.55 *18.89 * 37.34 ***30.37 **
DW65 2112.74 *31.45 *** 39.07 ***15.08 *
1 Grain yield; 2 thousand-grain weight; 3 seed number/main spike; 4 seed weight/main spike; 5 seed width; 6 NDVI at Z45; 7 NDVI at Z65; 8 NDVI value at Z83; 9 SPAD at Z45; 10 SPAD at Z65; 11 SPAD value at Z77; 12 SPAD value at Z83; 13 SPAD value at Z85; 14 plant height up to the flag leaf collar; 15 plant height up to the base of the ear; 16 plant height up to the tip of the ear; 17 ear length; 18 fertile tiller number; 19 flag leaf area; 20 fresh biomass at Z65; 21 dry biomass at Z65.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bányai, J.; Maccaferri, M.; Láng, L.; Mayer, M.; Tóth, V.; Cséplő, M.; Pál, M.; Mészáros, K.; Vida, G. Abiotic Stress Response of Near-Isogenic Spring Durum Wheat Lines under Different Sowing Densities. Int. J. Mol. Sci. 2021, 22, 2053. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22042053

AMA Style

Bányai J, Maccaferri M, Láng L, Mayer M, Tóth V, Cséplő M, Pál M, Mészáros K, Vida G. Abiotic Stress Response of Near-Isogenic Spring Durum Wheat Lines under Different Sowing Densities. International Journal of Molecular Sciences. 2021; 22(4):2053. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22042053

Chicago/Turabian Style

Bányai, Judit, Marco Maccaferri, László Láng, Marianna Mayer, Viola Tóth, Mónika Cséplő, Magda Pál, Klára Mészáros, and Gyula Vida. 2021. "Abiotic Stress Response of Near-Isogenic Spring Durum Wheat Lines under Different Sowing Densities" International Journal of Molecular Sciences 22, no. 4: 2053. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22042053

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop