Crop Science Grow Your Career with CSSA
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (3)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Yu, J.
Right arrow Articles by Bernardo, R.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Yu, J.
Right arrow Articles by Bernardo, R.
Agricola
Right arrow Articles by Yu, J.
Right arrow Articles by Bernardo, R.
Related Collections
Right arrow Crop Genetics
Right arrow Maize
Published in Crop Sci. 44:405-410 (2004).
© 2004 Crop Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA

CROP BREEDING, GENETICS & CYTOLOGY

Changes in Genetic Variance during Advanced Cycle Breeding in Maize

Jianming Yu and Rex Bernardo*

Dep. of Agronomy and Plant Genetics, Univ. of Minnesota, 411 Borlaug Hall, 1991 Upper Buford Circle, St. Paul, MN 55108-6026

* Corresponding author (berna022{at}umn.edu).


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Advanced cycle breeding in crops aims to develop newer, improved versions of inbreds from crosses among current, elite inbreds. Although genetic gains in many crops have been achieved through advanced cycle breeding, this breeding process systematically leads to a narrow germplasm base. Our objectives were to estimate the loss, if any, in genetic variance during advanced cycle breeding in maize (Zea mays L.), and to compare these genetic variance estimates with those expected under a testcross additive model. We evaluated the Iowa Stiff Stalk Synthetic (BSSS), B73 x B37, and Lo904 x Lo916 populations, which have a common genetic background but have different levels of genetic variance expected from the testcross additive model; that is, 100% in BSSS, 82% in B73 x B37, and 45% in Lo904 x Lo916. Contrary to expectations, testcross additive variances (VTA) for grain yield did not significantly differ among the three populations. As expected, VTA for grain moisture, stalk lodging, plant height, and ear height were significantly lower in B73 x B37 than in BSSS. The decreases in VTA for these traits were not significant from B73 x B37 to Lo904 x Lo916. Fitting a testcross epistatic model did not result in positive estimates of both VTA and testcross additive-by-additive variance (VTAA) for any trait. In summary, genetic variance, except for grain yield, was lost during advanced cycle breeding and the amount of loss was greater than expected.

Abbreviations: BSSS, Iowa Stiff Stalk Synthetic • VTA, testcross additive variance • VTAA, testcross additive-by-additive variance • VTC, testcross variance


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
IN ADVANCED cycle breeding, two elite inbreds are crossed to form a base population from which new inbreds are developed (Allard, 1960, p. 115; Bernardo, 2002, p. 70). Advanced cycle breeding has become the most common method for developing inbreds in major crops such as barley (Hordeum vulgare L.; Rasmusson and Phillips, 1997), maize (Hallauer, 1990), rice (Oryza sativa L.; Peng et al., 1999), soybean [Glycine max (L.) Merr.; Cornelious and Sneller, 2002], and wheat (Triticum aestivum L.; Knott, 1987). In maize, for example, open-pollinated cultivars were the base populations used to develop nearly 100% of the inbreds (i.e., first-cycle inbreds) in the USA in the 1930s (Jenkins, 1978). The frequency of second-cycle inbreds, developed from crosses among the first-cycle inbreds, reached 50% by 1960 (Jenkins, 1978). Advanced-cycle breeding in maize is currently in about the fifth cycle (J. Geadelmann, 2001, personal communication).

The average yield in different crops continues to increase by about 1 to 2% each year and these linear increases across time demonstrate that advanced cycle breeding has not reached a yield plateau (Poehlman and Sleper, 1995, p. 4). Advanced cycle breeding, however, systematically reduces genetic variation (Rasmusson and Phillips, 1997; Tanksley and McCouch, 1997). From a population genetics standpoint, crossing two inbreds to form a breeding population at each cycle creates a bottleneck, which is defined as a severe reduction in the number of parents that are mated to form the next generation. Assuming no dominance and no epistasis, the additive genetic variance in a bottlenecked population decreases as a linear function of the inbreeding coefficient (Falconer and Mackay, 1996, p. 264). Bottlenecks would therefore lead to a continuous loss of genetic variation in breeding populations from cycle to cycle.

