Multi-Trait Selection Index for Superior Agronomic and Tuber Quality Traits in Bush Yam ( Dioscorea praehensilis Benth.)

: Dioscorea praehensilis Benth. is a semi-wild yam species and a valuable source of resistance trait genes. To access the agronomic and tuber quality performance, eleven quantitative phenotypic traits were used to discriminate and identify promising accessions among 162 accessions of D. praehensilis collected in Ghana. Signiﬁcant and high genetic variability ( p < 0.001) for all eleven quantitative traits was found among the evaluated accessions. Moderate broad-sense heritability (H 2 ) (30–60%) was observed for all the evaluated quantitative traits except the response to YMV and tuber hardness. The accessions were clustered into three groups; each cluster displayed genotypes with good potentiality for the different traits evaluated. Path coefﬁcient analysis revealed positive contributions ( p < 0.01) of the number of tubers per plant, tuber length, tuber width, stem internode length, number of internodes, and tuber ﬂesh hardness to the total tuber weight per plant. Through the multi-trait genotype–ideotype distance index (MGIDI), 24 accessions were identiﬁed from the 162 evaluated accessions as top-ranking and could be used as progenitors for trait introgression. The results of this study provide insight for future yam breeding and improvement programs in West Africa.


Introduction
Yam (Dioscorea species) is a popular tropical and subtropical food crop.It is widely grown in West Africa, contributing significantly to food security and poverty reduction [1,2].West Africa accounts for more than 95% of global yam production, with Nigeria, Ghana, Côte d'Ivoire, and Benin as the leading producers [3].
D. praehensilis is an edible semi-cultivated wild yam species used chiefly to ease food insecurity among local farmers in forest zones of West African nations such as Nigeria, Ghana, Benin, and Togo during lean seasons [4,5].D. praehensilis has a high yield potential, pest and disease tolerance, in-soil storage ability, and the capacity to blossom and fruit profusely, making it ideal for hybridization in breeding programs [5].
Despite these significances, the economic values of D. praehensilis have not been fully realized due to the oxidative browning and hardening of tuber flesh a few days after harvesting, resulting in poor utilization and under-exploitation of its potential [5].The germplasm collection and the estimation of morphological divergence in D. praehensilis, compared to widely cultivated and utilized yam species such as D. rotundata and D. alata, are partial and not comprehensive, especially in Ghana.These factors have resulted in rapid genetic erosion and the risk of extinction of this valuable yam species [5].The information Agronomy 2023, 13, 682 2 of 16 on tuber culinary quality traits, yield performance in a natural environment, and disease resistance of D. praehensilis is highly deficient and not fully documented.Identifying accessions with high-yielding attributes, tolerance to yam mosaic virus (YMV), and good post-harvest tuber qualities that farmers and consumers prefer will contribute to unlocking the genetic potential of bush yam.
Several studies have been carried out on yam using morphological descriptors [6][7][8][9].Many authors [10][11][12][13] have reported the effectiveness of morphological markers in assessing genetic diversity in D. alata.In Guinea yams (D. rotundata and D. cayenensis), morphological traits have been successfully employed in evaluating genetic variability [14,15].The application of morphological descriptors has also been used in D. bulbifera [16].In addition, the genetic diversity among 140 accessions of D. trifida was successfully determined in the municipality of Caapiranga, in the central Amazon region of Brazil, using 64 morphological descriptors [17].The genetic diversity in D. praehensilis germplasm in Ghana has been less explored, and there is also uncertainty about the applicability of morphological traits to evaluate the genetic diversity in the crop.Hence, a comprehensive analysis of the phenotypic diversity of D. praehensilis germplasm indigenous to Ghana may be critical for identifying and developing cultivars with economically valuable traits and for the conservation and sustainable utilization of the species.
In any crop improvement initiative, breeders often keep in mind a combination of attributes that, when combined into a genotype, would result in excellent performance; this genotype is referred to as an ideotype [18].The goal of ideotype design is to improve crop performance by selecting genotypes based on many attributes at the same time [19].The Smith-Hazel (SH) index is a linear selection index frequently used by breeders for multi-trait selection [20,21].However, in the case of the SH index, the presence of multicollinearity and the difficulties in assigning economic weightage to the qualities under evaluation can have an impact on genetic gain [19].To address these shortcomings, a multivariate selection index, the multi-trait genotype-ideotype distance index (MGIDI), has been created [19].This index accounts for multi-collinearity and favorably selects all variables under consideration, resulting in significant genetic gain [19].The application of MGIDI as a multi-trait selection index has been reported for white Guinea yam (D. rotundata) [22], but limited information is available on D. praehensilis.
The objectives of the present study were: (i) to quantify the agronomic and tuber quality trait performance of D. praehensilis accessions in Ghana, and (ii) to identify D. praehensilis accessions with higher agronomic and yield-related traits for future genetic improvement initiatives.

