Crop Science Journal of Natural Resources and Life Sciences Education
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Published online 1 January 2005
Published in Crop Sci 45:18-26 (2005)
© 2005 Crop Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA
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 Related articles in Crop Science
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 (7)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Casanoves, F.
Right arrow Articles by Balzarini, M.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Casanoves, F.
Right arrow Articles by Balzarini, M.
Agricola
Right arrow Articles by Casanoves, F.
Right arrow Articles by Balzarini, M.
Related Collections
Right arrow Statistics
Right arrow Biometrics

Evaluation of Multienvironment Trials of Peanut Cultivars

F. Casanovesa,*, J. Baldessarib and M. Balzarinic

a Centro Agronómico Tropical de Investigación y Enseñanza, 7170 Turrialba, Costa Rica
b EEA-Manfredi, Instituto Nacional de Tecnología Agropecuaria, Manfredi, Córdoba, Argentina
c Facultad de Ciencias Agropecuarias, Universidad Nacional de Córdoba, cc 509, (5000) Córdoba, Argentina

* Corresponding author (casanoves{at}catie.ac.cr)


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Multienvironment yield trials (MET) for advanced peanut lines are conducted each year at the EEA-Manfredi Peanut Breeding Program, the main INTA program for developing new peanut (Arachis hypogaea L.) cultivars for cultivation in the Argentinean crop area. The main objective of this work was the simultaneous analysis of several multienvironment yield tests first to identify superior cultivars for the peanut crop area in Argentina, and second to investigate if different megaenvironments exist. The simultaneous evaluation of several years of MET provides information that allows researchers to better guide breeding strategies. We analyze a 6-yr series of grain yield data from MET, involving 18 genotypes and five test locations using six by-year analyses of complete yield data sets and an Additive Main Effect and Multiplicative Interaction (AMMI) mixed model analysis combining all 6 yr of MET. AMMI models in a mixed model framework were used for exploring genotype–environment (GE) interaction since the lists of genotypes annually tested in multienvironment trials vary from year to year since new genotypes are introduced every year and others are withdrawn. The results allowed us to identify mf484 and mf505 as superior cultivars and confirm the existence of a unique megaenvironment for identifying high yield cultivars in the peanut crop area of Argentina. The mixed model approach of MET data was successfully implemented to analyze highly unbalanced GE data sets.

Abbreviations: AIC, Akaike Information Criterion • AMMI, additive main effect and multiplicative interaction • BIC, Schwarz Bayesian Criteria • COI, Cross-over interaction • E, environment main effect • EEA, Estación Experimental Agropecuaria • FA, Factor Analytic • G, genotypic main effect • GE, genotype by environment interaction effect • GGE, G plus GE • GL, genotype x location interaction effect • INTA, Instituto Nacional de Tecnología Agropecuaria • L, location main effect • MET, multienvironment trials • PBP, peanut breeding program • PC, principal component(s) • SREG, sites regression


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
THE GERMPLASM EVALUATION is a crucial activity in plant breeding (Stroup, 2000). MET are commonly conducted annually to obtain information that supports recommendations of superior cultivars for cultivation. MET are used to evaluate several genotypes in multiple environments (locations and/or years), and they are essential because of the presence of GE interaction, i.e., differential genotypic responses in different environments. The main objective in the evaluation of a series of MET is to identify superior cultivars for a target region and to determine if this region can be subdivided into different megaenvironments to better guide breeding strategies (Kang, 2002). Important concepts such as ecological regions, ecotypes, megaenvironments, specific adaptations, and stability originated from the GE interaction analysis (Yan and Hunt, 2002). The identification of megaenvironments is associated with the exploration of the annually repeatable GE interaction patterns. For a particular megaenvironment, genotypes are evaluated on the basis of mean yield and stability of yields across environments.

The EEA-Manfredi, INTA, Argentina, conducts MET in a Peanut Breeding Program (PBP-INTA) in different environments (locations and years). The selected locations of MET are representative of the environmental characteristics of the northern, central and the southern zones of the Argentinean peanut crop area (Córdoba province) that extends from 32° to 33°50' S latitude and from 63° to 64°35' W longitude. This region produces 95% of the Argentinean peanuts since it shows high homogeneity of climate and soil characteristics appropriated for the cropping. The genotypes evaluated during all 6 yr of MET vary from year to year because new genotypes are introduced every year and others are withdrawn. Therefore, the MET databases through the years are incomplete, i.e., all of the genotypes are not present in all of the combinations of locations and years.