Concern about the loss of genetic variation in crops has increased (Tanksley and McCouch, 1997; Hoisington et al., 1999). No research to systematically monitor the changes in genetic variance during advanced cycle breeding has been reported. Our objectives were to estimate the loss, if any, in genetic variance during advanced cycle breeding in maize, and to compare these genetic variance estimates with those expected under a testcross additive model.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Plant Populations
The maize population BSSS was created in the 1930s by intermating 16 inbreds of above-average stalk quality (Sprague and Jenkins, 1943). Inbred B37 was developed in 1958 directly from BSSS (Henderson, 1976), and inbred B73 was developed in 1972 from the fifth generation of recurrent selection in BSSS (Russell, 1972). Inbred Lo904 was developed in 1983 from a (B73 x B37)B73 population, and inbred Lo916 was developed in 1984 from a B73 x B37 population (Bertolini et al., 1991). BSSS, B73 x B37, and Lo904 x Lo916 represent three maize breeding populations at different breeding cycles and have a common genetic background.

Seeds of BSSS were obtained from Iowa State University. Seeds of B73 and B37 were obtained from stocks maintained at Purdue University. Seeds of Lo904 and Lo916 were obtained from the Instituto Sperimentale per la Cerealicoltura, Bergamo, Italy. Six populations (BSSS, B73 x B37, and Lo904 x Lo916, each at two selfing generations, S0 and S2) were synthesized. The BSSS-S0 was formed by random mating 180 S0 plants, and the B73 x B37-S0 and Lo904 x Lo916-S0 populations by random mating 180 F2 plants. Random plants from the BSSS-S0, B73 x B37-S0, and Lo904 x Lo916-S0 populations were selfed for two generations to create S2 families. The resulting populations were denoted by BSSS-S2, B73 x B37-S2, and Lo904 x Lo916-S2. Having populations at two selfing generations allowed us to calculate a pooled variance estimate with the assumption of no epistasis, and to test an epistatic model by considering the two estimates separately.

Three non-BSSS inbreds (LH185, LH295, and a proprietary Syngenta inbred) were used as testers (Table 1). These three testers were chosen so that the resulting testcrosses could be evaluated in the maturity zone of the testing environments. An individual S0 plant was crossed as the male parent to two plants of a tester, and the resulting seeds were bulked as an S0 testcross. A single S2 plant from each S0 family was likewise crossed as the male parent to two plants of a tester to obtain an S2 testcross.


View this table:
[in this window]
[in a new window]
 
Table 1. Number of families for different maize testcross populations evaluated.

 
Field Evaluation
The testcross populations with LH185 as the tester (i.e., LH185 testcrosses) were evaluated at North Vernon, Romney, and West Lafayette, IN, in 2000. The testcross populations with the Syngenta inbred as the tester (i.e., Syngenta testcrosses) were evaluated at North Vernon, Romney, West Lafayette-1, and West Lafayette-2 (planted 2 wk later than West Lafayette-1), IN, in 2001. The testcross populations with LH295 as the tester (i.e., LH295 testcrosses) were evaluated at Lamberton-1, Lamberton-2 (planted 2 wk later than Lamberton-1), Rosemount, and Waseca, MN, in 2002. The testcrosses were randomly divided into different sets and each set was evaluated in a randomized complete block design (Table 1). A single replication was used for the LH185 and Syngenta testcrosses, and two replications were used for the LH295 testcrosses.

The plants were grown in two-row plots, each row 6.7 m long and spaced 0.76 m apart, at a plant population density of 68678 plants ha–1. Data recorded on plot were grain yield (Mg ha–1, adjusted to 155 g H2O kg–1), grain moisture (g kg–1), stalk lodging (percentage of stalk lodged plant per plot), plant height (cm, distance from the soil surface to the tip of the tassel), and ear height (cm, distance from the soil surface to the highest ear-bearing node). Grain yield and grain moisture were recorded for the LH185, Syngenta, and LH295 testcrosses. Stalk lodging was recorded for the LH185 and Syngenta testcrosses, and plant height and ear height were recorded for the LH295 testcrosses.