Experimental Site
The study was conducted at the Teaching and Research Farm of the School of Agriculture, the University of Cape Coast, Ghana (5 • 07 7.6" N, 1 • 17 18.9"W; 15 m above sea level).This farm is located in the Central region with semi-deciduous forests and coastal savannah ecozones.The trial was conducted under field conditions during the 2020 and 2021 growing seasons.The annual rainfalls for the experiment period were 1246.2mm for 2020 and 1170.2 mm for 2021; the average maximum and minimum temperatures for 2020 were 27.9 and 26.9 • C, while those for 2021 were 28.6 and 25 • C, respectively.The average relative humidity for 2020 was 75.7%, and it was 81.2% for 2021.The soil on this experimental site is sandy loam with a pH of 6.72, 1.31% organic carbon, 754.6 µg/g of available phosphorus, and potassium content estimated at 0.081 cmol/kg.

Plant Materials and Experimental Design
A total of 162 D. praehensilis accessions were used, of which 71 were collected from the Central region, 25 from the Eastern region, and 66 from the Western North region of Ghana (Figure S1 and Table S1).These accessions were collected from farmers during the 2019 harvest season.The details, including the accession codes and regions of the collection, are presented in Table S1.The experiment was laid out in a 15-by-11 simple lattice design with two replicates.Each plot size consisted of 3 m long ridges containing three plants at 1 m intra-and inter-row spacing.The field layout was generated using Agricolae package [23].The recommended cultural practices, such as ridging, weeding, staking, etc., were implemented during the growing seasons.

Data Collection
Data were collected for 11 quantitative traits (Table S2) according to the standard operating protocol for yam performance evaluation trials [24] and the yam trait ontology available at YamBase (www.yambase.org,accessed on 7 January 2023).The area under the disease progression curve (AUDPC) for YMV severity, dry matter content, tuber flesh oxidation intensity, and tuber flesh hardness were evaluated as described below.
The AUDPC, a valuable quantitative summary of disease intensity or severity for YMV over time, was estimated using the trapezoidal method [25].This method discretizes the time variable and calculates the average disease intensity or severity between each pair of adjacent time points: where N is the number of observations, y i is the disease severity at the ith observation, and t i is the time at the ith observation.The dry matter content was determined by chopping 100 g of fresh tuber flesh into small pieces and then oven-drying it at 105 • C for 24 h until a constant weight was achieved.The percentage of dry matter content was then estimated as follows: % dry matter content (DMC) = Dry tuber flesh weight (g) Wet tuber flesh weight (g) × 100 (2) The intensity of tuber flesh oxidation (color change or browning of cut tuber flesh) was assessed 60 min after cutting using a Chroma (colorimeter) meter (CR-400, Konica Minolta, Japan).The lightness (L*), red/green coordinate (a*), and yellow/blue coordinate (b*) values were recorded.A reference of white and black porcelain tiles was used to calibrate the Chroma meter before each reading.The delta (color difference) (∆E*) among all three coordinates was calculated using the following formula: where F∆E* is the final delta and I∆E* is the initial delta.Tuber flesh hardness was assessed with a 6.00 mm probe digital penetrometer.Tuber samples of 1 cm thickness and ~5 cm diameter were prepared from each genotype/accession, and the probe was pressed into the tuber.The force necessary for its penetration into the tuber was considered an indicator of the hardness of the tuber.Three measurements were taken per accession, the average was calculated, and the data were expressed in Newtons.