The unbalanced data sets obtained from compiling several years of MET data was handled by modeling the response under the mixed model theory (Balzarini, 2000). This modeling strategy allows inferring the genotype performance across environments even when the GE tables are incomplete. Additionally, the stability and interaction measures can be obtained as certain mixed model parameters (Piepho, 1998; Balzarini, 2002).

In this study, the genotype effects were considered fixed because they refer to genotype sets in advanced selection stages. The environmental effects were considered random, except for an annual MET analysis where differences among environments refer to differences among locations and they could be associated with predictable differences among locations (Allard and Bradshaw, 1964). For the combined MET analysis (through 6 yr), the environments were defined as combinations of location and year effects, and treated as random. For both by-year and all-year MET analyses, the GE interaction was modeled by an AMMI model, but with the fixed or the mixed approach to AMMI according to the balancedness of the data sets. The evaluation of several years of multienvironment yield trials provides useful information to help researchers to better guide breeding strategies.

Up until now, studies have not been done to investigate better cultivars and the possible existence of different megaenvironments using yield data from several MET. The main objective of this work was the simultaneous analysis of several multienvironment yield tests to identify superior cultivars for the peanut crop area in Argentine, and secondarily to determine if different megaenvironments exist.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Database
The information sources used include the MET of PBP-INTA conducted from year 1996–1997 to 2001–2002, which correspond to the first 6 yr of multienvironmental evaluation in the program. The trials were conducted in 10 sites, some of which because of their proximity and similarity in soil and weather conditions were regarded as the same location for the purposes of the analysis (Table 1). Altogether 18 genotypes, four short cycle and 14 long cycle, were evaluated (Table 2). The set of evaluated genotypes was the same for all the locations within a year, but the genotype list in the study period varied from year to year. The MET of the PBP-INTA for advanced lines began in 1996 with 11 genotypes without structural changes (genotypes and locations) in the first 2 yr. During the next 4 yr, seven new genotypes were incorporated, six genotypes were dropped and new locations were added. Florman, mf484, mf485, mf487, and mf489 were the only cultivars evaluated during all 6 yr. The genotype x location x year database is highly unbalanced. Table 2 contains a list of the evaluated genotypes. In all of the sites, the trials were laid out as a randomized complete block design, with four replications. The plots were made of two 10-m-long furrows with 70 cm between furrows. Recommended seeding rates (15 seeds/m2) and cultural practices were used at all locations. Each plot was manually harvested. The analyzed yield values correspond to kilograms of peanuts per plot at constant moisture (80 g kg–1). In Table 3, the main events that characterized the years involved in this evaluation of MET data are shown.


View this table:
[in this window]
[in a new window]

 
Table 1. Participating locations in multienvironment trials of the INTA-EEA-Manfredi Peanut Breeding Program of Argentina from 1996–1997 to 2001–2002.

 

View this table:
[in this window]
[in a new window]

 
Table 2. Genotypes evaluated in multienvironment trials in the EEA-Manfredi Peanut Breeding Program, INTA, from 1996–1997 to 2001–2002.

 

View this table:
[in this window]
[in a new window]

 
Table 3. Main characteristics of the agricultural years 1996–1997 to 2001–2002 in the EEA-Manfredi Peanut Breeding Program, INTA.

 
Statistical Analysis of MET
Considering the locations as the basic megaenvironment units and that the genotype x year and location x year data tables are highly unbalanced, genotype x location (GL) interaction was explored via individual yearly analyses (complete data tables). In this first modeling attempt, we run an analysis of variance model which included fixed effects for location (L), genotypes (G), and genotype x location interaction (GL). The adjusted model was:

Formula 1[1]
where Yijk is the yield of Genotype i, in the Location j, Block k; µ is the overall mean; Lj is the effect of Location j with j = 1,...,s; B(L)k(j) is the effect of Block k within Location j with k = 1,...,n; Gi is the effect of Genotype i with i = 1,.. g; GL(ij) is the interaction effect between Genotype i with Location j and {epsilon}ijk is the random error term associated with observation Yijk. Two models were developed, one model assumed homogenous residual variances and the other allowed heterogeneous residual variances across locations. The latter model was evaluated because the trials conducted in different locations could have different precision (residual variance).

For each year, the comparison of genotype performance across environments (broad inference) was based on the least squares means. Since a fixed effects model was considered, the standard errors of mean differences only depended on the residual variance. The narrow inferences (environment-specific inference) about a genotype performance were obtained from genotype means in each environment. Symmetrical biplots were used (Gabriel, 1971) to analyze GL interaction patterns associated with the AMMI model (Gauch, 1988) and with the site regression model (SREG) (Cornelius et al., 1996). The biplot algorithm is based on the single value decomposition of a residual matrix. The AMMI model biplots are constructed from the residual matrix of the additive model, i.e., Eq. [1] without the GL interaction effect. The SREG model biplots, known as GGE biplot (Yan et al., 2000), are constructed from the residual matrix corresponding to an adjustment Eq. [1], which omits G and GL.

Scatter plots were generated for each year using mean yield (x axis) and stability measurement (y axis) for each genotype. Four stability statistics were calculated: (i) the CV (Francis and Kannenberg, 1978) because it is the statistic traditionally used in the PBP-INTA, (ii) the stability variances of Shukla (1972), (iii) the first principal component (PC1) of the AMMI model analysis, and (iv) the first and/or second principal component (PC1 and PC2) of the SREG model.

Combined MET analysis
A model with environment, genotype, and genotype x environment effects was used. Environmental effects were defined as the combination of the location and years. All the effects used in the model were considered random, except for genotype. The equation for the response of Genotype i in Block k in Environment j is:

Formula 2[2]
where Yijk is the response of Genotype i, in Environment j and Block k, µ is the overall mean, Gi is the fixed effect of Genotype i, Ej is the random effect of the Environment j, GEij is the random effect of the interaction of Genotype i and Environment j, B(E)k(j) is the effect of Block k within the Environment j, {epsilon}ijk is the random error term associated with Observation Yijk.

The overall genotype performance (broad inference) was based on the genotype means, but the standard errors of the mean differences depended on the mixed model used for the variance of genotype means that included variance components associated with the GE interaction. Thus, the statistical comparison of the genotype performance, both yield differences and differences in yield stability were considered. The variance components were estimated via restricted maximum likelihood (REML) using SAS (SAS Institute, 1997) PROC MIXED.

Different models were developed for the variance and covariance matrix of the random interaction terms for a given environment: (i) traditional mixed model, i.e., Eq. [2] assuming homoscedastic variances for the random GE terms and (ii) mixed stability variance model, i.e., model analogous to Eq. [2] considering heteroscedastic variances across genotypes for the GE interaction terms, and 3) AMMI mixed model. The AMMI mixed model is expressed as

Formula 2
where yijk is the response of Genotype i, in Environment j and Block k, µ is the overall mean, Gi is the fixed effect of Genotype i, Ej is the random effect of the Environment j, and Formula 2{lambda}m{gamma}im{alpha}mj models the GE interaction as the sum of M multiplicative components that explain the GE interaction in M orthogonal directions [in each direction, fixed genotype scores ({gamma}im) and random environment ({alpha}mj) intervene]; {lambda}m is a weight factor associated with the mth direction; and {epsilon}ijk represents the portion of the GEij interaction not explained by the model plus the random error term associated with yijk.

