Published in Crop Sci. 43:2018-2027 (2003).
© 2003 Crop Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA
CROP BREEDING, GENETICS & CYTOLOGY
Genetic Components of Yield Stability in Maize Breeding Populations
E. A. Lee*,
T. K. Doerksen and
L. W. Kannenberg
University of Guelph, Department of Plant Agriculture, Crop Science Building, Guelph, ON, N1G 2W1 Canada
* Corresponding author (lizlee{at}uoguelph.ca).
 |
ABSTRACT
|
|---|
Phenotypic stability has long been recognized as an important target in plant breeding. Stability is influenced in part by the genetic structure, i.e., level of heterogeneity and heterozygosity, of the cultivar. Yet, very little is known about the genetic components underlying stability, and how population improvement strategies influence stability. We examined 12 maize (Zea mays L.) breeding populations selected via reciprocal recurrent selection (RRS), selfed progeny recurrent selection (S), or a method combining RRS and S (COM), to examine changes in the genetic structure of the phenotypic stability of three traits (grain yield, grain moisture, and broken stalks), and two associated selection indices. Partitioning of the genotype x environment sums of squares from diallel matings of the original (C0) and advanced (CA) cycle populations into linear trends indicated that only grain yield and the unadjusted performance index (UPI) followed a predictable linear response. Grain yield and UPI linear trends were further partitioned by Gardner and Eberhart Analysis III to examine the genetic components of stability. We found that recurrent selection (RS) improved grain yield stability, and that this trait is heritable, predictable, and mostly controlled through additive gene action. Improvement in grain yield stability was observed both in cross and per se performance and was accompanied by significant improvement in the mean performance of the populations. However, the improvement in grain yield stability did not result in substantial changes in the general combining ability (gi) estimates of most populations. Our results indicate that grain yield stability can be improved through RS by selecting solely for mean performance across multiple environments.
Abbreviations: API, adjusted performance index C0, original cycle CA, advanced cycle COM, combined recurrent selection GCA, general combining ability RS, recurrent selection RRS, reciprocal recurrent selection S, selfed-progeny recurrent selection SCA, specific combining ability UPI, unadjusted performance index P, phenotype G, genotype E, environment ANOVA, analysis of variance CG, Canada-Guelph CCGP, Cross Canada Gene Pool RCBD, randomized complete block design LR, linear regression AMMI, additive main effects and multiplicative interaction IPCA, interaction principal component axis
 |
INTRODUCTION
|
|---|
QUANTITATIVE GENETIC THEORY states that an individual's phenotype (P) is the product of the genetics (G) of the individual, the environment (E) that the individual is exposed to, and the interaction that occurs between the genotype of the individual and the environment (G x E). Large G x E effects tend to be viewed as problematic in breeding because the lack of a predictable response hinders progress from selection (Dudley and Moll, 1969). This idealized predictable response across multiple environments is generally referred to as stability (Cannon, 1932). When stability concepts are applied to breedingselection for increased grain yield, the idealized genotype is one that is capable of utilizing the resources available in higher yielding environments and has a mean performance that is above average in all environments.
Yield stability is influenced in part by the genetic structure of the variety. More heterozygous varieties and more heterogeneous varieties are less affected by environmental differences (Lerner, 1954; Lewontin, 1957; Allard and Bradshaw, 1964). For example, maize double-cross hybrids have smaller G x E interactions and are more stable than maize single-cross hybrids (Sprague and Federer, 1951). More homogeneous generations, inbred and F1, have larger G x E interactions than more heterogeneous generations, F2 and BC1, and a more homozygous generation, inbred, has a larger G x E interaction than a more heterozygous generation, F1 (Valdivia-Bernal and Hallauer, 1991). To reduce the complications that the G x E interaction creates when selecting superior genotypes, many attempts have been made to (i) understand the environmental components causing the G x E interaction (Epinat-Le Signor et al., 2001; Jeutong et al., 2000), (ii) examine the G x E interaction biometrically (Lin et al., 1986; Zobel et al. 1988; Finlay and Wilkinson, 1963; Yan et al., 2001), and (iii) develop selection strategies that involve a stability parameter (Magari and Kang, 1993). However, very little effort has been focused on examining the genetic components underlying, in particular, yield stability. In maize, the limited information regarding the genetics of yield stability comes from diallel studies involving inbred lines. These results indicate that stability (both positive and negative) is heritable and controlled by additive gene action (Eberhart and Russell, 1966, 1969).
One approach to examining stability is to further partition the G x E interaction from a traditional analysis of variance (ANOVA) into linear trends and a departure from linear (residual) (Finlay and Wilkinson, 1963; Eberhart and Russell, 1966). Linear trends can then be further partitioned by Analysis III of Gardner and Eberhart (1966) to arrive at estimates of genetic effects. It is then possible to examine the genetic structure underlying yield stability in terms of heterosis, general combining ability (GCA), and specific combining ability (SCA). The linear portion of the G x E interaction presumably has an additive and a nonadditive genetic component. Understanding the relative importance and genetic nature of these components and how selection influences them will greatly enhance our knowledge of yield stability and will help determine if selection strategies designed to enhance stability are important in breeding programs. In this paper, we present the results from two diallels involving, respectively, the originating gene pool (C0) and the most advanced cycle (CA) of 12 breeding populations. Specifically, we were interested in examining (i) if genetic progress for yield stability has been achieved indirectly through RS, (ii) if RS produced changes in the genetic structure of the populations for yield stability, and (iii) the underlying genetic components that influence yield stability.
 |