Data Analysis
Lme4 [23], an R package, was used to perform the analysis of variance (ANOVA) using a mixed linear model (MLM) fitted across cropping seasons, as shown below.
where Y ijk is the value of the observed quantitative trait; µ is the population mean; G h is the effect of the hth accession; S i is the effect of the ith growing season; (G h * S i ) is the accessions and season interaction associated with accession h and season i; R ij is the effect of the jth replicate (superblock) in the season ith; B k is the effect of the kth incomplete block within the jth replicate; and ε hijk is the experimental error.Accessions were considered fixed effects, while growing seasons, replicates, and blocks were considered random effects.The means between growing seasons were compared using the least significant difference (LSD) test at a p-value threshold of 0.05.Variations in the quantitative traits among D. praehensilis accessions were assessed using descriptive statistics such as means, standard deviations, minimum and maximum values, and coefficients of variation.Pearson's correlation coefficients (phenotypic and genotypic) among the quantitative traits were estimated using the Corrplot R package [26].FactoMineR [27] and Factoextra [28] in R were used to evaluate the contributions of the quantitative traits using principal component analysis (PCA).Path coefficient analysis was conducted using lavaan and semPlot in the R package [23], considering the tuber weight per plant and dry matter content as response variables.A path diagram was constructed to depict the direct effect of key agronomic and tuber quality traits on the tuber yield and dry matter content to determine which traits can be adopted for indirect selection.A cluster dendrogram for estimating the genetic relationship among the 162 D. prehensilis was visualized using the dendextend [29] and circlize [30] packages.The variance components for each quantitative trait were estimated from the expected mean square (EMS) in the analysis of variance.The broad-sense heritability and genotypic and phenotypic coefficients of variation were calculated based on the estimated variance components as follows: Genotypic coefficient of variance where δ 2 g is the genotypic variance and δ 2 p is the phenotypic variance.Shabanimofrad et al. [31] categorized the estimated values of PCV and GCV as low for 0-10%, intermediate for 10-20%, and high for greater than or equal to 20%.Broad-sense heritability (h 2 ) was categorized as low for 0-29%, intermediate for 30-60%, and high for greater than 60%.
The multi-trait genotype-ideotype distance index (MGIDI) index theory is based on four key steps: (i) rescaling the traits so that they all have a 0-100 range, (ii) using factor analysis to account for the correlation structure and data dimensionality reduction, (iii) planning an ideotype based on known/desired trait values, and (iv) computing the distance between each genotype and the planned ideotype.
i.The following formula was used to rescale traits: where η nj and ϕ nj are the new maximum and minimum values for trait j after rescaling, respectively; ϕ oj and ϕ oj are the original maximum and minimum values for trait j, respectively, and hij is the original value for the jth trait of the ith genotype/treatment.The values for η nj and ϕ nj were chosen as follows.For the traits in which negative gains are desired, η nj = 0 and ϕ nj = 100 should be used.For the traits in which positive gains are desired, η nj = 100 and ϕ nj = 0 [32,33].In the rescaled two-way table (rX ij ), each column has a 0-100 range that considers the desired sense of selection (increase or decrease) and maintains the correlation structure of the original set of variables.
The factor analysis and ideotype design index (MGIDI) were calculated for ranking the genotypes based on multiple traits free from multi-collinearity [34].The radar chart was generated using the radar chart function of the fmsb package [35].The predicted genetic gain SG (%), was computed using the MGIDI index for each trait considering the α% selection intensity, as follows: where s is the mean of the selected genotypes, o is the mean of the original population, and h 2 is the heritability.

