Departments of Agronomy (Avis and Taliaferro) and Statistics (Holbert), Oklahoma State University, Stillwater 74078
Gentoype × environment (GE) interactions were studied using forage yield data from 11 bermudagrass [Cynodon dactylon (L.) Pers.] genotypes agronomically evaluated for 3 years (1974 through 1976) at each of four Oklahoma locations. Mean squares, estimated variance components, and regression analysis were used to draw conclusions.
The first-order genotype × year (GY) interaction had a statistically nonsignificant (P > .05) mean square and a relatively small variance component (a2GY^= 0.11). The mean square and the variance component attributable to the genotype × location (GL) interaction effect were respectively, significant (P < .05), and of intermediate magnitude (a2GL^= 0.33). The genotype × replication in location GR(L) and the second-order genotype × location × year (GLY) interactions had highly significant (P < .01 ) mean squares and large variance components (a2GR(L)^= 0.73, a2GLY^= 0.64).
The results emphasized the considerable magnitude of GE interactions and the necessity for multiple environment testing through time and space to characterize relative genotypic differences. The nonsignificant mean square and the relatively small variance component attributable to the effect of years indicates that the duration of testing is less important than the number of locations at which the testing is conducted.
Regression analysis gave some indication of, but did not elucidate, relative genotypic yield stability across environments. It was very useful, however, in ascertaining specific instances of genotypic instability.
The failure of two or more genotypes (cultivars, lines, strains, etc.) to respond similarly to an environment (genotype × environment interaction) complicates their agronomic evaluation and characterization with respect to relative performance potential and usefulness. Such genotype × environment (GE) interactions prevent the extrapolation of results of agronomic evaluations from one location to another, thus requiring expensive trials at multiple locations. Knowledge of the magnitudes of GE interactions and of the various sources of variation in GE interactions, such as the genotype × year or the genotype × location, is useful in decision-making regarding the amount of testing necessary through time (years) and space (locations) to accurately measure relative genotypic performance (1-5). Additionally, regression analyses methods have been devised to measure the stability of performance, or general adaptation, of individual genotypes across varying environments in an effort to provide meaningful biological explanation of GE interactions (6, 7).
The Oklahoma Agricultural Experiment Station supports an extensive bermudagrass [Cynodon dactylon (L.) Pers.] breeding improvement program having as its principal goal the transfer of higher nutritive value from unadapted tropical species to adapted types with high yield potential. In such a program selected genotypes must be thoroughly compared to existing cultivars and their geographic areas of adaption should be described prior to release for public use. In Oklahoma, bermudagrass is used over a wide range of climatic and edaphic conditions ranging from the semi-arid western regions to the more humid regions of the east. The temperature differential across the state is of major importance since bermudagrass genotypes often differ in winterhardiness and it is primarily this trait which dictates the latitude to which they can be successfully utilized.
The studies reported herein assess the importance of GE interactions in bermuda-
grass forage yields in Oklahoma, estimate the relative magnitudes of the components of GE interactions, and draw conclusions regarding evaluation procedures in the state.
The data used in the yield analyses are from four cultivars and seven experimental genotypes of bermudagrasses grown for 3 years (1974-1976) at each of four Oklahoma locations near Altus, Chickasha, Mangum, and Muskogee. The experimental strains were from the Oklahoma Agricultural Experiment Station bermudagrass breeding project and carried the designations S-13, S-29, S-54, S-24, SS-16, SS-27, and S-78. The cultivars were 'Midland', 'Hardie', 'Oklan', and 'Alicia'. All represent single plants that have been clonally propagated and will henceforth be referred to as genotypes except where a distinction is made between cultivars and experimental lines.
The test locations were selected to sample the climatic and edaphic conditions likely to be encountered in growing bermudagrass throughout the state. Mangum, in the southwestern portion of Oklahoma, is characterized by semiarid conditions and has a Meno sandy loam soil, a member of the loamy, mixed, thermic Aquic Arenic Haplustalfs. Altus, also in the southwest, is likewise characterized by semi-arid conditions but some supplemental irrigation water was applied during dry periods. That location has a complex soil series of Tillman clay loam, a member of the fine, mixed, thermic Typic Paleustolls and Hollister clay loams, a member of the fine, mixed, thermic Pachic Paleustolls. Chickasha, in central Oklahoma, represents a moderate rainfall area and has a Reinach silt loam soil, a member of the coarse-silty, mixed, thermic Pachic Haplustolls. Muskogee, in eastern Oklahoma, represents a relatively high rainfall area and has a Taloka silt loam soil, a member of the fine, mixed, thermic Mollic Albaqualfs.
All test plats were planted in the spring of 1973. The field plot design at each location was a randomized complete block with four replications. The grass in individual plots was started by transplanting ten greenhouse-grown plants approximately 0.62 m apart in a 6-m long row. The laterally spreading, sod-forming grasses were restricted to plots measuring approximately 1.83 × 6.10 m. All plots were fully covered by grass by the end of the 1973 growing season. During each year of the 3-year test period all plots received two applications of nitrogen fertilizer. The first application was made in April and the second in late June or July. The yearly total amounts of actual elemental nitrogen applied to the respective tests varied between locations from approximately 134 to 179 kg/ha. Forage yields were measured by harvesting and weighing material from a 0.92 × 5.49 m swath from the center of each plot. Samples of the harvested forages were taken for determination of moisture content, and these data were then used to convert yields to a dry matter basis. Yield data within a year were summed over cuts to give total seasonal production figures on which the analyses herein were conducted.
A conventional analysis of variance based on that outlined by Comstock and Moll (2) and modified to fit a perennial crop by Taliaferro et al. (8) was used to separate the components of variance. A random-effects model was used assuming a random sampling of genotypes, locations, years, etc. Unbiased estimates of the genetic and environmental components of variance were obtained by equating the expected mean squares with those calculated from the experiment. Cochran's (9) method for deciding upon the appropriate F-test and calculation of degrees of freedom was used.
Because the analysis of variance indicated the presence of significant GE interactions, a joint regression analysis as outlined by Eberhart and Russell (6) was used to determine whether the interactions were a linear function of the additive environmental component. An environmental index was calculated by taking the means of all genotypes grown at a particular location in a particular year and subtracting from that quantity the mean of all genotypes over all environments. The regression of mean yield on environmental index was obtained for each genotype within and over years. The regression coefficient (b) measures the increase in mean yield of a genotype per unit increase in environmental index. The standard deviation (sd)¯of residuals about the regression line measures how well predicted response agrees with that actually observed and includes GE
interactions. A stable genotype, as described by Eberhart and Russell (6), has a regression coefficient of unity (b = 1.0) and deviations from regression as small as possible (sd¯= 0.0).
Yield, Analyses of Variance, and Estimates of Variance Components. Mean genotypic yields, averaged over locations, differed significantly (P < .01) during each test year and over test years. There was only slight variation in the relative rankings of individual genotypes from year to year (Table 1). The mean yields of geno-\
types at the respective locations, averaged over years, also differed significantly (P < .05), but there was considerably more variation in the relative rankings of individual genotypes from location to location, as compared to that from year to year (Table 2).
Mean squares and estimates of variance components relevant to the study of GE interactions are given in Tables 3 and 4, respectively. Mean squares associated with the GR(L) and GLY interactions were highly significant (P < .01). The GL and GY interactions had respectively, significant (P < .05) and nonsignificant (P > .05) mean squares.
The variance component estimate for a2GY^was relatively small (0.11) compared to estimates for the other components. The estimate for the a2GL^component was of intermediate size (0.33) while the estimates for the a2G,^ a2GR(L)^and a2GLY^components were large (0.62, 0.73 and 0.64, respectively) and were exceeded only by the error component a2e^ (0.99).
Regression coefficients and their standard errors for the 11 genotypes are listed in Table 5. With yields averaged over locations and years, all genotypes had regression coefficients within two standard errors of unity (1.0). For location yields, averaged over years, several genotypes had regression coefficients differing from unity by more than two standard errors.
The standard deviation of residuals about the regression line vary considerably within and over locations (Table 6). The genotypes could be arbitrarily divided into low, intermediate, and high groups based on their standard deviation of residuals. The genotypes S-29, S-54, Hardie, and SS-27 had relatively low standard deviations. Genotypes S-24, S-78, Midland, Alicia, and SS-16 had standard deviations of intermediate magnitude while genotypes Oklan and S-13 had relatively high standard deviations. Specific instances of genotypic instability, based on significant yield deviations from a 95% confidence interval set about the
regression line for each genotype, could be detected at individual locations in individual years.
The highly significant mean square and the large variance component estimate for the GLY interaction indicate the necessity for multiple environmental testing if the relative performance of bermudagrass genotypes is to be accurately assessed for a large geographic region. This gives no clues, however, as to how such testing should be distributed through time and space. The magnitudes of the mean squares and variance components for the GY and GL interaction effects suggest that years have considerably less effect than locations on relative yield performance. However, bermudagrass genotypes differ in their ability to persistently maintain good stands and yield potential. Some genotypes, particularly those with intermediate levels of winterhardiness, may perform satisfactorily for a period of time, perhaps a year or two, then begin to deteriorate in stand and yield. This is due apparently to an interaction of the stress brought upon the genotypes by multiple seasonal harvests and natural environmental hazards including freeze damage, drought, and attacking insect or disease pests or both. Genotypes have a threshold tolerance level to these stresses which when surpassed results in stand decline and yield reduction. Consequently, it is necessary to test genotypes for a minimum of 2 years, with any additional years of testing giving greater reliability to the results. However, if a choice has to be made between increasing the number of years of testing beyond 2 (or perhaps 3) and increasing the number of testing locations, our data indicate the latter to be the better choice.
The highly significant mean square and the large variance component for the GR (L) interaction mean that the magnitude of the yield differences between genotypes varied from replication to replication in the individual tests. In the testing of perennial crops such as bermudagrass, replications remain static for the duration of the test and there is consequently some correlation of individual plot yields from year to year. Such correlations tend to cause underestimation of true GE interactions and an overestimation of the true genotypic component. The large GR(L) interaction and error effects are not unusual in tests of this type but are nonetheless disturbing by virtue of their graphic illustration of the amount of variability generated at the submacro- to microenvironmental levels. Since clonal propagules of single genotypes are used to establish plot replicates, the yield variability from such plots should be due entirely to environment or experimental errors (measurement error) or both. The relative size of these effects should conceivably be reduced by using more uniform test environments or more replications or both. From a practical standpoint, however, these measures are frequently either not possible or not feasible. The use of more elaborate experimental designs, e.g. lattices and latin squares, may sometimes prove beneficial in reducing experimental error.
The regression analysis did not elucidate the relative performance stability over environments of genotypes but did pinpoint specific instances of yield instability. On the basis of the standard deviations of residuals about the regression line, the genotypes S-54 and S-29 were the most stable. Though their performance was relatively poor, they were the most consistent over locations and years. Several genotypes had regression coefficients larger than one though they did not differ from one by more than two standard errors. Joppa et al. (7) point out that regression coefficients larger than one are indicative of either better than average performance in high-yield environments or worse than average performance in lowyield environments.
From the data presented, we conclude that the large GE interactions for bermudagrass forage yields in Oklahoma mandate multiple environmental testing to elucidate true genotypic yield performance across the state. These tests must be conducted at several locations representing as many of the major climatic and edaphic regions as possible. We believe that the duration of the tests should be at least 3 years including the year of establishment.
1. E. L. BREESE, Heredity 24: 27-44 (1969).
2. R. E. COMSTOCK, and R. H. MOLL. In W. D. HANSON and H. F. ROBINSON (ed.), Statistical genetics and plant breeding. Natl. Acad. Sci.-Natl. Res. Coun., Washington, D. C., 1963, pp. 164-196.
3. J. HILL, J. Agric. Sci., Camb. 85: 477-493 (1975).
4. R. R. KARWASRA, H. R. TADAV, and R. S. PARODA, Euphytica 24: 261-267 (1975).
5. J. C. MURRAY and L. M. VERHALEN, Crop Sci. 10: 197-199 (1970).
6. S. A. EBERHART and W. A. RUSSELL, Crop Sci. 6: 36-40 (1966).
7. L. R. JOPPA, K. L. LEBSOCK, and R. H. BUSCH, Crop Sci. 11: 238-241 (1971).
8. C. M. TALIAFERRO, C. E. DENMAN, R. D. MORRISON, and D. HOLBERT, Crop Sci. 13: 619-622 (1973).
9. W. G. COCHRAN, Biometrics 7: 17-32 (1951).