MATERIALS AND METHODS
|
|---|
Genetic Materials
The 11 C0 and 12 CA (C2C6) early (25502900 Ontario Crop Heat Units [OCHU] [Brown and Bootsma, 1993]) University of Guelph breeding populations in this study included CG-SynA (S), CG-SynB (S), CG-Wigor (RRS) and CG-SynA (RRS), CG-HOPE Elite A (RRS) and CG-HOPE Elite B (RRS), CG-Stiff Stalk (COM) and CG-Lancaster (COM), CG-CBI (RRS) and CG-CBII (RRS), and CG-Cross Canada Gene Pool A (RRS) and CG-CCGP B (RRS). The two CCGP populations were selected from a common C0, resulting in 11 C0 populations and 12 CA populations. All other populations were unrelated except that CG-SynA (RRS) and CG-SynA (S) originated from different cycles of selection in the CG SynA gene pool. Four to six cycles of RS had been completed in each population except for the CCGP populations, which had undergone two cycles of RS. For each cycle, families within a population were selected visually during development and subsequently in replicated trials in two environments, i.e., locations. Further information regarding the populations can be found elsewhere (Lee, 2002). Unadjusted (UPI) and adjusted (API) performance indices were used as the selection criteria for advancement until the populations had attained a level of maturity comparable to adapted germplasm for 2600 OCHU. The UPI was the ratio of grain yield to the percentage of grain moisture at harvest, while the API was the same ratio, but based on standing plants only, i.e., whole-plot grain yield x the percentage of plants neither lodged nor with broken stalks. Progenies were ranked by each index and selected based on the sum of ranks to identify those progenies for intercrossing to form the next cycle. When satisfactory grain moisture at harvest was reached, selection was based only on the sum of ranks for grain yield and grain yield of standing plants. In general, 20 selected families were intercrossed to form a new cycle.
In 1999, separate half-diallel matings, not including reciprocals, were made with the C0 populations (n[n - 1]/2 = 55 crosses) and CA (n[n - 1]/2 = 66 crosses). For each cross, six ears were pollinated with bulked pollen from six plants from the other population. Reciprocal crosses were then bulked, resulting in a total of 12 ears per cross. Individual plants were used as either males or females in the crosses, but not both, for a total of 24 S0 plants, representing 24 gametes from each population. This approach samples predominant alleles and linkage blocks; rare alleles and linkage blocks are likely not represented in the population crosses and is similar to other studies utilizing population crosses (Moll et al. 1977; de Leon and Lonnquist 1978; Douglas et al. 1961; Gouesnard et al. 1991). Seed was increased for each of the 11 C0 and 12 CA parental populations via chain sibbing and at least 70 individual ears were harvested from each population. Equal amounts of seed from each ear of the populations per se were bulked for use in the experimental trials.
Experimental Procedure and Data Collection
The experiment included 144 entries, 66 crosses from the CA diallel mating, 55 crosses from the C0 diallel mating, 12 advanced populations per se, and 11 C0 populations per se. Yield trials were grown in both 2000 and 2001 at three southwestern Ontario locations (Alma [2500 OCHUs], Elora [2600 OCHUs], and Woodstock [2850 OCHUs]) in a 12 x 12 partially balanced lattice design with three replications at each location. The soil type at all locations is Guelph loam (Typic Hapludalf). Fertilizer (N, P2O5, and K2O, respectively) was applied on the basis of soil tests at the rate of Alma (2000) 155, 58, and 74 kg ha-1; Alma (2001) 122, 75, and 57 kg ha-1 supplemented with 33 690 L ha-1 liquid swine manure; Elora (2000) 140, 50, and 50 kg ha-1 and 135, 32, and 76 kg ha-1; Elora (2001); Woodstock (2000) 150, 66, and 108 kg ha-1; and Woodstock (2001) 167, 59, and 75 kg ha-1. Weeds were controlled by conventional herbicides. Experimental units were two-row plots, 5.78 m long, with a spacing of 0.76 m between rows. Trials were overplanted and thinned at the six-leaf stage to uniform stands of 68 000 plants ha-1 (60 plants plot-1). All trials were machine planted (New Idea four-row planter New Idea Company, Cold Water, OH) in May with Almaco planting cones and machine harvested (Almaco SPC40-2 two row combine [ALMACO, Allan Machine Company, Nevada, IA)] in October or November. Three traits and two indices were measured, respectively: machine-harvested grain yield (kg ha-1) adjusted to 155 g kg-1 grain moisture grain moisture at harvest, percentage of broken stalks (plants broken below the ear or inclined more than 45° from the vertical), UPI, and API. All five parameters were analyzed and will be referred to as traits.
Statistical Analysis
Data were combined across locations and analyzed as a partially balanced lattice by means of PROC MIXED and a randomized complete block design (RCBD) by means of PROC GLM (SAS, 1996). The efficiency of the lattice design was tested by determining the ratio of the lattice MSerror x 100/RCBD MSerror for each analyzed trait. Analysis of variance indicated only slight gains in efficiency when a lattice was used rather than a RCBD, and thus the RCBD analysis and means were used. Data across locations were pooled and analyzed as a randomized complete block design (RCBD) by PROC GLM for the plot-based linear regression (LR) and PROC GLM and PROC PRINCOMP for the mean based AMMI models using SAS (1996). The LR on environmental mean index model is
and the additive main effects and multiplicative interaction model (AMMI) is
where Ygre is the yield of genotype g in replication r at environment e; µ is the grand mean;
g and ße are the genotype and environment deviations from the grand mean;
r(ße) is the replicate effect nested within an environment;
g is the yield sensitivity (regression coefficient) of genotype g to the change in the environmental mean;
n is the singular value for interaction principal component axis (IPCA) n;
gn and
en are the genotype and environment eigenvectors for axis n, respectively;
ge is the deviation from linear trends (LR model) or the residual not accounted for by retained IPCAs (AMMI model) for the G x E interaction, respectively; and
rge (or
ge) is the random experimental error. Both genotype (
n0.5
gn) and environment (
n0.5
en) IPCA scores are scaled so that yield is transformed to the square root of yield. When IPCAs are graphed together, they have the units of the square root of yield and the product of the scaled genotype and environment IPCAs, thus giving the value directly in terms of the model (Gauch, 1992). The
gße component (linear trends) in the LR model was further partitioned to the linear response due to GCA and SCA combining ability effects (Eberhart and Russell, 1966, 1969), according to Analysis III (Gardner and Eberhart, 1966). Genotypes in the LR model were considered stable when
g = 1. All effects partitioned out of the G x E interaction in the LR model were tested for significance by
gre mean squares (Finlay and Wilkinson, 1963), not pooled deviations mean squares (Eberhart and Russell, 1966). Thus, standard errors of the regression parameter estimates are from the
gre mean squares. In addition, an F test of the ratio of GCA to SCA mean squares was used to test the significance of additive to nonadditive variation or error, if the SCA is assumed randomly and normally distributed with constant variance (Lin and Binns, 1991). The
ge in the AMMI model was taken from the plot-based LR model divided by the average number of replications as suggested by Cornelius (1993) to test the significance of IPCAs by means of the FGH2 test.
Least significant differences (LSDs) from zero and between estimates were calculated according to Griffing (1956) and Zhang and Kang (1997) for Analysis III. The calculation of LSDs for gi estimates from zero, and between estimates were, respectively
where p was the number of parents in the diallel cross, n was the error degrees freedom, e was the number of environments, r was the number of replications, and MSerror was the respective source x environment mean square. Rank correlations were calculated between C0 and CA gi estimates. All hypotheses were tested with a Type I error (
) of 0.05.
 |
