Twenty cassava genotypes were arranged in a randomised complete block design with three replications and assessed for growth and yield stability using the additive main effect and multiplicative interaction (AMMI) analysis. Highly significant (P<0.001) effects of genotype, environment, and genotype ⁎ environment interaction were observed for all traits studied. The AMMI analysis of variance indicated that genotype accounted for 51% of the total sum of squares for height at first branching followed by environment (33%) and interaction (15%). For fresh root yield, environment effects accounted for 37% of the total sum of squares, whilst genotype and interaction accounted for 32% and 29%, respectively. Genotypic variances for harvest index (HI), plant height, storage root yield, and dry matter content contributed a greater proportion of the phenotypic variance indicating stronger genetic control. This suggests better chance of progress in the genetic improvement of these traits. Genotype MM96/1751 combined high yield with stability according to the yield stability index ranking across environments. On the other hand genotypes UCC 2001/449 and 96/1708 though high yielding were unstable according to AMMI stability value scores. However they can be tested further in more environments to ascertain their specific adaptability for release to farmers for cultivation to boost cassava production and ensure food security.
Generation Challenge ProgrammeAlliance for Green Revolution in Africa1. Introduction
Cassava (Manihot esculenta Crantz) is grown in a wide range of environments but severe yield losses occur in traditional cultivars grown by poor farmers [1–3]. Root yield losses in poor drought-prone environments can be as high as 80% when compared with root yield in optimum environments [4]. This is because genotypes exhibit different levels of phenotypic expression under different environmental conditions resulting in crossover performances [5]. Crossover performances by genotypes in different environments result from differential responses under varying environmental conditions [6, 7]. This was defined as genotype∗environment (G∗E) interaction [8]. Studies on cassava’s performance in different environments have indicated significant G∗E effect for root yield and yield components [9–12].
Genotype∗environment interaction also results from differences in the sensitivities of genotypes to the conditions in the target environment [13] leading to inconsistent performances of genotypes across environments. This limits the efficiency of selection of superior genotypes [14]. In G∗E interaction, the magnitude of the observed genetic variances changes from one environment to the next and tends to be larger in better environments than poorer environments [15]. Breeding cassava for yield involves different attributes of genotypes that are subject to influences by variation in environmental conditions [16]. The objective of most cassava yield improvement programmes is to identify and select high yielding and stable genotypes across several environments and seasons [11]. The efficiency and success of such selection depend on the consistency of performances of genotypes in varying environments [17]. For this reason, genotypes are tested in diverse environments to assess their adaptability and stability [5]. Genotypes whose G∗E effects are not significant are said to be stable [14].
Several methods have been used to assess the G∗E effect and stability in crop performances. Eberhart and Russell [18] proposed joint regression analysis to estimate the average performance of a genotype in different environments relative to the mean performance of all genotypes in the same environment. Multiplicative models which include the additive main effect and multiplicative interaction (AMMI) model have also been used [19, 20]. The AMMI model fits the sum of several multiplicative terms rather than only one multiplicative term in assessing the performance of genotypes in different environments [21]. AMMI analysis can be used to determine stability of the genotypes across locations using the PCA (principal component axis) scores and AMMI stability value (ASV) [22]. The ASV is based on the AMMI model’s IPCA1 and IPCA2 (interaction principal components axes 1 and 2, resp.) scores for each genotype [23]. Genotypes having the least ASV are considered as widely adapted genotypes. Similarly, IPCA2 score near zero indicates more stable genotypes whilst large values represent more responsive and less stable genotypes.
Stability per se does not give much information about the level of yield so Farshadfar et al. [24] and Tumuhimbise et al. [25] used yield stability index (YSI) and genotype stability index (GSI) which combined high yield performance with stability. Both the YSI and the GSI are based on the sum of the ranking due to ASV scores and yield or performance ranking. Lower YSI and GSI values indicate genotypes that combine high yield or performance with stability [24, 26, 27]. The study was conducted to assess the performance and stability of cassava genotypes for fresh root yield and other yield components. Specific objectives were to identify high yielding stable cassava genotypes and estimate the variance components as well as broad sense heritability of root yield and other traits.
2. Materials and Methods
The experiments were conducted during the 2012/2013 and 2013/2014 growing seasons. The trials were conducted at the research fields of the CSIR-Savanna Agricultural Research Institute located at Nyankpala (9°25′N, 0°58′W) in the Guinea Savannah Agroecological Zone and the Crops Research Institute located at Fumesua (6°41N 1°28W) in the Forest Agroecological Zone of Ghana. The Guinea Savannah Zone which covers over 40% of the land area of Ghana is characterised by high temperatures and low humidity for most parts of the year [28]. The rainfall pattern is monomodal and erratic with an annual mean of 1100 mm. The rains usually begin in April-May, end in October, and are followed by a long dry season (5–7 months) which normally begins in November and lasts till April. Intermittent dry spells usually occur even during the rainy season. This experimental site represented the low potential environment.
The site at the CSIR-Crops Research Institute (CSIR-CRI) was located at Fumesua within the Forest Agroecology which is characterised by a bimodal rainfall pattern, major and minor rainy seasons. The major rainy season begins in February and ends in August and then is followed by the minor season which begins in September and ends in November. This site represented the high potential environment for assessing the genotypes for high yield potential.
Twenty cassava genotypes were used for the experiment (Table 1). These included six local landraces collected from farmers’ fields, seven drought tolerant genotypes obtained from the International Institute of Tropical Agriculture (IITA), Nigeria, and seven genotypes from the International Center for Tropical Agriculture (CIAT), Columbia. The genotypes were arranged in a randomised complete block design with three replications under each growing condition. A plot consisted of four rows with ten plants in a row giving a population of 40 plants in each plot.
List of cassava genotypes used in the study.
Genotype
Source
Genotype
Source
ATR 002
Local
00/0203
IITA
ATR 007
Local
96/1708
IITA
NWA 004
Local
I91934
IITA
Pontisange∗
Local
CTSIA 45
CIAT
Biabasse∗
Local
CTSIA 48
CIAT
UCC2001/449
Local
CTSIA 65
CIAT
TME 419
IITA
CTSIA 72
CIAT
TME 435
IITA
CTSIA 112
CIAT
96/409
IITA
CTSIA 230
CIAT
MM96/1751
IITA
CTSIA 110
CIAT
Local = genotypes from local sources, IITA = genotypes obtained from the International Institute for Tropical Agriculture, and CIAT = genotypes obtained from the International Center for Tropical Agriculture. ∗Local farmer preferred varieties.
At each site the land was ploughed, harrowed, and laid out in a randomised complete block with three replications. Ridges were manually raised for the trials at CSIR-SARI, Nyankpala, but planting was done on the flat at Fumesua based on the regular practice. The 20 cassava genotypes represent the treatments. Each plot was made up of four rows with 10 plants in a row. Mature cassava cuttings measuring 25–30 cm were planted using a spacing of 1 m × 1 m on the top of ridge or flat as the case may be at an angle of 45° to the ground surface [29]. No fertilizer was applied to the trials but weeds were manually controlled using cutlasses and hoes as and when necessary.
Monthly data were recorded on growth parameters including plant height (cm), stem girth (cm), and plant height at first branching (cm) which were collected at three months after planting and continued till harvest at 12 months. Harvesting was done at 12 months after planting. Eight plants from the two inner rows were sampled for yield determination. Storage root yield was estimated per plot and extrapolated to tons per hectare. For storage root dry matter content determination, roots were first chopped into pieces (about 1 cm thick) and mixed thoroughly. Afterwards three 100 g subsamples were taken and oven-dried at 80°C for 48 hours and weighed. Dry matter content for each sample was obtained by expressing the dry weight as a percentage of the original fresh weight of the sample taken. An average of the three samples was calculated to obtain the dry matter content of the respective genotypes. Cassava mosaic disease (CMD) severity was scored using a scale of 1–5 where 1 denotes no symptom and 5 denotes highly diseased plants [30].
The data were subjected to analysis of variance using GenStat software version 12.1 [31] to determine the significance of the main effects and interactions. Combined analysis of all growth parameters as well as root yield and yield components from the different growing environments was done. Two growing environments (Nyankpala and Fumesua) and two growing seasons (2013 and 2014) were used.
Variance components (genotypic, phenotypic, and environment as well as genotype∗environment variances) were also estimated from their respective mean squares obtained from the analysis of variance table according to Ntawuruhunga and Dixon [16]. Broad sense heritability (H2b) was estimated according to Padi [32] as follows:(1)H2=σ2gσ2p,σ2p=σ2g+σ2g∗eE+σ2eER,where
σ2g is genotypic variance;
σ2p is phenotypic variance;
σ2g∗e is genotype∗environment variance;
σ2e is pooled error;
E is number of environments;
R is number of replications.
Genotypic coefficient of variation (GCV) and phenotypic coefficient of variation (PCV) were estimated according to Singh and Chaudhary [33] as follows: (2)GCV%=√σ2gX¨∗100%,PCV%=√σ2pX¨∗100%,where X¨ is grand mean.
Expected genetic advance, GA, was calculated as (3)GA=KσAH2,where GA is expected genetic advance.
K is Selection differential (2.06 at 5% selection intensity).
σA is phenotypic standard deviation.
Genetic advance as a percentage of mean (GAM) was also estimated using the following formula:(4)GAM=GAX¨∗100.
2.1. Stability Analysis
Additive main effect and multiplicative interaction (AMMI) model in GenStat 12.1 [31] was used to determine the stability of the genotypes across environments. The AMMI model first fits the additive effects for the genotypes and the growing environments (two growing environments and two seasons) and multiplicative term for genotype∗environment interactions.
The AMMI model according to Farshadfar et al. [24] is presented as(5)Yij=μ+gi+ej+∑k=1nλkαikγjk+eij,where Yij is the yield of the ith genotype in the jth environment, gi is the mean of the ith genotype minus the grand mean, λk is the square root of the eigenvalue of the PCA axis k, αik and γjk are the principal component scores for PCA axis k of the ith genotype and the jth environment, respectively, and eij is the residual. The environment and genotypic PCA scores are expressed as unit vector times the square root of λk; i.e., environment PCA score=λk0.5Yik; 0.5; genotype PCA score=λk0.5αik.
AMMI stability value (ASV) was calculated for each genotype according to the relative contributions of the principal component axis scores (IPCA1 and IPCA2) to the interaction sum of squares.
The AMMI stability value (ASV) as described by Purchase et al. [23] was calculated as follows:(6)ASV=IPCA1Sum of squaresIPCA2Sum of squaresIPCA1score2+IPCA2score2,where IPCA1Sum of squares/IPCA2Sum of squares is the weight given to the IPCA1-value by dividing the IPCA1 sum of squares (from the AMMI analysis of variance table) by the IPCA2 sum of squares. The larger the IPCA score is, either negative or positive, the more adapted a genotype is to a certain environment. Smaller ASV scores indicate a more stable genotype across environments [24].
Yield stability index was also calculated using the sum of the ranking based on yield and ranking based on the AMMI stability value. (7)YSI=RASV+RY,where
RASV is the rank of the genotypes based on the AMMI stability value;
RY is the rank of the genotypes based on yield across environments (RY).
YSI incorporates both mean yield and stability in a single criterion. Low values of both parameters show desirable genotypes with high mean yield and stability [25, 34].
3. Results
The environments differed in total rainfall amounts, mean temperature, and relative humidity (Table 2). The total rainfall recorded in 2014 was higher than that recorded in 2013 for both locations. For Fumesua in 2014, 1620.00 mm was recorded which was higher than in 2013 (1147.08 mm). Higher amount of rainfall was also recorded at Nyankpala in 2014 (1181.04 mm) than Nyankpala in 2013 (991.68 mm). Mean monthly temperature recordings at Nyankpala were higher than at Fumesua in both years. The environments at Fumesua (FUM13; 75.4% and FUM14; 81.3%) also had higher relative humidity than Nyankpala in both years (Nyank2013; 68.68% and Nyank2014; 68.5%).
Total rainfall, mean temperature, and average relative humidity during the period of the experiment.
Environment
Total rainfall (mm)
Mean temperature (°C)
Relative humidity (%)
Nyankpala 2013
991.68
28.29
68.68
Nyankpala 2014
1181.04
28.71
68.5
Fumesua 2013
1147.08
27.3
75.4
Fumesua 2014
1620
25.8
81.3
There were highly significant (P<0.001) effect of genotypes and environment for all the traits studied (Table 3). Genotype∗environment effect was significant (P<0.05) for storage root dry matter content, very significant (P<0.01) for plant height, and highly significant (P<0.001) for stem diameter, plant height at first branching, storage root yield, harvest index, and cassava mosaic disease score.
Combined analyses of variance for seven traits evaluated on 20 cassava genotypes in two ecologies of Ghana.
Source
df
Plt_ht (cm)
Stem diam (cm)
Ht_branch
Rt_yld (t/ha)
DM%
HI
CMD score
Rep
2
741.60
0.74
1739.90
178.22
11.20
0.02
0.84
Geno (G)
19
9548.50∗∗∗
0.46∗∗∗
9375.50∗∗∗
604.61∗∗∗
91.87∗∗∗
0.05∗∗∗
7.07∗∗∗
Env (E)
3
298047.90∗∗∗
2.77∗∗∗
39224.30∗∗∗
4385.50∗∗∗
105.56∗∗∗
0.16∗∗∗
17.57∗∗∗
G∗E
57
1376.30∗∗
0.18∗∗∗
993.50∗∗∗
182.23∗∗∗
20.47∗
0.03∗∗∗
1.08∗∗∗
Residual
158
750.70
0.09
481.50
61.58
13.95
0.01
0.22
Total
239
∗,∗∗,∗∗∗=significant at P<0.05, P<0.01, and P<0.001, respectively, ns = not significant (P>0.05), df = degree of freedom, Plt_ht = plant height (cm), Stem diam = stem diameter (cm), Ht_branch = height at first branching (cm), Rt_yld = storage root yield (t/ha), DM% = dry matter content (%), HI = harvest index, and CMD_Score (1–5) = cassava mosaic disease score (based on scoring chart where 1 denotes no symptom and 5 denotes highly diseased plants).
The AMMI analysis of variance indicated highly significant (P<0.001) effects of genotype, environment, and interaction for all traits (Table 4). Genotypic factors accounted for larger proportion of the treatment sum of squares for plant height at first branching, dry matter content, harvest index, and cassava mosaic disease score environment and genotype∗environment interactions. Genotype effect accounted for 15.72% of the treatment sum of squares (SS) for plant height (cm) whilst environment and interaction accounted for 77.48% and 6.80%, respectively. The first two interaction principal component axes (IPCA1 and IPCA2) accounted for 85.10% of the interaction sum of squares. For stem diameter, 32.10% of the treatment sum of squares was due to genotype effect whilst environment and interaction effect accounted for 30.38% and 37.52%, respectively. The IPCA1 accounted for 74.63% with IPCA2 accounting for 13.95%. A greater proportion of the treatment sum of squares for plant height at first branching was due to genotype effect (51.34%) followed by environment effect (33.38%) and interaction effect (15.27%). The first two interaction principal component axes accounted for a total of 79.42%.
AMMI analyses of variance for 20 cassava genotypes evaluated at 12 months after planting in two ecologies of Ghana for growth, storage root yield, and related traits.
Source
Mean squares
df
Plt_ht (cm)
St_diam (cm)
Ht_branch (cm)
Rt_yld (t/ha)
DM%
HI
CMD score (1–5)
Treatments
79
14608∗∗∗
0.35∗∗∗
4407∗∗∗
443.40∗∗∗
40.79∗∗∗
0.024∗∗∗
1.19∗∗∗
Genotypes (G)
19
9548∗∗∗
0.46∗∗∗
9408∗∗∗
604.10∗∗∗
91.76∗∗∗
0.047∗∗∗
3.15∗∗∗
Environments (E)
3
298048∗∗∗
2.77∗∗∗
38745∗∗∗
4384.70∗∗∗
105.08∗∗
0.163∗∗∗
7.07∗∗∗
Block
8
1037.00
0.31∗∗∗
1734∗∗∗
60.70ns
26.46ns
0.014∗∗∗
17.57∗∗∗
G∗E
57
1376∗∗
0.18∗∗∗
1003∗∗∗
182.40∗∗∗
20.41∗
0.010∗∗∗
0.80∗∗∗
IPCA1
21
1875∗∗
0.36∗∗∗
1308∗∗∗
324.50∗∗∗
32.74∗∗∗
0.017∗∗∗
1.08∗∗∗
IPCA2
19
1441∗
0.08ns
777∗
170.30∗∗∗
13.32ns
0.007∗
2.52∗∗∗
Error
152
736.00
0.08
436.00
63.20
13.25
0.010
0.12ns
Source
Sum of squares
df
Plt_ht (cm)
St_diam (cm)
Ht_branch
Ry_yld (t/ha)
DM (%)
HI
CMD score
Treatments
79
1154016.00
27.32
420944.00
35030.00
3222.00
1.92
285.00
Genotypes (G)
19
181421.00
8.77
348166.00
11479.00
1743.00
0.89
248.67
Environments (E)
3
894144.00
8.30
178757.00
13154.00
315.00
0.49
134.25
Block
8
8294.00
2.45
116235.00
486.00
212.00
0.11
52.72
G∗E
57
78452.00
10.25
13874.00
10397.00
1163.00
0.54
6.42
IPCA1
21
39381.00
7.65
53174.00
6815.00
688.00
0.37
61.69
IPCA2
19
27380.00
1.43
27470.00
3236.00
253.00
0.13
52.93
Error
152
111801.00
12.59
10944.00
9477.00
2001.00
0.60
2.03
% treatment SS due to G
19
15.72
32.10
51.34
32.77
54.10
46.22
53.99
% treatment SS due to E
3
77.48
30.38
33.38
37.55
9.78
25.39
21.20
% treatment SS due to G∗E
57
6.80
37.52
15.27
29.68
36.10
28.38
24.81
% G∗E SS due to IPCA1
21
50.20
74.63
51.66
65.55
59.16
67.36
85.80
% G∗E SS due to IPCA2
19
34.90
13.95
27.76
31.12
21.75
23.77
10.93
∗,∗∗,∗∗∗=significant at P<0.05, P<0.01, and P<0.001, respectively, ns = not significant (P>0.05), df = degree of freedom, Plt_ht = plant height (cm), Stem diam = stem diameter (cm), Ht_branch = height at first branching (cm), Rt_yld = storage root yield (t/ha), DM% = dry matter content (%), HI = harvest index, and CMD_Score (1–5) = cassava mosaic disease score (based on scoring chart where 1 denotes no symptom and 5 denotes highly diseased plants).
Genotype effects accounted for 53.99% of the treatment sum of squares for cassava mosaic disease severity score whilst environment and interaction effects accounted for 21.20% and 24.81%, respectively. The IPCA1 accounted for 85.80% of the interaction sum of squares with IPCA2 accounting for 10.93%. For storage root yield, environment effect (37.55%) contributed a greater proportion of the treatment sum of squares compared with genotype effect (32.77%) and genotype∗environment interaction (29.68%). The first two interaction principal component axes (IPCA1 and IPCA2) accounted for 65.55% and 31.12% of the interaction sum of squares respectively. The influence of genotype (54.10%) on root dry matter content was greater than both genotype∗environment interaction (36.10%) and environment effect (9.78%). Genotype effect (46.22%) on harvest index was greater than the effect of interaction (28.38%) and environment (25.39%). The first two principal component axes contributed a total of 91.13% of the interaction sum of squares.
A large proportion of the phenotypic variance (σ2p) for plant height was accounted for by the genotypic variance (σ2g) (Table 5). Genotypic variances (σ2g) for plant height and height at first branching were higher than the G∗E variance (σ2g∗e). Storage root yield and dry matter content were mostly influenced by their σ2g compared with σ2g∗e. Error variances (σ2e) for all traits were higher than their phenotypic variances (σ2p) except for plant height, height at branching, storage root yield, and harvest index.
Estimates of variance components, broad sense heritability, PCV, and GCV for six traits in 20 cassava genotypes from four environments.
Parameter
Traits
Plt_ht (cm)
Ht_branch (cm)
St_diam (cm)
Rt_yld (t/ha)
HI
DM (%)
σ2g
746.98
630.94
0.013
24.85
0.003
4.90
σ2g∗e
72.47
111.68
0.012
6.07
0.001
2.38
σ2e
843.90
387.10
0.091
45.08
0.004
13.27
σ2p
901.67
805.59
0.031
47.66
0.004
6.30
H2b
82.84
78.32
43.08
43.46
71.66
77.7
GCV (P)
14.75
36.19
5.14
18.36
9.79
6.53
PCV (G)
15.70
38.89
7.06
21.93
10.81
7.46
P-G
0.95
2.70
1.92
3.57
1.02
0.93
GA
51.25
45.79
0.16
6.77
0.10
4.02
GAM
27.66
65.98
6.94
27.31
17.06
11.86
σ2g=genotypic variance, σ2g∗e=genotype∗ environment variance, σ2e=error variance, σ2p=phenotypic variance, H2b=broad sense heritability, GCV = genotypic coefficient of variation (%), PCV = phenotypic coefficient of variation (%), GA = genetic advance, GAM = genetic advance as percentage of the mean, Plt_ht = plant height (cm), Ht_branch = plant height at first branching (cm), St_diam = stem diameter (cm), Rt_yld = root yield (t/ha), DM% = storage root dry matter content (%), and HI = harvest index.
The estimates of broad sense heritability varied for all the traits and were especially low for stem diameter (43.08%) and root yield (43.46%) (Table 5). Relatively high broad sense heritability estimates were observed for plant height (82.84%), height at first branching (78.32%), and storage root dry matter content (77.70%). Phenotypic coefficient of variation for all traits was higher than the corresponding genotypic coefficient of variation (GCV). Wide differences between PCV and GCV were observed for height at branching and storage root yield. PCV ranged from 7.06% to 38.89% for stem diameter and height at first branching respectively. GCV varied from 5.14% (stem diameter) to 36.19% (height at first branching). Moderate PCVs (10–20) were observed for plant height and harvest index whereas high PCV (>20) was recorded for height at first branching. Stem diameter and root dry matter content had the lowest PCVs (<10). The analysis of the expected genetic advance as percentage of the mean (GAM) indicated that height at first branching could be improved by 65.98% whilst only 6.49% progress could be made in the improvement of stem diameter. Root yield could also be improved by 27.31% whilst progress of 17.06 and 11.86% could be made in harvest index and dry matter content, respectively.
The mean performance of all genotypes across the six environments showed significant (P<0.05) genetic variability for plant height, height at branching, and cassava mosaic disease severity score (Table 6). CTSIA 65 had the highest plant height (249.67 cm) with genotype 00/0203 having the shortest plants across all environments. The lowest height at branching genotype was I91934 (36.38 cm) whilst CTSIA 65 had the highest (136.65 cm). Average stem diameter ranged from 1.84 cm (CTSIA 72) to 2.61 cm (96/409). Cassava mosaic disease scores varied significantly among the different genotypes (Table 6). MM96/1751 and 00/0203 had the lowest overall mean CMD severity score (1.00) whilst CTSIA 230 had the highest. (3.13). Average root yield for all genotypes across the four environments was 24.79 t/ha (Table 5). UCC2001/449 had the highest overall root yield of 35.97 t/ha with CTSIA 48 having the lowest root yield (11.87 t/ha). Three genotypes (UCC2001/449, MM96/1751, and 96/1708) had significantly higher root yields than the best farmer preferred variety, Pontisange, whilst four genotypes had significantly better root yield than the second best farmer preferred variety, Biabasse. Harvest index also varied from 0.43 (CTSIA 110) to 0.66 (00/0203 and MM96/1751). The average harvest index for all genotypes across environments was 0.56. Storage root dry matter ranged between 28.24% (I91934) and 38.99% (CTSIA 110) with an average of 32.38%.
Mean performance of 20 cassava genotypes evaluated in four environments.
Genotypes
Plt_ht (cm)
Ht_branch (cm)
St_diam (cm)
CMD (1–5)
Rt_yld (t/ha)
HI
DM (%)
00/0203
137.73
35.38
2.33
1.00
31.29
0.66
31.71
96/1708
142.52
75.27
2.22
1.08
32.95
0.65
30.78
96/409
180.88
58.58
2.61
1.17
27.69
0.53
29.80
ATR 002
160.53
62.17
2.36
1.15
22.36
0.61
36.05
ATR 007
221.98
77.95
2.16
1.48
24.29
0.53
35.06
Biabasse
158.42
69.87
2.14
2.58
24.32
0.59
29.87
CTSIA 110
214.08
46.58
2.39
2.33
16.41
0.43
38.99
CTSIA 112
203.43
55.30
2.40
2.54
21.25
0.54
33.32
CTSIA 230
191.73
45.72
2.14
3.13
13.52
0.47
31.01
CTSIA 45
174.40
48.95
2.30
2.17
21.15
0.55
36.05
CTSIA 48
157.05
57.93
2.21
2.75
11.87
0.52
33.37
CTSIA 65
249.67
136.65
2.01
2.75
12.38
0.47
31.78
CTSIA 72
211.65
110.83
1.84
2.75
14.17
0.59
34.40
I91934
146.37
36.38
2.13
1.25
27.94
0.60
28.24
MM96/1751
190.92
75.80
2.32
1.00
32.43
0.66
29.25
NWA 004
200.22
69.47
2.15
1.23
22.97
0.54
34.49
Pontisange
202.22
115.20
2.19
2.58
25.77
0.57
30.94
TME 419
199.12
89.42
2.34
1.10
26.53
0.56
31.38
TME 435
188.97
79.95
2.33
1.10
26.60
0.59
31.69
UCC 2001/449
173.37
41.03
2.52
1.33
35.97
0.59
29.45
Mean
185.26
69.42
2.25
1.82
23.59
0.56
32.38
SED
9.68
6.56
0.10
0.38
3.20
0.02
3.05
Plt_ht (cm) = plant height (cm), Ht_branch (cm) = height at first branching, CMD (1–5) = cassava mosaic disease score (1–5 scale) where 1 denotes no symptom and 5 denotes highly diseased plant, and HI = harvest index.
Additive main effect and multiplicative interaction (AMMI) stability value (ASV) ranked the genotypes based on the least score (Table 7). Low scores represent the most stable genotypes. Based on the ASV, the most stable genotype for root yield was ATR 007 since it had the highest ASV ranking. 96/1708 was ranked the least stable because it had the highest ASV. In terms of root yield UCC2001/449 had the highest rank with CTSIA 48 being ranked the least. The sum of the yield and stability rankings (YSI) ranked MMM96/1751 as the genotype that combined high yield with stability. 96/1708 though high yielding was unstable because of its low rank according to the YSI. Three genotypes MMM96/1751, 00/0203, and UCC2001/449 can be considered as high yielding and stable.
Ranking of 20 cassava genotypes based on storage root yield, AMMI stability value (ASV), and yield stability index (YSI) at 12 months after planting.
Fresh storage root yield (t/ha)
Genotype
Mean
Rank (A)
IPCA[1]
IPCA[2]
ASV
ASV rank (B)
YSI (A+B)
YSI rank
MM96/1751
32.43
3
0.45
-2.26
0.54
4
7
1
Pontisange
25.76
9
-0.31
-1.43
0.37
2
11
2
00/0203
31.29
4
-1.10
-1.07
1.33
8
12
3
ATR 007
24.29
11
0.25
-1.54
0.31
1
12
4
UCC 2001/449
35.97
1
-1.48
-1.90
1.78
11
12
5
Biabasse
24.32
10
0.42
-1.88
0.50
3
13
6
96/409
27.69
6
-1.18
-0.03
1.42
9
15
7
TME 435
26.60
7
-1.19
0.66
1.44
10
17
8
I91934
27.94
5
-1.71
1.50
2.06
13
18
9
ATR 002
22.36
13
-0.49
0.96
0.59
6
19
10
NWA 004
22.79
12
-0.80
-0.18
0.96
7
19
11
CTSIA 110
16.41
16
0.45
2.10
0.54
5
21
12
96/1708
32.99
2
-2.91
0.95
3.51
20
22
13
TME 419
26.53
8
-2.21
1.12
2.66
19
27
14
CTSIA 45
21.15
15
2.01
-0.74
2.42
14
29
15
CTSIA 112
21.24
14
2.06
-0.18
2.48
16
30
16
CTSIA 65
12.38
19
1.48
0.88
1.78
12
31
17
CTSIA 230
13.52
18
2.04
1.54
2.45
15
33
18
CTSIA 72
14.17
17
2.09
0.66
2.52
17
34
19
CTSIA 48
11.87
20
2.13
0.83
2.57
18
38
20
Rt_yld = storage root yield (t/ha), IPCA [1] score = interaction principal component axis one score, IPCA [2] score = interaction principal component axis two score, ASV = AMMI stability value, and YSI = yield stability index.
Additive main effect and multiplicative interaction analysis identified four highest yielding genotypes in each of the four environments (Table 8). Four genotypes 96/1708, UCC2001/449, TME 419, and I91934 were selected for the high potential environment, Fumesua 2014. In environment FUM13, MM96/1751, UCC2001/449, 00/0203, and Biabasse were selected as the top genotypes. The top four genotypes in Nyankpala 2013 were 96/1708, 00/0203, MM96/1751, and TME 435. In the least favourable environment (Nyankpala 2014), genotypes I91934, UCC2001/449, 96/1708, and MM96/1751 were selected.
First four AMMI selections based on best root yielding genotype in each environment.
Environment
Mean storageroot yield (t/ha)
Effect
Rank
1
2
3
4
FUM14
34.82
6.283
96/1708
UCC2001/449
TME419
I91934
FUM13
25.62
-0.732
MM96/1751
UCC2001/449
00/0203
Biabasse
NIR13
17.56
-1.481
96/1708
00/0203
MM96/1751
TME 435
NIR14
16.34
-1.989
I91934
UCC2001/449
96/1708
MM96/1751
4. Discussion
The performance of cassava is subject to strong influence of genotype, environment, and genotype∗environment interaction [10–12, 16]. Significant genotypic variations were observed for growth parameters such as plant height, stem diameter, and height at first branching indicating opportunity for selection. Similar variations in response to cassava to different environments have been reported [9, 12, 35]. Fresh root yield was significantly affected by genotype and the growing environments with much higher root yields being obtained at Fumesua than Nyankpala in both years probably due to the favourable environmental conditions typical of forest agroecologies compared with the guinea savannah ecologies. This indicates that cassava responds favourably to better environments. Similar observations were made by Turyagyenda et al. [4] and Okogbenin and Fregene [9] who observed higher yields under optimum conditions than under limited moisture. Genotype ∗ environment interaction was observed over the different environments as indicated by crossover performances for some of the genotypes. This led to variations in the mean ranks of the genotypes in the different environments [8, 36]. This implies different adaptation by the different genotypes suggesting the need to identify and select location specific genotypes for different environments. Alternatively yield stability analysis can be performed to identify genotypes whose performance remains stable over several years and environments [35]. In this study, genotypes performed better in the optimum environments than under stress conditions.
Harvest index was higher in the stress environments than in the humid environments (Fumesua in 2013 and 2014). Okogbenin and Fregene [9] reported higher harvest index under limited moisture conditions than under optimum conditions. This was possibly due to the tendency of cassava to grow more vigorously and produce more above ground biomass relative to their root growth in such environments [3, 37]. This could be attributed to differential partitioning of dry matter into the above ground parts and storage roots [1, 37]. Genotype CTSIA 65 had the tallest plants but had the lowest root yield suggesting that dry matter partitioning was probably directed to top growth instead of storage root bulking. Okogbenin and Fregene [9] stated that certain cassava genotypes allocate initial dry matter production for shoot growth before storage root development commences. Genotypes that are able to balance initial dry matter accumulation between shoot and root would be ideal for moisture limited environments.
Genotypic variance accounted for a large proportion of the observed phenotypic variance for HI and plant height indicating inherent genetic variation for these traits [38]. On the contrary, Tewodros et al. [39] found lower genotypic variances for these traits indicating strong environmental influences. Though estimates of genetic components indicated large genotypic variance (σ2g) for a number of traits, their error variances were higher than the phenotypic variances in most cases implying strong environmental influence. Progress in their improvement requires evaluation in several environments to minimise the effect of G∗E interactions that may influence progress through selection [17]. Root yield and dry matter content were influenced more by their genotypic variances than the G∗E interaction variances (σ2g∗e). This indicates the presence of a significant genetic potential and therefore provides opportunity remarkable improvement with selection. Similar observations were made by Ntawuruhunga and Dixon [16] who observed large environmental variance for fresh root yield in cassava.
Broad sense heritability estimates for dry matter content and plant height were similar to the findings of Ntawuruhunga and Dixon [16]. Heritability estimates influence the amount of progress that can be made in selection for a trait of interest [40]. However, broad sense heritability alone does not always give a full indication of genetic gain that can be made through selection since it includes both additive and nonadditive components of the variation [41, 42]. It is necessary to estimate the genetic advance as a percentage of the mean to determine the actual progress. When the genetic advance as a percentage of the mean was estimated, it was found that though storage root yield had a relatively low heritability (43%) compared with plant height (82%), they had the same almost the same GAM. This suggests that plant height had greater nonadditive variance and environmental influence compared with storage root yield [41].
Stability analysis methods are often used by breeders to identify genotypes that have stable performance and respond positively to improvements in environmental conditions [24]. AMMI stability value (ASV) indicates the stability of genotypes. Genotypes having low ASV are considered more stable whilst those with high values are less stable genotypes [22]. CTSIA 110, MM96/1751, and TME 435 were the most stable for root yield. Stability alone for yield performance does not warrant selection since a consistently low yielding genotype can still be stable [43]. In some cases the most stable genotypes do not always have the best yield performance [44]. Therefore, high root yield is considered with stability in the estimation of yield stability index (YSI). The YSI which is similar to genotype stability index (GSI) proposed by Fardshadfar [26] integrates both yield and stability across environments into a single index, to select varieties. The YSI sums the rank of mean yield across environments with the rank of the ASV of genotypes [25, 27]. Genotypes with lower YSI are desirable since they combine high mean yield performance with stability [24, 25, 27, 34]. Based on the YSI, genotypes MM96/1751, Pontisange, 00/0203, ATR 007, and UCC2001/449 were selected as combining high yield performance with stability. Genotypes such as 96/409, TME 435, and I91934 had intermediate yields with moderate stability. However genotypes 96/1708 and TME 419, although high yielding, had high ASV scores resulting in high YSI scores. However they can be recommended for specific environments where they performed well. The range of variation in the favourable environments (Fumesua 2013 and 2014) was larger than in the poor environments indicating that genotypes were better able to exploit their full potential yield in the good environments [15].
5. Conclusion
Genotype∗environment interaction was significant for most of the traits indicating the need to test the genotypes in multiple environments before effective selection can be made. The high broad sense heritability estimates obtained for plant height and height at first branching suggest that significant progress can be made through selection for improvement of genotypes for these traits. Six genotypes (UCC2001/449, 96/1708, MM96/1751, 00/0203, 96/409, and I91934) had significantly higher root yield than the best local farmer preferred genotype, Pontisange. These genotypes can be evaluated in more environments to assess their adaptability and possible recommendation for release to farmers for cultivation.
Conflicts of Interest
The authors declare no conflicts of interest regarding the publication of this paper.
Acknowledgments
The authors are grateful to the Generation Challenge Programme and the Alliance for Green Revolution in Africa for sponsoring this study. The staff of the Root and Tuber Crops Improvement Sections of the CSIR-Savanna Agricultural Research Institute, Nyankpala, and the CSIR-Crops Research Institute are also deeply acknowledged for their efforts and assistance during the study period.
El-SharkawyM. A.Cassava biology and physiology200456448150110.1007/s11103-005-2270-72-s2.0-14744273154AinaO. O.DixonA. G. O.AkinrindeE. A.Effect of soil moisture stress on growth and yield of cassava in Nigeria20071018308530902-s2.0-3534886857810.3923/pjbs.2007.3085.3090El-SharkawyM. A.Physiological characteristics of cassava tolerance to prolonged drought in the tropics: Implications for breeding cultivars adapted to seasonally dry and semiarid environments20071942572862-s2.0-4624910628610.1590/S1677-04202007000400003TuryagyendaL. F.KizitoE. B.FergusonM.BagumaY.AgabaM.HarveyJ. J.OsiruD. S.Evaluation of Ugandan cassava germplasm for drought tolerance20135321222610.1093/aobpla/plt007HaldavankarP. C.JoshiG. D.BhaveS. G.KlandekarR. G.SawantS. S.Stability of yield and yield attributing phenotypic characters in sweet potato20093512835CrossaJ.CorneliusP. L.Sites regression and shifted multiplicative model clustering of cultivar trial sites under heterogeneity of error variances19973724064152-s2.0-003089237410.2135/cropsci1997.0011183X003700020017xMkumbiraJ.MahunguN. M.GullbergU.Grouping locations for efficient cassava evaluation in Malawi20033921671792-s2.0-034451168310.1017/S0014479702001199DixonA. G.NukenineE. N.Genotype * environment interaction and optimum resource allocation for yield and yield components of cassava20008111010.4314/acsj.v8i1.27711OkogbeninE.FregeneM.Genetic analysis and QTL mapping of early root bulking in an F1 population of non-inbred parents in cassava (Manihot esculenta Crantz)20021061586610.1007/s00122-002-1068-02-s2.0-0036984424AinaO. O.DixonA. G. O.AkinrindeE. A.Trait association and path analysis for cassava genotypes in four agroecological zones of Nigeria2007757597642-s2.0-3494889159110.3923/jbs.2007.759.764EgesiC. N.IlonaP.OgbeF. O.AkorodaM.DixonA.Genetic variation and genotype * environment interaction for yield and other agronomic traits in cassava in Nigeria20079941137114210.2134/agronj2006.02912-s2.0-34447555895AinaO. O.DixonA. G. O.IllonaP.AkinrindeE. A.G * E interaction effects on yield and yield components of cassava (landraces and improved) genotypes in the savannah regions of Nigeria200981949334945FalconerD. S.MackayT. F. C.1996New York, NY, USALongmanSsemakulaG.DixonA. G. O.Genotype * environment interaction, stability and agronomic performance of carotenoid-rich cassava clones200729390399PrzystalskiM.OsmanA.ThiemtE. M.RollandB.EricsonL.ØstergårdH.LevyL.WolfeM.BüchseA.PiephoH.-P.KrajewskiP.Comparing the performance of cereal varieties in organic and non-organic cropping systems in different European countries200816334174332-s2.0-5024916369810.1007/s10681-008-9715-4NtawuruhungaP.DixonA. G. O.Quantitative variation and interrelationship between factors influencing cassava yield20102615941602OtooE.OkonkwoC. R.AsieduR.Stability studies of hybrid yam genotypes in Ghana200641234238EberhartS. A.RussellW. A.Stability parameters for comparing varieties196663640GauchH. G.Model selection and validation for yield trials with interaction198844370571510.2307/2531585Zbl0707.62236van EeuwijkF. A.Linear and bilinear models for the analysis of multi-environment trials: I. An inventory of models1995841172-s2.0-002883403510.1007/BF01677551R. Bernardo2010Woodbury, Minn, USAStemma Press10.2135/cropsci2003.1578HagosH. G.AbayF.AMMI and GGE biplot analysis of bread wheat genotypes in the northern part of Ethiopia201311218PurchaseJ. L.HattingH.van DeventerC. S.Genotype * environment interaction of wheat in South Africa: stability analysis of yield performance200017310110710.1080/02571862.2000.10634878FarshadfarE.MahmodiN.YaghotipoorA.AMMI stability value and simultaneous estimation of yield and yield stability in bread wheat (Triticum aestivum L.)2011513183718442-s2.0-84855403733TumuhimbiseR.MelisR.ShanahanP.KawukiR.Genotype * environment interaction effects on early fresh storage root yield and related traits in cassava20142532933710.1016/j.cj.2014.04.008FarshadfarE.Incorporation of AMMI stability value and grain yield in a single non parametric index (GSI) in bread wheat200811141791179610.3923/pjbs.2008.1791.17962-s2.0-54049118084BarakiF.TsehayeY.AbayF.AMMI analysis of genotype * environment interaction and stability of sesame genotypes in Northern Ethiopia201413417818310.3923/ajps.2014.1781832-s2.0-84925683767EPANational action programme to combat drought and desertification2003EkanayakeI. J.Cassava procedures for growth analysis199624MsikitaW.JamesB.NnoduE.LeggJ.WydraK.OgbeF.2000IITA publicationPayneR. W.MurrayD. A.HardingS. A.BairdD. B.SoutarD. M.200912thHemel Hempstead, UKVSN International10.1080/02674649366780421PadiF. K.Genotype * environment interaction and yield stability in a cowpea-based cropping system20071581-2112510.1007/s10681-007-9420-82-s2.0-35348838992SinghK.ChaudharyB.1979Kalyani publishers/Lyall Bk Depot10.2307/2530404BoseL. K.JambhulkarN. N.PandeK.SinghO. N.Use of AMMI and other stability statistics in the simultaneous selection of rice genotypes for yield and stability under direct-seeded conditions2014741192-s2.0-8490165689410.4067/S0718-58392014000100001MutegiR. W.2009Georg-August University GottingenMalosettiM.RibautJ.-M.van EeuwijkF. A.The statistical analysis of multi-environment data: Modeling genotype-by-environment interaction and its genetic basis20134, article 4410.3389/fphys.2013.000442-s2.0-84879467512AlvesA. A. C.HillocksR. J.ThreshJ. M.BellottiA.Cassava Botany and Physiology2002Wallingford, UKCABI679010.1079/9780851995243.0000TsegayeE.DechassaN.SastryE. V. D.Genetic variability for yield and other agronomic traits in sweet potato20076194992-s2.0-3424714694010.3923/ja.2007.94.99TewodrosM.GetachewW.TadesseB.TesfayeW.Genetic diversity of cassava (Manihot esculenta Crantz) genotypes in Ethiopia20133963664210.15580/GJAS.2013.3.061013658SabesanT.SureshR.SaravananK.Genetic variability and correlation for yield and grain quality characters of rice grown in coastal saline low land of Tamil Nadu200915659PradeepkumarT.BastianD.JoyM.RadhakrishnanN. V.AipeK. C.Genetic variation in tomato for yield and resistance to bacterial wilt200139157158SinghA. K.SharmaP.SinghP. K.Studies on genetic characteristic of upland rice (Oryza sativa L.)20136451552010.5958/j.2230-732X.6.4.025YanW.TinkerN. A.Biplot analysis of multi-environment trial data: Principles and applications20068636236452-s2.0-3424724787310.4141/P05-169OliveiraE.GodoyI.Pod yield stability analysis of runner peanut lines using AMMI20066431031710.12702/1984-7033.v06n04a09