Heat and Drought Stress Impact on Phenology, Grain Yield, and Nutritional Quality of Lentil (Lens culinaris Medikus)

Lentil (Lens culinaris Medikus) is a protein-rich cool-season food legume with an excellent source of protein, prebiotic carbohydrates, minerals, and vitamins. With climate change, heat, and drought stresses have become more frequent and intense in lentil growing areas with a strong influence on phenology, grain yield, and nutritional quality. This study aimed to assess the impact of heat and drought stresses on phenology, grain yield, and nutritional quality of lentil. For this purpose, 100 lentil genotypes from the global collection were evaluated under normal, heat, and combined heat-drought conditions. Analysis of variance revealed significant differences (p < 0.001) among lentil genotypes for phenological traits, yield components, and grain quality traits. Under no stress conditions, mineral concentrations among lentil genotypes varied from 48 to 109 mg kg−1 for iron (Fe) and from 31 to 65 mg kg−1 for zinc (Zn), while crude protein content ranged from 22.5 to 32.0%. Iron, zinc, and crude protein content were significantly reduced under stress conditions, and the effect of combined heat-drought stress was more severe than heat stress alone. A significant positive correlation was observed between iron and zinc concentrations under both no stress and stress conditions. Based on grain yield, crude protein, and iron and zinc concentrations, lentil genotypes were grouped into three clusters following the hierarchical cluster analysis. Promising lentil genotypes with high micronutrient contents, crude protein, and grain yield with the least effect of heat and drought stress were identified as the potential donors for biofortification in the lentil breeding program.


