Abstract

Understanding the genetic relationships and diversity of Ethiopian lentil in relation to lentil from other countries is important in attempting to widen the genetic base of germplasm in the country. The objectives of this study were to generate information on agromorphological variability, to estimate PCV, GCV, heritability, and expected genetic advance of quantitative traits of lentil. 228 genotypes with different population types were studied for 11 agromorphological traits and rust disease severity score for two seasons (2011-2012) over three locations. The analysis of variance showed highly significant variations ( ) among genotypes for all characters studied. As per genetic parameter values, four groups of character were deduced. It is inferred that the exotic genotypes introduced from ICARDA showed rich genetic bases for 100-seed weight, number of seeds per plant, seed weight per plant, resistance source for rust, and high yielder in high yielding environment, where rainfall is not a major problem. Use the Ethiopian accessions for developing cultivars that could be used in double cropping and drought prone areas.

1. Introduction

Lentil (Lens culinaris Medik. subsp. culinaris) is a short, slender annual cool-season food legume. The center of origin of lentil is the Near East where it was first domesticated in the Fertile Crescent around 7000 BC [1]. According to Barulina [2], two varietal types based on seed size were recognized: the small-seeded microsperma and the macrosperma with relatively large flattened pods. The largest lentil producer is India, followed by Canada and Turkey, accounting for 68% of the global production [3]. In Africa, Ethiopia, Morocco, and Tunisia are the leading producers, but 61% of the areas and 68% of the production belong to Ethiopia [4]. Lentil is an important dietary source of macro- and micronutrients for both rural and urban dwellers. It provides sufficient amounts of the most essential amino acids to meet nutrient requirements. It is a cash crop fetching the highest price in domestic market compared to all other food legumes and major cereal crops [5]. Global, African, and Ethiopian lentil productivity is about 887 kg/ha, 644 kg/ha [3], and 1168 kg/ha, respectively [6]. Biotic and abiotic factors limit lentil productivity and seed quality. The major yield gap contributing factors are insect pests, diseases, low yielding landraces grown by farmers, and the narrow genetic base [7, 8].

Prior knowledge of genetic variability and characterization of genetic resources within the germplasm available at a gene bank has an important implication for future utilization and collection activities, to identify areas of major priority for conservation and improvement programs [9]. Quantitative traits provide an estimate of genetic diversity, and various numerical taxonomic techniques have been successfully used to classify and measure the pattern of genetic diversity in germplasm, as in lentil [10, 11]. Morphological characterization is the first step in the classification and description of any crop germplasm [12]. The assessment of genetic variability present in a given crop population can be determined by using the biometrical components such as range, variance, coefficient of variation, standard error, and heritability. The relative magnitude of these components determines the genetic properties of the population, in particular, the degree of resemblance between germplasm, that is, various determinants of the phenotype [13]. Reference [14] emphasized that the uses of the genetic coefficient of variability together with heritability and genetic advance values would determine the best picture of the amount of progress to be expected from a selection and determine the selection method to improve a character.

Several researchers investigated the diversity of economically important traits of lentil germplasm and found considerable variations [10, 11, 1520]. Some quantitative evaluations of the Ethiopian landraces revealed potential variability for morphological traits [5, 7, 8, 2125]. Understanding the genetic relationships and diversity of the Ethiopian lentil in relation to the germplasm from other countries is important to widen the genetic base of a germplasm. However, information is not exhaustive for different population structures of local and exotic genotypes on the genetic variability of important traits, which have economic benefit. Hence, the objectives of this study were to generate information on morphological variability and to estimate the phenotypic coefficient of variation, the genetic coefficient of variation, the heritability, and genetic advance of quantitative traits in lentil.

2. Materials and Methods

2.1. Description of the Study Sites

The field experiments were conducted on a hot spot area for rust. The areas included Sirinka Agricultural Research Center (SIARC) in the northeastern part of Ethiopia for two seasons (2011 and 2012), Chefe Donsa in the central part of Ethiopia, and Sinana Agricultural Research Center (SARC) in southeastern part of Ethiopia during the 2011/12 cropping season (Table 1 and Figure 1).

2.2. Plant Materials

