|
|
||||||||
a Inst. of Plant Breeding, Seed Sci., and Population Genetics, Univ. of Hohenheim, 70593 Stuttgart, Germany
b Int. Maize and Wheat Improvement Center (CIMMYT), Apdo. Postal 6-641 06600 Mexico D.F., Mexico
* Corresponding author (melchinger{at}uni-hohenheim.de)
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: ANOVA, analysis of variance AMOVA, analysis of molecular variance CIMMYT, International Maize and Wheat Improvement Center GCA, general combining ability MRD, modified Roger's distance PC, principal coordinate PCoA, principal coordinate analysis PMPH, panmictic midparent heterosis Pop, population SCA, specific combining ability SSR, simple sequence repeat
| INTRODUCTION |
|---|
|
|
|---|
Most evidence in maize suggests that the genetic basis of heterosis is partial to complete dominance (Hallauer et al., 1988; Stuber et al., 1992). Overdominance has long been discussed as the basis of heterosis (East, 1936; Crow, 1948). However, many data supporting overdominance presumably resulted from pseudooverdominance, arising from dominant alleles in repulsion phase linkage (Stuber et al., 1992; Crow, 1999). Epistasis, particularly between linked loci, may also be an explanation for heterosis in maize (Cockerham and Zeng, 1996). No data exclude the possibility of all three mechanisms contributing to heterosis, albeit in different proportions.
Lamkey and Edwards (1999) coined the term panmictic midparent heterosis to describe the deviation in performance between a population cross and the mean of its two parent populations in Hardy-Weinberg equilibrium. Quantitative genetic theory shows that in the absence of epistasis and two alleles per locus, PMPH is a function of the product of the dominance effect and the square of the difference in gene frequencies at the respective locus (Falconer and Mackay, 1996, p. 255), which corresponds to the square of the MRD (Melchinger, 1999). In fact, a linear increase in PMPH with increasing genetic distance (Hypothesis 1) was hypothesized in a diallel of U.S. maize populations (Moll et al., 1962).
In contrast, experimental data reported by Moll et al. (1965) in a study with tropical maize populations of diverse geographic origin suggest that PMPH increases with increasing genetic distance only up to an optimum level but thereafter decreases in extremely wide crosses (Hypothesis 2). The authors explained this by fertility distortion in wide crosses and epistatic interactions of genes. While Moll et al. (1962)(1965) inferred the genetic distance from the geographic origin of the populations, to our knowledge no attempts have been made to verify or falsify the above hypotheses with more reliable data based on molecular markers.
The choice of heterotic groups is fundamental in hybrid breeding of maize (Melchinger and Gumber, 1998). While heterotic patterns in temperate maize have been established more than 50 yr ago, a clearly defined heterotic pattern does not exist in the tropical maize of the CIMMYT germplasm. Therefore, before embarking on a hybrid breeding program, CIMMYT conducted several diallel studies for identifying populations showing not only good per se performance but also high heterosis in their crosses (Beck et al., 1990; Crossa et al., 1990; Vasal et al., 1992a,b,c). Genetic distances based on molecular markers have been suggested as a tool for grouping of similar germplasm as a first step in identifying promising heterotic patterns (Melchinger, 1999).
The major goal of this study was to investigate the relationship between heterosis and genetic distance determined with SSR markers. The objectives of our research were to (i) compare the genetic diversity within and between seven tropical maize populations, (ii) test alternative hypotheses on the relationship between PMPH and genetic distances determined with SSR markers, and (iii) evaluate the use of SSR markers for grouping of germplasm and establishing heterotic patterns for hybrid breeding of tropical maize.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
53 333 plants ha-1. All rows were hand-harvested and grain yield was calculated from dry ear weight at harvest assuming 80% shelling and adjusted to 155 g kg-1 grain moisture.
Simple Sequence Repeat Analyses
From each of the seven populations, 48 randomly chosen individuals were analyzed separately. The seeds used for extracting DNA were from the same selection cycle as the populations tested in the field trials; however, the populations were multiplied repeatedly by CIMMYT's maize genebank since 1985.
DNA was extracted employing the CTAB procedure (Clarke et al., 1989). The 85 SSR markers were chosen from the MaizeDB database (http://nucleus.agron.missouri.edu/cgi-bin/ssr_bin.pl) based on repeat unit and bin location to provide uniform coverage of the entire maize genome. The SSRs were multiplexed for maximum efficiency. Fragments were separated using acrylamide gels run on an ABI 377 automatic DNA sequencer. Fragment sizes were calculated with GeneScan 3.1 (Perkin Elmer/Applied Biosystems) using the Local Southern sizing method; allele identity was assigned using Genotyper 2.1 (Perkin Elmer/Applied Biosystems) and the two inbred lines CML51 and CML292 as control. Data have been stored in the MaizeDB database (http://nucleus.agron.missouri.edu/cgi-bin/ssr_bin.pl).
Statistical Analyses
Analyses of variance (ANOVA) were computed for the three plant traits. A mixed linear model was used with the assumption that effects of entries were fixed and all other effects were considered random. Following Analysis III of Gardner and Eberhart (1966), the sums of squares and degrees of freedom (27 df) for entries were orthogonally partioned into the contrast between parents vs. crosses (1 df), the variation among populations (6 df), and the variation among crosses (20 df) with a further subdivision into general combining ability (GCA) and specific combining ability (SCA) effects. A corresponding subdivision was made on the entry x environment interaction sums of squares. Entry mean squares were tested by F tests for significance by using the corresponding entry x environment mean squares. Entry x environment mean squares were tested for significance by using the pooled error mean square. The PMPH of each cross was calculated as the difference between the F1 mean and the respective midparent mean across all environments.
The gene diversity (D) based on SSR data was calculated for each population according to Weir (1996)(p. 151):
![]() | [1] |
We calculated the MRD between two populations or individuals (Wright, 1978, p. 91; Goodman and Stuber, 1983) as:
![]() | [2] |
Associations among the populations were revealed with principal coordinate analysis (PCoA) (Gower, 1966) based on MRD estimates. Multiple regression analysis was used to study the relationship between PMPH and squared modified Roger's distance (MRD2). The PCoA was performed with the statistical software R (Ihaka and Gentleman, 1996) and multiple regression analysis with the statistical software SAS (SAS Institute, 1988).
| RESULTS |
|---|
|
|
|---|
|
Simple Sequence Marker Data
The 85 SSR primers generated a total of 532 alleles in the 336 genotypes (7 populations x 48 individuals) analyzed. The number of alleles per marker across all seven populations was on average 6.3 and ranged from 2 to 16 (Table 3). Gene diversity D within the seven populations ranged from 0.503 to 0.580 with a mean of 0.539 (Table 3). Values of MRD between pairs of populations averaged 0.258 and ranged from 0.203 (Pool24 x Pop22) to 0.318 (Pop32 x Pop43) with significant (P < 0.01) differences between MRD estimates (Table 4). The AMOVA revealed that 89.8% of the molecular genetic variance was found within populations and 10.2% between populations (Table 5).
|
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
Genetic Diversity among and within the Populations
In this study, we found on average across the seven populations 6.3 alleles per marker. Lu and Bernardo (2001) detected for 40 U.S. inbred lines an average of 4.9 alleles using 83 SSR markers. Senior et al. (1998) reported an average of five alleles in their study with 94 elite maize inbreds, representative of the diversity in the U.S. maize germplasm, and 70 SSR markers. Hence, the total number of alleles per marker was higher in our study than previously reported in the literature. This and the high average number of alleles per population (Table 3) in our study suggests a broad genetic base of the seven populations.
Pop29 had the highest gene diversity D followed by Pool24 and Pop21 (Table 3). This is consistent with pedigree information (Table 1) because the populations have been established using a wide range of germplasm. The lowest D value observed for Pop43 is also in accordance with its pedigree, because it was generated from 16 S1 lines including only Tuxpeño germplasm. Ranking of the populations based on D was almost identical with their ranking based on the average number of alleles per marker (rank correlation rs = 0.93; P < 0.01). Altogether, the high percentage (89.8%) of the molecular variance revealed by the AMOVA (Table 5) within populations is in harmony with the broad genetic base of the materials used for their synthesis (Table 1). Since related germplasm such as various sources from Tuxpeño or ETO entered different populations, it was also not surprising to find only a minor variance between populations (Table 5). A more detailed analysis of the population subdivision with test statistics of the AMOVA was not possible, because this would require knowledge of the gametic phase for linked loci (Michalakis and Excoffier, 1996), which cannot be determined from SSR analyses of heterozygous individuals.
Correlation between MRD2 and Panmictic Midparent Heterosis
We investigated the correlation between PMPH and MRD2 because quantitative genetic theory suggests a linear relationship between both measures under certain assumptions (Falconer and Mackay, 1996, p. 255). This is in harmony with related studies on midparent heterosis in crosses of inbred lines (see Melchinger et al., 1991; Boppenmaier et al., 1993), where the commonly employed Roger's distance (1972) is equal to MRD2 (Melchinger, 1993). A low correlation between PMPH and MRD2 can be attributable to several causes: (i) a poor association between heterozygosity estimated from marker data and heterozygosity at quantitative trait loci affecting the trait examined, (ii) a poor association between heterozygosity and heterosis at quantitative trait loci in the crosses examined (Charcosset et al., 1991), (iii) existence of multiple alleles (Cress, 1966), and (iv) epistasis (Moll et al., 1965).
The low correlations between MRD2 and PMPH for plant height and days to silking were mostly due to small PMPH estimates for these traits (Table 2). By comparison, the corresponding correlation for grain yield was surprisingly high (r = 0.63; P < 0.01). This is consistent with the relative large contribution of SCA effects to the total sums of squares, which accounted for 33% of the genetic variation among crosses for this trait (Vasal et al., 1992a). In accordance with quantitative genetic theory (Melchinger, 1999) the correlation of MRD2 was lower with hybrid performance (r = 0.41; P < 0.05) than with PMPH for grain yield (r = 0.63; P < 0.01). On the basis of a literature survey with single crosses produced from inbreds, Melchinger (1999) pointed out that only intragroup crosses show a correlation between parental genetic distance and midparent heterosis, but not intergroup crosses. However, a closer examination of the graph between MRD2 and PMPH (Fig. 3) did not provide any clue in this direction.
While Hypothesis 1 postulates a linear relation between MRD2 and PMPH, under Hypothesis 2 a quadratic or cubic regression is expected to fit the data better than linear regression. However, in our study neither a quadratic nor a cubic regression model gave a significantly better fit to the data than linear regression. This is in accordance with the graphs shown in Fig. 3. Consequently, our results confirm Hypothesis 1 for the tropical maize germplasm investigated here.
A decrease in PMPH of genetically very distant populations is generally attributed to the lack of coadaption between both allelic and nonallelic combinations of genes from the two parental haploid genomes, resulting in reduced or negative dominance and negative epistatic effects, respectively (Falconer and Mackay, 1996, p. 255). A major reason for the absence of an optimum in the relationship between genetic distance and PMPH in our study could be that all populations (Table 1) were more or less well adapted to the test environments. In addition, we did not include extremely wide crosses, as was the case in the experiment of Moll et al. (1965).
For hybrid breeding, Melchinger and Gumber (1998) recommended the following criteria for the choice of heterotic patterns: (i) high mean performance and large genetic variance in the hybrid population; (ii) high per se performance and good adaption of the parent populations to the target region(s); (iii) low inbreeding depression, if hybrids are produced from inbreds. Under Hypothesis 1 (PMPH increases with increasing genetic distance), genetic distance could be used as a further criterion for the identification of heterotic patterns. Considering all four criteria, the following promising heterotic patterns can be suggested: (i) Heterotic Group A with Heterotic Group B; (ii) Pop43 with Heterotic Group A or B; (iii) Pop29 with Heterotic Group B or Pop43.
Grouping of Germplasm
We chose the MRD as genetic distance measure because of its mathematical and genetic properties. In particular, it is an Euclidean distance, which is an often-overlooked prerequisite for most multivariate analysis methods (Jacquard, 1974, p. 465). Furthermore, in the absence of epistasis and two alleles per locus, PMPH is a linear function of the product of the dominance effect and the square of the MRD (Melchinger, 1999).
Principal coordinate analysis based on MRD revealed very clearly a major split between the populations from Heterotic Group A and Pop32 (Fig. 1). Pop25 is separated from the other populations by PC3 and had an average MRD at the population level to Heterotic Group A of 0.24 and to Pop32 of 0.26. The assignment of Pop25 to Heterotic Group B together with Pop32 originally based on testcross data was not supported by our molecular data. This could be interpreted as an indicator that Pop25 should have been established as a separate Heterotic Group C. The values of PMPH (Table 2) support this hypothesis in that Pop25 had a low average PMPH with Heterotic Group A. In addition, PCoA accurately portrayed the relationship of Pop43 to Heterotic Group A and B. It is closer to Heterotic Group A (
= 0.26) than to Heterotic Group B (
= 0.29), but the distance from Pop43 to Heterotic Group A was higher than the average distance between Heterotic Groups A and B. This together with the diallel analysis suggests classification of Pop43 as a separate Heterotic Group D. According to the PCoAs (Fig. 1, 2), Pop29 could be assigned to Heterotic Group A, because it had a smaller average MRD to Heterotic Group A (0.22) than to B (0.26). The diallel analysis supports this suggestion.
In conclusion, classification of the seven populations based on SSR data mostly confirmed the results from the diallel data except the assignment of Pop25 to Heterotic Group B. Furthermore, it was possible to assign Pop29 to the established Heterotic Groups A and to propose new heterotic groups (Pop25, Pop43). When a large number of germplasm exists but no established heterotic groups are available, genetically similar germplasm can be identified with molecular markers. On basis of this information, field trials can be planned more efficiently. Thus, by using molecular data to focus the search for heterotic groups on a smaller number of promising heterotic patterns and evaluating these intensively, breeders should arrive at a more economic and solid approach for making this important decision at the beginning of a hybrid breeding program.
| CONCLUSIONS |
|---|
|
|
|---|
Simple sequence repeat markers provide a valuable tool for grouping of germplasm and are a good complementation to field trials for identifying groups of genetically similar germplasm. Consequently, field trials for identification of promising heterotic patterns can be planned more efficiently based on prior information obtained from SSR analyses.
| ACKNOWLEDGMENTS |
|---|
Received for publication June 7, 2002.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
E. A. Lee, M. J. Ash, and B. Good Re-examining the Relationship between Degree of Relatedness, Genetic Effects, and Heterosis in Maize Crop Sci., March 1, 2007; 47(2): 629 - 635. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. C. Xia, J. C. Reif, A. E. Melchinger, M. Frisch, D. A. Hoisington, D. Beck, K. Pixley, and M. L. Warburton Genetic Diversity among CIMMYT Maize Inbred Lines Investigated with SSR Markers: II. Subtropical, Tropical Midaltitude, and Highland Maize Inbred Lines and their Relationships with Elite U.S. and European Maize Crop Sci., October 27, 2005; 45(6): 2573 - 2582. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Chen and R. L. Nelson Relationship between Origin and Genetic Diversity in Chinese Soybean Germplasm Crop Sci., June 24, 2005; 45(4): 1645 - 1652. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. C. Reif, A. E. Melchinger, and M. Frisch Genetical and Mathematical Properties of Similarity and Dissimilarity Coefficients Applied in Plant Breeding and Seed Bank Management Crop Sci., January 1, 2005; 45(1): 1 - 7. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. C. Xia, J. C. Reif, D. A. Hoisington, A. E. Melchinger, M. Frisch, and M. L. Warburton Genetic Diversity among CIMMYT Maize Inbred Lines Investigated with SSR Markers: I. Lowland Tropical Maize Crop Sci., November 1, 2004; 44(6): 2230 - 2237. [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 | |||