INTRODUCTION
Micronutrient malnutrition, hidden hunger, is a serious health problem affecting more than two billion people globally (1). Among the 17 essential minerals required by the human body, iron and zinc serve a multitude of biological functions (2). An inadequate bioavailability of these two elements, in commonly eaten foods, leads to serious physiological disorders, resulting in a range of health hazards including stunting and learning disabilities in children, anemia, increased morbidity, and mortality rate (3). In Sub-Saharan Africa and South Asia regions, populations living below the poverty line are at high risk of iron (Fe) and zinc (Zn) deficiencies mainly due to the predominance in their everyday diets, of cereals-based foods which are inherently low in essential micronutrients. Iron deficiency is estimated to be prevalent among 60% of the world population (4,5), whereas zinc deficiency is widespread among 31% population (6). There are effective interventions, including diet diversification, food fortification, and nutrients supplementation, to combat micronutrient malnutrition. However, these interventions are beyond the reach of poor households due to poverty, poor distribution network, and limited access to commercially marketed fortified foods (7,8). As a result, sustainable food approaches are essential to mitigate micronutrient malnutrition among developing countries.
Biofortification through conventional plant breeding techniques emerges as the most promising and sustainable strategy to enhance both the level of micronutrient and their bioavailability in the edible part of food crops (9). Therefore, nutritionally improved varieties hold a great promise for meeting the challenge of addressing malnutrition and attaining nutritional security worldwide. However, achieving food security will pose a serious challenge due to the rising global population, from 7.7 billion currently to 9.7 billion in 2050 (10). In addition, climate change and variability have started showing its consequences on agricultural productivity and the nutritional quantity of foods, especially in rainfed dry regions.
Lentil (Lens culinaris Medikus) is one of the most important cool season food legume crops in the world with its cultivation in as many as 52 countries. Globally, it is grown on 6.1 million hectares, with an annual production of 6.3 million tons and productivity of 1,038 kg ha −1 . The major lentil producing countries are Canada, India, USA, Turkey, Australia, Kazakhstan, Nepal, Russian Federation, Bangladesh, China, and Ethiopia, contributing more than 93% to the global output (11). Lentil is known for its high nutritive value. It contains a relatively high amount of protein (20-36%) and provides an affordable source of micronutrients in highly bioavailable forms. It also provides a range of important dietary compounds like carbohydrates, prebiotics, fiber, vitamins, amino acids, and antioxidants (12). Due to its high nutrient density, lentil has emerged as an excellent candidate for micronutrient biofortification (13). Several studies have shown significant genetic variation for nutritional parameters including iron, zinc, and protein content in lentil genotypes, including cultivated and wild species (14)(15)(16). Under the HarvestPlus program, ICARDA, in collaboration with NARS partners, has developed and released many biofortified lentil varieties for cultivation in South Asia, Sub-Saharan Africa, West Asia, and North Africa (https://www.icarda.org/media/ drywire/lentil-biofortification-research-fight-hidden-hunger). With climate change and variability, lentil crop encounters water and heat stresses at different growth stages, causing severe economic damage (17). The severity of damage caused by these stresses depends on timing and the intensity of the stress (18). Lentil is extremely sensitive to heat stress, especially during the flowering and seed filling stages. Temperature exceeding 32 • C can restrict photosynthesis, metabolic pathways, electron flow, and respiration rate (18), which cause flower abortion, pollen infertility, and reduction in the number of pods in lentil (19,20). These changes lead to significant losses in grain yield (20,21). Water stress also affects plants at different growth stages, including vegetative (intermittent drought) and reproductive (terminal drought) stages. Terminal drought can suppress nearly all the processes of lentil growth and metabolism, causing heavy yield losses (19), as it reduces flower production, pod number, and seed number (22). However, very little information is available about the negative effect of heat and drought, during seed filling, on the nutritional quality of lentil. A recent study revealed that heat and drought stress adversely impacted the nutritional quality of lentil grains (23). Therefore, this study was undertaken to (a) investigate the effect of heat and drought stress on traits associated with phenology, grain yield, and nutritional quality under field conditions, and (b) assess the genetic variation for agro-morphological and quality traits in lentil genotypes.

Description of the Study Area
The present study was conducted at ICARDA's Marchouch experimental station (33.56 • N, 6.63 • W, 392 m altitude) in Morocco. This station represents the Mediterranean semiarid environment with high precipitation during November and December. The mean annual precipitation is around 400 mm, and the mean annual temperature is 18 • C. The soil is characterized by the dominance of Vertisols with a clayey texture. The pH of the soil is neutral and ranges from 6.3 to 7.9.

Plant Material and Experimental Design
A set of 100 lentil genotypes comprising 93 germplasm accessions originating from 46 countries and seven elite lines imanating from ICARDA breeding program were evaluated in an alpha lattice design with two replications during 2016-17 cropping season. Each genotype was grown in a 2-row plot of 1 m length, with a spacing of 30 cm between rows. In each row, seeds were sown by hand at a 2-cm depth maintaining 10-cm space between plants. In total, three experiments involving the same set of genotypes were conducted side by side by manipulating the planting date and water supply, in order to impose heat and water stress at the reproductive phase. These three experiments were considered to represent three treatments, namely the normal date of planting (treatment A), late planting with irrigation at field capacity throughout the crop period (treatment B), and late planting without irrigation during the reproductive phase (treatment C), which were referred to as treatments. Treatment A resulted in optimal growing conditions (∼150 mm welldistributed rainfall and below 27 • C temperature) without any heat and water stress to the plants. Treatment B (planted 65 days after normal planting date with irrigation at field capacity throughout the crop duration) imposed heat stress, as the plants were frequently exposed under field conditions to a temperature above 35 • C during the reproductive phase.
In contrast, regular irrigation at field capacity avoided any water stress to the plants. Treatment C (planted 65 days after normal planting date without irrigation during the reproductive phase) imposed a combined heat and water stress. Treatment A was planted on 27 December 2016, and no irrigation was applied during the crop period as the crop received well-distributed enough rainfall at Marchouch (149.8 mm). Treatments B and C were planted on 1 March 2017. Irrigation was applied to maintain water supply at field capacity using sprinkler system throughout the crop duration in treatment B, whereas irrigation was stopped from the flowering initiation stage onward in treatment C, to impose water stress (<5 mm rainfall during the reproductive stage) in addition to the heat stress. All recommended agronomic packages of practices were followed to raise a successful crop, including nutrient and weed management throughout the growing season.

Weather Data
Daily minimum and maximum temperatures ( • C), precipitation (mm), and relative humidity (%) were recorded throughout the crop season at the experimental site. Treatment-A experienced 9.2-23.3 • C maximum and −2.4-10.7 • C minimum temperatures during the vegetative stage and 11.4-30.6 • C and 0.3-11.7 • C during the reproductive stage (Figure 1). For Treatments B and C, maximum temperatures ranged from 13.3 to 30.6 • C and minimum temperatures from 0.3 to 10.5 • C during the vegetative stage and from 22.4 to 43 • C and 3.5 to 18.1 • C during the reproductive stage. The maximum temperatures during the reproductive stage reached the threshold level of 35 • C. Therefore, late planting with irrigation at field capacity was successful in imposing heat stress during the reproductive stage of the test genotypes in treatment B, and late planting without irrigation in imposing combined stress of drought and heat in treatment C.

Observations on Phenological and Yield Traits
Phenological observations were recorded on days to 50% flowering, days to the first pod, and days to physiological maturity on whole plot basis. At the time of harvesting, five plants were randomly selected from each plot for recording observations on yield traits namely, plant height (cm), numbers of primary, secondary and tertiary branches, numbers of total pods, filled pods, and unfilled pods, 100-seed weight (g), biological yield (g) and grain yield (g). The mean value of five plants for each trait was used for statistical analysis.

Determination of Grain Iron and Zinc Concentrations
Seeds of each sample were ground by using a Cyclone mill (Twister, 10 mm−250 µm, Retsch). Mineral (Fe and Zn) concentrations were determined following a modified diacid protocol (24,25). In the digestion block (QBlock series, Horiba), 0.5 g of each ground sample was placed in individual tubes and digested with 6 mL of HNO 3 at 90 • C for 60 min, followed by adding 3 mL of 30% hydrogen peroxide (H 2 O 2 ) to each tube with another period of digestion of 15 min at 90 • C. Then, 3 mL of 6 M hydrochloric acid (HCl) was added to each digestion tube. Finally, after the sample solutions were cooled down, the volume was adjusted to 10 mL, and then filtered. The Fe and Zn concentrations were estimated using inductively coupled plasma-optical emission spectroscopy (ICP-OES); (iCAP-7000 Duo, Thermo Fisher Scientific) at the Cereal and Legume Quality Laboratory, ICARDA, Morocco. Calibration curves for Fe and Zn were made using serial dilution from 0.1 to 10 mg L −1 .

Crude Protein Analysis
The nitrogen content of the samples was determined using the Kjeldahl procedure (26). In this method, the weighted seed samples of each genotype were digested with sulfuric acid followed by alkalizing the solution by adding of sodium hydroxide, and then determining the resulting ammonia by distillation into a measured volume of standard acid, the excess of which is determined by titration. Seed protein contents then were calculated using nitrogen values multiply by 6.25. For this, triplicate analyses were carried out on each sample.

Statistical Analysis
Analysis of variance was performed using the General Linear Model (GLM). Tukey's post-hoc test was applied to compare the differences between the mean values. Associations among traits were assessed using the Pearson correlation coefficient (r). The principal component analysis was conducted to determine which traits explain most of the variation under normal, water stress and heat stress treatments and also to determine different groups of genotypes associated with major traits. Highly correlated traits with grain yield were discarded from the analysis to avoid redundant results.
Hierarchical Cluster Analysis (HCA) using Ward's method based on Euclidian distance was performed to group lentil genotypes based on mineral composition, crude protein, and grain yield. Euclidean distances were calculated after the standardization of the data. The result of the hierarchical clustering procedure was displayed graphically using a tree diagram (Dendrogram). All statistical analysis were conducted using IBM SPSS statistics version 23.

RESULTS
The analysis of variance revealed significant (p < 0.001) variation among lentil genotypes for all the traits under stressed and non-stressed conditions ( Table 1). The treatment effect was also significant for the traits under investigation. The highly significant interaction effect between genotypes and treatments was detected for day to flowering, podding and maturity, grain yield, hundred-seed weight, and a number of filled and unfilled pods per plant. As for seed quality traits, the genotype × treatment interaction effect was highly significant for crude protein, iron, and zinc concentrations.

Effect Heat and Drought on Phenology and Morphological Traits
Based on each agro-morphological trait, the response of genotypes at each condition differed. The mean value of time to flower was 69 days in normal sowing while in late conditions,  Genotype ** ** ** ** ** ** ** ** ** ** ** ** ** ** ** ** ** Treatment ** ** ** ** ** ** ** ** ** ** ** ** ** ** ** ** **   Table 2). Similarly, days to the first pod declined by 33.4% in heat stress and 33.7% in combined heat and drought stresses. Day to maturity ranged from 100 to 120 days in normal planting, which was declined significantly in late planting with more reduction in combined heat and drought stress (23.1%) compared with heat stress (22.9%). The mean plant height was 31 cm, with a range from 21.6 to 38.3 cm in normal planting. In late sowing, plant height ranged from 16.2 to 31.2 cm in heat stress and 13.4 to 29.4 cm in heatdrought stress. The mean was 24.9 and 22.1 cm, respectively. The number of primary, secondary, and tertiary branches decreased markedly in late planting and found to decrease more in heatdrought stress. The mean of primary branches was 2.4 in normal planting. In late planting, the mean was 2.3 under heat stress condition and 2.1 under heat-drought stress condition. The secondary branches showed a marked reduction in late planting; the reduction was more significant in heat-drought stress (41.3%) than heat stress (38.1%). Similar results were found for tertiary branches, which was reduced by 53.1% in heat and by 56.8% in the combined stress.

