The study material consisted of 2,324 cultivated lentil germplasm accessions comprising 1,796 indigenous collections (ICs) and 528 exotic collections (ECs). The indigenous accessions were collected from different lentil-growing areas in India after exploring several crops from 1976 to 2017. Geo-coordinates of each collection site were mapped using geographic information system (GIS) tools (Figure 1). Geo-referenced maps were prepared with WGS84 datum and geographic projection system using GIS tools (Semwal et al., 2014). The germplasm was characterized during the winter season of 2017–2018 and 2018–2019 at Research Farm of ICAR-NBPGR located at a latitude of 28°38′N and longitude of 77°10′ E and an altitude of 228.61 m above the mean sea level. The soil of the location was sandy loam type having an optimal pH range. Lentil germplasm accessions were sown under natural field conditions. The recommended package of practices for growing lentil was followed. To meet the crop’s balanced nutrient demand, 20 kg ha^–1 nitrogen and 40 kg ha^–1 phosphorus was applied as a basal dose before sowing. One presowing irrigation was given to ensure proper germination in the field, while irrigation was also given at the pod formation stage. The seeds were treated with a mixure of thiram (2g) and carbendazim (1 g) per kg of seed. Manual weeding was completed at 25–30 and 45–50 days after sowing. The experimental design consisted of an augmented block (Federer and Raghavarao, 1975) with four popular checks, namely, DPL62, IPL316, IPL526, and L4717. The experiment was replicated in 24 blocks. A hundred accessions were sown per block uniformly in 23 blocks, whereas, the 24^th block was incomplete with 24 accessions. The accessions were sown in paired rows of two m length, with a row-to-row distance of 30 cm. The detailed timeline of the experiment is shown in Figure 2.

Twenty-six agro-morphological descriptors (nine quantitative and 17 qualitative) were selected for phenotypic characterization using descriptor lists of lentil (IBPGR, 1985; Mahajan et al., 2000) (Table 1). The qualitative characters were (1) early plant vigor (EPV), (2) seedling stem pigmentation (SSP), (3) growth habit (GH), (4) leaf color (LC), (5) leaf pubescence (LP), (6) leaflet size (LS), (7) flower ground color (FGC), (8) tendril length (TL), (9) biomass score (BS), (10) lodging score (LoD), (11) pod pigmentation (PP), (12) pod dehiscence (PD), (13) pod shedding (PS), (14) seed shape (SS), (15) seed coat color (SCC), (16) pattern of seed testa (PST), and (17) cotyledon color (CC). The quantitative traits included (1) days to 50% flowering (DF), (2) days to 80% maturity (DM), (3) plant height (PH), (4) number secondary branches per plant (SBP), (5) seed diameter (SD), (6) seed thickness (ST), (7) pods per plant (PPP), (8) seeds per pod (SPP), and (9) 100-seed weight (SW).

A statistical software, PowerCore, was used to constitute a core collection based on the characterization data generated using 26 different descriptor traits (Kim et al., 2007). Continuous variables were classified into different categories based on Sturges’ rule (Sturges, 1926) (K = 1 + log₂n), where n and K mean the number of observed accessions and number of classes, respectively. The number of classes can be modified. The modified heuristic algorithm was implemented to identify and select accessions of core collections using the strategy of advanced maximization or allelic richness. Subsets created using PowerCore represented all observation classes with the least allelic redundancy, ultimately ensuring a highly reproducible entry list. Core sets were extracted with different combinations of variables to assess their efficacy using PowerCore and using qualitative traits, quantitative traits, a combination of quantitative and qualitative traits, and a combination of quantitative and qualitative traits and passport data. Another method, known as principal component score strategy (PCSS), was used to extract the core (Noirot et al., 1996). This method uses principal component analysis (PCA) to eliminate collinearity between variables and sample individuals based on their cumulative relative contribution. The Newman–Keuls procedure was followed to compare the means of the entire collection and developed core sets (Newman, 1939; Keuls, 1952) for descriptor traits. The similarity of distribution frequencies in the core and entire sets was evaluated using the χ² test. The Shannon–Weaver diversity index (SDI) was used to compare the representativeness of the whole set and the germplasm selected as core entries.