Out of the 228 genotypes considered for this study, 158 were planted for morphological evaluation at SRARC in the 2010/2011 cropping season. Of these, 104 genotypes were from the Ethiopian lentil gene pool collected from six major lentil production regions (Figure 2) and were kept by the Institute of Biodiversity Conservation (IBC). The accessions represent over 17% of the 618 lentil accession holdings of IBC. Rust susceptible check/EL-142/and wilt susceptible check/ILL-590/NEL 590/, as well as six commercial national released varieties (Alemaya, Adaa, Teshale, Alem Tena, Chekol and Derash), were used as a check. Ten elite breeding lines from the DebreZeit Agricultural Research Station (DZARC), Ethiopia, were included. Thirty-six parental lines and elite lines based on superiority for agronomic traits introduced from the International Center for Agricultural Research in Dry Areas (ICARDA) were included (Table 2). In the 2011/12 cropping season, 228 genotypes plus RILs were included in the study across three locations: Sirinka, Chefe Donsa, and Sinana.

2.3. Experimental Layout and Design

A randomized complete block design (RCBD) was used with three replications, in the 2010/2011 cropping season at SIARC. A unit plot comprised 2-meter length row with a plot size of (0.8 m2). Row-to row distance was 20 cm. The distance between two plots was 50 cm and the distance between two blocks was 100 cm. The genotypes were planted in the first week of July. In the 2011/12 cropping season, the experiment was carried out using the augmented design with a single row of 2 m length for each genotype planted in five blocks. Eight checks were replicated within each block. Planting was done in the first week of July at SIARC, in mid-August at Chefe Donsa and in mid-September at SARC. The recommended agronomic packages were applied for raising a successful crop.

2.4. Data Collection

Data were recorded on 10 randomly selected plants on plant basis for plant height (PH), number of pods/plant (NP), number of seeds/pod (NS), number of seeds/plant (NSPP), and seed weight/plant (SWPP), whereas days to flowering (DFF), days to maturity (DM), above ground biomass (BI), 100-seed weight in gram (SW), seed yield/plant (SY), pod set score (PSS) (1–9 rating scale where 1 = very poor pod setting and 9 = excellent pod setting), and rust severity scored on a 1–9 scale (1–9 rating scale where 1 = highly resistant and 9 = highly susceptible) before flowering (RSBF) and after flowering (RSAF) were recorded on plot basis according to [26].

2.5. Data Analysis

Plot mean values were calculated for all traits and used for the analysis of variance (ANOVA). The estimation of genetic parameters was analyzed using GenStat Release 13.3. [27]. Phenotypic and genotypic variances for the augmented design were computed for all traits based on the methods of [28]. The pooled analysis for augmented design was carried out according to [29, 30] using GenStat Release 15.1 [31]. The broad sense of heritability ( ) was estimated using the genotype mean according to [14]. Consider

Genetic advance (GA) for selection intensity ( ) at 5% is estimated according to [14] using the formula: .

The genetic advance as the percentage of the mean (GAM) was estimated as per [32] formula: .

3. Results

3.1. Analysis of Variance

Highly significant genotypic ( ) differences were observed for all parameters measured in the 2010/11 cropping season at Sirinka indicating that the variation was genetic. The mean square, CV, and mean values for traits of the genotype were presented in Table 3.

Moderate variabilities with an intermediate coefficient of variation (CV) among the genotypes were observed for days to 50% flowering, days to 90% maturity, plant height, and pod set score, whereas considerable variations were observed for the rest of the traits among the genotypes (Table 3).

The analysis of variance for the 2011/12 cropping season within the check showed significant variation at Sirinka, except for two traits: number of seeds/pods and number of seeds/plant (Table 4). At Chefe Donsa, days to 50% flowering, number of pods/plant, number of seeds/pods, number of seeds/plant, and seed weight/plant had nonsignificant variation within the checks (Table 5). At Sinana, all traits showed significant variability, except days to 50% flowering, days to 90% maturity, number of seeds/pods, and above ground biomass within the checks (Table 6). Consistent significant variations were noted across locations within the checks for plant height, 100-seed weight, and seed yield (Tables 4, 5, and 6).