Effect of Heat and Drought on Grain Yield and Yield Traits
Under normal planting conditions, the mean value of biomass was 11.1 g per plant and ranged from 5.8 to 20.4 g ( Table 3).
On the other hand, genotypes under late planting conditions produced less biomass, which was reduced by 42.2% in heat stress and by 46.3% in combined heat-drought stress. Compared with no stress condition, the number of filled pods reduced by 55.6% under heat-drought stress and 50.7% under heat stress. Similarly,

Association Between Seed Yield and Other Agro-Morphological Traits
The coefficients of correlation between seed yield and other studied traits under normal and late conditions are shown in Tables 5, 6 Under stress conditions, grain yield showed a significant positive correlation with plant height, biomass, and harvest index in both heat and combined heat-drought stresses ( Table 6). Grain yield was positively correlated with secondary and tertiary branches in both treatments (p < 0.01). A significant positive correlation was obtained between grain yield and the number of filled pods per plant in both    heat and combined heat-drought stresses; the coefficient of correlation was r = 0.87 and 0.82, respectively. Similar to the normal planting, the correlation coefficient of grain yield and 100-seed yield was highly significant and positive under both heat and heat-drought treatment (r = 0.65 and 0.49; p < 0.01).
In both normal and late planting conditions, seed yield showed no significant correlation between seed mineral contents and crude protein contents. A significant positive correlation was demonstrated between iron and zinc concentrations in normal planting. The correlation was relatively low (r = 0.28; p < 0.01). Similarly, a significant positive correlation was found between iron and zinc concentrations under both heat and combined heat-drought treatment (r = 0.26; p < 0.05 and r = 0.36; p < 0.01, respectively). On the other hand, the correlation coefficient between zinc and crude protein was significant and negative (r = −0.26; p < 0.05) in normal condition (Figure 3).

