Trait Associations in Diversity Panels of the Two Common Bean (Phaseolus vulgaris L.) Gene Pools Grown under Well-watered and Water-Stress Conditions

Common beans are a warm-season, food legume cultivated in areas prone to water limitation throughout their growing season. This study assessed the magnitude and pattern of trait associations for a total of 202 common bean genotypes divided into panels of 81 Andean and 121 Mesoamerican gene pool accessions grown under contrasting treatments of well-watered, non-stress, and water-limited, terminal drought-stress conditions. Linear correlation, complex path coefficient, and genetic divergence analyses were used to dissect the relationship dynamics between traits and the relative contribution of adaptive traits to differentiation among gene pools and genotypes based on drought stress. Drought severity level for the trial was high and created the ideal condition to reveal genotypic differences, as seen by the differential response of the genotypes for the various traits measured. The value for phenotypic coefficients of variation for all traits was higher than the corresponding genotypic values. Seed yield had positive and strong genotypic and phenotypic correlation with pods per plant across gene pools and stress levels. The overall amount of genetic correlation was greater than the corresponding phenotypic correlation matrix for all the traits within the gene pool and across stress levels. Moreover, the results depicted the phenotypic correlation as equal or better than its genotypic counterpart in estimating drought tolerance in common bean plants. Clustering analysis with Mahanalobis's coefficient of generalized distance grouped genotypes with a differential level of drought adaptation into different classes within each panel. This indicates drought tolerance involves different mechanisms of plant response and is present separately in each gene pool panel. Pods per plant, seed weight, pod partitioning index, and harvest index are useful selection objectives to improve drought adaptation in common bean, but must be differentially weighted in each gene pool. The analysis of genetic variation and association between adaptive traits on the two panels provided useful insights on which traits could be used to improve common bean adaptation to low water availability during the growth season.


INTRODUCTION
Common bean (Phaseolus vulgaris L.) is one of the important grain legumes in food and agricultural systems of Africa and Latin America. It is widely cultivated throughout the tropics for its edible green leaves, green pods, mature or immature seeds for human food and straw as fodder for animal feed (Broughton et al., 2003;De Luque et al., 2014). Despite its important role in agri-food systems of the tropics, common bean production in many regions is being challenged by drought stress that is part of the many calamities of climate change and a frequent abiotic constraint in East Africa (Wortmann et al., 1998;Asfaw and Blair, 2014). Drought stress, which refers to inadequate water availability in quantity and distribution during the life cycle of the crop, is the most important production risk for common bean worldwide . Drought stress potentially affects over 60% of the common bean production area globally in any given year, including vast areas in eastern and southern Africa, Latin America and the Caribbean (White and Singh, 1991;Thung and Rao, 1999). Drought is the most significant contributor to yield reduction and seed insecurity in the rain-fed common bean production systems of most developing countries, where there are few investments in irrigation to raise a good common bean crop and where bean growing is characteristic of poverty hot spots (Asfaw et al., 2013).
Drought during a planting season has a multi-faceted effect on common bean growth and performance whenever and wherever it occurs. This abiotic stress causes poor root growth and development, reduced P uptake and N fixation, poor rates of photosynthesis, low biomass production, and inefficient partitioning (Acosta-Gallegos and Kohashi-Shibata, 1989;Singh, 1995;Serraj and Sinclair, 1998;. Drought also results in flower abortion, pod drop, reduced seed filling, lower grain weight and ultimately in lower seed yield Asfaw and Blair, 2014;Mukeshimana et al., 2014). Drought stress can cause 20-90% common bean yield reduction in farmers' fields which in the worst scenario could be up to 100% yield loss (Ambachew et al., 2015). Fortunately, the drought problem in common bean production can be mitigated.
Development and deployment of climate resilient varieties of common bean that can stabilize or increase seed yield under drought stress while meeting the producer and consumer demands would position this legume to have a greater role in improved livelihoods in the tropics and elsewhere. These crop improvements require a better understanding of the dynamics of the drought tolerance mechanisms functioning during the exposure of common bean to differential water availability and their relative contributions to stabilized or increased seed yield (Ambachew et al., 2015). Improved seed yield under variable stress factors is the result of not only a single trait but rather a cumulative function of many interdependent plant traits . The genetic and phenotypic relationships of traits are of practical interest in breeding programs as selection for one trait may cause improvement or deterioration in associated trait (Baker, 1986) and simultaneous selection of multiple traits would be desirable if possible.
Accurate assessment of the interconnectedness of plant traits under exposure to variable and low water availability and their contribution to seed yield formation is a key factor for making breeding progress in drought stress. Different methods have been extensively employed in many crops including common beans to understand and exploit the genetic basis of relationship between overall and component production traits in the breeding of new varieties. These methods include simple correlation and causation assessment using linear correlations and regression or complex path coefficients and principal component analysis (Wright, 1923;Li, 1956;Sharma, 1995;Falconer and Mackay, 1996;Waitt and Levin, 1998;Yan and Rajcan, 2000;Önder et al., 2013). Despite the potential of these methodologies for dissection of the genetic basis of trait associations and determining the selection strategies in breeding programs, reports in the literature on its application to common bean under drought conditions are scanty. Previous trait association studies on drought affected common bean plants are primarily based on breeding populations or advanced lines and do not provide both genotypic and phentoypic correlation of traits together. Therefore, few studies in common beans report the extent to which phenotypic correlations between traits reflect its genotypic counterpart over stress factors.
While most drought studies in common bean have used limited genetic diversity found in one or two races of the crop, in this study we develop a germplasm panel for the dissection of drought tolerance traits in all types of the grain legume. Our overall goal was to assess the inter-relationships of multiple drought adaptive traits in two major gene pools of common bean using two reference diversity panels for Andean and Mesoamerican beans from two previous studies grown in contrasting water availability: one with high and one with low total rainfall during the growing season in southern Ethiopia. The reference panels were expected to be informative based on the fact that the different gene pools are the result of separate domestication process in common bean (Blair et al., 2009) and therefore should not be expected to have similar genotypic and phenotypic trait correlations. This study compared and contrasted the pattern and structure of trait correlations and differentiation in the two gene pools using correlation, path and genetic variability analyses. Genetic correlation matrices were compared with their phenotypic counterparts to assess how close the correlation estimates were for traits in the two different gene pools and across stress levels. Such analysis could yield significant additional insights on the organization of trait interrelationships in the common bean plants and would inform a breeding program on how to consolidate the main traits together into new varieties for tropical environments like the one used in the study.