The analysis of variance showed a highly significant variation ( ) for days to 50% flowering, days to 90% maturity, number of pods/plant, and number of seeds/pods within the test genotypes at Sirinka (Table 4). At Chefe Donsa, analysis depicted a highly significant variation ( ) within test genotypes for days to 50% flowering, plant height, 100-seed weight, and rust disease severity score after flowering (Table 5), whereas at Sinana all traits revealed significant variation ( ) except days to 90% maturity, plant height, above ground biomass, and rust disease severity score after flowering within test genotypes (Table 6).

The analysis of variance between check and test genotypes for the 2011/12 cropping season showed highly significant differences except for days to 50% flowering, days to 90% maturity, number of pods/plant, and plant height at Sirinka (Table 4). At Chefe Donsa, significant ( ) variations were observed for days to 50% flowering, number of pods/plant, seed weight/plant, above ground biomass, seed yield, and rust disease severity score between checks and test genotypes (Table 5), whereas, at Sinana, significant variations were recorded between checks and test genotypes for number of pods/plant, plant height, 100-seed weight, seed yield, and rust severity score before and after flowering (Table 6). The test for homogeneity of error variance was significant across locations because the genotypes perform differently across the three locations for yield and yield component traits. Thus, it could not make further valid inferences for the pooled ANOVA (Table 7). Besides, rust infection was severe at Sinana, mild at Chefe Donsa, and absent at Sirinka and then no combined analysis was made for this character.

The minimum and maximum values of 11 agromorphological traits analyzed for the 2010/2011 cropping season at Sirinka were presented in (Table 8). The range for days to 50% flowering was recorded from 42 to 78 days for Acc. number 36120 and breeding line (FLIP-04-26L), respectively (Table 12). The shortest maturity period was recorded (77 days) for Acc. number 36025 and the longest (138 days) for FLIP-04-26L. Parent X2003S 223 produced the minimum (1.6) number of pods/plant while the maximum (73.3 pods/plant) was noted for breeding line (R-186XFLIP-86-38L-24). The shortest plant height of 20.6 cm was recorded for breeding line (FLIP-97-61L) and the highest was 44.7 cm for Acc. number 36001. The minimum number of seeds/pod was recorded (0.5 seed) for landrace (Acc. number 36098) and the maximum (2 seeds) for parent line (X2003S 233/8009/). Seed weight/plant ranged from 0.1 g for parent line X2002S 219/shehor-74/7554 to 3.4 g for breeding line R-186XFLIP-86-38L-24 (Table 12). The number of seeds/plant ranged from 5.5 for landrace (Acc. number 36098) to 103.1 for lines R-186XFLIP-86-38L-24. The largest seed size was observed for parent line FLIP-2004-7L (4.9 g) while the parent line X2003S 233/8009/had a minimum seed weight of 1.4 g. The lowest biomass/plant was recorded for the parent line X2003S 223 (35 g), and the highest biomass/plant was recorded for the parent line X2006S 128/5480/(767.8 g). The lowest seed yield (0 g) was recorded for landrace X2003S 223 whereas breeding line ILL-590 showed the highest seed yield (347.2 g) (Table 12).

3.2. Phenotypic and Genotypic Coefficients of Variation

Though variability in population is an indispensable prerequisite for any improvement, it is not the only criterion for deciding as to which trait is showing the highest degree of variability. Phenotype coefficient of variation (PCV) and genotype coefficient of variation (GCV) can help in this regard. Estimates of genetic parameter are shown in Tables 8, 9, 10, and 11. PCV and GCV values of approximately more than 20% are regarded as high, whereas values less than 10% are considered low and values in between are considered as medium [33]. Based on this delineation, high PCV and GCV values were recorded for number of pods/plant, number of seeds/pod, pod weight/plant, seed weight/plant, number of seeds/plant, 100 seed weight, above ground biomass, and seed yield. However, days to 50% flowering, days to 90% maturity, and plant height showed lower GCV and PCV values in the 2010/2011 cropping season. In the 2011/12 cropping season, similar GCV and PCV patterns were recorded across the three locations for number of pods/plant, seed weight/plant, number of seeds/plant, above ground biomass, seed yield, and rust disease score, whereas low PCV and GCV were noted for number of pods/plant, number of seeds/pods, seed weight/plant, and number of seeds/plant at Chefe Donsa.