Principal Component Analysis
Principal component analysis (PCA) was used to determine which traits would differentiate the genotypes under three different treatments. Eigenvalues of the studied traits and cumulative variance are reported for each treatment in Table 7. Under stress conditions, the first two principal components explained 88% of the total variability. PCA1 and PCA2 accounted for 74.4 and 13.6% of the total variation. The PCA1 showed high association with filled pods per plant and low association with harvest index, while PCA2 was highly associated with iron content. Under heat stress treatment, the first two axes of the PCA explained 79.8% of the total variation. PCA1 accounted for 62.2% of the total variation, and PCA2 explained 17.6% of the total variation. The PCA1 was mainly associated with filled pods per plant, while PCA2 was negatively associated with unfilled pods per plant and positively with HI. Under the combined heat-drought conditions, the first two axes of the PCA explained 84.3% of the total variation with PCA1 and PCA 2, explaining 73% and 11.3% of the total variation. The PCA1 was mainly associated with filled pods per plant, while PCA2 was negatively associated with unfilled pods (−0.31) and iron content (−0.33) and positively associated with harvest index (0.88). These results revealed that the iron content contributed significantly to the variation among genotypes under normal conditions while its contribution was limited under heat and combined droughtheat stress. Biplot of PCA (Figure 4) showed that the evaluated genotypes were clustered according to iron content and number of filled pods. Under no stress conditions, the PCA1 showed clear differences among genotypes based on the number of filled pods. The right hand of this axis was characterized by genotypes with a high number of filled pods, which was highly correlated with grain yield. The PCA2 showed clear differences among genotypes based on iron content. In the upper part of the graph, all lines had higher iron concentration than the average. The upper and right part of the biplot, 17 genotypes, namely ILL 82, ILL 890, ILL 6870, ILL 6848, ILL 6281 ILL 6272, ILL 5553, ILL 5505, ILL 5416, ILL 5384, ILL 5151, ILL 4830, ILL 4606, ILL 4164, ILL 2297, ILL 1861, and ILL 224 contained higher iron and showed a higher number of filled pods compared to tested genotypes. Given the fact that iron content was highly correlated with zinc and number of filled pods with yield, the genotypes mentioned above should also be high yielding and high zinc concentration.