Quantitative Traits Variation
The interactions among seasons and accessions were significant (p < 0.05) for only the tuber flesh hardness.Significant differences (p < 0.001) among the accessions were observed across all evaluated quantitative traits.The season main effect was significant (p < 0.05) for the tuber length and tuber flesh hardness.The overall mean, range, and coefficients of variation of the evaluated traits are presented in Table 1.The responses to YMV severity based on the AUDPC score varied from 135.00 to 320.00, with an average of 149.44.The average dry matter content was 34.01%, ranging from 17.84 to 50.49%.The tuber flesh oxidation varied from −46.71 to 7.77, with an average of −13.34.The tuber flesh hardness ranged from 48.43 to 55.26 N, averaging 50.86 N. The coefficients of variation varied from 2.39% for the tuber flesh hardness to 76.70% for the tuber flesh oxidation.

Principal Component Analysis of Evaluated Quantitative
The first three principal components (PCs) accounted for ~55% of the total variation (Table 2; Figure 1 (left)).PC1 explained ~23% of the total variation, with traits such as the number of tubers per plant, tuber length, and tuber width having larger contributions to the explained variation by the PC.PC2 contributed 16.55% of the total variation, with traits such as the stem diameter, shoot internode, and dry matter content having larger contributions to the explained variation by this PC.Approximately 15% of the total variation was detected in PC3, with the number of internodes to the first branch, tuber oxidation, and tuber hardness having larger contributions to the explained variation by this PC (Table 2).The influence of the traits on the principal components and the levels of correlation among them are presented in Figure 1 (right).The PCA biplot indicates higher and positive correlations between the tuber weight per plant and the number of tubers per plant and between the dry matter content and the tuber flesh oxidation, but negative correlations between the tuber flesh hardness and the dry matter content and between the response to YMV severity and the dry matter content.

Phenotypic and Genotypic Correlation Coefficients of Quantitative Traits
The results of the phenotypic and genotypic relationships among the quantitative traits are presented in Figure 2. Significant correlation coefficients (genotypic and phenotypic) were observed between the tuber weight per plant and the number of tubers per plant, tuber length, tuber width, stem diameter, and stem internode length (p < 0.001).The dry matter content was positively correlated with the tuber flesh oxidation (p < 0.05), the number of tuber per plant (p <0.01), tuber length (p < 0.001), and the number of internodes to the first branching (p < 0.01), but significantly and negatively correlated with the tuber flesh hardness (p < 0.001), the response to YMV severity (p < 0.01), and stem diameter and stem internode length (p < 0.001).Tuber flesh oxidation was significantly and negatively correlated with the tuber flesh hardness (p < 0.001), the response to YMV severity (p <0.01), and tuber length (p < 0.001).Tuber flesh hardness revealed significant and positive correlations with the response to YMV severity (p < 0.001) and tuber length (p <0.01) (Figure 2).

Phenotypic and Genotypic Correlation Coefficients of Quantitative Traits
The results of the phenotypic and genotypic relationships among the quantitative traits are presented in Figure 2. Significant correlation coefficients (genotypic and phenotypic) were observed between the tuber weight per plant and the number of tubers per plant, tuber length, tuber width, stem diameter, and stem internode length (p < 0.001).The dry matter content was positively correlated with the tuber flesh oxidation (p < 0.05), the number of tuber per plant (p <0.01), tuber length (p < 0.001), and the number of internodes to the first branching (p < 0.01), but significantly and negatively correlated with the tuber flesh hardness (p < 0.001), the response to YMV severity (p < 0.01), and stem diameter and stem internode length (p < 0.001).Tuber flesh oxidation was significantly and negatively correlated with the tuber flesh hardness (p < 0.001), the response to YMV severity (p <0.01), and tuber length (p < 0.001).Tuber flesh hardness revealed significant and positive correlations with the response to YMV severity (p < 0.001) and tuber length (p <0.01) (Figure 2).