Expected Genetic Variance
A zero inbreeding coefficient was assumed for the BSSS population. The inbreeding coefficient of B73 x B37 (FB73xB37), corresponding to the coefficient of coancestry between B73 and B37 (fB73, B37), was calculated from the average of three studies (Godshalk et al., 1990; Messmer et al., 1991; Bernardo, 1993) as summarized by Bernardo (1993). The parental contribution of B73 and B37 to Lo904 and Lo916 was based on a study with 195 SSR markers (Bernardo et al., 2000). The inbreeding coefficient of Lo904 x Lo916 (FLo904xLo916), corresponding to the coefficient of coancestry of Lo904 and Lo916 (fLo904, Lo916), was further calculated as described by Bernardo et al. (2000).

Because testcross genotypic values are additive (Hallauer and Miranda Filho, 1988, p. 28; Bernardo, 2002, p. 79), the expected testcross variance can be calculated as

where VTCnew is the expected testcross variance of the bottlenecked population, VTCoriginal is the testcross variance of the original population, and Fnew (i.e., FB73xB37 or FLo904xLo916) is the inbreeding coefficient of the bottlenecked population with reference to the original population (Falconer and Mackay, 1996, p. 264).

Statistical Analysis and Observed Genetic Variance
The linear model for the LH185 and Syngenta testcrosses at each selfing generation was

whereas the linear model for the LH295 testcrosses was

where yijm (yijkm) is the observed value for testcross m from set j at environment i, (replication k); µ is the overall mean; li is the effect of environment i; sj is the effect of set j; slij is the effect of interaction between set j and environment i; rk(ij) is the effect of replication k; gm(j) is the effect of testcross m in set j; glim(j) is the effect of the interaction between testcross m and environment i; and {epsilon}ijm ({epsilon}ijkm) is the random error. The environments and testcrosses were considered as random effects. The mean squares from the ANOVA were equated to their expectations to solve for the testcross variance (VTC) (SAS Institute, 1990).

The VTC is a function of VTA and VTAA. We considered a testcross additive model and a testcross epistatic model. For the testcross additive model, VTA was solved from

where n is the number of selfing generations, that is, n = 0 for S0 and n = 2 for S2. A pooled estimate (across selfing generations) of VTA was obtained for each population. A pooled VTA estimate (across testers) was also calculated when a trait was evaluated for more than one tester. For the testcross epistatic model, VTA and VTAA were solved from

A 90% confidence interval on the pooled VTA was obtained based on the method of Ting et al. (1990) and Burdick and Graybill (1992)(p. 56). The variance estimates were regarded as significantly different from zero if their confidence intervals did not include zero. A significant difference between two estimates was declared if their confidence intervals did not overlap.


    RESULTS AND DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Expected versus Observed Genetic Variance
The inbreeding coefficients calculated from published molecular data (Bernardo, 1993; Bernardo et al., 2000) increased from FBSSS = 0 in BSSS, to FB73xB37 = 0.18 in B73 x B37, and to FLo904xLo916 = 0.55 in Lo904 x Lo916. This result indicated that, relative to the original VTA in BSSS (100%), the VTA was expected to be reduced to 82% in B73 x B37, and to 45% in Lo904 x Lo916, assuming no epistasis (Falconer and Mackay, 1996, p. 264).

For grain yield, the pooled estimates of VTA across two selfing generations were all positive (Table 2). The sizes of VTA in BSSS, B73 x B37, and Lo904 x Lo916 showed a similar pattern for the LH185 and Syngenta testcrosses: a decrease from BSSS to B73 x B37, followed by an increase from B73 x B37 to Lo904 x Lo916. The pooled estimates across three testers also followed the same pattern (Table 2, Fig. 1) . For the LH295 testcrosses, the sizes of VTA showed a different pattern: an increase from BSSS to B73 x B37, followed by a decrease from B73 x B37 to Lo904 x Lo916. The differences, however, were not significant among the pooled estimates of VTA across two selfing generations for each tester and across the three testers.


View this table:
[in this window]
[in a new window]
 
Table 2. Pooled estimates of testcross additive variances (VTA) in the BSSS, B73 x B37, and Lo904 x Lo916 maize populations, across selfing generations and testers.

 


View larger version (23K):
[in this window]
[in a new window]
 
