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 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 Similar articles in ISI Web of Science
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 (6)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Dudley, J. W.
Right arrow Articles by Rocheford, T. R.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Dudley, J. W.
Right arrow Articles by Rocheford, T. R.
Agricola
Right arrow Articles by Dudley, J. W.
Right arrow Articles by Rocheford, T. R.
Related Collections
Right arrow Crop Genetics
Right arrow Maize
Right arrow Biometrics
Published in Crop Sci. 44:1419-1428 (2004).
© 2004 Crop Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA

GENOMICS, MOLECULAR GENETICS & BIOTECHNOLOGY

Effects of Random Mating on Marker–QTL Associations in the Cross of the Illinois High Protein x Illinois Low Protein Maize Strains

J. W. Dudleya,*, A. Dijkhuizenb, C. Paula, S. T. Coatesc and T. R. Rocheforda

a Dep. of Crop Sci., Univ. of Illinois, 1102 S. Goodwin Ave., Urbana, IL 61801
b Monsanto SAS, Route D'Epincy, Louville, La Chenard 28150, France
c 311 N. Osborn Lane, Aberdeen, MD 21001

* Corresponding author (jdudley{at}uiuc.edu).


    ABSTRACT
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
The effects of random mating on population performance and identification of marker–quantitative trait loci (QTL) associations were studied in the cross of the Illinois High Protein (IHP) and Illinois Low Protein (ILP) maize (Zea mays L.) strains. The F1 was randomly mated to produce the Syn0 and Syn4 generations. Two hundred Syn0 S1 and 200 Syn4 S1 lines were crossed to two inbred testers. The S1 lines per se and the testcrosses (TCs) were evaluated for grain protein, starch, and oil concentrations in four environments. Genetic variance for protein and starch was reduced from the Syn0 to the Syn4 and genetic variance for the S1 lines per se was greater than for the TCs. Heritability was similar in the Syn0 and Syn4. With single factor (SF) analysis, the number of significant marker–QTL associations for protein and starch in the Syn4 was drastically reduced from the number found in the Syn0. Multiple regression (MR) analysis supported this result. Very few significant marker–QTL associations for oil were found in either the Syn0 or Syn4. Agreement between the testers and the S1 lines per se for chromosomal regions identified as associated with QTL for starch and protein was high. Genotypic correlations of starch with protein ranged from –0.96 to –0.98. However, three markers had the same sign and significant additive effects for both starch and protein in the S1 progeny and both TCs. Such markers should be useful for simultaneous selection for high starch and high protein.

Abbreviations: BLUP, best linear unbiased predictor • IHO, Illinois High Oil • IHP, Illinois High Protein • ILO, Illinois Low Oil • ILP, Illinois Low Protein • MR, multiple regression • QTL, quantitative trait loci • RFLP, restriction fragment length polymorphism • SF, single factor • SSR, simple sequence repeat • TC, testcross


    INTRODUCTION
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
THE PRESENCE of linkage disequilibrium between molecular markers and QTL is necessary for identification of marker–QTL associations (Liu, 1998). Maximum linkage disequilibrium is obtained in the F2 or backcross one generation. However, use of such early generations makes discrimination among closely linked QTL difficult. Random mating within a mapping population breaks up linkage blocks and reduces the average size of an unbroken linkage segment (Hanson, 1959). Linkage disequilibrium decays, leaving only close linkages detectable after a few generations of random mating (Wright, 1969). Thus, markers significantly associated with QTL in a randomly mating population should be those most closely linked to the QTL. Random mating also should allow breakup of linked QTL. In the case of repulsion phase linkages, this could allow identification of QTL, which might be missed in an early generation because of canceling effects caused by linkage of QTL of opposite sign. In the case of coupling phase linkages, random mating should allow for more precise estimation of QTL effects because the presence of two QTL linked in coupling would lead to overestimation of the effect of a single QTL and a misleading estimate of the location of the QTL.

Dudley (1993) suggested the use of randomly mated progeny to increase recombination between markers, between QTL, and between markers and QTL to increase resolution. Liu et al. (1996) showed empirically that much higher resolution could be obtained among tightly linked markers after four generations of randomly mating an F2. Darvasi and Soller (1995), using simulation, demonstrated a reduction in the confidence interval surrounding a QTL by successive generations of random mating. Lee et al. (2002), with 190 RFLP loci, showed a nearly four-fold increase in genetic map distance following five generations of intermating the cross of maize inbreds B73 and Mo17. Recombinant inbreds were derived from this population after four generations of randomly mating the F2 population. Lee et al. (2002) suggested this resource should connect research on dense genetic maps, physical mapping, gene isolation, comparative genomics, and analysis of QTL. Winkler et al. (2003) showed the relationship between crossover rates in an F2 and expected recombination in a population of recombinant inbreds developed from a randomly mated population tracing to the F2. Lu et al. (2003) used random mating to break up repulsion linkages in an attempt to separate pseudooverdominance from true overdominance for grain yield in maize.