Phenotypic and Genotypic Correlation Coefficients of Quantitative Traits
The results of the phenotypic and genotypic relationships among the quantitative traits are presented in Figure 2. Significant correlation coefficients (genotypic and phenotypic) were observed between the tuber weight per plant and the number of tubers per plant, tuber length, tuber width, stem diameter, and stem internode length (p < 0.001).The dry matter content was positively correlated with the tuber flesh oxidation (p < 0.05), the number of tuber per plant (p <0.01), tuber length (p < 0.001), and the number of internodes to the first branching (p < 0.01), but significantly and negatively correlated with the tuber flesh hardness (p < 0.001), the response to YMV severity (p < 0.01), and stem diameter and stem internode length (p < 0.001).Tuber flesh oxidation was significantly and negatively correlated with the tuber flesh hardness (p < 0.001), the response to YMV severity (p <0.01), and tuber length (p < 0.001).Tuber flesh hardness revealed significant and positive correlations with the response to YMV severity (p < 0.001) and tuber length (p <0.01) (Figure 2).

Path Coefficient Analysis
Direct path coefficient analysis revealed that the tuber weight per plant gained significant positive contributions from the number of tubers per plant (r = 0.26; p ≤ 0.05), tuber length (r = 0.05; p ≤ 0.05), and tuber width (r = 0.03; p ≤ 0.05), while the dry matter content gained significant positive contributions from the number of tubers per plant (r = 0.80; p ≤ 0.05) and tuber flesh oxidation (r = 0.08; p ≤ 0.05) (Figure 3).

Path Coefficient Analysis
Direct path coefficient analysis revealed that the tuber weight per plant gained significant positive contributions from the number of tubers per plant (r = 0.26; p ≤ 0.05), tuber length (r = 0.05; p ≤ 0.05), and tuber width (r = 0.03; p ≤ 0.05), while the dry matter content gained significant positive contributions from the number of tubers per plant (r = 0.80; p ≤ 0.05) and tuber flesh oxidation (r = 0.08; p ≤ 0.05) (Figure 3).

Hierarchical Clustering
Hierarchical clustering based on the Gower dissimilarity matrix grouped the 162 D. praehensilis accessions into three clusters, accounting for 51, 69, and 42 accessions for clusters 1, 2, and 3, respectively (Figure 4).Cluster 1 was characterized by accessions with increased tolerance and resistance to the YMV severity, high dry matter content, and moderately low tuber flesh oxidation and hardness but low tuber weight per plant and a low number of tubers per plant (Table 3).Cluster 2 consisted of accessions that were characterized by a significantly high yield (tuber weight per plant, long, and many tubers), resistance to YMV severity, high dry matter content, moderately low tuber flesh oxidation, and low tuber flesh hardness (Table 3).The accessions in Cluster 3 were associated with high tuber yield component attributes (tuber weight per plant, tuber length, and width) and low tuber flesh oxidation.

Hierarchical Clustering
Hierarchical clustering based on the Gower dissimilarity matrix grouped the 162 D. praehensilis accessions into three clusters, accounting for 51, 69, and 42 accessions for clusters 1, 2, and 3, respectively (Figure 4).Cluster 1 was characterized by accessions with increased tolerance and resistance to the YMV severity, high dry matter content, and moderately low tuber flesh oxidation and hardness but low tuber weight per plant and a low number of tubers per plant (Table 3).Cluster 2 consisted of accessions that were characterized by a significantly high yield (tuber weight per plant, long, and many tubers), resistance to YMV severity, high dry matter content, moderately low tuber flesh oxidation, and low tuber flesh hardness (Table 3).The accessions in Cluster 3 were associated with high tuber yield component attributes (tuber weight per plant, tuber length, and width) and low tuber flesh oxidation.

