Genotype × Environment Interaction Studies of Promising Teppi Coffee (Coffea arabica L.) Genotypes in Southwestern Ethiopia

Coffee is the stimulant crop plant that belongs to the genus Coffea, in the family Rubiaceae. (e yield variation and pattern of the crop varied within short geographic distance and thus attributed to low productivity and unpredictable production. (e objective of this study was to determine the extent of genotype× environment interaction on yields of promising Teppi coffee genotype. In this study, seventeen Arabica coffee genotypes, representing coffee growing areas of Teppi, southwestern Ethiopia were evaluated at six environments. (is study was conducted by using completely randomized block design with two replications. (e extent of genotype× environment interaction was analyzed by using additive main effect andmultiplicative interactionmodel and genotype main effect and genotype× environment interaction biplot analysis model. (e analysis of variance for yield revealed the presence of highly significant differences (P< 0.01) within genotypes, environments and genotype by environment interaction. From the total sum of squares of genotype× environment interaction of yield, environment accounted for 51% while genotypic and GEI effects were 20.5% and 25.32%, respectively. Generally, genotypes with superior performance and suitable environments for coffee were identified. However, further researches have to target hybridization of those superior genotypes with elite lines in feature Arabica coffee breeding program.


Introduction
Coffee is a useful stimulant beverage plant that belongs to the genus Coffea, in the family Rubiaceae. e genus Coffea contains around 124 species [1]. e crop is commonly grown in tropical and subtropical agroecological zones [2]. Arabica coffee (Coffea arabica) originated in Ethiopia where it grows wild in the natural forests of southwestern parts of the country, and there is a great genetic diversity in the region [3].
Arabica coffee is an essential commercial crop in the national economy of Ethiopia. It is still Ethiopia's number one export item [4]. It accounts for 25 to 30 percent of total export earnings. However, at present the crops contribution to total export earnings has increasingly declined due to increased exports of other agricultural and industrial commodities like gold, flowers, khat (Catha edulis), textiles, and leather products [5].
Eastern Africa have suitable lands for Arabica coffee production. Ethiopia has potential to produce huge amounts of differentiated high-quality green coffee, which is liked for its unique flavor and taste. According to Central Statistical Agency [6], the estimated area of land covered by coffee in Ethiopia is about 700474.69 ha. On the other hand, many growing locations in the country did not benefit from this huge potential, as they should have. In the same way, despite the consequences of coffee plays central role in the national economy and in spite of the fact that the country is home of C. arabica, the country's coffee production is characterized by low productivity about 619 kg ha −1 [7], which is below global productivity of 790.7 kg ha −1 [8]. is could be because of a number of factors, among which lack of recommended genotypes to specific and wider ranges of production areas is one.
In cultivar development plan, multi-location performance trial across diversified environments is undertaken to understand the nature of genotype × environment interaction (GEI) in that specific set of genotypes and environments under study. By doing so, genotypes that have affinity to doing well in a wider range of environment and specific adapting genotypes would be identified so that they can be recommended based on situations. [9]. From earlier study in Ethiopia, indigenous coffee cultivar has its own specific adaptation region. So that the coffee growing environment is subdivided into different coffee growing regions to reduce the effect of GEI. e extent of GEI helps to know about the probable area of adaptation of specific genotypes. It is also useful in determining efficient methods of using time and resources in a breeding program. Previous studies of GEI on C. arabica have illustrated significant interaction of genotypes with environment for yield and yield-related traits [10][11][12]. In addition, Montagnon et al. [13] reported the presence of strong GEI in C. canephora with some stable genotypes for bean yield. Mesfin and Bayetta [14] also reported that C. arabica varieties that show superior performance at one site of coffee growing region did not show good performance at other locations of a dissimilar geographic region. Montagnon et al. [13] also identified that the existence of significant interaction affects genotypes performance across different coffee growing agroecologies.
Although some studies were done in the Ethiopia, there is no information about the GEI of promising Teppi coffee genotypes in diverse agroecologies of Southwestern Ethiopia. So, this study was done with the objective of determining the extent of GEI on yield of some promising Teppi coffee genotypes in Southwestern Ethiopia.

Description of the Study Area.
e study was done at three locations, namely, Teppi, Gemadro, and Godare, for two cropping seasons (Table 1).

Genetic Materials Used and Experimental Design.
Seventeen coffee genotypes were used in this study (Table 2). e genotypes consisted fourteen selections screened from previous trials, which originally were collected from coffee growing areas of Teppi and its surroundings were planted along with three checks (Geisha, Catimor J-19, and Catimor J-21). e genotypes were planted in June 2013. e experiments were conducted by superimposing on first and second bearing trees, which were grown under the shade trees of Sesbania sesban and Albizia schimperiana simultaneously at recommended space. It was arranged in randomized complete bock design with two replications. A plot consisted of a double row with ten trees per row and had border on both sides of the block to reduce border effects. Spacing between rows and plants was 2 m × 2 m and spacing between blocks was maintained at 4 m. All necessary field supervision was done as per the recommendation [15].

Yield of Clean Coffee in Kg Ha
e total fresh cherry yield was harvested from all the trees in a plot and measured in gm. at was used to compute mean yield per trees. e clean coffee yield in kg ha −1 was obtained by multiplying the yield of the fresh cherry by percent out-turn [16].

Statistical Analysis
e joint analysis of variance over environments was analyzed by the PROC GLM method in SAS [17] version 9.3 software. Prior to conduct combined analysis, homogeneity of variances between environment was checked by Bartlett's test. e effects of the genotypes, locations, and years as well as their interactions were determined by using SAS version 9.3. e following statistical model was used for combined analysis of variance over locations: where μ � general mean, τ i � effects of genotype i, ε j � effects of environment j, c ij � effects of genotype × environment interaction, ρ k(j) � the k th block effect within location j; and ε ijk � residual variation or error assumed to be normally distributed with mean 0 and variance σ 2 (i.e., ε ijk ∼ N (0, σ 2 )).

AMMI Model.
Additive main effect and multiplicative interaction (AMMI) model, which combines standard analysis of variance with PCA [18], was used to investigate GEI. In AMMI model, the contribution of each genotype and each environment to the GEI was assessed by use of the biplot graph display in which the means of genotypes are plotted against the scores of the IPCA1 [18]. It was done by using GEA-R [19] version 4.0. e AMMI model is as follows: where Y ge � the yield for genotype (g) in environment (e); µ � the grand mean; α g denotes genotype deviation; β e � environment deviation; λ n � the singular value for component n; c gn is the eigenvector value for g; and δ en � the eigenvector value for e and the residual term is ρ ge .

GGE Biplot Analysis.
e GGE biplot method is composed of two concepts, the biplot concept and the GGE concept, which were used to visually analyze the genotypes at different environments. e GGE biplot analysis was done by using Genotype x Environment Analysis with-R [19] version 4.0. e general model for GGE biplot is as follows: where Y ij � the performance of the i th genotype in the j th environment; μ � the grand mean; β j � the main effect of the environment j: λl and λ2 � singular value for IPCA1 and IPCA2, respectively: ε i 1 and ε i 2 � eigenvectors of genotype i IPCA1 and IPCA2, respectively: η j 1 and η j 2 � eigenvectors of environment j for IPCA1 and IPCA2, respectively; and ε ij � residual associated with genotype i and environment j.  e mean performance of genotypes was different from one environment to another.

Result and Discussion
is indicates that different genotypes react differently to a particular environment. Similarly, different researchers [10,11,21] reported that different coffee genotypes responded differently to diverse agroecology of southwestern Ethiopia.

Combined Analysis of Variance over Locations.
e analysis of variance for coffee bean yield revealed the presence of highly significant difference (P < 0.01) between genotypes, environments, and their interaction (Table 3). e existence of considerable GEI indicated that the genotypes showed inconsistent performances across the tested environments. is could be due to unpredictable variations of weather and soil characteristics. e presence of considerable GEI showed the differential performance of coffee genotypes across environments.
is means different genotypes react differently to a particular environment. Related results were reported by Yonas et al., Afework, and Lemi et al. [10][11][12]. e total variations explained by environment genotype and GEI were 51%, 20.5%, and 25.32% (Table 3). e high proportion of the environment sum of squares means the  International Journal of Agronomy major factor that manipulates yield performance of coffee genotypes is the environment. e GEI is highly significant (P < 0.01) scoring 25.32% of the total sum of squares implying the call for investigating the nature of degree of difference reaction of the genotypes to environments. e presence of GEI shows that the performance of one genotype might be superior to the other genotypes in one environment but inferior in another environment. It complicates interpretation of the results. As a result, it is difficult to identify consistently superior genotypes across environments. Similarly, Afework [11] reported that the highest variation explained was due to environment and GEI. In general, from the combined ANOVA (Table 3), superior genotypes across environments cannot be identified by considering their mean yield performance because GEI is highly significant. Because of the interactions between genotypes and environments, yield of genotypes tested across environments fluctuated. erefore, it was a problem to identify varieties consistently giving high yields in locations with diverse environmental conditions. e existence of highly significant GEI decreases the usefulness of genotypes. e suitable statistical analyses are required for quantifying GEI. Furthermore, the conventional analysis of variance determines the values of each variation source and the significance of the contribution of each component, but it does not partition the interaction into a number of components and hence other types of analyses have to perform. For this reason, such multi-environment trial data along with considerable GEI require measures of stability analysis.
at helps to get more information on the GEI as well as to evaluate the adaptation regions of the genotypes according to their favorable interaction. Similarly, previous researchers [12,22] reported the significant effect of GEI on C. arabica in Ethiopia.

AMMI Analysis.
e AMMI analysis of variance showed highly significant differences (P < 0.01) for environments, genotypes, and their interaction (Table 4). e F-test was highly significant (P < 0.01) for the first five interaction principal component analyses (IPCAs). e IPCAs were put in order according to their decreasing significance. e AMMI analysis of variance for bean yield showed that most of the total sum of squares was explained by environment (52.65%) followed by genotypic (21.19%) and then by GEI (26.16%) illustrated in Table 4. is indicated that the environments used were diverse, with large differences that caused most of the variation in yield, confirming the significant influence of environments on the yield performance of coffee. Similarly, Afework, Lemi et al., and Demissie et al. [11,12,22] reported the related results. (Table 5) e Gollob test discovered that the first five IPCAs were significant (P < 0.01), indicating that the total information contained in GEI can be explained using these IPCAs. Out of the total IPCAs, the first five IPCAs axes explained 100% of the GEI sum of squares. e first IPCA1 captured about 44.11% of the interaction sum of squares, while the second IPCA explained 27.03% of the GEI sum of squares. is is related to the findings of Yonas et al., Lemi et al., and Demissie et al. [10][11][12]22] who reported significance of the 4 to 6 IPCAs in C. arabica for bean yield evaluation. e AMMI model performed by using the first two IPCAs and the rest of the IPCAs mostly capturing noise showed that IPCA 1 against IPCA2 is generally informative. e pattern in GEI of the given coffee data set was also predicted by using the first two PCAs of genotypes and environments, since IPCA1 and IPCA2 cumulatively accounted for 71.14%, which was greater than half of a total GEI.

AMMI 2 Biplot Analysis.
e AMMI 2 biplot is performed by using the genotype and environment score of the first two successive AMMI components [23]. e first IPCA captured 44.11%, while the second IPCA captured 27.03%. e first two successive IPCAs cumulatively captured 71.14% of sum of square of the GEI of studied genotypes. From earlier yield experiment of GEI in C. arabica, [22] identified that the first two IPCAs explained 74% of AMMI total interaction sum of squares, [11] reported 58.55%, and [12] reported 63.3%. e genotypes and environments located far away from the center are more reactive or unstable, while genotypes that are positioned closer to the biplot center have superior stability achievement [24]. Hence, genotypes 32/82, 45/82, 235/71A, 17/79, 20/79, and Catimor J-19 were plotted relatively close to the origin in the AMMI 2 biplot indicating their related yielding potential to all environments. erefore, genotypes 32/82, 17/79, and Catimor J-19 were considered as a high yielding and widely adopted genotypes indicating their minimum contribution to the total GEI variance. On the other hand, genotypes like Geisha, Catimor were relatively distant from the origin and have considerable contribution to the GEI variance considered as specifically adopted to their individual favorable environments ( Figure 1). e AMMI 2 biplot showed that E1 (Teppi 1), E4 (Teppi 2), E5 (Gemadro 2), and E6 (Godare 2) were far from the origin indicating that these environments contributed higher amount of variation to the total GEI. However, because of their longest distance between their marker and the origin, genotypes inconsistency at this environment may not exactly reflect the average genotypes performance across environments. On the other hand, E3 (Godare 1) and E2 (Gemadro 1) are located close to the origin indicating their lower contribution to the GEI variance. It shows that they are the least discriminating and stable environments ( Figure 1).

GGE Biplot Analysis
3.5.1. e "Which-Won-Where" Patterns. e partitioning of GEI through GGE biplot analysis showed that IPCA1 and IPCA2 accounted for 52.51% and 17.59% of sum of squares, respectively, with 70.1% variation for yield ( Figure 2). e six environments fell into two sectors with different winner genotypes and the biplot showed five vertex genotypes, Geisha, Catimor J-19, 42/82, 39/82, and 37/82. ere are rays, which divided the biplot into five sections. e genotypes fell into five sections but all the tested environments fell into two sections. e genotypes that fell on vertex of all sectors were the ones that performed the highest yield in those environments, which fell within that sector. In this study, GGE-biplot analysis identified two different coffee growing mega-environments. e first environment contains E1(Teppi 1), E2(Gemadro1), E4 (Teppi 2), E3 (Godare 1), and E6 (Godare 2) with a vertex genotypes Geisha and Catimor J-19; the second environment contains the environment E5 (Gemadro 2) with 32/82 the winner genotype. It had also been observed that no environments fell into sectors where genotypes 37/82, 42/82, and 39/82 were the vertex genotypes, indicating that these genotypes were not adapted in any of the test environments. Some authors also reported related findings in the country; Lemi [21] identified four different C. arabica growing mega-environments and Afework [11] identified four different C. arabica growing mega-environments.

Mean Performance and Yield Stability of Genotypes.
Yield performance and stability of the genotypes were evaluated by an average environment coordination method.
In this method, genotypes were ranked beside the average environment coordinate axis through an arrow demonstrating the uppermost value based on their mean performance across every part of environments. e average environment coordinate ordinate separates genotypes, which are above-average means and below-average means. Accordingly, genotypes with above-average means were 45/82, 3/82, 32/82, Catimor J-21, 17/ 79, 45/82, 48/82, Catimor J-19, and Geisha, while genotypes with below-average means were 42/82, 20/79, 37/82, 29/82, 44/ 82, 235/71A, 22/79, 28/82, and 39/82 (Figure 3). e line, which passes through the center and that is perpendicular to the average environment coordinate by having double arrows, represents the consistency of genotypes. e direction of genotypes away from the biplot, basis, on the axis, shows greater GEI and instability. e better genotypes are those that perform high mean yield and similar consistency across the environment. In the biplot, they are close to the origin and have the shorter vector from the average environment coordinate. e longer projection to the average environment coordinate, without considering direction, represents a greater tendency of the GEI of a genotype that means less constancy across environments. Accordingly, genotypes 17/79, 22/79, 29/82, Catimor J-19, 45/82, 48/82, and 44/82 were the most stable. On the other hand, the genotypes Geisha, 3/82, 42/82, 39/82, and 37/82 were far from AEC (long vector) indicating their least constancy (Figure 3). ese results are in accordance with the report of Lemi and Afework [11,21] on Limmu and Ilu Abbabora coffee genotypes, respectively. ey found the ideal genotypes by using this mean performance and yield stability measuring parameters.

Rank Genotype and Environments
Relatively. An ideal genotype is one that shows superior performance in mean yield and that has high performance in constancy. e interior of the concentric circles ( Figure 4) represents the situation of an ideal genotype, which is defined by a projection onto the mean-environment axis. For the reason that the units of IPCA1 and IPCA2 for the genotypes are the novel component of yield in the genotype-focused scaling, the units of the AEC abscissa and ordinate should also be the original unit of yield. Accordingly, genotypes Catimor J-19, 17/79, 32/82, and 3/82 followed by Geisha, which fell closer to the center of concentric circles, were desirable genotypes in terms of higher yield potential and constancy, compared with the rest of the genotypes. Genotypes 39/82, 37/82, 42/ 82, 44/82, 45/82, 235/71A, 48/82, and 29/82, followed by 20/ 79, were unattractive genotypes located distant from the first concentric circle of the ideal genotype. e result of the study was in accordance with Lemi and Afework [24,11) that they had found an ideal genotype in C. arabica growing domain of Ethiopia.
In the ways of ranking environments relative to the ideal environment method, the average IPCA1 and IPCA2 scores of all environments, represented by a small circle (Figure 4), describe an average environment. Like that of ideal genotype, an ideal environment is more advantageous if it is located closer to the ideal environment. e model environment, represented by a symbol of a small circle with an arrow pointing to it, is mainly discriminating the genotypes and representativeness of the other test environments. Figure 4 indicates that E4 (Teppi2), which fell in the interior part of circles, was an ideal test environment in terms of being the most representative of the all environments and being mainly powerful to discriminate genotypes, where Gemadro2 was far from the ideal environment and considered as less powerful to discriminate genotypes.
Discriminating ability and representativeness are the important properties of an experimental location. Ideal locations have to be highly differentiating for the evaluated genotypes and simultaneously representative of the target location. Like ideal genotype, an ideal location is defined and shown by the small circle with an arrow pointing to it. It shows that experimental site is more desirable and discriminating when situated closer to the origin of concentric circle or to an ideal environment. is result is related to the findings of Lemi and Afework [11,21] that they identified an ideal genotypes and environments as well as undesirable genotypes by using ranking genotypes and environment method.

Conclusion
e oscillating agroecology within a short distance in Ethiopia leads to fluctuation of coffee yield due to GEI and  thus contributes to low productivity. erefore, testing genotypes across different environments is important to identify genotypes with low GEI, which can increase C. arabica productivity in the country. e large sum of squares and highly significant performance mean square of environments showed that the environments were diverse, by means of large differences among environmental means causing difference in yield. is shows the main influence of environments on yield performance of the crop. e existence of important genotypes, environments, and GEI indicates that the considerable influence of environments on C. arabica and the genotypes showed instability across environments, and high percentage variation of yield was due to environmental effect.
e Gollob test revealed that the first five IPCAs were significant; and the total information contained in GEI explained IPCAs. e patterns in GEI of the given coffee data set were predicted by the first two IPCAs of genotypes and environments, since IPCA1 and IPCA2 cumulatively accounted for 71.14% of the total GEI. Among the tested genotype, Catimor J-19 showed lower GEI across environments. GGE biplot identified two different mega-environments and two winner genotypes in each megaenvironment. e first mega-environments contained E4 (Teppi 2), E6 (Godare 2), E3 (Godare 1), and E1 (Teppi 1) with a vertex genotype Geisha and Catimor J-19; the second environment contained E2 (Gamadro 1) and E5 (Gemadro 2). It has also been identified that genotypes 37/ 82 and 39/82 were not adopted in any of the test environments.
Even though no genotypes showed superior performance across the experimental site, some genotypes with better mean performance coupled with lower GEI were identified. erefore, priority should be given to the genotypes that showed desirable performance. In this case, Catimor J-19 is the most desirable genotype for its low GEI and high yield potential. Generally, in the study, genotypes were adapted and an ideal location was identified. However, further researches have to target hybridization of those superior genotypes with elite lines in feature Arabica coffee breeding program.
Data Availability e data used during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.