The data recorded for nine quantitative traits were first analyzed separately for each year, after which a combined analysis was performed using the general linear model (PROC GLM) method of SAS 9.3 for augmented design. The accessions were considered random and year as a fixed effect. Levene’s test was performed to assess the homogeneity of error variances before performing the combined analysis (Levene, 1960). Descriptive statistics were estimated for the complete set and the core collections separately. Frequency distribution for different classes of qualitative traits was obtained using Excel, and histograms for quantitative traits were drawn using IBM statistical package for the social sciences (SPSS) statistics (version 20.0). Variability in descriptor traits between ICs and ECs was compared using boxplots developed using Statistical Analysis Software-JMP 14 software. Hierarchical cluster analysis was performed using the Euclidean distance matrix following Ward’s minimum variance method for accession number grouping and estimating genetic relationships. Core sets were extracted using PowerCore and the PCS strategy. The principles laid by Oliveira et al. (2010) were used to calculate the range retention percentage. In addition, the distribution of homogeneity for each of the descriptor traits was analyzed using the χ² test. For the quantitative traits, classes were formed based on Sturges’ formula. The observed number of accessions in the core set in each class was determined and was tested against the expected number of accessions using the χ² test. Further, the SDI (H’) was computed as a measure of phenotypic variability using the phenotypic frequencies of quantitative and qualitative traits (Shannon and Weaver, 1949) in the entire and different core sets extracted using the following formula.

where p_i is the proportion of germplasm accessions in the i^th class of an n-class characteristic and n is the number of phenotypic categories for a character. Phenotypic correlation coefficients (r) among descriptors in the core collection (Core^d) and entire set were estimated. Moreover, the contribution of different descriptor traits to multivariate polymorphism and conservation of contribution in the core collection were assessed based on PCA. The SAS 9.3 program was used to estimate the correlation and PCA (SAS Institute, 2012).

The actual germplasm collection sites of lentils in India (Figure 1) were presented using geo-referenced map. These maps revealed that most of the lentil germplasm were collected from Uttar Pradesh, Bihar, Uttarakhand, and Madhya Pradesh. Geo-referenced map showed that lentils were collected from Western Himalayas (Kashmir, Himachal Pradesh, and Uttarakhand) to Eastern Himalayas (Assam and Arunachal Pradesh). It was collected from high-altitude areas in north (Leh) to low-lying areas of eastern part of India (West Bengal). Collections from dry areas of Rajasthan to wet areas of North-east Himalayas indicated that India has sufficient diversity in its collections. Indian lentil collections showed the unique distribution pattern having representation from four global biodiversity hotspots, the Himalayas, the Western Ghats, the Indo-Burma region, and the Sundaland (Includes Nicobar group of Islands). All 2,324 accessions and four checks were characterized under the agro-climatic region in the north-western Indian conditions for two consecutive years (2017–2018 and 2018–2019). The experimental field view and seed coat color diversity are presented in Figure 3.

The frequency distribution of the entire set was not normal and varied among traits. In certain cases, a particular trait was predominant, whereas it was dispersed in other cases (Figure 4). The majority of the accessions (2,061) showed purple pigmentation on the stem at the seedling stage, whereas leaf color was recorded as green. A slight to medium leaf pubescence was prominent, whereas 48 accessions were glabrous. EPV of most of the accessions (2,044) was good, whereas 275 accessions were highly vigorous. Prominent tendril was recorded in 1,360 accessions, whereas 954 accessions showed rudimentary tendril. High biomass was observed at the mid pod filling stage in 1,438 accessions, whereas 72 accessions showed low biomass. Semi-erect plant growth was more prominent (1,648 accessions), and horizontal growth was recorded in 525 accessions, whereas the growth in 151 accessions was erect. Around 67% of the accessions showed lodging at the time of pod filling. The maximum number of accessions (2,159) were of purple flower ground color. After a week of maturity, medium pod shedding was recorded in 180 accessions, whereas most of the accessions (1,671) showed no pod shedding. Pod dehiscence was absent in 1,765 accessions, whereas low to medium pod dehiscence was noted in 549 accessions. Pod pigmentation was absent in the majority of the accessions, whereas brown seed coat color and dotted testa pattern were prominent. A flattened seed shape was observed in 2,060 accessions, whereas 254 accessions showed the globose seed shape. There was substantial variation in the cotyledon color with a maximum of orange type.

