Evaluation of 130 Eggplant ( Solanum melongena L.) Genotypes for Future Breeding Program Based on Qualitative and Quantitative Traits, and Various Genetic Parameters

: Eggplant is an essential widespread year ‐ round fruit vegetable. This study was conducted using 130 local germplasm of brinjal to select diverse parents based on the multiple traits selection index for the future breeding program. This selection was performed focusing on 14 qualitative and 10 quantitative traits variation and genetic parameters namely, phenotypic and genotypic variance (PV and GV) and genotypic and phenotypic coefficients of variation (GCV and PCV), broad ‐ sense heritability (hBS), genetic advance, traits association, genotype by trait biplot (G × T), heatmap analysis and multi ‐ trait index based on factor analysis and genotype ‐ ideotype distance (MGIDI). Descriptive statistics and analysis of variance revealed a wide range of variability for morpho ‐ physiological traits. Estimated hBS for all the measured traits ranged from 10.6% to 93%, indicating that all the traits were highly inheritable. Genetic variances were low to high for most morpho ‐ physiological traits, indicating complex genetic architecture. Yield per plant was significantly correlated with fruit diameter, fruits per plant, percent fruits infestation by brinjal shoot and fruit borer, and fruit weight traits indicating that direct selection based on fruit number and fruit weight might be sufficient for improvement of other traits. The first two principal components (PCs) explained about 81.27% of the total variation among lines for 38 brinjal morpho ‐ physiological traits. Genotype by trait (G×T) biplot revealed superior genotypes with combinations of favorable traits. The average genetic distance was 3.53, ranging from 0.25 to 20.01, indicating high levels of variability among the germplasm. The heat map was also used to know the relationship matrix among all the brinjal genotypes. MGIDI is an appropriate method of selection based on multiple trait information. Based on the fourteen qualitative and ten quantitative traits and evaluation of various genetic parameters, the germplasm G80, G54, G66, and G120 might be considered as best parents for the future breeding program for eggplant improvement. and Quantitative


