Introduction

With few exceptions (Richards 2006), breeding for drought adaptation has been largely empirical or based on drought escape through manipulation of development, so that sensitive development stages do not coincide with the greatest risk and effect of stress (Ludlow and Muchow 1990). Although conventional breeding has had substantial impact in marginal wheat growing environments (Lantican et al. 2003), future genetic gains will require a more systematic use of physiological and genetic approaches, facilitated by the rapid increase in genome knowledge and understanding. An example of physiological application is how the ability to use water more slowly and efficiently has improved performance of wheat cultivars under late-season drought conditions of Australia, via the crossing and selecting for high transpiration efficiency (TE) (Condon et al. 2004). In situations where water is available in deep soil profiles, the ability to extract it under water-limited conditions has been identified as the principal mechanism for drought adaptation as shown by a robust association of canopy temperature (CT) with wheat yields (Olivares-Villegas et al. 2007). Traits such as these that are rapid and integrative in nature are increasingly recognized not only as useful selection tools in breeding (Araus et al. 2008) but as valuable screens for high throughput phenotyping of mapping populations (Montes et al. 2007; Olivares-Villegas et al. 2008). In addition to CT, several other physiological traits can be either rapidly phenotyped and/or informative about how adaptation to drought and heat can arise. Normalized difference vegetative index (NDVI) is an integrated measure of both ground cover (leaf area) and the nitrogen content of the canopy. Both NDVI and chlorophyll meters provide an indirect estimation of leaf health for photosynthesis as associated with leaf nitrogen content and early senescence (Araus et al. 2008). Given similar canopy leaf area, the maintenance of healthy green leaves under stress can be interpreted as stress tolerance for either drought or heat conditions (Olivares-Villegas et al. 2007; Tao et al. 2000).

Mapping quantitative trait loci (QTL) in bi-parental populations allows the detection of chromosome segments controlling traits of agronomic interest with the opportunity to dissect complex traits (those traits that integrate the effect of a number of factors occurring during the crop cycle) into component loci (Marza et al. 2006). Detected regions can be used in breeding and pre-breeding to increase its efficiency through marker-assisted selection (Campbell et al. 2003). In both maize and soybean, new cultivars have been released to markets following the rapid recombination of favorable yield QTL through marker-assisted recurrent selection strategies based on multi-environment assessments of bi-parental populations, e.g. Eathington et al. (2007). While effective in generally favorable production environments, this approach is potentially confounded for traits under stress due to the existence of interactions of multiple physiological processes and environmental influences over the period of the crop cycle.

In recent years, large numbers of QTL have been reported in diverse cereals for a range of agronomic traits: for example, in barley, QTL have been reported for yield under drought environments (Comadran et al. 2008; Talamè et al. 2004), and in wheat, QTL for seed size, seed shape (Breseghello and Sorrells 2007), plant height, maturity, and grain yield (Kato et al. 2000; Kuchel et al. 2007; Marza et al. 2006; McCartney et al. 2005; Snape et al. 2007). However, relatively few studies have examined QTL for physiological traits and their co-location with effects on crop yield and quality. Exceptions are for concentration of water soluble carbohydrates (WSC) and carbon isotope discrimination (CID) (Rebetzke et al. 2008a; Rebetzke et al. 2008b; Yang et al. 2007) although these studies were partially confounded by variation in flowering time and height in the populations. Public research on QTL in cereals has often utilized greatly contrasting parents, including land-races. However, QTL associated with stress responses in such populations may be confounded by other differences, or have already been fixed in elite breeding programs.

In this paper, the population studied is an elite cross (Seri/Babax), and does not segregate for known major height (Rht genes), vernalization (Vrn1:Vrn-A1c, Vrn-B1, and Vrn-D1) or photoperiod response genes (Ppd1). Consequently, the lines have shown low variation for both flowering time and height (Olivares-Villegas et al. 2007). This potentially increases the opportunity to identify QTL that may otherwise be masked by the confounding effects of crop development. By phenotyping in both drought and hot/irrigated environments, the study aimed to confirm the existence of common genetic bases between adaptation to moisture and heat stress that have been indicated by physiological studies in the same environments (Reynolds et al. 2007). Specific objectives of this study were: (1) to identify QTL associated with agronomic and physiological traits in six environments encompassing drought, heat, and well-irrigated conditions; (2) to demonstrate the value of using progeny from an elite by elite cross expressing a restricted range of height and phenology for improving the QTL detection; (3) to determine genomic regions consistently associated with adaptation to both drought and hot, irrigated environments; and (4) to interpret ‘genetic dissection’ of yield in terms of QTL simultaneously associated with yield, yield components, and stress-adaptive traits contributing to performance.

Materials and methods

Plant material and field trials