Using the homocedastic factor analytic (FA1) structure (Jennrich and Schluchter, 1986) to model the variance and covariance matrix of the GE interaction terms within environment, the estimated covariance parameters can be used as genotype scores (genotypic sensitivity) to explain the GE interaction. The general factor analytic (FA) model is QQ' + D, where Q is a scores matrix of qs and D is a diagonal matrix with the possibility of nonnegative parameters on the diagonal. If matrix D is omitted, that is to say D = 0, then the model is denoted as FA0. If matrix D has all its diagonal elements equal (D = {sigma}2I), then the model is denoted as FA1. In this work structures FA1, of order 1, 2, and 3 corresponding to AMMI models with M = 1, M = 2 and M = 3 multiplicative terms, respectively, were used. The selection of the most suitable parameterization to model the GE interaction was made through maximum likelihood ratio tests, the Akaike information criterion (AIC) and the Schwarz Bayesian criteria (BIC) (Littell et al., 1996). The graphical biplots associated with a AMMI mixed model (AMMI biplot) were obtained by graphing the covariance (standardized) parameters associated with each genotype in each multiplicative term versus the empirical best linear unbiased estimator (BLUP) of the environmental effect on the same term using the statistical software Info-Gen (Info-Gen, 2003).


    RESULTS AND DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Analysis of MET by Year
The model with heterogeneous residual variances across locations fit better than the model with homogenous variances in 4 of the 6 yr of MET data analyzed (Table 4). The Akaike information criteria (AIC) and Schwarz (BIC) (in which a greater value implies a better fit) coincided with the likelihood ratio test in identifying the years with different precision from the trials conducted in different locations.


View this table:
[in this window]
[in a new window]

 
Table 4. Fitting criteria for the fixed effects multienvironment trials model with and without heterogeneous residual variance across the locations

 
Table 5 contains the analysis of variance for each year of MET data processed with the models with better fits as it is suggested in Table 4. In all the years analyzed, except 1996–1997, the GL interaction was statistically significant (P < 0.01). The PBP-INTA MET began in 1996, the year in which the residual variance was highest. This low accuracy of the tests could mask the presence of GL interaction. The relative magnitude of the GL interaction with respect to the variability explained by (G+L+GL) for each year is given in Table 5. The variation caused by the GL interaction was smaller than the variation among genotypes in 5 of the 6 yr of the MET analyzed. The variation between locations was always the most important one, explaining between 74.9 and a 96.0% of the total variation, which justifies the selection of the biplots based on the SREG model for MET analysis (Yan et al., 2000).


View this table:
[in this window]
[in a new window]

 
Table 5. Variability of yields between genotypes (G), locations (L), block between site [B(L)], and G x S interaction for PBP-INTA multienvironment trials, for the 1996–1997 to 2001–2002 years.

 
The first two principal components (PC1 and PC2) obtained by singular value decomposition of the centered data (SREG model) explained more than 85% of the total variability caused by (G+GL) for all the years except in 1998–1999 and 2000–2001 where these percentages were 79 and 83%, respectively. Figure 1 represents the GGE biplots obtained for each year of MET.


Figure 1
View larger version (39K):
[in this window]
[in a new window]

 
Fig. 1. GGE biplots based on 6 yr of PBP-INTA multienvironment trials. The dark points represent genotypes and the light points represent locations.

 
In the year 1996–1997 (Fig. 1.a), there was a high correlation of the PC1 (r = 0.98, P < 0.0001) and low correlation of the PC2 (r = 0.15, P = 0.641) with yields. The extreme genotypes on the PC1 axis, mf478 (positive PC1) and Florman (negative PC1), do not show a crossover interaction (COI); these were genotypes of higher and lower rankings, respectively. The GL interaction was not significant and all the locations behaved as a unique megaenvironment where some advantages were demonstrated for the short-cycle cultivars over those of long cycle. This is possibly due to the low precipitation registered during this year. The mf487 and mf489 genotypes (positive PC2) and the mf480 genotype (negative PC2) were the only ones that showed a differential response, COI across environments.

In the year 1997–1998 (Fig. 1.b), the yields of the long-cycle check cultivars were relatively high and homogenous across environments, showing advantages with respect to the other genotypes in Location 2 where the harvest was late. The mf485 and mf484 genotypes showed COI, with a relatively superior performance in Location 4 and with a lower ranking in Location 1. The mf480 and mf457 genotypes showed the opposite behavior with relatively high performances in Location 1, where the crop was subjected to stress by a Sclerotinia minor Jaegger infection. In this year, the yield did not correlate with the PC1 (r = 0.062; P = 0.857). Therefore, the genotypes with greater projections of the PC1 (mf485, mf484, mf457 and mf480) showed COI, whereas the Tegua, Florman, and mf480 genotypes, with greater positive projections of the PC2, yielded relatively better than the rest during this year.