Introduction
Eggplant or brinjal (Solanum melongena L.; 2n = 2x = 24) is considered a rich member of the species Solanaceae, which contains approximately 1300 species. It can be grown in diversified climatic conditions of various ecological regions. It possesses high species richness with considerable flexibility of phenotypic adaptability that made the species the most important vegetable economically. Eggplant is a general term for various Solanum species cultivated for their fruits, including the East Asian aubergine (S. melongena L.) and the two African native eggplants, Scarlet (S. aethiopicum L.) and Gboma (S. macrocarpon L.) [1]. Eggplant has become prominent due to its health-promoting properties. Therefore, eggplant and its relatives have numerous medicinal applications, with 77 distinct medicinal properties [2]. Eggplant is a widespread vegetable that grown from the subtropics to the Mediterranean region, popularly in Asia, Africa, and the southern part of the USA, with significant production in 2019 globally (55.15 million tons). Asia produces more than 90% of global eggplant production with 87% of the growing area coverage [3,4]. It ranks second most-produced vegetable after potato in Bangladesh.
It occupies roughly 15% of total vegetable farmland and produces about 8% of total vegetable production [5]. However, wild forms grow in sympatry with landraces and cultivars throughout their distinct areas of origin and domestication. Natural gene flow between wild and cultivated materials, followed by natural and human selection, has resulted in intermediate phenotypes that correlate with many wild features. The contribution of features of breeding importance to diversity is unequal. The variables that contributed the most to the divergence between accessions in Indian landraces of S. melongena were yield per plant, fruit width, number of long-styled flowers per plant, flowering earliness, total phenolic content, and ascorbic acid content [6]. Nevertheless, the results depend on the sample size used. To date, there has not been any large-scale study of a representative sample containing the complete phenotypic diversity of each cultivated eggplant. In Bangladesh, eggplant is grown throughout the country, however, the yield is not sufficient due to the lack of improved and desired variety and remarkable infestation of insect pests. However, morphological characterization has been useful in studying the relationship and diversity of various eggplant varieties. The European Eggplant Genetic Resources Network (EGGNET) defined the morphological characterization for eggplant [7], which has been validated and used in the characterization of eggplant breeding materials in numerous studies [8][9][10]. Therefore, creating variation through mutation, hybridization, and biotechnology approaches is an expensive and time-consuming method [8]. Consequently, characterizing collected germplasm (populations) is required to identify lines suitable for new variety development [11]. Plant breeders are interested in genetic diversity studies based on qualitative and quantitative traits because such traits can be scored quickly and easily using low-cost methods.
The phenotypic variation of fruits, plants, and other interesting traits has been demonstrated in many articles on Solanum or S. aethiopicum or two or more eggplant species [12][13][14][15][16][17][18][19]. Moreover, summarizing the phenotypic diversity of eggplants following the Mendelian or quantitative heredity patterns of traits of interest have been widely studied in many reports [20][21][22]. Different scientists [23][24][25] studied in-depth genetic diversity, heritability and genetic advance in eggplant genotypes. Consequently, breeders face a challenge in selecting genotypes that combine high yields in multiple attributes, which requires a reliable decision support tool. In plant breeding studies, a strong selection approach can save a lot of time and resources.
The Smith-Hazel Index (SH index) is widely utilized in plant breeding as a multitrait selection index [26]. Reversing a phenotypic covariance matrix and a vector of economic weights is required to calculate the SH index. As a result of the presence of multicollinearity, poorly conditioned matrices and biased index coefficients would occur, affecting genetic gain estimations [26]. To account for the multicollinearity issue in multitrait indexes, a combination of multivariate approaches is effective in overcoming their limits. FAI-BLUP is a factor analysis-based model in which each ideotype's factorial scores are created based on desirable and undesirable elements [26]. Then, depending on the genotype-ideotype distance, a geographic probability is calculated allowing genotype ranking. Olivoto and Nardino [27] offered a new multi-trait genotype-132 ideotype distance index (MGIDI) and the entire current index is combined with the exercise calculations in the R-Metan package, which contains all the functions required for genotype selection in plant breeding programs.
In Bangladesh, however, numerous genotypes of eggplant are available. The Bangladesh Agricultural Research Institute (BARI Plant)'s Genetic Resource Center (PGRC) collects and preserves several types of eggplant germplasm from all around Bangladesh. Studying the level of accessible diversity in a crop development program is a crucial stage in crop improvement, which can be accomplished through the collection and evaluation of germplasm. Therefore, the present study was carried out to determine the inherent variation of local eggplant germplasm to identify the promising germplasm that exhibits genetic diversity for crop improvement programs through advanced multidisciplinary analysis.

Experimental Site
The experiment was executed at PGRC of BARI at Gazipur, Bangladesh during winter (Rabi season) 2019-2020 at 23.988929 N latitude, 90.412393 E longitudes and 8.40 m above sea level. The soil in the test field was silty clay with a pH of 6.

