|
|
||||||||
a Dep. of Crop Sci., Univ. of Illinois, Urbana, IL 61801
b Dep. of Food Sci. and Human Nutrition, Univ. of Illinois, Urbana, IL 61801
c USDA-ARS, Soybean/Maize Germplasm, Pathology, and Genetics Research Unit, Dep. of Crop Sci., Univ. of Illinois, Urbana, IL 61801
d Dep. of Agronomy, Iowa State Univ., Ames, IA 50011
e Dep. of Agronomy, Purdue Univ., West LaFayette, IN 47907
* Corresponding author (ekabelka{at}uiuc.edu).
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: LOD, likelihood of odds MG, maturity group PCR, polymerase chain reaction QTL, quantitative trait locus/loci R1, days to flower R8, days to maturity RIL, recombinant inbred line SSR, simple sequence repeat
| INTRODUCTION |
|---|
|
|
|---|
Introgressing exotic germplasm into cultivars to increase genetic diversity within domesticated crops has been used to enhance complex traits such as yield (Tanksley and McCouch, 1997). In soybean, Thompson and Nelson (1998) tested experimental lines derived from crossing North American cultivars with several plant introductions. Lines that yielded significantly more than their North American parent were identified which indicated that soybean PIs may have alleles for enhanced seed yield.
Molecular marker analysis can aid in the discovery and mapping of QTL associated with beneficial and novel alleles from exotic parents. For example, in soybean, two populations of recombinant inbred lines (RILs) derived from Minsoy and Archer, and Noir 1 and Archer (Orf et al., 1999) have been used extensively for the identification of QTL governing agronomic and seed composition traits utilizing molecular markers (Orf et al., 1999; Terry et al., 2000; VanToai et al., 2001). In one study comparing these two RIL populations, Orf et al. (1999) identified QTL for plant height, lodging, days to flowering and maturity, reproductive period, seed yield, seed weight, and seed oil and protein concentration. At many of the QTL identified in this study, the beneficial or novel alleles were those provided by Minsoy and Noir 1. Minsoy and Noir 1 are divergent in lineage from each other and from North American cultivars (Mansur et al., 1993a, 1993b; Orf et al., 1999), whereas Archer is a North American cultivar (Cianzio et al., 1991).
Yield alleles from exotic sources are not well documented in soybean and those identified have often been related to other agronomic traits such as maturity. In this study, a soybean population consisting of F5derived lines was developed from a cross between BSR 101 and the experimental line LG82-8379. The pedigree of BSR 101 has nine of 10 ancestral lines that have made the greatest contribution to the commercially used gene pool of North America (Bernard et al., 1988; Gizlice et al., 1994) whereas LG82-8379 was selected from a cross between two soybean PIs, PI 68508 and FC 04007B. Using random amplified polymorphic DNA markers, the two PIs had been shown to be genetically distinct from all major ancestral lines of modern North American soybean cultivars (Thompson et al., 1998). From this cross, high-yielding and genetically diverse experimental lines have been developed (Thompson and Nelson, 1998) and released (Thompson et al., 1999). The objective of this study was to evaluate the two PIs as sources of alleles for the enhancement of seed yield in North American cultivars.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Days to maturity, plant height, lodging, and seed yield were recorded in all trials. Days to maturity was recorded as the number of days after planting when approximately 95% of the pods had reached mature pod color (R8; Fehr et al., 1971). Plant height (mm) was measured at maturity as the average distance from soil surface to the apex of the main stem. Lodging was scored at maturity on a scale of 1 to 5, with 1 designated as plants standing erect and 5 as plants prostrate. The two center rows of each subplot were harvested and seed yield, expressed as kg ha1, was adjusted to 130 g kg1 moisture. At three locations each year, days to flower was recorded as the number of days after planting when approximately 50% of the plants had at least one flower (R1; Fehr et al., 1971). Reproductive period was calculated by subtracting R1 from R8. Seed protein and oil concentration were measured on whole beans from four locations each year using near infrared transmission (NIT) grain analyzers, models ZX800 and ZX880 (Zeltex Inc., Hagerstown, MD). All samples were measured in duplicates in order to improve reliability. The seed protein and oil concentrations are reported on a moisture-free basis as grams per kilogram.
Leaf tissue DNA was extracted from eight greenhouse-grown seedlings per F5:8 line according to Keim et al. (1988) with modifications as described by Kisha et al. (1997). The SSR markers used in this study were developed by P.B. Cregan (USDA-ARS, Beltsville, MD). Nonlabeled and fluorescently-labeled primers were obtained from Applied Biosystems (Foster City, CA) and Research Genetics (Huntsville, AL). Polymerase chain reactions (PCRs) were performed according to Cregan and Quigley (1997). The nonlabeled PCR products were analyzed by electrophoresis in 6% nondenaturing polyacrylamide gels (Sambrook et al., 1989; Wang et al., 2003) and stained with 1 µg mL1 ethidium bromide. The fluorescently-labeled PCR products were analyzed using an ABI Prism 310 Genetic Analyzer (Applied Biosystems, Foster City, CA).
More than 500 SSR markers covering the 20 chromosomes of the soybean genome were screened against the parent BSR 101 and a set of seven DNA bulks that each consisted of a mixture of 22 different F5:8 lines to identify markers segregating within this population. The bulk DNA set of F5:8 lines was used instead of LG82-8379 to identify polymorphisms because LG82-8379 is heterogeneous and the exact parental genotypes are unknown. A total of 145 polymorphic SSR markers were identified. For each polymorphic marker, only two alleles were present within the population and allelic frequencies at each locus were as expected within and across sets. Since only two alleles were found at each locus, this suggests that the three F1 plants used to create this population were homogeneous at the loci evaluated in this study. To determine their contribution, PI 68508 and FC 04007B were each genotyped with SSR markers associated with yield enhancement.
Agronomic traits were subjected to ANOVA by the PROC GLM and VARCOMP functions of SAS (Statistical Analysis System version 8.0, SAS Institute, Cary, NC). Sets, based on maturity, were considered as a fixed effect whereas environment, replications within environment, lines within sets, environment by sets, and environment by lines within set interactions were considered as random effects. Each locationyear combination was considered an environment in the analysis. Significance was determined by an F test as described by McIntosh (1983), with the numerator and denominator degrees of freedom approximated as described by Satterthwaite (1946). Means, ranges, and standard deviations across environments within each set were computed for each trait. Least significant differences between set means were determined at a 5% significance level. Broad-sense heritabilities of each trait were calculated on a line-mean basis according to Fehr (1987). Pearson correlation coefficients among all traits within and across sets were calculated from the means of lines across environments using PROC CORR in SAS. The experiment-wide error rate in this study was calculated according to Lander and Botstein (1989).
Single marker-trait and composite interval mapping methods were used to identify marker association with each agronomic trait. For both methods, sets were analyzed separately and combined. Single marker-trait analysis was performed by PROC GLM in SAS. Marker genotypes, used as class variables, were considered as fixed effects whereas environments, replications, and lines were considered random effects. In the combined analyses, means of all lines from the three sets were analyzed together to detect QTL. Pairwise comparisons of least square means were used to test the significance (P < 0.05) of each marker class. Lines heterozygous for the genetic markers were included in all analyses, but their means were not included in the summary tables. Composite interval mapping was performed using the multiple-regression based computer program PLABQTL (Utz and Melchinger, 1996), with the required genetic linkage map constructed with the computer program JoinMap v. 3.0 (Van Ooijen and Voorrips, 2001). A minimum LOD of 3.0 and a maximum distance of 50 cM was used for testing linkages among markers. A minimum LOD of 2.5 was chosen to declare a marker association significant. The R2 value was used to describe the phenotypic variance explained by individual markers. The total phenotypic variance explained by two or more QTL for a given trait within a single set or across sets was determined by using a multifactor ANOVA.
| RESULTS |
|---|
|
|
|---|
|
Mean squares were estimated for the 167 F5derived lines grouped into sets according to maturity and grown across 12 environments. Mean square estimates for environment, sets, lines within sets, environment by set, and environment by lines within set interactions were significant (P < 0.01) for all traits measured except for protein concentration and seed yield where sets were not significant. Variance component estimates revealed which components accounted for the variance of each trait measured (Table 2). For yield, the environment accounted for more of the variance compared with lines within sets, environment by sets, and environment by lines within set interactions.
|
|
QTL Analysis
In this study, we report QTL that were significantly associated (P < 0.05, LOD > 2.5) with agronomic and seed composition traits across 12 environments. Fifteen QTL were found significantly associated with seed yield (Table 4). At nine of the 15 QTL, the LG82-8379 alleles increased seed yield (Table 4A). Two QTL, located on linkage groups C2 and O, were found significant (P < 0.05) in combined analysis of all lines across the three sets. The alleles contributing to the increase in yield at these loci were both introgressed from FC 04007B. The linkage group C2 QTL explained 10% of the phenotypic variation and increased seed yield 60 kg ha1, while the linkage group O QTL explained 14% of the phenotypic variation and increased seed yield 47 kg ha1. Collectively, the two QTL accounted for 18% of the phenotypic variation for seed yield. The remaining seven QTL with the LG82-8379 alleles increasing seed yield were detected only within specific sets. The alleles contributing to the increase in yield on linkage groups A1, G, K, D2, and M were from PI 68508, whereas FC 04007B contributed the alleles on linkage groups B2 and H, in addition to alleles on linkage groups C2 and O. Within specific sets, these QTL explained 14 to 38% of the phenotypic variation and increased seed yield 60 to 148 kg ha1. Including the QTL that were significant across sets, the QTL collectively explained 47% of the phenotypic variation for seed yield within Set 1, 21% within Set 2, and 52% within Set 3.
|
No QTL were found significant for R8 in combined analysis of all lines across the three sets, but four QTL were detected among lines of Set 1 (Table 5). The QTL located on linkage groups A1, G, K, and N, explained 6 to 12% of the phenotypic variation for R8 and accounted for a 1 to 2 d difference between the BSR 101 and LG82-8379 alleles at these loci. Collectively, the four QTL explained 24% of the phenotypic variation for R8. Twelve QTL were detected for plant height, one was significant across sets, four each within Sets 1 and 3, and three within Set 2 (Table 5). Each QTL explained 2 to 21% of the phenotypic variation for plant height and accounted for 20- to 60-mm differences between the BSR 101 and LG82-8379 alleles at these loci. The plant height QTL, including the QTL that was significant across sets, collectively accounted for 30% of the phenotypic variation within Set 1, 42% within Set 2, and 55% within Set 3. One QTL was detected for lodging within Set 3 on linkage group D1b. This QTL explained 6% of the phenotypic variation (Table 5) and the LG82-8379 allele at this locus increased lodging slightly. No QTL for R1 or reproductive period were identified.
|
Three QTL were significantly (P < 0.05, LOD > 2.5) associated with oil concentration (Table 5). Two QTL, located on linkage group C1 and D1b, were significant in combined analysis of all lines across the three sets and explained 5 to 10% of the phenotypic variation for oil concentration. Collectively, the two QTL accounted for 10% of the phenotypic variation for oil concentration. The third QTL located on linkage group J was detected only within Set 1 and explained 16% of the phenotypic variation for oil concentration. Including the QTL that were significant across sets, the oil QTL collectively explained 26% of the phenotypic variation for oil concentration within Set 1. The difference in oil concentration between the BSR 101 and LG82-8379 alleles at these loci ranged from 2 to 3 g kg1.
| DISCUSSION |
|---|
|
|
|---|
To determine the impact of the significant environment by set and environment by lines within set interactions had on our QTL analysis for yield, we examined the effects of the 15 yield QTL in their respective sets in each of the 12 environments. The effect of the yield QTL was positive and statistically significant in 21% of the 180 (i.e., 15 x 12) yield QTL and environment combinations. In 60% of the yield QTL and environment combinations, the QTL effect was positive but not statistically significant. In only five of the yield QTL and environment combinations were the effects on yield negative, and none of these negative effects were statistically significant. Three of the five negative effects were associated with Satt142 on linkage group H. Although this QTL seems to be highly variable, the mean increase in yield among Set 3 lines associated with this QTL was the largest recorded for any of the QTL identified in this research.
Comparison of our results with those of others (Orf et al., 1999; Specht et al., 2001; Yuan et al., 2002) revealed QTL mapping to similar positions. Specht et al. (2001) identified a QTL, marked by Satt277, near our QTL on linkage group C2 that accounted for phenotypic variances of 13 to 15% for maturity, 5 to 15% for plant height, 6% for lodging, and 7 to 13% for seed yield in a RIL population derived from crossing Minsoy and Noir 1. In their study, the impact on seed yield from the beneficial Minsoy allele was attributed to the maturity locus E1 that is located 4 cM below Satt277 on linkage group C2 (Cregan et al., 1999). Orf et al. (1999) also discovered a QTL on linkage group C2 in a Noir 1 x Archer RIL population where the Noir 1 allele, marked by Satt277, increased seed yield 172 kg ha1 and accounted for 11% of the phenotypic variance. In this case, Satt277 was not significantly associated with maturity. In our population, the LG82-8379 allele for Satt363 on linkage group C2 was associated with increased seed yield of 60 kg ha1 and accounted for 10% of the phenotypic variance for yield across all sets. The LG82-8379 alleles in this region also were associated with increased plant height of 60 mm among Set 2 lines (R2 = 20%) and 40 mm among Set 3 lines (R2 = 21%) but were not associated with lodging. While Satt363 maps 15.2 cM from the maturity locus E1 on linkage group C2 (Cregan et al., 1999), no maturity differences were associated with Satt363 in our study.
Specht et al. (2001) identified a QTL, marked by Satt317, near our QTL on linkage group H, that accounted for phenotypic variances of 2% for plant height, 5% for lodging, and 3% for seed yield in a Minsoy x Noir 1 RIL population. The QTL association in their population was attributed to the Noir 1 semisparse pubescence (Ps-s) allele located near Satt317. Lines homozygous for the Noir 1 Ps-s allele tended to be spindly, lodging-prone, and generally yielded less than lines homozygous for the Minsoy ps allele. In Set 3 lines of our population, the LG82-8379 allele on linkage group H, marked by Satt142, increased seed yield 148 kg ha1 (R2 = 35%) and decreased plant height 20 mm (R2 = 4%). This genetic region was not associated with lodging and there was no obvious segregation for pubescence density in our population. On the University of Utah genetic linkage map (Cregan et al., 1999), Satt317 is 4.2 cM from Satt142 on linkage group H.
Yuan et al. (2002) identified yield QTL on linkage group K in two RIL populations. Essex provided the beneficial allele (R2 = 10%, 110 kg ha1 yield increase) marked by Satt337 in their Essex x Forrest RIL population, whereas Flyer provided the beneficial allele (R2 = 15%, 220 kg ha1 yield increase) marked by Satt326 in their Flyer x Hartwig RIL population. A QTL on linkage group K, marked by Satt326, also was identified by Specht et al. (2001) in a Minsoy x Noir 1 RIL population and it accounted for phenotypic variances of 2% for maturity and 2% for seed yield in only 1 yr of a 2 yr study. In our Set 1 lines, Satt544 marked the LG82-8379 allele on linkage group K that was associated with increasing seed yield 114 kg ha1 (R2 = 38%), increasing days to maturity by 1 d (R2 = 7%), and increasing plant height by 30 mm (R2 = 6%). Markers Satt326 and Satt337 map approximately 11.3 cM from Satt544 on linkage group K of the USDA/Iowa State University genetic linkage map (Cregan et al., 1999).
Specht et al. (2001) reported that 4% of the phenotypic variance for seed yield within a Minsoy x Noir 1 RIL population was accounted for by segregation at the Rpg4 locus on linkage group N for resistance/susceptibility to bacterial blight. Bacterial blight occurred within the test the year this association was made and lines homozygous for the susceptible Noir 1 allele yielded 181 kg ha1 less. In our population, the allele from BSR 101 on linkage group N, marked by Satt387 (R2 = 15%) increased seed yield 128 kg ha1 among Set 1 lines while Satt339 (R2 = 12%) increased seed yield 67 kg ha1 among Set 2 lines. The BSR 101 allele at Satt387 also was associated with increasing days to maturity by 2 d (R2 = 11%) and increasing plant height 30 mm (R2 = 5%). At Satt339, the BSR 101 allele also was associated with increasing plant height 20 mm (R2 = 2%). On the University of Utah genetic linkage map (Cregan et al., 1999), the Rpg4 locus on linkage group N is 7.5 cM from Satt387 and 23 cM from Satt339 that marks the QTL in our population.
Several yield QTL identified in this population also were associated with other agronomic traits. For example, the marker alleles associated with an increase in seed yield on linkage groups A1, G, K, and N also were associated with an increase in days to maturity and with taller plants. In contrast, marker alleles that were associated with greater yield on linkage groups H and D1a were associated with shorter plants. Although maturity and plant height cosegregated with several of the yield QTL, the variation observed was small and it is uncertain whether these differences would account for the change in seed yield associated with these regions. The correlation coefficients between seed yield and these traits were generally small and even though both plant height and maturity were significantly different between sets, seed yield was not. Further genetic dissection of the regions containing these QTL would be needed to distinguish between pleiotropy or gene linkage to define these relationships.
Identifying QTL with positive effects on seed yield while simultaneously improving or at least not decreasing protein concentration would be beneficial. Eight seed yield QTL and four protein concentration QTL were found to segregate independently of each other; however, three of the 15 yield QTL identified also were associated with increasing protein concentration. The yield-increasing LG82-8379 alleles, located on linkage groups B2, M, and O, also were associated with increased protein concentration of 2 to 3 g kg1. Unique seed yield and protein QTL may be segregating within this population.
Despite the observation of significant variation among lines within each set of this population and the high heritability estimates calculated for each trait measured (Table 1), many QTL identified were found only within lines of specific sets. It is not clear whether this is the result of a particular QTL identified only having an effect in certain maturities or whether there was insufficient resolution to detect a QTL within a particular set. The number of lines evaluated within each set was low, which can limit precision of QTL localization, cause an underestimation of the number of QTL controlling the trait, overestimate the effects of QTL, and reduce the power of QTL detection (Beavis, 1998). Confirmation studies will determine if the QTL identified in this study, within and among sets, are truly real or not. Three set-specific yield QTL that we identified in our population were found by other researchers to be significant within lines of different maturities in their populations (Specht et al., 2001; Yuan et al., 2002). This indicates that these particular yield QTL may not be specific to maturity groups or may be specific to maturity groups within certain genetic backgrounds. We have also determined that the early MG IV experimental line, LG91-7350R, previously developed from BSR 101 x LG82-8379 and used as a check in this study, carries the yield enhancing, Set 2 QTL marked by Satt168 on linkage group B2 from FC 04007B. Further examination of the environmental or developmental conditions the set-specific QTL respond to will be necessary.
Results from this study indicate that soybean PIs contain QTL alleles that may be of benefit to North American soybean cultivars. Of particular interest is the yield QTL marked by Satt358 on linkage group O. The QTL was significant in the combined analysis across sets of this population. For this QTL, the FC 04007B allele enhanced yield 47 kg ha1 and increased seed protein concentration 3 g kg1. To date, no other studies have detected a yield QTL at this chromosomal location. Confirmation of the yield-enhancing QTL from LG82-8379 is currently being performed in the genetic background of the population used in the current study and in different genetic backgrounds.
| ACKNOWLEDGMENTS |
|---|
Received for publication December 13, 2002.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
D. Li, T. W. Pfeiffer, and P. L. Cornelius Soybean QTL for Yield and Yield Components Associated with Glycine soja Alleles Crop Sci., March 19, 2008; 48(2): 571 - 581. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. S. Guzman, B. W. Diers, D. J. Neece, S. K. St. Martin, A. R. LeRoy, C. R. Grau, T. J. Hughes, and R. L. Nelson QTL Associated with Yield in Three Backcross-Derived Populations of Soybean Crop Sci., January 22, 2007; 47(1): 111 - 122. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. R. Walker, A. M. Scaboo, V. R. Pantalone, J. R. Wilcox, and H. R. Boerma Genetic Mapping of Loci Associated with Seed Phytic Acid Content in CX1834-1-2 Soybean Crop Sci., January 24, 2006; 46(1): 390 - 397. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. R. Panthee, V. R. Pantalone, D. R. West, A. M. Saxton, and C. E. Sams Quantitative Trait Loci for Seed Protein and Oil Concentration, and Seed Size in Soybean Crop Sci., August 26, 2005; 45(5): 2015 - 2022. [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 | |||