The year 1997–1998 was the only year of MET data where the GL interaction explained a greater percentage of (G+L+GL) than G. Although this result could suggest the possible existence of different megaenvironments in the PBP-INTA target region, it is important to note that the observed pattern was not consistent through the years. Besides that, short- and long-cycle cultivars were evaluated, and the harvests in all of the MET, mainly in Location 2, were delayed during this year.

In the year 1998–1999 (Fig. 1.c), the correlation of the PC1 with yield was 0.89 (P < 0.0001). The genotypes with greater PC1 values were mf505 and mf484 (long cycle), both with the highest ranking in most of the locations, except Location 2. The short-cycle cultivars, mf447 and mf480 (negative PC1), showed relatively poor performances except in Location 2. The precipitation was low during the cultivation cycle and high during the harvest, which could explain the differential performance of the short-cycle cultivars. The differences in the genotype response across environments were less than the differences in the mean genotype response, thus the data did not suggest the presence of megaenvironments. Location 4 was more favorable than Location 2 in this year for genotype mf505. Nevertheless, in the following year (Fig. 1.d), this genotype showed advantages relative to the other genotypes in Location 4 and 2. The remaining genotypes are grouped into a unique megaenvironment where the short cycle mf447, mf480, and manf393 genotypes had the worst performance.

In the year 2000–2001 (Fig. 1.e), the contribution of the GL interaction was approximately equal to the variability among genotypes and the yield correlated significantly (P < 0.05) with PC1 and PC2, which suggested that the projections of PC2 did not necessarily indicate COI. The extreme genotypes in their projections on the PC1, mf508, Florman, and mf487, showed greater COI. The Florman and mf487 rankings were relatively higher in Location 1 than in the rest, which showed a behavior opposite to mf508. The mf484 and mf496 genotypes with opposite PC2 values responded proportionally to the differences between high- and low-ranking locations, respectively.

In the year 2001–2002 (Fig. 1.f), the variability between genotypes was much greater than the GL interaction and the yields correlated significantly with PC1 (r = 0.95, P < 0.0001). The mf484, mf505 and mf489 genotypes had the best performance in Locations 2, 3, and 4, whereas Florman and mf487 showed an advantage in Location 1. The genotypes mf508 and mf496 were the worst performers. The separation of Location 1 with respect to the rest, where Florman and mf487 showed relatively high yields, could be due to factors such as scarce rainfall during the grain filling and fruit loss due to the high precipitation registered at harvest in that location. In Location 3, the precipitation was also low, but the differences between the yields of the genotypes were small since a low yield in all of the genotypes was registered.

Figure 2 represents a comparison of the performances of mf480 (short-cycle genotype) and mf489 genotypes (long-cycle) across three locations using the GGE biplot for year 1997–1998, where the GL contribution was greater than that of G. The comparison was made by connecting the point that represents the mf480 genotype with the point that represents mf489 genotype with a straight line and drawing a perpendicular line that passes through the origin (Yan and Hunt, 2002). The perpendicular line separates the locations into two groups, which shows that mf480 genotype had greater relative yields than mf489 genotype in Location 1. Meanwhile, for the mf489 genotype the favorable location was Location 4 followed by Location 2, although yields at Location 2 were barely above average. If genotypes mf480 (short cycle) and mf485, mf484, mf487, and Florman (long-cycle) are compared this way, one concludes that Locations 4 and 2 are favored by the long-cycle cultivars.


Figure 2
View larger version (29K):
[in this window]
[in a new window]

 
Fig. 2. Comparison of short and long cycle genotype performance in all of the intervening multienvironment trials locations during the 1997-1998 year. The dark points represent genotypes and the light points represent locations.

 
A collective analysis of the six GGE biplots, plus the result in Fig. 2, suggest that the peanut region in which cultivars were tested does not show repeatable megaenvironments for breeding purposes. From the results obtained with the analysis of MET data conducted in 1997–1998 where the COI was significant because of differences in natural infections across locations, experiments specially designed to create different environmental conditions appear as appropriated for breeding with respect to pathogen resistance.