Fig. 1. Expected (dashed line) and observed (solid line) testcross additive variance (VTA) for grain yield, grain moisture, plant height, ear height, and stalk lodging in the three maize populations. Estimates were pooled across selfing generations and testers.

 
For grain moisture, the pooled estimates of VTA across two selfing generations were significantly smaller in B73 x B37 than in BSSS. For all three testers, VTA consistently decreased from BSSS to B73 x B37 and then increased in Lo904 x Lo916. This pattern was similar to that of VTA for grain yield (Fig. 1, Table 2). The pooled estimate of VTA across testers increased significantly from B73 x B37 to Lo904 x Lo916 (Table 2).

For stalk lodging, plant height, and ear height, the pooled estimates of VTA across the two selfing generations followed a common pattern: a decrease in B73 x B37 and a further decrease in Lo904 x Lo916 (Table 2, Fig. 1). For stalk lodging, the pooled estimates of VTA across two testers decreased significantly from BSSS to B73 x B37 (Table 2). For plant height and ear height, the pooled estimates of VTA across the two selfing generations decreased significantly from BSSS to B73 x B37 (Table 2).

Across all five traits, the behavior of VTA can be summarized as follows. As expected, VTA decreased from BSSS to B73 x B37 for all traits. But VTA in Lo904 x Lo916 was not always less than VTA in B73 x B37. The VTA in B73 x B37 was lower than expected for all traits (Fig. 1). The VTA in Lo904 x Lo916 was lower than expected for grain moisture, plant height, ear height, and stalk lodging, but not for grain yield.

The classical additive model predicts a decrease in genetic variance after a population bottleneck. Several studies, however, showed that additive genetic variance after a population bottleneck was significantly higher than expected from the additive model (Bryant et al., 1986; Wade et al., 1996; Cheverud et al., 1999). Our results showed that traits react differently to population bottlenecks. At a given bottleneck, VTA may decrease for one trait but increase for another. Meanwhile, for a given trait, VTA may decrease after one bottleneck but increase after the next bottleneck. Reviews of experimental and theoretical studies suggested that bottlenecks affect traits in different ways depending on the relative amount of additive and nonadditive genetic variances (Wade et al., 1996; Meffert, 1999). Fertility in flour beetle [Tribolium castaneum (Herbst)], for example, showed increased genetic variance after the bottleneck, but pupal weight showed decreased genetic variance as predicted by an additive model (Wade et al., 1996). Conversion of additive-by-additive epistatic variance into additive variance was proposed as a mechanism for maintaining additive variance in a bottlenecked population (Cockerham and Tachida, 1988; Goodnight, 1988; Whitlock et al., 1993; Cheverud and Routman, 1996). For example, for two loci with duplicate dominant epistasis, random fixation of either locus would lead to an average increase of 300% in additive variance (Bernardo, 2002, p. 115). The changes in VTA for grain yield and grain moisture in our study deviated prominently from a continuously decreasing trend after two population bottlenecks, which was expected from the testcross additive model. This deviation suggested that testcross nonadditive effects, which can be converted into additive genetic variance during a bottleneck, might have been involved.

Testcross Additive versus Testcross Epistatic Models
For all traits examined, the testcross epistatic model did not fit the data for any tester because it failed to give positive estimates of both VTA and VTAA (data not shown). Theoretical studies suggested that even a small epistatic variance could have substantial effect in maintaining additive variance in a bottlenecked population (Goodnight, 1988; Bernardo, 2002, p. 116). The negative estimates of VTA or VTAA in the testcross epistatic model did not allow us to verify the maintenance of additive variance through the conversion of epistatic variance. Although physiological epistasis undoubtedly exists, the general conclusion drawn from empirical studies is that realistic estimates of epistatic variance are difficult to obtain (Hallauer and Miranda Filho, 1988, p. 140).

General reasons that lead to the difficulties in detecting a significant epistatic variance were summarized by Bernardo (2002)(p. 141). First, a complex mating design is needed to obtain multiple covariances between relatives. Second, epistatic variances are by definition expected to be smaller than additive variance and dominance variance. Third, epistatic variances are difficult to separate from additive variance and dominance variance. The large sampling error variance associated with estimates of epistatic variance is also another reason for nonsignificance of epistatic variance (Mackay, 2001); the nonsignificance of epistatic variance does not suggest an absence of epistasis.