The results revealed significant variations as evident from the mean, range, and coefficient of variations (Supplementary Table 1 and Figure 5). The days to 50% flowering and maturity ranged from 51 to 123 and 93 to 140 days, respectively, whereas 100-seed weight varied from 0.75 to 7.6 g. The average seed thickness and seed diameters were 2.42 and 4.18 cm, respectively. The mean plant height and the number of secondary branches were 35.34 and 28.64 cm, respectively. The coefficient of variation ranged from low (<10.17%; days to flowering, seed thickness, days to maturity, and number of seeds per pod) to high (>17.34%; the number of secondary branches, 100-seed weight, number of pods per plant, and seed diameter) for key descriptor traits. In addition, substantial variation in the range for the number of pods per plant, number of secondary branches per plant, 100-seed weight, and earliness was observed. The least variation was observed in the number of seeds per pod, with the majority of the accessions with two seeds per pod in studied accessions. ICs showed earliness, small seeds, and low biomass as compared to the ECs. The days to 50% flowering, days to 80% maturity, and seeds per pod were negatively skewed, indicating that a majority of the accessions had higher mean values for these traits and days to 80% maturity and seeds per pod showed high kurtosis values with a relatively higher peak. Boxplot analysis revealed a comparison of trait distribution between exotic and indigenous accessions. The average performance for plant height, days to 50% flowering, pods per plant, and number of secondary branches was lower in ECs than in ICs. However, the mean value of seed weight was higher in ECs than in ICs (Figure 6).