Table 6 contains for each one of the analyzed years of MET data, the mean yield for the genotypes, the AMMI model PC1, and the coefficient of variation (CV) as a stability statistics. Although the simultaneous interpretation of the CV and the yield mean constitute a common strategy of analysis implemented in this breeding program, our results show that use of the CV is associated with recommendations for low performance genotypes. For example, in the year 1996–1997, the mf480 genotype has a low CV, and nevertheless, it is a major contributor to the GL interaction. In Location 1, this genotype had a ranking of 11 whereas in Location 4, it ranked third. Something similar happened for the 1997–1998 data where although the ranking of mf480 changed from 11 to 1 across locations, this genotype showed the lowest CV. On the contrary, the use of the PC1 from the AMMI model as a stability measure allowed for better identification of the cultivars with COI. Genotype mf484 was among the superior ones in 5 of the 6 yr and mf505 genotype in 3 of the 4 yr of MET data.


View this table:
[in this window]
[in a new window]

 
Table 6. Yield means and stability measures for each genotype during the 6 yr (1996–1997 to 2001–2002).

 
Combined Analysis of MET
The environments (combination of years and locations) constituted a source of important variation (90.5% of the total variation). The high variations due to environmental differences is expected in MET conducted through several years (Yan and Kang, 2003). The GE interaction was considerably higher than the variability attributable to G, and more than 50% of GE interaction variation is due to the G x L x Y interaction (Table 7). The GL interaction only represents 23.9% of the GE interaction.


View this table:
[in this window]
[in a new window]

 
Table 7. Variance components for genotypes (G), environments (E) and G x E interaction for PBP-INTA multienvironment trials. Joint analysis of the 1996–1997 to 2001–2002 years.

 
Table 8 includes the statistical values that were used to select between the best mixed models: (i) homogenous variance model for GE, (ii) heterogeneous across genotypes variance model for GE, and (iii) AMMI mixed model for GE. The best model for the joint analysis across the 6 yr of MET data was AMMI mixed model with two multiplicative terms (AMMI [2]).


View this table:
[in this window]
[in a new window]

 
Table 8. Fitting criteria for 3 mixed models adjusted to the PBP-INTA multienvironment trials data. Joint analysis of the 1996–1997 to 2001–2002 years.

 
Figure 3 , obtained from the AMMI mixed model, was used to explore the GE interaction in the joint analysis across 6 yr. The mf480 and Florman genotypes made high contributions to the GE interaction; the first had important COI (Environments 1 and 4). Florman responded proportionally to the differences between Environment 1 and 4 environments in the majority of the years. The dispersion of the environmental scores, although they demonstrate a high interaction between years and locations, also show that Location 1 was one of the greater contributors to the GE interaction, especially in the 1997–1998 year because the mf480 and mf457 genotypes yielded better in this location than in others.


Figure 3
View larger version (15K):
[in this window]
[in a new window]

 
Fig. 3. Biplot based on 6-yr multienvironment trials data for PBP-INTA. Dark points represent genotypes and light points represent environments coded according to the harvest year-location.

 
The GE interaction pattern could be difficult to repeat since groups are not visible between the environmental scores from the same location. The GL interaction is small in relation to the GLY interaction. Therefore, the GE interaction is not expected to be repeatable. The results suggest that this set of locations should be considered a unique megaenvironment for breeding objectives.

Figure 4 is a dispersion graph of the product of the genotype scores in the PC1 and PC2 of the AMMI model (2) versus the yield means for each genotype. It shows that Florman, mf480 and mf457 genotypes had lower stability. The mf484 and mf505 genotypes had the best performance across all environments.


Figure 4
View larger version (9K):
[in this window]
[in a new window]

 
Fig. 4. Yield stability versus yield mean.

 
Yield means across all of the environments and the standard errors to evaluate mean differences, obtained from the AMMI[2] mixed model are shown in Table 9. These errors depend on the amount of information available for each genotype through the collection of analyzed MET and on the contribution of the genotype to the GE interaction (Balzarini, 2000). In addition, the table contains the coefficients of variation (CV) across environments and the stability variances (SV) (Shukla, 1972) obtained as covariance parameters of model 2. The stability variances for genotypes mf480, mf457, mf484, Florman, mf447, and manf393 were significantly different from zero (P < 0.05).