RESULTS AND DISCUSSION
|
|---|
Analysis III
Environments, genotypes, and G x E were significant sources of variation for all traits (Doerksen et al., 2003). Genotype x environment was further partitioned into linear trends and deviations from linear. The linear trend was not significant for grain moisture, broken stalks, or API. Grain yield and UPI were the only traits where both linear trends and deviations from linear trends across environments were significant (Table 1) . However, as a percentage of the G x E sums of squares, linear trends for grain yield and UPI accounted for only 25 and 26% of the variance, respectively. Linear trends were further partitioned into CA vs. C0, E x C0linear, and E x CAlinear. The CA vs. C0 contrast was significant for both grain yield and UPI, indicating significant differences in phenotypic stability between C0 and CA. Furthermore, there were significant differences among regression coefficients within the E x C0linear and within the E x CAlinear, indicating that differences in phenotypic stability exist between populations.
View this table:
[in this window]
[in a new window]
|
Table 1. Mean squares from analyses of variance for grain yield and unadjusted performance index (UPI) of an original and advanced cycle maize population diallels, conducted over six location-years, Ontario, 2000-2001. Differences in linear trends were partitioned according to Gardner-Eberhart's (1966) Analysis III.
|
|
The E x CAlinear and E x C0linear variances for grain yield and UPI were further partitioned into parents, heterosis, and crosses by Analysis III (Gardner and Eberhart, 1966). In both CA and C0, differences in linear trends were largely due to crosses (88 and 83%, respectively). Surprisingly, no significant heterotic effect among regression coefficients was detected, implying that heterosis in this set of material is not conferring increased phenotypic stability. This is contrary to other studies that have suggested that heterosis is important for conferring stability (Valdivia-Bernal and Hallauer, 1991). Heterosis is important for grain yield in this set of materials. When the entry variance was partitioned in a similar manner, a significant heterosis effect was observed for grain yield in both C0 and CA (Doerksen et al., 2003). But, perhaps there is sufficient heterogeneity and heterozygosity in C0, which is maintained in CA, that additional heterozygosity does not confer additional stability. Crosses were further partitioned into GCA and SCA for both E x CAlinear and E x C0linear. Differences in linear trends due to GCA were significant in both C0 and CA, while differences due to SCA were significant only in CA. F tests of the ratio of GCA to SCA mean squares were significant in C0 and CA for grain yield (3.03 and 2.57, respectively), but not for UPI (1.99 and 1.87, respectively), indicating that grain yield stability is mostly controlled by additive genetic effects while UPI stability is controlled both by additive and nonadditive effects. Yield stability is heritable and predictable as indicated by significant GCA effects in both the C0 and CA. Yield stability has been improved through recurrent selection indicated by the significant CA vs. C0 contrast. In both gene pools, i.e., C0 and CA, there appears to be sufficient heterozygosity and heterogeneity that heterosis does not confer any additional stability. And finally, mostly additive genetic affects are controlling grain yield stability.
Stability Analysis
Yield stability was examined initially by AMMI analysis on the C0 and CA separately and combined. The first two principal components accounted for 41 and 20% of the variation, respectively. However, large environmental scores and small genotype scores were obtained, resulting in a tight clustering of the genotypes at the center of the plot (data not shown). Because of the lack of informativeness of the AMMI results, the LR approaches of Finlay and Wilkinson (1963) and Eberhart and Russell (1966) were used to examine the mean response of a population both in crosses and per se. Linear regression, a Type 2 stability statistic (Lin et al., 1986), was used because we were only interested in comparisons among the set of genotypes, and the LR approach allowed us to exploit the genetic structure that exists among the entries.
Because of the additive nature (significant GCA and lack of heterosis) of grain yield stability observed with Analysis III, we used, as a measure of a population's stability, the mean bi values for a population averaged over all crosses involving that population (Table 2)
. As implied by Analysis III, all CA populations in crosses had significantly higher bi values than their corresponding C0 populations, showing that RS had improved the responsiveness of the populations to higher yielding environments. Accompanying this increase in responsiveness was also a significant increase in mean yield in crosses for all populations. Taken together, it appears that the responsiveness did not come at the expense of average performance. This trend is consistent across all populations and can be easily observed graphically (Fig. 111)
. The deviations
associated with the average bi values were not very informative in this set of material (Table 2). Also, when per se yield stability was examined, no pattern of improvement was apparent. For some populations, no significant change in bi values was observed, while other populations did exhibit an improvement in responsiveness. Somewhat surprisingly, the SynB(S) population actually showed a loss of responsiveness from C0 to CA (1.23 to 0.82). Even though the responsiveness of the populations per se was not always improved, all populations exhibited significant improvement in grain yield on a per se basis. In this set of 12 RS populations, RS has improved cross yield stability while simultaneously improving mean performance of the populations.
View this table:
[in this window]
[in a new window]
|
Table 2. Mean yield stability values (bi), deviations (s2d) and grain yield both in maize crosses and per se, and gi estimates for bi values and grain yield for C0 and CA populations.
|
|
Relationship between C0 and CA Genetic Effects for Yield Stability
In C0, the differences among linear trends for gi estimates were largely due to a large negative contribution from Wigor (-0.31) and a large positive contribution from SynA(S) (0.19) (Table 2). This indicated that Wigor and SynA(S) were contributing lower and higher grain yield to their crosses, respectively, in higher yielding environments. In CA, the differences among linear trends for gi estimates were largely due to large negative contributions from Wigor (-0.24) and SynA(S) (-0.19) (Table 2). The significant rank correlation (r = 0.59, P > 0.05) between C0 and CA gi estimates suggested that the stability parameter, while improved on through recurrent selection, was essentially a characteristic present in the original population.
Consequences of Recurrent Selection in the 12 Breeding Populations
We have now examined two important aspects of the 12 RS populations. First, we were interested in what changes in genetic structure had occurred that influenced grain yield (Doerksen et al., 2003). We found that genetic improvement occurred in both the per se and cross performance of most populations. Accompanying the favorable changes in population performance were less favorable shifts from predominantly additive genetic effects in C0 to greater nonadditive genetic effects in CA. This shift did not substantially change gi estimates of most populations. However, in the case of grain yield, the underlying components of gi effects were altered in their relative importance. General combining ability effects in C0 were caused primarily by population per se effects (vi), while in CA, the GCA effects were caused predominately by parental heterotic effects (hi) (Doerksen et al., 2003).
The other aspect of these populations, addressed in this paper, is phenotypic stability for grain yield. Yield stability is heritable and predictable. It has also been improved through recurrent selection. In both gene pools, i.e., C0 and CA, there appears to be sufficient heterozygosity and heterogeneity that heterosis (i.e., population crosses) does not confer any additional stability. Mostly additive genetic affects are controlling grain yield stability, with no apparent shift toward more nonadditive effects in CA. Recurrent selection for only grain yield across multiple environments has resulted in improved cross yield stability while simultaneously improving mean performance of the populations. It is important to note that the number of test environments during RS was only two per cycle, i.e., four cycles = eight performance tests, six cycles = 12 performance tests. Apparently, the two environments per cycle, plus the different years of testing each cycle, even though a relatively small number (eight to 12 for most populations), has been sufficient to select (indirectly) for stability of grain yield. Or, interpreted another way, the reduction of deleterious gene frequencies for grain yield, which would be especially effective in the early cycles of RS, develops a more refined heterogeneous, heterozygous population, and therefore a more stable one. In effect, a good yield gene is also a good stability gene.
We also can draw some preliminary conclusions regarding the effectiveness of the various methods of improvement, S, RRS, and COM. One of the most surprising observations from the diallels is that in many cases in CA, the crosses involving the RRS pairs were not superior to crosses involving other populations (Table 3) . And in some instances, Wigor-SynA(R) and Stiff Stalk-Lancaster, most of the crosses involving other populations were superior to the RRS pair cross. This is contrary to RS theory. Given this set of data, it appears that RRS will not result in population pairs that are consistently superior to populations improved independently from one another.
View this table:
[in this window]
[in a new window]
|
Table 3. Mean grain yield for C0 and CA maize population crosses. In each grouping, the cross involving the RRS pair is italicized. LSD (0.05) = 608 kg ha-1.
|
|
 |
