|
|
||||||||
a Embrapa, C.P. 179, Santo Antonio de Goias, GO, 75375, Brazil
b Dep. of Plant Breeding and Genetics, 240 Emerson Hall, Cornell Univ., Ithaca, NY 14853
c Roman Meal Co., 2101 S. Tacoma Way, Tacoma, WA 98409
d USDA, Soft Wheat Quality Lab., Williams Hall, 1680 Madison Ave. Wooster, OH 44691
e NeoVentures Biotechnology Inc., 69 Mary Street, Guelph, ON, Canada, N1G 2A9F
* Corresponding author (mes12{at}cornell.edu)
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: AWRC, Alkaline water retention capacity CIM, Composite interval mapping SMR, Single marker regression analysis
| INTRODUCTION |
|---|
|
|
|---|
Hybridizations between hard and soft wheat types could expand the genetic base of wheat breeding and create new possibilities for combinations of desirable alleles from both germplasm subgroups. However, this type of cross is not common practice in wheat breeding because the two classes have distinct quality goals. Carver (1996) compared interclass hybrids, backcrosses and progeny from a hard x hard cross, and concluded that the interclass crosses resulted in progenies with higher grain yield but lower flour yield and larger variability for quality traits, and that recovering the quality profile of the hard type through intensive selection would be feasible. Identification of quantitative trait loci (QTL) related to quality differences between classes could help in planning complementary crosses and backcrosses, and in designing selection schemes to recover the quality characteristics needed in either class.
Major genes controlling the difference in kernel texture between soft and hard wheat were mapped on the short arm of chromosome 5D (Sourdille et al., 1996), and as a group are named Ha or Hardness locus (Baker, 1977). The product of that locus is called friabilin, which is a composite of related proteins that include puroindoline a and puroindoline b (Giroux and Morris, 1998). Other QTLs, with minor effects on kernel texture within the hard wheat type, were detected on chromosomes 1A and 6D (Perretant et al., 2000).
Populations derived from crosses between hard and soft genotypes have higher expected marker polymorphism because of the divergence between the two breeding groups, which can facilitate building a linkage map, compared with other elite x elite crosses. Additionally, there is potential to discover QTLs that are fixed for different alleles within each texture class, which could not be detected in a hard x hard or soft x soft cross. A population derived from the cross of the soft elite line NY18 and the hard variety Clark's Cream has been used to map QTLs of grain quality traits by Campbell et al. (2001). These authors found protein QTLs on groups 1, 2, and 7, and a flour yield QTL on group 3. Markers related to high molecular weight glutenins had major influence on mixogram traits.
In the study of Campbell et al. (2001), the RFLP marker PinB, linked to Ha-5DS, was highly significantly associated with flour yield, softness equivalent, damaged starch, AWRC, and cookie diameter. However, the effect of the Hardness gene was disproportionately large compared with the other loci influencing grain quality, resulting in bimodal distribution of traits affected by that gene. Many statistical methods used in QTL analysis assume normal distribution (Lynch and Walsh, 1998). One approach to meet this assumption when analyzing quantitative variation mixed with qualitative variation is to regress the quantitative trait of interest on the categorical trait and use residuals as phenotypes in the analysis. That approach has been used in genetic analysis of humans and animals to account for effects of sex or race (Grosz and MacNeil, 2001; Friedlander et al., 2003). The effect of a major gene segregating in a mapping population is analogous to those cases and could be controlled by using similar means, provided that categorical variation can be clearly identified.
This paper reports the results of a QTL analysis of several milling and baking quality traits, conducted with the objective of identifying genomic regions related to the variation in kernel quality in a population derived from a cross between a hard and a soft wheat variety. Our focus was on loci other than the known Hardness locus; hence, the effect of the segregation of this locus on quality traits was quantified and controlled by regression analysis.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Genotypic Data and Linkage Map Construction
The genetic map was constructed with 340 molecular markers, including 222 AFLP, 42 RFLP, 75 SSR, and one STS (Lox, Hessler et al., 2002). AFLP markers were evaluated at Agriculture Canada by a modified method based on Vos et al. (1995). All the other markers were evaluated at Cornell University. Hybridization probes included 22 CDO (oat cDNA) and 9 BCD (barley cDNA, Heun et al., 1991), 7 RZ (rice cDNA, Causse et al., 1994), and markers of known genes: PinA (puroindoline a, Giroux and Morris, 1998), MTA9 (esterase, Jouve and Diaz, 1990), and A1tgh (dihydroflavanol reductase, R.A. Graybosch, pers. comm.). SSR markers included 31 WMC (Gupta et al., 2002), 30 GWM (Röder et al., 1998), and 14 BARC (P. Cregan, Q. Song and coll., unpublished). PCR reactions were modified from Röder et al. (1998), and PCR runs consisted of 5 min at 94°C, 35 cycles of 45 s at 94°C, 45 s at specific marker annealing temperature (www.graingenes.org; verified 15 April 2005), and 90 s at 72°C, followed by a 10 min final extension at 72°C.
The linkage map was constructed with Map Manager QTX13 (Manly and Olson, 1999) using the Kosambi mapping function and type-I error probability of 0.001, followed by "ripple" command to check the ordering of markers within each linkage group. Genetic markers presenting high segregation distortion by the
2 test (P < 0.001) were not used in the first stage of map construction, when the core map was defined. In a second round, those markers were allowed to link to the existing groups using the "distribute" command, so that they could be used for QTL analysis by composite interval mapping. The resulting linkage groups were compared to the Synthetic W7984 x Opata 85 map (Marino et al., 1996; Nelson et al., 1995; Van Deynze et al., 1995) and tentatively assigned to wheat chromosomes on the basis of common markers. Mapchart (Voorips, 2002) was used to draw the linkage map with QTLs.
Phenotypic Evaluations
Seed was thoroughly air-aspirated to remove shriveled kernels and nonwheat materials before evaluation for milling and baking quality parameters at the USDA-ARS, Soft Wheat Quality Laboratory, at Wooster, OH. Four traits were evaluated on samples from three locations: Quadrumat mill flour yield and softness equivalent, according to Finney and Andrews (1986), flour protein content and alkaline water retention capacity (AWRC), following AACC Approved Methods 46-12 and 56-10, respectively (www.aaccnet.org; verified 15 April 2005). A sample of 500 g of grain from Tulelake was Allis-Chalmers milled, yielding Allis-Chalmers flour yield, break-flour yield, endosperm separation index (Yamazaki and Andrews, 1977), and friability. Friability was the ratio of the weight of flour produced divided by the sum of the weight of the mill stock fed to each of the break and reduction rolls. The consistency of the Allis-Chalmers milling data were evaluated by analysis of variance of two independent samples from each of seven checks [AC Reed, Grandin, Serra (CI 87738), Yolo (CI 17961), Katepwa (AFRC 5621), Anza (CI 15284), and UC896 (University of California-Davis experimental line)]. The check samples were produced in the same year and location as the mapping population. Ten-gram mixograph assay (Finney and Shogren, 1972) resulted in mixing time and peak height, plus mixogram curve area and area under the curve determined by image analysis of the mixogram. Additional tests included: kernel volume (AACC Method 55-10), sucrose retention capacity and lactic acid retention capacity (AACC Method 56-11), cookie diameter, and top grain (AACC Method 10-52).
QTL Analysis
The DH lines were classified as hard or soft type on the basis of the mean softness equivalent score (Yamazaki and Andrews, 1977). Phenotypic data were regressed on this classification, and residual variance was used for QTL analysis. The purpose of this adjustment was to eliminate the qualitative-like variation caused by segregation of the Ha-5DS locus (Sourdille et al., 1996), and achieve normal distribution of data, as checked by Shapiro-Wilk test at
= 0.01 (SAS System V.8). QTL analysis was performed in QTL Cartographer V2.0 (Wang et al., 2001). All traits were analyzed by single marker regression analysis (SMR) using all markers.
Replicated traits (Quadrumat flour yield, softness equivalent, protein content, and AWRC) were analyzed by composite interval mapping (Zeng, 1994), using a reduced version of the map, including only linkage groups containing at least one locus significant for any trait at
= 0.01 in the previous SMR analysis. The parameter settings for CIM were model 6, stepwise selection of cofactors with SLE = SLS = 0.01, window size 10 cM and testing step 2 cM. Experiment-wise type I error rate of detected QTLs was estimated by one thousand permutations with the same settings (Churchill and Doerge, 1994).
Multiple-trait analysis method was used to jointly analyze QTL over locations using each location as a "correlated trait" (Jiang and Zeng, 1995). A threshold of LOD = 3 was arbitrarily applied for those QTLs. Approximate 95% confidence intervals were built by concatenating adjacent positions within a 1-LOD difference from the value of the peak, although those intervals tend to have an actual probability level lower than 95% (Mangin et al., 1994), and the interval probability is unknown when it is truncated by the end of the linkage group.
Analysis of variance of phenotypes for locations and lines was done, and significance was calculated using the variance of location x line interaction as error. Multiple regression models were built, beginning with markers within the confidence interval of joint analysis plus the significant markers by SMR at
= 0.01. Predictors were selected by backward elimination with SLS = 0.05 for main effects and SLS = 0.01 for interaction effects (SAS System V.8).
| RESULTS |
|---|
|
|
|---|
2 p < 0.001). Thirty-seven loci were skewed toward AC Reed, most of them mapped to the linkage group 1AD. The largest skewed region toward Grandin was on 1BL-2, including 6 consecutive markers.
|
|
|
|
|
|
Protein Content
Three QTLs for protein content were detected by CIM at Tulelake, on linkage groups 2AS, 4B-1 and 6B, and two QTLs were detected in Swift Current, on 2BL-1 and 4B-2 (Table 2). All of those regions were significant by SMR (Fig. 3). No QTL was found for Lethbridge by CIM, and the single markers were only marginally significant. Joint analysis detected QTLs on linkage groups 1BL-2, 2AS, 2DS-2, 4B-1, 5BL and in 6B, which was the most significant (LOD = 4.9, Fig. 1). The main effect of location in the ANOVA of protein content was very small (r2 = 0.022), but the location x line interaction was the highest among the traits studied (r2 = 0.296). A multiple regression model included seven markers but no interaction effects. Together, those markers explained proportions of the variance varying from 36.5% at Lethbridge to 65.6% at Tulelake (Table 3).
Alkaline Water Retention Capacity
AWRC was correlated with markers on 1BL-2 and 4B-3, plus other minor effects (Fig. 3). CIM detected the QTL on 4B-3 for Swift Current, which had the highest coefficient of determination in our study (r2 = 0.403, Table 2). AC Reed, the lower parent for AWRC, contributed the increasing allele at that QTL. Joint analysis of environments detected the QTL on 4B-3 with LOD = 6.1, plus three other significant regions on 1BL-2, 4B-2, and 6AL-2 (Fig. 1). In the ANOVA, main effect of location and the location x line interaction on AWRC was intermediate compared with the other traits. A regression model composed of three markers explained a modest proportion of the phenotypic variance (Table 3).
Allis-Chalmers Milling Traits
The Allis-Chalmers milling data for the population were unreplicated. An approximate error associated with those determinations was derived from two replicates of each of seven checks, grown in the same field and year. The error variance of all of the Allis-Chalmers traits was very small (Table 4), with cultivars accounting for more than 95% of the variation for all traits except kernel volume. The standard deviations among lines for the same traits (Table 3) were two- to four-fold greater than the error among checks, implying that a large proportion of the phenotypic variance was due to genetic differences among lines. The locus Xgcat6 on 1AD was highly significant (P < 0.001) for Allis-Chalmers milling flour yield and friability and marginally significant (P < 0.05) for break flour yield and endosperm separation index (Fig. 3). AC Reed had the favorable allele at that locus, increasing flour yield and softness. Another region influencing friability was found on 1BL-2, with increasing effect from Grandin. A QTL for break flour yield was detected by Xbcd808 on 1AS,L, in a region apparently involved with softness. Three highly significant markers for kernel volume were found on 2ABS, 2BS,L, and 4B-3. Grandin contributed the increasing alleles for those three loci (Fig. 3).
|
Baking Assay Traits
Lactic acid retention capacity correlated positively with AC Reed alleles on 1AS (Fig. 3). A region affecting sucrose retention capacity was detected on 4B-2, with the increasing effect from Grandin. Cookie diameter was influenced by a QTL on 6B, with positive effect from Grandin. Top grain was related to markers on 2ABS and 4B-3, both with favorable effect from AC Reed.
| DISCUSSION |
|---|
|
|
|---|
The D-genome was populated with fewer markers as a consequence of lower polymorphism compared with A and B genomes, as observed in previous studies (e.g., Cadalen et al., 1997; Marino et al., 1996). More complete linkage maps normally require populations derived from wide crosses; however, this type of population frequently presents limitations to QTL analysis because of lack of agronomic adaptation. Adapted genotypes and their progeny are more likely to express their genetic potential under normal growing conditions and to produce fully developed grain for quality assessment. In this study, both parents were highly adapted to the growing conditions of the tested locations.
Hard and soft wheat varieties are normally developed by separate breeding programs; hence, the polymorphism of the AC Reed x Grandin cross was probably higher than that of a cross within texture class. However, the segregation of the hardness locus Ha-5DS introduced a "nuisance" variance in the population. This was a similar situation to human or animal populations where the effect of sex or race can override the expression of genes of interest (Grosz and MacNeil, 2001; Friedlander et al., 2003). Removing the effect of Ha by regression facilitated the detection of minor QTLs.
Single marker regression is the simplest method of QTL analysis. Although this method has some disadvantages, such as low power and biased estimation of QTL effects (Lander and Botstein, 1989), its results allowed a thorough comparison of locus-wise significance across traits. QTLs detected through CIM above thresholds that are corrected for multiple testing normally do not allow such straightforward comparison because many moderate but nonetheless real effects are not detected. A conservative threshold may result in different QTLs being detected in each location and could overestimate the importance of the QTL x environment interaction. The use of lower thresholds would probably result in better agreement among different environments and studies, but in this case, the protection of the experiment-wise type I error rate would require larger population sizes. On the other hand, CIM allows for a more precise location of important QTLs within linkage groups and better estimation of the QTL effects (Zeng, 1994). Therefore, the results of the two methods complemented each other in this study.
Rapid microtests for milling quality that identify the same genes as commercial scale milling evaluations are critical for wheat millers and breeders. Previous reports described a modified Quadrumat mill that produced flour yield data that ranked cultivars similarly to the longer-flow Allis-Chalmers mill (Finney and Andrews, 1986) but noted that grain moisture content and texture strongly influenced the results. Gaines et al. (2000) developed algorithms to adjust Quadrumat flour yield to 150 g kg1 moisture and a constant softness equivalent based on particle size, which raised the correlation between the results from the two mills from 55 to 90%, without having to temper the wheat before Quadrumat milling. However, phenotypic correlations often do not accurately reflect a common genetic basis for two traits. In this study, we compared the QTL for Quadrumat flour yield with the Allis-Chalmers flour yield (Fig. 3). All markers that were significant for Quadrumat flour yield were also significant for either Allis-Chalmers flour yield or friability or both and Allis-Chalmers flour yield identified only one QTL that was not detected by the Quadrumat (Xbarc61-1BL-1). These results validate the use of the modified Quadrumat mill and the algorithms of Gaines et al. (2000) for predicting flour yield and friability of Allis-Chalmers mill for selecting wheat genotypes with superior milling characteristics.
Some of the QTLs reported here may have a common genetic base with previously reported QTLs. The locus Xbcd1431 on 4DL-2 was related to AWRC and damaged starch in the study of Campbell et al. (2001) and may be homoeologous to the QTL for AWRC (Swift Current) on 4B-3 in this study. Our QTL for kernel volume on 2ABS agreed with a major QTL reported by Campbell et al. (1999). The largest QTL for Quadrumat flour yield on 3ABS may be the same reported on 3A in bread wheat by Parker et al. (1999). QTLs for softness equivalent were detected here on linkage group 1AS,L, which may match a QTL for kernel hardness previously reported on 1A (Perretant et al., 2000).
The highly significant QTL for protein content that we found on 2AS could refer to the same gene as one reported on chromosome 2A of other bread wheat populations by Groos et al. (2003) and Prasad et al. (2003). A QTL for protein content was mapped close to Xgwm193-6B by Khan et al. (2000). In our map, that marker is between two non-overlapping confidence intervals of protein QTLs on 6B. However, a more conservative criterion (e.g., 2-LOD drop) would merge the two confidence intervals (result not shown), hence the possibility of a single QTL cannot be eliminated. We found QTLs for protein content on 4B-1 (Tulelake) and 4B-2 (Swift Current) and confirmed by joint analysis on 4B-1, which could match the QTL on 4B reported to be the most stable across environments by Blanco et al. (2002).
Some mixogram QTLs agreed with previous QTLs or known genes. On the basis of the linked loci Xrz166 and Xcdo1173, we infer that the effect detected by Xwmc216a-1BS,L was probably caused by the gliadin Gli-B3 gene (Galili and Feldman, 1984). The gene underlying the QTL for mixogram time and curve area near Xbcd738-1BL-1 was probably the glutenin Glu-B1 gene (Payne et al., 1982), which is near the locus Xbarc61 (5 cM from Xbcd738) in the W7984 x Opata 85 map. The locus Xagat12, near the QTL for mixogram height on 3ABS, was closely linked to Xcdo718, which was related to flour viscosity in the population NY18 x Clark's Cream (Udall et al., 1999).
Several QTLs had effects opposite to the phenotypes of the parents. This was expected because parental phenotypes were largely defined by their alleles for Ha-5DS. Grandin, the hard wheat parent, contributed the allele for softness at a QTL on 1AS,L (Table 2) and for larger cookie diameter on 6B. Similarly, the Grandin allele at the QTL near Xwmc44-1BL-2 was favorable for soft wheat products, increasing friability, flour yield, cookie diameter, and top grain, while decreasing protein content, AWRC, sucrose retention capacity, and mixogram height (Fig. 3). On the other hand, AC Reed alleles increased protein content at the QTLs on 2BL-1, 4B-1, and 6B, and flour yield at QTLs on 2BL-2 and 4B-2 (Table 2). Those sources of variation represent opportunities to improve quality traits through hard x soft hybridization. QTLs reported here and in previous papers could orient marker-assisted selection strategies to accelerating the recovery of quality characteristics required for each class following hybridizations while retaining favorable alleles from both parents.
| NOTES |
|---|
|
|
|---|
Received for publication May 19, 2004.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
E.J. Souza, C. Griffey, M. Kweon, and M. Guttieri Sources of Variation for Long-Flow Experimental Milling Crop Sci., July 1, 2008; 48(4): 1432 - 1440. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. Li, L. Huang, and B. S. Gill Recurrent Deletions of Puroindoline Genes at the Grain Hardness Locus in Four Independent Lineages of Polyploid Wheat1,[W],[OA] Plant Physiology, January 1, 2008; 146(1): 200 - 212. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Agronomy Journal | Vadose Zone Journal | |||
| Journal of Natural Resources and Life Sciences Education |
Soil Science Society of America Journal | ||||
| Journal of Plant Registrations | Journal of Environmental Quality |
The Plant Genome | |||