View this table:
[in this window]
[in a new window]

 
Table 9. Genotypic performance across 23 environments.

 
The main differences in yield stability across environments were observed in mf480, mf457, mf484, and Florman genotypes according to stability variance, which was consistent with the results obtained from the mixed AMMI biplot. The coefficient of variation fails to identify unstable genotypes that showed relatively high interaction with environments, e.g., mf480, mf478, manf393, and mf447, which are all short-cycle cultivars. Taking into account variability across environments, the cultivar with the best mean response was mf484.


    CONCLUSIONS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
The by-year analysis of MET data showed that the variability due to the GL interaction in the PBP-INTA is relatively small in relation to the variability among genotypes. Although the GGE biplots obtained for each year of MET data allowed us to identify superior genotypes in some locations, the locations that were favorable in a particular year for some genotypes were not favorable through out the years. In 5 of the 6 yr of MET, the GGE PC1 biplot significantly correlated with yield means demonstrating proportional genotype responses across the locations; on occasion, it implied rank changes in genotype order. Using additional information of registered climatic contingencies in each year could explain the interaction. The results obtained from mixed AMMI models using data from 6 yr confirmed the random nature and the relatively low magnitude of the GE interaction in the PBP-INTA. Both strategies to evaluate MET data, the by-year analyses and the joint analysis indicate the existence of a unique megaenvironment for breeding purposes in the peanut crop area of Argentina. The results suggest that it is not necessary to partition the region into subregions to make cultivar recommendations. It also suggests that instead of increasing the number of locations where MET are conducted, it would be better to redirect the resources available toward implementing more efficient experimental designs.


    ACKNOWLEDGMENTS
 
We thank INTA for providing the data. This work was supported by the National Agency of Science and Technology of Argentina (FONCYT PICT 2000/08-08-302).

Received for publication October 16, 2003.


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


Related articles in Crop Science:

THIS ISSUE IN CROP SCIENCE

Crop Science 2005 45: xi. [Full Text]  



This article has been cited by other articles:


Home page
Crop Sci.Home page
B. Glaz and M. S. Kang
Location Contributions Determined via GGE Biplot Analysis of Multienvironment Sugarcane Genotype-Performance Trials
Crop Sci., May 1, 2008; 48(3): 941 - 950.
[Abstract] [Full Text] [PDF]


Home page
Crop Sci.Home page
W. Putto, A. Patanothai, S. Jogloy, and G. Hoogenboom
Determination of Mega-Environments for Peanut Breeding Using the CSM-CROPGRO-Peanut Model
Crop Sci., May 1, 2008; 48(3): 973 - 982.
[Abstract] [Full Text] [PDF]


Home page
Crop Sci.Home page
W. Yan, M. S. Kang, B. Ma, S. Woods, and P. L. Cornelius
GGE Biplot vs. AMMI Analysis of Genotype-by-Environment Data
Crop Sci., March 1, 2007; 47(2): 643 - 653.
[Abstract] [Full Text] [PDF]


Home page
Crop Sci.Home page
H. G. Gauch Jr.
Statistical Analysis of Yield Trials by AMMI and GGE
Crop Sci., May 18, 2006; 46(4): 1488 - 1500.
[Abstract] [Full Text] [PDF]


Home page
Crop Sci.Home page
F. Casanoves, R. Macchiavelli, and M. Balzarini
Error Variation in Multienvironment Peanut Trials: Within-Trial Spatial Correlation and Between-Trial Heterogeneity
Crop Sci., August 26, 2005; 45(5): 1927 - 1933.
[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 Related articles in Crop Science
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 (7)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Casanoves, F.
Right arrow Articles by Balzarini, M.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Casanoves, F.
Right arrow Articles by Balzarini, M.
Agricola
Right arrow Articles by Casanoves, F.
Right arrow Articles by Balzarini, M.
Related Collections
Right arrow Statistics
Right arrow Biometrics


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