Grouping several accessions into certain homogenous clusters assists in selecting diverse parents. This allowed the accurate comparison among all probable pairs of individuals and brought together gene constellations, yielding desired progenies through crossing between different parents. Hierarchical cluster analysis with Ward’s method of minimum variance was applied to assess genetic diversity. Based on nine quantitative traits, phenotypic relationships among lentil accessions were determined using Euclidean distances. The Indian national gene bank (INGB) lentil accessions were grouped into 12 clusters with a varied number of accessions in each cluster (Figure 7 and Supplementary Table 2). A list of core set accessions with cluster number is given in Supplementary Table 3. The cluster also showed heatmap distribution of accessions, showing trait variability in different clusters represented by different colors. Cluster II was the largest with 393 accessions followed by Clusters III and IV with 382 and 253 accessions, respectively. Cluster II grouped 320 ICs and 73 ECs accessions; whereas Cluster XI was the smallest with only 38 accessions (24 ECs and 14 ICs). The cluster mean values of studied traits showed that Clusters XI and XII included extra-early accessions, such as IC241529, IC241531, and IC241532, whereas early to medium maturing accessions were grouped into Clusters VI and VII. Late flowering accessions such as IC73690, IC73692, IC381127, and IC398044 grouped in cluster I. Proportionately, exotic accessions were relatively dominant in Clusters X, XI, and XII, and indigenous were dominant in Clusters II to IX. Cluster XII comprised all exotic accessions from Syria. Cluster V contained accessions having very tall plant types (IC329110, IC398793, and IC59038, a high number of pods per plant (IC78387 and IC78398), and a high number of secondary branches (IC148333, IC16453). The accessions such as EC223214 and IC241543 grouped in cluster XII had highest mean value for seed weight (3.239 g) whereas cluster III had accessions with mean value of 1.523 g.

Potential accession numbers for unique traits were identified from characterization datasets and validated over the years and locations at NBPGR, New Delhi, The International Center for Agricultural Research in the Dry Areas (ICARDA), Amlaha, and NBPGR, Regional Station, Ranchi. Some of them are depicted in Figure 8. The accession IC317520 was identified with a unique seed morphotype having extended funiculus. A total of 50% of plants with accessions such as IC241532 and IC241529 exhibited very early flowering (51 days) relative to all checks used in the study. The accession EC267615 displayed a novel type of erect plant architecture, indicating its suitability for mechanical harvesting. Pod setting in this genotype was observed at more than 15 cm above the ground. The accessions IC199461 and IC201716 had a high number of secondary branches (>50) and pods per plant (>400) relative to all checks used in the study. The accession EC499760 was of bold seeded type (100-seed weight, 7.1–7.83 g) and was imported from the United States of America (USA). A novel and unique multiflower (MF) germplasm accession, IC241473, was identified the first time that formed up to 16 flowers per peduncle (FPP) at multiple flowering nodes. In addition, this accession showed fasciation of the main stem.

The core sets were extracted using combinations of different datasets using the PowerCore software. Additionally, the PCSS was used to develop the core sets. These core sets were designated as Core^a, Core^b, Core^c, Core^d, and Core^e when generated using qualitative traits, quantitative traits, a combination of both traits, a combination of both traits and passport data, and PCSS. A total of 170 accessions were selected based on the trio of qualitative and quantitative traits and passport data using PowerCore (Table 2 and Supplementary Table 3). A total of 67 accessions were selected when the data of only qualitative variables were used, whereas 71 accessions and 79 accessions were selected, respectively, when data of only quantitative characteristics and both variables were used for the analysis. Accessions with the maximum relative contribution corresponding to 10% representation were used as the criteria to select the core; thus, 221 accessions were selected based on the PCS strategy. Relative SDI was estimated for both quantitative and qualitative traits in entire lentil accessions and in core sets developed using a different strategy that revealed adequate diversity. The lentil germplasm was more diverse (H’ > 0.5) for all descriptors, except for EPV, SSP, FGC, PP, SS, CC, and SP. The average SDI for quantitative traits was higher than qualitative traits. The SDIs were consistently higher in the core sets derived from different datasets than in the entire set, signifying an enhanced representation of the available diversity in each core set. PowerCore yielded a better representation of diversity than PCSS. The average SDIs for all traits revealed that Core^d (1.361, 170 accessions) was more diverse, followed by Core^b (1.340, 71 accessions) and Core^a (1.302, 67 accessions) (Table 2).

Qualitative traits in the INGB collections were represented in the core sets developed by all strategies, except the core developed by quantitative method and PCSS strategy, indicating that the Core^a, Core^c, and Core^d captured the allelic richness of the entire INGB collections. The range of quantitative descriptors was classified into 12 groups for frequency distribution and χ² analysis based on PowerCore. Frequency distribution analysis indicated the homogeneity of distribution of several traits by Core^a(10), Core^b(13), Core^c(11), Core^d (18), and Core^e(9) (Table 3). Thus, the core set developed by combining qualitative and quantitative traits and passport data well-represented the structure of the entire lentil germplasm. The results showed that the variability existing in the entire set was well-symbolized in the core collection. Differences in the means of the entire collection and Core^d were insignificant for quantitative descriptors, except days to flowering and seed weight. In addition, the range retention percentage was 100, capturing a complete range of diversity, including all extremes (Table 4).

The lentil germplasm accessions recorded significant and positive correlations between seed thickness and seed diameter, equally in both entire set and core set (r = 0.539) and between days to 50% flowering and days to 80% maturity (r = 0.582 in the entire set, r = 0.731 in the core set). Plant height was positively correlated with days to 50% flowering (i = 0.378), days to 80% maturity (r = 0.287), pods per plant (r = 0.496), seeds per pod (r = 0.521), and number of secondary branches (r = 0.542) and negatively with seed thickness (r = –0.034) and seed diameter (r = –0.234) in the core set. Similar trends were recorded in the entire lentil collection. The number of secondary branches was positively correlated with days to 50% flowering, days to 80% maturity, plant height, pods per plant, and seeds per pod negatively with 100-seed weight, both in the core collection and entire set (Figures 9, 10).

Principal component analysis revealed the association among different descriptors and their contribution toward variability. The first three PCA components provided a realistic summary of the data and explained 70.69% of the total variation in the core set. The first principal component (PC1) accounted for 39.36% of the total variation, whereas PC2 and PC3 accounted for 17.73% and 13.60% variance, respectively. These values were comparable with those of the entire lentil collection, explaining 68.03% of the variation (PC1:37.49%; PC2:16.62%, and PC3:13.91%). The negative variation in PC1 was mainly contributed (>50%) by days to 50% flowering, days to 80% maturity, plant height, number of secondary branches per plant, and pods per plant and 100-seed weight in the entire and lentil core set (Supplementary Table 4 and Supplementary Figure 1).

Geo-referenced maps of India revealed that most of the Indian accessions were from the Gangetic plains and eastern part of India belonging to Uttar Pradesh (254), Bihar (203), Uttarakhand (181), and Madhya Pradesh (96). These are the most populated regions of India, with the massive scope of pulse production in their rice fallow areas (Chowdhury et al., 2020). The majority of the germplasm was collected from hot subhumid (dry and moist) to hot subhumid agro-ecological zones of India. In contrast, the ECs, conserved in INGB, represented all six continents, primarily Asia.

The major steps in crop improvement are assessing the variability for desired traits and their utilization in breeding programs (Kumar et al., 2015). Previous researchers adopted the piecemeal approach to characterize the agro-morphological traits in lentil germplasm collections (Gautam et al., 2013; Jha et al., 2014; Kumar and Solanki, 2014). However, in this study, the entire genebank collections representing all agro-ecological zones were characterized in one go. Significant differences among lentil accessions based on agro-morphological characterization indicated their potential for use in current and future lentil breeding programs. The results showed significant variations in certain traits such as days to flowering, seed thickness, days to maturity, number of secondary branches, 100-seed weight, number of pods per plant, and seed diameter. The wide range of variability obtained was attributable to the diverse collection assessed, involving native and introduced germplasm originating from diverse geographical origins with different genetic makeup (Phogat et al., 2020).

Lentil is a less vigorous grain legume that faces competition from weeds during the early growth stage. Lentil germplasm with increased vigor can dominate the weeds, ultimately enhancing the yield (Sharma et al., 2018). The previous researchers reported that EPV and rapid canopy development as key traits for escaping drought at the terminal stage because these traits lead to the early onset of maturity (Sarker et al., 2005; Kumar et al., 2012). The majority of exotic accessions had light-green leaves, whereas indigenous varieties were dominated by green to dark green canopy. In addition, the purple pigmentation was observed to be disappeared with the increase in the ambient temperature. Similarly, in subterranean clover, Nichols et al. (1996) reported prominence of purple pigmentation during winter, which is appeared with the arrival of spring. Wide variations were recorded in leaf pubescence, ranging from the absence of pubescence to dense pubescence. Low or moderate leaf pubescence has promoted tolerance to aphids (Kumari et al., 2009). However, Tripathi et al. (2020) indicated that pubescence on the pod is a unique feature of cowpea, a marker-trait for insect resistance. Variation was also reported in leaflet size, tendril length, flower ground color, pod shedding, pod dehiscence, testa color, and seed shape.

Elias et al. (1979) suggested that the protein quality is affected by seed coat color in beans. Diversity in seed color among accessions may provide valuable genetic resources for biofortified lentils. Orange cotyledons are preferred by Indians, which is also reflected in Indian genebank collections (Singh et al., 2014). These agro-morphological descriptors can facilitate the differentiation of distinct phenotypic classes and are used as diagnostic keys for taxonomic delineation (Pundir et al., 1985; Gore et al., 2019). The mode of inheritance of these traits can be investigated using the principles of classical genetics. This mega characterization program showed that most of the Mediterranean germplasm were found with low biomass and poor yield. The direct use of these accessions appeared difficult in Indian conditions, leading to the ubiquitous selection of Indian germplasm with high biomass and better yield traits.

The entire lentil germplasm was more diverse for quantitative traits such as number of secondary branches per plant, pods per plant, plant height, days to 50% flowering, seed thickness and qualitative traits, leaf pubescence, biomass score, pod shedding, growth habit, leaflet size, and pod dehiscence. Plant height (cm) in the INGB lentil core set ranged from 9.70 to 52.13 with mean value of 34.04, while it was ranged from 15 to 40 cm with mean value of 23 cm in the USDA core set sown at Pullman, Washington, United States, in 1996 (Tullu et al., 2001). An average value of days to flowering and maturity is lower in the USDA lentil core set than the INGB lentil core set (Tullu et al., 2001). The grouping of germplasm in different clusters allowed us to study their genetic diversity, whereas the heatmap depicted the trait diversity. Germplasm accessions belonging to diverse clusters can be used in cross-breeding programs and recombination studies; for example, Clusters XI comprised early maturing accessions such as IC241529, IC241531, and IC241532 (less than 100 days), and cluster X comprised late maturing accessions such as EC267554, EC267615, IC258265, and EC440747 > 125 days. Similarly, Clusters XI had accessions with low plant height IC241488, IC241501, IC241529, and IC241531 (<20 cm) and cluster V had accessions with larger plant height, IC329110, IC398793, and IC59038 (>50 cm). The accessions such as EC223214, IC241543 (with high 100-seed weight) grouped in cluster XII and may be crossed with accessions with low seed weight IC361467, IC521438 (<1.4 g) from cluster III for genetics and mapping purposes.

The ultimate aim of germplasm characterization is to support crop-breeding programs (He and Li, 2020; Guerra-García et al., 2021). Similarly, the present work unfolded the treasure of lentil collections conserved in the NGB. The accession IC317520 collected from Rajasthan was identified as a unique seed morphotype with intact funiculus at maturity (Tripathi et al., 2019). This was reported for the first time and is registered as a genetic stock (INGR19072) for its utilization in lentil breeding programs. An Intact funiculus may provide resilience to drought stress, resulting in high yields. This trait may help to gain area in rice-fallow areas of eastern India, where the crop is sown at conditions of conserved moisture (Lamichaney et al., 2021; Tripathi et al., 2021). Rice–wheat and rice–rice systems in the Indo-Gangetic plains dominate South Asia’s cereal systems (Kumar et al., 2013). A substantial reduction in the acreage under pulses is attributed to the expansion of high-yielding wheat and boro rice in northern India and eastern to north-eastern India, respectively. This trend of pushing out pulses is likely to continue unless extra-early varieties that fit into these cropping systems are developed (Erskine et al., 1993; Shrestha et al., 2006). Similarly, in the rice-fallow areas of South Asia, the top soil layer generally dries at the harvesting stage, making it impossible to sow the next crop. Under such conditions, extra-early lentil germplasm can be used as donors to breed early maturing varieties to fit into the cropping system and convert these mono-cropped areas into double-cropped areas, thereby increasing the production of lentils under rice-based systems. In a previous study, a lentil accession with 59 days to flower was reported (Kumar and Solanki, 2014). However, the present experiment discovered accessions, IC241529 and IC241532, flowered in 51 days. The number of pods per plant and the number of secondary branches are the primary quantitative descriptor traits that can be used as kickstarting points for selecting high-yielding germplasm. Further, the seed size and shape of lentils are important traits because they affect the market class, cooking time, and the quality and yield of milled lentils (Fedoruk, 2013). Indigenous and exotic accessions are broadly classified into small- and large-seeded types, respectively. Descriptor traits, such as tall, erect, top pod-bearing habits, thick stems, synchronous maturity, lodging tolerance, and non-shattering reported contributing to the development of varieties with machine harvestability (Kumar et al., 2013). Genetic variability for these traits was recorded in the INGB lentil germplasm. In an earlier study, Idleb-2 variety was released for machine harvesting (El-Ashkar et al., 2004). Identified germplasm, EC267615, may provide avenues for breeding machine harvestable varieties in India. Multiflowering accessions can be identified with the expression of three or more flowers at one or multiple flowering nodes (Gaur and Gour, 2002; Benlloch et al., 2015; Sanwal et al., 2016; Devi et al., 2018). In several legumes, previous researchers have documented multiflowering clusters for a few flowering nodes in Cicer (up to nine FPP) (Gaur and Gour, 2002), Pisum (up to five FPP) (Devi et al., 2018), and Lens (up to seven FPP) (Mishra et al., 2020). However, this study found a unique multiflowering germplasm accession, IC241473, in cultivated lentils, formed up to 16 FPP at multiple flowering nodes. This trait can fulfill the aim of genebanks to work in the direction of conservation and utilization. Trait specific germplasm (TSG) developed in this study were shared with the National Agricultural Research System (NARS) partners under the material transfer agreement for utilization in the varietal development programs.

Globally, the first and foremost reason behind the limited use of conserved germplasm by plant breeders is huge and uncharacterized germplasm collections in certain genebanks (Ortiz et al., 2008). Characterized germplasm and smaller subsets can enhance germplasm use in trait discovery (Dutta et al., 2015). Core sets are subsets of large and diffused germplasm collections, containing selected accessions that represent the genetic variability of the entire set and acts as the entry point to the entire collection for trait-specific evaluation (Ortiz et al., 1998, 1999; Huamán et al., 1999; Upadhyaya et al., 2003). PowerCore is a fast and reliable tool for extracting a core set with a higher percentage of diversity than other methods used (Kim et al., 2007; Zhang et al., 2011; Dutta et al., 2015).

A subsampling procedure based on 26 agro-morphological descriptors and geographical origin data resulted in the selection of 170 germplasm accessions as the lentil core set, whereas only 79 accessions were sampled without passport information. Therefore, the inclusion of stratification into the geographical origin and places should be a choice for the sampling procedure. The INGB lentil core set amounted to 7.31% of the entire collection (2,324 accessions) conserved in the INGB, New Delhi. Similarly, Archak et al. (2016) selected 1,103 accessions (7.5%) as the chickpea core set from 14,651 accessions conserved in the INGB. As expected, the INGB lentil core set comprised 80.58% accessions belonging to Indian origin. However, only 12.89% (37 accessions) of the Indian accessions contributed to a lentil core set developed for agro-morphological descriptor traits by United States Department of Agriculture (USDA) (Tullu et al., 2001). The USDA lentil core consists of 12% of the 2,390 lentils (Simon and Hannan, 1995; Tullu et al., 2001). The maximum representation of germplasm accessions in INGB lentil core collections was from India (137), followed by Syria (22), indicating the rich diversity of cultivated lentils in India. A total of 95.88% of the core set belonged to Asian countries, including 15.29% from other Asian countries (Syria [mainly from ICARDA], Bangladesh, Israel, and Pakistan). As lentil is the main crop of Asia (Laskar et al., 2019), the INGB core set can be used as a diversity hub to meet the present and future challenges of the Asian lentil breeding program. Simon and Hannan (1995) also reported the utility of the legume cores in chickpea, lentil, and pea developed by USDA for directing users toward desirable germplasm from targeted geographic regions, and facilitating users at the preliminary stages of germplasm evaluation.

While sampling for a core set, the phenotypic associations in the core set need to be preserved for maintaining coadapted genetic complexes (Ortiz et al., 1998) and efficient utilization of germplasm (Upadhyaya et al., 2001). Nine traits were studied and 36 comparisons were performed between the pairs of characteristics. All trait combinations maintained similarity in correlation coefficient value between the entire INGB collection and the INGB core set. Earlier studies on the development of core in various crops such as chickpea (Upadhyaya et al., 2001; Archak et al., 2016), eggplant (Gangopadhyay et al., 2010), and wheat (Phogat et al., 2020) also reported the preservation of trait association in the core set. The strong correlation among some of the traits, such as days to 50% flowering and days to maturity (r = 0.52 in the entire collection, r = 0.731 in core set), and secondary branches per plant and pods per plant (r = 0.63 in the entire collection, r = 0.65 in core set), showed that future germplasm characterization might use days to 50% flowering and secondary branches per plant during preliminary evaluation. This trait is easier to measure than other traits, such as days to maturity and pods per plant. Similar observations were reported by Upadhyaya et al. (2003). Examination of the entries’ spatial distribution and explaining the variance through PCA were also reported as an exploratory criterion for evaluating the core set in the study of Bisht et al. (1998) and Wang et al. (2008). In the present study, the first three PCs collectively explained 70.39% of the total variation in the core sets and 68.03% in the entire set. Germplasm accessions were majorly distributed in one large group and one small group in both the complete set and the core collection of INGB lentil germplasm.