The Study Material and Trial Design
The study was conducted at Hawassa, in South Nations, Nationalities and Peoples Regional State (SNNPRS), Ethiopia. Hawassa is located at the 7 • 03 ′ N latitude and 38 • 30 ′ E longitude at an elevation of 1700 meter above sea level. The soil at this site is a sandy loam (Flavisol, FAO classification) with pH 7.0 which is characterized by low water retention capacity and which is quick to drain under normal rainfall and drainage system. Hawassa has a bimodal rainfall pattern with an extended period of wet season from March to October with mean monthly rainfall varying from 18 to 128 mm and annual maximum and minimum temperatures ranging from 24 to 29 • C and 9.5 to 14.1 • C, respectively, based on the long-term weather record from Hawassa station. A quick onset of permanent plant wilting in cases of rainfall deficit during the season due to intense evaporation and percolation is the main feature of the experimental site.
The study material consisted of 202 common bean accessions, of which 81 genotypes were from the Andean gene pool and 121 genotypes were from the Mesoamerican gene pool. Both diversity panels represented accessions held at the International Centre for Tropical Agriculture (CIAT) as described in Pérez et al. (2011) andSimbarashe (2013). The genotypes from the two gene pools were evaluated in separate trials planted in a lattice design with three replications. The trials were planted from July-October 2009 season using two different sowing dates: one that was early to avoid drought and one that was late to expose the crop to terminal drought. These differing conditions created two treatments, one with high and one with low total seasonal rainfall during the growing season. The trials were planted in four rows of 3 m long plots using 60 cm between-row and 10 cm within-row spacing. Across all trials, diammonium phosphate fertilizer was side drilled into each row in a plot at the rate of 100 kg ha −1 at planting. A seed dressing of fungicide (Benomyl, a systemic benzimidazole) and insecticide (Thiamethoxam, a systemic insecticide in the class of neonicotinoids) was applied to reduce the effect of diseases and insect pests prevalent during both the rainy and dry season planting for valid testing of the drought resistance of the germplasm. Therefore, disease and insect pressure was insignificant during both trials and for the season overall. The experimental plots were hand-weeded before flowering and again when needed. Total rainfall was recorded on a daily basis with a rain catchment system. Drought-stressed and non-stressed conditions were considered separate experiments for analysis.