The lowest PCV (10%) was recorded for days to 50% flowering and the highest (60.7%) for pod weight/plant at Sirinka 2010/2011. At Sirinka 2010/2011, GCV values ranged from 7.6% for days to 50% flowering to 58.4% for pod weight/plant (Table 8). The maximum PCV and GCV recorded for number of pods/plant were 89.8 and 75, respectively, at Sirinka 2012, while the lowest PCV and GCV recorded for plant height were 14.5 and 8.6, respectively (Table 9). Days to 90% maturity and 100-seed weight showed moderate GCV. At Chefe Donsa and Sinana, values recorded for both PCV and GCV were similar with those of Sirinka (Tables 10 and 11). Traits such as number of pods/plant, seed weight/plant, 100-seed weight, above ground biomass, seed yield, and rust disease score showed relatively high PCV and GCV values across the four environments (Tables 811).

3.3. Heritability and Genetic Advance

The broad sense of heritability ( ) estimates of the traits ranged from 0.4% (number of seeds/plant) to 92.6% (pod weight/plant). According to [34], heritability estimate in cultivated plants can be placed in the following categories: heritability estimate ranging from 5 to 10% low, values ranging from 10 to 30% medium, and values 30% and above as high heritability. Based on the above classification, most of the characters have shown high heritability over locations except at Chefe Donsa. At Chefe Donsa, low heritability estimates were recorded for number of seeds/plant (0.4%), number of seeds/pods (0.7%), and number of pods/plant (2.5%) (Table 10). At Sinana, heritability estimates ranged from 18% to 88% for biomass and 100-seed weight, respectively. Days to 90% maturity, number of seeds/plant, plant height, and above ground biomass had a moderate heritability (17.7–26.1%) (Table 11), whereas seed yield, 100-seed weight, rust disease score, number of pods/plant, and number of seeds/plant scored higher heritability.

Consistently higher value of heritability was revealed for days to 50% flowering, 100-seed weight, days to 90% maturity, and plant height over locations. However, all these four characters associated with low genetic advance value (GA) because of low GCV, whereas 100-seed weight, rust severity score, and seed yield showed high heritability values coupled with high genetic advance values across locations because of high of GCV. Numbers of seeds/plant, seed weight/plant, and number of seeds/pod have intermediate GCV and heritability with a mild response to selection.

4. Discussion

Morphological diversity has been used to characterize germplasm from a range of plant species and allows the assessment of genetic diversity. The purpose of this study was to compare the genetic diversity within and among Ethiopian and exotic lentil genotypes of the Ethiopian possessions. In this paper, the characterizations of morphological variation of lentil germplasm using biometrical model was described.

As per ANOVA, highly significant phenotypic ( ) differences were observed for all traits measured indicating the existence of variation. The presence of highly significant variation among the genotypes for the morphological traits was indicative of the presence of high degree of genetic variation for future breeding programs through selection. Similarly, from follow-up study with the simple sequence repeat (SSR), genetic diversity analysis was a supportive evidence of the presence of a considerable variability among both local and exotic lentil genotypes of the Ethiopia germplasm (unpublished data). In agreement with this report, [21, 3439] reported significant genetic variations among lentil landraces. On the other hand, some studies showed that no significant variation was recorded for some traits, like seed weight/plant and days to flowering [21, 38, 40].

The results recorded for coefficient of variation were comparable with [38] for days to 50% flowering (8.8), plant height (9.1), number of seeds/plant (43.2), and above ground biomass (31.5). However, [38] reported higher values vis-a-vis with this study for days to 90% maturity (19.2), pod weight/plant (47), and seed weight/plant (47.8). In agreement with this finding, [41] also noted comparable range for plant height on 3974 accessions of lentil. However, they recorded a wider range of values for 100-seed weight with an overall mean than with our records on lentil accessions. In another study, [37] reported a narrow range for yield component traits among 46 genotypes from South East Anatolia region of Turkey as opposed to our findings.

