|
|
||||||||
a Advanta Semillas S.A.I.C., Ruta Nac. 33 Km 636, CC 559, (2600) Venado Tuerto, Argentina
b CSIRO Plant Industry, Queensland Bioscience Precinct, 306 Carmody Rd., St. Lucia, Qld 4067, Australia
* Corresponding author (avega{at}waycom.com.ar)
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: AMMI, additive main effects and multiplicative interaction BLUE, best linear unbiased estimate E, environment F, female G, genotype GCA, general combining ability M, male NCII, North Carolina Experiment II PC, principal component PCA, principal component analysis SCA, specific combining ability REML, restricted maximum likelihood SREG, site regression
| INTRODUCTION |
|---|
|
|
|---|
Yan and Hunt (2002) used the site regression (SREG2) model to analyze a diallel mating dataset and identify an ideal tester. The biplot of the first two principal components (PC) of the SREG2 model displays the GCA of lines or testers (depending on the type of centering of the two-way lines x testers dataset) as well as the SCA of the line x tester interaction. In the SREG2 biplot, the ideal tester should be the one that has the longest PC1 vector (high power to discriminate lines) and zero projection on the average tester ordinate (more representative of all testers) (Yan and Hunt, 2002). Narro et al. (2003) pointed out that, since the SREG2 model contains in its multiplicative components both the main effects of the lines as well as the line x tester interaction, the visualization of the SCA of each line with each tester is often not clear and is imprecise. Hence, Narro et al. (2003) proposed the use of the additive main effects and multiplicative interaction (AMMI2) biplot, which uses only the line x tester interaction (SCA), to estimate the interaction parameters and to obtain a precise and clear ranking of the lines for their SCA with each tester.
Both the SREG2 (Yan and Hunt, 2002) and AMMI2 (Narro el al., 2003) biplots of lines x testers display the GCA and SCA of the lines averaged across environments. However, hybrid relative performance is not only a consequence of gamete segregation and recombination but also of the environmental effect where cultivars were evaluated (Narro et al., 2003). The G x E interactions result when there is a change in the relative performance of genotypes when they are tested in different environments and have the potential to influence the nature and magnitude of the selection response achieved by a breeding program (Cooper and DeLacy, 1994; Kang, 1998). Some understanding of the nature of the G x E interactions is needed to accommodate their effects through appropriate selection strategies aimed at exploiting broad and/or specific genotype-adaptation patterns (Basford and Cooper, 1998, and references therein). Similarly, consideration of the impact of G x E interactions on GCA and SCA is central to assess the relative merits of different inbred lines potentially used as testers in a breeding program targeting broad and/or specific adaptation.
The variance components estimated from multienvironment trials can be used to judge the relative magnitude of G and G x E interaction and to predict the response to selection (Cooper and DeLacy, 1994). In hybrid crops, the G component of variance can be partitioned into its components, namely female line (F), male line (M) and F x M interaction, allowing the determination of the relative size of the F x E, M x E and F x M x E interaction components of variance as well.
If F x E, M x E and/or F x M x E interactions explain a large portion of the variance in a half-sib mating design study (i.e., GCAF x E, GCAM x E and/or SCA x E interactions are important), then it would be useful to examine the three-way matrix of F x M x E means. Three-mode PCA (Tucker, 1966; Kroonenberg, 1983), an extension of the standard PCA to handle such three-way datasets, has been used for handling genotypes, environments, and attributes simultaneously (Kroonenberg and Basford, 1989; Basford et al., 1990; Chapman et al., 1997; de la Vega et al., 2002; Bertero et al., 2004), allowing an examination of the relationships between genotypes and attributes associated with specific patterns of environmental variability. In the same way, we propose that three-mode PCA can be used to display the relationships between female and male inbred lines associated with predictable or unpredictable environmental contrasts. If the F, M, F x M interaction and E effects are removed from the three-way dataset and F x E, M x E and F x M x E interactions are retained, then the graphical output of three-mode PCA will display the nature of the pattern-rich GCA x E and SCA x E interactions, complementing the analytical tools demonstrated by Yan and Hunt (2002) and Narro et al. (2003). In this analytical strategy (Strategy 1), SGREG2 (Yan and Hunt, 2002) could be used to display GCAF + GCAM + SCA (i.e., F, M and F x M interaction effects, respectively) and three-mode PCA to display GCAF x E + GCAM x E + SCA x E.
In Strategy 1, two- and three-mode multivariate analyses are used to analyze the sources of variation across environments and their interactions with the environments, respectively. Alternatively (Strategy 2), two-mode PCA on environment-centered and normalized (Fox and Rosielle, 1982; Cooper and DeLacy, 1994) two-way F (M) x E tables for hybrid yield can be used to display GCAF + GCAF x E + GCAM + GCAM x E effects, while three-mode PCA is used to display the SCA + SCA x E interaction effect. In this case, the F, M, E, F x E, and M x E interaction effects should be removed from the three-way dataset and F x M and F x M x E interactions should be retained. Together, the two strategies account for all components of the G and G x E interaction variances in a half-sib mating design analysis.
The main objective of this paper is to critically examine the use of two-mode PCA on environment-centered and normalized F (M) x E tables and three-mode PCA on F x M x E tables to graphically display the pattern rich variability associated with GCA, SCA and their interactions with the trial environments in a half-sib mating design study. The proposed methods (including Strategies 1 and 2) are evaluated together with the across-environments SREG2 (Yan and Hunt, 2002) and AMMI2 (Narro et al., 2003) biplots on the basis of their explanation of the total G + G x E variation of the system under study. The paper describes the methods completely and compares the methods for their suitability to detect elite tester lines for a breeding program that aims to target both broad and specific (i.e., regional) adaptation. Although the case study consists in a sunflower mating design, the methods described, as well as those proposed by Yan and Hunt (2002) and Narro et al. (2003), are applicable to any hybrid crop.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
|
Planting took place within the normal sowing window at each location (i.e., about October to mid-November in central and southern locations and about August to mid-September in the northern locations). In each environment, a 4 x 4 square lattice design with three replicates was used to test the 16 hybrids. The trials were over-planted and thinned to 47, 600 plants ha1. A plot size of four rows x 6 m and interrow spacing of 0.70 m was used. All trials were rain-fed and nutrient deficiencies were prevented with fertilization when necessary. Weeds and insect pests were controlled chemically. Plot data of grain yield were determined by hand harvesting of 3.99 m2 (one central row, discarding the border plants). Grain oil concentration was determined on a 10 g oven-dried achene sample by nuclear magnetic resonance (Granlund and Zimmerman, 1975). Oil yield was calculated as the product of grain yield and grain oil concentration.
Analysis of Variance Components
A mixed model analysis was conducted to examine partitions of the G and G x E interaction components of variance for grain yield, oil content, and oil yield. The phenotypic observation yijkmn on hybrid derived from the cross between female inbred line i and male inbred line j in incomplete block n of replicate m of environment k was modeled as:
![]() | [1] |
NID(0,
2r); (b/r/e)kmn the random effect of the incomplete block n nested within the replicate m of environment k and is
NID(0,
2b); fi the random effect of female inbred line i and is
NID(0,
2f); mj the random effect of the male inbred line j and is
NID(0,
2m); (fm)ij the random effect of the interaction between female i and male j and is
NID(0,
2fm); (fe)ik the random effect of the interaction between the female i and environment k and is
NID(0,
2fe); (me)jk the random effect of the interaction between male j and environment k and is
NID(0,
2me); (fme)ijk the random interaction effect for female i, male j, and environment k and is
NID(0,
2fme); and
ijkm is the random residual effect for the combination between female i and male j in replicate m of environment k (experimental error) and is
NID(0,
2
). Restricted maximum likelihood (REML) methods (Patterson and Thompson, 1975) in GenStat v8.1 was used to estimate the variance components of the random terms in the model and their standard errors for grain yield, oil content, and oil yield. Given that oil yield is the most important economic trait for sunflower, we used this attribute to compare different multivariate approaches aimed at understanding patterns of GCA, SCA, and their interactions with the trial environments. The best linear unbiased predictors (BLUPs) (Robinson, 1991) of the F x M, F x E, M x E, and F x M x E means for oil yield were calculated from REML by the mixed model of Eq. [1] and used further to construct the matrices for two- and three-mode multivariate analyses.
Two-Mode PCA-Derived Biplots of GCA and SCA across Environments
The 4 x 4 (F x M) two-way table of across-environments oil yield BLUPs was either female- and male-centered (i.e., centered by subtraction of female (male) mean and normalized by division of the remainder by the within-female (male) standard deviation (Fox and Rosielle, 1982; Cooper and DeLacy, 1994; Yan and Hunt, 2002) and double-centered [i.e., centered by subtracting the across-females male means and the across-males female means, and adding the overall mean, as in the AMMI model (residuals from additivity; Gauch, 1988; Narro et al., 2003)] to graphically display the across-environments GCAF + SCA, GCAM + SCA, and SCA, respectively (Yan and Hunt, 2002; Narro et al., 2003). The PCs of the squared Euclidean distance matrices of oil yield were estimated by a singular value decomposition procedure and biplots of the first two PCs were constructed from these analyses (Gabriel, 1971) by GenStat 8.1.
Two-Mode PCA-Derived Biplots of GCA + GCA x E Interaction
The two-way tables of F x E (4 x 11) and M x E (4 x 11) oil yield BLUPs were centered by subtraction of environment mean from columns and normalized by division of the remainder by the within-environment standard deviation (Fox and Rosielle, 1982; Cooper and DeLacy, 1994). The PCs of the squared Euclidean distance matrices of oil yield were estimated by a singular value decomposition procedure and biplots of the first two PCs were constructed from these analyses (Gabriel, 1971).
Three-Mode PCA
This procedure derives components, i.e., linear combinations of the levels of the modes, for each of the three modes (say, P, Q, and R components for female lines, male lines, and environments, respectively). It can be assumed that these components together contain the only relevant systematic variation of the three-way array dataset. The components of the three modes can be labeled on the basis of the patterns shown by the levels with high loadings on such components. In this model, each mode is allowed to have a different number of components. The number of components for each mode needs to be simultaneously determined for all modes. Therefore, several solutions have to be inspected to come to an adequate description of a dataset (Kroonenberg, 1983, chap. 2; Kroonenberg and Basford, 1989).
A three-way array of order P by Q by R (the core array) contains the weights assigned to each of the combinations of the components for the three modes. The model is written as (Kroonenberg, 1983, chap.2; Kroonenberg and Basford, 1989):
![]() | [2] |
The component loadings of the female lines can be investigated jointly with the component loadings of the male lines, by projecting them together in one space, as it then becomes possible to display the interaction between female and male inbred lines. The plot of the common space is called the joint biplot, a variant of Gabriel's (1971) biplot (Kroonenberg, 1983; Basford et al., 1990) and is constructed from the core matrix as follows. For each component r of the environment mode, the female components and the male components are scaled by dividing the core slice associated with that component r between them (using singular value decomposition) and weighting the scaled female and male components by the relative number of elements in the modes to make the distances comparable. For the rationale behind this construction and more detailed discussion, see Kroonenberg and De Leeuw (1977).
The treatment of the raw data, i.e., centering and normalization of the F x M x E input data arrays, is a central decision to obtain the "best" analysis for a particular dataset (Kroonenberg, 1983, chap. 6). In three-way data, different types of centering and normalization lead to different solutions. There are also many more ways of centering and normalizing to choose from, compared with two-way data. Kroonenberg (1983, chap. 6) and Harshman and Lundy (1984, p. 225253) give detailed discussions on this issue.
Centering and Normalization to Retain GCA x E + SCA x E Interactions in the Three-Way Table (3-Mode PCA1)
The two-way 16 x 11 (hybrids x environments) array of BLUPs for oil yield was double-centered by subtracting the across-hybrids environment means and the across-environment hybrid means, and adding the overall mean (residual from additivity, Gabriel, 1978), as in the AMMI model (Gauch, 1988) for two-way tables. The purpose of centering was to eliminate from the analysis those means that should not be modeled multiplicatively. To understand its implications, three-way data can be written as if they were generated by a three-way analysis of variance as that detailed in Eq. [1]. The centering applied leads to
![]() | [3] |
Normalization, which is also named scaling or standardization, is the process of equalizing sums of squares. This process is necessary since environments showed different oil yield means. Normalization within-environments (i.e., rijk = yijk''/S..k) was used because it causes each environment to have a mean of zero and a standard deviation of one, this being the most appropriate treatment for reducing the influence of environmental main effects (Cooper and DeLacy, 1994). The standardized residuals from additivity were further arranged in a three-way 4 x 4 x 11 (F x M x E) matrix and three-mode PCA was applied by the program TUCKALS3 (Kroonenberg, 1994).
Centering and Normalization to Retain SCA + SCA x E Interaction in the Three-Way Table (3-Mode PCA2)
The three-way (F x M x E) table of oil yield BLUPs was arranged with the females as rows, the males as columns, and the environments as slices. The three-way dataset was centered within environment by subtracting both the across-females male means and the across-males female means and adding the overall environment mean (Kroonenberg, 1983, chap. 6), which leads to
![]() | [4] |
With the aim of eliminating the environmental effects, normalization across slices was used (Kroonenberg, 1983, chap. 6), which in this type of analysis implies the division of the double-centered data by the within-environment standard deviation over all females and males (i.e., rijk = yijk''/S..k; Kroonenberg and Basford, 1989; Basford et al., 1990). Three-mode PCA was applied to the 4 x 4 x 11 (F x M x E) double-centered and normalized BLUPs residual array for oil yield using the program TUCKALS3 (Kroonenberg, 1994).
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
|
In Fig. 1b, which retains the GCAM and the SCA, the ideal female tester is apparently A3, which is parallel to the PC1 and shows an almost zero ordinate on PC1 (Yan and Hunt, 2002). The PC1 is related to the GCAM across environments, with the highest yielding males at the right hand side of the diagram. The PC2 is dominated by the strong SCA between A2 and R3.
AMMI2 Model
The AMMI2 biplot of the first two components for oil yield explains 99% of the SCA between males and females across environments (Fig. 1c) and can be interpreted as follows (Narro et al., 2003). A positive perpendicular projection of a male point on an individual female vector indicates that both lines combine well (positive SCA). The greater the distance from the origin to the intersection of a male projection on a female vector, the better this male combines with this female line. A male projection around the origin of a female vector indicates nearly null SCA, while a male negative projection on an individual female vector indicates that these lines do not combine well in relative terms (negative SCA).
The patterns of SCA across environments revealed by Fig. 1c indicate for female A2: a strong positive SCA with the male R3, a positive SCA with R2, a negative SCA with R4, and a strong negative SCA with R1. Since the main effects of the lines (GCA) have been removed, this pattern does not directly reflect yield rankings. In this particular case, the hybrid A2/R4 (negative SCA) still had a greater average oil yield across environments than did hybrid A2/R2 (positive SCA) (Table 1). However, this analysis allows the detection of SCA across environments, which is highly useful information for a breeder, especially when selecting parents for new breeding populations. For example, male lines of similar SCA pattern for oil yield with an individual elite female (e.g., R1 and R4 show positive SCA with A3, Fig. 1c), but that are genetically distant, could be combined in experimental crosses to deliver inbred lines that bring better opportunities for yield improvement in hybrid combinations with this or similar female. In the genotypeenvironment system under study, the females A1 and A2 showed a positive SCA with the males R3 and R2, while the females A3 and A4 showed a positive SCA with the males R1 and R4 across environments (Fig. 1c).
GCA x E and SCA x E Interactions (Three-Mode PCA1)
The three-mode model with 2 x 2 x 3 components for F x M x E, respectively, was considered adequate for fitting the data (r2 = 0.51, Table 4), on the basis of informal judgments of the small increases in r2 when more dimensions were used and on the increased difficulty of interpretation (Kroonenberg, 1983). In this model, the two components of the female mode accounted for 40 and 11% of the variation; the two components of the male mode accounted for 31 and 20%; and the three components for the environments accounted for 35, 10, and 5% of the variation (Table 4). Not all females, males, and environments were fitted equally well by the model. For example, the model only explained 13% of the interactions of the female line A4 compared to the overall fit of 51% (Table 4). While environments TA and VM were less well fitted by the model, the selected model accounted for more than 50% of the contribution of BA, GP, QQ, SB, VO, and VT to the total GCA x E and SCA x E interactions (Table 4).
|
Northern and central regions of Argentina are different subregions for sunflower, i.e., their differences in terms of genotypic discrimination for oil yield are large and repeatable over years (de la Vega et al., 2001). The northern region shows a higher variation across years than the central region in terms of the environmental factors affecting sunflower growth and development (Chapman and de la Vega, 2002). Consistent with this observation, the ratio of the genotype x year and genotype x year x location interaction effects to the genotypic effect within-region is higher in the northern region (de la Vega and Chapman, 2006). In the central region, 1 yr of trials can adequately represent the target environment if sufficient locations are sampled (de la Vega and Chapman, 2006). Conversely, in the northern region different years can show opposite discrimination effects on northern-adapted hybrids differing in maturity. According to the results of the sunflower private trials of Advanta Semillas SAIC (unpublished data), the 20042005 season could be considered representative of the northern environments that promote a relative improvement of early-maturity hybrids.
The joint biplots of female and male inbred lines for the first and the second environment components will be used to investigate the relationships between females and males (GCA and SCA) associated with the two specific patterns of genotype discrimination, namely "Central versus North and South" and "North versus Central and South." The third environmental component reflects unpredictable G x E interactions (Table 4). The proportion of variability explained by this component (5%) is too small to justify discussion.
First Environment Component: Interpreting the Contrasting Effects of Central versus Northern and Southern Regions on GCA and SCA
The component weights for the first and second axes of the joint biplot of female and male inbred lines in the first environmental component were 0.29 and 0.07, respectively (Fig. 2a
). This bidimensional joint biplot displays those aspects of the relationships between female and male lines that are influenced by the differences between northern and southern versus central environments, after the genotype and environment effects have been removed. These effects can be selected in breeding for specific adaptation to the central subregion or discarded in breeding for broad adaptation to the undivided target region.
|
In this study (Fig. 2a), the male inbred lines showed a wider pattern of interactions than the female lines, since they occupy a wider range of the Euclidean space. While all females can interact positively (negatively) with some males (e.g., all females have a positive G x E interaction effect with R2 in VO and a negative effect with R4 in the same environment), all males together cannot interact in the same direction with any female.
A practical application of the above observation is related to the definition of selection strategies for broad or specific adaptation. If division of the breeding program's target population of environments in subregions is appropriate and sources of specific adaptation to those subregions have to be identified, they are more likely to be found within the male germplasm base of the program than within the female one. This is assuming that such germplasm bases are adequately represented by the sampled lines. Regarding the selection of testers, for example, a male interacting positively with all females in the southern region (e.g., R2) could be used as a tester for this subregion, while a male interacting positively with all females in the central region (e.g., R4) could be used as a tester for this contrasting subregion. The male R3 showed positive and negative G x E interaction effects in the southern region with A3 and A1, respectively, while R1 is located around the origin of the joint biplot, showing a weak pattern of responses in this environmental component.
Second Environment Component: Interpreting the Contrasting Effects of Northern versus Central and Southern Regions on GCA and SCA
The component weights for the first and the second axes of the joint biplot of female and male inbred lines in the second environment component were 0.10 and 0.01, respectively. Therefore, major effects of this particular environmental contrast can be described in a single dimension, corresponding to the first component of the joint biplot (Fig. 2b). When one dimension effectively explains most of variability, the joint biplots collapse into a single line, in which it is possible to include the component loadings of the environment mode as well. In such a case, a product term to compare scores may be calculated as a product of any combination of the scores of the three modes, e.g., Basford et al. (1990) and Chapman et al. (1997), for genotypes, environment, and attributes. From Fig. 2b, female A2 (score 1.37), for example, will show a negative G x E interaction effect with male R1 (score 0.15) in CM (score 0.44) (i.e., 1.37 x 0.15 x 0.44 = 0.09) and a positive G x E interaction effect in VO (score 0.74).
All females have a negative (or near zero) score in this joint biplot, while the males clearly split into two groups, showing, as in the case of the first environment component, a higher variability in terms of their G x E interaction effects (Fig. 2b). This confirms the notion that a search for specific adaptation traits in this particular breeding program should focus on the male germplasm base. Male R1 (negative score) showed a positive G x E interaction effect with all females (negative score) in the northern environments (positive score). Males R2, R3, and R4 (positive score), on the other hand, showed a positive G x E interaction effect with the same females in the southern and central environments (negative score).
In summary, a joint analysis of the two joint biplots (Fig. 2), which together account for the pattern-rich G x E interaction variability of the system under study, demonstrates that male R4 showed a positive G x E interaction effect with all females in the central environments; males R2 and R3 showed a positive G x E interaction effect with most females in the southern environments (except A1 for R3 in environment component 1); and R1 showed a positive G x E interaction effect with all females in the northern environments associated with environment component 2 and a weak response in environment component 1. This information is highly relevant for the definition of the use of tester lines in the breeding program. Regarding the female lines, they showed similar patterns of combining ability by environment interactions. However, line A4 was poorly explained by the model when comparing with the other females, suggesting that it does not constitute an appropriate representation of the interaction patterns of the female germplasm base of the program and should not be used to test the combining ability of new male inbred lines.
GCA and GCA x E Interaction Using 2-Mode PCA
The environment-centered and normalized biplot of the first two components for female oil yield explains 88% of the GCAF + GCAF x E interaction (Fig. 3a
). The environment-centered and normalized biplot of the first two components for male oil yield explains 86% of the GCAM + GCAM x E interaction (Fig. 3b) and can be interpreted in analogous way to the biplot of Fig. 2a.
|
However, in the two biplots (Fig. 3), the PC2 adds the environmental dimension to the analysis, reflecting the main pattern of GCA x E interaction for females and males and discriminating between the lines that showed high GCA in both the SREG2 model (Fig. 1a and 1b) and in the direction of the PC1 of Fig. 3a and 3b. In the GCAF + GCAF x E interaction biplot (Fig. 3a), the three environments of the northern region (i.e., VO, MA and SB) had negative scores for the PC2, while most central environments had positive scores for that PC. The female A3 had, on average, a higher GCA than A2 in the northern environments and a lower GCA in the central environments. Thus, a program aimed at exploiting specific adaptations to the subregions should select A3 and A2 as female testers for the northern and central subregions, respectively. Similar analysis could be done for the GCAM + GCSM x E interaction biplot (Fig. 3b), where, on average, the male R1 had a higher GCA than R4 in the northern environments and a lower GCA in the central environments.
SCA and SCA x E Interactions (Three-Mode PCA2)
The three-mode model with 2 x 2 x 2 components for F x M x E, respectively, was considered adequate (r2 = 0.76, Table 4). The first environmental component displays the overall variability across environments, with positive scores for all environments. The second environment component reflects the portion of the G x E interaction that separates CM, SB, and VM, with negative scores, from GP, MA, TA, and VO, with positive scores. This environmental contrast is considered as unpredictable as there is no obvious regional association. The proportion of variability explained by this component (5%) is too small to justify discussion.
The joint biplot of female and male inbred lines for the first environment component can be described in one dimension, which retains 63% of the SCA + SCA x E interaction variability (Fig. 4 ). The interpretation rules of this one-dimensional joint biplot are the same as those for the second environment component of three-mode PCA1 analysis (Fig. 2b).
|
Comparison of the Models to Display GCA, SCA, and Their Interactions with Environments
Table 5 summarizes the main features of all methods utilized in this study to partially explain GCA, SCA, and their interactions in a half-sib mating (NCII) multienvironment trial. Two combinations of methods were explored to account for all sources of variation. In Strategy 1, two-mode PCA was used to display GCA + SCA across environments and three-mode PCA to display GCA x E + SCA x E interactions. In Strategy 2, two-mode PCA was used to analyze GCA + GCA x E interaction and three-mode PCA to display SCA + SCA x E interaction. In this example, Strategy 2 explained almost 18% more of the total G + G x E interaction variation of the system than did Strategy 1 (Table 3), proving to be an adequate analytical approach to understand the combining ability of the core germplasm of a hybrid crop breeding program within its target population of environments. However, in another dataset, Strategy 1 may be superior to Strategy 2, so both strategies need to be considered when dealing with a F x M x E dataset.
|
| CONCLUSIONS |
|---|
|
|
|---|
The SREG2 and AMMI2 models are able to identify the lines showing the higher GCA and the main patterns of SCA across environments. However, the compression of the environmental dimension implicit in these analytical approaches do not permit display of the variation of those patterns across the population of environments targeted by a breeding program. Two-mode PCA of environment-centered and normalized F (M) x E tables and three-mode (F x M x E) PCA provide effective tools to complement these methods in visualizing and studying GCA x E and SCA x E interactions, allowing the selection of the best tester for each selection strategy (broad or specific adaptation) and displaying the variability of the tested lines for adaptation and combining ability.
| ACKNOWLEDGMENTS |
|---|
Received for publication August 31, 2005.
| REFERENCES |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME< |