Identification of Potential Genotypes for Biofortification
Given the cluster analysis of genotypes based on grain yield, crude protein, and nutrient content, three distinct clusters were generated. The mean values of each trait in each of the three  clusters are presented in Table 8. Under normal conditions, the first cluster grouped 34 lentil genotypes; this group is characterized by a high value of Fe and Zn concentrations and a moderate value of crude protein and grain yield. The second cluster contained 34 genotypes, which are characterized by a high value of iron concentration, crude protein and grain yield, and a moderate level of Zn concentration. On the other hand, the third cluster had 32 lentil genotypes. Most of the genotypes remaining in this group showed a high level of protein and a moderate level of zinc, iron, and grain yield (Figure 5).  Under the late condition, cluster analysis revealed the existence of three groups for both heat and combined heat and drought stresses. For heat stress, genotypes in the first cluster had a high level of Fe, moderate level of Zn concentration and grain yield, and low crude protein level. The second cluster was characterized by a moderate level of Zn and protein content despite reduced grain yield. Whereas, genotypes belonging to the third cluster showed high grain yield and a moderate level of Fe and crude protein. Likewise, under combined stress, genotypes in the first cluster showed moderate grain yield and Fe concentration and low crude protein. While genotype existing in the second cluster showed poor yield and low crude protein but moderate Fe and Zn concentrations. Genotypes existing in the third cluster exhibited low grain yield, crude protein, Fe, and Zn concentrations.
The top 10 lentil genotypes with a combination of high yield, high crude protein, and high micronutrient contents under normal conditions are listed in Table 9. Among the 10 top genotypes, ILL 6281, ILL 223, ILL 5151, and LSI88 are characterized by the high level of iron concentration. In comparison, the highest amount of zinc concentration was recorded for LSI88, ILL 6281, and ILL 5505. The highest crude protein was obtained for LSI88, followed by ILL 7084. As for grain yield, the greatest amount was observed for ILL 5505.