Crosses between very divergent parents are desirable to determine the genetic architecture of a trait and identify many of the QTL affecting the trait. The IHP and ILP strains of maize are ideally suited for this purpose. These strains were divergently selected for high and low protein concentration in the grain beginning in 1896 (Dudley and Lambert, 2004). At the time of initiation of the experiment reported in this paper, ILP had reached a lower limit of protein concentration, although progress was still being made in IHP. In a Design III study with the cross of generations 70 of IHP and ILP, Dudley (1994) demonstrated a significant reduction in additive genetic variance for protein concentration following four generations of random mating. This result is what is expected with coupling phase linkages (Comstock and Robinson, 1948).

Increasing protein concentration in maize grain usually results in reduced starch concentration (Dudley and Lambert, 2004). Thus, in generations 70, the generations from which the parents of this study came, IHP (440 g kg–1 starch, 261 g kg–1 protein) and ILP (745 g kg–1 starch, 58 g kg–1 protein) (Dudley and Lambert, 1992) were low and high starch strains as well as being high and low protein strains. In a previous study, Goldman et al. (1993) found 19 significant marker–QTL associations for starch, of which 16 were also significant for protein. Whether these shared associations were the result of linkage among trait-specific QTL or single QTL with pleiotropic effects could not be determined. Information on the genetic control of starch and protein is desirable because of the importance of increased starch concentration, with accompanying reduced protein concentration, for increasing efficiency of ethanol production, and because of the importance of increased protein concentration in grain for cattle feeding.

The objectives of this study were (i) to determine the effects of random mating on the ability to identify marker–QTL associations in the cross of generations 70 of IHP and ILP; (ii) to compare S1 per se performance with performance of TCs to two different testers for ability to identify marker–QTL associations; and (iii) to use this information to identify chromosomal regions containing QTL controlling starch, protein, and oil concentration in maize grain.


    MATERIALS AND METHODS
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Genetic Materials
The parental materials for this study were generations 70 of IHP and ILP. Illinois High Protein and ILP resulted from 70 generations of selection for high protein concentration (IHP) and 70 generations of selection for low protein concentration in the grain (ILP), respectively. Selection was initiated in the open-pollinated maize cultivar ‘Burr’s white' in 1896. Selection procedures have been described (Dudley, 1974; Dudley and Lambert, 1992, 2004). Production of the F1 and randomly mated generations was described by Dudley (1994). Briefly, five to seven plants from IHP were crossed to five to seven plants from ILP to produce F1 seed. The F1 plants were randomly mated to produce the F2, which was designated the Syn0 generation. The Syn0 generation was randomly mated four times to produce the Syn4. Approximately 200 plants were used for random mating each generation. Progenies for the current study were produced by selfing 200 plants in the Syn0 and 200 plants in the Syn4 to produce Syn0 S1 and Syn4 S1 lines. The sample of plants selfed in the Syn4 was independent of the sample selfed in the Syn0 because the plants selfed in the Syn0 were not those used to randomly mate the Syn0. Testcross progenies were produced by crossing each of the S1 lines to two commercial inbreds, FR1064 (a Stiff Stalk Synthetic derivative) and FR616 (a non-Stiff Stalk line). Seeds of FR1064 and FR616 were provided by Illinois Foundation Seeds, Inc., Champaign, IL.

Evaluation
In 1994 and 1995, the 200 Syn0 S1 and 200 Syn4 S1 lines were evaluated as lines per se in adjacent experiments at the Crop Sciences Research and Education Center near Urbana, IL. Each year, each set of lines was grown at two different planting dates approximately 3 wk apart to provide two environments in each year. To control environmental variation within a year–environment combination, families were arranged in a generalized lattice [{alpha}(0,1)] design with two replications of 20 incomplete blocks of 10 entries each (Carmer, 1985). Plots of the Syn0 experiment were grown adjacent to those of the Syn4 experiment. Plots consisted of a single 4.5-m row with 0.76 m between rows. Plots were overplanted and thinned to 15 plants per row. In each row, six plant-to-plant crosses were made to produce seed for chemical analysis. Ears from the hand-pollinated plants were harvested by hand.

In 1996 and 1997, TCs to FR1064 and FR616 of the 200 Syn0 and Syn4 lines were grown. Testcrosses were grown at two locations near Urbana, IL, each year and the Syn0 and Syn4 TCs to a particular tester were grown in adjacent experiments. Each year–location combination was arranged in a generalized lattice [{alpha}(0,1)] design with two replications of 20 incomplete blocks of 10 entries each. Plots consisted of 2 rows 5.3 m long with 0.76 m between rows. Plots were overplanted and thinned to 23 plants per row. Plots were machine-harvested and a sample of grain collected from the combine for chemical analysis. Grain was weighed and moisture data collected on the combine.