Experimental Materials, Treatments, Desing and Procedures
A total of 130 genotypes (126 locally collected germplasm and Bangladesh has developed and released four Bt brinjal varieties expressing Cry1Ac gene (Bt brinjal), viz., BARI Bt brinjal-1,2,3 and 4 variety as a check) were used in the experiment (Table S1).
All 130 genotypes were arranged in an augmented randomized complete block design (augmented RCBD) with four check varieties and seven blocks was followed in this study. All check varieties received seven replications, giving a total of 154 experimental plots. The plot size was 3 m × 2.1 m. Each genotype was implanted in three rows per plot. The spacing was 70 × 60 cm.
Direct seeding was completed with inside the well-organized seedbeds on 12 November 2019. Thirty-day-old seedlings were transplanted in the organized pits of the predominant experimental plot on 22 December 2019.
Fertilizer doses were 10 tons ha −1 Cowdung, 210 kg ha −1 Urea, 33 kg ha −1 triple superphosphate (TSP), 200 kg ha −1 , MP and 5 kg ha −1 Borax ( [28]. The total requirement of Cow-dung, TSP and Borax was applied during final land preparation about one week before transplanting. Urea and MP were supplied in the three equal splits [29]. Four times weeding and mulching were carried out in the first 25 days of mid-December. Sumithion 60 EC at 2.5 ml L −1 , Sevin 75 WP at 0.1 g pit −1 and Vertimac 18 EC at 1.2 ml L −1 were sprayed for controlling insect and mite, respectively. The data was noted as per the descriptor developed by IBPGR, 1990.

Data Recorded
Fourteen qualitative traits and ten quantitative traits (Table 1) were considered during the morphological characterization based on EGGNET [7] and IBPGR descriptors [30]. The normalized difference vegetation index (NDVI) was measured by hand green seeker (Trimble) and green seeker RT100 (Agri Optics). Quantification of two natural light sources (near infrared-NIR and red light) is measured by NDVI. These two natural lights have individual mechanisms on vegetation. For example, NIR is reflected by the vegetation, whereas red light is absorbed by the vegetation.
The NDVI formula is: The chlorophyll content in plant leaves was determined with a SPAD meter (Model: SPAD-502). The SPAD value was carried from the middle portion of the leaf of the tagged main shoot at the first flowering stage through using a self-calibrating Minolta chlorophyll meter. Measurements at each experimental plot consisted of an average of five readings [31].
Leaf area was calculated with an Automatic Leaf Area Meter (Leaf area meter-LICOR-3300, USA) at the first flowering stage of a single leaf. The leaf area index (LAI) of the crop at different growth stages was calculated using the equation as described by [32].
Data on fruit length, fruit diameter, mean number of edible fruits, and mean weight of edible fruits were measures from 10 randomly selected fruits.
The multi-trait index based on factor analysis and genotype-ideotype distance (MGIDI) proposed by [27] was used to select the novel donors with high performing under optimum and low nitrogen conditions. We also compared the result of the MGIDI index with the result of the Smith-Hazel (SH) index proposed by Smith (29) and Hazel (30) multiple trait index based on factor analysis and ideotype-design (FAI-BLUP) index proposed by [33]. The MGIDI [33] was computed as follows: where MGIDIi is the distance index of multi-trait genotype-ideotype for the ith genotype, γij is the score of the ith genotype in the jth factor I = 1, 2,..., g; j = 1, 2,..., f), g and f are the number of genotypes and factors, respectively, and Fj is the ideotype's jth score. This means that the genotype with the lowest MGIDI is closest to the ideotype, and as a result, it has the ideal values for all the analyzed features.

Statistical Analysis
Statistical analyses were performed under the R-statistics platform (software version 4.0.2) [34]. Analyses of variance (ANOVA) for each trait were assessed by using the R package 'augmented RCBD' [35]. Phenotypic and genetic variance (PV and GV), along with genotypic and phenotypic coefficients of variation (GCV and PCV), were calculated using the formula provided by [36]. Broad sense heritability (hBS), genetic advance (GA) calculation as formula elucidated in [37]. Components of the phenotypic variance of each trait were estimated using restricted maximum likelihood methods. For the estimation of variance components of linear mixed-effect ''lmer'', lme4 package was employed. R package ggplot2, scales and GGally were used for heatmap analysis.
The hierarchical clustering was performed using Spearman's rank correlation algorithm. Principal component analysis (PCA) was performed using R package ggplot2, ggfortify, usethis, devtools, plyr, scales and grid. Using a two-way matrix of 10 characteristics and 130 genotypes, a G×T biplot was constructed. The first two PCs were plotted. Genotypes were schemed according to scores on each PC, and traits were plotted based on the eigenvectors on each PC. The genotypic, phenotypic variance and broadsense heritability were estimated using Agricola R-package [38] with "metan" package. Mathematic figures were plotted using the ggplot2 package [39].

Qualitative Traits
To determine the variability of the examined germplasm, different qualitative traits were evaluated (Tables 2 and 3).

The Analysis of Variance and Frequency Distribution of Quantitative Traits
The germplasm panel, consisting of 126 accessions with four check varieties, was tested for characterized through different morphological traits including 14 qualitative and 10 quantitative traits. Analysis of variance (ANOVA) showed high significant variation among the accessions with a check for all the investigated traits. ANOVA for all the traits also revealed highly significant differences among the check, accession by check interaction as well as among the accession (Table 4)