Factor Analysis and Selection-Based Multi-Trait Genotype-Ideotype Distance Index (MGIDI)
The MGIDI index identified three factors based on the eleven quantitative traits.Factor analysis 1 was associated with the tuber weight per plant, the number of tubers per plant, and tuber length and width.Factor analysis 2 was correlated with the dry matter content, tuber flesh oxidation, and tuber flesh hardness, while Factor analysis 3 was associated with the stem diameter per plant, number of internodes, and stem internode length (Table 4).The average communality and uniqueness accounted for 55% and 46% of all the genetic variability in the dataset, respectively ( .57TWPL: tuber weight per plant; NTP: number of tubers per plant; TBL: tuber length; TBW: tuber width; SDMP: stem diameter per plant; NIFB: number of internodes before first branching; SINL: stem internode length; YMV: yam mosaic virus; DMC: dry matter content; TBOXI: tuber flesh oxidation; TBHard: tuber flesh hardness.Of the 11 traits evaluated, nine had desired genetic gains using the MGIDI index (Table 4).The traits with undesired selection gain using the MGIDI index were the response to YMV severity (1.92%) and tuber flesh oxidation.The MGIDI index provided a total genetic gain of 12.04% for the assessed multi-traits for which increases are desired and 3.69% for those for which decreases are desired (Table 4).

Selection of Genotypes
Of the 162 D. praehensilis accessions evaluated, the MGIDI identified 24 accessions as high-performing accessions for multiple traits (Figure 5, Table S3).These accessions show the greatest potential for the simultaneous improvement of the measured traits in a yam breeding program.
Accessions associated with FA1 (WNDpr13, WNDpr18, WNDpr4, WNDpr46, WNDpr84, WNDpr88, Asamankese002, CDpr1, CDpr28, CDpr29, CDpr4, CDpr68, and CDpr72) showed strength for traits such as the dry matter content, yam mosaic virus severity, tuber flesh hardness, and stem internode length (Figure 6; Table 5).Accessions related to FA2 (WNDpr63, WNDpr76, Asamankese009, CDpr24, and CDpr8) had strength for the tuber weight per plant, the number of tubers per plant, tuber length, and the number of internodes to the first branching (Figure 6; Table 5), while accessions associated with FA3 (WNDpr71, CDpr46, CDpr48, CDpr57, EDpr23, and PGR002) revealed strength for traits such as the tuber width, tuber flesh oxidation, and stem diameter.The strengths and weaknesses of the selected clones Figure 6.The strengths and weaknesses of the selected genotypes are shown as the proportion of each factor on the computed multi-trait genotype-ideotype index (MGIDI).The smaller the proportion explained by a factor (closer to the external edge), the closer the traits within that factor are to the ideotype.The black broken circle at the center shows the theoretical value if all the factors contributed equally.A002: Asamankese002; A009: Asamankese009.

Discussion
In this study, we employed 11 quantitative traits to assess the variations among 162 accessions of bush yam collected from three different regions of Ghana.The results show significant variations among the accessions of D. praehensilis in all the evaluated traits, indicating a high level of diversity among the D. praehensilis accessions studied.The high coefficients of variation (CV > 20) that were observed in some of the quantitative traits, especially the yield component traits, indicate huge and readily available genetic differentiation in D. praehensilis.Kouam et al. [16] reported high coefficients of variation for tuber yield components in a study conducted on D. bulbifera accessions.Similar observations of high genetic variability using quantitative traits have been reported in other yam species such as D. alata [20], D. rotundata [24], and D. dumetorum [36,37].These high variations in quantitative traits are an indication that these traits could be used as the basis for the selection of accessions with high genetic merit.The knowledge of existing variability and the degrees of association among quantitative traits are paramount for selecting superior accessions for breeding programs.
The degree of genetic variability and heritability determines the response to selection in any crop improvement initiative [38].The high GCV (>20%) and moderate broad-sense heritability (30-60%) observed for all the yield-related traits suggest high selection pressure, which could be enforced on these traits for future breeding activities.The GCV coupled with heritability estimates offers the best information about the extent of progress that can be expected from selection [39].In contrast, the low GCVs recorded for the dry matter content and YMV severity response imply that these traits can only be improved using selection methods that are not under the influence of environmental factors.Norman et al. [40] and Asfaw et al. [41] have reported high GCVs and PCVs for tuber yield-related traits and low GCVs and PCVs for the dry matter content in studies conducted on advanced and early-generation breeding populations of D. rotundata, respectively.Padhan et al. [42] also reported high GCVs and PCVs for tuber yield in a study conducted on India's wild and cultivated yam species.High H 2 observed in the YMV severity response is an indication that this trait could be improved through natural selection for superior accessions.The results from this study corroborate the findings by Agre et al. [43], who reported high H 2 in the YMV severity response in their study conducted on elite populations of D. rotundata.
The correlation analysis in this study revealed that improvement in yield and superior tuber quality traits is possible through the selection of attributing traits, such as an increased size of tubers and resistance and tolerance to YMV.The positive relationships observed among yield-related traits in this study imply that indirect selection could be adopted for significantly correlated traits.Positive correlations observed among yield-related traits have been reported by Asfaw et al. [41] and Padhan et al. [42].The significant positive correlation observed in this study between the dry matter and the number of tubers per plant agrees with the findings of Satie et al. [44] in a study on eight white yam landraces.
The cluster analysis showed that Cluster 2 was the most promising group for superior tuber yield attributes, high resistance to YMV severity, high dry matter content, and low tuber flesh hardness.Cluster 1 had some promising accessions for resistance to YMV severity and Cluster 3 had accessions with potential for high tuber yield attributes and low tuber flesh oxidation.Hybridization within each cluster may result in less genetic gain due to the close relatedness of the accessions within each cluster [45].Similarly, hybridization between accessions belonging to different clusters will result in the generation of different breeding materials.
Breeders frequently attempt to blend numerous desirable features into a new genotype in order to create high performance.It is frequently difficult to choose a genotype from the ideotype when assessing many attributes.The MGIDI was used to rank the D. praehensilis accessions based on the data from the multiple traits that were measured.The MGIDI index selected accessions (Asamankese 002, Asamankese 009, CDpr1, CDpr24, CDpr28, CDpr29, CDpr4, CDpr46, CDpr48, CDpr57, CDpr68, CDpr72, CDpr8, EDpr23, PGR002, WNDpr13, WNDpr18, WNDpr4, WNDpr46, WNDpr63, WNDpr71, WNDpr76, WNDpr84, and WNDpr88) as promising D. praehensilis accessions for a yam improvement program.The MGIDI model has also been used to assess ideal yield and yield-related variables in white Guinea yam genotypes [22], wheat genotypes [46], and eggplant genotypes [47].Multiple trait selection using the MGIDI was found to be beneficial in identifying highperforming bush yam accessions and estimating the expected genetic gains of the selected

Figure 4 .
Figure 4. Clustering showing the grouping patterns of D. praehensilis accessions using 11 quantitative traits based on the Gower dissimilarity matrix.Cluster 1-Red, Cluster 2-Green, and Cluster 3-Blue.

Figure 4 .
Figure 4. Clustering showing the grouping patterns of D. praehensilis accessions using 11 quantitative traits based on the Gower dissimilarity matrix.Cluster 1-Red, Cluster 2-Green, and Cluster 3-Blue.

Figure 5 .
Figure 5. D. praehensilis accession rankings showing selected accessions using the multi-trait genotypeideotype index (MGIDI).The selected accessions are shown as red dots, while the unselected accessions are shown as black dots.The red circle represents the cut point according to the selection pressure.

Table 1 .
Mean squares, mean performance, and variance components and broad-sense heritability of quantitative traits of D. praehensilis accessions across two seasons.

Table 2 .
Principal components and correlations among quantitative traits of D. praehensilis accessions (Bold values represent traits with high contribution to each component).

Table 3 .
Description of clusters for 11 quantitative traits of D. praehensilis accessions.

Table 4 .
Factorial loadings, communalities, uniqueness, and predicted genetic gains (PSG) based on the multi-trait genotype-ideotype distance index (Bold values represent traits with high contribution to each component).

Table 5 .
Factorial loadings, communalities, and uniqueness of twenty-four selected genotypes based on the multi-trait genotype-ideotype index (MGDI).(Bold values represent traits with high contribution to each component).