Next Article in Journal
Molecular T-Cell Repertoire Analysis as Source of Prognostic and Predictive Biomarkers for Checkpoint Blockade Immunotherapy
Next Article in Special Issue
Conventional and Molecular Techniques from Simple Breeding to Speed Breeding in Crop Plants: Recent Advances and Future Outlook
Previous Article in Journal
Nanogels of a Succinylated Glycol Chitosan-Succinyl Prednisolone Conjugate: Release Behavior, Gastrointestinal Distribution, and Systemic Absorption
Previous Article in Special Issue
OsINV3 and Its Homolog, OsINV2, Control Grain Size in Rice
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mapping of QTLs Associated with Yield and Yield Related Traits in Durum Wheat (Triticum durum Desf.) Under Irrigated and Drought Conditions

1
Nuclear Institute for Agriculture and Biology, Jhang Road, Faisalabad, 38000, Pakistan
2
Department of Plant and Environmental Sciences, University of Copenhagen, Højbakkegård Allé 13, DK-2630 Taastrup, Denmark
3
Leibniz Institute of Plant Genetics and Crop Plant Research, Corrensstr. 3, Seeland OT, 06466 Gatersleben, Germany
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(7), 2372; https://doi.org/10.3390/ijms21072372
Submission received: 18 February 2020 / Revised: 10 March 2020 / Accepted: 11 March 2020 / Published: 30 March 2020
(This article belongs to the Collection Genetics and Molecular Breeding in Plants)

Abstract

:
Global durum wheat consumption (Triticum durum Desf.) is ahead of its production. One reason for this is abiotic stress, e.g., drought. Breeding for resistance to drought is complicated by the lack of fast, reproducible screening techniques and the inability to routinely create defined and repeatable water stress conditions. Here, we report the first analysis of dissection of yield and yield-related traits in durum wheat in Pakistan, seeking to elucidate the genetic components of yield and agronomic traits. Analysis of several traits revealed a total of 221 (160 with logarithm of odds (LOD) > 2 ≤ 3 and 61 with LOD > 3) quantitative trait loci (QTLs) distributed on all fourteen durum wheat chromosomes, of which 109 (78 with LOD > 2 ≤ 3 and 31 with LOD > 3) were observed in 2016-17 (S1) and 112 (82 with LOD > 2 ≤ 3 and 30 with LOD > 3) were observed in 2017-18 (S2). Allelic profiles of yield QTLs on chromosome 2A and 7B indicate that allele A of Xgwm895 and allele B of Xbarc276 can enhance the Yd up to 6.16% in control and 5.27% under drought. Moreover, if combined, a yield gain of up to 11% would be possible.

Graphical Abstract

1. Introduction