Descriptive Statistic of the Traits
For morpho-physiological traits evaluated in this study, their descriptive statistics including means and standard error (std.Error), standard deviations (std.Dev), minimum (Min), maximum (Max), skewness and kurtosis are summarized in Table 5  Skewness is a measure of the asymmetry and kurtosis is a measure of 'peakedness' of a distribution. The skewness and kurtosis were non-significant for all the traits except PH, TF and YPP indicating all the traits fitted with a normal distribution ( Figure 2). The traits SLA, DFF, FD, TF and YPP were significant, and the distribution is positively skewed, which means that more accessions are below the mean than expected in a normal distribution. Only the trait SPAD was non-significant, and the distribution is negatively skewed, which means that more accessions above the mean than expected in a normal distribution (Figure 2). The traits PH, TF and YPP were significant and positive for kurtosis which means heavily leptokurtic distributions (Figure 2).
The results pertaining to genetic parameters viz., phenotypic coefficient of variation (PCV), genotypic coefficient of variation (GCV), broad-sense heritability (h 2 BS), and genetic advance as percent of the mean (GAM) for all the 10 traits are summarized in Table  6. TF had the highest PCV (54.76%) and GCV (53.64%), followed by FD, FL, FW, and YPP. High PCV and high GCV suggesting that these traits were under the influence of genetic control. The traits PH, NDVI, and SPAD were noted for moderate magnitudes of both PCV and GCV respectively. The traits DFF recorded for low magnitudes of both PCV and GCV, respectively.

Variability in Fruit Traits
All the traits studied related to the fruit of eggplant showed distinct variation among the germplasm except fruit cross-section and fruit position ( Table 6). The maximum variation was found in 'fruit Color at ripening stage'. Six categories of fruit Color at ripening such as milky white (1.54%), lilac grey (5.38%), Purple (24.62%), green with mottled at the distal end (39.23%), and Green with yellowish stripe (28.46%) and purple with light green at the distal end (2.31%) were observed. The next higher variation was found in 'fruit flesh density'. Very loose (spongy), loose (crumbly), average density, dense, and very dense type of fruit flesh density was found where the majority of the germplasm exhibited average density type. On the other hand, fruit calyx prickles were found as very few (1.6%), few (13.9%), intermediate (36.5%), many (23.0%), and very many (25.0%). Fruit Color distribution was exhibited as uniform (50.8%), mottled (13.1%), netted (16.8%), and striped (27.9%). No fruits were found curved in 80.56% of the germplasm. Only 8.33 % germplasm showed curved and 11.11% showed slightly curved. Fruit apexes were exhibited as two categories as rounded (40.28%) and depressed (59.72%).

Analysis of Correlation Matrix
The phenotypic correlation analysis is being used to explore a linear relationship between various traits, which was visualized in the correlation matrix ( Figure 3). In this analysis, DFF displayed a significant positive correlation with FL, while NDVI was correlated negatively. FW exhibited a strong negative correlation with TF and TF showed a negative correlation with YPP. FW and PH showed a moderate positive significant correlation with YPP ( Figure 3).

Multivariate Analysis
Multivariate analysis is a tool to find patterns and relationships between several variables simultaneously. To understand the relationship among 130 eggplant genotypes with various morpho-physiological traits, principal component, biplot, and heatmap analysis were done which revealed different clusters of genotypes that performed better in different aspects. The genotypes by traits biplot were constructed from a two-way matrix of 10 morpho-physiological traits and 130 eggplant genotypes using the relative value of the trait (Figure 4). Again, biplot analysis showed the trait profiles of the genotypes, especially, those genotypes positioned far away from the origin and the results indicated a correlation between traits with genotypes. Again, traits on opposite sides of the origin are negatively correlated and traits near each other are positively correlated. Moreover, traits at 90° to each other are not correlated, concerning the origin. The principal component (PC) analysis identified a total of 10 principal components (PCs) for the morpho-physiological traits. Among them, the first two PC explained 82.26% of the entire morpho-physiological variations (Figure 4). This biplot revealed superior genotypes with higher levels of expression of favorable trait combinations. The total outcome proposed that TF, FD, SLA, FW, PH, and YPP could help to detect superior genotypes in elite germplasm.