A recombinant inbred line (RIL) population consisting of 167 sister lines was studied in managed environments in north-west Mexico. The population was derived from a reciprocal cross between the related elite lines: semi-dwarf spring wheat variety Seri M82 from the “Veery” cross (KVZ/BUHO//KAL/BB) and a fixed line (Babax) derived from the “Babax” cross (BOW/NAC//VEE/3/BJY/COC). Seri M82 carries the T1BL.1RS (rye) translocation, and is characterized by moderate tolerance to drought conditions and high yield potential. The Babax parental line has a coefficient of parentage of 0.3316 with Seri M82 (Mathews et al. 2008). It is a sister line of the elite variety Baviacora M92 (recognized for drought tolerance and also has a high yield potential), but was selected so as not to carry the T1BL.1RS rye translocation that is in present in Baviacora M92 and Seri M82 (see Mathews et al. 2008 for details). The parental lines, termed here as Seri and Babax, have been screened for known phenology alleles, and both parents have the photoperiod-insensitive allele at Ppd-D1, as well as spring-type alleles for at least two vernalization (Vrn-B1 and Vrn-D1) loci (K Cane, Department of Primary Industries, Victoria, Australia). Hence, the population is characterized by its narrow range of height and flowering time (ca. 10–15 days), and was developed primarily for genetic mapping and screening for the physiological basis of stress tolerance in drought and heat environments (Olivares-Villegas et al. 2007).

Six field trials were sown between 2002 and 2006, under three different environments: two under terminal drought (D02 and D05; total crop water supply ≤300 mm); two under high environmental temperatures due to a delayed planting date, but irrigated throughout the crop cycle (H05 and H06; average T max > 32°C; total crop water supply >700 mm); and two controls under well-irrigated conditions (I02 and I06; total crop water supply >700 mm). The trial naming protocol was: letters D, H, and I to indicate drought (DRT), heat (HOT) or irrigated (IRR) conditions, and the two-digit numbers indicate the year of harvest. During the two DRT cycles, available water (including rainfall) was estimated at ca. 195 and 260 mm for D02 and D06, respectively. For HOT and IRR treatments, irrigation was applied when approximately 50% of available soil moisture (in the 0–1 m profile) was depleted, so water was assumed not to limit growth. Available soil water for the DRT trials was estimated from gravimetric sampling of neighboring plots after irrigation. All trials were sown in two-replicate alpha-lattice designs in the Yaqui Valley, Mexico at CIMMYT’s Obregon Experimental Station, located in north-western Mexico (27°25′N 109°54′W, 38 m above sea level). The site is a high radiation, irrigated environment (Table 1). The soil is a Typic Calciorthid, low in organic matter (0.76%) and slightly alkaline (pH 7.7) with a plant-available water holding capacity of about 200 mm. Soil analyses conducted previously at various profiles did neither indicate problems associated with mineral deficiencies or toxicities, nor salinity problems (Olivares-Villegas et al. 2007). Appropriate fertilization, weed, disease, and pest control were implemented to minimize other yield limitations. Plots comprised two 80-cm raised beds spaced at 60 cm between centers and comprising two rows per bed. There were slight variations in plot lengths (l) and seeding rate (d): in D02 l = 6 m, d = 12.5 g m−2; in D05 l = 4.5 m, d = 12.5 g m−2; H05 l = 4 m, d = 11.7 g m−2; H06 l = 4 m, d = 15.6 g m−2; I02 l = 5 m, d = 15 g m−2; in I06 l = 5 m d = 10 g m−2. Olivares-Villegas et al. (2007) described the 2002 trials in detail.

Table 1 Trial mean yield, anthesis range, thousand grain weight (TGW), and heritability (h 2G ) for six trials grown in three environments: drought, heat, and irrigated

Field methods and trait calculations

The recorded traits were classified into three groups: (1) agronomic traits, (2) phenological traits and height, and (3) physiological traits. Agronomic traits included final grain yield (yield, g m−2), grain number (GM−2), thousand grain weight (TGW, g), and kernel weight index (KWI). According to the population average grain size to calculate KWI, the harvested grains from each line were sampled, and using a 19 × 3 mm sieve (to avoid a group containing only broken grains and rubbish), the sample was classified into two groups (large grains and small grains). KWI was calculated as:

$$ \left( {{\frac{\text{Lg}}{{{\text{Lg}} + {\text{Sg}}}}}} \right) \times \left( {{\frac{{{\text{TGW}}({\text{I}})_{\text{ave}} }}{{{\text{TGW}}({\text{I}})_{\text{line}} }}}} \right), $$

where Lg is the mass of large grains in the sample, Sg is the mass of small grains in the sample, TGW(I)ave is the average thousand grain weight for both irrigated trials (I02 and I06), and TGW(I)line is the average thousand grain weight for individual lines in the irrigated trials (a larger KWI indicating an increased proportion of larger grains normalized to irrigated conditions). Anthesis (ANTH), maturity (MAT), and height were included in the group of agronomic and morphological traits. ANTH was recorded as the number of days for 50% of plants to exhibit burst anthers, and MAT was recorded when 50% of the spikes in a plot showed total loss of green color (MAT). Physiological traits included canopy temperature in the vegetative stage (CTv, °C), canopy temperature in the grainfill stage (CTg, °C), normalized difference vegetative index in the vegetative stage (NDVIv), normalized difference vegetative index in the grainfill stage (NDVIg), chlorophyll content in the grainfill stage (CHLg), and water soluble carbohydrate content of stems at anthesis (WSC, %). Canopy temperature was measured using a portable infrared thermometer (Mikron M90 series) twice in the day if possible, in the morning (10:00–12:00 h) and in the afternoon (12:00–15:00) of fine windless days during vegetative (CTv) and grainfill (CTg) stages (Olivares-Villegas et al. 2007). CT was recorded on an average of four occasions during each growth stage by pointing to the canopy from the edge of the plot for approximately 5 s and located with the sun behind the observer, avoiding shadows. Each CT measurement (stage × by time of day) was statistically modeled and QTL mapped individually. Adjusted means were also averaged for each growth stage, and one data point for CTv and one for CTg was reported. NDVI was obtained using a GreenSeeker sensor (Optical Sensor Unit, 2002 Ntech Industries, Inc., Ukiah, CA, USA) with repeated measurements during the vegetative (NDVIv) and grainfill (NDVIg) stages (Araus et al. 2008; Marti et al. 2007). The instrument records the reflectance in one bed per plot at speeds of 10–20 times/plot; all zero values resulting from measuring the soil were deleted, and the average per plot was reported; the statistical analysis was also done individually for each measurement. CHLg was recorded using a portable chlorophyll meter (SPAD-502 Minolta, Spectrum Technologies Inc., Plainfield, IL, USA). During the grainfilling stage, six leaves were sampled per plot, and the average was reported; WSC was measured by sampling and drying 12 stems cut between 6 and 8 days after 50% anthesis for each plot, and estimating the percentage with Near-infrared Reflectance Spectroscopy (Foss NIRSystems mod. 6500, Uruguay) previously calibrated against chemical measurement (Ruuska et al. 2006). For logistical reasons, not all traits were collected in all experiments: in D02 observations were missing for KWI and WSC; in D05 CTv; in I02 KWI, CTv, NDVIv, NDVIg, and WSC. A detailed description of the traits measured in year 2002 can be found in Olivares-Villegas et al. (2007).

Statistical analysis of the traits data

Adjusted means (Best Linear Unbiased Estimates; BLUEs) were calculated from a spatial model using the REML procedure in Genstat Discovery Edition 3 (http://www.vsni.co.uk/). Although the Seri/Babax population has a narrow range of anthesis dates, in order to detect QTL independent from phenology, anthesis was tested as a covariate for all the traits resulting significant in the major of the physiological traits and also KWI but generally not significant for yield and yield components. The anthesis date was included as a covariate in the models of all the physiological traits plus KWI under all the environments × year combinations. To obtain variance components, first, the best spatial model was determined individually for each trait × environment × year combination, assuming random genotype effects and an autoregressive process of order 1 in both the row and columns directions (AR1 × AR1). The models included the experimental design factors for an α-lattice design: replicates and replicates × subblock. Global spatial trends in the row and column directions and extraneous spatial effects, like the direction of recording data, were also modeled following Gilmour et al. (1997). Genotype was then fitted as a fixed effect using these best spatial models to obtain the BLUEs and weights (Smith et al. 2001), which were used in the subsequent QTL analysis. A multiple linear regression was performed for yield, TGW, and GM−2 using Statgraphics Plus V4.0. Genetic correlations (\( r_{{{\text{g}}ij}} \)) between traits i and j were calculated using the genotypic variance and covariance component estimates:

$$ r_{{{\text{g}}ij}} = {\frac{{\sigma_{{{\text{G}}ij}} }}{{\sigma_{\text{G}} \sigma_{{{\text{G}}j}} }}}, $$

where σ Gij is the genotypic covariance between traits i and j, and σ Gi and σ Gj are the estimated genotypic standard deviations for trait i and j, respectively (Holland 2006).

For each trait, in each trial, the broad sense heritability was calculated based on Falconer (1981) and modified for spatial adjustments according to Cullis et al. (2006) and Oakey et al. (2006). The components of the following formula were derived from the random genotype effects model. The formula is:

$$ h_{\text{G}}^{2} = 1 - {\frac{{{\text{avesed}}^{2} }}{{2{\text{Vg}}}}}, $$

where avesed is the average standard error of the difference for the set of BLUPs and Vg is the genetic variance for the trait.

QTL mapping of the agronomic, phenological, and physiological traits

Prior to this study, a molecular map was constructed from 74 SSR (Single Sequence Repeat), 249 AFLP (Amplified Fragment Length Polymorphisms), and 264 DArT (Diversity Array Technology) markers (McIntyre et al. 2010). The final map used 401 markers to define 29 linkage groups with only chromosomes 3D and 7D missing (see Mathews et al. 2008 for further details on the map construction). Prior to QTL analysis, variance analyses were performed using Genstat Discovery Edition 3 to determine the existence of significant genetic variability between RILs. Those traits without a significant genotypic variance were not included in the QTL mapping. Composite interval mapping (CIM) was undertaken using Windows QTL Cartographer v2.0 (Wang et al. 2007) with the program’s default values: a genome 10.0 cM wide, significance level of 0.05, and a walk speed of 2 cM in forward regression with five cofactor markers. The BLUEs generated from the spatial analyses were used and the LOD (logarithm of odds ratio) scores obtained. A QTL was declared if two or more close markers (≤30 cM) were linked to a locus of significant LOD. Significant profiles were determined if the LOD score of the locus was greater than 2 (McIntyre et al. 2010). Other authors have used LOD as low as 2.5 (Sun et al. 2008; Yang et al. 2007), but the current study has an exploratory focus and hence indicative QTL signals (2–2.5) were not excluded. Significant QTL (LOD > 2) were further classified according to their appearance in more than one environment as follows: (a) repeatable QTL when it appeared only in both years of any environment; (b) stress QTL when it appeared in both years of one stress environment and at least 1 year of the other stressed environment; and (c) robust when the QTL was detected in at least 1 year of irrigated plus 2 years of an stressed environment, or vice versa. Weaker effects that appeared only in 1 year of any environment were classified as suggestive QTL; isolated LOD significant signals in only one marker but not in nearby markers (≤30 cM) were not considered as a QTL even when it appeared in more than one environment.

Results

Weather conditions under which the six trials were undertaken

Minimum and maximum temperature and total rainfall in the crop were recorded for each trial (Table 1). The temperatures were summarized into four crop stages: emergence to vegetative (EV), late vegetative (LV), anthesis (AN), and grainfill (GF) (Table 1). For IRR treatments, conditions were typical of a high radiation, irrigated, spring wheat environment with solar radiation and temperatures increasing from March through grain-filling (Table 1). The rainfall for all six trials was negligible (Table 1). Consequently, the DRT trials resulted in gradually intensifying moisture stress levels as water from seeding irrigations was depleted. For the late sown trials (HOT), gradually increasing above-optimal temperatures were experienced throughout the cycle (Table 1). Maximum temperatures in all stages were recorded in the HOT trials, as expected. In both years of the HOT environment, all the stages but early vegetative, experienced maximum temperatures greater than 30°C on more than 90% of days.

Agronomic and physiological trait screening

Mean, minimum, maximum, and standard deviation of the agronomic and physiological traits in all six trials are presented in Table 2. Yield under DRT and HOT were reduced, on average, by 50 and 59%, respectively, compared with the IRR environment. The average IRR yield was approximately 620 g m−2 with the best lines yielding almost 800 g m−2 in the best year and the worst lines 475 g m−2 in the worst year (data not shown). Compared with the IRR treatments, stress treatments generally reduced grain number (GM−2) more than grain weight (TGW), i.e. in D02, the GM−2 reduction was almost 51% while TGW was reduced by 21%. Kernel weight index (KWI) values were lower in the stressed environments, indicating a greater proportion of small and shriveled grains (Table 2).

Table 2 Summary statistics (mean, maximum, minimum, and standard deviation) for all traits measured on the Seri/Babax population in all trials

The growth cycle averaged 126 days from crop emergence to physiological maturity under IRR conditions, but was reduced, on average, by 12 days (10%) under DRT and 45 days (40%) under HOT conditions (Table 2). The relative duration of grainfilling (MAT-ANTH)/MAT, as a fraction of crop cycle length, was 0.29, 0.26, and 0.30 for DRT, HOT, and IRR, respectively. Mean canopy temperatures varied predictably according to air temperature and soil moisture while the range among genotypes within trials was 1–4°C (from the coolest to the warmest genotypes), with the range typically greater under stress compared with IRR. The spectral index NDVIv measured in the vegetative stages was generally consistent with high levels of ground cover for all sites and genotypes, with the exception of the HOT site in 2005, where incomplete ground cover resulted in low values. Stem water soluble carbohydrate (WSC) data were not collected in all environments, but the mean values were greatest in the HOT environments, and genetic variation was substantial in all environments.

The lowest yielding trials were H05 and D02 (Table 1). However, the ten best genotypes within these trials recorded 73% (H05) and 39% (D02) greater yield than the parents (Table 1); in comparison, in each of the IRR trials, the ten best genotypes yielded only 9% more than the parents (Table 1). When compared with the parents (Table 1), the ten best yielding genotypes across trials showed their highest yields under D02 and H05 (68% and 30% more than the parental mean, respectively). The maximum range of anthesis recorded was 18 days in I06 (Table 1). Heritability across trials for yield, anthesis, and TGW is presented in Table 1, and show similar values for the three traits in all of the trials.

Associations between agronomic and physiological traits

The phenotypic correlations across years of HOT, DRT, and IRR yield were 0.68, 0.66, and 0.56, respectively, while correlations for grain yield across environments varied from 0.48 to 0.74 (Fig. 1). Phenotypic and genotypic correlations between yield and the remaining traits were calculated, and the significant correlations are reported (Table 3). In terms of yield components, GM−2 showed relatively large and consistent correlations with yield in all environments, while KWI was weakly associated with yield (Table 3). There was a general trend for earliness to be associated with better performance in most environments, even within the narrow range observed here. Consistent negative correlations with yield were observed for CTv and CTg (high temperatures) in the stressed environments but not in IRR (Table 3). The spectral index (NDVIv) showed an association with yield in all environments being stronger under both stress environments. CHLg showed moderate to high association with yield and was strongly related under DRT. The association of WSC with yield was significant but showed some direction inconsistencies across environments.

Fig. 1
figure 1

Trait association for yield across all six trials of the Seri/Babax population grown between 2002 and 2006. The diagonal contains the yield histogram for each trial, the lower diagonal a scatter plot and loess smoothing line between all trials, and the upper diagonal contains the phenotypic correlations

Table 3 Phenotypic (P) and genotypic (G) correlations of yield with agronomic, phenological, and physiological traits, and heritability across all environments (h 2G ALL ) and across heat and drought environments (h 2G HD )

For each trait, its heritability calculated across all environments (h 2G ALL ) or across stress environments (h 2G DH ) did not change substantially (Table 3). The most heritable traits were anthesis and maturity (>0.9). Comparing only the agronomic traits, TGW had the highest heritability, grain number (GM−2), and KWI reported heritabilities >0.5. From the physiological traits, WSC and NDVIv reported the lowest heritabilities (Table 3), and all the rest showed moderate to high values.

QTL mapping results for agronomic, phenological, and physiological traits

In total, 115 trait × environment × year combinations were analyzed for QTL. Of these, 109 combinations were associated with at least one significant QTL (see Supplementary Table 1 for the complete table of results). The QTL related to agronomic, phenological, and physiological traits averaged LOD scores of 2.92, 3.5, and 3.56, respectively. The average LOD for all significant QTL (LOD > 2) was 3.6 with maximum of 15 (see Supplementary Fig. 2 for LOD score distribution). A summary of the results is presented as a matrix of linkage groups by traits in Table 4; shading indicates repeatable (across years), stress, and robust (across environments) QTL (different font formats differentiate between these three groups); unshaded cells indicate QTL that did not appear in both years of any single environment. Table 5 shows the markers (and flanking markers) related to the largest effects in each repeatable, robust or stress QTL, and also other traits with QTL at the same marker. Except for D02, the greatest phenotypic variances (R 2) in all the environments × year combinations were associated with QTL of traits other than yield; maximum %R 2 values explained by markers associated with repeatable, robust, and stress QTL are presented in Supplementary Table 2. Not all markers of each QTL are listed, the criteria being to include only consistent QTL with markers showing an R 2 > 10% for at least one trait × environment × year combination.

Table 4 Summary of significant QTL found for all traits across environments showing the marker related to the maximum phenotypic variance (% R 2) for each trait and its corresponding linkage group and position
Table 5 Markers (and flanking markers) associated with the highest additive effect across all trials where a consistent QTL was detected for each trait

The 109 combinations were related to 104 significant QTL (across the 28 linkage groups), of which 14, 16, and 10 QTL were associated exclusively with DRT, HOT, and IRR environments, respectively, and 94 were associated with stress although not necessarily exclusively (Table 4). From the 104 QTL detected, 7 were repeatable, 17 were stress QTL, and 14 robust QTL. In addition, 9 linkage groups appear to contain two QTL for the same trait (Table 4; Supplementary Fig. 1). Including suggestive QTL, 1B-a contained the most QTL (10) related to yield, GM−2, CT, NDVI, and CHLg. All traits were related to multiple genomic regions with six linkage groups 1B-a, 2B-a, 3B-b, 4A-a, 4A-b, and 5A-a (Table 4; Supplementary Fig. 1) being identified as being most significant in controlling the traits studied based on the number of QTL signals detected. Grain yield was associated with QTL on eight different linkage groups; all of them were identified under DRT while 6 and 7 were detected under HOT and IRR environments, respectively (Supplementary Fig. 1). Six genomic regions were related mainly to yield, CT, and GM−2 with positive contributions mainly from Babax alleles (linkage groups 1B-a, 3B-b, and 4A-a), but some effects also from Seri on linkage groups 2B-a and 6B-a (Fig. 2). All QTL are represented in Supplementary Fig. 1; however, Fig. 2 shows that yield, GM−2, and CT were associated with QTL on at least five of the same linkage groups: 1B-a, 2B-a, 3B-b, 4A-a, 4B-b, and 6B-a.

Fig. 2
figure 2

Main genomic regions associated with yield, grain number (GM−2), and canopy temperature (CT) under drought (D), heat (H), and irrigated (I) environments. Blue and red colors are used to distinguish between Babax (red) and Seri (blue) LOD scores. On both sides of the figure is indicated the corresponding trial × year, canopy temperature (CT) is additionally labeled with the date of measurement given in number of days after emergence

QTL related to agronomic traits

Of the eight QTL related to yield, four were robust (i.e. across stressed and irrigated environments) appearing in linkage groups 1B-a, 3B-b, 4A-a, and 4A-b (Tables 4, 5; Supplementary Table 1). However, considering that distances greater than 30 cM between markers with LOD > 2 may indicate different QTL it is possible that linkage groups 3B-b and 4A-b contained more than one QTL for yield (Table 4; Supplementary Table 1). No repeatable yield QTL were detected as unique for any given environment. The QTL for yield in 1B-a appeared in all environments and were associated with markers mainly located in the region of the rye translocation (see below). The additive effects of all significant QTL expressed under DRT explained the largest proportion of variation in yield when compared with the other environments (Fig. 3). The strongest effect for yield (47 g m−2) of any QTL was found in linkage group 4A-a under DRT, and the largest favorable effects on yield were contributed by the Babax allele.

Fig. 3
figure 3

QTL effects for agronomic traits in the six trials. Positive (+) and negative (−) values are used to distinguish between additive effects of Babax and Seri alleles in each linkage group. Only effects where LOD > 2 are shown

On linkage group 1B-a, there was co-location of yield QTL with QTL for GM−2, NDVIv, CTv, CTg, and WSC (Table 4). Except for WSC, all these traits were found to be significantly correlated with yield (Table 3). In the 1B-a region, yield increases as well as the increases in GM−2 and NDVIv were usually controlled by the Babax allele (Figs. 3, 5), but warmer (unfavorable) canopy temperatures were generally dominated by the Seri allele while higher WSC contents were associated with both parental alleles (Fig. 5). Repeatable, stress, and robust QTL for GM−2 that coincided with those for yield were found on 1B-a, 3B-b, and 4A-a (Table 4). The strongest effect for GM−2 (654 grains/m2) was located on 1B-a, with a %R 2 of 11.8%; other traits associated with this marker were yield, height, CTv, CTg, NDVIv, CHLg, and WSC (Table 5). This QTL was robust and associated with the presence of the Babax (non-rye) allele, being detected under DRT, IRR and both years of HOT (Supplementary Table 1). It was also associated with the location of the third largest QTL effect for yield. On the 5B-a linkage group, the GM−2 QTL detected was repeatable for both years of the HOT treatment. Favorable effects from Seri and from Babax alleles (Fig. 3) were found for GM−2, but the expression of QTL showed the strongest additive effects were dominated by the Babax allele under the HOT environments. The strongest QTL effects for TGW were found on linkage groups 3B-a, 4A-b, and 4B-b (Table 4; Supplementary Table 1); the greatest repeatable effects for TGW were found in IRR environment (Fig. 3) mainly associated with the Seri allele in the 4B-b region (Table 5). The largest effects for KWI were found in the HOT environment, and the trait was increased by Seri alleles on 2B-a and by Babax alleles on 1A-a, 3B-a, and 6B-a (Fig. 3).

QTL related to phenological traits and height

Diverse genomic regions were found to contain QTL related to the phenological traits, and some of them were common between anthesis and maturity, which is consistent with the high genetic correlation between these traits averaged across experiments (data not shown, R 2 = 0.79). Three main regions affected anthesis date: 1D-b, 4D-a, and 5A-a (Table 4). The QTL located in the first two regions were classified as robust, while the QTL on 5A-a was found to be specific for stress. Linkage group 1D-b showed the strongest additive effects under DRT and IRR (0.5–1.3 days) contributed by the Babax allele (Fig. 4); conversely, a similar size QTL effect on 4D-a was contributed by a Seri allele in HOT environments (Table 5). The greatest QTL effects (ca. 1 day) were found in the IRR environment (Fig. 4). QTL for maturity were co-located with those for anthesis on 1A-b, 1D-b, 4D-a, 5A-a, and 7B-a (Table 4). The QTL at 1A-b and 4D-a were detected as repeatable under DRT and HOT, respectively; the QTL on 1D-b was robust, and the QTL at 5A-a and 7B-a were specifically associated with stress environments. The greatest effects for maturity (1 day) were found on 4D-a and were contributed by the Seri allele (Table 5; Supplementary Table 1) and in the HOT environments (Fig. 4; Supplementary Table 1).

Fig. 4
figure 4

QTL effects for the phenological traits. Positive (+) and negative (−) values are used to distinguish between additive effects of Babax and Seri alleles in each linkage group. Only effects where LOD > 2 are shown

Most of the repeatable, robust, and stress QTL related to other traits were found in locations independent of the main regions controlling phenology. The exceptions were two: a robust QTL for CTg on 5A-a linkage group and one stress QTL for NDVIv on 7B-a, both co-located with QTL for maturity. Those regions that had the highest number of QTL identified (1B-a, 3B-b, and 4A-a) were not related to anthesis or maturity.

One robust QTL was detected for plant height at 4B-b (Table 4) with the Seri allele (Fig. 4) increasing plant height by 1.6 cm (Table 5). The largest effect for this trait (4.4 cm) was a suggestive QTL detected under DRT on linkage group 4A-b with the increased height alleles again from Seri (Fig. 4).

QTL related to physiological traits

The QTL for CTv were identified on linkage groups 1B-a, 2B-a, 3B-b, and 7A-a as exclusively related to stress environments (Table 4) and on 4A-a as robust across environments. Irrespective of the time of the measurement (am/pm), the greatest effects for CTv were found in the HOT environment (Fig. 5). Unfavorable effects (high CT) were dominated by the Seri allele, and the largest increases in CT (0.46°C) were associated with the QTL on 4A-a; other negative effects of Seri existed on 3B-b and 7A-a. Unfavorable effects from the Babax allele of similar magnitude were also detected, i.e. at 1B-a (0.39°C) and 2B-b, but they were less frequent than the Seri effects. Four QTL for CTg (three stress and one robust) were co-located with four stress QTL for CTv (Table 4). The QTL identified at 3B-b and 5A-a were robust, while at 1B-a, 2B-a, and 7A-a, the QTL were specific to stress environments (Table 4). Two additional QTL were detected on 5A-a (robust) and UA-b (repeatable under DRT). Effects from alleles of both parents were found to be associated with CTg in both the morning and afternoon (Fig. 5), but the Seri allele was typically associated with the unfavorable increasing effect; the highest effects were detected in the DRT environment across diverse linkage group, i.e. increases of 1.5°C by a QTL on 4A-a (Fig. 5). NDVIv-related QTL showed the most favorable increasing effects (0.027) in the 1B-a region, controlled by the Babax alleles (Fig. 5; Table 5). This trait showed one robust QTL in 1B-a and one repeatable QTL under HOT on 4A-b (Table 4); the HOT environment showed the strongest and most frequent significant effects for the NDVIv (Fig. 5). During grainfilling, NDVIg was associated with specific environments, and the QTL effects varied according to the environment and linkage groups. One repeatable QTL under DRT was detected at 2B-a and two stress QTL on 4A-a and 7B-a (Table 4). The highest effects were found in the IRR environment (Fig. 5) contributed by alleles of both parents. One stress QTL for CHLg was detected on 1B-a, with positive effects associated with the Seri allele (Table 4; Fig. 5); this QTL at 1B-a was co-located with a suggestive QTL associated with NDVIg; other coincident QTL were found on linkage groups 1D-b and 5A-a. The suggestive QTL found on 1D-b under DRT and IRR showed the highest effect (0.78) for the CHLg (Supplementary Table 1). Only suggestive QTL related to WSC were detected (Table 4); the most and the largest WSC QTL were identified under the DRT environment and were mainly contributed by the Babax allele (Fig. 5).

Fig. 5
figure 5

QTL effects for the physiological traits. Positive (+) and negative (−) values are used to distinguish between additive effects of Babax and Seri alleles in each linkage group. Only effects where LOD > 2 are shown

Common QTL associated with drought and heat adaptive traits

Co-locations of DRT and HOT QTL were frequently detected in this study, pointing to the potential existence of common adaptation mechanisms useful under both stress conditions. Common QTL were found for all traits in these two stress environments (17 stress QTL), with eight QTL co-located with QTL in IRR environments. The QTL exclusively related to both DRT and HOT environments (Table 4) were specific by trait and classified as stress QTL, as previously indicated. These stress QTL were found on linkage groups 3B-b, 4A-a, and 6B-a for GM−2; on 4A-b for TGW; on 5A-a for anthesis; on 5A-a and 7B-a for maturity; on 1B-a, 2B-a, 3B-b, and 7A-a for CT; on 4A-a and 7B-a for NDVIg; and on 1B-a for CHLg (Table 4; Supplementary Table 1). A number of regions contained QTL in at least 1 year of both DRT and HOT environments, namely linkage groups 4B-b for yield; 3B-a GM−2; 1A-a for KWI; 3B-b and 4A-b for height; 4A-a for CTg; 3B-b NDVIv; 4B-b, 6A-a for NDVIg; 1D-a, 3B-b for CHLg and 1A-b for WSC. An exploratory analysis for stress indexes (calculated as IRR-Stress/IRR) was conducted in order to detect QTL directly related to stress performance traits. The regions associated with the stress indexes confirmed the previously identified QTL (data no shown), but no novel locations were identified.

Genetic dissection of yield into agronomic and physiological components

One approach for explaining yield through the genetic control of component traits is to consider genomic regions/QTL associated with yield components starting with TGW and GM−2 (Fig. 6). A more complete arithmetic analysis would include all yield components, namely number of plants per unit area, number of spikes per plant, number of spikelets per spike, number of fertile florets, and grain weight. In general, grain yield can be assessed through the combination of two components: grain weight and grain number. A multiple linear regression model across environments showed that 66% of the yield variance was explained by TGW and GM−2 (data not shown), and QTL for diverse traits have been found to be co-located with yield QTL (Fig. 2). The genetic dissection of yield QTL with respect to yield components is expressed herein as a Venn diagram (Fig. 6). Of the four robust QTL for yield, three (1B-a, 3B-b, and 4A-a) were co-located with QTL for GM−2 and one (4A-b) with TGW (Fig. 6a). Similarly, with respect to physiological traits, the yield QTL were also co-located with QTL for the following physiological traits: CT (three QTL), NDVI (three QTL), and CHL (one QTL), while suggestive QTL for WSC were common to all four QTL (Fig. 6b).

Fig. 6
figure 6

Dissection of yield QTL into agronomy- and physiology-related QTL. a QTL for yield (single lined curve) dissected into QTL related to agronomic traits, QTL for GM−2 (dotted circle) QTL for TGW (double lined circle); b QTL for yield dissected into physiological traits, QTL for CT (double lined circle); QTL for NDVI (dotted circle); QTL for CHL(single lined curve); Big fonts are robust QTL for yield while small fonts are suggestive QTL for yield

Discussion

QTL mapping for agronomic and physiological performance traits

Although several physiological traits have been shown in wheat to be associated with genetic variation or improvement under drought (Olivares-Villegas et al. 2007; Richards 2006), relatively few QTL for physiological traits have been identified to date (Rebetzke et al. 2008a, b; Reynolds and Tuberosa 2008), fewer still have been utilized in breeding, and none have been cloned (Collins et al. 2008). Table 6 summarizes the results from 23 studies that have identified chromosomal regions in wheat associated with the expression of agronomic, and in a few cases, physiological traits. Results of the current study show broad agreement with these previous studies (shaded cells in Table 6) with the main difference being fewer co-locations with phenology-related QTL. While the nature of QTL mapping restricts the comparison of specific markers across populations, chromosomal regions are still indicative. For the Seri/Babax population all traits were controlled by multiple genomic regions (Figs. 3, 4, 5), especially under stress, highlighting the complex nature of the adaptation process. In this study, the most important chromosomal regions were 1B-a, 2B-a, 3B-b, 4A-a, 4A-b, and 5A-a (Table 4; Supplementary Fig. 1) which were associated with 10, 9, 8, 7, 7, and 7 traits, respectively.

Table 6 Summary of recent research in QTL mapping for wheat performance-associated traits and their relevance to the current study

In other populations, QTL on the short arm of chromosome 1B have been related to performance traits such as those related to transpiration efficiency (Rebetzke et al. 2008a) and to yield (Marza et al. 2006). The increasing yield effects in this study associated with the short arm of chromosome 1B were from the Babax parent, i.e. the allele from Seri associated with the T1BL.1RS (rye) translocation resulted in decreased yield. This effect was also reported by Mathews et al. (2008) for Australian drought and irrigated trials of the Seri/Babax population. In previous studies of isolines of Seri M82, the T1BL.1RS rye translocation was shown to increase yield by 11% under terminal stress (Villareal et al. 1998). The 1BL.1BS donor for the Villareal study was an older CIMMYT line, Pavon F76, which may well have had weaker alleles for yield on 1BS compared to those present in the Babax parent in this study. In common with Mathews et al. (2008), Peake (2003), examining a cross between Seri M82 and Hartog (a sister line of Pavon F76), also found no advantage of the T1BL.1RS rye translocation in Australian environments. The major use of the T1BL.1RS rye translocation has been to improve disease resistance although it has also frequently been reported to decrease grain quality, e.g. Gobaa et al. (2008). In conjunction with other reports such as Rebetzke et al. (2008a), it seems that there is alleleic variation in the effects of the 1BS chromosome segment, and that the impact on yield is not always inferior to that of 1RS, at least in environments where the disease pressures are low.

The co-location of two or more QTL is indicative of pleiotropic or linkage effects (Huang et al. 2004), for example where yield QTL co-located with other QTL for CTg, CTv and GM−2 on 3B-b. The 3B-b region was previously reported by Miura and Worland (1994) to affect vernalization, but in these relatively warm (spring wheat) environments it was not associated with any variation in flowering time (Table 1). The yield and CTv increments controlled on region 4A-a by contrasting alleles confirms their strong negative correlations (Olivares-Villegas et al. 2007) and suggests that cooler canopies, due to the absence of the Seri allele on the 4A-a region, would result in higher yields, especially under irrigated conditions exposed to continuous high aerial temperature. The robust QTL found on 4A-a (identified as the yield QTL of highest effect in Table 5) was previously reported by Kuchel et al. (2007) and Marza et al. (2006) as being related to yield, but those pedigrees were quite different to the Seri-Babax population (coefficient of parentage with Seri M 82 and Babax <0.077). This region has been identified as affecting a number of important traits, including yield, flowering time, amylase content, and grain weight (Araki et al. 1999; Kirigwi et al. 2007; McCartney et al. 2005). Coincident location of diverse QTL for yield, yield components, and some physiological traits in the 1B-a, 3B-b, 4A-a, and 4A-b (Table 4) suggest the possible existence of a group of genes involved in important mechanisms for performance of wheat under DRT (Kirigwi et al. 2007), as well as HOT, and IRR environments.

The repeatable QTL for GM−2 on 5B-a (Table 4) was also reported as co-locating with a yield QTL by Kuchel et al. (2007). Two robust QTL for TGW and height on 4B-b (Table 4) have been previously reported by McCartney et al. (2005) and Huang et al. (2004) in a similar region. These authors had suggested that this region may be associated with the Rht-B1 gene, and later Marza et al. (2006) confirmed that the 4B QTL had a strong influence on height in his population. However, Seri and Babax do not segregate for any known Rht genes, and given their common pedigree likely share the same Rht allele regardless.

Value of using progeny expressing a restricted range of phenology

Variation in wheat phenology is determined by genes of major effect (Ppd and Vrn), as well as minor effects associated with “Earliness per se” (Eps) (Worland 1996). Given that alleles for major Ppd and Vrn are known, their effects on drought adaptation can be studied relatively easily using near isogenic lines generated by backcrossing into any genetic backgrounds (Dyck et al. 2004). When studying minor genetic effects, a population should comprise individuals with synchronized phenological cycles, such that all lines will pass through the same developmental stages at the same time, and hence, assure that all the genotypes will experience the same environmental conditions in each stage. By avoiding segregation for genes of major phenology effects in mapping populations, the probability of identifying genes of minor effect (i.e. those affecting complex traits) is expected to increase (Reynolds and Tuberosa 2008) as recently demonstrated by subdividing a large DH wheat mapping population (Kukri/RAC875) into early and late sub populations for QTL analysis (Reynolds et al. 2009). The relatively narrow range of phenology of the elite Seri/Babax population helps avoid confounding of environment patterns (Olivares-Villegas et al. 2007) as indicated by the following observations.

  • In contrast to previously reported studies (Table 6), there were minimum coincidence of any repeatable, stress or robust QTL for yield, yield components, CHL, NDVI or WSC with phenology-related QTL (Table 4); many other studies found common regions for phenology and yield (Araki et al. 1999; Kirigwi et al. 2007; Kuchel et al. 2007; Marza et al. 2006; McCartney et al. 2005), CID (Rebetzke et al. 2008a), TGW (Börner et al. 2002; Huang et al. 2004), and GM−2 (Shah et al. 1999).

  • In the current study, the HOT environments with the lowest variation in anthesis date also had more significant QTL effects for all physiological traits (yield, NDVI, WSC, CT, and CHL) compared with IRR environments which had greater variation in flowering time (Figs. 3, 5).

  • Herein, a group of QTL “clusters” were clearly identified for physiological traits, i.e. CT, CHL, WSC, and NDVI, indicating the existence of consistent genomic regions that contain important groups of genes controlling or affecting their expression. The QTL appeared repeatedly across stressed and non-stressed environments at 1B-a, 2B-a, 3B-b, 4A-a, and 7A-a linkage groups.

The fact that most reported studies in wheat mapping populations do not even report phenological range, and in all cases where it is the range is relatively large, up to 3 weeks (Table 4), confirms that control of phenology is generally an experimental bottleneck to gene discovery for complex traits, such as those associated with drought and heat adaptation.

Common QTL associated with drought and heat adaptive traits

Different conditions induce gene expression at different loci, but some genomic regions may have value under multiple environments. QTL that are significant under different kinds of abiotic stress, such as those identified in common to both DRT and HOT environments in the current study (Fig. 2), could be useful in breeding for target environments that encompass a range of abiotic stresses. The overlapping of intervals for CTv and CTg further suggested that the ability to maintain low CT may be associated with the same genes at different growth stages as well as under different environments (Fig. 2); favorable expression of CT under drought is linked to increased root depth in drought adapted Seri/Babax sisters (Lopes and Reynolds 2010) while improved root capacity is also likely to permit a larger transpiration rate and therefore cooler CT under hot conditions. Other regions have suggested common drought and heat associated genes for anthesis (4D-a, 5A-a), maturity (5A-a, 7B-a), KWI (1A-a, 2B-a), and TGW (4A-a) (Table 4).

Genetic dissection of yield

Grain yield is an integrative trait including the net effect of all physiological processes during the crop cycle. However, significant interactions of yield and QTL effects with environment, especially under stress, make it difficult to use yield QTL directly in selection. The dissection of yield into associated traits could assist breeders to identify specific desirable characteristics with higher heritability and tagged by DNA markers for more reliable selection; cf. selection for yield per se. As shown in Fig. 6, all of the repeatable, robust, and stress QTL for yield were common to GM−2 or TGW QTL, while considering also the suggestive QTL, more than 50% of the QTL detected for yield were common to the yield components. Common genomic regions for yield and yield components have also been reported in previous studies, i.e. co-locations of QTL related to yield and QTL for grain weight (Araki et al. 1999; Kato et al. 2000; Marza et al. 2006; McCartney et al. 2005) and grain number (Campbell et al. 2003; Kirigwi et al. 2007). The co-location of QTL associated with different traits can be the result of: (1) Two strongly linked genes affecting different traits, (2) One single gene that produces a series of effects in related traits, (3) One gene affecting two or more independent traits, (4) Two linked genes with effects in the same traits (Yang et al. 2007). Groos et al. (2003) have suggested that QTL for traits like TGW could be successfully used for yield improvement due to its accurate detection and repeatability across environments in comparison to the QTL for yield, and it could be complemented by adding correlated grain number QTL. However, in practice, this could be difficult to achieve to the extent that grain number and potential grain size could be affected by competing mechanisms during early stages of ovule and grain development. In this study, GM−2 and TGW showed five and three repeatable, robust or stress QTL, respectively (which could therefore be selected for), in comparison to yield which showed four QTL.

Physiological dissection showed that all of the four yield QTL were accounted by CT, NDVI, CHL, and WSC together, the first three being traits that may be more closely associated with biomass production per se (via maintenance of transpiration and leaf nitrogen content) rather than grain set processes. In contrast to yield QTL, selection for the respective QTL associated with physiological traits can be guided by conceptual models of trait expression, in accordance with the specific characteristics of the target environments (Reynolds et al. 2007). For example, QTL for cooler canopies associated with access to water at deeper soil profiles, (Lopes and Reynolds 2010) would only be a priority for target environments with suitable soil characteristics and QTL associated with high ground cover (NDVI) may not be needed in farming systems where soil cover is achieved through crop residue retention and/or crops that rely on stored soil water.

Conclusions

This population did not segregate for major flowering time genes, with a consequent restriction of phenology. Unlike most other similar studies, QTL for anthesis were generally not co-located with performance QTL. Furthermore, QTL were identified that explained up to 27% variation in yield and 28% variation in canopy temperature, independently of confounding effects of phenology.

Common genomic regions for drought and heat suggest a generic value across stresses for several of the QTL identified. Chlorophyll (CHLg), NDVI, and canopy temperature (CT) are examples of traits previously reported as being related to yield under both drought and heat. According to the results presented here, their genetic basis may be associated with the same genomic regions and controlled by the same parent allele irrespective of the environment.

The co-location of QTL for diverse agronomic and physiological traits with QTL for yield directly support the genetic dissection of the crop performance in order to facilitate a more strategic approach to breeding for adaptation. Those regions identified across environments are candidates that can be used in MAS or gene cloning, especially if they show moderate to high broad sense heritabilities.