DISCUSSION
In this study, the late planting technique was used to expose flowering and pod filling period to high temperatures. However, the confounding effects of some other environmental variables, such as relative humidity, cannot be completely controlled. The fact that heat stress occurs together with drought stress due to rapid water loss from the soil and plants (27), makes it difficult to investigate their unique impact under field conditions. Thus, to maintain heat stress conditions, high-temperature treatment was accompanied by frequent irrigation.
Lentil is particularly sensitive to high temperature at the time of flowering and seed-filling period. It may also face the combined effects of heat and drought stress, which are more severe than the individual effect (28). The temperature for optimal growth and development of lentil ranges from 18 to 30 • C (29), whereas in the current study, plants were exposed to a temperature of 43 • C during the reproductive stage. Moreover, exposure to high daytime temperatures (30-35 • C) at this stage can have a damaging effect on many processes, including photosynthesis, metabolic pathways, electron flow, and respiration rate (19).
Our results indicate that heat stress and combined heatdrought stress, during flowering and pods filling stage, had severe impacts on growth, seed yield components, and, especially, seed nutritional quality in studied lentil genotypes. Heat stress and combined heat-drought stress accelerated all phenological growth stages. Therefore, both stresses reduced the number of days to flowering, number of days to podding, and number of days to maturity, resulting in rapid transition from vegetative to reproductive stages. Compared to the normal sowing, flowering and maturity duration were drastically reduced under both stresses, which is in accordance with earlier findings in lentil (28), chickpea (30), and common bean (31). Plant height and total plant biomass were significantly reduced by heat stress and combined stress; similar reduction has been observed in chickpea (30), faba bean (32), and common bean (33). This reduction is related to the inhibited expression of growth-related metabolism (34). A significant reduction in the number of branches per plant was noticed under both stresses. This trait was positively correlated with grain yield under all treatments, which agrees with previous observation in lentil (35)(36)(37)(38). Similar results were found under heat and combined heat-drought stress. Moreover, it has been shown that branch production is more important in achieving higher seed yield in lentil (39), which indicates the direct effect of the number of branches on seed yield.
Heat and combined heat-drought stress considerably decreased hundred-seed weight, which was consistent with earlier findings in lentil (23) and chickpea (30). Seed weight was a positive correlation with grain yield under both normal and late conditions; a similar kind of association was observed previously in lentil (40). Harvest index was significantly reduced under both stresses, which is consistent with the result found by Bourgault et al. (41), who reported a 16% decrease in harvest index of lentil due to the heatwave. Harvest index was also found to decrease under drought (42,43). This index was a positive correlation with grain yield under stress and no stress condition which was in accordance with previous studies on lentil (44,45), chickpea (46,47). Moreover, Berny Mier Y Teran et al. (48) recently reported the existing genetic relationship between harvest index and yield in common bean and confirm the role of harvest index in the selection of both additive and epistatic effects controlling drought tolerance.
A considerable reduction in the total number of filled pods in response to heat stress and combined heat-drought stress is consistent with previous studies, where a reduction in a number of filled pods of lentil was caused by increased flower and seed abortion (49). Similar results were found in chickpea, where water deficit during pod filling period increased pod abortion (50). Several reports showed a reduction in pod set under a high-temperature regime in chickpea, which was attributed to a reduction in pollen viability and germination (51,52).
Heat stress and heat-drought stress significantly decreased seed yield and seed weight due to the limited number of filled pods and pods losses. Reductions in grain yield varied among tested genotypes, indicating their different responses to drought and heat stress conditions. Grain yield and seed weight declined more under combined heat-drought stress than heat stress only. Our findings are supported by those found earlier (28,53), where the impact of combined heat-drought on grain yield and grain weight was attributed to physiological and metabolic impairment of the photosynthetic components and water relations. Similar results suggested that heat alone and, in combination with drought stress, severely impaired photosynthetic function in chickpea (30).
Several studies reported that seed compound concentrations in legumes varied in response to genetic and environmental factors (54). Seed nutritional quality in lentil, including protein and mineral contents, is also significantly influenced by environmental factors and might be very sensitive to extreme weather conditions. Our results showed that iron and zinc concentrations significantly decreased under heat stress conditions, this might be attributed to the reduction of root nutrient uptake, by decreasing root biomass and metabolic rate (55) or by direct damage to roots (56). Combined heatdrought stress further decreased iron and zinc level in lentil seeds, which might be explained by the fact that decreasing water availability under drought conditions results in low micronutrient uptake, including iron and zinc (57). Furthermore, a reduced transpiration rate due to water deficit may also reduce nutrient absorption and efficiency of their utilization by the plant (58). Supporting results were recently published (23), where combined stresses intensified the impact on seed nutrients, by decreasing leaf water relations to inhibit translocation into the developing seeds. On the other hand, Smith et al. (59), reported that drought stress under field conditions had no negative impact on the accumulation of Fe and Zn in common been seeds.
Heat stress-induced significant reduction in the amount of nitrogen and, consequently, in protein content of lentil seeds. Whereas, crude protein content was dramatically reduced in response to combined heat-drought stress. Several researchers reported that drought stress is associated with changes in physiological and biochemical processes, including inhibition of protein synthesis (54). Furthermore, drought stress reduces nitrogen partitioning and fixation, resulting in a reduction in the rate of protein accumulation in the seeds (60). Our findings are consistent with early findings in lentil (23), chickpea (50), and bean (60)(61)(62), where drought stress reduced seed nitrogen and consequently seed protein content. In contrast to our findings, Hummel et al. (63) revealed that the protein level increased under drought stress conditions.
The correlation coefficient between seed mineral contents under both normal and late conditions indicates that Fe is positively correlated with Zn. In the same way, many researchers reported similar observation for Fe and Zn association in lentil seeds (12,14,64), chickpea (65), bean (66), and faba bean (67). In contrast, Kumar et al. (15), Thavarajah et al. (68), and Darai et al. (69) revealed no correlation between Fe and Zn contents. The positive correlations among Fe and Zn contents suggest common uptake pathways or transporters of these minerals (70,71). These results imply that high Fe concentration can be accompanied by high Zn concentration. Under normal condition, crude protein content exhibit a low but positive correlation (r = 0.03) with iron content and a significant negative correlation (r = −0.26; p < 0.05) with zinc content. Our results were consistent with earlier findings in lentil (72) and chickpea (73). The negative association of zinc content with crude protein content indicates that high crude protein of lentil seeds may result in decreasing zinc concentration. This can be explained by the fact that zinc plays an important role in protein synthesis in plants (74). Moreover, it has been documented that zinc is required as a cofactor for the activity of various enzymes (75). In our study, micronutrient contents and crude protein showed a nonsignificant correlation with yield, making it possible to develop cultivars with high micronutrient concentrations in combination with high yield. In contrast, several studies suggested positive correlation between iron and grain yield in lentil and chickpea (65,76).
In the present study, the multivariate analysis provides a useful mechanism for pinpointing the components that determine the variation when considering several traits simultaneously. Principal Component Analysis (PCA) indicated that most of the variation among tested genotypes was attributed to the number of filled pods and nutritional value. The identified genotypes coincided with those obtained when using the hierarchical cluster analysis (HCA). These methods were efficient in classifying lentil genotypes based on similar seed nutritional traits and grain yield, resulting in 3 clusters under both normal and stress conditions. Among the three clusters, the second cluster of the normal condition conciliates all qualities, lentil genotypes in this group are characterized by high productivity, high levels of iron concentration and crude protein, and by a moderate level of zinc concentration. Even though seed yield was markedly reduced under stress conditions, the plants were able to defend themselves to some extent from heat stress. Genotypes from clusters I and II were able to perform better and showed high grain yield per plant and a moderate level of micronutrient (Fe and Zn) and crude protein contents. However, these defenses failed under the combined heat-drought stress, resulting in a very low amount of crude protein and poor yield. Regardless, cluster I genotypes showed moderate yield and micronutrient contents.
The challenge of biofortification programs is to produce grains with high yields and superior nutritional quality, including high amounts of protein and several essential mineral elements (77). A Positive association between Fe and Zn concentrations provides the opportunity to combine these two traits in the same genotype. Furthermore, the broad-sense heritability estimate was 64 and 68% for Fe and Zn concentrations, respectively, indicating a high heritability of Fe and Zn accumulation in lentil (68,76). Our result revealed that Fe and Zn levels seem to be yieldindependent, which indicates that an increase could follow an increase in the concentration of simultaneous nutrients in grain yield, or it should at least be unaffected.
Analysis of variance results showed a high genetic variation for yield, crude protein, and micronutrient concentrations. The presence of genetic variability can be due to the diverse nature of tested genotypes. Carrying out a further analysis like cluster analysis led to the distinction of three groups of genotypes, in which yield, crude protein and iron and zinc concentrations were most responsible for their discrimination. Thus, it was possible to identify 10 lentil genotypes with a combination of high grain yield and high micronutrient concentrations, and hence, have the potential for biofortification. These genotypes can be used in genotype x environment studies to investigate whether different environments influence mineral concentrations before being included in the breeding program.