Significant variations were noted over locations among the test genotypes and within the checks with a few exceptions. In line with this study, [24] recorded the wide range of genetic variability among yield related characters in exotic lentil lines in Ethiopia. Comparable to our study, [15, 18] noted the presence of variation in characters, such as grain and straw yield, 100-seed weight, days to 50% flowering, days to 90% maturity, number of seeds/pods, plant height, pod number/peduncle, and resistance to various biotic and a biotic stresses in the lentil. Besides, the variation among environments was significant for all morphological traits in view of diverse agro-climatic features of the test sites. Similarly, [42] also observed the greater portion of total variance was due to environment. Reference [5, 7, 23] also reported consistent regional differences among lentil landraces of Ethiopia for time to flowering and maturity, 100-seed weight, seed/pod, and plant height.

ANOVA over locations showed consistent significant variations were recorded within the checks for 100-seed weight, above ground biomass, seed yield, and rust disease score. These implied that yield increment was achieved for the last three decades, largely for 100-seed weight, above ground biomass, and seed yield of lentil genotypes. A consistent significant variation was also recorded within test genotypes over three locations for days to 50% flowering and 100-seed weight. Besides, reliable variation was noted over all locations between checks and test genotypes for above ground biomass and seed yield. Similarly, according to this (unpublished data) follow-up study, using the multivariate analysis reported that seed yield, above ground biomass, seed weight/plant, 100-seed weight, rust severity score, and plant height were the major agromorphological traits contributing to variations among the genotypes.

Estimate of PCV and GVC for days to 50% flowering, days to 90% maturity, and 100-seed weight showed narrow differences, indicating a relatively low influence of the environment on these characters. Similar observation of considerable closeness between GCV and PCV for these characters made by [43] indicated low influence of environmental factors in lentil. In general, PCV values for most of the traits at Sirinka and Sinana were higher than that record noted at Chefe Donsa. In addition, the PCV was generally higher than the GCV for all characters. Consistently higher PCV values were observed across four environments for rust disease scores, seed yield, pod weight/plant, seed weight/plant, number of pods/plant, 100-seed weight, above ground plant biomass, and number of seeds/plant. Other studies showed similar high PCV and GCV for seed yield, above ground biomass, and 100-seed weight [44, 45]. In contrast, days to 50% flowering, days to 90% maturity, and plant height showed comparatively low PCV values (less than 20%) for all locations. It was observed that GCV for number of pods/plant, seed yield, rust disease score, 100-seed weight, number of seeds/plant, and biological yield/plant was high in both seasons. However, traits such as days to 50% flowering, days to 90% maturity, and plant height consistently showed lower GCV values of between 2.4% and 17%. The remaining traits lack consistency and showed low to moderate GCV for number of seeds/pods, seed weight/plant, and pod weight/plant. Similar significant and considerable variation for (phenotypic and genotypic) results was reported for biological yield/plant, seed yield/plant, and 100-seed weight by [15, 35, 43, 4548]. However, as opposed to this finding, [49, 50] reported higher GCV and PCV for days to 50% flowering but they reported the same record for days to 90% maturity and plant height. Similarly, [35] reported the minimum PCV and GCV for days to maturity and number of seeds/pod.

According to [51], most characters have showed high heritability estimate over locations except at Chefe Donsa due to severe wilt incidence. Reference [52] reported comparable high heritability estimates for all traits except for plant height and biological yield on lentil. Reference [15] similarly recorded high heritability for 100-seed weight, days to 50% flowering, and days to 90% maturity. Reference [7] also reported high heritability for seed weight and days to 50% flowering. In line with our finding, [53] reported that those days to 50% flowering showed higher heritability estimates; however, they differed in reporting with higher genetic advance. 100-seed weight, rust severity score, and seed yield showed high heritability values coupled with high genetic advance values across locations because of high GCV. Numbers of seeds/plant, seed weight/plant, and number of seeds/pods have intermediate GCV and heritability with a mild response to selection. Similar to our findings, [45] reported higher degree of PCV, GCV, and genetic advance for seed yield/plant on lentil. However, in contrast to these reports, higher PCV, GCV, and genetic advance were recorded for biological yield/plant.

