Application of a library of near isogenic lines to understand context dependent expression of QTL for grain yield and adaptive traits in bread wheat
BMC Plant Biology volume 16, Article number: 161 (2016)
Previous quantitative trait loci (QTLs) studies using the Avalon × Cadenza doubled haploid (DH) population identified eleven QTLs determining plant height, heading date and grain yield. The objectives of this study were: (i) to provide insight into the effects of these QTLs using reciprocal multiple near isogenic lines (NILs) with each pair of alleles compared in both parental backgrounds (Avalon or Cadenza), (ii) quantifying epistasis by looking at the background effects and (iii) predict favourable allelic combinations to develop superior genotypes adapted to a target environment.
To this aim, a library of 553 BC2 NILs and their recurrent parents were tested over two growing seasons (2012/2013 and 2013/2014). The results obtained in the present study validated the plant height, heading date and grain yield QTLs previously identified. Epistatic interactions were detected for the 6B QTL for plant height and heading date, 3A QTL for heading date and grain yield and 2A QTL for grain yield.
The marker assisted backcrossing strategy used provided an efficient method of resolving QTL for key agronomic traits in wheat as Mendelian factors determining possible epistatic interactions. The study shows that these QTLs are amenable to marker assisted selection, fine mapping, future positional cloning, and physiological trait dissection.
Many of the key traits for wheat such as yield, flowering time, and canopy architecture are quantitatively inherited, displaying continuous variation due to the contribution of multiple genes, each with relatively small effects, called quantitative trait loci (QTL). Greater understanding of the physiological and molecular basis of individual QTL is a major goal in the identification of fundamental biological processes controlling them and guiding improved manipulation of higher level quantitative traits. One way to achieve this is to convert quantitative to qualitative traits through the production of Near Isogenic Lines (NILs). For bread wheat QTLs this has been achieved in just a few instances eg [1, 2]. This study represents a process of systematic QTL validation by the production of NILs in a multi stream backcrossing strategy that not only serves as a validation platform for the majority of QTLs identified in Avalon × Cadenza UK reference segregating population, but also as a unique tool for the further understanding of context dependent QTL expression.
The major target in wheat breeding is grain yield (GRYLD). Increasing GRYLD is challenging for breeders due to its complexity, environmental and genetic interactions and low heritability. To reduce the complexity of GRYLD it can be broken down into numerical components: grain number (GRpsqm) and thousand grain weight (TGRWT) which are often negatively associated [3, 4]. GRpsqm is the strongest historical determinant of GRYLD increases (usually due to the number of grains per spike) [4–6] and currently the yield component used to improve GRYLD potential in wheat breeding . GRpsqm is also one of the most plastic  and least heritable  of the yield components, so precise genetic stocks, which allow the precise quantification of genetic and environmental interactions are essential for deeper genetic analysis of this trait. Potential improvement in GRYLD driven by GRpsqm can only be fully realised if increases in GRpsqm are complemented with improvements in TGRWT, while avoiding compensatory effects (trade-off) between both yield components [10, 11]. TGRWT can be further dissected into individual components including grain length (GRL), width and area. GRL was found to be positively associated with TGRWT . The identification of QTLs in the same position for TGRWT and GRL support the idea that GRL is a key trait for TGRWT and it is possible to achieve higher yields with large grains . Several studies have identified QTLs for GRYLD and grain yield–related traits [12–28].
Lying beneath the numerical components of grain yield are a complex web of developmental traits, the genetic control of which has just as an important influence on the final expression of GRYLD, as well as crop adaptation, and must be considered in wheat breeding. These include crop height (PH) and heading date (HD) which are more stably inherited. Short stature and optimal flowering are breeding pre-requisites for achieving high yield varieties. Reduced PH has been exploited by the introduction of Rht dwarfing genes (Rht-B1b, Rht-D1b) which reduce lodging, increase harvest index and yield potential [29, 30]. Rht-B1 (on 4BS) and Rht-D1 (on 4DS) reduce PH by reducing the response to gibberellin with pleiotropic effects on grain number and yield . The increase in grain yield and grain number was driven by the increase in grains per spike [29, 30, 32]. Another dwarfing gene, Rht8, on 2DS has also been used to reduce stature and improve lodging resistance, without yield penalty [33, 34]. Gasperini et al.  determined that the reduced stature is a consequence of shorter internodes with the largest differences obtained in the peduncle and the first internode. In addition to the major crop height genes, other QTLs for PH have been identified in a number of studies .
Wheat is the most widely grown grain crop in the world because of its broad geographical adaptation. Part of this adaptation is due to the manipulation of genes that control HD (Photoperiod (Ppd), Vernalization (Vrn) and ‘Earliness per se’ (Eps) genes). Ppd genes control day length response, Vrn genes determine the requirement for a period of cold, Eps genes are important for the fine-tuning of flowering time and account for the variation in HD after photoperiod and vernalization requirements are fully satisfied. QTLs for HD and PH have been found in the same position, with early alleles associated with increased height in the Avalon × Cadenza DH population [36, 37]; this could therefore, be a trade-off regarding yield manipulation. QTLs for HD in wheat have been identified in several studies [13, 38–43].
Using the Avalon × Cadenza DH population as a platform, Ma et al.  suggested that the ideal UK target genotype should be: PH and HD similar to Avalon; Rht-D1b and Vrn-A1b alleles; high GRpsqm and TGRWT; long and wide grains and maximum GRYLD. Breeding for superior genotypes can be achieved through the combination of favorable alleles using marker-assisted selection (MAS). To this aim, NILs have to be developed to validate the QTLs identified in the Avalon × Cadenza DH population for PH , HD  and GRYLD . The objectives of this study were: (i) to validate these QTLs using the NILs developed by MAS using two recurrent parents (Avalon and Cadenza), (ii) determine possible epistatic interactions and (iii) define the expected genotype if the target regions were introgressed in the same background in order to to develop superior genotypes adapted to a target environment.
Plant material and genotyping
The Avalon × Cadenza DH population was developed as a tool for the genetic discovering of wheat, as part of the Wheat Genetic Improvement Network (WGIN) (www.wgin.org.uk) and represents a broad spectrum of elite UK winter germplasm produced in different UK wheat breeding programmes. The parents were chosen for contrasting canopy architecture traits. Avalon and Cadenza carry recessive photoperiod sensitive alleles of Ppd-D1 and Ppd-B1. Avalon (UK winter wheat) carries recessive alleles of Vrn-A1, Vrn-B1 and Vrn-D1 whereas Cadenza (UK alternative wheat) carries the dominant Vrn-A1a allele. The parents also differ for Rht-D1, Avalon carries the dwarf allele Rht-D1b whereas Cadenza carries the wild type allele (Rht-D1a).
QTLs controlling flowering time, plant height and grain size and shape had been previously identified in the Avalon × Cadenza DH population using meta-QTL analysis [15, 36, 38]. Eleven QTLs, on chromosomes 1B, 1D, 2A, 2D, 3A, 3B, 5A, 6A, 6B, 7B and 7D were chosen as target regions for introgression in our marker assisted backcrossing scheme (Fig. 1). A total of 31 microsatellite (SSR) markers linked to the target QTLs were used for foreground selection (Table 1).
Donor Avalon × Cadenza DH lines were crossed with recurrent parents (either Avalon or Cadenza) to produce F1 plants. These plants were grown and backcrossed twice with the recurrent parent to produce two independent BC2 streams. After each backcross, progenies were genotyped at the QTL of interest, and those heterozygous at the QTL were carried forward to the next cross. BC2 individuals were then self fertilised and homozygous (carrying Avalon or Cadenza alleles in the QTL region) were selected with flanking SSR markers, self fertilised and provided seeds for multiplication and then large scale experiments. The BC2 derived NILs contain an estimated recurrent parent genomic background of at least 93.75 %. A total of 553 BC2 NILs were generated (250 and 303 NILs with Avalon and Cadenza background, respectively; Table 2). More details of the development of the Avalon × Cadenza NIL population and streams pedigree are available in WGIN website (www.wgin.org.uk).
Experimental design and phenotyping
Field trials were conducted at Church Farm, Bawburgh, Norfolk (UK) in two growing seasons (2012/2013 and 2013/2014). Sowing date differed between trials, the 2012–2013 trial was sown the first week of November 2012 whereas the second trial was sown the second week of October 2013. Church Farm has a sandy clay loam soil type. Fungicides and fertiliser products were applied according to standard agronomic practice. The experimental design for both trials consisted of a randomized complete block design with three replicates (some NILs had more replicates) and the two parents (Avalon and Cadenza). The experiment contained a total of 1800 and 1868 plots in 2012/2013 and 2013/2014, respectively. The additional plots in 2013/2014 were due to an increase number of replicates for the parents. Each plot consisted of 7 rows 20 cm apart and 4 m long. A sowing rate of 260 seeds.m−2 was used in all cases. An additional file shows the details of experimental conditions in more detail (see Additional file 1). Heading date was scored as the day when 50 % of the ears in a plot were more than 50 % emerged from the flag leaf and was assessed in thermal time (HD) (°C days, using a base temperature of 0 °C). Plant Height (PH) was measured from soil level to the top of each wheat ear. Grain Yield (GRYLD) was recorded per plot in both years. Prior to harvest, 10 main shoots were sampled from each plot for assessments of yield components including thousand grain weight (TGRWT), grain width (GRW), grain length (GRL) and grain area (GRA) using a Marvin digital seed analyser (GTA Sensorik GmbH). Grain number (GRpsqm) was estimated from GRYLD and TGRWT. PH components (ear length (EL), peduncle (PL) and internode lengths (from the first to the fifth counted from the top, abbreviated as 1stITL, 2ndITL, 3rdITL, 4thITL, 5thITL)) and GRpsqm components (spikes/m2 (S), spikelet/spike (s/S) and grains/spikelet (G/S)) were only measured in 2013. Harvest Index (HI) and Biomass was only measured on 3A NILs. This was done by selecting an internal row from the plot and selecting a 60 cm section of this row that was judged to be representative of the whole plot. Plants within this strip were then cut to ground level and threshed so that the ration of grain yield to total above ground biomass could be calculated.
In order to remove spatial effects and to estimate heritabilities a mixed model was fitted for all traits. The model used included environment, replicates within environment and checks within environment as fixed factors, and row (of plots) within environment, column within environment, genotype (NIL lines) and the interaction between genotype x environment as the random factors. Narrow-sense heritabilities were estimated using the variance components from this model, as follows: h2 = σ2 g / [σ2 g + (σ2 ge/e) + (σ2 res/re)]; where σ2 g is the variance component of genotype; σ2 ge is the variance component of the interaction genotype x environment; σ2 res the error; e, the number of environments, and r, the number of replicates per trial.
Best Linear Unbiased Estimators (BLUEs) were estimated for each genotype within individual trials considering genotype as fixed, and replicate, row and column as random factors. The estimated BLUEs were used to determine significant difference between Avalon and Cadenza alleles in the target QTL regions using linear mixed models. The first model included as fixed factors the environment, background, allele nested within QTL and the 2, 3 and 4-level interactions, while the random factor was the remaining genotypic variance. All analysis were conducted using Genstat 16th .
Eleven QTLs for three agronomic traits (HD, PH and GRYLD) in the two genetic backgrounds of the bi-parental Avalon × Cadenza DH population, which were used for their discovery, were investigated using 553 BC2 NILs and their recurrent parents in the field for two years. The mean values of the two parents (Avalon and Cadenza) and the NILs in each environment are shown in Table 3. In all cases, Cadenza was taller and flowered earlier than Avalon (with a difference of 3.14 cm and 14 °C days to heading in 2013 and 6.12 cm and 24.90° days to heading in 2014). In 2013 Avalon and Cadenza had similar grain yield (8.27 T/Ha and 8.11 T/Ha, respectively). However, in 2014, Avalon had lower grain yield (10.45 T/Ha) than Cadenza (11.56 T/Ha). These differences can be explained dissecting GRYLD into its components (GRpsqm and TGRWT). Avalon grown in 2013 had greater GRpsqm and TGRWT compared to 2014. The mean values of GRA and grain shape (L/W) for the parents were consistent in both years; Avalon had lower GRA and higher L/W ratio than Cadenza. There were no significant differences between the parental lines for most of the traits with the exception of PH in 2013. However, in 2014, all the traits except GRA (P = 0.301) and GRpsqm (P = 0.004) showed significant differences at P < 0.001 (Table 3).
Means of the NILs for GRYLD values ranged from 5.01 to 9.70 and 8.20 to 14.13 T/Ha in 2013 and 2014, respectively (Table 3; Fig. 2). For PH, the values ranged from 58.65 to 82.35 cm in 2013 and 70.11 to 106.26 cm in 2014. For HD, the values ranged from 1287 to 1430 and 1717 to 1848° days to heading in 2013 and 2014, respectively. This large genetic variability reflected transgressive segregation since differences between parents were smaller than the range of NILs lines indicating the polygenic inheritance of the trait alleles. Distributions of NILs lines were approximately normal in all cases with the exception of GRA and GRW in 2013. NILs can be classified in two backgrounds depending on the recurrent parent (Avalon or Cadenza) used in the marker assisted backcrossing scheme. Significant differences of the mean of NILs with the Avalon versus Cadenza background were detected for all the traits in 2014; however, only HD and PH, were significant in 2013 (Table 3). NILs with the Avalon background were shorter, flowered later and had lower grain yield than NILs with the Cadenza background. Overall, the mean values of the NILs were similar compared to the recurrent parent used in each case for all the traits in both years with the exception of GRYLD, GRpsqm and TGRWT in 2013.
These results showed a significant background × year interaction that can be partially explain by the difference between the two growing seasons (2012/2013 and 2013/2014) in terms of temperature and precipitation during the growing cycle (see Additional file 1). To meteorologically characterize each environment, we arbitrarily defined three consecutive ‘phenophases’: vegetative phase (from sowing until March), reproductive phase (from April until flowering) and grain filling phase (one month after flowering). Overall, the rainfall during 2014 (591 mm) was higher than in 2013 (471 mm) throughout the growing cycle. For each phenophase, the biggest difference was observed in the grain filling phase (14.7 and 30.1 mm in 2013 and 2014, respectively), followed by the reproductive phase (98.6 mm in 2013 and 153.4 mm in 2014) and finally the vegetative phase (357.8 mm in 2013 and 407.5 mm in 2014). In terms of temperature, the biggest difference was observed during the vegetative phase (604.15 and 1240.45° days in 2013 and 2014, respectively), followed by the grain filling phase (943.10° days in 2013 and 678.00° days in 2014) and finally the reproductive phase (520.30 and 415.85° days in 2013 and 2014, respectively). The duration of each phenophase is generally determined by interactions between genetical and environmental factors such as temperature and rainfall (these effects differ depending on when the factors act). Increasing temperatures, in general, accelerate phenological development resulting in shorter developmental phases and reduce GRYLD in the absence of terminal heat and drought . Another important factor to consider, which differed between trials, is the sowing date as this can also affect the duration of the phenophases. In our study, the 2012–2013 trial was sown in November 2012 whereas the second trial was sown in October 2013. The combination of early sowing with favourable meteorological conditions (wetter and warmer) resulted in taller plants, earlier flowering and increased GRYLD and yield components in 2014. Our data suggested that GRYLD is largely determined by GRpsqm which is established by anthesis and can be reduced as a consequence of water availability and/or high temperatures during the critical period (20 days before anthesis and 10 days late. The incresased GRpsqm observed in 2014 was therefore partially due to high rainfall and slightly lower temperatures during this phase. Lower night-time temperatures would also reduce rates of respiration making more photosynthetic products available for growth and development. TGRW is determined by the rate and duration of grain filling which is driven by genetic factors (cultivar) and environmental factors (temperature) [46–48]. The temperature during grain filling was slightly lower in 2014 but higher rainfall was observed. In addition, GRpsqm and TGRW can also be positively influenced by the duration of the vegetative phase which has a positive effect on GRYLD.
Correlations among traits were analysed for each environment and summarized in Table 4. Correlation coefficients were low to moderate among GRYLD, PH and HD. There was no significant correlation of HD, with GRYLD, or PH in 2013 (0.07 and −0.04, respectively). A statistically significant correlation was found between PH and GRYLD (0.48). However in 2014, the correlations between traits were 0.25, −0.16 and −0.37 for GRYLD – PH, GRYLD – HD and PH – HD, respectively, indicating that early flowering genotypes were associated with higher grain yield and taller plants. GRYLD was highly positively correlated with GRpsqm (0.80 in 2013 and 0.89 in 2014) and displayed lower positive correlation with TGRWT (0.31 and 0.22 in 2013 and 2014, respectively). As expected, a significant negative correlation was found between TGRWT and GRpsqm (−0.31 in 2013 and −0.24 in 2014). In a high yielding year, like 2014, TGRWT showed high significant correlation with GRA (0.87) and GRW (0.86) and was moderately correlated with GRL (0.38). Measures of grain size (GRA and TGRWT) were not independent of grain shape (L/W) as significant negative correlations were obtained between GRA and L/W (−0.79 in 2013 and −0.33 in 2014) and TGRWT and L/W (−0.36 and −0.54 in 2013 and 2014, respectively). Broad sense heritabilities ranged from 0.16 to 0.82 with PH and HD showing the highest heritability across years (Table 3).
In order to determine if there was a significant difference between Avalon and Cadenza alleles in the target QTL region, based on the background and year, a mixed model was carried out for all the traits. The results are showed in Tables 5 and 6.
For PH, the selected target regions for the development of NILs were on chromosomes 2A, 2D, 3A, 3B, 6A and 6B. Significant differences between Avalon and Cadenza alleles were found for the QTLs on chromosomes 2D, 3A, 3B, 6A and 7D (Tables 5 and 6). For the 3B and 6A QTLs, NILs carrying the Avalon allele were taller than NILs carrying the Cadenza allele. However, on 2D, 3A and 7D QTLs, NILs carrying the Avalon allele were shorter. In all cases, the results agree with previous QTL analysis [36, 37]. Compensatory effects were found between PH and GRYLD such as for the QTLs 3B and 6B in which the reduction of PH was not associated with GRYLD penalty.
The complexity of PH can be reduced by studying its components including EL, PL and internode lengths (1stITL, 2ndITL, 3rdITL, 4thITL and 5thITL) which were measured in 2013 only. For those QTLs that showed significant differences for PH (2D, 3A, 3B, 6A and 7D), the ratio between each PH component and the total PH was analysed. PL had the strongest influence for PH on 3B (1.27 cm), 6A (0.72 cm) and 7D (2.04 cm). On 2D, 1stITL showed the strongest effect (0.72 cm). On 3A, 3rdITL had the strongest influence on PH (1.18 cm). Additional file 2 supporting these results is included.
For HD, the chromosomes specifically selected for NIL development were 1B, 1D, however QTLs for HD were also captured on 3A, 6A and 6B as these were located in regions that also influenced PH and GRYLD. Significant differences were observed for 1B, 1D, 3A, 6A and 6B QTLs (Tables 5 and 6). NILs carrying the Avalon allele flowered later at the 1D QTL whereas those on 1B, 3A, 6A and 6B QTLs flowered earlier. These results are consistent with previous studies [37, 38]. The 1D QTL allow the manipulation of heading date for the purposes of better crop adaptation, with minimal ramifications for GRYLD.
Grain yield and yield components
For GRYLD, the chromosomes selected for the development of NILs were 2D, 3A, 3B, 5A, 7B and 7D. Significant differences between Avalon and Cadenza alleles were observed for the 1D, 3A and 7B QTLs (Tables 5 and 6). NILs carrying the Cadenza allele in the introgressed region on chromosomes 7B had reduced GRYLD. However, the Cadenza allele on the 1B and 3A QTLs conferred the increased GRYLD.
GRYLD was broken down into its components in order to elucidate which explained the significant differences detected for GRYLD and to identify instances in which variation in components is not translated into variation for GRYLD. For GRpsqm, significant differences were found for the 1D, 3A, 6A, 6B and 7B QTLs (Tables 5 and 6). On 1D and 7B QTLs, NILs carrying the Avalon allele had more grains/m2 whereas those on 3A, 6A and 6B QTLs had less grains/m2. GRpsqm is the result of the spikes/m2 (S), spikelet/spike (s/S) and grains/spikelet (G/S) which were only measured in 2013. For S, the Avalon allele significantly increased the spikes/m2 on the 7B QTL. For s/S, significant differences were detected on the 3A, 6A and 7B QTLs. NILs carrying the Cadenza allele had more spikelets/spike. For G/S, significant differences between Avalon and Cadenza alleles were observed for the 2A, 3A and 6A QTLs. These had effects in opposite directions; NILs carrying Cadenza allele on 3A and 6A had more grains per spikelet whereas those on 2A had less grains. An additional file supporting these results is included (see Additional file 3).
For TGRWT, significant differences were found for the 2A, 5A and 6A QTLs (Tables 5 and 6). The Avalon allele increased grain weight on 6A, whereas Cadenza provided the increased grain weight on 2A and 5A. TGRWT can be further dissected into individual components including GRL, GRW, GRA and L/W. Three QTLs for GRA (2A, 5A and 6A) and GRL (3A, 5A and 6A) and two for GRW (5A and 6A) and L/W ratio (3A and 6A) showed significant differences.
In our study, GRYLD was driven by GRpsqm as the QTLs that showed significant difference between alleles for GRYLD were also significant for GRpsqm. Only the 6A QTL showed significant TGRWT effects linked to GRpsqm. Avalon allele increased TGRWT through an increasing in GRA, GRL and GRW and reduced GRpsqm by reducing spikelet/spike (s/S) and grains/spikelet (G/S).
The mixed model fitted was used to determine possible epistasis by the analysis of the Background × Chromosome (Allele) interaction (B × C(A)) which indicates possible QTL × background or QTL × donor segment interactions (Tables 5, 6 and 7). HD (p-value = 0.004) and PH (p-value = 0.050) showed significant (B × C(A)) interaction. For HD, epistatic interactions were detected for the 3A and 6B QTLs. In both cases, the magnitude of the differences was greater in the Avalon background (34 °C and 21 °C on 3A and 6B, respectively) compared to the Cadenza background (16 °C on 3A and 9 °C on 6B). These QTLs also showed epistatic interactions for other traits. On 3A QTL, background effects were found for GRYLD whereas on 6B were detected for PH, an increase of GRYLD and PH was greater in the Cadenza background (3.42 cm on 6B and 0.88 T/Ha on 3A). On 3A, the greater GRYLD effects on the Cadenza background can be explained by the reduction of PH (7.31 and 4.65 cm in the Avalon background and the Cadenza background, respectively) and the increase of Harvest Index (HI) (0.01 in the Avalon background and 0.02 in the Cadenza background) through the decrease in Biomass (3.47 in the Avalon background and 3.10 in the Cadenza background) in 2013. In wheat breeding, HI and GRYLD have been increased by reducing PH providing crop advantage as not only the risk of lodging is reduced also the distribution of resources in the plant is better . For GRYLD, epistatic interactions were also detected on 2A QTL. The magnitude of the differences was greater in the Avalon background (0.53 T/Ha) compared to the Cadenza background (0.03 T/Ha).
Previous QTL studies using the Avalon × Cadenza DH population identified the location of QTL determining HD, PH and GRYLD [15, 36–38]. The development of NILs using marker-assisted selection helped to validate these QTLs and will define favourable allelic combinations to develop superior genotypes adapted to a target environment. In the present study, a set of NILs were developed using two recurrent parents (Avalon or Cadenza). This reciprocal approach allowed us to confirm not only the effect of the QTLs, but also determine if the expression of these QTLs were specific for the recurrent parent used in the backcrossing scheme or year.
The 553 BC2 NILs are divided in two groups depending on the background (Avalon and Cadenza). NILs carrying the Cadenza background were taller, flowered earlier and had increased GRYLD and yield components in both years. Part of these differences can be explained by Vrn-A1 (on 5A) and Rht-D1 (on 4D). Although the NILs were not genotyped for Vrn-A1 and Rht-D1 we expected them to be segregating as the recurrent parents differed for Vrn-A1 and Rht-D1. Avalon carries the recessive vrn-A1 and dwarf Rht-D1b alleles whereas Cadenza carries the dominant Vrn-A1a and the wild type Rht-D1a allele. Rht-D1 had no effect on flowering time . Although plants experienced warmer temperatures during the vegetative phase in 2014, there was still enough vernalizing cold as the difference for HD between NIL groups were similar for both years (27.62 and 26.28° days in 2013 and 2014, respectively). However, warmer temperatures and early sowing allowed NILs carrying the Cadenza background (spring types) to develop more quickly than NILs carrying the Avalon background (winter types), then slowed growth and development due to lower temperatures during winter without affecting plant survival. This could explain the inconsistency in the magnitude of the effects between NIL groups for PH (1.22 cm in 2013 and 4.24 cm in 2014) which can also be associated with higher GRYLD and yield components.
Comparison with previous studies
The present study validated the HD, PH and GRYLD QTLs identified by [15, 36–38]. However, some differences were observed with the direction of the QTL effect in the NILs when compared to those expected from the original A × C DH population. These differences are probably due to G × E interactions. For HD, the Avalon allele delayed flowering for the 7B QTL whereas in the original DH population the delay was detected for the Cadenza allele. For the 6A QTL in the NILs, the Cadenza allele increased the grain length whereas the Cadenza allele decreased it in the original population. In addition, this study allowed the identification of QTLs that showed significant differences which were not detected in the original DH population such as 2A for HD, 7D for PH, 1B and 7B for GRYLD.
An additional point to consider is the magnitude of the effects, so the effect in the ‘new’ background is the same as the effect estimated in the original DH population. An additional file shows the comparison between the additive effects of the Avalon × Cadenza DH population and the NILs for PH and HD (see Additional file 4). Similar QTL effects were observed in most of the cases, with the exception of the 2D QTL for PH and 1D and 3A QTLs HD. These discrepancies can be explained by QTL × E interaction, QTL× background interaction and weak phenotypic effects [51, 52]. These results agree with  who detected QTL × Env interactions on 1D and 3A for HD. Evidence of epistatic interactions were found in the original DH population (unpublished data) which can also explain these discrepancies. Further genotyping of the NILs should provide more evidence to clarify these interactions, observed as background effects in this study.
Common effects among traits
In our study, GRYLD was driven by GRpsqm as all the QTLs that showed significant difference between alleles for GRYLD were also significant for GRpsqm (1B, 3A and 7B). This was consistent with the highly positive correlation coefficient between both traits (0.80 in 2013 and 0.89 in 2014). Same effect direction was observed for GRYLD and GRpsqm. For GRYLD, differences were almost significant (dropping below significant threshold) on 1D QTL. It was interesting that the increasing allele on 1B and 3A was Cadenza whereas on 1D and 7B the increasing allele was Avalon. For the 1B QTL, the Avalon allele reduced GRYLD and HD with neutral effects on PH. The 1D QTL was also associated with HD suggesting that the early allele (Cadenza) would increase GRYLD with neutral effects on PH. These results are consistent with the effects detected in the Avalon × Cadenza DH population . For the 3A QTL, GRYLD difference was associated with HD and PH; the late allele (Cadenza) conferred tallness, GRYLD and GRpsqm increased with neutral effects on TGRWT. The combination of the Cadenza allele on 1D and 3A will result in an increase of GRYLD with neutral effects on HD. It has been shown  the segregation of both HD QTLs could confer a crop performance advantages. The direction of the effects for all traits agreed with the results obtained by Griffiths et al. [36, 38] and Ma et al. , suggesting possible tight linkage or pleiotropic effects. In addition, the Cadenza allele increased the roundness of the grain by decreasing grain length and the L/W ratio. QTLs for GRL and the L/W ratio in the Avalon × Cadenza DH population were found by Gegas et al. ; the Avalon allele increased grain length and the L/W ratio. The increased GRpsqm from the Cadenza allele can be explained by the increase in the grains/spikelet and spikelet/spike. On chromosome 3A, clusters of QTLs for GRYLD, yield components, PH and HD have been found in previous studies [24, 53–56]. It is not surprising that our results agreed with the previous findings as the interval defined by the SSR markers used for the NILs development on 3A is 48 cM based on the Avalon × Cadenza reference map (http://www.cerealsdb.uk.net). Ma et al.  found two closely linked GRYLD QTLs, qGY-psr-3A.1 and qGY-psr-3A.2 at around 73 cM and 87 cM, respectively. The qGY-psr-3A.1 co-located with PH (qPH-psr-3A) and HD (qEM-psr-3A) QTLs and the qGY-psr-3A.2 with GRpsqm (qGN-psr-3A) and GRL (qGRL-psr-3A) QTLs with the Cadenza alleles increasing all the traits except GRL. Therefore, selecting for the Avalon allele for the qGY-psr-3A.1 and for the Cadenza allele for the qGY-psr-3A.2 should provide an elite line with early flowering, shorter stature and increased in GRYLD and GRpsqm with round grains. However, high-resolution mapping together with the development of recombinants will be required to break this linkage between these QTLs or confirm pleiotropic effects. For the 7B QTL, the Avalon allele increased GRYLD through an increased in GRpsqm with neutral effects on PH and HD.
In our study, no common QTLs were found between GRYLD and TGRWT. Only the 6A QTL showed significant TGRWT effects linked to GRpsqm. Our data suggested that TGRWT is mainly explained by GRA and GRW on the 6A QTL. The early allele (Avalon) associated with TGWRT via increasing GRA, GRL and GRW and reduced GRpsqm by reducing spikelet/spike (s/S) and grains/spikelet (G/S) with some height penalty. The SSR markers (barc23 to gwm570) used for the foreground selection on 6A include the TGWRT QTL reported by Simmonds et al.  in the Spark × Rialto population; however, in this study it was expressed as a GRYLD QTL. Griffiths et al.  detected common QTLs for PH and HD on 6A with effects in opposite directions which support our results. For the 5A QTL, the Cadenza allele increased TGWRT by increasing GRA, GRL and GRW with neutral effects for PH and HD penalty. In the Avalon × Cadenza DH population, QTLs controlling TGRWT, GRpsqm, GRW, GRL, PH and HD were detected on chromosome 5A .
In the Avalon background, the introgression of the five target regions on chromosomes 1D, 3A, 5A, 6A and 7B should provide an elite line with increased GRYLD (1.87 T/Ha), TGRWT (8.56 g), Grpsqm (5364 grains/m2), GRA (2.49 mm2), GRL (0.36 mm), PH (14.73 cm) and HD (94° days to heading). In the Cadenza background, the introgression of 3A and 6A target regions will increase GRYLD (1.30 T/Ha), TGRWT (3.43 g), Grpsqm (2830 grains/m2), GRA (1.04 mm2), GRL (0.21 mm), PH (9.99 cm) and HD (30° days to heading). An alternative approach to take advantage of QTL information was proposed by Ma et al. (2015). They used a ‘breeding by design’ simulated approach to define the best two parental DH lines and selection methods in order to develop a target genotype that performed well to UK conditions. These crosses are now being conducted and will verify the simulation results.
In the present study we validate eleven QTLs determining PH, HD and GRYLD previously identified in the UK reference DH population Avalon × Cadenza in two different backgrounds. The development of 553 NILs provides a unique opportunity to assess the ‘real’ QTL effects and epistatic interactions estimated by looking at the background effects which will allow breeders to define favourable allelic combinations to develop superior genotypes adapted to a target environment. In addition, these results provides an important first step to further fine map of these traits.
DH, doubled haploid; EL, ear length; Eps, earliness per se; G/S, grains/spikelet; GRA, grain area; GRL, grain length; GRpsqm, grain number; GRW, grain width; GRYLD, grain yield; HD, heading date; HI, harvest index; L/W, grain shape; MAS, marker-assisted selection; NILs, near isogenic lines; PCA, principal component analysis; PH, crop height; PL, peduncle; Ppd, photoperiod; QTLs, quantitative trait loci; s/S, spikelet/spike; S, spikes/m2; SSR, microsatellite; TGRWT, thousand grain weight; Vrn, vernalization; WGIN, Wheat Genetic Improvement Network
Zikhali M, Leverington-Waite M, Fish L, Simmonds J, Orford S, Wingen LU, Goram R, Gosman N, Bentley A, Griffiths S. Validation of a 1DL earliness per se (eps) flowering QTL in bread wheat (Triticum aestivum). Mol Breeding. 2014;34:1023–33.
Simmonds J, Scott P, Leverington-Waite M, Turner AS, Brinton J, Korzun V, Snape J, Uauy C. Identification and independent validation of a stable yield and thousand grain weight QTL on chromosome 6A of hexaploid wheat (Triticum aestivum L.). Bmc Plant Biol. 2014;14:191.
Slafer GA, Calderini DF, Miralles DJ. Yield components and compensation in wheat: opportunities for further increasing yield potential. In: Reynolds MP, Rajaram S, McNab A, (Eds). Increasing Yield Potential in Wheat: Breaking the Barriers. Mexico DF: CIMMYT. 1996:101–34.
Miralles DJ, Slafer GA. Sink limitations to yield in wheat: how could it be reduced? J Agric Sci. 2007;145:139–49.
Borras L, Slafer GA, Otegui ME. Seed dry weight response to source-sink manipulations in wheat, maize and soybean: a quantitative reappraisal. Field Crop Res. 2004;86:131–46.
Fischer RA. Understanding the physiological basis of yield potential in wheat. J Agric Sci. 2007;145:99–113.
Peltonen-Sainio P, Kangas A, Salo Y, Jauhiainen L. Grain number dominates grain weight in temperate cereal yield determination: Evidence based on 30 years of multi-location trials. Field Crop Res. 2007;100:179–88.
Slafer GA, Savin R, Sadras VO. Coarse and fine regulation of wheat yield components in response to genotype and environment. Field Crop Res. 2014;157:71–83.
Sadras VO, Slafer GA. Environmental modulation of yield components in cereals: Heritabilities reveal a hierarchy of phenotypic plasticities. Field Crop Res. 2012;127:215–24.
Calderini D, Quintero A, Hasan A, Wingen L, Griffiths S, Girousse C, Le Gouis J. Improving grain filling and potential grain size. Proceedings of the international workshop of the wheat yield consortium. 2013; 3 (CENEB, CIMMYT. Cd. Obregon, Sonora (Mexico)): 39.
Griffiths S, Wingen L, Pietragalla J, Garcia G, Hasan A, Miralles D, Calderini D, Ankleshwaria J, Leverington-Waite M, Simmonds J, Snape J, Reynolds M. Genetic Dissection of Grain Size and Grain Number Trade-Offs in CIMMYT Wheat Germplasm. Plos One. 2015;10:118847.
Bennett D, Izanloo A, Reynolds M, Kuchel H, Langridge P, Schnurbusch T. Genetic dissection of grain yield and physical grain quality in bread wheat (Triticum aestivum L.) under water-limited environments. Theor Appl Genet. 2012;125:255–71.
Bennett D, Reynolds M, Mullan D, Izanloo A, Kuchel H, Langridge P, Schnurbusch T. Detection of two major grain yield QTL in bread wheat (Triticum aestivum L.) under heat, drought and high yield potential environments. Theor Appl Genet. 2012;125:1473–85.
Börner A, Buck-Sorlin GH, Hayes PM, Malyshev S, Korzun V. Molecular mapping of major genes and quantitative trait loci determining flowering time in response to photoperiod in barley. Plant Breeding. 2002;121:129–32.
Gegas VC, Buck-Sorlin GH, Hayes PM, Malyshev S, Korzun V. A Genetic Framework for Grain Size and Shape Variation in Wheat. Plant Cell. 2010;22:1046–56.
Groos C, Robert N, Bervas E, Charmet G. Genetic analysis of grain protein-content, grain yield and thousand-kernel weight in bread wheat. Theor Appl Genet. 2003;106:1032–40.
Huang XQ, Cloutier S, Lycar L, Radovanovic N, Humphreys DG. Molecular detection of QTLs for agronomic and quality traits in a doubled haploid population derived from two Canadian wheats (Triticum aestivum L.). Theor Appl Genet. 2006;113:753–66.
Huang XQ, Kempf H, Ganal MW, Röder MS. Advanced backcross QTL analysis in progenies derived from a cross between a German elite winter wheat variety and a synthetic wheat (Triticum aestivumL.). Theor Appl Genet. 2004;109:933–43.
Kirigwi FM, Ginkel MV, Brown-Guedira G, Gill BS, Paulsen GM, Fritz AK. Markers associated with a QTL for grain yield in wheat under drought. Mol Breeding. 2007;20:401–13.
Kuchel H, Williams KJ, Langridge P, Eagles HA, Jefferies SP. Genetic dissection of grain yield in bread wheat. I. QTL analysis. Theor Appl Genet. 2007;115:1029–41.
Maphosa L, Langridge P, Taylor H, Chalmers KJ, Bennet D, Kuchel H, Mather DE. Genetic control of grain yield and grain physical characteristics in a bread wheat population grown under a range of environmental conditions. Theor Appl Genet. 2014;127:1607–24.
Marza F, Bai G, Carver B, Zhou W. Quantitative trait loci for yield and related traits in the wheat population Ning7840 x Clark. Theor Appl Genet. 2006;112:688–98.
McIntyre CL, Mathew KL, Rattey A, Chapman SC, Drenth J, Ghaderi M, Reynolds M, Shorter R. Molecular detection of genomic regions associated with grain yield and yield-related components in an elite bread wheat cross evaluated under irrigated and rainfed conditions. Theor Appl Genet. 2010;120:527–41.
Rustgi S, Shafqat MN, Kumar N, Baenziger PS, Ali ML, Dweikat I, Campbell BT, Gill KS. Genetic Dissection of Yield and Its Component Traits Using High-Density Composite Map of Wheat Chromosome 3A: Bridging Gaps between QTLs and Underlying Genes. Plos One. 2013;8, e70526.
Sun XY, Wu K, Zhao Y, Kong FM, Han GZ, Jiang HM, Huang XJ, Li RJ, Wang HG, Li SS. QTL analysis of kernel shape and weight using recombinant inbred lines in wheat. Euphytica. 2009;165:615–24.
Tsilo TJ, Hareland GA, Simsek S, Chao SM, Anderson JA. Genome mapping of kernel characteristics in hard red spring wheat breeding lines. Theor Appl Genet. 2010;121:717–30.
Varshney RK, Prasad M, Roy JK, Kumar N, Harjit S, Dhaliwal HS, Balyan HS, Prasad HS, Gupta PK. Identification of eight chromosomes and a microsatellite marker on 1AS associated with QTL for grain weight in bread wheat. Theor Appl Genet. 2000;100:1290–4.
Xiang-Zheng L, Jin W, Rong-Hua Z, Zheng-Long R, Ji-Zeng J. Mining favorable alleles of QTLs conferring thousand-grain weight from synthetic wheat. Acta Agrnomica Sinica. 2008;34:1877–84.
Flintham JE, Borner A, Worland AJ, Gale MD. Optimizing wheat grain yield: Effects of Rht (gibberellin-insensitive) dwarfing genes. J Agric Sci. 1997;128:11–25.
Miralles DJ, Katz SD, Colloca A, Slafer GA. Floret development in near isogenic wheat lines differing in plant height. Field Crop Res. 1998;59:21–30.
Zanke CD, Ling J, Plieske J, Kollers S, Ebmeyer E, Korzun V, Argillier O, Stiewe G, Hinze M, Neumann K, Ganal MW, Röder MS. Whole Genome Association Mapping of Plant Height in Winter Wheat (Triticum aestivum L.). Plos One. 2014;9:113287.
Brandle JE, Knott DR. The effect of a gene for semidwarfism (Rht1) on various characters in a spring wheat cross. Can J Plant Sci. 1986;66:529–33.
Gasperini D, Greenland A, Hedden P, Dreos R, Harwood W, Griffiths S. Genetic and physiological analysis of Rht8 in bread wheat: an alternative source of semi-dwarfism with a reduced sensitivity to brassinosteroids. J Exp Bot. 2012;63:4419–36.
Worland AJ, Law CN. Genetic-analysis of chromosome 2D of wheat.1. The location of genes affecting height, day-length insensitivity, hybrid dwarfism and yellow-rust resistance. Zeitschrift Fur Pflanzenzuchtung-Journal of Plant Breeding. 1986;96:331–45.
Bai C, Liang Y, Hawkesford MJ. Identification of QTLs associated with seedling root traits and their correlation with plant height in wheat. J Exp Bot. 2013;64:1745–53.
Griffiths S, Simmonds J, Leverington M, Wang Y, Fish L, Sayers L, Alibert L, Orford S, Wingen L, Snape J. Meta-QTL analysis of the genetic control of crop height in elite European winter wheat germplasm. Mol Breeding. 2012;29:159–71.
Ma J, Wingen LU, Orford S, Fenwick P, Wang J, Griffiths S. Using the UK reference population Avalon x Cadenza as a platform to compare breeding strategies in elite Western European bread wheat. Mol Breeding. 2015;35:70.
Griffiths S, Simmonds J, Leverington M, Wang Y, Fish L, Sayers L, Alibert L, Orford S, Wingen L, Herry L, Faure S, Laurie D, Bilham L, Snape J. Meta-QTL analysis of the genetic control of ear emergence in elite European winter wheat germplasm. Theor Appl Genet. 2009;119:383–95.
Hanocq E, Niarquin M, Heumez E, Rousset M, Le Gouis J. Detection and mapping of QTL for earliness components in a bread wheat recombinant inbred lines population. Theor Appl Genet. 2004;110:106–15.
Kamran A, Iqbal M, Spaner D. Flowering time in wheat (Triticum aestivum L.): a key factor for global adaptability. Euphytica. 2014;197:1–26.
Le Gouis J, Bordes J, Ravel C, Heumez E, Faure S, Praud S, Galic N, Remoue C, Balfourier F, Allard V, Rousset M. Genome-wide association analysis to identify chromosomal regions determining components of earliness in wheat. Theor Appl Genet. 2012;124:597–611.
Reif JC, Maurer HP, Korzun V, Ebmeyer E, Miedaner T, Würschum T. Mapping QTLs with main and epistatic effects underlying grain yield and heading time in soft winter wheat. Theor Appl Genet. 2011;123:283–92.
Rousset M, Bonnin I, Remoué C, Falque M, Rhoné B, Veyrieras JB, Madur D, Murigneux A, Balfourier F, Le Gouis J, Santoni S, Goldringer I. Deciphering the genetics of flowering time by an association study on candidate genes in bread wheat (Triticum aestivum L.). Theor Appl Genet. 2011;123:907–26.
Payne RW, Murray DA, Harding SA, Baird DB, Soutar DM. Genstat for windows (16th Edition) Introduction. VSN International, Hemel Hempstead. 2013.
Asseng S, Foster I, Tuner NC. The impact of temperature variability on wheat yields. Glob Chang Biol. 2011;17:997–1012.
Sofield I, Evans LT, Cook MG, Wardlaw IF. Factors influencing rate and duration of grain filling in wheat. Aust J Plant Physiol. 1977;4:785–97.
Wardlaw IF, Sofield I, Cartwright PM. Factors limiting the rate of dry-matter accumulation in the grain of wheat grown at high-temperature. Aust J Plant Physiol. 1980;7:387–400.
Zahedi M, Jenner CF. Analysis of effects in wheat of high temperature on grain filling attributes estimated from mathematical models of grain filling. J Agric Sci. 2003;141:203–12.
Fisher RA. Selection traits for improving yiled potential. In: Reynolds MP, Ortiz-Monasterio JI, McNab A, editors. Application of physiology in wheat breeding. Mexico, DF: CIMMYT; 2001.
Langer SM, Longinand CFH, Wuerschum T. Flowering time control in European winter wheat. Front Plant Sci. 2014; doi:10.3389/fpls.2014.00537.
Bouchez A, Hospital F, Causse M, Gallais A, Charcosset A. Marker-assisted introgression of favorable alleles at quantitative trait loci between maize elite lines. Genetics. 2002;162:1945–59.
Thabuis A, Palloix A, Servin B, Daubeze AM, Signoret P, Hospital F, Lefebvre V. Marker-assisted introgression of 4 Phytophthora capsici resistance QTL alleles into a bell pepper line: validation of additive and epistatic effects. Mol Breeding. 2004;14:9–20.
Ali ML, Baenziger PS, Ajlouni ZA, Campbell BT, Gill KS, Eskridge KM, Mujeeb-Kazi A, Dweikat I. Mapping QTL for Agronomic Traits on Wheat Chromosome 3A and a Comparison of Recombinant Inbred Chromosome Line Populations. Crop Sci. 2011;51:553–66.
Campbell BT, Baenziger PS, Gill KS, Eskridge KM, Budak H, Erayman M, Dweikat I, Yen Y. Identification of QTLs and environmental interactions associated with agronomic traits on chromosome 3A of wheat. Crop Sci. 2003;43:1493–505.
Dilbirligi M, Erayman M, Campbell BT, Randhawa HS, Baenziger PS, Dweikat I, Gill KS. High-density mapping and comparative analysis of agronomically important traits on wheat chromosome 3A. Genomics. 2006;88:74–87.
Shah MM, Gill KS, Baenzinger PS, Yen Y, Kaeppler S, Ariyarathne. Molecular mapping of loci for agronomic traits on chromosome 3A of bread wheat. Crop Sci. 1999; 39: 1728–1732.
All contributors are acknowledged by authorship.
This work was funded by the Department for Environment, Food and rural Affairs as part of the UK Wheat Genetic Improvement Network (WGIN).
A.F. wrote the manuscript and performed the statistical analysis; L.S. carried out the first two rounds of backcrossing; M.LW. led MAS work; R.G. supported MAS work; L.W. carried out preliminary statistics analysis; S.O. conducted phenotyping; C.M. conducted phenotyping; S.G. supervised, designed the study and received the funding. All autors have read and approved the manuscrit.
The autors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Means of minimum (blue line) and maximum (red line) temperatures and cumulated rainfall (blue bars) during the months of the growth cycle for both trials. (PDF 95 kb)
Difference between Avalon and Cadenza alleles for the ratio between the lenghts of each PH component and the total PH for the chromosome that showed significant differences for plant height (values above zero indicates that the Avalon allele increased height whereas below zero the Cadenza allele reduced height; ear length (EL), peduncle (PL) and internode lengths (from the first to the fifth counted from the top, abbreviated as 1stITL, 2ndITL, 3rdITL, 4thITL, 5thITL); * P-value < 0.05; ** P-value <0.01 and *** P-value < 0.001). (PDF 165 kb)
(a) Summary of the mixed model analysis performed for the grain number components studied comprising 553 NILs carrying the Avalon or Cadenza alleles in the introgressed region in 2013. (b) Average values for the two groups (carrying the Avalon or Cadenza alleles in the QTL region) based on the chromosome and background in 2013. Significant difference between Avalon and Cadenza alleles are highlighted in bold (spikes/m2 (S), spikelet/spike (s/S) and grains/spikelet (G/S). (PDF 111 kb)
Comparions between the magnitude of the effects of the A × C DH population (red) and the NILs for heading date (a) and plant heigth (b) (2013 and 2014 are indicated in bakc and grey color, respectively; the Avalon background is showed in dashed line whereas the Cadenza background is showed in straigth line). (PDF 91 kb)
About this article
Cite this article
Farré, A., Sayers, L., Leverington-Waite, M. et al. Application of a library of near isogenic lines to understand context dependent expression of QTL for grain yield and adaptive traits in bread wheat. BMC Plant Biol 16, 161 (2016). https://doi.org/10.1186/s12870-016-0849-6
- Near isogenic lines