Genetic Variance and Advanced Cycle Breeding
Advanced cycle breeding is generally believed to improve mean performance at the cost of reducing genetic variance. Our results confirmed the improvement in mean performance: higher grain yield and fewer plants lodged in Lo904 x Lo916 than in BSSS, and acceptable grain moisture, plant height, and ear height in Lo904 x Lo916 (Table 3). The testcross means (pooled across selfing generations and testers) indicated a 30% increase in grain yield from BSSS to B73 x B37, and a 9% increase in grain yield from B73 x B37 to Lo904 x Lo916. The results in testcross means for grain yield and grain moisture were in agreement with the improvement documented for BSSS, B73, B37 (Troyer, 1999), Lo904, and Lo916 (Bertolini et al., 1991).


View this table:
[in this window]
[in a new window]
 
Table 3. Pooled estimates of testcross means of the BSSS, B73 x B37, and Lo904 x Lo916 maize populations, across selfing generations and testers.

 
With regard to reduction in genetic variance, our study showed that, except for grain yield, genetic variance was lost during advanced cycle breeding (Table 2) and the amount of loss was greater than expected (Fig. 1). For stalk lodging, the negative estimates of VTA in Lo904 x Lo916 indicated the difficulty for further genetic improvement beyond the current level. But for grain yield, the lack of significant differences among VTA in BSSS, B73 x B37, and Lo904 x Lo916 may explain the increases in yield that have been achieved so far through advanced cycle breeding.

Molecular marker studies of progenitors and elite lines from BSSS (Messmer et al., 1991; Bernardo et al., 2000) confirmed that the number of alleles per locus has been reduced in B73 x B37 and further reduced in Lo904 x Lo916 compared with the original BSSS. In practice, advanced cycle breeding involves crossing many pairs of elite inbreds, for example, B14 and B84, in addition to B73 and B37. This would result in a slower rate of fixation of alleles than in BSSS to B73 x B37, and B73 x B37 to Lo904 x Lo916. But, if the inbreds all belong to the same genetic background, those at advanced cycles would have an increasingly higher proportion of loci with fixed alleles, which could eventually lead to a lower genetic variance even if additive genetic variance is maintained during the first few cycles. For any crop species, it is perhaps premature to predict an exhaustion of genetic variance, but broadening the genetic base beyond the current level would help to sustain long-term genetic gains (Hallauer, 1990; Hoisington et al., 1999). Brown (1975) suggested that, on a global basis, only about 5% of the available maize germplasm is used commercially. The amounts of available but unused germplasm for other major crops are also huge (Hoisington et al., 1999). A strategy of maintaining favorable gene combinations through advanced cycle breeding, while increasing genetic variation through introgressing new germplasm into elite germplasm pools, has led to promising results in wheat (Carver and Rayburn, 1994), tomato (Lycopersicon esculentum Mill.; Bernacchi et al., 1998), and barley (Peel and Rasmusson, 2000). But because of its proven success, breeders, particularly in private companies, will probably continue to practice advanced cycle breeding until it becomes obvious that genetic variance has been exhausted and that new germplasm sources are needed.

Received for publication April 16, 2003.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 




This article has been cited by other articles:


Home page
Crop Sci.Home page
D. A. Tabanao and R. Bernardo
Genetic Variation in Maize Breeding Populations with Different Numbers of Parents
Crop Sci., September 23, 2005; 45(6): 2301 - 2306.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (3)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Yu, J.
Right arrow Articles by Bernardo, R.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Yu, J.
Right arrow Articles by Bernardo, R.
Agricola
Right arrow Articles by Yu, J.
Right arrow Articles by Bernardo, R.
Related Collections
Right arrow Crop Genetics
Right arrow Maize


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
The SCI Journals Agronomy Journal Vadose Zone Journal
Journal of Plant Registrations Soil Science Society of America Journal
Journal of Natural Resources
and Life Sciences Education
Journal of
Environmental Quality