Heatmap Analysis
The heatmap represented the overall performance of 10 observable traits among the 130 germplasm. A heatmap is a two-dimensional data visualization technique that uses color to show the scope of a phenomenon. Color variation by hue or intensity provides the reader with a visual representation of how the phenomenon is grouped or varies over space. It depicts the relative patterns of highly abundant features against a background of mostly low-abundance features ( Figure 5). A heatmap analysis of characteristics was performed to demonstrate a chromatic examination of the genotypes. The heatmap analysis produced two dendrograms: one in the vertical direction, representing the germplasm, and one in the horizontal direction, representing the traits that caused the diffusion. The red diagonal depicts each accession's perfect connection to itself. The relationship measurements for pairs of germplasm are represented by the symmetric off-diagonal elements. Based on the morpho-physiological properties of the germplasm studied, six clusters emerged through hierarchical clustering ( Figure 5). Group (a) included 43 germplasm from clusters v (30 germplasm) and vi (13 germplasm), while group (b) included the remaining 87 germplasm from four clusters. Out of a total of 87 germplasm, cluster-i received 33, cluster-ii received 13, cluster-iii received 30, and cluster-iv received 11 ( Figure 5). Dendrogram two also revealed two significant groups: group (a) is associated with seven traits (FD, FL, FW, YPP, NDVI, SPAD, and TF), while group (b) is associated with three traits (DFF, PH, and SLA) (b). Surprisingly, the dendrogram two groups and sub-groups revealed the disparity effects of different eggplants.

Discussion
In the present study, numerous qualitative and quantitative features were assessed in 130 brinjal genotypes in order to discover superior genotypes. Genetic diversity research is critical for the successful evaluation, preservation, and use of germplasm resources [40]. The breeding strategy is mostly determined by the degree of genetic variation, and morphological characteristics are viewed as a critical initial step in characterizing and identifying plant genetic resources [41]. Screening for qualitative and quantitative features is crucial for determining a plant's socioeconomic preferences.

Qualitative Traits
All qualitative traits are found to have significant variations except leaf hairs, corolla color, fruit cross-section, and fruit position (Tables 2 and 3). Some other previous reports also published a similar type of fruit curvature [42][43][44]. Eggplant is a herbaceous plant, mostly upright in nature [45]. Our results demonstrated that 17.69% upright, 49.23% intermediate, and 33.08% prostrate growth habit among all the studied genotypes ( Table  2). Shekar et al. [46] clustered the eggplant plants into upright and intermediate. Islam et al. [42] saw 48% intermediate, 45% upright, and 7% prostrate growth habit at vegetative stage among the studied genotypes. They suggested that the plant growth habit is used for the identification of eggplant varieties.
As per leaf blade length, eggplant genotypes were divided into two groups, viz., intermediate and short [47]. They reported that 30 genotypes had intermediate leaf blade length while five genotypes showed short leaf blade length. Similarly, based on leaf blade width, the genotypes were divided into three types, viz., wide, intermediate, and narrow [45]. Osei et al. [13] found significant variation in leaf blade length and width in eggplant. In the current study, we distributed eggplant genotypes into three types, viz., very strong (5 genotypes), strong, intermediate (43 genotypes), and weak (59 genotypes), based on leaf blade lobing (Table 2). Sunseri et al. [48] reported 60% weak, 37% intermediate, and 3% strong leaf blade lobbing in eggplant. In our investigation, 6 genotypes had very acute, 28 had acute leaf blade tip angle followed by intermediate leaf blade tip angle, which includes 64 genotypes, while 32 genotypes each had obtuse. Sunseri et al. [48] observed 62% acute, 33% very acute, and 5% intermediate leaf blade tip angle in eggplant. Dash et al. [49] noticed acute leaf blade tip angle in many eggplant genotypes. A much higher variation was observed in leaf prickle (Table 2), where most of the genotype (72 genotypes) had intermediate prickle on the leaf. Sunseri et al. [48] noted that leaf prickles were absent in 25% of genotypes, whereas the remaining genotypes produced 40% very few, 25% few, and 5% both intermediate and many prickles. Tiwari et al. [50] noted prickles in stem, petiole, calyx including peduncle, and leaf including veins in eggplant. Many leaf hairs (98.46%) are found in most of the genotypes. Light violet (31.54%), pale violet (64.62%), and white (3.85%) colors were noted among the studied genotypes for corolla Color (Table 2).
In the present study, we described eight distinctive traits of eggplant fruits (Table 3). For traits related to the fruit calyx prickles, fruit color distribution, fruit color at ripening, and fruit flesh density much higher variation were observed in the cultivated eggplant (Table 3). Similar kinds of fruit calyx prickle distribution were also reported by [45]. It was noticed that 68 genotypes had no curvature on fruit; 37 genotypes showed slightly curved fruit and 25 genotypes had curved fruit. Sunseri et al. [48,51] reported sickle-shaped, snaked-shaped, curved, and U-shaped eggplant fruits along with fruits with no curvature. Eggplant genotypes were distributed into two types: rounded (39 genotypes) and depressed (91 genotypes), on the basis of fruit apex shape. Sunseri et al. [48] reported that 38, 34, and 28% of eggplant genotypes had depressed, rounded, and protruded types of fruit apex shape, respectively. Eggplant fruits traits variations are important in protection from UV irradiation, insect attack in plants as well as socio-economic value. More fruit color variation was noted by [52] and [53], which support our present study. Solaimana et al. [51] found uniform and stripped fruit color distribution. The variations in fruits color were also described by [51,54,55]. Tiwari et al. [50] divided the eggplant fruits into six color groups i.e., green (37.27%), purple (25.45%), milky white (13.62%), purple-black (12.72%), light purple or lilac grey (9.09%), and scarlet red (1.08%). Fruit flesh density is an important characteristic for the determination of fruit volume and weight. A wide range of variation was observed in fruit flesh density, the highest number of genotypes had very loose flesh density (74 genotypes).