Protein, starch, and oil concentrations in the grain were measured with a DICKEY-john near-infrared analyzer (Hymowitz et al., 1974), as described by Dudley and Lambert (1992). Thus, the phenotypic data set consisted of 14400 data points.

Molecular Marker Data
Molecular marker genotypes for the 200 S1 lines from the Syn0 and the 200 lines from the Syn4 were obtained for 44 RFLP and 20 SSR markers. DNA was isolated from bulks of tissue from 25 random plants from each S1 line according to established procedures (Mikkilineni, 1997). The RFLP procedures (Berke and Rocheford, 1995; Dijkhuizen et al., 1998) and SSR procedures (Wong et al., 2003) used have been described previously and are based on protocols reported by Hoisington et al. (1989) and Senior et al. (1996). The RFLP and SSR probes are part of larger sets listed in the Maize Genomics and Genetics Database (http://www.maizegdb.org/; verified 4 Mar. 2004). In addition, the cDNA clone of shrunken2 (sh2) (Bhave et al., 1990) was used as a probe.

Markers were chosen on the basis of polymorphism between IHP and ILP (based on differences in allelic frequencies between IHP and ILP as measured on a sample of plants from generations 70), previous association with QTL for kernel composition in the cross of generations 76 of IHP and ILP (Goldman et al., 1993), and coverage of all chromosome arms in the maize genome. Because the parents were not homozygous and the allelic frequencies in the parents were not known, it was not possible to assign alleles on the basis of their origin as being from IHP or ILP. Consequently, the fragments were simply scored by relative molecular weight and given a letter designation of A, B, C from highest molecular weight to lowest. For some markers, different single alleles were fixed in IHP and ILP and designations were simply A and B, scored in the same manner.

Statistical Analysis
Analyses of variance were done and relevant variance components estimated for each generation (Syn0 and Syn4) of each progeny type (S1 progeny, TC to FR1064, TC to FR616). PROC MIXED in SAS (SAS Institute, 1999) was used for statistical analysis assuming years and environments (planting dates for S1's; locations for TCs) (E) were fixed effects. Blocks (B), reps (R), and lines (L) were considered random. The interactions of lines with years (LY), lines with environment (LE), and lines x years x environments (LYE) were considered random. Variance components and their 95% confidence intervals were estimated from PROC MIXED by the {alpha} = 0.05 option (SAS Institute, 1999).

Heritabilities (H2), for each generation of each progeny type, were calculated from the equation H2 = {sigma}2g/ (Bernardo, 2002). Confidence intervals (90%) for H2 were calculated as described by Knapp et al. (1984). Best linear unbiased predictor (BLUP) estimates for protein, starch, and oil concentration were obtained for each line for each progeny type by PROC MIXED.

A procedure described by Zehr et al. (1992) was used to test for significance of marker–QTL associations because many of the markers were segregating for multiple alleles and allelic frequencies were not 0.5 in the F1 generation. Briefly, each allele for each marker was analyzed separately. Homozygotes for allele A were designated AA. All families heterozygous for allele A were designated AX regardless of the designation of the alternative allele. Genotypes not AA or AX were designated XX. The BLUP estimates were then used in a SF ANOVA to identify significant associations of markers with QTL for protein, oil, or starch. This procedure tests for a QTL linked to allele A. A similar procedure was used for alleles B, C, etc. for each marker. When one or more alleles showed a significant effect at the 0.05 significance level for a marker, the allele with the lowest probability for significance was selected to represent that marker with the restriction that there had to be at least 10 individuals in each marker genotype class. For markers with only two alleles, the probability for each allele was the same and the A allele was selected in every case. Because the parents were not homozygous and the allelic frequencies in the parents were not known, it was not possible to assign alleles on the basis of their origin as being from IHP or ILP.

Multiple regression models were developed only for starch and protein because of the low genetic variance and small number of significant effects from the SF analysis for oil. All markers with at least one allele significant at the 0.05 level for starch in at least one of the progeny types were selected as the set to be used in MR analysis. The same set of markers was used for analysis of starch and protein because the markers selected based on starch included all the markers with significant effects for protein. This procedure reduced the set of markers included in the regression analysis to 42.

The stepwise option in PROC REG in SAS (SAS Institute, 1999), with a probability of entering or dropping a marker as 0.15, was used in MR analysis. Because missing marker data reduced the number of families available for MR analysis, a sequential procedure was used in that analysis. First, an analysis was done for each chromosome separately. The markers selected from each chromosome in this step were then entered into a similar analysis across all chromosomes to produce the final model. This procedure reduced the number of markers entered into any one analysis, thus maximizing the number of families available for QTL identification.

Effects of random mating were measured by comparing means, genetic variances, heritabilities, and proportion of genetic variance accounted for by the best multiple regression model (p) measured in the Syn0 and Syn4. Estimates of p were obtained by dividing R2 (proportion of phenotypic variability accounted for by the multiple regression model) by H2. Effects of random mating on marker–QTL associations were measured by comparing the number of marker–QTL associations significant at the 0.01 probability level in the Syn0 with the number significant in the Syn4 for each trait and progeny type. Significance of changes in number of significant marker–QTL associations between the Syn0 and Syn4 was measured by a contingency table {chi}2 test (P < 0.05) adjusted for continuity (Snedecor and Cochran, 1989). Marker allelic frequencies in the Syn0 and Syn4 were compared with standard errors of allelic frequencies calculated as described by Weir (1996). Chi-square analysis (Weir, 1996) was used to test for agreement of genotypic frequencies with random-mating expectations in the Syn0 and Syn4. For tests of changes in allelic frequency between generations, the marker loci and alleles used in the MR analysis were used and the genotypes associated with those alleles were used in the {chi}2 analysis. Linkage disequilibrium among markers in each generation was measured by the correlations among markers within each generation. This procedure is equivalent to the method described by Weir (1996)(p. 137). Correlations were calculated for the set of 42 markers used in MR analysis.

Phenotypic correlations among traits were calculated for all possible combinations of traits for each progeny type and generation by PROC CORR in SAS. BLUP estimates of entry values were used in the calculations. Genotypic correlations were calculated with a SAS program developed by Holland (personal communication, 2003, available at http://www4.ncsu.edu/%7Ejholland/correlation/correlation.html). This procedure uses PROC MIXED with the ASYCOV option to calculate genetic variances and covariances between traits. The output is then used in an IML subroutine to calculate genetic correlations and their standard errors by the procedure of Mode and Robinson (1959). Markers showing significant associations for both starch and protein in either the SF or MR analysis were identified to determine the extent of common genetic control of these traits.


    RESULTS AND DISCUSSION
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Means, Correlations, Genetic Variances, and Heritabilities
Although least squares means for the Syn0 and Syn4 generations were significantly (P < 0.05) different, as measured by a t test, for many of the trait progeny combinations (Table 1), the maximum difference between Syn0 and Syn4 means was only 13 g kg –1 for any trait. Starch and oil concentrations were significantly (P < 0.05) higher and protein concentrations lower for the TCs than for the S1's. Protein was significantly (P < 0.05) lower and starch significantly higher for FR1064 TCs than for FR616 TCs. Entry BLUP estimates for starch and protein were highly negatively correlated with phenotypic correlations ranging from –0.78 to –0.95 (Table 2) and genotypic correlations ranging from –0.96 to –0.98. From these results, marker-associated effects for starch and protein are expected to have opposite signs; that is, a QTL allele for high starch is expected to also be an allele for low protein. In contrast to the correlations between starch and protein, both the phenotypic and genotypic correlations of starch with oil were low, even though the phenotypic correlations were significant (P < 0.05). For FR616 TCs, (in both generations) and for the S1 correlations in the Syn4, correlations of oil with protein were not significant. No significant changes in genetic correlations between the Syn0 and Syn4 were found.


View this table:
[in this window]
[in a new window]
 
Table 1. Means across years and planting dates for starch, protein, and oil percentages by generation and tester.

 

View this table:
[in this window]
[in a new window]
 
Table 2. Phenotypic and genotypic (±SE) correlations among starch, protein, and oil by generations and progeny types.

 
The correlations among the S1, FR1064, and FR616 progeny types were all significant (Table 3). They ranged from 0.69 to 0.81 for protein and starch in the Syn0. The correlations in the Syn4 were generally lower ranging from 0.49 to 0.66. For oil, the correlations were also low, ranging from 0.33 to 0.61. The lower correlations involving oil are likely due to the low genetic variance for oil.


View this table:
[in this window]
[in a new window]
 
Table 3. Correlations among progeny types for starch, protein, and oil by generation.

 
Genetic variances for the S1 progenies per se were approximately four times those for either tester for protein and six times those for either tester for starch (Table 4). Confidence intervals for the tester genetic variance estimates did not overlap those for the S1's. Starch and protein genetic variance estimates for FR1064 and FR616 were very similar with overlapping confidence intervals. The results for protein are those expected by theory for additive genetic variance; that is, genetic variance among S1 family means is expected to be four times the genetic variance among TC means. The results observed for starch suggest dominant alleles for some loci involved in the inheritance of starch, because dominant alleles in the tester would reduce TC variance more than expected for an additive model. In contrast to the starch and protein results, genetic variance for oil was reduced from Syn0 to Syn4 only for the S1 progenies (Table 4). In addition, genetic variance for FR616 was significantly less than for FR1064 in both the Syn0 and Syn4. These results suggest FR616 may contain dominant alleles for oil.


View this table:
[in this window]
[in a new window]
 
Table 4. Genetic variance and heritability (H2) estimates for protein, starch, and oil.{dagger}

 
Theory (Comstock and Robinson, 1948; Moreno-Gonzalez et al., 1975; Dudley, 1994) suggests divergent selection should build up coupling linkages for protein in IHP and ILP. With coupling phase linkages in the F1, random mating is expected to reduce the additive genetic variance (Comstock and Robinson, 1948). Although the genetic variance among S1 lines is expected to include some dominance variance as well as additive variance, Dudley (1994) found no significant dominance variance for protein in a Design III study of the cross used in this experiment. In addition, variance among TC means does not include dominance variance. Thus, estimates of genetic variance in the Syn4 are expected to be lower than those from the Syn0 for all progeny types. For protein and starch, in agreement with theory and with the results of Dudley (1994), none of the Syn4 genetic variance estimates based on per se and TC evaluations had 95% confidence intervals which overlapped the corresponding Syn0 estimate. Moreno-Gonzalez et al. (1975) found similar reductions in genetic variance for oil from a cross of Illinois High Oil (IHO) x Illinois Low Oil (ILO). In our study, the difference between Syn0 and Syn4 genetic variance estimates for oil, based on 95% confidence intervals, was only significant for the S1 lines per se. However, indirect selection had likely not built up coupling phase linkages for oil given the low genetic correlations between oil and protein and the relatively small difference in oil between generations 70 of IHP (54 g kg –1 oil) and ILP (42 g kg –1 oil). This difference is much less than that between the generations 70 IHO (160 g kg–1 oil) and ILO (<10 g kg–1 oil) strains used by Moreno-Gonzalez et al. (1975).

Heritabilities for protein and starch ranged from 0.81 to 0.93 with no significant difference between the Syn0 and Syn4 and no difference among progeny types (Table 4) based on overlapping 90% confidence intervals. For oil, Syn0 and Syn4 heritabilities did not differ significantly. Thus, even though the genetic variance was reduced from the Syn0 to the Syn4 for starch and protein in the S1 and both TCs and for oil in the S1 progenies, heritabilities did not change because the interaction components in the denominator of the heritability equation were also reduced (data not shown).

Allele and Genotypic Frequencies, Linkage Disequilibrium
Significant (P < 0.05) differences in allele frequency between the Syn0 and Syn4 were found for 18 of the 42 markers (Table 5) used in MR analysis. Significant (P < 0.05) deviations from random mating as measured by {chi}2 were found for 14 markers in the Syn0 and five markers in the Syn4. For 11 markers with significant {chi}2 values in the Syn0, the deviation from random mating was not significant in the Syn4. This result might be expected with random mating because markers not in equilibrium in the Syn0 would be expected to go to equilibrium with random mating. On the other hand, there was only one marker with a significant deviation from random-mating equilibrium in the Syn4 which had a nonsignificant deviation in the Syn0.


View this table:
[in this window]
[in a new window]
 
Table 5. Tests for differences in allele frequency and for deviations from random-mating equilibrium for markers used in multiple regression analysis.{dagger}

 
Linkage disequilibrium among markers is expected to be reduced with random mating. Of 88 significant (P < 0.05) correlations among markers in the Syn0, only 28 (32%) were significant in the Syn4, indicating breakup of linkage disequilibrium. However, on chromosome 3 where marker density was increased because significant associations had been found in previous studies, 18 of 24 correlations significant in the Syn0 were also significant in the Syn4, indicating continued disequilibrium among markers.

QTL Analysis Results
Syn0 vs. Syn4—Number of Significant Marker–QTL Associations
A major objective of this study was to compare results of QTL analysis from the Syn0 and Syn4 generations. The Syn4 was expected to have fewer significant marker trait associations than the Syn0. The reasons for this can be illustrated by three examples. In the first example, assume marker density is low. Where marker density is low, a marker may be loosely linked with a QTL having a large effect. In the Syn0, linkage disequilibrium is high and this marker will show a significant effect. Because random mating reduces linkage disequilibrium and the size of unbroken linkage blocks (Hanson, 1959), disequilibrium between a loosely linked marker and a QTL may be reduced to the point that the marker QTL association is no longer detectable. In the second example, assume a region with 10 markers in, for example, a 20-cM distance. Within that region, assume there are two QTL linked in coupling, each with a relatively small effect. In the Syn0 all 10 markers may show significant marker QTL effects. After random mating, perhaps only four of the markers (two flanking each QTL) may show significant effects. Thus, the number of significant marker effects detected is reduced. As a third example, QTL effects measured in the Syn0 could be the result of blocks of genes linked in coupling with a marker located near this block. The effect of each individual gene may be too small to be detectable. When the block of genes is broken up by random mating, the marker QTL association is no longer detectable. The second and third examples are based on coupling phase linkages as is expected in the IHP x ILP cross. If two QTL are linked in repulsion and their effects cancel, random mating may break up the linkage and allow identification of two QTL where none were identified before random mating.

For comparisons of the Syn0 with the Syn4, the number of marker trait associations significant at P < 0.01 was used. This P value was chosen arbitrarily as a compromise between a very stringent P value and a more relaxed P value. The expectation of fewer significant associations in the Syn4 was realized because for protein there were 38, 32, and 32 markers significantly (0.01 probability level) associated with QTL in the Syn0 for the S1, FR1064, and FR616 progenies, respectively, but only 13, 7, and 3 in the Syn4 (Table 6). Similarly, for starch there were 35, 30, and 29 significant associations in the Syn0, and 9, 3, and 6 in the Syn4 for the S1, FR1064, and FR616 progenies, respectively. For each progeny type for starch and protein, the contingency table {chi}2 value (Table 6) comparing the Syn0 and Syn4 was highly significant (P < 0.01), indicating a significant reduction in the number of significant marker QTL associations. For oil, only 6, 3, and 8 markers were significant in the S1, FR1064, and FR616 progenies, respectively, in the Syn0, and only 3, 2, and 0 for the Syn4. Chi-square values for oil were not significant for the S1 and FR1064 progenies. In evaluating these results, it is important to remember that the 200 S1 lines used in the Syn4 are independent of those used in the Syn0 because the sample of parents used to continue random mating from the Syn0 was not the set of lines used to measure effects in the Syn0. In addition, the fact that the standard errors of the means of the Syn0 and Syn4 were very similar (Table 1) and the number of progenies used was the same in both generations suggests that the smaller number of QTL detected in the Syn4 are not the result of a difference in precision between the Syn0 and Syn4.


View this table:
[in this window]
[in a new window]
 
Table 6. Number of marker QTL associations significant (P < 0.01) and nonsignificant for each trait, progeny type, and generation and heterogeneity {chi}2 values comparing the Syn0 and Syn4.

 
The average marker-associated effect is expected to be lower in the Syn4 than the Syn0. To measure this, ratios of Syn4 marker-associated effects to Syn0 marker-associated effects for markers showing significant associations at the 0.01 P level in the Syn0 and at the 0.05 level in the Syn4 were calculated for protein and starch. The 0.05 level was chosen for the Syn4 because there were too few associations significant at the 0.01 level to give meaningful results. For S1 progenies, the average ratios were 0.47 and 0.50 for protein and starch, respectively. For FR616 TC progenies, they were 0.44 and 0.59. For FR1064, the ratio for protein was 0.44 but there were too few markers meeting the criteria for starch to obtain a meaningful value. Thus, where there were significant associations in the Syn0 and Syn4, the marker-associated effects in the Syn4 were lower than in the Syn0 in agreement with theoretical expectations.

Syn0 vs. Syn4—Variability Accounted for in Multiple Regression Analysis
Results from the MR analysis agreed with the SF results. For protein, the final genetic models in the Syn0 accounted for 68 to 73% of the genetic variance. For starch, 52 to 75% of the genotypic variance in Syn0 was accounted for (Table 7). However, for the Syn4 the final regression models accounted for only 12 to 18% of the genetic variance for protein and 3 to 20% of the genetic variability for starch. Thus, with the limited number of genetic markers available in this study, no marker model accounted for a high proportion of the genetic variance in the Syn4. In addition, the numbers of markers included in the final models for starch and protein in the Syn4 were less than half the numbers included in the models for the Syn0 (Tables 8 and 9).


View this table:
[in this window]
[in a new window]
 
Table 7. Proportion of phenotypic (R2) and genotypic variability (p) accounted for by multiple regression models for starch and protein.

 

View this table:
[in this window]
[in a new window]
 
Table 8. Markers significant at the 0.01 level for at least one progeny type in one generation from the single factor (SF) ANOVA or included in one of the final multiple regression (MR) models for protein.

 

View this table:
[in this window]
[in a new window]
 
Table 9. Markers significant at the 0.01 level for at least one progeny type in one generation from the single factor (SF) ANOVA or included in one of the final multiple regression (MR) models for starch.

 
Agreement of Progeny Types
Marker-associated effects in the TCs are expected to be only half the magnitude of those in the S1 progenies assuming an additive model. Thus, a reduction in the number of significant effects between the S1 and TC progenies might be expected. However, there were no significant differences in number of significant marker trait associations between the S1 and either the FR1064 or FR616 progenies (Table 6) in either generation based on {chi}2 tests or between the FR1064 and FR616 progenies (results not shown). The ability of the TC progenies to detect a similar number of significant effects as the S1 progenies, when the magnitude of the effects is expected to be half that of the S1's, resulted from greater precision in the TC data as evidenced by the standard error of means of the TCs being approximately half the standard errors of S1 means (Table 1).

In the Syn0, agreement among progenies as to the markers found significantly associated with traits was high. For both protein and starch, 30 of the 38 markers significant (P < 0.01) in the Syn0 for S1 progenies were also significant in the FR1064 and FR616 progenies (Tables 8 and 9). Only two markers significant in either the FR1064 or FR616 progenies were not significant in the S1 progenies. Twenty-four markers were significant for all three sets of progenies in the Syn0. Agreement was not as good for the Syn4 progenies likely because of the small number of significant markers in the Syn4 for any progeny type.

Identification of QTL
Criteria for declaring a marker associated with a QTL included significance (P < 0.01) of marker trait associations in the SF analysis of the Syn0 for all three progeny types, retention in the Syn0 multiple regression model, significance in the SF analysis in the Syn4 of markers found significant in the Syn0, or inclusion in the Syn4 multiple regression model. The number of markers used was too small to make the Syn4 results definitive, thus the Syn4 results were used to reinforce the interpretation of the Syn0 results. On the basis of these criteria, markers npi262 (bin 1.04), bmc1643 (bin1.08), npi337 (bin 2.06/07), bmc1144 (bin3.03), bnl5.37b (bin 3.06), umc63a (bin 3.09), php20020 (bin 7.06), npi260b (bin 8.03), umc2c (bin8.05), php20075b (bin 9.03), umc130 (bin 10.03), and bmc1185 (bin 10.07) are considered highly likely to be linked to QTL for protein. In addition to being significant for all three progeny types in the Syn0, markers npi262, bmc1643, bnl5.37b, php20020, and umc2c were included in the final Syn0 multiple regression models for all three progeny types (Table 8). For bmc1643, the Syn4 data were also significant (P < 0.01) in the SF analysis for the S1 progenies. Marker npi337 was in the Syn0 final multiple regression model for the S1. Marker umc63a was included in the final Syn0 regression models for the S1 and FR616 progenies and was significant in the Syn4 SF analysis for FR1064. Markers umc130 and bmc1185 were also in the final multiple regression models for FR616 and FR1064 progenies, but not in the S1 progenies. Marker npi260b was included in the list of marker QTL associations highly likely to be linked to QTL for protein, even though it was not significant in all three progeny types in the Syn0 because it was significant in all three progeny types in the Syn4.

With the same criteria as for protein, markers bmc1953 (bin 1.02), npi262 (bin 1.04), umc67a (bin 1.06), bmc1643 (bin 1.08), umc6a (bin 2.03), npi298 (bin 2.08), bnl17.14 (bin2.10), bmc1144 (bin 3.02), umc154 (bin 3.04), bnl5.37b (bin 3.06), php20726 (bin 3.09), umc63a (bin 3.09), php20020 (bin 7.06), bnl9.08 (bin 8.03/04), npi260b (bin 8.03), umc2c (bin 8.05), umc130 (bin 10.03), and bmc1185 (bin 10.07) were considered highly likely to be linked to QTL for starch. Markers bmc1953, npi262, bmc1643, and php20020 were included in the final multiple regression models in the Syn0 for all three progeny types. Marker npi260b was significant in the Syn4 SF analysis for all three progeny types.

Another indication of whether a marker is highly likely to be linked to a QTL is agreement with results of other studies. Eighteen of the RFLP markers used in this study were also used in a study of the cross of generations 76 of IHP and ILP (Goldman et al., 1993). Five of those markers were significantly associated with protein concentration in the Goldman et al. (1993) study. Of those five, three (php10080, sh2, and umc63a) were significantly (P < 0.01) associated with protein in the S1 and both TC progenies in the Syn0 of this study. All three markers are in bin 3.09 (Table 8). Of three markers significantly associated with starch, in the Goldman et al. (1993) study, two (sh2 and umc63a) were also significant for the S1 and both TCs in the Syn0 (Table 9) of our study.

In a study to compare the relationship between a micro wet-milling procedure and the near-infrared reflectance method used in this study, Dijkhuizen et al. (1998) presented results for 26 RFLP markers used in this study. The phenotypic data were the 1995 S1 progeny data. Thus, the Dijkhuizen et al. (1998) results were from a subsample of the data reported herein. Twenty-six RFLP markers were used in both studies. Of the 26 markers in common, 20 showed significant effects in both studies, four were not significant in either, and one was significant in this study and not in Dijkhuizen et al. (1998), while one was significant in the Dijkhuizen et al. (1998) study and not here. Thus results from the one year of data were consistent with the total set of two year's data reported in this study.

Agreement between Traits
The same markers are expected to be associated with significant effects for both starch and protein because of the high negative correlations between starch and protein. In the Syn0, 23 markers were significant at the 0.01 level for both traits for all three progeny types (Tables 8 and 9). To select for increased starch and protein simultaneously, the signs of the additive effects should be in the same direction. As expected, based on the genetic correlations between traits, most of the signs of the additive effects for protein were in the opposite direction from the signs for starch and this result was consistent across all three types of progeny. However, for umc67a in bin 1.06, npi238 in bin 1.11, and bnl17.14 in bin 2.10, the signs were the same for the additive effects for both protein and starch for all three progeny types. Thus, these regions may be of interest for developing high-protein, high-starch genotypes.


    CONCLUSIONS
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
The agreement between S1 progeny results and results from the FR616 and FR1064 TC progeny data in the Syn0 for number of markers identified as significantly associated with QTL for starch or protein was high. This agreement was likely because of the reduced error associated with progeny means in the TC data because the magnitude of the additive effects in the TCs are expected to be half of those in the S1's, based on an additive model. This result, in conjunction with the agreement of the reduction of genetic variance variances in the TCs with theoretical expectations, agrees with the results of Dudley (1994) in suggesting the lack of dominance variance for protein. The agreement between S1 and TC results is particularly encouraging in that the S1 progenies were grown in a different set of years from the TC progeny. The agreement between results for the two testers is also encouraging in that the two testers were from different heterotic groups. Even though the genetic correlations between starch and protein were highly negative, three markers were identified for which the signs of the additive effects were the same in both the Syn0 and Syn4, suggesting it should be possible to make at least small increases in starch and protein simultaneously.

To the authors' knowledge, this is the first report of the effects of random mating on identification of QTL in plants. From a cross for which the parents differed drastically for protein and starch, large numbers of QTL should have been segregating. Use of such a cross assured the presence of coupling phase linkages for protein and, given the high negative correlation between protein and starch, for starch as well. Random mating drastically reduced the number of significant marker–QTL associations. This result was expected because of the breakup of marker QTL linkages and of linkages between QTL. Accompanying the reduction in number of significant effects was a drastic reduction in the proportion of the genetic variability accounted for by multiple regression models. This reduction is likely because of the small number of markers available. This observation points to the value of emerging common resource populations such as the Intermated B73 x Mo17 (IBM) population that has more than 1200 molecular markers. This type of new resource with extensive marker coverage of the genome will enable more precise identification of QTL in randomly mated populations.


    NOTES
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
This research was supported by the Illinois Agricultural Experiment, a grant from the Consortium for Plant Biotechnology Research, Monsanto Corp., Cargill, Inc., and Golden Harvest Seeds.

Received for publication September 2, 2003.


    REFERENCES
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 




This article has been cited by other articles:


Home page
Crop Sci.Home page
J.W. Dudley
From Means to QTL: The Illinois Long-Term Selection Experiment as a Case Study in Quantitative Genetics
Crop Sci., December 18, 2007; 47(Supplement_3): S-20 - S-31.
[Abstract] [Full Text] [PDF]


Home page
Crop Sci.Home page
J. W. Dudley, D. Clark, T. R. Rocheford, and J. R. LeDeaux
Genetic Analysis of Corn Kernel Chemical Composition in the Random Mated 7 Generation of the Cross of Generations 70 of IHP x ILP
Crop Sci., January 22, 2007; 47(1): 45 - 57.
[Abstract] [Full Text] [PDF]


Home page
Crop Sci.Home page
D. Clark, J. W. Dudley, T. R. Rocheford, and J. R. LeDeaux
Genetic Analysis of Corn Kernel Chemical Composition in the Random Mated 10 Generation of the Cross of Generations 70 of IHO x ILO
Crop Sci., February 24, 2006; 46(2): 807 - 819.
[Abstract] [Full Text] [PDF]


Home page
Crop Sci.Home page
B. A. Barrett, I. J. Baird, and D. R. Woodfield
A QTL Analysis of White Clover Seed Production
Crop Sci., August 1, 2005; 45(5): 1844 - 1850.
[Abstract] [Full Text] [PDF]


Home page
Plant CellHome page
L. M. Wilson, S. R. Whitt, A. M. Ibanez, T. R. Rocheford, M. M. Goodman, and E. S. Buckler IV
Dissection of Maize Kernel Composition and Starch Production by Candidate Gene Association
PLANT CELL, October 1, 2004; 16(10): 2719 - 2733.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
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 Similar articles in ISI Web of Science
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 (6)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Dudley, J. W.
Right arrow Articles by Rocheford, T. R.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Dudley, J. W.
Right arrow Articles by Rocheford, T. R.
Agricola
Right arrow Articles by Dudley, J. W.
Right arrow Articles by Rocheford, T. R.
Related Collections
Right arrow Crop Genetics
Right arrow Maize
Right arrow Biometrics


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