NOTES
|
|---|
Part of a thesis submitted by T.K. Doerksen in partial fulfillment of the requirements for the M.Sc. degree. Financial support from the Ontario Ministry of Agriculture and Food, and the National Sciences and Engineering Research Council of Canada is acknowledged.
Received for publication May 27, 2003.
 |
REFERENCES
|
|---|
- Allard, R.W., and A.D. Bradshaw. 1964. Implications of genotype-environmental interactions in applied plant breeding. Crop Sci. 4:503508.[Free Full Text]
- Brown, D.M., and A. Bootsma. 1993. Crop heat units for corn and other warm season crops in Southern Ontario. Ont. Minist. Agric. Food Factsheet, Agdex 111/31.Ontario Ministry of Agriculture and Food, Toronto, ON.
- Cannon, W.B. 1932. The wisdom of the body. Norton, New York.
- Cornelius, P.L. 1993. Statistical tests and retention of terms in the additive main effects and multiplicative model for cultivar trials. Crop Sci. 33:11861193.[Abstract/Free Full Text]
- de Leon, J.L., and J.H. Lonnquist. 1978. Heterosis in full-sibs within and between half-sib families in open-pollinated varieties of maize. Crop Sci. 18:2628.[Abstract/Free Full Text]
- Doerksen, T.K., L.W. Kannenberg, and E.A. Lee. 2003. Effect of recurrent selection on combining ability in maize breeding populations. Crop. Sci. 43:16521658.[Abstract/Free Full Text]
- Douglas, A.G., J.W. Collier, M.F. El-Ebrashy, and J.S. Rogers. 1961. An evaluation of three cycles of reciprocal recurrent selection in a corn improvement program. Crop Sci. 3:157161.
- Dudley, J.W., and R.H. Moll. 1969. Interpretation and use of estimates of heritability and genetic variances in plant breeding. Crop Sci. 9:257262.[Abstract/Free Full Text]
- Eberhart, S.A., and W.A. Russell. 1966. Stability parameters for comparing varieties. Crop Sci. 6:3640.[Abstract/Free Full Text]
- Eberhart, S.A., and W.A. Russell. 1969. Yield and stability for a 10-line diallel of single-cross and double-cross maize hybrids. Crop Sci. 9:357361.[Abstract/Free Full Text]
- Epinat-Le Signor, C., S. Dousse, J. Lorgeou, J.-B. Denis, R. Bonhomme, P. Carolo, and A. Charcosset. 2001. Interpretation of genotype x environment interactions for early maize hybrids over 12 years. Crop Sci. 41:663669.[Abstract/Free Full Text]
- Finlay, K.W., and G.N. Wilkinson. 1963. The analysis of adaptation in a plant breeding programme. Aust. J. Agric. Res. 14:742754.[ISI]
- Gardner, C.O., and S.A. Eberhart. 1966. Analysis and interpretation of the variety cross diallel and related populations. Biometrics 22:439452.[ISI][Medline]
- Gauch, H.G. 1992. Statistical analysis of regional yield trials: AMMI analysis of factorial designs. Elsevier, Amsterdam.
- Gouesnard, B., A. Gallais, M. Lefort-Buson, and J.P. Sampoux. 1991. Comparison of direct and correlated expected responses to intrapopluation and interpopulation selections for two forage maize populations. Maydica 36:4956.
- Griffing, B. 1956. Concept of general and specific combining ability in relation to diallel crossing systems. Aust. J. Biol. Sci. 9:463493.
- Jeutong, F., K.M. Eskridge, W.J. Waltman, and O.S. Smith. 2000. Comparison of bioclimatic indices for prediction of maize yields. Crop Sci. 40:16121617.[Abstract/Free Full Text]
- Lee, E.A. 2002. Corn breeding and genetics germplasm page. Available at http://www.plant.uoguelph.ca/research/corn_breeding/Populations2.htm (verified 4 June 2003).
- Lerner, I.M. 1954. Genetic homeostasis. Oliver and Boyd, London.
- Lewontin, R.C. 1957. The adaptation of populations to varying environments. Cold Spring Harb. Symp. Quant. Biol. 22:395408.[Medline]
- Lin, C.S., and M.R. Binns. 1991. Genetic properties of four types of stability parameters. Theor. Appl. Genet. 82:505509.[ISI]
- Lin, C.S., M.R. Binns, and L.P. Lefkovitch. 1986. Stability analysis: Where do we stand? Crop Sci. 26:894900.[Abstract/Free Full Text]
- Magari, R., and M.S. Kang. 1993. Genotype selection via a new yield-stability statistic in maize yield trials. Euphytica 70:105111.
- Moll, R.H., A. Bari, and C.W. Stuber. 1977. Frequency distribution of maize yield before and after reciprocal recurrent selection. Crop Sci. 17:794796.[Abstract/Free Full Text]
- SAS Institute. 1996. The SAS system for windows. Version 6.12. SAS Inst., Cary, NC.
- Sprague, G.F., and W.T. Federer. 1951. A comparison of variance components in corn yield trials. II. Error, year x variety, location x variety components. Agron. J. 43:535541.[Free Full Text]
- Valdivia-Bernal, R., and A.R. Hallauer. 1991. Estimates of genetic homeostasis in maize. Bras. J. Genet. 14:483499.
- Yan, W., P.L. Cornelius, J. Crossa, and L.A. Hunt. 2001. Two types of GGE biplots for analyzing multi-environment trial data. Crop Sci. 41:656663.[Abstract/Free Full Text]
- Zhang, Y., and M.S. Kang. 1997. DIALLEL-SAS: A SAS program for Griffing's diallel analyses. Agron. J. 89:176182.[Abstract/Free Full Text]
- Zobel, R.W., J.J. Wright, and H.G. Gauch. 1988. Statistical analysis of yield trial. Agron. J. 80:388393.[Abstract/Free Full Text]
This article has been cited by other articles:

|
 |

|
 |
 
F. J. Kutka and M. E. Smith
How Many Parents Give the Highest Yield in Predicted Synthetic and Composite Populations of Maize?
Crop Sci.,
September 1, 2007;
47(5):
1905 - 1913.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
E.A. Lee, R. Chakravarty, B. Good, M.J. Ash, and L.W. Kannenberg
Registration of 38 Maize (Zea mays L.) Breeding Populations Adapted to Short-Season Environments
Crop Sci.,
November 21, 2006;
46(6):
2728 - 2733.
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
Y. Zhang, M. S. Kang, and K. R. Lamkey
DIALLEL-SAS05: A Comprehensive Program for Griffing's and Gardner-Eberhart Analyses
Agron. J.,
June 17, 2005;
97(4):
1097 - 1106.
[Abstract]
[Full Text]
[PDF]
|
 |
|