Dear Label_Name,

Greetings! I hope you are well.

I am writing on behalf of an international publishing house, Science and Education Publishing, USA (http://www.sciepub.com).

Science and Education Publishing (SciEP) is a multidisciplinary open access publisher of Journals covering the fields of Science, Technology and Medicine. All SciEP publications are subject to high-quality peer review, editorial and production processes and are freely available online.

Welcome to the April, 2013 edition of SciEP newsletter, providing monthly summaries of articles on science & education related topics of interest to scientist and professionals. The mission of our journals is to provide an open access forum for publishing new findings on issues pertinent to scientific and academic field. We also invite you to submit original research papers, review articles, short reports, or case studies are considered for publication in our Multidisciplinary Journals.

American Journal of Public Health Research
DOI: 10.12691/ajphr-1-1-3

The Implementation of a Tele-Homecare System with Service Oriented Architecture and HL7 Message Transmission Standard
Ching-Sung Wang

American Journal of Medical and Biological Research
DOI: 10.12691/ajmbr-1-1-2

The Influence of Ethanol on Pyruvate Kinases Activity in Vivo, in Vitro, in Silico
Lelevich Sergey Vladimirovich, Khrustalev Vladislav Victorovich, Barkovsky Eugene Victorovich, Shedogubova Tatyana Aleksandrovna

American Journal of Medical Sciences and Medicine
DOI: 10.12691/ajmsm-1-1-1

Clinical and Histological Features in Wegener Granulomatosis
Li Cavoli Gioacchino, Passantino Rita, Bono Luisa, Tortorici Calogera, Ferrantelli Angelo, Giammarresi Carlo, Rotolo Ugo

American Journal of Infectious Diseases and Microbiology
DOI: 10.12691/ajidm-1-1-2

The Curse of Lymphatic Filariasis: Would the Continual Use of Diethylcarbamazine Eliminate this Scourge in Papua New Guinea?
Anne Lanham, Lillian Mwanri

We have released issues of some journals on February 2013.
American Journal of Public Health Research
American Journal of Medical Sciences and Medicine
American Journal of Medical and Biological Research
American Journal of Pharmacological Sciences
American Journal of Clinical Medicine Research
American Journal of Nursing Research
American Journal of Cancer Prevention
American Journal of Medicine Studies
American Journal of Infectious Diseases and Microbiology
American Journal of Cardiovascular Disease Research
American Journal of Biomedical Research
American Journal of Mechanical Engineering
American Journal of Microbiological Research
Journal of Preventive Medicine
American Journal of Food and Nutrition
American Journal of Hypertension Research
American Journal of Sports Science and Medicine
Journal of Cancer Research and Treatment
American Journal of Zoological Research
Journal of Biomedical Engineering and Technology
American Journal of Microbiological Research
Biomedical Science and Engineering

We are pleased to announce that submissions are now being accepted for peer review for all the journals across all our subject areas in SciEP.

Please submit your manuscript via our online submission system:
http://www.sciepub.com/author/signin

To unsubscribe from our mailing list, please click on the following link: unsubscribe

I am looking forward to hearing from you.

Kind regards,

David Rain
Email:join@sciepub.com
Science and Education Publishing (http://www.sciepub.com)

Genotype x Environment Interactions on Seed Yield of Inter-racial Common Bean Lines in Kenya

Jean M. Mondo*, Paul M. Kimani, Rama D. Narla

Department of Plant Science and Crop Protection, University of Nairobi, P.O. Box 29053-00625, Nairobi, Kenya

Abstract: Determination of yield stability is critical in identifying new common bean cultivars with either specific or broad adaptation in target environments. This study aimed to assess genotype by environment (G x E) effects on agronomic performance of 78 F1.7 lines selected with molecular markers for multiple disease resistance from 16 inter-racial bean populations. Field trials were conducted in low-, medium- and high altitude conditions in Kenya. Data collected on seed yield were subjected to additive main-effects and multiplicative interaction (AMMI) model to separate additive variance from the G x E interaction and to determine the stability of genotypes across locations. Results showed that G x E effects were highly significant (P<0.001), implying that tested lines behaved differently across the three locations. Better yields were recorded from high altitude Tigoni site while the lowest were from low altitude Mwea site. Yield across sites ranged from 1,518 to 2,748; 1,324 to 3,860; 1,537 to 3,722 and 1,010 to 3,718 kg ha-1 for pinto, red mottled, red kidney and mixed color bean lines, respectively. Number of pods plant⁻¹ was the most strongly correlated to seed yield and could be, therefore, used as an indirect selection criterion for seed yield. The environment was responsible for the largest part of yield variability (86.4%, 84.8%, 82.3% and 49.5% for pinto, red kidney, red mottled and mixed color bean lines, respectively). KMA13-22-21 and KMA13-29-21 were the most stable high yielding lines across locations. Higher yielding lines were the most unstable across sites. Two pinto, four red kidney, 15 red mottled, and two mixed color lines did better than their corresponding checks with yield advantages of 7.6, 14.3, 71.5, and 34.9%, respectively. These lines should, therefore, be selected for further testing and release.

Keywords: Phaseolus vulgaris, inter-racial crosses, gamete selection, market class, AMMI model

Received March 01, 2019; Revised April 14, 2019; Accepted May 05, 2019

Citation: Jean M. Mondo, Paul M. Kimani, Rama D. Narla, Genotype x Environment Interactions on Seed Yield of Inter-racial Common Bean Lines in Kenya. World Journal of Agricultural Research, 7 (3). 76-87. January 2019. doi:10.12691/wjar-7-3-1

1. Introduction

Common bean (Phaseolus vulgaris L.) is the most important grain legume consumed worldwide. Six races of common bean have been distinguished based on morphological, agronomic, adaptive, and molecular characteristics [1]. Three of these races (Durango, Jalisco, and Mesoamerica) belong to the Middle American also referred to as Mesoamerican gene pool. The other three races (Chile, Nueva Granada, and Peru) belong to the Andean South American gene pool [1,2,3]. Small-seeded beans (<25 g 100-seed mass) belong to Mesoamerica race and are well adapted to relatively warmer tropical lowlands. Medium-seeded beans have a 100-seed mass of 25 to 40 g and belong to Durango race for the semi-climbers and Jalisco race for the climbers. They are as well adapted to tropical and subtropical environments [4]. Small- and medium-seeded beans often have indeterminate growth habit and out-yield their large-seeded counterparts (>40 g 100-seed mass) from Chile and Nueva Granada races by as much as 500 to 2000 kg ha⁻¹ [1,5]. In addition to their high yield potential, small- and medium-seeded beans are resistant to several diseases devastating the large-seeded beans such as angular leaf spot, anthracnose, rust, bean golden yellow mosaic virus and bean common mosaic virus and possess genes and high level of resistance to drought stress [6]. However, large-seeded Andean beans are the most widely grown in Eastern Africa because they are preferred by farmers and consumers for their seed quality and often fetch higher prices [5,7]. Andean genotypes possess a narrow to moderate genetic base especially for disease resistance and yield potential and thus, threatening progress toward improvements for those traits [8]. The genetic base of common bean varieties grown in Eastern Africa needs to be broadened to enhance yield potential and resistance to diseases. Inter-racial and inter-gene pool crosses provide an important opportunity to create useful genetic variation for maximizing gains from selection, broadening the genetic base of commercial cultivars and making efficient use of available resources [4,9]. Despite of the limitations faced in developing inter-racial and inter-gene pool cultivars (notably the F₁ hybrid dwarfism, weakness, or incompatibility, problems in recovering desirable seed quality and adaption characteristics, and cripples or virus-like foliage symptoms), breeding programs across the world have succeeded through inter-racial and inter-gene crosses to develop genotypes combining desirable traits such as tolerance to production limiting factors (especially diseases, drought), seed quality and high yield potential [10].

Common bean yields in Eastern and Central Africa are among the lowest in the world (0.5 t ha⁻¹) while the potential is between 1.5 to 3 t ha⁻¹ for bush beans, and up to 5 t ha⁻¹ for climbing beans [11,12,13]. Diseases are among the major factors constraining bean production in the region, and therefore, the improvement of bean productivity requires effective and efficient selection for yield traits along with controlling major diseases [8,14,15]. Plant breeding is the most cost-effective and sustainable approach to cope with bean diseases since no additional investment is required from farmers. However, as several pathogens co-infect beans at the farmer level, breeding for one trait would not result in a significant change, and thus, a multiple disease resistance approach should be promoted for a more and durable impact on yield and farmers’ livelihood.

In bean breeding programs, a large number of genotypes are tested for many generations within contrasting environments before release for seed multiplication and distribution to growers. Because environmental conditions for testing are distinct, the genotype by environment interaction (G x E) affects the agronomic performance, and thus, it is necessary to analyze the stability of genotypes across environments [16,17]. This allows the assessment of the real impact of selection and ensures high reliability in the genotype recommendation for a specific place or environment groups [18,19]. Another key reason for the G x E analyses in bean breeding in Africa is that lines adapted to an African bean environment (AFBE) can be grown in similar areas in other parts of Africa [20]. Due to differences among growing regions, breeding might be more effective if it was AFBE based. Therefore, we hope that lines developed through the current breeding program in the AFBE of Kenya could be adapted and disseminated in African areas with similar agro-ecological conditions.

The specific objective of this study was to assess the G x E effects on the agronomic performance of 78 advanced F1.7 lines previously selected for multiple disease resistance using molecular markers. These lines originated from 16 small- and medium-seeded inter-racial bean populations, which were subsequently grouped in four market classes.

2. Material and Methods

2.1. Experimental Sites

This study was conducted in three different agro-ecological zones, representing major bean growing environments in Kenya. Experiments were conducted during 2017 short rain season at KALRO (Kenya Agricultural and Livestock Research Organization)-Mwea representing low altitude conditions, Kabete Field Station of the University of Nairobi, the medium altitude, and KALRO-Tigoni for the high altitude environments. KALRO-Mwea is located on coordinates 0°38’S (latitude); 37°22’E (longitude) and at approximately 1150 masl. This research station receives mean precipitation of 850 mm per year with a bimodal distribution. The long rain season starts in March and ends in May. The short rain season usually starts in October to end in late December. Mean annual temperatures range from 15.6°C to 28.6°C. Soils at this station are vertisols with an acidic pH of about 5.1 [21]. KALRO-Tigoni is located at coordinates 01°08’S; 036°40’E and at approximately 2130 masl. It receives bimodal rainfall of 1100 mm per year. Temperatures range from 12°C to 24°C. Soils at Tigoni are humic nitisols with soil pH of approximately 4.6 [22]. Kabete Field Station is located at 01°15’S; 036°44’E and 1820 masl. The station experiences mean bimodal precipitation of 1059 mm per year. Mean monthly temperatures range between 12.3°C and 22.5°C [23]. Soils are humic nitisols, very deep, well-drained, friable clay with acid humic topsoil, and dark reddish brown in color. The pH is about 5.0 to 5.4 and a mean sunshine of 6.6 hours per day. Following Wortmann and Allen (1994) classification, Kabete and Tigoni are located in the African Bean Environment I (AFBE 1) while Mwea is in the AFBE 8.

2.2. Plant Materials

Study materials were 78 lines including 73 F1.7 bean lines selected from 16 inter-racial populations, and five check varieties (Mex54, AND1062, BRB191, GLP92, and KATB1). The 73 F1.7 lines were grouped in four market classes on the basis of their seed color, shape and size. The 73 F1.7 bean lines comprised of, 14 red kidney, 16 red mottled, 12 pinto and 31 mixed color bean lines. The four market classes were evaluated in separate trials and compared with appropriate checks selected among parental genotypes. In these trials, AND1062 and Mex54 were used as checks for red kidney, BRB191 for red mottled, GLP92 for pinto, and KATB1 for mixed color lines. These plant materials were developed following the gamete selection procedure as summarized in Table 1.

2.3. Experimental Design and Crop Management

A simple lattice experimental design with four replicates was used for each market class depending on the number of tested lines; a 4 x 4 lattice design for red kidney, red mottled and pinto and a 6 x 6 lattice for mixed color market class. A plot consisted of three 4m rows. Seed rate was 10 seeds m⁻¹ spaced by 0.2 m within rows and 0.5 m between rows. Two guard rows were erected to avoid competition and interference between genotypes. All the field experiments were planted in October 2017 during the short rain season. Diammonium phosphate (DAP) at a rate of 80 kg ha⁻¹ was applied at planting. Weeding at all sites were carried out three times: two weeks after seedling emergence, before flowering and after podding. The pesticide Confidor (200 g l⁻¹ Imidacloprid) was used to control whiteflies and leafminers.

Table 1. Breeding scheme for the development of bean lines evaluated in three agro-ecological conditions in Kenya

2.4. Data Collection and Analysis

Data were collected on seedling emergence rate, plant vigor, days to flowering, growth habit, days to maturity, number of pods per plant, number of seeds per pod, 100-seed mass, grain yield, and harvest index using the standard system for the evaluation of bean germplasm as described by [24]. Statistical analyses were performed using GenStat 17th edition [25] and Statistix 8.0 version [26]. Combined analysis of variance (ANOVA) was conducted to determine the magnitude of variation associated with each source (environment, genotype and their interaction). Fisher’s least significant difference (LSD) test was used for separation of means at 5% probability level. ANOVA is an additive model in which the G x E interaction is a source of variation, but its intrinsic effects are not analyzed. The additive main effect and multiplicative interaction (AMMI) model was, therefore, necessary to separate the additive variance from the G x E interaction [27,28]. In fact, AMMI uses ANOVA to test the main effects of genotypes and environments, and principal component analysis (PCA) to analyze the residual multiplicative interaction between genotypes and environments to determine the sum of squares of the G × E interaction, with a minimum number of degrees of freedom [29]. The AMMI model used was:

[picture] | (11)

Where: Yger is the yield of genotype g in the environment e for replicate r; μ is the grand mean; αg is the genotype mean deviations; βe is the environment mean deviation; n is the number of PCA axes retained in the model, λn is singular value for PCA axis n; ygn is the genotype eigenvector values for PCA axis n; δen is the environment eigenvector values for PCA axis n; ρge represents the residuals and εger is for error.

AMMI analysis was also used to determine the stability of the genotypes across locations using the PCA scores (IPCA1 and IPCA2). The IPCA score near zero reveals more stable genotypes, while large values indicate more responsive and less stable genotypes. AMMI stability value for the grain yield was estimated as shown as follows [30]:

Where: ASV is the AMMI stability value, SS IPCA 1 and SS IPCA 2 are the sum of squares of IPCA 1 and 2, respectively and IPCA is the interaction principal component analysis. Thus, lowest ASV indicates a wide adaptation of specific genotypes for certain environments and vice-versa.

Genotype main effect plus genotype by environment interaction (GGE) biplots were subsequently constructed to determine adaptation and stability of genotypes across test environments. From this analysis, genotypes located near the biplot origin were considered as widely adapted, while genotypes located far were specifically adapted. All the genotypes with positive IPCA1 scores responded positively to the environment having positive IPCA1 scores, and were, therefore, adapted to that particular environment [31,32].

3. Results

3.1. ANOVA of Main Effects and Multiplicative Interaction (AMMI)

Analysis of the main effects and multiplicative interaction (AMMI) for inter-racial bean lines showed that the effects on seed yield due to genotypes (G), environments (E) and interactions between genotypes and environments (G x E) were significant (P<0.01) regardless of the market class. Treatments (G, E, and G x E) contributed up to 83.7% to the total variability for pinto, 90.8% for red kidney, 91.3% for red mottled, and 93.3% for mixed color bean lines. When partitioning the treatment variability, the environment contributed most to the variance (86.4% for pinto, 84.8% for red kidney, 82.3% for the red mottled, and 49.5% for mixed color lines). The variability due to interaction between genotypes and environments was high for mixed color bean lines (26.7%). This high contribution of variability due to the interaction between genotypes and environments suggests that test lines were not stable and thus responded differently across locations and should, therefore, be selected and recommended to specific environments. IPCA1 contributed the most to the G x E effects accounting for more than 80% of the variability regardless of the market class, suggesting a high contribution of the genotypes in the interaction (Table 2, Table 3, Table 4, and Table 5).

Table 2. Summary of ANOVA for Additive Main effects and Multiplicative Interaction (AMMI) for seed yield (kg ha⁻¹) of pinto bean lines grown at three locations in Kenya

Table 3. Summary of ANOVA for Additive Main effects and Multiplicative Interaction (AMMI) for seed yield (kg ha⁻¹) of red kidney bean lines grown at three locations in Kenya

Table 4. Summary of ANOVA for Additive Main effects and Multiplicative Interaction (AMMI) for seed yield (kg ha⁻¹) of red mottled bean lines grown at three locations in Kenya

Table 5. Summary of ANOVA for Additive Main effects and Multiplicative Interaction (AMMI) for seed yield (kg ha⁻¹) of mixed color bean lines grown at three locations in Kenya

3.2. Stability Analysis

3.2.1. Pinto Bean Lines

The AMMI model showed that the highest seed yields of pinto bean lines across sites were recorded at Tigoni in the high altitude (4,347 kg ha⁻¹), followed by Kabete in medium altitude (1,388 kg ha⁻¹) whereas the lowest yields were from Mwea located in the low altitude (585.6 kg ha⁻¹) (Table 6). Across sites, the genotypes KMA13-22-21 (P5) and KMA13-22-30 (P6) were the best yielding with 2,748 kg ha-1 and 2,726 kg ha⁻¹, respectively, but not significantly different from the check variety GLP92 which yielded 2,543 kg ha⁻¹. All the other lines were either statistically equal or inferior to the check variety.

The AMMI stability value (ASV) of pinto bean lines showed that the check variety GLP92 was the most stable across sites (ASV=3.5). Among advanced lines, KMA13-24-6 (P11) and KMA13-21-10 (P1) were the most stable genotypes across sites with ASV of 15.6 and 45.8, respectively. KMA13-22-30 (P6) was the least stable across sites (ASV=816.7). The first four AMMI selections per environment were KMA13-21-10 (P1), GLP92, KMA13-21-19 (P2) and KMA13-22-21 (P5) for low altitudes; GLP92, KMA13-22-21 (P5), KMA13-23-18 (P9) and KMA13-23-13 (P8) for medium altitudes and KMA13-22-30 (P6), KMA13-22-21 (P5), KMA13-23-22 (P10) and KMA13-24-7 (P12) for high altitudes (Table 6).

3.2.2. Red Kidney Bean Lines

For the red kidney bean lines (Table 7), the highest seed yields across sites were recorded at Tigoni (4,642 kg ha⁻¹), much higher than Kabete (1,238 kg ha⁻¹) and Mwea (954 kg ha⁻¹). Mex54 with a mean of 3,722 kg ha⁻¹ out-yielded all the advanced bean lines and the other check variety AND1062 which yielded 2,266 kg ha⁻¹. The best genotype among the advanced red kidney bean lines was KMA13-30-22 (RK13) with a seed yield of 3,226 kg ha⁻¹ which was higher than all test lines and one of the check varieties, AND1062. It was not, however, significantly different from the best check variety (Mex54). Most of the red kidney bean lines were bush (Type I and II). This could explain why the cultivar Mex54 which is Type III growth habit (i.e semi-climber) had a significantly higher yield than most of the red kidney bean lines.

Table 6. Seed yield (kg ha⁻¹), ranking (in parenthesis), IPCA scores and AMMI stability values (ASV) of advanced pinto bean lines grown at three locations in Kenya

Table 7. Seed yield (kg ha⁻¹), ranking (in parenthesis), IPCA scores and AMMI stability values (ASV) of red kidney bean lines grown at three locations in Kenya

KMA13-19-16 (RK3), KMA13-25-20 (RK7), KMA13-20-3 (RK4) and AND1062 were the most stable genotypes across sites with ASV of 1.5, 1.8, 5.6 and 9.5, respectively. The high yielding genotypes KMA13-30-22 (RK13) and Mex54 were the least stable across sites. The first four AMMI selections per environment were KMA13-21-11 (RK5), KMA13-30-22 (RK13), Mex54, and AND1062 for medium altitude areas such as Kabete; Mex54, KMA13-26-32 (RK8), KMA13-27-31 (RK9), KMA13-25-3 (RK6) for low altitude agro-ecological zones such as Mwea, and Mex54, KMA13-30-22 (RK13), KMA13-28-2 (RK10), and KMA13-19-16 (RK3) for high altitude zones such as Tigoni.

3.2.3. Red Mottled Bean Lines

Across environments, the highest yields for the red mottled lines were recorded from Tigoni (4,703 kg ha⁻¹), followed by Kabete (1,358 kg ha⁻¹). The lowest means were from Mwea (551 kg ha⁻¹). KMA13-29-21 (RM13) with a mean seed yield of 3,860 kg ha⁻¹ out-yielded all the advanced red mottled lines and the check variety BRB191 which recorded a mean yield of 1,352 kg ha⁻¹. Among the red mottled bean lines, only KMA13-24-11 (RM6) yielded lower than the check variety but the difference was not significant (1,324.0 kg ha⁻¹) (Table 8).

Table 8. Seed yield (kg ha⁻¹), ranking (in parenthesis), IPCA scores and AMMI stability values (ASV) of red mottled bean lines grown at three locations in Kenya

Table 9. Seed yield (kg ha⁻¹), ranking (in parenthesis), IPCA scores and AMMI stability values (ASV) of mixed color bean lines grown at three locations in Kenya

KMA13-20-14 (RM3) and KMA13-24-16 (RM3) were the most stable genotypes across environments with ASV scores of 5.8 and 6.0, respectively. However, the best yielding line, KMA13-29-21 (RM13), was also the least stable across environments. The first four AMMI selections per environment were KMA13-24-11 (RM16), KMA13-27-25 (RM10), KMA13-29-21 (RM13) and KMA13-29-24 (RM14) for the low altitude agro-ecological zones; KMA13-29-21 (RM13), KMA13-17-17 (RM17), KMA13-32-24 (RM15), KMA13-17-25 (RM1) for the medium altitudes and KMA13-29-21 (RM13), KMA13-24-17 (RM8), KMA13-24-5 (RM5) and KMA13-29-24 (RM14) for the high altitudes (Table 8).

3.2.4. Mixed Color Bean Lines

The best yields for mixed color bean lines were obtained from Tigoni (2,550 kg ha⁻¹), followed by Kabete (1,797 kg ha⁻¹). Mwea with a mean grain yield of only 742 kg ha⁻¹ was the least productive site. Among test lines, KMA13-28-21 (MC28), a black-seeded line out-yielded all the other lines and the check varieties with a mean seed yield of 3,718 kg ha⁻¹. The other high performing lines included KMA13-27-27 (MC10) with a yield of 2,845 kg ha⁻¹, KMA13-21-20 (MC32) (2,329 kg ha⁻¹), KMA13-27-12 (MC27) (2,044 kg ha⁻¹) and KMA13-23-20 (MC5) (2,017 kg ha⁻¹). The lowest yielding line was KMA13-27-1 (MC31). This line characterized by greyish green seeds had a mean yield of 1,010 kg ha⁻¹, which was lower than the greyish green-seeded check variety KATB1 which yielded 1,144 kg ha⁻¹.

The most stable lines across sites were KMA13-23-20 (MC5) (ASV score of 5.7) and KMA13-22-322 (MC15) (ASV score of 8.3). KMA13-27-27 (MC10) (ASV score of 152.2) was the least stable across environments. The first four AMMI selections per environment were KMA13-28-5 (MC11), KMA13-31-62 (MC18), KMA13-23-9 (MC4) and KMA13-27-27 (MC10) for low altitudes; KMA13-28-21 (MC28), KMA13-21-20 (MC32), KMA13-27-12 (MC27), KMA13-28-13 (MC12) for medium altitudes, and KMA13-28-21 (MC28), KMA13-27-27 (MC10), KMA13-22-23 (MC21) and KMA13-21-20 (MC32) for high altitude bean growing environments (Table 9).

Figure 1. Polygon view of GGE biplots for best pinto (A), red kidney (B), red mottled (C), and mixed color (D) bean genotypes for seed yield across three environments in Kenya

3.3. GGE Biplots for G x E Analysis “Which Won Where” of the Inter-racial Bean Lines across Three Locations in Kenya

Most of the variability was explained by the 2 PCs regardless of the market class (97.6%, 94.5%, 95.9%, and 92% for pinto, red kidney, red mottled and mixed color bean lines, respectively). PC1 contributed the most to that variability (92%, 86.4%, 87.4%, and 81.8% for pinto, red kidney, red mottled and mixed color genotypes, respectively). Tigoni, the high altitude site, was the best environment for most of the genotypes regardless of the market class. The variability across environment was high for the red mottled genotypes for which there are three distinct mega-environments; genotypes having performed differently in each site. The variability across environments was low for mixed color genotypes for which there is only one mega-environment suggesting that better yielding genotypes in one site were the better in the other two environments. From graphs, pinto bean lines KMA13-22-21 (P5) and GLP92 (P13) performed best at Kabete and Mwea while KMA13-22-30 (P6) was best for Tigoni. Among red kidney lines, Mex54 (RK15) was best for Mwea while KMA13-21-11 (RK5) and KMA13-30-22 (RK13) did better at Kabete. Red mottled bean lines KMA13-17-17 (RM17) and KMA13-24-5 (RM5) were suited for Kabete and KMA13-24-17 (RM8) for Mwea. KMA13-27-27 (MC10) and KMA13-28-21 (MC28) were the best mixed color bean lines for Mwea and Tigoni whereas KMA13-21-20 (MC32) and KMA13-27-12 (MC27) won at Kabete (Figure 1). From the GGE biplots, Tigoni was the most discriminative as it was far from the origin of the biplot graph regardless of the market classes. All genotypes inside the polygon, mainly those located close to the plot origin were less responsive than the vertex genotypes and not the best in any environment.

3.4. Pearson’s Correlation Coefficients among Yield and Yield Components of Inter-racial Bean Lines Grown in Three Locations in Kenya

Table 10, Table 11, Table 12, and Table 13 present the Pearson’s correlation coefficients among seed yield and yield-related parameters for pinto, red kidney, red mottled and mixed color bean market classes, respectively. Regardless of the market class, seed yield was positively correlated with days to flowering, days to maturity, number of pods per plant, number of seeds per pod, 100-seed mass and harvest index (P<0.05). It was, however, negatively correlated with seedling emergence rate and plant vigor (except the mixed color market class). This would imply that the higher the number of pods per plant and the higher the number of seeds per pod, the higher the yield was. Better yielding plants were late to reach the 50% flowering stage as they contained a large number of flowers which appeared progressively. This had also an impact on the days to maturity which was delayed compared to plant developing fewer flowers and fewer pods. As the plant vigor score varies from 1 to 9 [24] from which 1 is the best score and 9 the worst, the more a plant was vigorous, the more it could carry more flowers and more pods and consequently, the more the yield was higher. The negative correlation observed in this study between seed yield and seedling emergence rate could be largely due to extrapolation as the yield ha⁻¹ was estimated on the single plant basis. If the yield per m² (or per plot) was considered in extrapolation regardless of the number of plants, the relationship may change. Regardless of the market class, the number of pods per plant was the most highly correlated with the seed yield, suggesting its usefulness as an indirect selection criterion for seed yield.

Looking at the correlation between growth habit and yield and yield components such as number of pods per plant, it has been observed heterogeneity among market classes. There were negative but not significant correlations between the growth habit and the seed yield (r=-0.05ns) and between the growth habit and the number of pods per plant (r=-0.06ns) for the pinto bean lines. However, the trend was different for other market classes for which the growth habit was positively correlated with seed yield and number of pods per plant. The growth habit was positively but not significantly correlated with seed yield (r= 0.04ns) for red mottled market class. However, the correlation between the growth habit and number of pods per plant on red mottled market class was positive and significant (r=0.14*). The trend was the same on red kidney bean lines for which the correlations were significant and positive between the growth habit and the seed yield (r=0.20*) and between the growth habit and the number of pods per plant (r=0.24**). This study reflected the general assumption that yield increases with growth habit such that Type IVs (climbers) are the best yielding.

Table 10. Pearson’s correlation coefficients among seed yield and yield components of pinto bean lines grown at three locations in Kenya

Table 11. Pearson’s correlation coefficients among seed yield and yield components for red kidney bean lines grown at three locations in Kenya

Table 12. Pearson’s correlation coefficients among seed yield and yield components for red mottled bean lines grown at three locations in Kenya

Table 13. Pearson’s correlation coefficients among seed yield and yield components of mixed color bean lines grown at three locations in Kenya

4. Discussion

4.1. Agronomic Performance of Inter-racial Bean Lines across Sites in Kenya

Effects due to interactions between the sites and the genotypes for all the traits and all the market classes were significant (P<0.05), implying that advanced bean lines responded differently to environmental conditions prevailing at test sites. As a result, their ranking varied significantly across the three sites. For all the traits, crops grown at Tigoni in high altitude recorded the highest means statistically superior to the other two sites namely Kabete and Mwea located in medium and low altitudes, respectively. The better performance recorded at Tigoni could be attributed to the relatively cooler conditions offered to crops; which led to slower plant growth and delayed maturity and, therefore, longer seed filling period which resulted in higher seed yields. Similar results were reported by [5]. The low yield recorded at Mwea in low altitude could be due to dry spells and erratic rainfall observed in that site during the experiment. In fact, the mean monthly temperature was 24.3°C with a total rainfall of approximately 311.4 mm for the period of September 2017 to February 2018. In addition, more than 85% of that rainfall was recorded during October and November, flooding young bean seedlings. The most critical phases (flowering and podding) experienced a dry period as no rain was recorded in January and February 2018 (0 mm), and thus, affecting negatively the grain yield. [33,34,35,36] observed that water stress during flowering, pod filling stages severely affects the harvest index and the seed yield. Seed yield losses might exceed 20% if the stress occurs during the early vegetative growth and could reach up to 50% in the early pod filling [32,37,38]. As most of test bean lines were of indeterminate growth habit, effects of water stress in low altitude Mwea site were more pronounced compared to dwarf cultivars as also reported in Malawi by [39]. [40] demonstrated that humid high altitude conditions are more conducive to indeterminate growth habit cultivars.

Seed yield and yield related components were varying significantly among genotypes and market classes. In all the market classes, there were promising genotypes for seed yield and which performed better than corresponding commercial check varieties, apart from the red kidney market class where the best yielding line was not significantly different from the best check variety (Mex54). This was probably an effect of growth habit as Mex54 is a semi-climber cultivar while most of test red kidney lines were bush lines (Type I and Type II growth habit). However, four of the 15 advanced red kidney lines were superior to the other check variety (AND1062) which is a bush cultivar. The presence of promising lines, regardless of the market class, demonstrated the effectiveness of inter-racial crosses to improve the seed yield of common bean. After studying the effects on seed yields of the Andean intra-gene pool and Andean-Middle America inter-gene pool crosses, [5] concluded that the utilization of high yielding genotypes from both gene pools which are diverse and with positive general combining ability could maximize gains from seed yield selection. [9] and [41] had previously demonstrated the superiority of the inter-racial lines over the intra-racial, suggesting the necessity to explore them as a mean to create useful genetic variations and to broaden the genetic base of commercial cultivars as well as maximizing gains from selections.

The seed yield was high for market classes with higher 100-seed mass compared to smaller seeds. While assessing effects of size of seed grown on the growth and yield of common bean, [42] concluded that sowing larger seeds improves the early-season plant growth which is advantageous for crop establishment in stressed environments. This could explain why red kidney and red mottled market classes had higher yields than pinto and mixed color market classes. The effects of seed size on yield were much more pronounced among the lines within the same market class than among market classes. This study which had both large- and small-/medium- seeded genotypes disagrees with the general observation (especially in Colombia/CIAT) that small-seeded lines yield better than large-seeded types [5]. Another key reason is that of [43] who presented evidence that large-seeded bean lines adapt better to cooler conditions from higher elevations than small-seeded counterparts.

4.2. Correlations between Seed Yield and Yield Components of Inter-racial Bean Lines grown in Three Locations in Kenya

Seed yield was significantly and positively correlated with days to flowering, days to maturity, number of pods per plant, number of seeds per pod, 100-seed mass and harvest index. The most important of these yield components regardless of the market class was the number of pods per plant, suggesting that it can be used by plant breeders as an additional and indirect selection method for seed yield. Similar results were found by [32,36,44]. This study reflected the general assumption that yield increases with growth habit such that Type IVs (climbers) are the best yielding. In fact, this study revealed a positive correlation between the growth habit and the number of pods per plant and between the growth habit and the seed yield regardless of the market class. However, the trend was opposite for the pinto bean lines for which correlations were negative but not significantly. [45] explained that in stressed environments, these climbing genotypes possess a yield compensation capacity to recover rapidly from stress.

Better yielding lines were late to reach the 50% flowering stage as they contained a large number of flowers which appeared progressively. This had also impacted the days to maturity which was delayed compared to plant developing fewer flowers and fewer pods. These findings are similar to those of [5,9,46]. However, opposite results were found in drought stress environments where higher yield was in negative correlation with days to maturity [47,48]. There were no significant correlations between the growth habit and the duration to flowering and to maturity for all market classes. This contrasts the general assumption that climbers take longer to flower and to mature compared to bush bean lines.

Significant negative correlations were detected between seed yield and plant vigor and between seed yield and seedling emergence rate. As the plant vigor score varies from 1 to 9 [24] from which 1 is the best score and 9 the worst, the more a plant was vigorous, the more it could carry more flowers and more pods and consequently, the more the yield was higher. The negative correlation existing between yield and seedling emergence rate could be attributed to extrapolation as seed yield ha⁻¹ was estimated on the basis of single plants. If the yield per m² (or per plot) was considered in extrapolation regardless of the number of plants, the relationship may change.

4.3. Yield Stability and Genotype-environment Interaction (G x E) Effects on Seed Yield

Variability among genotypes across sites was highly significant regardless of the market class. Treatments (G, E, and G x E) contributed the most to the variability for up to 80% regardless of the market class. This showed the diversity of sites and the existence of significant genetic differences among the advanced lines for seed yield as also reported by [16] and [48]. By partitioning treatments’ contribution for every market class, the environment was responsible for the largest part of the variability. Similar results were found on common bean by [39] in Malawi and [16] and [17,50] in Ethiopia. Although the environment is a very broad term and includes many factors (predictable and unpredictable); it was the temperature and the amount and distribution of rainfall that had mainly contributed to observed results. Tigoni in high altitude experienced cooler conditions (15.8°C) with a relatively well-distributed rainfall along the growing season (506 mm). Kabete experienced mean monthly temperatures of 18.2°C and an amount of rainfall of 372 mm. Mwea in low altitude was warmer (24°C) with erratic rainfall as described previously (311 mm). Other key environmental factors (e.g. soil type, nutrients, pH, etc.) were not significantly different among the three sites.

Interaction between genotype and environment was high for mixed color market class (26.7%), suggesting that test lines were not stable and thus responded differently across locations. These genotypes should, therefore, be selected and recommended to specific environments. From ASV, higher yielding lines were also the most unstable across sites. This is supporting results found by [17,50,51] showing that the stable lines are not always the better yielding. In fact, [52] demonstrated that a satisfactory Type I stability parameter (i.e., CV) is often linked with reduced yield performance.

5. Conclusion

Promising genotypes combining high seed yield potential and high stability across environments were identified from all market classes. The environment contributed the most to the variability among lines. The high altitude Tigoni site was the best environment for bean cultivation regardless of the market classes. Although the best yielding lines were not the most stable across sites, KMA13-22-21 a pinto bean line and KMA13-29-21 a red mottled line combined high yield potential and wider adaptation across the three agro-ecological conditions. Two pinto, four red kidney, 15 red mottled, and two mixed color bean lines did better than their corresponding checks with yield advantages of 7.6, 14.3, 71.5, and 34.9%, respectively. These lines should, therefore, be selected for further testing and release.

Acknowledgments

Financial support received from the Université Evangélique en Afrique to conduct this research is gratefully acknowledged. Authors are thankful to the Center Directors of KALRO Mwea (Dr. John Kimani) and KALRO Tigoni (Dr. Moses W. Nyongesa) and the Field Station Manager of the University of Nairobi for providing lands and facilities for this research to be successfully conducted.

Statement of Competing Interests

The authors do not have any competing interests.

References

[1]   | Singh, S.P., P. Gepts & D. Debouck (1991). Races of common bean (Phaseolus vulgaris Fabaceae). Economic Botany 45(3): 379-396.
[2]   | Beebe, S., A.V. Gonzalez & J. Rengifo (2000). Research on trace minerals in the common bean. Food and Nutrition Bulletin 21(4): 387-391.
[3]   | Kwak, M., O. Toro, D. Debouck & P. Gepts (2012). Multiple origins of the determinate growth habit in domesticated common bean (Phaseolus vulgaris L.). Annals of Botany 110(8):1573-1580.
[4]   | Singh, S.P. (2001). Broadening the genetic base of common bean cultivars: A Review. Crop Science 41(6): 1659-1675.
[5]   | Singh, S.P, H. Teran, C.G. Muñoz & J.M. Osorno (2002). Selection for seed yield in Andean intra-gene pool and Andean × Middle American inter-gene pool populations of common bean. Euphytica 127(3): 437-444.
[6]   | Terán, H. & S.P. Singh (2002). Comparison of sources and lines selected for drought resistance in common bean. Crop Science 42(1): 64-70.
[7]   | Sichilima, T., L. Mapemba & G. Tembo (2016). Drivers of dry common beans trade in Lusaka, Zambia: A trader’s perspective. Sustainable Agriculture Research 5(2): 15.

[8]  
Kimani, P.M., R. Buruchara, K. Ampofo, M. Pyndji, R. Chirwa & R. Kirkby (2005). Breeding beans for smallholder farmers in Eastern, Central and Southern Africa: Constraints, achievements and potential. In: Pan-African Bean Research Network (PABRA) Millennium Workshop, 28 May - 1 June, Arusha, Tanzania. pp. 11-28.

[9]   | Welsh, W., W. Bushuk, W. Roca & S.P. Singh (1995). Characterization of agronomic traits and markers of recombinant inbred lines from intra- and interracial populations of Phaseolus vulgaris L. Theoretical and applied genetics 91(1): 169-177.
[10]   | Kelly, J.D. & M.W. Adams (1987). Phenotypic recurrent selection in ideotype breeding of pinto beans. Euphytica 36(1): 69-80.
[11]   | Kaizzi, K.C., J. Byalebeka, O. Semalulu, I.N. Alou, W. Zimwanguyizza, A. Nansamba, E. Odama, P. Musinguzi, P. Ebanyat, T. Hyuha, A.K. Kasharu & C.S. Wortmann (2012). Optimizing smallholder returns to fertilizer use: Bean, soybean and groundnut. Field Crops Research 127: 109-119.
[12]   | Ronner, E., K. Descheemaeker, C.J.M. Almekinders, P. Ebanyat & K.E. Giller (2017). Farmers’ use and adaptation of improved climbing bean production practices in the highlands of Uganda. Agriculture, Ecosystems and Environment 261: 186-200.
[13]   | FAO (2018). FAOSTAT: FAO Statistical Databases. Available online at: http://faostat.fao.org/
[14]   | Wortmann, C.S., R.A. Kirkby, C.A. Eledu & D.J. Allen (1998). Atlas of common bean (Phaseolus vulgaris L.) production in Africa. No 297. CIAT, Cali, Colombia.

[15]  
Okii, D., P. Tukamuhabwa, G. Tusiime, H. Talwana, T. Odong, C. Mukankusi, A. Male, W. Amongi, S. Sebuliba, P. Paparu, S. Nkalubo, M. Ugen, S. Buah & P. Gepts (2017). Agronomic qualities of genetic pyramids of common bean developed for multiple-disease-resistance. African Crop Science Journal 25(4): 457-472.

[16]   | Ashango, Z., B. Amsalu, K. Tumisa, K. Negash & A. Fikre (2016). Seed Yield Stability and Genotype x Environment Interaction of Common Bean (Phaseolus vulgaris L.) Lines in Ethiopia. International Journal of Plant Breeding and Crop Science 3(2): 135-144.
[17]   | Tadesse, T., A. Tekalign, B. Mulugeta & G. Sefera (2017). Identification of Stability and Adaptability of Small Red Bean Cultivars Using AMMI Analysis. Plant 5(6): 99-103.
[18]   | Corrêa, A.M., A.R.S. Lima, D.C. Braga, G. Ceccon, P.E. Teodoro, A.C. Silva Junior & F.A. Silva (2015). Agronomic Performance and Genetic Variability among Common Bean Genotypes in Savanna/Pantanal Ecotone. Journal of Agronomy 14(3): 175-179.
[19]   | Corrêa, A.M., M.C. Gonçalves & P.E. Teodoro (2016). Pattern analysis of multi-environment trials in common bean genotypes. Bioscience Journal 32(2): 328-336.
[20]   | Wortmann, C.S. & D.J. Allen (1994). African bean production environments: their definition, characteristics and constraints. Network on Bean Research in Africa, Occasional Paper Series No. 11, Dar es Salaam, Tanzania.
[21]   | Wahome, S.W., P.M. Kimani, J.W. Muthomi, R.D. Narla & R. Buruchara (2011). Multiple disease resistance in snap bean genotypes in Kenya. African Crop Science Journal 19(4): 289-302.
[22]   | Njoki, N.W.B. (2013). Breeding for durable resistance to angular leaf spot (Pseudocercospora griseola) in common bean (Phaseolus vulgaris) in Kenya. Ph.D Thesis, University of Kwa Zulu-Natal, Republic of South Africa, p.145.
[23]   | Jaetzold, R., H. Schmidt, B. Hornetz, and C. Shisanya. 2006. Farm Management Handbook of Kenya. Vol II, Natural conditions and farm management information, 2nd Edition Part B Central Kenya. Subpart B2. Central Province.
[24]   | Schoonhoven, A. & M.A. Pastor-Corrales (1987). Standard System for the Evaluation of Bean Germplasm. Centro Internacional de Agricultura Tropical, CIAT Apartado Areo 6713 Cali, Colombia, p.56.
[25]   | VSN International (2014). GenStat reference manual (17th edition). VSN International, Hemel Hempstead, UK.
[26]   | USDA and NRCS (2007). Statistix 8 User Guide for the Plant Materials Program, USA, p.80.
[27]   | Gauch, G.H. & R.W. Zobel (1997). Interpreting mega-environments and targeting genotypes. Journal of Crop Science 37(2):311-326.
[28]   | Gauch, H.G., H.P. Piepho & P. Annicchiarico (2008). Statistical analysis of yield trials by AMMI and GGE: Further considerations. Crop science 48(3): 866-889.
[29]   | Zobel, R.W., M.J. Wright & H.G. Gauch (1988). Statistical analysis of a yield trial. Agronomy Journal 80(3): 388-393.
[30]   | Purchase, J.L. (1997). Parametric analysis to describe genotype x environment interaction and yield stability in winter wheat. PhD. Thesis. University of the Orange Free State.
[31]   | Samonte, S.O.P., L.T. Wilson, A.M. McClung & J.C. Medley (2005). Targeting cultivars onto rice growing environments using AMMI and SREG GGE biplot analyses. Crop Science 45(6): 2414-2424.
[32]   | Assefa, T., I.M. Rao, S.B. Cannon, J. Wu, Z. Gutema, M. Blair, P. Otyama, F. Alemayehu & B. Dagne (2017). Improving adaptation to drought stress in white pea bean (Phaseolus vulgaris L.): Genotypic effects on grain yield, yield components and pod harvest index. Plant Breeding 136(4):548-561.
[33]   | Mwale, V.M., J.M. Bokosi, C.M. Masangano, M.B. Kwapata, V.H. Kabambe & C. Miles (2008). Yield performance of dwarf bean (Phaseolus vulgaris L.) lines under Researcher Designed Farmer Managed (RDFM) system in three bean agro-ecological zones of Malawi. African Journal of Biotechnology 7(16).
[34]   | Beebe, S.E., I.M. Rao, M.W. Blair & J.A. Acosta-Gallegos (2013). Phenotyping common beans for adaptation to drought. Frontiers in Plant Physiology 4:35.
[35]   | Rao, I., S. Beebe, J. Polania, J. Ricaurte, C. Cajiao, R. García & M. Rivera (2013). Can tepary bean be a model for improvement of drought resistance in common bean? African Crop Science Journal 21(4): 265-281.
[36]   | Rao, I.M., S.E. Beebe, J. Polania, M. Grajales, C. Cajiao, J. Ricaurte, R. García & M. Rivera (2017). Evidence for genotypic differences among elite lines of common bean in the ability to remobilize photosynthate to increase yield under drought. Journal of Agricultural Science 155(6): 857-875.
[37]   | White, J.W. & S.P. Singh (1991). Breeding for adaptation to drought. pp. 501-560. In Schoonhoven, A.V., and O. Voysest (Eds.). Common Beans: Research for Crop Improvement. C.A.B. International, CIAT, Cali, Colombia.

[38]  
Blair, M.W., C.H. Galeano, E. Tovar, M.C.M. Torres, A.V. Castrillón, S.E. Beebe & I.M. Rao (2012). Development of a Mesoamerican intra-gene pool genetic map for quantitative trait loci detection in a drought tolerant × susceptible common bean (Phaseolus vulgaris L.) cross. Molecular Breeding 29(1): 71-88.

[39]  
Mwale, V.M., J.M. Bokosi, C.M. Masangano, M.B. Kwapata, V.H. Kabambe & C. Miles (2009). Performance of climber common bean (Phaseolus vulgaris L.) lines under Researcher Designed Farmer Managed (RDFM) system in three bean agro-ecological zones of Malawi. African Journal of Biotechnology 8(11): 2460-2468.

[40]   | Singh, S.P. (1989). Patterns of variation in cultivated common bean (Phaseolus vulgaris, Fabaceae). Economic Botany 43(1): 39-57.
[41]   | Singh, S.P. & C.A. Urrea (1995). Inter- and intra-racial hybridization and selection for seed yield in early generations of common bean, Phaseolus vulgaris L. Euphytica 81(2):131-137.
[42]   | Lima, E.R., A.S. Santiago, A.P. Araújo & M.G. Teixeira (2005). Effects of the size of sown seed on growth and yield of common bean cultivars of different seed sizes. Brazilian Journal of Plant Physiology 17(3): 273-281.
[43]   | Debouck, D.G., O. Toro, O.M. Paredes, W.C. Johnson & P. Gepts (1993). Genetic diversity and ecological distribution of Phaseolus vulgaris in northwestern South Africa. Economic Botany 47(4): 408-423.
[44]   | Darkwa, K., D. Ambachewa, H. Mohammed, A. Asfawa & M.W. Blair (2016). Evaluation of common bean (Phaseolus vulgaris L.) genotypes for drought stress adaptation in Ethiopia. The Crop Journal 4(5): 367-376.
[45]   | Mekbib, F. (2003). Yield stability in common bean (Phaseolus vulgaris L.) genotypes. Euphytica 130(2): 147-153.
[46]   | Lad, D.B., N. Longmei & U.M. Borle (2017). Studies on Genetic Variability, Association of Characters and Path Analysis in French Bean (Phaseolus vulgaris L.). International Journal of Pure and Applied Bioscience 5(6): 1065-1069.
[47]   | Polania, J.A., C. Poschenrieder, S. Beebe & I.M. Rao (2016). Effective Use of Water and Increased Dry Matter Partitioned to Grain Contribute to Yield of Common Bean Improved for Drought Resistance. Frontiers in Plant Science 7: 660.
[48]   | Gereziher, T., E. Seid & G. Bisrat (2017). Performance evaluation of common bean (Phaseolus vulgaris L.) varieties in Raya Valley, Northern Ethiopia. African Journal of Plant Science 11(1):1-5.
[49]   | Tamene, T.T & S.G. Tadese (2014). Sites Regression GGE Biplot Analysis of Haricot Bean (Phaseolus vulgaris L.) Genotypes in three Contrasting Environments. World Journal of Agricultural Research 2(5): 228-236.

[50]  
Tadesse, T., A. Tekalign, B. Mulugeta & G. Sefera (2018). Evaluation of the effect of genotype, environment and genotype x environment interaction on white common bean varieties using additive main effect and multiplicative interaction (AMMI) analysis in the mid altitude of Bale zone, Southeastern Ethiopia. African Journal of Agricultural Research 13(7):338-344.

[51]   | Swegarden, H.R., C.C. Sheaffer & T.E. Michaels (2016). Yield stability of heirloom dry bean (Phaseolus vulgaris L.) cultivars in midwest organic production. HortScience 51(1):8-14.
[52]   | Lin, C.S., M.R. Binns & L.P. Lefkovitch (1986). Stability analysis: where do we stand? Crop Science 26(5): 894-900.