Genetic Components
The proposed index of 0-10% for low, 10-20% for moderate, and 20% for high variation was used to characterize the projected GCV and PCV values. In our experiment, closer PCV and GCV values were estimated in most of the traits which possibly were less influenced by the environment suggesting the reliability of selection based on these traits. The assessed 10 quantitative traits exhibited a wide range of variation and more or less similar results were observed by [48,53,56]. The selection procedure considers the differences between traits based on the degree of heredity. In order to understand the predicted selection benefits, assessing genetic progress may be an important method for improving crops. Several studies have shown that selection can effectively use available genetic variations with a certain degree of heritability to improve specific traits [26,57]. The concern of both heritability and genetic advance is more effective over the distinctively use of heritability. We found that the phenotype variance values of all traits were higher than the genotype variance, indicating that the environment regulates the expression of traits. The same kind of outcome also was gained from several studies for various traits in eggplant [51,[58][59][60][61].
The evidence of the traits having essential potential in the assortment process due to low environmental impacts was calculated using strong GA with hBS for all yield-related traits except DFF. Heterosis breeding has the potential to improve traits with poor heritability and genetic advance [28]. Regarding the hBS and GA index [62], which were greater than 60% for high, 30-60% for moderate, and 0-30% for low, we discovered that all traits were highly heritable together with a high genetic advance mean value except DFF, implying that direct selection can be effective for eggplant crop improvement based on these traits with the effect of additive genes; a more or less similar result was obtained by many researchers [63,64]. Due to the strong influence of the environment on genetic effects, low to moderate heritability and genetic advanced values will inhibit the improvement of traits. Therefore, effective selection can only be achieved by selecting higher values of GCV, PCV, hBS, and GA, which means that the influence of additive genes is more stable than the influence of the environment. High GAM was also observed for all traits except DFF. This revealed that if the selection was carried out for the next generation for these features, a greater improvement in the population mean may be seen.

Correlation Matrix
The plant breeding correlation matrix is an outstanding method to assess the relationship between two or more variables. For higher genotype selection procedures, considering the correlation matrix can be a scaling measure [65]. In the correlation matrix, DFF displayed a strong positive correlation with FL, while NDVI was correlated negatively, FW and PH showed a moderate positive correlation with YPP ( Figure 6). A moderate to positive significant association can be proposed; selection based on these traits can help increase the yield of this crop. More or less similar results were consistent with the study of [28,64,66,67].

Multivariate Statistical
Compared with univariate and bivariate statistical methods, multivariate statistical methods can analyze more than one relationship at the same time. There are many multivariate data analysis methods, each of which has a different purpose, such as Regression analysis, factor analysis, cluster analysis, analysis of variance, discriminant analysis, etc. [68]. Biplot analysis is usually employed to assess the component effects creating the genotypic variations.
The highest values indicate the highest influence of the trait on the total variation. Biplot analysis determines varietal stability in the multi-environmental trial [69]. It describes the relationship between different genotype traits. The association between morpho-physiological traits among the 130 genotypes was observed by the biplot analysis [70]. Again, the biplot analysis showed the trait profiles of the genotypes, especially, those genotypes positioned far away from the origin and the results indicated a correlation between traits with genotypes ( Figure 4). An acute angle between two elements indicates a positive correlation, and an obtuse angle between two elements indicates a negative correlation. As a result, principal components (PC) analysis provides a good screening of available genotypes and aids in the selection of possible parents for crop breeding initiatives. In our data, the first two PC accounted for 82.36% of the overall variation (Figure 4). The yield potential of accessions was represented in PC 1; thus, the accessions contributing to this component are likely to undergo direct selection, or selected parents can be used in hybridization operations. These PC1 results are consistent with the results of the correlation analysis. The figure summarizes the information of the matrix in principal components, where the cosine of the angle between the vectors connecting the objects to the origin is proportional to the correlation coefficient between these objects.
The heatmap shows the highest and lowest values of each genotype in different colors against all the traits comparing. The intensity of the color indicates the degree of high or low of the traits. Hierarchical clustering based on the morpho-physiological traits of the studied germplasm revealed six clusters ( Figure 5). The heatmap analysis depicted the degree of correspondence among the morphological traits assessed in brinjal genotypes, and this result was consistently supported by [71,72].

Multi-Trait Index Based on Factor Analysis and Genotype-Ideotype Distance (MGIDI)
Experienced breeders often try to combine several desired traits into a new genotype to produce high performance. When measuring multiple traits, it is often difficult to select a genotype from the ideotypes. In this regard, various multivariate methods are widely used, such as principal component analysis, factor analysis, cluster analysis, and different samples to group measured traits or select test genotypes [73]. We used a two-way heat map clustering pattern and PCA to connect test genotypes and measured attributes in this study (Figures 6 and 7), however, we could not pick specific genotypes. To make the selection of genotypes with several features easier, [27] the recently introduced MGIDI (multi-trait genotype-ideotypes distance index) is a new method for genotype selection based on multiple trait information. The eggplant genotypes were ranked based on information on measured multiple traits (Figure 7). The MGIDI index selected genotypes G80, G54, G66, G120, G46, G61, G65, G108, G4, G79, G42, G77, G47, G50, G51, G43, G44, G48, and G49 as promising eggplant genotypes. Apart from these genotypes, G80 was very close to the cut point, which recommends that this genotype can exist desirable features. Hence, the researcher should pay particular attention to assessing genotypes that are very close to the cut point [27]. The application of the MGIDI index to plant crop research is predicted to grow rapidly. Similarly, This index was used to find the best strawberry genotype [74].

Conclusions
The current study clearly established that the improvement of brinjal yield and related traits can be obtained through selection with the valuation of different genetic parameters analysis. Most of the qualitative traits showed distinct variations among the germplasm. Qualitatively, the maximum variation was observed in leaf prickles, fruit calyx prickles, fruit color at ripening, and fruit flesh density. Quantitatively, the highest variation was observed in fruit yield per plant which was followed by fruit weight. Nevertheless, the present study has shown that selection with the evaluation of various analyzes of genetic parameters such as GCV, PCV, hBS, and GA can achieve an improved eggplant yield and related traits. It can be observed that practically all of the agronomic traits tested in this study exhibit significant variability based on the recorded data and additional analyses (heatmap analysis, correlation matrix, PCA, MGIDI analysis). Therefore, the germplasm G80, G54, G66, and G120 might be considered as best parents based on the qualitative and quantitative characters for the future breeding program. The present findings have a great genetic potential for the studied germplasm. Concurrently, the promising germplasm identified in the present study might be used in future breeding programs for eggplant improvement.