The GCV along with heritability estimates provide a reliable estimate of the amount of genetic advance expected through phenotypic selection [54]. Based on the underlying facts, four groups of traits were categorized as per the analysis: the first group included plant height, days to 50% flowering, and days to 90% maturity which have low GCV with a high heritability but they end up with low genetic advance. [42, 44] made a similar report that high heritability estimate for days to 50% flowering and days to 90% maturity with low genetic advance values on lentil. The second group of characters was above ground biomass, number of seeds/plant, and number of pods/plant with intermediate GCV, heritability, and genetic advance. However, deviated from our findings, studies by [42, 44, 45] reported high heritability and expected genetic gain estimate for above ground biomass and number of pods/plant.

The third group of characters, rust severity score, seed yield, and 100-seed weight, has higher GCV with high heritability coupled with high genetic advance as percentage of the mean. Reference [42, 44] reported high heritability estimate for seed yield and 100-seed weight coupled with higher genetic gain. High heritability for the characters indicated that these traits were less affected by environmental factors. This is indicative of the fact that these traits are mostly controlled by genetic factors and expected to respond to direct selection for traits improvement. The fourth group of characters, pod weight/plant, seed weight/plant, and number of seeds/pods, showed inconsistent GCV, heritability, and genetic advance values across locations. These findings indicated that measuring yield components in breeding program would be difficult. From our follow-up study by (unpublished data), the second and the third groups of characters responded more to selection and the most important traits that contributed to the genetic divergence of lentil. From our follow-up association study (unpublished data), these traits attributed a strong association with seed yield. Therefore, for selecting high yielding lentil genotypes, the breeder should give emphasis to higher seed weight/plant, more number of pods/plant, and plants with short to intermediate height. This observation is in accordance with [55, 56] selection based on number of pods/plant and seed weight/plant which were the most important characters that contributed to seed yield. However, days to 90% maturity, rust disease severity score, and days to 50% flowering which otherwise had a negative phenotype correlation with seed yield had substantially negative direct effect on seed yield (unpublished data).

5. Conclusion

This study has described a high phenotypic diversity for important agromorphological characteristics of landraces and exotic genotypes of Ethiopia gene pool. The incidence of highly significant variation between landraces and exotic line for the majority of traits considered is a sign of the presence of a high degree of genetic variation implying the great potential of the Ethiopian germplasm in future breeding programs through selection. As per the analysis, four groups of traits were categorized. Some promising landraces were identified with superior plant characteristics, such as early flowering and early maturing accessions, wide range of genetic base for plant height, number of seeds/plant, resistance to rust, and both low and high yielding lines. Seed yield, 100-seed weight in gram, rust disease score, and total number of seeds/plant consistently showed higher GCV, heritability, and genetic advance in the lentil genotypes. All these important agromorphological characteristic traits with their associated genotypes could be successfully utilized in breeding programs that are aimed at improving the yield and the yield components of lentil. The wide genetic base that is possessed by the Ethiopian germplasm would be integrated in a breeding program and used as a parent in a recombination breeding.

In conclusion, this study has demonstrated that the exotic accessions introduced from ICARDA have a rich genetic potential for 100-seed weight, number of seeds/plant, seed weight/plant, short stature, and resistant source for rust. The Ethiopian landraces were potentials for earliness, higher number of pods/plant, and long plant height. This study has showed a high phenotypic diversity of important plant characteristics of the Ethiopian lentil germplasm. However, additional search is needed to develop a specific plant idotype for major agroecological zone to provide as a platform for local adaptation and then breed for location relevant traits, such as a biotic stresses and identification source of resistance for other major important diseases.

Conflict of Interests

The authors declare that they have no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors are most grateful for Gonder University, Haramaya University, Federal Ministry of Education, ICARDA Legume Research Program, ICARDA Ethiopia country office, Ethiopian Institute of Agricultural Research Center, National Chickpea and Lentil Research Project in Debrezeit Agricultural Research Center, Sirinka Agricultural Research Center and Sinana Agricultural Research Center, and Institute Biodiversity Conservation, for research funding supports of this study.