CONCLUSION
The present study revealed the existence of adequate genetic variation for Fe (48-109 mg kg −1 ), Zn (31-65 mg kg −1 ), and crude protein content (22.5-32.0%) in lentil genotypes. Our results also indicate that rising temperatures alone and in combination with drought stress at the time of flowering and seed filling period severely reduced grain yield and seed quality components. We report a significant reduction in Fe (30%), Zn (35%), and crude protein (57%) levels in lentil seeds under combined heat-drought stress environments compared to 18, 22, and 14% reductions in Fe, Zn, and crude protein, respectively, under heat stress alone. The impact of combined heat-drought stress was more severe on nutritional quality compared to the mild effect of heat stress alone. This difference was more pronounced for protein compared to micronutrients. This study showed a differential response of lentil genotypes to water and heat stress and identified lentil genotypes with the least effect of stress on their grain yield and nutritional quality. These genotypes require detailed study to understand the mechanism behind no impact of heat and drought on their nutritional quality. Non-significant correlations of micronutrient and crude protein with grain yield further make it possible to develop high yielding cultivars with better micronutrient concentrations and protein. With climate change and variability, the crop is expected to encounter frequent and intense water and heat stress in rainfed dry areas. Therefore, it is essential to develop high yielding nutrient-rich varieties with the least effect of heat and water stress for sustainable lentil production.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

AUTHOR CONTRIBUTIONS
HC performed the research work and prepared first draft of the manuscript. HC and KH raised experiments in fields and recorded phenotyping data. HC and AE-B analysized the seed samples for grain nutrient analysis in quality lab. NE and FM performed statistical analysis of data, SK planned and supervised research activity, and AS, FM, DT, and SK contributed in the final draft of the manuscript. All authors reviewed and approved the final version of the manuscript.

FUNDING
The partial funding support provided by the CGIAR Research Program on Grain Legumes and Dryland Cereals (GLDC) and Government of India is duly acknowledged.