Plant Trait Measurements
Multiple plant traits were measured using non-destructive and destructive sampling at different crop growth stages based on accepted common bean drought trait ontology (http://www. cropontology.org/ontology/CO_335/Common%20bean). Traits were: (1) days to flowering (DF) based on number of days from sowing to when 50% of plants in a plot opened at least one flower; (2) days to physiological maturity (DPM) based on number of days from sowing to when the first pod begins to discolor in 50% of the plants; (3) non-destructive SPAD chlorophyll meter reading (SCMR) measured at mid-pod filling stage, about 1 month after flowering on 10 fully expanded young leaves of five comparable plants in each plot using a portable SPAD-502 chlorophyll meter (Minolta Camera Co., Ltd., Japan). Other measurements were recorded on five comparable plants per plot sampled at mid pod-filling and at harvest using sampling described in  or Beebe et al. (2013) including; (4) pod partitioning index (PPI, %) determined as the ratio of dry weight of pods at harvest over dry weight of total biomass at mid-pod fill multiplied by 100; (5) pod harvest index (PHI, %) determined as the ratio of dry weight of seed over dry weight of pod at harvest multiplied by 100; and (6) harvest index (HI, %) determined as the ratio of seed dry weight at harvest over dry weight of total biomass at the mid-pod fill stage multiplied by 100. The biomass sampled from each plot at mid-pod fill and harvest were hand-separated into leaves, stems, pods and seeds and were oven-dried at 80 • C for 48 h. The dried samples were weighted for determination of the indices. Meanwhile, we also measured four additional traits, namely trait (7) pods per plant (PDPL); trait (8) seeds per pod (SDPD); trait (9) seed weight (100SW, g/100seed); and trait (10) seed yield per hectare (YLDH, kg/ha) all recorded at harvest. PDPL was recorded by counting the average number of pods on five randomly selected plants in a plot. SDPD was recorded by counting the average number of seeds on 10 randomly selected pods from five plants in a plot. Seed weight was measured as weight in grams of 100 seeds and seed yield per hectare as weight (kg) of seed harvest divided by the effective plot in m 2 and multiplied by 10. Seed yield was on a plot basis after the grain was sun-dried to maintain the seed moisture at 12 % and recorded using a sensitive digital balance (A&D FX3000i, 3,200 × 0.01 g). Simple sampling was also used to determine 100SW as the weight of 100 seeds randomly sampled from the plot yield. In addition to the individual genotype measurements, drought intensity index (DII%) calculated as DII = [1 − (Y ds /Y ns )] 100, where Y ds and Y ns are the mean seed yield of all genotypes under drought stress and non-stress treatments, respectively (Fischer and Maurer, 1978), was used to quantify the severity of drought stress on seed yielding potential of the genotypes. Drought intensity index was also calculated for other traits measured in the trials to quantify the drought stress effect on trait expression compared with their non-stress counterpart. We also estimated the drought response index (DRI) for individual genotypes as suggested by Bidinger et al. (1987) by adjusting genotype capacity for seed yield potential and growth duration. For the drought response estimation, we regressed stress seed yield on potential (non-stress) yield and phenology (days to flowering and days to physiological maturity under nonstress) and calculated standardized residual from the regression estimation as the difference between measured seed yield under drought stress and predicted stress seed yield divided by the standard error. The standardized residual value was used to define the genotype as drought resistance (a positive deviation from the regression line) or susceptible (a negative deviation), independently of the effect of potential yield or growth duration on its seed yield under stress. The response of a genotype to drought was assumed zero (no response) if the difference between predicted and measured seed yield does not exceeds the experimental error.

Statistical Analysis
The data recorded for each trait per trial were initially subjected to mixed model analysis using the program Genstat v. 12.1 (Payne et al., 2007) to test the genotype differences in respective trials. We fit the data to a linear mixed model following Gilmour et al. (1997) using Genstat, where replication and genotypes were fixed and block was a random factor. We also used flowering and maturity time as a covariable in the model to correct their effect on seed yield. The restricted maximum likelihood (REML) procedure was used to estimate the variance components and residual variances. The model used was: where Y ilm was plot mean performance of a certain genotype i, in replication l and block m, µ the overall mean, g i the effect of genotype i, r l the effect of replicate l, b m the effect of block m and ε ilm the residual associated with plot for single trial analysis. The variance component estimates from the mixed model analysis were used calculate the phenotypic (PCV) and genotypic (GCV) coefficients of variation as where, σ p , σ g , and X were the phenotypic, genotypic standard deviation and grand mean of the traits, respectively (Singh and Chaudhary, 1985). Heritability in the broad sense (H 2 ) was estimated on genotypic mean as: where σ 2 g and σ 2 p were genotypic and phenotypic variance, respectively (Allard, 1999;Galwey, 2014). Expected genetic advance (GA) and percentage of GA were calculated according to Shukla et al. (2006) GA = iσpH 2 (5) where i was standardized selection differential, a constant (2.08) and σ p was the standard deviation of the phenotype. A linear correlation analysis was applied in a pair-wise manner between traits to all the characteristics measured. SAS Proc MIXED model was employed for REML estimation of genetic and phenotypic correlations among traits and their standard errors following Holland (2006). The correlation matrixes between traits were visualized graphically with RCricos package (Krzywinski et al., 2009). Pairs of genetic and phenotypic correlation matrices from the REML estimation of each trial were used to calculate the overall magnitude of correlation within a matrix, average disparity and pattern of similarity between corresponding genetic and phenotypic correlation matrices. The overall magnitude of correlation within a matrix was calculated as: where r i,j refers to the correlation between traits i and j, and n was the number of off-diagonal non-redundant elements in the matrix (Waitt and Levin, 1998). The overall magnitude of correlation within a matrix was implemented to assess the extent to which phenotypic correlation reflected their genetic counterparts. The average disparity between corresponding genotypic and phenotypic correlation matrices was determined as: where r G,i,j and r P,i,j referred to genetic and phenotypic correlations between trait i and j, and n was the number of offdiagonal non-redundant elements in the matrix (Willis et al., 1991). The average disparity between corresponding genetic and phenotypic correlation matrices was implemented to assess the proximity, on average, of the estimates of genotypic and phenotypic correlation. A pattern similarity of corresponding genetic and phenotypic correlation matrices was assessed using a matrix correlation based on: (9) where a ij and b ij were the i th and j th elements of correlation matrices being compared (A and B), a ij and b ij were mean of the off-diagonal elements in matrices A and B, respectively, and n was the number of traits sampled (cf. Waitt and Levin, 1998). The large positive matrix correlation indicated that genetic and phenotypic correlation vary in similar directions but not that the magnitudes of individual correlations were identical. Hence, the matrix correlation and average disparity were considered together for interpretation of the results. High matrix correlations and low average disparity indicated that the genetic and phenotypic correlation were similar in magnitude and tended to occupy the same position in their respective matrices (Waitt and Levin, 1998). The statistical significance of matrix correlation was tested using Mantel's randomization test (Mantel, 1967) in Genstat. In addition to assessing the magnitude, disparity and pattern similarity of corresponding genetic and phenotypic correlation matrixes, the causal relationship involved for the seed yield and all other component traits measured for the trials were assessed with path coefficient analysis as per the procedure described in Wright (1923Wright ( , 1934. Level of genetic divergence of the diversity panels on DS attributes and their relative contribution to the differentiation of genotypes was assessed with Mahalanobis D 2 -statistics as described in Sharma (1998). The correlations analyses were done using SAS v 9.1software. Clustering and cluster statistics analyses were performed in R statistical software using various packages (eclodist, stats and clv). Mahalanobis's distance matrix in Ward2 hierarchical agglomerative clustering method (Murtagh and Legendre, 2014) was implemented to group the genotypes. The path coefficient analyses were performed using Microsoft Excel.

Variation on Seed Yield and Other Traits
The weather parameters (maximum and minimum daily temperatures and the daily total rainfall) during the cropgrowing period are presented in Figure 1. The non-stress trial received a total rainfall of 235 mm during the growth period through evenly spaced rains ( Figure 1A). Meanwhile the drought stress trial was planted later and only received 149 mm, with terminal drought (Figure 1B). The rainfall variation during the crop's life cycle caused a seed yield penalty under drought stress compared to the non-stress treatment and poor performance of plant traits in the terminal drought ( Table 1).
Drought stress resulted in a reduction in expression of all common bean traits assessed in this study except pods per plant in the Andean accessions in which a slight increase was observed with drought exposure. The overall seed yield reduction (drought severity index) was 59 and 63% under drought stress for the Andean and Mesoamerican gene pool accessions, respectively. Drought also caused early flowering and maturity in the bean diversity panels. Stress trials showed 9% (in Andean accessions) to 17% (in Mesoamerican accessions) reduction in leaf chlorophyll content. Exposure to drought also resulted in 37 and 13 % reduction in seeds per pod in Andean and Mesoamerican accessions, respectively. The stress on 100 seed weight and pod harvest index were generally mild irrespective of the gene pool origins as drought caused less than 30% reduction in these traits' mean performance in the accessions. Drought stress effects were lower for the Andean accessions in pod partitioning and harvest index compared to the Mesoamerican accessions as percent reduction in mean trait values were less than 30% Andean and above 40% in the Mesoamerican gene pool accessions.
In general, Mesoamerican diversity panel accessions showed slightly higher sensitivity to drought stress for expression of seed yield, days to flowering and maturity, pod partitioning index and harvest index compared to the Andean accessions for this environment. Andean accessions expressed more sensitivity to drought for expression of seed weight compared to the Mesoamerican accessions. Table 1 also presents different statistics for comparing the variation in various common bean plant traits between gene pools and across stress and nonstress conditions. Highly significant genotypic differences were observed for the expression of all the different traits assessed within the Andean and Mesoamerican gene pool accessions and across the drought stress and non-stress treatments.
The values of phenotypic coefficient variability (PCV) were higher than the corresponding genotypic coefficient of variability (GCV) values for all the traits within gene pools and across the water regimes. PCV values ranged from 3.90 to 60.63 % for Andean accessions and from 3.91 to 43.38 % for the Mesoamerican accessions under drought stress with the highest values for the trait of seed yield per hectare. Under nonstress treatment, the PCV values ranged from 3.52 (days to physiological to maturity) to 36.3 3% (seeds per pod) for the Andean gene pool accessions and 3.16 (days to physiological maturity) to 28.85% (pods per plant) for the Mesoamerican gene pool accessions.
Broad sense heritability estimates were in general greater than 0.50 for the majority of the traits except SCMR in the Andeans or DF and DPM in the Mesoamericans under drought stress, and PDPL in the Andeans or SDPD in the Mesoamericans under non-stress treatment. Genetic advance as percent of the mean was highest for yield (46.75-122.31%), 100 seed weight (39.82-48.00%) and pod partitioning index (45.86-71.16%) while the remaining traits showed a moderate to very low amounts of genetic advance over gene pools and across stress levels. The genetic advance values were observed to be higher for seed yield, pods per plant, pod partitioning index and harvest index on both gene pool accessions under drought stress compared to the nonstress treatment. This indicates that selection will be beneficial for improvement of these traits. Figure 2 presents drought response of the genotypes in the trials based on the simultaneous regression of measured stress FIGURE 1 | Rainfall distribution (bars), as well as maximum and minimum temperatures (red and blue broken lines, respectively) during the crop growth period of two separately planted experiments (A) non-stress with early sowing to avoid drought and (B) drought-stress due to late sowing to expose the crop to low overall rainfall and terminal drought. The x-axis shows days after planting at the bottom of the subfigures and the y-axes to the right and left of the graphs show the daily rainfall and temperature records, respectively. TABLE 1 | Probability values for testing genotypic effect, mean traits values and REML estimates of variance components (σ 2 g, σ 2 e, σ 2 p), broad sense heritability (H 2 ), genotypic (GCV) and phenotypic (PCV) coefficients of variation, genetic advance (GA%) and drought intensity index (DII%) for various adaptive and seed yield traits in Andean and Mesoamerican common bean gene pool diversity panels grown under contrasting water availability. seed yield on the non-stress potential seed yield and growth duration (flowering and maturity time). The magnitude of deviation of a genotype from the regression line which is referred as a drought response index (DRI) was used to assess the response of genotypes to drought which is not explained by the yield potential and/or growth duration. We regressed measured stress seed yield on yield potential, growth duration (flowering and physiological maturity time) and DRI to assess the individual contribution of these four factors to explaining the variation in seed yield under stress. Yield potential and times to flowering were major factors contributing to variation in the seed yield in both gene pools. Yield potential explained an average of 23 and 15 % of the variation in seed yield under terminal drought in Andean and Mesoamerican accessions, respectively (data not presented). Flowering time explained an average of 11% of the variation in seed yield in both gene pools while times to physiological maturity was low in explaining the variation in seed yield under stress. Deviation from the regression or DRI explained 65% yield variation in Andean and 71% of the variation in seed yield under stress in Mesoamerican genepool accessions. More than half of the genotypes in the trial (both Andean and Mesoamerican genepool) expressed no specific response to the drought stress indicating their measured seed yield under stress was adequately estimated by their potential yield and growth duration (flowering time and maturity). These genotypes are those with a drought response index value of near zero or within the limits of experimental error hence are average in their susceptibility to drought. Few genotypes expressed a high deviation from the regression line or high DRI values. Genotypes such as AND1005, SEQ1027, PVA111, G6873, G4493, G1836, and AFR619 from the Andean genepool and G801, G10945, G10982, G18454, G2093, G2277, G2379, G2635, G5733, G7932, and PINTOVILLA from the Mesoamerican gene pool with a high and positive deviation from the regression line showed a resistance reaction to the drought stress treatment. Whereas, G2686, G17070, G12529, G11787, and DRK47 from the Andean accessions and DOR364, MASSAYRED, G16400, G1815, G19941, G22787, G2445, G3595, G7742, and G7863 from the Mesoamerican accessions with a negative and higher deviations from the regression line expressed a susceptible response to drought. Application of the drought response index in this study identified drought resistant and susceptible genotypes of common beans. Such an attempt to dissect the genotype's drought response status should consider also the absolute yield under non-stress and drought stress along with other traits defining drought adaptation and cultivar acceptability in breeding of common beans for tropical environments.

Correlations among Adaptive Traits and between the Seed Yield
Correlations between seed yield and adaptive traits determined for accessions of the two common bean gene pools evaluated under both drought stress and non-stress treatments are presented in Figure 3 and Table 2. Seed yield had positive and high genotypic and phenotypic correlation with PDPL, 100SW, PPI, and HI under both drought stress and non-stress treatments in Andean gene pool accessions (Figure 3, Supplementary Table  1). While for the Mesoamerican gene pool accessions, the correlation with seed yield was positive and significant only for PDPL both at genotypic and phenotypic levels across the water regimes (Table 2). Similarly, days to flowering had a significant and positive genotypic and phenotypic correlation with days to physiological maturity over the gene pools and across the water regimes. SPAD chlorophyll meter reading had significant but negative genotypic and phenotypic correlation with days to flowering and days to physiological maturity over the gene pool and across water regimes. The association of SAPD chlorophyll reading with SDPD was significant but negative both at genotypic and phenotypic correlation under drought stress while it was significant and negative at phenotypic level and nonsignificant genotypic level under non-stress treatment for the Andean accessions. Correlations among yield parameters were more complex than for the phenological and photosynthetic traits. For example, the correlation between PDPL and SDPD was negligible both at genotypic and phenotypic levels across the drought stress and non-stress treatments for the Mesoamerican accessions, but was significant at 0.05 p value for the Andean accessions. Hundred seed weight had negative and highly Under the non-stress treatment, 100 seed weight showed negative and significant correlation with days to flowering both at genotypic and phenotypic level for the Mesoamerican accessions while it was negative and high at phenotypic level and weak at the genotypic level for the Andean accessions. Pod partitioning index showed strong positive correlation (both at genetic and phenotypic level) with harvest index over the gene pools and across stress factors. The positive and significant associations were observed between the two indices HI and PHI, both at genotypic and phenotypic levels under the drought stress on both gene pool accessions while the same correlation was positive and significant at genotypic and phenotypic levels for the Andean accessions and positive but non-significant for the Mesoamerican accessions under non-stress. The correlations of pods per plant with seed yield, harvest index with pod harvest index, and days to flowering with days to physiological maturity suggest the potential for simultaneous improvement of these traits through breeding. Furthermore, selections for increased levels of traits such as pods per plant, seed weight, pod partitioning index and harvest index in Andean and pods per plant in Mesoamerican gene pool accessions would improve the seed yielding potential under both drought stress and non-stress environments.

Magnitude, Disparity and Pattern Similarity of Genetic and Phenotypic Correlations
The corresponding pairs of genetic and phenotypic correlations of traits in the diversity panels representing the Andean and Mesoamerican gene pools of common bean grown under contrasting water regimes (well-watered non-stress and water stressed terminal drought treatments) were compared with average absolute values of genetic and phenotypic correlation matrices, average disparity and matrix correlation ( Table 3). The overall magnitude of correlations as measured by absolute value had a mean value of 0.313 for the genotypic correlation compared with 0.206 for the corresponding phenotypic correlation. The overall magnitude of genetic correlation matrices was higher than the corresponding phenotypic correlation in 100% of the comparisons ( Table 3). The average disparity as a measure of how close, on average, the corresponding correlation matrices was ranged from a low of 0.080 to a high of 0.282. The correlations of seed yield with other traits in Andean accessions under non-stress had little disparity whereas genetic correlation among all traits in Andean accessions evaluated under drought stress with its corresponding phenotypic correlation matrices under non-stress showed a high disparity.
The matrix correlation that measured the pattern similarity of the magnitude of the corresponding correlation matrices ranged from 0.524 to 0.982 with a mean value of 0.804 ( Table 3). The phenotypic correlation in the Mesoamerican accessions under drought stress with the corresponding phenotypic correlation under non-stress showed little pattern similarity while the corresponding genotypic and phenotypic correlation between seed yield and other traits under drought stress in Andean accessions showed high pattern similarity. Eighty-three percent of the corresponding matrix correlations were more similar than would be anticipated by chance alone with the table-wide significance level of 0.01.

Causal Relationship of Seed Yield with Other Traits
Path coefficients between individual traits and YLDH were calculated separately for the Andean and Mesoamerican gene pool accessions and were useful for dissection of the interrelationships between the adaptive traits and their effects on seed yield at the genotypic level based on the correlation matrices (with direct effects on the main diagonal and indirect effects on both off-diagonal parts of the table). The genotypic correlation matrices of adaptive traits with the seed yield in each gene pool were asymmetrical, consisting of nine components. Meanwhile the direct genotypic effect and the indirect genotypic effects of adaptive traits on seed yield had relationships with the other eight traits.
For the Andean gene pool accessions (Table 4), the direct genotypic effects of pods per plant and seeds per pod on seed yield were positive, and those of SPAD leaf chlorophyll meter reading and pod partitioning index were negative over the stress factors. Days to flowering, 100 seed weight and harvest index had positive and negative direct genotypic effect on seed yield under drought stress and non-stress, respectively. Days to physiological maturity and pod partitioning index had negative direct genotypic effects on seed yield under drought stress and positive direct genotypic effects under non-stress treatment. Days to physiological maturity and pod harvest index had direct 3 | Comparison of the corresponding pairs of genetic and phenotypic correlations common bean traits grown under contrasting water regimes, where /RG/ and /Rp/ are the average absolute values of genetic and phenotypic correlation, respectively; D is the average disparity; rM is the matrix correlation between genetic and phenotypic correlation matrices; PrM is the proportion of permutation matrix correlation exceeding the observed matrix correction.    (Table 5), days to physiological maturity and 100 seed weight had a positive direct effect, and pod harvest index had a negative direct effect on seed yield over stress factor. Days to flowering, SPAD leaf chlorophyll meter reading, pod partitioning index and harvest index had positive direct genetic effects on seed yield under drought stress and negative under non-stress. Pods per plant, seeds per pod and harvest index had negative direct genetic effects on seed yield under drought stress and highest positive direct genetic effect under non-stress. Harvest index had positive indirect under non-stress. Under drought stress, pod-partitioning index had the highest positive direct genetic effect on seed yield.
Meanwhile for this same gene pool, harvest index had genetic effects on seed yield through SPAD leaf chlorophyll meter reading, pods per plant, seeds per pod, pod partitioning index, and pod harvest index under non-stress whereas the same route had negative indirect genetic effect on seed yield under drought stress. Harvest index had the larger positive indirect genetic effect on seed yield through pod partitioning index under drought stress in the Mesoamerican gene pool accessions.
The path analysis results in Tables 4, 5 depicted different trait responses involve in modulating the seed yielding potential of common bean plants in the two gene pools. Pods per plant and seeds per pod in Andean accessions and seed weight (grain filling) in Mesoamerican accessions are potential traits to practice indirect selection improving seed yield both under drought stressed and non-stressed environments. In addition, harvest index, seed weight (better grain filling) and flowering time in Andean accessions and pod partitioning index, leaf chlorophyll content, and flowering time in Mesoamerican accessions are relevant traits to attempt indirect selection improving the seed yielding potential under drought stress.

Differentiation on the Basis of Phenotypic Performance
Clustering analysis of the phenotypic performance under drought stress of the total number of 202 accessions representing the two common bean gene pools using the Mahalanobis's distance in Ward's hierarchical agglomerative clustering method grouped the genotypes in six clusters (Figure 4). The clustering analysis with phenotypic traits did not generate a meaningful differentiation of the accessions into Andean and Mesoamerican gene pools as well as the distinctness in degree of drought adaptation. Accessions from the two gene pools and with varied of level of drought tolerance tended to cluster together. This indicates different mechanisms of plant responses to drought adaptation involve in common bean. The Andean accessions with a high degree of drought tolerance such as CAL143, G19882, AND1005, and SEQ1027 were grouped under clusters I and II while those from the Mesoamerican gene pools such as G10982, G7932, G10945, and PINTOVILLA were grouped under cluster I. Accessions G10982, G10945, and G2402 which belongs to the race Durango of highland Mexico were grouped with the drought bred lines SER16, SEA15, SXB418, and SER109 in cluster I. The genotypes in cluster I will be good sources of genes for drought  adaptation in future breeding. Cluster size varied between the groups identified: with a larger number of genotypes in cluster I (33%) and cluster VI (29% of genotypes); whereas clusters V and II were small groups containing 11 (5%) and 19 genotypes (9%), respectively. Clusters III and V together had 24% of all the genotypes. The Mahalanobis generalized distance for the 202 genotypes ranged from 0.46 to 97.52. The highest inter-cluster divergence appeared between clusters II and V (D 2 = 30.89) followed by clusters III and V (D 2 = 39.10) ( Table 6). The lowest intercluster distance was found between clusters I and VI (D 2 = 17.74) followed by that between clusters III and VI (D 2 = 18.78) and IV and VI. The relative contribution of the multiple traits to the differentiation of the diversity panels at genotypic and the inter-cluster level is demonstrated by the coefficient of variation (CV) values (Figure 5). Seed yield, harvest index, pod partitioning index and 100 seed weight are potential contributors to phenotypic differentiation both at the genotype and inter-cluster level in the common bean gene pool accessions considered for this study. Days to flowering, days to maturity and pod partitioning index were least contributors to phenotypic divergence in the reference collections. Seeds per pod and leaf chlorophyll content were intermediate in contribution to the differentiation of accessions. This indicates no trait is equally important to phenotypic differentiation and drought adaptation in common beans. Methods to recognize relevant traits to drought adaptation in common beans should consider the relevance of different traits in different environmental situations.

DISCUSSION
The drought stress imposed by deficient water availability during the post-rainy season had a significant impact on the seed yield potential, overall performance and plant attributes of the genotypes. The total rainfall received during the drought stress treatment was 37% lower than that of the non-stress treatment. The drought experiment received a low amount of rainfall during flowering that coupled with a total cut-off of rain during pod filling exposed the common bean plant to water limitation over a long critical growth period for the plants (Figure 1). In contrast, the amount and distribution of rainfall over the growth period for the non-stress, rainy-season trial, was near to the average rainfall Values in parenthesis are number of accessions forming cluster (group of accessions resembling each other for their phenotypic performance under drought stress for the ten traits and hence with low intra-cluster divergence). Bold values refer to intra-cluster distance. The χ2 significance was calculated using CHIINV functions in Microsoft excel with n-2 degree of freedom, where n is the number of traits considered for clustering which in this case were 10 traits.
FIGURE 5 | The relative contribution of the multiple traits to the phenotypic differentiation in common bean diversity panels evaluated under drought stress at a genotypic and inter-cluster level as demonstrated by coefficient of variation (CV%) values. Y-axis is CV value for the trait and x-axis is multiple traits contributing to the group formation. Refer Table 1 for trait abbreviations. required by the common bean crop for the optimum growth and development (350-500 mm). In summary, the pattern of rainfall in the Ethiopian highlands allows for drought tolerance testing using staggered alternate sowing dates as shown by our results.
The degree of stress on seed yield potential imposed by alternating the sowing dates as expressed by drought intensity index values were 59 and 63% for the Andean and Mesoamerican gene pool panel accessions, respectively. This drought level was high and adequate to reveal genotypic differences, as seen by the differential response of the genotypes for the various traits measured (Table 1). Previous research on common bean also successfully created consistent stress between experiments by alternating sowing dates to expose the crop to high and low rainfall during the growing period (Acosta-Gallego and White, 1995;Abebe et al., 1998;Rezene et al., 2013;Asfaw and Blair, 2014). Due to the environment used, the kind of abiotic stress that the common bean plants experienced under this study was terminal drought stress, which is a serious problem for common bean production in the tropics Beebe et al., 2014). This study reflected the typical stress condition that common bean growers often face in Ethiopia and also was in accordance with other drought studies in tropical regions where the crop is often grown during the short rains, when precipitation totals decline early Ambachew et al., 2015;Darkwa et al., 2016).
Persistent or intermittent water limitation at any time during the growing season poses a great challenge for common bean production worldwide (Terán and Singh, 2002). This problem is further exacerbated in the rain-fed common bean production system in the tropics where the application of irrigation water is not very feasible, due to deep-rooted poverty and lack of infrastructure. Use of climate resilient varieties is among series of action steps to mitigate the drought problem and position the common bean crop to have a greater role for the livelihood improvement in drought prone production systems of the tropics.
Our study assessed the dynamics of trait linkage and the variability in two major common bean gene pool reference diversity panels grown under contrasting water availability, which provided additional insights instrumental for the making of new resilient varieties. The high variability for and strong positive associations among some adaptive traits over the gene pools and across stress levels were critical aspects of our study which could be useful to common bean breeding programs in best parent selection and trait targets necessary for progress in drought adaptation. The magnitude of drought resistance index for seed yield highlighted the response of some accessions over others (Figure 2) but our results also indicated that drought tolerance was not a function of a single plant trait but rather a cumulative function of many interdependent traits, as observed before (Rao, 2001;Beebe et al., 2013). Traits like pods per plant, pod partitioning index and harvest index showed high heritability combined with good genetic advance across both gene pools and under both treatments hence are important for crop improvement under drought stress. Understanding the pattern of association of such traits at genotypic and phenotypic level is of practical interest in a breeding program as selection for one trait may cause either improvement or deterioration in an associated trait (Baker, 1986;Rao et al., 2013;Ambachew et al., 2015).
In addition to increasing seeds and pods per plant for higher yield, seed weight (grain filling), pod partitioning index and harvest index can be used to improve seed yield for both drought stress and non-stress conditions in Andean gene pool accessions. However, for the Mesoamerican gene pool accessions, increases in pods per plant were associated with seed yield only under nonstress conditions while many pods were unfilled in drought stress. Hence, pods per plant, seed weight, pod partitioning index and harvest index are useful selection objectives to improve drought adaptation in common bean, but must be differentially weighted in each gene pool.
The extent to which the phenotypic correlations between traits reflect its genotypic counterpart is another aspect a breeding program has to take into account while designing a selection strategy for the crop. In this study, the correlation analyses for multiple traits indicated the higher values for genetic correlation matrices compared to the corresponding phenotypic matrices for all the traits, gene pools and stress levels. A review on genotypic and phenotypic correlations in 27 different species revealed that phenotypic correlation is a good reflection of genotypic correlations in plants (Waitt and Levin, 1998). As previously reported by Asfaw and Blair (2014) and Ambachew et al. (2015), the greater contribution of genotypic factors in the development of trait linkages and the strong positive correction between phenotypic traits may be an indication of pleiotropism, genetic coupling and / or linkage disequilibrium with population structure effects. In any case these genetic phenomena provide the opportunity to select superior genotypes that indeed offer increased seed yield and other desirable traits simultaneously (Dilday et al., 1990;Sharma, 1998).
Drought tolerance in common bean is made up of multiple sub-traits and is always quantitative in the nature of its inheritance. Considering trait association analysis results along with genetic diversity assessment could assist bean breeders to pyramid the different plant traits linked with drought tolerance. This study dissected level of difference among study materials for the plant attributes on drought stress treatment. The genetic differentiation assessment with Mahanalobis D 2 statistics indicated a high degree of phenotypic divergence among common bean accessions, between genepools and therefore within and between the reference Andean and Mesoamerican diversity panels. Genotype classification with phenotypic traits did not reveal clear distinction of the least and most highly drought-adapted genotypes belonging to different groups or genepools. This may be explained by the complementarity of drought tolerance traits found crossing the genepool divisions of Andean and Mesoamerican common beans. Accessions from different gene pools clustered together irrespective of their degree of phenotypic resemblance for drought response traits measured in this study and their genepool status, allowing us to conclude that drought tolerance involves different mechanisms of plant response and is present separately in each gene pool panel. Therefore, different phenotyping methods needed to identify relevant genotypes and traits. The genotypes presented here could be an excellent resource for parental selection and the two gene pools could provide complementary traits to improve the drought adaptation in common bean.

AUTHOR CONTRIBUTIONS
AA and MB conceived the study, AA conducted experiments, AA, DA, and TS analyzed data. AA and MB wrote the manuscript, DA contributed to writing. All authors read and approved the manuscript.