Wheat has been the world’s most important cereal from 10,000 years ago since the dawn of agriculture [1] up to the present day. With a production of ~620 million tons annually worldwide, wheat provides about 1/5th of the calories consumed by humans [2]. The current annual wheat yield increase is below 1%, but it should be around 1.6% in order to meet the wheat demand (which is anticipated to rise by 60% with the world population, which will be 9 billion by 2050 [3].
Durum wheat (Triticum durum or Triticum turgidum subsp. durum), which is also referred to as pasta or macaroni wheat, is the hardest of all wheats. Its density, combined with its high protein content and gluten strength, make durum the wheat of choice for producing quality products including bread, couscous, frekeh, bulgur, and most importantly pasta [4]. A total of 9.3 and 10.5 million tons of pasta was produced in 2001 and 2003, respectively. By 2013, its production reached 13.5 million tons [5], providing an insight into the growing global demand for durum wheat.
Durum wheat global consumption is ahead of its production. T. durum accounts for around 6% of total wheat production (37.7 MT in 2013; International Grain Council, October 2014) [6], occupying approximately 20 million hectares worldwide. High yielding cultivars endowed with drought tolerance and disease resistance, in addition to their high commercial and technological value, are therefore highly desirable. Durum wheat is primarily grown under rain-fed conditions where the frequent drought combined with heat stress is the major aspect of grain yield reduction [7,8,9]. Drought is one of the most severe stresses and is the main cause of significant losses in growth and productivity of crop plants [10]. About 45% of the 120 million hectares of land allocated in developing countries to wheat production is prone to drought [11]. Crop yield is seriously affected by water deficit in many growing areas of the world. It has been assessed that as much as 50% of the wheat cultivation area is regularly affected by drought [12].
In response to water deficit, morphological and biochemical changes are witnessed in plants which follow the functional damage and loss of important plant parts as the shortfall increases [13]. Various tolerance mechanisms (pigment content and stability and high relative water content) have been ascribed to the physiological changes in plants to water scarcity [14]. There are various ways in which drought-tolerant wheat species can be differentiated. These include growth response in water-limited conditions, alteration in water relations in tissues subjected to water stress, amassing of various ions, and early flowering [15]. As our comprehension of the physiological processes by which plants respond to water scarcity increases, we will be able to devise screening methods to select drought-tolerant accessions for further usage in breeding systems [16].
Generally, all yield components are disturbed under severe drought stress, most notably including fertile ears per unit area (reduction ≤ 60%), seeds per spike (reduction ≤ 60%), harvest index; and nitrogen harvest index (reduction ≤ 24%), whereas lower grain weight is witnessed under mild stress [17]. Different cultivars differ in their responses to drought. Several factors (CO2 exchange rate, stomatal conductance, relative water contents and osmotic changes) regulate the fate of drought resistance in a genotype and natural variation in these factors can be used as selection criteria to breed drought-tolerant cultivars [18]. Since the 1970s, the attained genetic gains in durum yield have been because of well-adjusted progress in fertility through distribution of relatively higher assimilates to the ear and the growing tillers which increased total biomass, whereas harvest index remained the same [19], probably because the ratio of improvement of grain yield to biological yield remained the same. Water usage of plants and minimization of the deleterious effects under water-scarce conditions can be modified by altering plant height, exposed leaf area and leaf index [20]. Moreover, varieties enriched with increased leaf turgidity and relative water contents in water-deficient environments have better yields in drought [21,22].
Several factors are responsible for complications in breeding drought-tolerant cultivars, such as fast repeatable screening techniques and the incapacity to create drought conditions in a repeatable manner to screen large amounts of germplasm [23]. Breeders look for yield performance under drought stress. However, there are other factors that lead to poor yield in drought. To nullify this, drought indices have been devised based on plants performance under normal and drought conditions [24]. The drought resistance or susceptibility of a given genotype forms the basis of these indices [25]. Resistance to drought of a given genotype can be defined as the relative yield of that genotype in comparison to other genotypes that have been grown under the exact same stress conditions [26].
Different mapping populations have been made to investigate both biotic (disease resistance) and abiotic (drought, salinity) stresses [27,28], and many QTLs have been identified for yield and related traits in field trails. For example, two QTLs for test weight on chromosomes 7A and 6B (explaining 30% variation) and five QTLs of thousand kernel weight (Tkw) (explaining 32% of total, including 25% genetic, variation) have been revealed [17]. Other investigations in durum wheat have reported QTLs for yield and its components and physiological developmental [29,30]. Although a plethora of information regarding genetic maps, populations and molecular markers exists, the transformation of this information into the breeding of drought-tolerant cultivars cannot yet be compared. This anomaly can be attributed to three reasons. The first is that yield, being a quantitative trait, has vague genetic architecture. The second is that the yield-linked loci are subject to a strong environmental influence. The third is the poor phenotyping of the concerned traits. Established QTLs in one environment often disappear in a different environment, and thus QTLs have been well-defined as constitutive versus adaptive. Nevertheless, robust QTLs of yield traits across 16 different environments on chromosomes 2B and 3B have been reported [29].
Here, we report the analysis of a bi-parental mapping population seeking to elucidate the genetic components of yield and agronomic traits under irrigated and drought conditions. Moreover, certain physiological traits were analyzed to check their response towards yield in both conditions. To add to this, various stress indices were calculated and their QTLs were sought. Thus, the current study was undertaken to achieve the following objectives: (1) to screen a panel of 114 recombinant inbred lines (RILs) of durum wheat for yield and yield-related traits under drought and irrigated conditions at the Nuclear Institute for Agriculture and Biology (NIAB), Pakistan; (2) to measure physiological traits and various stress indices of the population under irrigated and drought conditions; and (3) to uncover the potential QTLs of measured traits and their comparison.

2. Results

2.1. Phenotypic Variation

A three-way ANOVA revealed that treatment, i.e., drought, significantly reduced all the traits (Table 1). Moreover, germination (Gr), transpiration rate (Tr), stomatal conductance (Sc), photosynthetic rate (Phr), yield (Yd), harvest index (HI), and thousand kernel weight (Tkw) differed significantly at the genotype level. Impact of season was significant for Gr, Sc, Phr and Sl, where an increase for Gr was observed in S2, while Sc Phr and Sl were significantly higher in S1. Furthermore, genotype x treatment interaction was significant for Tr, Sc, Ph and HI. Season x treatment had a significant effect on Gr, Sc, Phr, Wue, (heading time) Hd and Sl. Finally, genotype x season x treatment were significant for Sc and Yd.

2.2. Agronomic Traits

Mean ± SD along with range of the traits are provided in Table 1, whereas box plots are given in Figure S1. Agronomic traits included Gr, Hd and Ph. S1 (2016-17) and S2 (2017-18) indicate different growing seasons. Mean germination under control in S1 reduced from 81.86 ± 9.63 (Gr_S1) to 72.55 ± 9.96 (GrD_S1), with relative germination being 88.2 4 ± 11.78 (GrD_S1). Similarly, mean germination under control in S2 decreased from 96.00 ± 3.72 (Gr_S2) to 91.24 ± 5.56 (GrD_S2) in drought, with relative germination being 95.11 ± 5.30. Heading time in control in S1 decreased from 100.47 ± 11.79 (Hd_S1) to 93.94 ± 12.41 9 (HdD_S1) in drought. The same decreased from 100.47 ± 11.79 (Hd_S1) to 92.32 ± 6.79 (HdD_S2) under drought in S2. The relative heading times in S1 and S2 were 94.30±7.38 and 94.18 ± 2.70, respectively. Reduction in plant height was also observed in drought as compared to control, as it was decreased from 114.78 ± 23.82 (Ph_S1) in control to 80.24 ± 12.45 (PhD_S1) in drought in S1. Likewise, it was decreased from 111.36 ± 26.31 (Ph_S2) to 83.65 ± 14.14 (PhD_S2) in drought in S2. (not mark in total)

2.3. Physiological Traits

These included Tr, Phr, Sc and Wue, which were higher in control than their corresponding values in drought in both seasons. In S1, Tr dropped form 37.38 ± 6.86 (Tr_S1) to 15.56 ± 4.65 (TrD_S1) under drought. On the same ground, in S2 it dropped from 34.17 ± 6.51 (Tr_S2) to 15.06 ± 1.90 (TrD_S2). However, RTr_S1 (42.98 ± 15.98) was lower than RTr_S2 (45.43 ± 9.78). Phr reduced from 3.97 ± 1.24 (Phr_S1) to 2.22 ± 0.44 (PhrD_S1) in drought. Likewise, a reduction from 3.49 ± 1.25 (Phr_S2) to 1.72 ± 0.67 (PhD_S2) in drought was observed in S2. RPhr_S1 (58.96 ± 15.51) was higher than RPhr_S2 (50.94 ± 15.76). In control, Wue in control treatment in S1 was 11.89 ± 4.03 (Wue_S1), which increased to 16.04 ± 7.93 in drought, whereas it was 11.14 ± 6.30 (Wue_S2) in control treatment in S2, which increased to 11.43 ± 5.14 (WueD_S2) in drought. This also raised the mean RWue_S1 (144.41 ± 87.62) and RWue_S2 (113.98 ± 41.44).

2.4. Spike-Related Traits

These included Sl, seeds per spike (Sps), spike weight (Sw) and Tkw. As per spike-related traits, the values of Sl, Sw and Tkw were higher in S1, whereas Sps was higher in S2. In drought, spike length was dropped to 7.13 ± 1.11 (SlD_S1) from 8.25 ± 0.96 (Sl_S1) in Sl, where RSl_S1 was 86.71 ± 11.40. On the other hand, in drought, Sl dropped to 6.86 ± 0.66 (SlD_S2) from 7.60 ± 0.71 (Sl_S2) with a relative value of 90.64 ± 6.81 in S2. Sps also witnessed a reduction in drought treatment which dropped from 43.35 ± 9.60 (Sps_S1) in S1 to 33.44 ± 7.75 (SpsD_S1) in drought. Sps_S2 was 47.15±8.18 in control, which decreased to 35.52 ± 5.73 (SpsD_S2) in drought. RSps_S1 was (78.75 ± 14.88), but was, however, higher than RSps_S2 (77.17 ± 13.78). Mean Sw_S1 was 3.05±0.65, whereas mean Sw_S2 was 2.92±0.65. Both of them decreased to 2.30 ± 0.57 (SwD_S1) and 2.25±0.41 (SwD_S2) in drought. The corresponding relative values of Sw were 77.17 ± 15.18 (RSw_S1) and 78.98 ± 14.60 (RSw_S2), respectively. Mean Tkw_S1 (48.88 ± 5.30) and mean Tkw_S2 (42.08 ± 6.26) were reduced to 45.07 ± 4.95 and 38.99 ± 5.57, respectively, in drought, whereas relative TKW values were comparable.

2.5. Yield and Related Traits

Yd, biomass (Bm) and HI were included in this section. The highest decrease among all traits from control to drought was observed for Bm, which decreased from 3535.80 ± 1,086.20 (Bm_S1) in S1 and 3551.20 ± 1,172.60 (Bm_S2) in S2 to 1070.50 ± 349.64 (BmD_S1) and 1191.40 ± 373.78 (BmD_S2), respectively. In line with the same trend, Yd in S1 dropped from 1,033.10 ± 361.77 (Yd_S1) in control treatment to 350.92 ± 109.75 (YdD_S1) in drought. Likewise, Yd in S2 dropped from 1,010.10 ± 335.93 (Yd_S2) in the control treatment to 332.46 ± 116.73 in the drought treatment in S2. An increasing trend in HI was observed from control to drought in S1. For example, HI in S1 in control was 30.50 ± 10.09, which increased to 33.59 ± 7.96 (HID_S1). In S2, HI in control and drought remained almost the same. Furthermore, there was very slight increase in relative HI during both seasons (RHI_S1 = 117.66 ± 36.37, RHI_S2 = 103.30 ± 30.84).

2.6. Stress Indices

Altogether, the five stress indices, that is, stress tolerance index (STI), mean productivity (MP), stress tolerance (ST), stress susceptibility index (SSI), and drought resistance index (DRI), were close to each other in both S1 and S2. However, there were differences in their ranges. For example, STI_S1 ranged between 0.06 and 1.04, whereas STI_S2 ranged between 0.06 and 1.41. Likewise, MP_S1 ranged between 266.51 and 1236.5, whereas MP_S2 ranged between 258 and 1306.5. Ranges of ST_S1 and ST_S2 were 158.19–1509 and 120–1401, respectively. Similarly, ranges of SSI_S1 and SSI_S2 were 0.44–1.32 and 0.49–1.31, respectively. Finally, DRI_S1 ranged between 49.91 and 1225.7 and DRI_S2 ranged between 49.97 and 1211.8.
There was a significant difference (α = 0.05) observed between S1 and S2 with respect to trait values in control for Gr, Tr, Sc, Phr, Sl, Sps and Tkw (data not shown). Moreover, there was also a significant difference (α = 0.05) between S1 and S2 with respect to trait values under drought of Gr, Sc, Phr, Wue, Ht, Sps, Tkw, Bm and HI. As far as relative values are concerned, significant differences (α = 0.05) were observed between S1 and S2 for Gr, Phr, Wue, Ht, Sl and HI. However, there was no significant difference observed between the S1 and S2 trait values of all five drought indices (ST, MP, STI, SSI and DRI).

2.7. Correlations

As before, the data were divided into five sections for correlation analysis (relative values were excluded). These included (1) common agronomic traits (Gr, Hd and Ph), (2) physiological traits (Tr, Sc, Phr and Wue), (3) spike-related traits (Sl, Sw, Sps, Tkw), (4) yield-related traits (Bm, Yd and HI) and (5) stress indices (STI, MP, ST, SSI and DRI). All sections were analyzed simultaneously for both S1 and S2. Data are shown in Tables S1–S5.

2.8. QTL Mapping

QTL mapping analysis of several traits revealed a total of 221 (160 with LOD > 2 ≤ 3 and 61 with LOD > 3) quantitative trait loci (QTLs) were distributed on all of the fourteen durum wheat chromosomes, of which 109 (78 with LOD > 2 ≤ 3 and 31 with LOD > 3) were observed in S1 and 112 (82 with LOD > 2 ≤ 3 and 30 with LOD > 3) were observed in S2. In S1, the donor of 54 QTLs was ‘Omrabi 5’ and the donor of other 55 QTLs was ‘Belikh 2‘. Likewise, in S2, the donor of 54 QTLs was ‘Omrabi 5’ and the donor of other 58 QTLs was ‘Belikh 2’. These QTLs explained 5%–33.73% phenotypic variation in S1 and 5%–29.9% in S2. All details are given in Table 2.
From a chromosome perspective, the highest number of QTLs was observed on chromosome 5B (31 QTLs), followed by chromosomes 7B and 5A with 29 and 27 QTLs each, and by chromosomes 3B and 4B with 21 and 20 QTLs, respectively. A total of 18 QTLs was detected on chromosome 3A, whereas 15 QTLs were detected on each of the chromosomes 1B and 7A, and 14 QTLs resided on chromosome 2A. There were 10 and nine QTLs discovered on chromosomes 6B and 2B, respectively, whereas 6 QTLs detected on chromosome 6A and the least numbers of QTLs were detected on chromosomes 1A and 4A (3 QTLs each).

2.8.1. Agronomic Traits

A total of 17 QTLs were detected for Gr on 9 of the 14 chromosomes (1B (2 QTLs), 3A, 3B, 4B, 5A (2 QTLs), 5B (6 QTLs), 6A, 6B (2 QTLs) and 7B) with five of them had an LOD > 3.0. There were two QTLs common, whereas the phenotypic variation explained ranged between 6% and 18%. As far as Hd is concerned, there were 8 QTLs discovered on chromosomes 2A, 3A (2 QTLs), 5A (2 QTLs), 6B and 7B (2 QTLs), with three of them being highly significant. In addition, there was no QTL observed for Hd_S2 and RHd_S2. Altogether, QTLs explained 9%–18% of the variation in Hd. For Ph, a total of 19 (4 highly significant) QTLs were revealed on chromosomes 2A, 2B (2 QTLs), 3A (2 QTLs), 4B (6 QTLs), 5A (4 QTLs), 6A, 7A and 7B (2 QTLs) responsible for 6%–22% variation.

2.8.2. Physiological Traits

There were 14 QTLs identified for Tr across both seasons on chromosomes 2A, 2B, 3A (2 QTLs), 3B, 5A, 5B (3 QTLs), 6A and 7B (4 QTLs) where two QTLs (QTr_S2.NIAB-5B and QRTr_S2.NIAB-5B) were co-located. These QTLs explained 6–14% variation. As regards to Sc, 9 QTLs were discovered on chromosomes 1B (2 QTLs), 2B, 3B, 5A, 5B, 6A (2 QTLs) and 7A responsible for 7–17% variation and one QTL (QSc_S1.NIAB-1B.2) was highly significant. For Phr, 17 QTLs were found on all the chromosomes except 1A and 4A. Among them, 7 QTLs were highly significant. The variation explained by these QTLs was 6–24%. There were 11 QTLs uncovered on chromosomes 2A, 2B, 4B (3 QTLs), 5A (2 QTLs), 5B, 7A and 7B (2 QTLs) for Wue from which three were highly significant, explaining 6%–22% of the phenotypic expression.

2.8.3. Spike Related Traits

Among spike related traits, there were 10 QTLs residing on chromosomes 2A, 2B, 3A, 4B (2 QTLs), 5B, 7A (2 QTLs) and 7B (2 QTLs) for Sl where 4 QTLs were highly significant. These QTLs explained 7-21% variation in Sl. Likewise, 18 QTLs were found for Sw that resided on chromosomes 1A (2 QTLs), 2A (2 QTLs), 3A (3 QTLs), 3B (2 QTLs), 5A (3 QTLs), 5B (3 QTLs), 7A (2 QTLs) and 7B. With respect to Sps, 19 QTLs were discovered, which explained 6%–11% of the variation. These QTLs were lying on chromosomes 1A, 1B (3 QTLs), 3A (2 QTLs), 3B (3 QTLs), 4B, 5A, 5B (7 QTLs) and 7B. The 16 QTLs linked to Tkw were distributed on 8 different chromosomes (1B (3 QTLs), 2A, 4B (3 QTLs), 5A (3 QTLs), 5B, 6B (2 QTLs), 7A (2 QTLs) and 7B, of which only one was highly significant. All together, they explained 5%–12% of the variation.

2.8.4. Yield and Related Traits

A total of 19, 14, 16 and 12 QTLs were discovered for the yield related traits, respectively. The 19 QTLs of Bm were present on chromosomes 3A, 3B (3 QTLs), 4B, 5A, 5B (3 QTLs), 7A (3 QTLs) and 7B (7 QTLs). The 14 Yd QTLs were distributed on chromosomes 1B, 2A (3 QTLs), 3A, 3B (2 QTLs), 5A (2 QTLs), 5A, 5B, 6B (2 QTLs) and 7B (2 QTLs) explaining 6–30% of phenotypic expression. There were three chromosomes, viz., 2A, 3B and 7B, where two Yd QTLs collocated. A total of seven chromosomes (1B, 2B, 3B (3 QTLs), 4B (2 QTLs), 5B (2 QTLs), 7A and 7B (2 QTLs) carried 12 HI QTLs which explained 7–11% phenotypic differences with only one highly significant QTL.

2.8.5. Stress Indices

A total of 18 QTLs were uncovered for various stress indices viz. STI (4 QTLs), MP (6 QTLs), ST (3 QTLs) and DRI (5 QTLs) on chromosomes 1B (2 QTLs), 2A, 2B, 3B (3 QTLs), 4A (3 QTLs), 5A (5 QTLs), 5B, 6B and 7A, which provided 8–26% of the variation among various stress indices.

3. Discussion

Climate change is causing more frequent and intense periods of drought as overall rainfall levels decline. Approximately 40% of the world’s land surface is characterized as dry land which inhabits 2.5 billion people, which makes up 33% of the global population [31]. Drought is currently one of the main constraints that prevent crop plants from expressing their full genetic potential. This makes the identification of drought-tolerant genotypes extremely important to secure productivity. In the present study, 114 RILs derived from a cross between Omrabi5 and Belikh2, where Omrabi5 was a drought tolerant but Bleikh2 a salt and high temperature tolerant cultivar were subjected to drought stress for two seasons to dissect drought tolerance in durum wheat.
The genetic basis of most of the important traits in cereals is complex [32]. The same was observed in this investigation where the evaluation of 94 parameters resulted in the discovery of 221 QTLs distributed all over the durum wheat genome (Table 3). There were thirteen traits (Hd_S2, RHd_S2, TrD_S2, RSc_S1, PhrD_S1, RPhr_S1, RWue_S1, RWue_S2, RSl_S2, Sw_S2, RHI_S1, SSI_S1 and SSI_S2) for which no QTL was discovered, possibly because of the limitation of bi-parental population which is derived from two parents following fewer recombination events.

3.1. Phenotypic Differences between Treatments and Seasons

The different phenotypic responses to drought treatment and seasons within our RILs empowered us to genetically dissect their performance through QTL analysis. From yield perspective, yield was considerably reduced in drought in both S1 and S2. However, reduction was slightly more pronounced in S2 (RYd_S2 = 34.60 ± 12.27) as compared to S1 (RYd_S1 = 36.39 ± 12.21). This superior yield in S1 could be attributed to Sl, Sw and Sps, as these traits were higher in drought treatments in S1 than in the same treatment in S2. Seasonal difference might also be due to change in mean minimum and maximum temperatures during the two seasons as the mean minimum temperature in S1 (16.81 °C) was 1.24 °C less than mean minimum temperature in S2 (17.57 °C) although there was almost no difference on average temperatures during S1 and S2 (Table S6). Besides, the relatively improved performance in drought in S1 of RILs could also be ascribed to 11 mm rain received in the month of Jan (when plants are heading which is sensitive to abiotic stress [33]) in S1. Superior performance in control over treatment in both S1 and S2 could also be due to prolonged growing season (~ 1 week) beside other advantages as treatment caused early flowering that reduced the assimilates time to move to plant development and growth [34].

3.2. Agronomic Traits

Germination plays an important role in the plant stand, which guarantees a successful harvest. In this study, the markers linked to QGr_S1.NIAB-5A.1 (Xbarc342) and QGr_S1.NIAB-6B.1 (Xbarc178) have been reported to be linked with initial and relative germinations, respectively [35,36], has reported one Hd QTL on the same chromosome at 58.2 ± 38.6 cM. Therefore, it can be concluded that the same QTL may be present in our study influencing Gr. Another Gr QTL (QGr_S1.NIAB.5A.2) is located with several other QTLs linked with Hd, Yd, ST and DRI. Six QTLs of Gr on chromosome 5B indicating that this is an important site for genes linked to plant survival.
As regards the heading date, [29] discovered an Hd QTL at 40 cM on chromosome 2A in a durum wheat population whereas the QTL QHdD_S1.NIAB.2A in our present study was at 5cM on chromosome 2A. The QTL QHd_S1.NIAB-3A is highly comparable with the one detected by [35] for relative germination. Moreover, it is also the site of many other QTLs linked to spike related traits in this study indicating that plant development influences spike related traits and in turn final yield. It should be mentioned that a flowering time locus has been reported by [37] on chromosome 3AL. Another QTL QHd_S1.NIAB-5A is linked to Yd, MP and ST QTLs in this study. In addition, this QTL and QHdD_S2.NIAB-5A are in close vicinity to each other. Thus, the Hd QTL on chromosome 5A can be a constitutive QTL whose selection may help in yield increase as QHdD_S2.NIAB-5A is also linked to Yd, ST and DRI QTL in this investigation. Two Hd QTLs on chromosome 5A have been reported by [32]. Moreover, [37] discovered an Hd QTL on chromosome 5A, affirming the finding of this study. On chromosome 7B, [29] have located an Hd QTL between 0 and 44.3 cM whereas our QTL (Q.RHd_S1.NIAB-7B) is located at 42.5 cM confirming the findings of [29].
Most of the Ph QTLs were concentrated on chromosomes 4B and 5A. Many Ph QTLs have been reported to be located on chromosomes 1A, 1B, 1D, 2D, 3A, 4A, 4B, 6A, 7A and 7B [32,38,39]. The marker linked to QPhD_S2.NIAB-4B and QRPh_S2.NIAB-4B (Xbarc193) has been reported to be linked with coleoptile length by [40]. The same area has also been reported to carry Ph QTLs in durum wheat [39]. Likewise QPh_S1.NIAB-6A is highly comparable with the QTL discovered by [32] on chromosome 6A. Ref. [38] reported a QTL for Ph under drought conditions between 38 and 50 cM and another QTL under well-watered conditions at the distal end of chromosome 7AL. In addition, another QTL for Ph under drought stress has been reported to be on chromosome 7B by [38], which is linked to marker Xbarc1073. Our QTL QPhD_S2.NIAB-7A linked to marker Xwmc517 can be highly comparable to the one reported by [38] because the marker Xwmc517 is closely linked to Xbarc1073 in [38]. Several Ph QTLs were co-located with QTLs of other traits. For example, QPhD_S2.NIAB-5A was co-located with Sc, Tr, Wue, Sw and Tkw loci. To add to this, QRPh_S1.NIAB-7B was co-located with a Tr QTL.

3.3. Physiological Traits

In our study, the QTL QTrD_S1.NIAB-2B for Tr may be the same as Qgs-2B.1 (for Tr) by [38], as both lie on the proximal ends of chromosome 2BS. Another marker (Xwmc50) of QRTr_S1.NIAB-3A is reported to be linked with relative germination after experimental ageing [35] indicating the presence of stress resistance genes in this area. Ref. [38] located a QTL linked to Tr under irrigated conditions at the distal end of chromosome 5AS, whereas the QTL QTr_S2.NIAB-5A in our investigation was at the far end of chromosome 5AL, indicating that QTr_S2.NIAB-5A is a novel QTL. Similarly, the two QTLs (QTr_S2.NIAB.5B and QRTr_S2.NIAB.5B) at 158 cM linked to marker Xgwm408 are highly comparable to one identified by [37] where the authors located a QTL linked to transpiration rate at early grain filling stage under irrigated conditions. Ref. [35] also located a Tr QTL in the distal region of chromosome 5AL. The QTL QTrD_S1.NIAB.7B.2 is connected to marker Xbarc276, which is the same marker reported by [37] to be linked with Tr under irrigated conditions on chromosome 7B.
In case of Sc, the QTL QScD_S2.NIAB.2B at 107.81 cM on 2B may be the same as observed by [38] at ~ 95 cM under drought stress. Ref. [37] detected several QTLs linked to leaf porosity (LP) on the distal end of chromosome 3BS. Our QTL QScD_S2.NIAB-3B at 13.91 cM on 3B might correspond to any of the LP QTLs of [37].
The 17 QTLs linked to Phr are distributed on all chromosomes except 1A, 4A and 5A, and most of them were discovered in S2 only. Nevertheless, QPhrD_S2.NIAB-3B.1 and QPhrD_S2.NIAB-3B.2 mirrored to a photosynthesis QTL identified by [37] under drought stress. The linked marker in both studies was Xgwm299. [35] located one Phr QTL under drought in the distal region of chromosome 3BS. In addition, [37] also located another Phr QTL under drought stress in the centromere region of chromosome 2B whereas in our study QRPhr_S2.NIAB.2B was located at the distal end of chromosome 2B. One drought and one well-watered Phr QTL by [35] and one LP QTL by [37] was located on chromosome 4B. Whether QPhr_S1.NIAB-4B matches any of them can only be speculated at this stage. To add to it, like in this study, [35] located Phr QTL on chromosomes 6B and 7A. However, they were at different positions. The QTL QPhrD_S2.NIAB-7B co-locates with normal seedling and total germination QTL of [40]. Interestingly, another QTL on 7B (QRPhr_S2.NIAB-7B) co-locates with Tkw, seed area and seed length by [40]. Furthermore, this QTL can also be mirrored with Phr QTLs by [41].
In the case of Wue, [38] located Wue QTLs at 94 cM on 2B and at ~ 238 cM on 7A, whereas QWueD_S1.NIAB-2B and QWue_S1.NIAB-7A in our study were located at the very distal ends of 2BS and 7AS, respectively, and hence are not comparable to the ones detected by [38]. However, the QTL QWueD_S2.NIAB-7B co-located with QRPhr_S2.NIAB-7B mirrored with Phr QTLs by [41].

3.4. Spike-Related Traits

With respect to Sl, QSl_S1.NIAB-2A mirrors with a grain protein content (GPC) QTL in [42], whereas QSlD_S1.NIAB-2B mirrors with Sps QTL in [43]. This QTL, however, is not the same as that described by [44], who detected consistent large effect QTLs for Sl in six environments. To add to this, Xwmc50 and Xgwm940a of QTLs QRSl_S1.NIAB-3A and QSlD_S2.NIAB-4B, respectively, showed connections with germination under normal and stress conditions [35]. The QTL QSl_S2.NIAB-5B is also at the location of numerous QTLs linked with yield related traits in [45]. Similarly, a QTL linked to Sl has been reported by [45] on the site of QSlD_S2.NIAB-7A. Ref. [38] also reported a QTL linked to Sl under drought stress at 160 cM on chromosome 7A. Furthermore, two QTLs at ~ 73 cM and 148 cM have been reported to be linked with Sl in irrigated and drought conditions, respectively. We speculate the 7B QTLs of Sl (QRSl_S1.NIAB-7B and QRSl_S2.NIAB-7B) might correspond to any of those mentioned in [38]. The site of QTLs QSwD_S1.NIAB-1A and QSw_S1.NIAB-1A is reported to carry Sps QTLs [38]. Likewise, the QTLs QSwD_S1.NIAB-2A and QRSw_S1.NIAB-2A are in the region of GPC QTLs [42]. The QTLs on 3A (QSwD_S1.NIAB-3A and QRSw_S1.NIAB-3A) can be compared with relative germination QTLs [35] and Sw QTLs in [45]. The site of QTL QSwD_S1.NIAB-3B has been reported for Sl QTL in different environments [43]. Likewise, other QTLs have been known to be in the region of QTLs linked to yield and related traits [38].
Many QTLs linked to yield and spike traits have been mapped to chromosomes where Sps QTLs of this study are mapped [32,37,43]. In particular, the marker Xwmc218 of QSpsD_S2.NIAB-7B has been reported to be in connection with seed area and length in [40], providing a clear connection towards yield. Moreover, QRSps_S1.NIAB in this study is located at 104 cM mirrored with grains per spike QTL of [44] who reported the Sps QTL at 102.5 cM.
Four QTLs among 16 of Tkw were co-located with other Tkw QTLs. The 1B loci (QRTkw_S.NIAB-1B, QTkwD_S1.NIAB-1B.1 and QTkw_S1.NIAB.2) can be compared with the Yd QTLs of [38] and meta-QTL of Yd by [43]. Furthermore, QRTkw_S.NIAB-1B, QTkw_S2.NIAB-4B, QRTkw_S2.NIAB-4B, QTkDw_S1.NIAB-5A, QRTkw_S1.NIAB.5A and QRTkw_S1.NIAB.7B mirror relative germination QTLs by [35].

3.5. Yield and Related Traits

Chromosome 1B is an important site of yield QTLs as several yield QTLs on this chromosome in hexaploid [29] or durum wheat [38,46]. Thus, the QTL QRYd_S2.NIAB-1B could be one reported previously. Interestingly, this QTL also co-locates with coleoptile length under drought [40] in the same population, pointing towards stress tolerance genes in this area. Similarly, there were three QTLs for yield (QYdD_S1.NIAB-2A, QYdD_S2.NIAB-2A and QRYd_S1.NIAB-2A) on chromosome 2A where QYdD_S2.NIAB-2A and QRYd_S1.NIAB-2A were linked to the same marker (Xgwm895). Consistent QTLs for Yd have been detected on chromosome 2A [43,45]. The QTL QYdD_S1.NIAB-2A can be compared with 2AS Yd QTLs of [42], whereas the other two can be compared with the QTLs on 2AL. Ref. [29] also detected a Yd QTL on 2A linked to marker Xwmc177, which is closely linked to Xgwm726 (related to QYdD_S1.NIAB-2A). The 3A QTL (QYd_S2.NIAB-3A) seems not to be comparable to those observed by [29] in durum wheat and by [36] in hexaploid wheat, because their QTLs were either in the centromere region or at the distal end of 3AL whereas our QTL was at 52.11 cM. However, this QTL could be compared Yd QTL described by [34], who described a cluster of QTLs including Yd QTL between 41.3 and 51.4 cM in durum wheat. Moreover, [46] have also located a Yd QTL under drought stress at 62.9 cM. To add to it, QYd_S2.NIAB-3A does co-locate with a grain number QTL of [32]. The QTLs viz. QYd_S2.NIAB-3B and QRYd.NIAB-3B linked to marker Xgwm547 are at the tip of long arm of 3B. Consistent QTLs for Yd under heat, drought and irrigated conditions have been described by [34] at 212 cM on chromosome 3BL indicating that our QTL could be the same as that of [34]. Likewise, both [29] and [45] have located QTLs for Yd at the same location, whereas [36] located a Tkw QTL in hexaploid wheat in this area indicating the connection of Yd with Tkw. One of the two 5A QTLs (QYd_S1.NIAB-5A and QRYd_S2.NIAB-5A) co-locates with a germination QTL in [35], whereas [32] located several Sl, Sw and Tkw QTL on the entire length of chromosome 5A, suggesting some connection between yield and these components. In addition, [32] have reported numerous QTLs linked to yield and related traits on chromosome 5A [47]. Furthermore, both [32] and [47] have reported Yd QTLs on the proximal end of chromosome 5BS, whereas QYdD_S1.NIAB.5B lies at 143.21 cM, and so cannot be compared with either of the above-mentioned investigations. Ref. [29] located a yield locus exactly in the same area as that of QRYd_S1.NIAB.6B. Nevertheless, [32] mapped numerous Yd related QTL on chromosome 6B. Chromosome 7B has been reported to be an important site for Yd and Yd related traits by [32], and [29] located a stable QTL for Yd between 19 and 33 cM and one QTL restricted to one environment between 106 and 122 cM. This study linked QYd_S1.NIAB.7B and QYdD_S2.NIAB.7B (216-226 cM) with Xbarc276, which is in between Xgwm333 and Xbarc176. This interval is linked with a Yd QTL in [29].
There are no investigations exclusively reporting Bm QTLs except [38] and [48]. No QTLs of Bm in the present study match with those reported by [38] and [48]. However, the QTL QBmD_S1.NIAB-7A of 7A is very near to GPC and Sl QTL as reported by [32] in hexaploid wheat. Likewise, the QRBm_S2.NIAB-7B of 7B is linked with QTLs of seed area and seed length [41]. Among other traits, the Bm QTLs co-locate with Gr, Ph, HI, Sw, Sps, Yd and HI QTLs in this population indicating that these are the sites of important genes.
The importance of HI in relation to drought has been indicated recently [36,44]. None of the 12 HI QTLs in this study matched with those of [38], who located HI QTLs in hexaploid wheat on chromosomes 1A, 2A, 2B, 3A, 4A, 4D and 6A. Numerous QTLs for HI have been identified on chromosomes 1B, 5B, 5D 6A, 7A, 7B and 7D using an association mapping approach and a different marker system [49]. The 1B QTL of HI (QRHI_S1.NIAB-1B) in this study was located at 93.11 cM whereas the QTL in [49] was located at 158 cM; hence, they seem not to be comparable with each other. QRHI_S1.NIAB-1B also co-locates with coleoptile length under drought stress in [41]. The QTL QHID_S.NIAB-5B on 5B can be the same as that of [49] who discovered an HI QTL under non stress conditions. The chromosome 7B carried QTLs for HI in non-stress, drought, heat and combined heat and drought stress most notable at 163 cM in [49], whereas our QTLs QRHI_S2.NIAB-7B and QHI_S2.NIAB-7B were at a distance of 184.41 and 336.01 cM, respectively.

3.6. Stress Indices

Several tolerance indices have been devised and used by breeders to identify the best drought-tolerant lines. It has been suggested that each index may be a potential indicator of different biological responses to drought [31]. STI, MP and geometric mean productivity have been concluded to be effective indices in discriminating drought tolerance genotypes in durum wheat [50]. Chromosome 2B has been an important site harboring STI QTLs under drought, heat and combined drought and heat stress [49]. However, our QTL QDRI_S2.NIAB-2B does not correspond to any of those QTLs because of a different position. Chromosome 4A also carried STI QTLs under combined stress at 67.71 and 71.47 cM in [49] indicating some similarity with our 4A QTLs which were at 44.51 (QSTI_S2.NIAB-4A.1), 64.31 (QDRI_S1.NIAB-4A) and 100.81 (QSTI_S2.NIAB-4A.2) cM. To add to it, [50] also located various QTLs linked with yield and yield related stress indices on chromosome 4A. The marker Xgwm299 on 3B in association with QMP_S2.NIAB-3B.1 and QST_S2.NIAB-3B co-locates with Phr, Sw and Sps QTLs. Furthermore, the same marker is linked with Tkw [51] and Phr at late growing stages [40] in other studies. Furthermore, the interval of QDRI_S1.NIAB-5B (Xwmc28-Xgwm1072) is linked with HI, spike HI and spike number QTLs in durum wheat [52].

3.7. Co-Location QTLs and MAS for Durum Wheat Breeding

From the total of 221 QTLs, 145 could be clustered (2 or more QTLs linked to the same marker) into a total of 53 clusters (Figure 1), covering all the durum wheat chromosomes except 4A. The number of QTLs clustered together ranged between two and six. Among them, there were at least 20 markers that were associated with any trait in both seasons, whereas 10 markers were associated with the same trait under various conditions in a single season. Moreover, there were four markers that were associated with the same trait under different conditions in both seasons. It is also evident that certain Yd QTLs match exactly with Gr (cluster 5A.5), Tr (cluster 7B.5), Hd (cluster 5A.4, 5A.5), Tkw (clusters 1B.1, 2A.1), HI (cluster 3B.5), Sps (cluster 5B.2) and stress indices (clusters 5A.4, 5A.5) QTLs. Clusters of QTLs have also been reported previously [32,34,38,45,47], indicating that multiple genes might be present at such sites which can either be constitutive or expressive. Plant breeders, however, tend to look for stable loci to improve yield. Although several QTLs can be targeted for yield enhancement, the most promising one could be the clusters 2A.4 and 7B.5, where markers Xgwm895 and Xbarc276 were associated with Yd in both S1 and S2. The allelic profile of yield QTLs in these clusters indicate that allele A of Xgwm895 and allele B of Xbarc276 can enhance the Yd up to 6.16% in control and 5.27% in drought (Table 4). Moreover, if combined, a yield gain of up to 11% is possible.

4. Materials and Methods

A set of 114 RILs was developed from a cross between a drought-tolerant cultivar Omrabi5 (P1) and a high temperature- and salinity-tolerant breeding line Belikh2 (P2) [17] at ICARDA, Syria, by repeated selfing of F1, using the single seed descent method [28]. Seeds from the 2015 season in Germany were multiplied in 2015-16 at NIAB, Pakistan, and were subsequently used for the investigation in 2016-17 and the seeds of the 2016-17 harvest were used in 2017-18.

4.1. Experimental Design and Traits Measured

A total of two experiments were conducted to investigate the population with respect to yield and other related traits. They were conducted in 2016-17 (S1) and 2017-18 (S2) seasons at NIAB, Pakistan, at exactly the same experimental site to avoid any soil interaction. Moreover, all the data were collected at almost the same dates in both experiments to minimize the effect of environment on phenotypic expression of the measured traits. The population was sown on 1st of November each season with one replicate in control and three under drought in both S1 and S2 cycles in the form of two rows per RIL. Seeds were sown manually with 5g seeds/row, in 2-meter rows using hand-driven drill and row to row distance was kept 20 cm. Control plots were irrigated four times (15th Nov., 1st Jan., 15th Feb. and 15th March), whereas drought plots were kept devoid of any irrigation. The drought and control experiments were in the same field, where the control experiment was separated from the drought by making a 2-ft-high ridge partition and there was a gap of at least 5 meters between drought and control treatments. Furthermore, there was also a gap of 5 meters between drought treatment and water channel to avoid any leach down effect that might have an effect on drought treatments. To add to it, we planted brassica in the leftover area so that the leached water from the water channel and control experiment is used up before it reached to the drought treatment. The total rainfall in S1 was 36.25 days, distributed over the period of 14 days (from Nov 2016-May 2017), where the maximum rain was received in April 2017 (16.58 mm). In S2, total rainfall was 42.48, distributed over a period of 30 days with maximum rainfall in April 2018 (Table S6). Furthermore, temperature data is provided in Table S7 and layout plan of the experiments is provided in Figure S1. A list of parameters recorded during the course of the two experiments is given in Table 4. Correspondingly, relative values were also calculated by dividing drought by the control times one hundred.

4.1.1. Agronomic Traits

Germination % (Gr/GrD) was recorded visually three weeks after sowing according to a scale ranging from 0%–100%, where 0% = no germination and 100% = complete germination. Heading time (Hd/HdD) was measured when 50% of spikes have left the flag leaf. Plant height (Ph/PhD) was recorded in cm from five random plants per replicate, with the help of meter rod, measured from the base of the plant to the spike tip excluding awns, before harvest.

4.1.2. Physiological Traits

Transpiration rate (Tr/TrD) was calculated using the handheld porometer (Model: LI-1600 steady state porometer) in µg cm−2s−1. Data were recorded separately on cloudless days for control (on 99th day of sowing) and drought (on 100th and 101st day of sowing) between 11:00 am and 2:00 pm. Each reading was a mean of three randomly selected flag leaves per replicate. Stomatal conductance (Sc/ScD) (mmol m−2s−1) was calculated using the formula: 1/DR × CF where DR is diffusible resistance calculated from porometer, CF = correction factor which is calculated using the formula: LT × constant where the value of constant at 25 °C for durum was 411.8, and LT was leaf temperature calculated from porometer. Photosynthetic rate (Phr/PhrD) was determined from Sc and Tr using the formula: Phr = Sc (mmol m−2s/Tr (µg cm−2s−1) × 10. Water use efficiency (Wue/WueD) was calculated by dividing the corresponding photosynthetic rate with transpiration rate times 100.

4.1.3. Spike-Related Traits

For spike-related traits, five spikes per row per replicate were collected randomly from the two experiments whereas for other traits, values from random plants were taken. Spike length (Sl/SlD) was calculated by taking mean length of five spikes/replicate using scale in centimeters. Spike weight (Sw/SwD) in grams was calculated by taking mean weight of five spikes using electronic weighing balance. The number of seed per spike (Sps/SpsD) was calculated by threshing and counting the seeds from five random spikes per replicate. Thousand kernel weight (Tkw/TkwD) in grams was recorded by weighing 1000 grains using electronic balance.

4.1.4. Yield and Related Traits

These included biomass (Bm), yield (Yd) and harvest index (HI). Biomass (Bm/BmD) (g) was measured by weighing the complete harvest of above ground part excluding roots per plot where harvesting was done manually to avoid any mixing of the straw. Yield (Yd/YdD) was recorded after manual threshing in grams per plot. Harvest index (HI/HID) was determined by dividing the corresponding Yd with respective Bm.

4.1.5. Stress Indices

Five stress indices were calculated using the yield data. These included stress tolerance index (STI), mean productivity (MP), stress tolerance (ST), stress susceptibility index (SSI), and drought resistance Index (DI), whose formulas are provided in Table 4.

4.2. Data Analysis

Statistica software was used for descriptive statistics of the measured traits and to construct box plots. ANOVA was performed using the RStudio-1.0.153.exe program.

4.3. Genotyping and Genetic Mapping

The genetic map of Omrabi5 × Belikh2 RIL population comprising 265 microsatellite loci, spanning 2864 cM (10.8 cM mean inter-marker separation) and covering all the 14 chromosomes of durum wheat, was available from [40]. The individual chromosomes ranged in genetic length from 37 cM (chromosome 1A) to 366 cM (chromosome 7B). The 265 microsatellites comprised 159 GWM, 62 BARC and 44 WMC markers. Further details about the map are available in [40]. QTLs were assigned using a composite interval mapping method (model 6 with forward stepwise regression) implemented in QTL Cartographer v2.5 [53]. To control the effect of genetic background, five markers, identified by forward regression, were used as a cofactor with a window size of 5cM. A LOD score of >2.0≤ 3 was applied to detect QTLs as significant and >3.0 as highly significant.

5. Conclusions

In the present study, we genetically dissected the performance of durum wheat RIL population for yield and yield-related traits and various morpho-physiological characteristics under irrigated and drought conditions through QTL mapping, which was mapped with the 265 assays. The experiments were conducted in two consecutive seasons at the exact same location to avoid variation due to soil properties and ground factor. A total of 221 QTLs were identified, where 145 QTLs were co-located with each other, and where 20 markers were associated with any trait in either S1 or S2. Furthermore, certain Yd QTLs match exactly with Gr, Tr, Hd, Tkw, HI, Sps and stress index QTLs, a finding consistent with previous reports. Several QTLs can be targeted for yield enhancement in durum wheat; the most promising ones from this study could be the clusters 2A.4 and 7B.5, where markers Xgwm895 and Xbarc276, if combined, might increase the yield up to 11%. In addition, an indirect improvement in yield is also possible by improving other traits (Tr, Tkw and HI), as we have found co-located QTLs for the said traits with Yd. Meanwhile, more studies are needed to validate our findings, our study still provides new insights into the genetic architecture of yield and related traits in durum wheat and the RILs with favorable loci could be used to improve yield in durum wheat.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/21/7/2372/s1.

Author Contributions

Conceptualization, M.A.R.A., S.S. and A.B.; Formal analysis, M.A.R.A.; Funding acquisition, M.A.R.A. and A.B.; Investigation, F.A., S.S., S.A. and M.Q.W.; Methodology, F.A., S.S., S.A., M.Q.W. and A.A.; Supervision, A.A.; Writing—original draft, M.A.R.A. and A.A.; Writing—review & editing, A.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by International Foundation for Science (IFS, Sweden), grant number C-5897-I.

Acknowledgments

The authors would also like to thank Sultan Mehmood Khan and Manzoor Ahmad for excellent technical assistance.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

QTLQuantitative trait loci
LODLog of odds
S1Season 1 (2016-17)
S2Season 2 (2017-18)
GrGermination
HdHeading time
PhPlant height
TrTranspiration rate
PhrPhotosynthetic rate
ScStomatal conductance
WueWater use efficiency
SlSpike length
SwSpike weight
SpsSeeds per spike
TkwThousand kernel weight
BmBiomass
YdYield
HIHarvest index
STIStress tolerance index
SSIStress susceptibility index
MPMean productivity
DRIDrought resistance index
STStress tolerance

References

  1. Lev-Yadun, S.; Gopher, A.; Abbo, S. The cradle of agriculture. Science 2000, 288, 1602–1603. [Google Scholar] [CrossRef] [PubMed]
  2. Tilman, D.; Balzer, C.; Hill, J.; Befort, B.L. Global food demand and the sustainable intensification of agriculture. Proc. Nat. Acad. Sci. USA 2011, 108, 20260–20264. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Breakout session P1.1 National Food Security–The Wheat Initiative–an International Research Initiative for Wheat Improvement. Available online: http://www.fao.org/docs/eims/upload/306175/Briefing%20Paper%20(3)-Wheat%20Initative%20-%20H%C3%A9l%C3%A8ne%20Lucas.pdf (accessed on 22 December 2019).
  4. Ames, N.P.; Clarke, J.M.; Marchylo, B.A.; Dexter, J.E.; Woods, S.M. Effect of environment and genotype on durum wheat gluten strength and pasta viscoelasticity. Cereal Chem. 1999, 76, 582–586. [Google Scholar] [CrossRef]
  5. World Pasta Market and Turkey. Available online: http://www.magazinebbm.com/english/world-pasta-market-and-turkey/.html (accessed on 22 December 2019).
  6. International Grains Council. Available online: http://www.igc.int/en/Default.aspx (accessed on 22 December 2019).
  7. Araus, J.L.; Bort, J.; Steduto, P.; Villegas, D.; Royo, C. Breeding cereals for Mediterranean conditions: ecophysiological clues for biotechnology application. Ann. Appl. Biol. 2003, 142, 129–141. [Google Scholar] [CrossRef]
  8. Araus, J.L.; Slafer, G.A.; Reynolds, M.P.; Royo, C. Plant breeding and drought in C3 cereals: what should we breed for? Ann. Bot. 2002, 89, 925–940. [Google Scholar] [CrossRef] [PubMed]
  9. Condon, A.G.; Richards, R.A.; Rebetzke, G.J.; Farquhar, G.D. Breeding for high water-use efficiency. J. Exp. Bot. 2004, 55, 2447–2460. [Google Scholar] [CrossRef] [Green Version]
  10. Ludlow, M.M.; Muchow, R.C. A critical evaluation of traits for improving crop yields in water-limited environments. Adv. Agron. 1990, 43, 107–153. [Google Scholar]
  11. Rajaram, S. Prospects and promise of wheat breeding in the 21st century. In Wheat in a Global Environment; Springer: Dordrecht, The Netherlands, 2001; pp. 37–52. [Google Scholar]
  12. Ashraf, M.; Harris, P. Abiotic Stresses: Plant Resistance through Breeding and Molecular Approaches; CRC Press: Boca Raton, FL, USA, 2005. [Google Scholar]
  13. Sangtarash, M.H. Responses of different wheat genotypes to drought stress applied at different growth stages. Pak. J. Biol. Sci. 2010, 13, 114–119. [Google Scholar] [CrossRef] [Green Version]
  14. Clarke, J.M.; Richards, R.A.; Condon, A.G. Effect of drought stress on residual transpiration and its relationship with water use of wheat. Can. J. Plant Sci. 1991, 71, 695–702. [Google Scholar] [CrossRef] [Green Version]
  15. Blum, A. Plant Breeding for Stress Environments; CRC Press: Boca Raton, FL, USA, 1998. [Google Scholar]
  16. Zaharieva, M.; Gaulin, E.; Havaux, M.; Acevedo, E.; Monneveux, P. Drought and heat responses in the wild wheat relative Aegilops geniculata Roth. Crop Sci. 2001, 41, 1321–1329. [Google Scholar] [CrossRef]
  17. Nachit, M.M.; Elouafi, I. Durum wheat adaptation in the Mediterranean dryland: Breeding, stress physiology, and molecular markers. In Challenges and Strategies of Dryland Agriculture, (Challengesandst); Crop Science Society of America and American Society of Agronomy: Madison, WI, USA, 2004; pp. 203–218, CSSA special Publication no. 32. [Google Scholar]
  18. Houshmand, S.; Arzani, A.; Mirmohammadi-Maibody, S.A.M. Effects of salinity and drought stress on grain quality of durum wheat. Commun. Soil Sci. Plant 2014, 45, 297–308. [Google Scholar] [CrossRef]
  19. Royo, C.; Di Fonzo, N.; Nachit, M.M.; Araus, J.L. Enhancing genetic grain yield potential and yield stability in durum wheat. In Proceedings of the Seminar Durum Wheat Improvement in the Mediterranean Region: New Challenges. Options Méditerranéennes, Zaragoza, Spain, 12–14 April 2000; Volume 40, pp. 88–93. [Google Scholar]
  20. Mitchell, J.; Siamhan, D.M.; Wamala, M.; Risimeri, J.E.; Chinyamakobvu, E.; Henderson, S.; Fukai, S. The use of seedling leaf death score for evaluation of drought resistance of rice. Field Crop Res. 1998, 55, 129–139. [Google Scholar] [CrossRef]
  21. Akram, M. Growth and yield components of wheat under water stress of different growth stages. Ban. J. Agric. Res. 2011, 36, 455–468. [Google Scholar] [CrossRef] [Green Version]
  22. Khakwani, A.A.; Dennett, M.; Munir, M. Drought tolerance screening of wheat varieties by inducing water stress conditions. Songklanakarin J. Sci. Technol. 2011, 33, 135–142. [Google Scholar]
  23. Ramirez-Vallejo, P.; Kelly, J.D. Traits related to drought resistance in common bean. Euphytica 1998, 99, 127–136. [Google Scholar] [CrossRef]
  24. Mitra, J. Genetics and genetic improvement of drought resistance in crop plants. Curr. Sci. 2001, 80, 758–763. [Google Scholar]
  25. Fernandez, G.C. Effective selection criteria for assessing plant stress tolerance. In Proceedings of the International Symposium on Adaptation of Vegetables and other Food Crops in Temperature and Water Stress, Shanhua, Taiwan, 13–16 August 1992; pp. 257–270. [Google Scholar]
  26. Hall, A. Is dehydration tolerance relevant to genotypic difference in leaf senescence and crop adaption to dry environments. In Plant Responses to Cellular Dehydration during Environmental Stress; Close, T.J., Bray, E.A., Eds.; The American Soc. Plant Pathologists: Rockville, MD, USA, 1993; pp. 1–10. [Google Scholar]
  27. Blanco, A.; Bellomo, M.; Cenci, A.; De Giovanni, C.; D’ovidio, R.; Iacono, E.; Laddomada, B.; Pagnotta, M.; Porceddu, E.; Sciancalepore, A. A genetic linkage map of durum wheat. Theor. Appl. Genet. 1998, 97, 721–728. [Google Scholar] [CrossRef]
  28. Nachit, M.M.; Elouafi, I.; Pagnotta, M.A.; El-Saleh, A.; Iacono, E.; Labhilili, M.; Asbati, A.; Azrak, M.; Hazzam, H.; Benscher, D.; et al. Molecular linkage map for an intraspecific recombinant inbred population of durum wheat (Triticum turgidum L. var. durum). Theor. Appl. Genet. 2001, 102, 177–186. [Google Scholar] [CrossRef]
  29. Maccaferri, M.; Sanguineti, M.C.; Corneti, S.; Ortega, J.L.A.; Salem, M.B.; Bort, J.; DeAmbrogio, E.; del Moral, L.F.G.; 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] [Green Version]
  30. Sourdille, P.; Perretant, M.; Charmet, G.; Leroy, P.; Gautier, M.; Joudrier, P.; Nelson, J.; Sorrells, M.E.; Bernard, M. Linkage between RFLP markers and genes affecting kernel hardness in wheat. Theor. Appl. Genet. 1996, 93, 580–586. [Google Scholar] [CrossRef]
  31. Nauraein, M.; Mohammadi, S.A.; Ahrizad, S.; Moghaddam, M.; Sadeghzadeh, B. Evaluation of drought tolerance indices in wheat recombinant inbred population. Ann. Biol. Res. 2013, 4, 113–122. [Google Scholar]
  32. Börner, A.; Schumann, E.; Fürste, A.; Cöster, H.; Leithold, B.; Röder, M.S.; Weber, W.E. Mapping of quantitative trait loci determining agronomic important characters in hexaploid wheat (Triticum aestivum L.). Theor. Appl. Genet. 2002, 105, 921–936. [Google Scholar] [CrossRef] [PubMed]
  33. Dolferus, R.; Ji, X.M.; Richards, R.A. Abiotic stress control of grain number in cereals. Plant Sci. 2011, 181, 331–347. [Google Scholar] [CrossRef] [PubMed]
  34. 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, 173, 1473–1485. [Google Scholar] [CrossRef] [PubMed]
  35. Rehman-Arif, M.A.; Börner, A. Mapping of QTL associated with seed longevity in durum wheat (Triticum durum Desf.). J. Appl. Genet. 2019, 60, 33–36. [Google Scholar] [CrossRef] [PubMed]
  36. Peleg, Z.V.I.; Fahima, T.; Krugman, T.; Abbo, S.; Yakir, D.A.N.; Korol, A.B.; Saranga, Y. Genomic dissection of drought resistance in durum wheat× wild emmer wheat recombinant inbreed line population. Plant Cell Environ. 2009, 32, 758–779. [Google Scholar] [CrossRef] [PubMed]
  37. Panio, G.; Motzo, R.; Mastrangelo, A.M.; Marone, D.; Cattivelli, L.; Giunta, F.; De Vita, P. Molecular mapping of stomatal-conductance-related traits in durum wheat (Triticum turgidum ssp. durum). Ann. Appl. Biol. 2013, 162, 258–270. [Google Scholar] [CrossRef]
  38. Xu, Y.F.; Li, S.S.; Li, L.H.; Ma, F.F.; Fu, X.Y.; Shi, Z.L.; Xu, H.X.; Ma, P.T.; An, D.G. QTL mapping for yield and photosynthetic related traits under different water regimes in wheat. Mol. Breed. 2017, 37, 34. [Google Scholar] [CrossRef]
  39. Hu, J.; Wang, X.; Zhang, G.; Jiang, P.; Chen, W.; Hao, Y.; Ma, X.; Xu, S.; Jia, J.; Kong, L.; et al. QTL mapping for yield-related traits in wheat based on four RIL populations. Theor. Appl. Genet. 2020, 133, 917–933. [Google Scholar] [CrossRef]
  40. Nagel, M.; Navakode, S.; Scheibal, V.; Baum, M.; Nachit, M.; Röder, M.S.; Börner, A. The genetic basis of durum wheat germination and seedling growth under osmotic stress. Biol. Plant. 2014, 58, 681–688. [Google Scholar] [CrossRef]
  41. Wang, S.G.; Jia, S.S.; Sun, D.Z.; Wang, H.Y.; Dong, F.F.; Ma, H.X.; Jing, R.L.; Ma, G. Genetic basis of traits related to stomatal conductance in wheat cultivars in response to drought stress. Photosynthetica 2015, 53, 299–305. [Google Scholar] [CrossRef]
  42. Blanco, A.; Simeone, R.; Gadaleta, A. Detection of QTL for grain protein content in durum wheat. Theor. Appl. Genet. 2006, 112, 1195–1204. [Google Scholar] [CrossRef] [PubMed]
  43. Patil, R.M.; Tamhankar, S.A.; Oak, M.D.; Raut, A.L.; Honrao, B.K.; Rao, V.S.; Misra, S.C. Mapping of QTL for agronomic traits and kernel characters in durum wheat (Triticum durum Desf.). Euphytica 2013, 190, 117–129. [Google Scholar] [CrossRef]
  44. Kumar, N.; Kulwal, P.L.; Balyan, H.S.; Gupta, P.K. QTL mapping for yield and yield contibuting traits in two mapping populations fo bread wheat. Mol. Breed. 2007, 19, 163–177. [Google Scholar] [CrossRef]
  45. Blanco, A.; Mangini, G.; Giancaspro, A.; Giove, S.; Colasuonno, P.; Simeone, R.; Signorile, A.; De Vita, P.; Mastrangelo, A.M.; Cattivelli, L.; et al. Relationships between grain protein content and grain yield components through quantitative trait locus analyses in a recombinant inbred line population derived from two elite durum wheat cultivars. Mol. Breed. 2012, 30, 79–92. [Google Scholar] [CrossRef]
  46. Roncallo, P.F.; Akkiraju, P.C.; Cervigni, G.L.; Echenique, V.C. QTL mapping and analysis of epistatic interactions for grain yield and yield-related traits in Triticum turgidum L. var. durum. Euphytica 2017, 213, 277. [Google Scholar] [CrossRef]
  47. Cuthbert, J.L.; Somers, D.J.; Brûlé-Babel, A.L.; Brown, P.D.; Crow, G.H. Molecular mapping of quantitative trait loci for yield and yield components in spring wheat (Triticum aestivum L.). Theo. Appl. Genet. 2008, 117, 595–608. [Google Scholar] [CrossRef]
  48. Ain, Q.; Rasheed, A.; Anwar, A.; Mahmood, T.; Imtiaz, M.; Xia, X.; He, Z.; Quraishi, U.M. Genome-wide association for grain yield under rainfed conditions in historical wheat cultivars from Pakistan. Front. Plant Sci. 2015, 6, 743. [Google Scholar] [CrossRef] [Green Version]
  49. Qaseem, M.F.; Qureshi, R.; Shaheen, H.; Shafqat, N. Genome-wide association analyses for yield and yield-related traits in bread wheat (Triticum aestivum L.) under pre-anthesis combined heat and drought stress in field conditions. PloS ONE 2019, 14, e213407. [Google Scholar] [CrossRef] [Green Version]
  50. Golabadi, M.A.; Arzani, A.; Mirmohammadi-Maibody, S.A.M. Assessment of drought tolerance in segregating populations in durum wheat. Afr. J. Agric. Res. 2006, 5, 162–171. [Google Scholar]
  51. Ballesta, P.; Freddy Mora, F.; Del-Pozo, A. Association mapping of drought tolerance indices in wheat: QTL-rich regions on chromosome 4A. Sci. Agricola 2020, 77, 2. [Google Scholar] [CrossRef] [Green Version]
  52. Golabadi, M.; Arzani, A.; Mirmohammadi-Maibody, S.A.M.; Sayed-Tabatabaei, B.E.; Mohammadi, S.A. Identification of microsatellite markers linked with yield components under drought stress at terminal growth stages in durum wheat. Euphytica 2011, 177, 207–221. [Google Scholar] [CrossRef]
  53. Wang, S.; Basten, C.J.; Zeng, Z.B. Windows QTL Cartographer 2.5. Department of Statistics, North Carolina State University, Raleigh. Available online: http://statgen.ncsu.edu/qtlcart/WQTLCart.html (accessed on 22 December 2019).
Figure 1. Distribution of QTL clusters in RIL population of durum wheat population. * = marker linked with a QTL in both seasons for any trait, ** = marker linked with the same trait in different conditions in a single season, ꙮ = marker linked with the same trait either in different conditions in both seasons. Markers linked with any yield QTL are under lined. For abbreviation, refer to Table 4.
Figure 1. Distribution of QTL clusters in RIL population of durum wheat population. * = marker linked with a QTL in both seasons for any trait, ** = marker linked with the same trait in different conditions in a single season, ꙮ = marker linked with the same trait either in different conditions in both seasons. Markers linked with any yield QTL are under lined. For abbreviation, refer to Table 4.
Ijms 21 02372 g001aIjms 21 02372 g001bIjms 21 02372 g001cIjms 21 02372 g001d
Table 1. Range and mean ± SD of traits measured in S1 and S2 including significant differences among genotypes, seasons, treatments, genotype x season, genotype x treatment, season x treatment and genotype x season x treatment interactions at p-values 0.05 (*), 0.01 (**) and 0.001 (***) in durum wheat RIL population. Gr = germination %, Tr = Transpiration rate, Sc = Stomatal conductance, Phr = Photosynthetic rate, Wue = water use efficiency, Hd = heading time, Ph = plant height, Sl = spike length, Sw = spike weight, Sps = seeds per spike, Bm = biomass, Yd = yield, HI = harvest index, Tkw = thousand kernel weight, STI = stress tolerance index, MP = mean productivity, ST = stress tolerance, SSI= stress susceptibility index, DRI = drought resistance index, S1 = season 1 (2016-17), S2 = (2017-18), D = drought, R = relative and NA = not applicable.
Table 1. Range and mean ± SD of traits measured in S1 and S2 including significant differences among genotypes, seasons, treatments, genotype x season, genotype x treatment, season x treatment and genotype x season x treatment interactions at p-values 0.05 (*), 0.01 (**) and 0.001 (***) in durum wheat RIL population. Gr = germination %, Tr = Transpiration rate, Sc = Stomatal conductance, Phr = Photosynthetic rate, Wue = water use efficiency, Hd = heading time, Ph = plant height, Sl = spike length, Sw = spike weight, Sps = seeds per spike, Bm = biomass, Yd = yield, HI = harvest index, Tkw = thousand kernel weight, STI = stress tolerance index, MP = mean productivity, ST = stress tolerance, SSI= stress susceptibility index, DRI = drought resistance index, S1 = season 1 (2016-17), S2 = (2017-18), D = drought, R = relative and NA = not applicable.
S1S2p-Value
TraitRangeMean ± SDTraitRangeMean ± SD
Agronomic Traits
Gr_S150–10081.86 ± 9.63Gr_S280–10096.00 ± 3.72PS = ***
GrD_S147.5–10072.55 ± 9.96GrD_S270–10091.24 ± 5.56PT = ***
RGr_S166.7–10088.24 ± 11.78RGr_S273.68–10095.11 ± 5.30PG = *
PS*T = *
Hd_S180–128100.47 ± 11.79Hd_S284–12098.31 ± 7.19PT = ***
HdD_S180–12693.94 ± 12.41HdD_S278.5–114.592.32 ± 6.79PS*T = *
RHd_S151–121.7894.30 ± 3.38RHd_S280–101.6694.18 ± 2.70
Ht_S161–148114.78 ± 23.82Ht_S262–185111.36 ± 26.31 PS = ***
HtD_S153–11080.24 ± 12.45HtD_S251–109.8383.65 ± 14.14PT = ***
RHt_S145.07–99.174.14 ± 10.05RHt_S256.15–98.0376.51 ± 8.94
Physiological Traits
Tr_S122.55–56.0237.38 ± 6.86Tr_S220.02–48.2534.17 ± 6.51PT = ***
TrD_S15.17–27.815.56 ± 4.65TrD_S28.59–19.6815.06 ± 1.90PG = **
RTr_S114.2–95.542.98 ± 15.98RTr_S221.28-73.2445.43 ± 9.78PT*G = *
Sc_S1675.08–3743.61492.70 ± 527.26Sc_S2526.72–2059.51143.90 ± 293.57PS = ***
ScD_S1118.33–675.08342.96 ± 114.03ScD_S299.85–486.87262.53 ± 93.35PT = ***
RSc_S13.16–61.225.29 ± 11.39RSc_S210.14–47.1923.50 ± 8.17PG = ***
PS*T = **
PT*G = **
PS*T = *
Phr_S12.5–10.23.97 ± 1.24Phr_S21.79–8.63.49 ± 1.25PS = ***
PhrD_S11.25–3.842.22 ± 0.44PhrD_S20.57–3.781.72 ± 0.67PT = ***
RPhr_S114.93–9358.96 ± 15.51RPhr_S215.95–97.8150.94 ± 15.76PG = *
PS*T = **
Wue_S17.32–30.5911.89 ± 4.03Wue_S24.41–36.1211.14 ± 6.30PS = **
WueD_S15.95–54.6316.04 ± 7.93WueD_S23.29–31.3311.43 ± 5.14PT = ***
RWue_S130.56–659.29144.41 ± 87.62RWue_S224.38–225.63113.98 ± 41.44PS*T = **
Spike Related Traits
Sl_S16.4–12.068.25 ± 0.96Sl_S25.8–9.967.60 ± 0.71PS = *
SlD_S14.43–9.587.13 ± 1.11SlD_S24.4–96.86 ± 0.66PT = ***
RSl_S156.51–99.4386.71 ± 11.40RSl_S274.48–10090.64 ± 6.81PS*T = ***
PT*G = *
Sps_S124.6–62.643.35 ± 9.60Sps_S232.8–91.347.15 ± 8.18PT = ***
SpsD_S116.6–56.633.44 ± 7.75SpsD_S215.8–48.135.52 ± 5.73
RSps_S143.24–99.4478.75 ± 14.88RSps_S234.35–99.7577.17 ± 13.78
Sw_S11.45–4.943.05 ± 0.65Sw_S21.58–4.962.92 ± 0.65PT = ***
SwD_S11.19–3.762.30 ± 0.57SwD_S21.28–3.662.25 ± 0.41
RSw_S142.41–10077.19 ± 15.48RSw_S243.2–99.6678.98 ± 14.60
Tkw_S127.52–61.9448.88 ± 5.30Tkw_S231.09–56.442.08 ± 6.26PT = ***
TkwD_S126.4–55.2145.07 ± 4.95TkwD_S226.53–54.8738.99 ± 5.57PG = **
RTkw_S167.96–99.9192.37 ± 6.66RTkw_S269.54–99.9392.20 ± 5.66
Yield and Related Traits
Bm_S11180.8–6266.53535.80 ± 1086.20Bm_S21120–90003551.20 ± 1172.60PT = ***
BmD_S1453.6–2282.61070.50 ± 349.64BmD_S2420–23201191.40 ± 373.78PS*G= *
RBm_S110.03–92.0132.15 ± 12.00RBm_S212–62.4234.60 ± 11.05
Yd_S1305.68–1882.71033.10 ± 361.77Yd_S2306–18211010.10 ± 335.93PT = ***
YdD_S1120.73–639.19350.92 ± 109.75YdD_S2138–792332.46 ± 116.73PG = *
RYd_S112.5–70.9336.39 ± 12.21RYd_S211.92–67.0834.89 ± 12.27PS*G= ***
PS*T*G=**
HI_S111.58–66.630.50 ± 10.09HI_S215.65–46.2728.56 ± 6.38PT = ***
HID_S114.16–60.4933.59 ± 7.96HID_S212.38–4428.47 ± 6.68PG = ***
RHI_S140.55–219.08117.66 ± 36.37RHI_S241.31–250.31103.30 ± 30.84PT*G =***
Stress Indices
STI_S10.06–1.040.36 ± 0.20STI_S20.06–1.410.35 ± 0.21NA
MP_S1266.51–1236.5693.77 ± 213.97MP_S2258–1306.5674.98 ± 198.37NA
ST_S1158.19–1509685.70 ± 322.29ST_S2120–1401684.24 ± 297.06NA
SSI_S10.44–1.320.96 ± 0.18SSI_S20.49–1.310.97 ± 0.18NA
DRI_S149.91–1255.7389.18 ± 208.56DRI_S249.97–1211.8371.11 ± 229.50NA
Table 2. Composite interval mapping analysis of the durum wheat RIL mapping population. Highly significant LODs are in bold. QTL with similar superscripts are likely identical loci.
Table 2. Composite interval mapping analysis of the durum wheat RIL mapping population. Highly significant LODs are in bold. QTL with similar superscripts are likely identical loci.
TraitQTL DesignationChrPosLODClosest MarkerAdd. EffectR2Donor
Agronomic Traits
Gr_S2QGr_S2.NIAB-1B.11B19.813.75Xgwm752−1.750.18B2
Gr_S2QGr_S2.NIAB-1B.21B31.015.41Xbarc8−1.770.18B2
Gr_S2QGr_S2.NIAB-3As3A96.912.23Xbarc3561.050.07O5
GrD_S2QGrD_S2.NIAB-3Bt3B11.912.51Xgwm2851.590.07O5
GrD_S1QGrD_S1.NIAB-4Baa4B36.512.04Xgwm710b−2.640.06B2
Gr_S1QGr_S1.NIAB-5A.15A73.712.51Xbarc3423.060.08O5
Gr_S1QGr_S1.NIAB-5A.2hh5A171.813.59Xgwm1171b−4.270.18B2
GrD_S1QGrD_S1.NIAB-5B5B48.712.81Xgwm1165−3.310.1B2
GrD_S2QGrD_S2.NIAB-5B.1nn5B207.512.6Xbarc266−1.710.08B2
RGr_S2QRGr_S2.NIAB-5B.1nn5B207.512.87Xbarc266−2.020.1B2
GrD_S2QGrD_S2.NIAB-5B.25B247.112.95Xbarc59−1.690.09B2
GrD_S2QGrD_S2.NIAB-5B.3pp5B256.012.03Xgwm790−1.40.06B2
RGr_S2QRGr_S2.NIAB-5B.2pp5B261.012.31Xgwm790−1.760.09B2
Gr_S2QGr_S2.NIAB-6B6A121.513.09Xgwm49151.440.14O5
Gr_S1QGr_S1.NIAB-6B.1rr6B141.412.9Xbarc178−3.640.12B2
Gr_S1QGr_S1.NIAB-6B.26B147.113.02Xbarc247−3.530.11B2
RGr_S1QRGr_S1.NIAB-7B7B216.212.25Xgwm983−3.690.08B2
HdD_S1QHdD_S1.NIAB-2A2A5.812.21Xgwm1244−4.230.08B2
Hd_S1QHd_S1.NIAB-3Aq3A73.514.54Xwmc50−5.110.17B2
HdD_S1QHdD_S1.NIAB-3A3A44.512.1Xgwm779−4.470.09B2
Hd_S1QHd_S1.NIAB-5Agg5A169.812.52Xbarc165−4.360.12B2
HdD_S2QHdD_S2.NIAB-5Ahh5A188.013.39Xgwm1171b−3.020.18B2
RHd_S1QRHd_S1.NIAB-6B6B123.62.35Xgwm963a2.50.1O5
RHd_S1QRHd_S1.NIAB-7B7B42.54.35Xgwm1954.020.16O5
RHd_S1QRHd_S1.NIAB-7Bccc7B223.82.05Xwmc396−2.30.08B2
PhD_S2QPhD_S2.NIAB-2Ah2A121.312.04Xgwm713.640.06O5
Ph_S1QPh_S1.NIAB-2Bl2B19.712.18Xwmc5976.880.07O5
PhD_S1QPhD_S1.NIAB-3A3A14.612.86Xbarc294−3.990.09B2
PhD_S1QPhD_S1.NIAB-2Bm2B52.812.29Xgwm972−3.410.06B2
Ph_S2QRPh_S2.NIAB.3A3A13.012.21Xbarc572.390.06B2
Ph_S2QRPh_S2.NIAB.4By4B9.912.8Xbarc193−2.660.08B2
PhD_S2QPhD_S2.NIAB-4B.1y4B14.912.92Xbarc1934.780.1O5
PhD_S1QPhD_S1.NIAB-4Bz4B15.915.2Xgwm9256.150.22O5
Ph_S2QPh_S2.NIAB-4B.1z4B24.012.21Xgwm9258.370.09O5
PhD_S2QPhD_S2.NIAB-4B.2z4B25.012.51Xgwm9254.690.1O5
Ph_S2QPh_S2.NIAB-4B.2aa4B35.712.99Xgwm710b8.580.1O5
Ph_S1QPh_S1.NIAB-5Aee5A46.912Xgwm1296.040.06O5
PhD_S2QPhD_S2.NIAB-5Aii5A277.712.3Xgwm126−4.360.09B2
Ph_S2QPh_S2.NIAB-5Ajj5A334.412.42Xbarc261−7.910.08B2
RPh_S1QRPh_S1.NIAB-5A5A311.713.3Xgwm9953.680.1O5
Ph_S1QPh_S1.NIAB-6A6A51.313.16Xgwm4608−8.310.1B2
PhD_S2QPhD_S2.NIAB-7A7A144.212.04Xwmc4883.680.06O5
RPh_S1QRPh_S1.NIAB-7Bbbb7B252.612.92Xgwm1123.840.09O5
RPh_S2QRPh_S2.NIAB.7B7B288.814.78Xwmc5174.730.16O5
Physiological Traits
Tr_S1QTr_S1.NIAB-2A2A189.712.8Xgwm382−2.410.11B2
TrD_S1QTrD_S1.NIAB-2Bk2B0.014.13Xbarc1831.90.12O5
RTr_S1QRTr_S1.NIAB-3Aq3A73.513.49Xwmc505.870.13O5
TrD_S1QTrD_S1.NIAB-3Ar3A77.513.48Xgwm51.740.13O5
RTr_S2QRTr_S2.NIAB.3B3B152.312.2Xgwm655−11.610.14B2
Tr_S2QTr_S2.NIAB.5Aii5A274.712.41Xgwm126−1.10.07B2
TrD_S1QTrD_S1.NIAB-5Bkk5B13.412.24Xgwm540a−1.240.06B2
RTr_S2QRTr_S2.NIAB.5Bmm5B158.912.03Xgwm408−8.410.06B2
Tr_S2QTr_S2.NIAB.5Bmm5B158.112.9Xgwm4081.340.09O5
Tr_S2QTr_S2.NIAB.6Aqq6A246.212.55Xwmc6211.450.12O5
TrD_S1QTrD_S1.NIAB-7B7B51.113.07Xgwm4001.680.11O5
TrD_S1QTrD_S1.NIAB-7B.2aaa7B219.112.48Xbarc276−1.310.07B2
RTr_S1QRTr_S1.NIAB-7Bccc7B238.912.19Xwmc3965.650.11O5
Tr_S1QTr_S1.NIAB-7Bbbb7B248.913.03Xgwm112−2.690.14B2
Sc_S1QSc_S1.NIAB-1B.11B168.612.47Xbarc61168.280.08O5
Sc_S1QSc_S1.NIAB-1B.21B192.415.03Xwmc134−245.340.17B2
ScD_S2QScD_S2.NIAB-2Bn2B107.812.17Xgwm471.810.07O5
ScD_S2QScD_S2.NIAB-3Bu3B13.912.96Xgwm685−2.180.09B2
RSc_S2QRSc_S2.NIAB-5Aii5A272.612.41Xgwm12610.030.1O5
Sc_S2QSc_S2.NIAB-5Bmm5B158.112.83Xgwm4084.290.09O5
Sc_S2QSc_S2.NIAB-6A.16A80.012.01Xbarc107−3.280.05B2
Sc_S2QSc_S2.NIAB-6A.2qq6A244.212.44Xwmc6214.880.13O5
ScD_S1QScD_S1.NIAB-7Att7A0.012.92Xgwm1171a37.930.09O5
RPhr_S2QRPhr_S2.NIAB-1Bd1B140.813.65Xwmc1566.260.15O5
PhrD_S2QPhrD_S2.NIAB-2A.1h2A121.313.15Xgwm710.140.1O5
PhrD_S2QPhrD_S2.NIAB-2A.2j2A175.215.41Xgwm1256−0.250.24B2
RPhr_S2QRPhr_S2.NIAB-2Bl2B18.712.25Xwmc5974.630.07O5
RPhr_S2QRPhr_S2.NIAB-3A.1o3A17.613.68Xgwm757−5.440.12B2
Phr_S2QPhr_S2.NIAB-3Ar3A77.514.5Xgwm50.120.17O5
RPhr_S2QRPhr_S2.NIAB-3A.23A108.012.36Xgwm666−4.580.07B2
PhrD_S2QPhrD_S2.NIAB-3B.1w3B181.813.75Xgwm299−0.150.13B2
PhrD_S2QPhrD_S2.NIAB-3B.2w3B186.812.82Xgwm299−0.140.1B2
Phr_S1QPhr_S1.NIAB-4Bbb4B43.312.34Xgwm1167−0.340.07B2
Phr_S2QPhr_S2.NIAB-5B5B6.414.5Xgwm1910.090.14O5
Phr_S1QPhr_S1.NIAB-6A6A2.012.15Xgwm459−0.310.06B2
Phr_S1QPhr_S1.NIAB-6B6B0.012.47Xgwm940b−0.370.08B2
RPhr_S2QRPhr_S2.NIAB-6B6B60.812.09Xbarc1368.310.28O5
Phr_S2QPhr_S2.NIAB-7A7A186.413.14Xgwm276−0.090.13B2
PhrD_S2QPhrD_S2.NIAB-7B7B106.812.06Xwmc1820.120.06O5
RPhr_S2QRPhr_S2.NIAB-7Bxx7B150.612.31Xgwm46−4.710.08B2
WueD_S2QWueD_S2.NIAB-2Aj2A166.212.33Xgwm1256−1.560.19B2
WueD_S1QWueD_S1.NIAB-2Bk2B0.013.71Xbarc183−0.030.11B2
WueD_S1QWueD_S1.NIAB-4B4B31.712.6Xbarc199−0.020.08B2
Wue_S1QWue_S1.NIAB-4B.1aa4B36.513.06Xgwm710b−0.010.1B2
Wue_S1QWue_S1.NIAB-4B.2bb4B43.312.12Xgwm1167−0.010.07B2
WueD_S2QWueD_S2.NIAB-5A5A259.014.5Xwmc7271.880.22O5
Wue_S2QWue_S2.NIAB-5Aii5A269.612Xgwm1260.930.08O5
Wue_S2QWue_S2.NIAB-5Bmm5B158.912.12Xgwm408−0.880.06B2
Wue_S1QWue_S1.NIAB-7A7A18.012.69Xwmc283−0.010.09B2
WueD_S2QWueD_S2.NIAB-7Bxx7B149.612Xgwm46−0.860.06B2
Wue_S1QWue_S1.NIAB-7B7B329.612.08Xbarc320.010.07O5
Spike Related Traits
Sl_S1QSl_S1.NIAB-2Ag2A76.416.89Xgwm1220.460.21O5
SlD_S1QSlD_S1.NIAB-2B2B33.214.23Xbarc55−0.570.16B2
RSl_S1QRSl_S1.NIAB-3Aq3A61.112.1Xwmc50−5.640.13B2
SlD_S2QSlD_S2.NIAB-4Bcc4B151.812.34Xgwm940a−0.180.08B2
Sl_S1QSl_S1.NIAB-4B4B154.813.8Xgwm935a−0.330.19B2
Sl_S2QSl_S2.NIAB-5B5B33.313.04Xgwm1180−0.240.1B2
SlD_S2QSlD_S2.NIAB-7Avv7A217.812.13Xgwm332−0.170.07B2
Sl_S2QSl_S2.NIAB-7A7A278.214.05Xgwm1061−0.270.14B2
RSl_S2QRSl_S2.NIAB-7Bww7B89.913.86Xgwm5733.110.18O5
RSl_S1QRSl_S1.NIAB-7B7B117.612.38Xgwm1184−5.050.1B2
SwD_S1QSwD_S1.NIAB-1Aa1A27.112.71Xwmc3120.180.09O5
Sw_S1QSw_S1.NIAB-1A1A30.115.25Xwmc93a0.290.19O5
SwD_S1QSwD_S1.NIAB-2Ag2A76.412.22Xgwm1220.140.06O5
RSw_S1QRSw_S1.NIAB-2Ag2A76.413.38Xgwm1225.10.1O5
SwD_S1QSwD_S1.NIAB-3Aq3A68.512.55Xwmc50−0.160.08B2
RSw_S1QRSw_S1.NIAB-3Aq3A72.513.02Xwmc50−5.640.12B2
Sw_S1QSw_S1.NIAB-3As3A97.812.62Xbarc3560.20.08O5
SwD_S1QSwD_S1.NIAB-3Bu3B13.912.82Xgwm6850.180.08O5
RSw_S1QRSw_S1.NIAB-3Bw3B194.812.28Xgwm299−5.650.12B2
SwD_S1QSwD_S1.NIAB-5Add5A0.012.94Xgwm1540.180.09O5
SwD_S2QSwD_S2.NIAB-5Aii5A277.712.14Xgwm126−0.120.09B2
RSw_S1QRSw_S1.NIAB-5Ajj5A337.413.1Xbarc261−5.740.13B2
RSw_S2QRSw_S2.NIAB-5B.1oo5B237.112.42Xgwm10724.680.09O5
RSw_S2QRSw_S2.NIAB-5B.2pp5B260.013.65Xgwm7905.870.15O5
SwD_S1QSwD_S1.NIAB-5B5B270.912.23Xgwm6050.140.06O5
SwD_S2QSwD_S2.NIAB-7Auu7A57.812Xwmc405−0.120.08B2
RSw_S1QRSw_S1.NIAB-7A7A98.812.59Xgwm710a−4.370.07B2
RSw_S1QRSw_S1.NIAB-7Bww7B92.213.09Xgwm5735.360.09O5
SpsD_S1QSpsD_S1.NIAB-1Aa1A25.112.41Xwmc3122.110.06O5
RSps_S1QRSps_S1.NIAB-1Bc1B104.612.12Xwmc419−4.030.06B2
Sps_S1QSps_S1.NIAB-1B.1d1B143.812.5Xwmc1563.220.11O5
Sps_S1QSps_S1.NIAB-1B.21B154.612.75Xwmc5482.910.08O5
Sps_S2QSps_S2.NIAB-3Ap3A56.112.5Xbarc452.630.1O5
SpsD_S1QSpsD_S1.NIAB-3Aq3A68.513.6Xwmc50−2.60.1B2
Sps_S1QSps_S1.NIAB-3B3B8.312Xgwm6252.40.06O5
SpsD_S1QSpsD_S1.NIAB-3Bu3B13.913.61Xgwm6852.670.1O5
RSps_S1QRSps_S1.NIAB-3Bw3B182.812.63Xgwm2994.510.08O5
Sps_S2QSps_S2.NIAB-4B4B155.912.94Xwmc428−3.290.09B2
SpsD_S1QSpsD_S1.NIAB-5Add5A1.012.84Xgwm1542.430.08O5
Sps_S1QSps_S1.NIAB-5B.1kk5B13.412.3Xgwm540a−3.340.07B2
RSps_S1QRSps_S1.NIAB-5Bll5B132.214.41Xgwm10438.720.33O5
RSps_S2QRSps_S2.NIAB-5B.2ll5B143.212.14Xgwm10433.90.07O5
SpsD_S2QSpsD_S2.NIAB-5B5B218.313.22Xbarc23220.11O5
RSps_S2QRSps_S2.NIAB-5B.2oo5B232.712.83Xgwm10724.410.09O5
Sps_S2QSps_S2.NIAB-5Boo5B233.112.27Xgwm1072−2.170.07B2
Sps_S1QSps_S1.NIAB-5B.25B273.912.85Xbarc2432.90.09O5
SpsD_S2QSpsD_S2.NIAB-7Byy7B172.612.1Xwmc218−2.360.07B2
RTkw_S1QRTkw_S1.NIAB-1Bc1B102.612.05Xwmc419−2.010.08B2
TkwD_S1QTkwD_S1.NIAB-1B.11B221.112.26Xgwm2681.560.09O5
TkwD_S1QTkwD_S1.NIAB-1B.21B230.312.21Xgwm1531.320.06O5
Tkw_S1QTkw_S1.NIAB-2Af2A17.512.11Xgwm7261.880.08O5
RTkw_S2QRTkw_S2.NIAB-4B.1aa4B35.713.98Xgwm710b−1.930.11B2
Tkw_S2QTkw_S2.NIAB-4Baa4B35.712.74Xgwm710b1.980.09O5
RTkw_S2QRTkw_S2.NIAB-4B.2bb4B42.514.24Xgwm1167−1.980.12B2
RTkw_S1QRTkw_S1.NIAB-5Aff5A124.713.22Xgwm1236−2.510.12B2
TkwD_S1QTkwD_S1.NIAB-5Aff5A126.712Xgwm1236−1.380.07B2
TkwD_S2QTkwD_S2.NIAB-5Aii5A277.712.62Xgwm126−1.850.1B2
RTkw_S2QRTkw_S2.NIAB-5B5B62.012.01Xbarc128−1.380.05B2
Tkw_S1QTkw_S1.NIAB-6Brr6B144.112.55Xbarc1781.860.08O5
TkwD_S1QTkwD_S1.NIAB-6Bss6B160.612.74Xgwm8891.610.09O5
RTkw_S2QRTkw_S2.NIAB-7A7A113.412.14Xwmc603−1.350.05B2
Tkw_S1QTkw_S1.NIAB-7Avv7A220.812.93Xgwm3322.030.1O5
RTkw_S1QRTkw_S1.NIAB-7B7B275.212.66Xbarc258−2.090.09B2
Yield and Related Traits
RBm_S1QRBm_S1.NIAB-3Ao3A16.612.32Xgwm757−1.290.06B2
RBm_S1QRBm_S1.NIAB-3Bt3B11.912.06Xgwm285−1.220.05B2
RBm_S2QRBm_S2.NIAB-3Bv3B88.112.85Xbarc3443.880.1O5
Bm_S2QBm_S2.NIAB-3Bv3B89.112.22Xbarc344−349.50.07B2
BmD_S1QBmD_S1.NIAB-4B4B88.313.56Xgwm998130.010.11O5
RBm_S1QRBm_S1.NIAB-5Aee5A46.912.4Xgwm129−1.390.07B2
Bm_S2QBm_S2.NIAB-5Boo5B232.712.65Xgwm1072−389.160.08B2
RBm_S2QRBm_S2.NIAB-5Boo5B232.712.03Xgwm10723.040.07O5
RBm_S1QRBm_S1.NIAB-5Bpp5B259.013.45Xgwm790−1.780.11B2
RBm_S1QRBm_S1.NIAB-7A7A213.912.83Xbarc253−1.460.08B2
RBm_S2QRBm_S2.NIAB-7A7A91.512.21Xbarc1025−3.380.08B2
BmD_S2QBmD_S2.NIAB-7Auu7A52.812.73Xwmc405−150.390.13B2
Bm_S1QBm_S1.NIAB-7B.17B137.513.06Xgwm540b−558.270.1B2
Bm_S1QBm_S1.NIAB-7B.27B164.413.16Xwmc476640.830.13O5
RBm_S2QRBm_S2.NIAB-7B.1yy7B172.613.55Xwmc218−4.240.12B2
RBm_S2QRBm_S2.NIAB-7B.27B179.413.39Xgwm963b−4.140.11B2
RBm_S2QRBm_S2.NIAB-7B.3zz7B184.412.43Xgwm1085−3.750.08B2
Bm_S1QBm_S1.NIAB-7B.3aaa7B218.212Xbarc276−334.540.07B2
BmD_S2QBmD_S2.NIAB-7Baaa7B221.112Xbarc276−113.280.07B2
RYd_S2QRYd_S2.NIAB-1Bb1B95.612.07Xgwm7843.610.06O5
YdD_S1QYdD_S1.NIAB-2Af2A26.012.65Xgwm72637.770.09O5
YdD_S2QYdD_S2.NIAB-2Ai2A124.612.51Xgwm89534.890.08O5
RYd_S1QRYd_S1.NIAB-2Ai2A125.212.44Xgwm895−4.070.07B2
Yd_S2QYd_S2.NIAB-3Ap3A52.112.38Xbarc4596.110.07O5
RYd_S2QRYd_S2.NIAB-3Bx3B205.713.58Xgwm547−4.910.11B2
Yd_S2QYd_S2.NIAB-3Bx3B205.712.2Xgwm54791.540.06O5
Yd_S1QYd_S1.NIAB-5Agg5A158.812.09Xbarc16599.370.07O5
RYd_S2QRYd_S2.NIAB-5Ahh5A193.014.04Xgwm1171b−7.680.3B2
YdD_S1QYdD_S1.NIAB-5Bll5B143.212.66Xgwm104338.450.09O5
YdD_S1QYdD_S1.NIAB-6B6B9.912.29Xgwm1199−35.550.08B2
RYd_S1QRYd_S1.NIAB-6B6B121.712.09Xgwm963a−3.850.07B2
Yd_S1QYd_S1.NIAB-7Baaa7B218.212.52Xbarc276−117.780.08B2
YdD_S2QYdD_S2.NIAB-7Baaa7B219.113Xbarc276−39.230.09B2
RHI_S2QRHI_S2.NIAB-1Bb1B93.112.22Xgwm7849.70.07O5
RHI_S2QRHI_S2.NIAB-2Bn2B108.812.4Xgwm47−10.10.08B2
RHI_S2QRHI_S2.NIAB-3B.1v3B89.812.91Xbarc344−10.710.09B2
HI_S2QHI_S2.NIAB-3Bx3B205.712.42Xgwm5472.030.08O5
RHI_S2QRHI_S2.NIAB-3B.2x3B205.712.8Xgwm547−10.170.09B2
HID_S2QHID_S2.NIAB-4B4B130.612.5Xbarc602.020.08O5
HID_S1QHID_S1.NIAB-4Bcc4B151.212.59Xgwm940a3.310.09O5
HI_S1QHI_S1.NIAB-5B5B228.612.77Xwmc283.190.1O5
HID_S2QHID_S2.NIAB-5Bll5B143.213.3Xgwm10432.240.11O5
HID_S2QHID_S2.NIAB-7A7A108.612.42Xbarc1081.80.07O5
RHI_S2QRHI_S2.NIAB-7Bzz7B184.412.91Xgwm108511.250.09O5
HI_S1QHI_S1.NIAB-7B7B336.012.42Xwmc2732.950.08O5
Stress Indices
STI_S2QSTI_S2.NIAB-4A.14A44.512.73Xbarc2460.060.09O5
STI_S2QSTI_S2.NIAB-4A.24A100.812.13Xgwm1234−0.060.06B2
STI_S1QSIT_S1.NIAB-6Bss6B159.912.14Xgwm8890.060.07O5
STI_S2QSTI_S2.NIAB-7Att7A0.013.36Xgwm1171a−0.070.11B2
MP_S1QMP_S1.NIAB-1B1B72.012.28Xgwm76287.360.08O5
MP_S2QMP_S2.NIAB-1B1B116.812.21Xbarc207−61.490.08B2
MP_S2QMP_S2.NIAB-2Aj2A158.212.04Xgwm125697.140.22O5
MP_S2QMP_S2.NIAB-3B.1w3B194.812.87Xgwm29982.740.16O5
MP_S2QMP_S2.NIAB-3B.2x3B205.713.72Xgwm54775.530.13O5
MP_S1QMP_S1.NIAB-5Agg5A158.812.4Xbarc16563.780.08O5
ST_S2QST_S2.NIAB-3Bw3B192.812.22Xgwm29999.60.11O5
ST_S1QST_S1.NIAB-5Agg5A165.812.12Xbarc165112.920.11O5
ST_S2QST_S2.NIAB-5Ahh5A179.813.43Xgwm1171b116.020.13O5
DRI_S2QDRI_S2.NIAB-2Bm2B52.212.02Xgwm97259.70.06O5
DRI_S2QDRI_S2.NIAB-4A4A64.312.47Xbarc34382.630.1O5
DRI_S2QDRI_S2.NIAB-5Ahh5A186.012.59Xgwm1171b−81.060.11B2
DRI_S1QDRI_S1.NIAB-5Ahh5A188.014.63Xgwm1171b−111.010.26B2
DRI_S1QDRI_S1.NIAB-5Boo5B232.712.44Xgwm107263.650.08O5
Table 3. Percentage increase or decrease in yield under various conditions with respect to Xgwm895 at 2A and Xbarc276 at 7B.
Table 3. Percentage increase or decrease in yield under various conditions with respect to Xgwm895 at 2A and Xbarc276 at 7B.
Trait Xgwm895, 2AXbarc276, 7B Combination
A alleleB alleleA alleleB alleleAABBAB
Yd_S15.5−9.315.56.6−4−1.892.16
YdD_S10.18−1.60.180.660.61−1.83−0.15
RYd_S1−6.656.82−6.65−4.661.562.26−14.1
Yd_S26.16−9.356.162.570.92−10.7810.2
YdD_S25.27−8.45.274.3922.98−10.4311.3
RYd_S2−1.591.85−1.592.61−3.044.74−0.39
Table 4. Absolute and relative traits measured for RIL population in 2016–17 (trait abbreviation followed by _S1) and 2017–18 (trait abbreviation followed by _S2) in control and drought and their time of measurement each season.
Table 4. Absolute and relative traits measured for RIL population in 2016–17 (trait abbreviation followed by _S1) and 2017–18 (trait abbreviation followed by _S2) in control and drought and their time of measurement each season.
TraitsControlDroughtRelativeTime
Agronomic Traits
1GerminationGr_S1, Gr_S2GrD_S1, GrD_S2RGr_S1, RGr_S224th Nov
2Heading dateHd_S1, Hd_S2HdD_S1, HdD_S2RHd_S1, RHd_S250% spike emerg
3Plant heightPh_S1, Ph_S2PhD_S1, PhD_S2RPh_S1, RPh_S2 Before harvest
Physiological Traits
4Transpiration rateTr_S1, Tr_S2TrD_S1, TrD_S2RTr_S1, RTr_S17th Feb
5Stomatal conductanceSc_S1, Sc_S2ScD_S1, ScD_S2RSc_S1, RSc_S27th Feb
6Photosynthetic ratePhr_S1, Phr_S2PhrD_S1, PhrD_S2RPhr_S1, RPhr_S27th Feb
7Water use efficiencyWue_S1, Wue_S2WueD_S1, WueD_S2RWue_S1, RWue_S2 7th Feb
Spike-Related Traits
8Spike lengthSl_S1, Sl_S2SlD_S1, SlD_S2RSl_S1, RSl_S2After harvest
9Spike weightSw_S1, Sw_S2SwD_S1, SwD_S2RSw_S1, RSw_S2After harvest
10Seeds per spikeSps_S1, Sps_S2SpsD_S1, SpsD_S2RSps_S1, RSps_S2After harvest
11Thousand kernel weightTkw_S1, Tkw_S2TkwD_S1, TkwD_S2RTkw_S1, RTkw_S2After harvest
TraitsControlDroughtTime
Yield and Related Traits
12BiomassBm_S1, Bm_S2BmD_S1, BmD_S2After harvest
13YieldYd_S1, Yd_S2YdD_S1, YdD_S2After harvest
14Harvest IndexHI_S1, HI_S2HID_S1, HID_S2After harvest
IndexAbbrev.FormulaTime
Stress Indices
15Stress Tolerance Index (STI)STI_S1, STI_S2 ( Y p × Y s ) ( Y p ¯ ) 2 After harvest
16Mean Productivity (MP)MP_S1, MP_S2 ( Y p + Y s ) 2 After harvest
17Stress Tolerance (ST)ST_S1, ST_S2 Y p Y s After harvest
18Stress Susceptibility Index (SSI)SSI_S1, SSI_S2 [ 1 ( Y s Y p ) ] / [ 1 ( Y s ¯ Y p ¯ ) ] After harvest
19Drought Resistance Index (DRI)DRI_S1, DR_S2 Y s × ( Y s Y p ) / Y s ¯ After harvest
where Ys = yield under stress, Y s ¯ = mean yield under stress, Yp = yield under control and Y p ¯ = mean yield under control

Share and Cite

MDPI and ACS Style

Rehman Arif, M.A.; Attaria, F.; Shokat, S.; Akram, S.; Waheed, M.Q.; Arif, A.; Börner, A. Mapping of QTLs Associated with Yield and Yield Related Traits in Durum Wheat (Triticum durum Desf.) Under Irrigated and Drought Conditions. Int. J. Mol. Sci. 2020, 21, 2372. https://doi.org/10.3390/ijms21072372

AMA Style

Rehman Arif MA, Attaria F, Shokat S, Akram S, Waheed MQ, Arif A, Börner A. Mapping of QTLs Associated with Yield and Yield Related Traits in Durum Wheat (Triticum durum Desf.) Under Irrigated and Drought Conditions. International Journal of Molecular Sciences. 2020; 21(7):2372. https://doi.org/10.3390/ijms21072372

Chicago/Turabian Style

Rehman Arif, Mian Abdur, Fauzia Attaria, Sajid Shokat, Saba Akram, Muhammad Qandeel Waheed, Anjuman Arif, and Andreas Börner. 2020. "Mapping of QTLs Associated with Yield and Yield Related Traits in Durum Wheat (Triticum durum Desf.) Under Irrigated and Drought Conditions" International Journal of Molecular Sciences 21, no. 7: 2372. https://doi.org/10.3390/ijms21072372

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