Next Article in Journal
Exogenous Application of Melatonin and Methyl Jasmonate as a Pre-Harvest Treatment Enhances Growth of Barhi Date Palm Trees, Prolongs Storability, and Maintains Quality of Their Fruits under Storage Conditions
Next Article in Special Issue
Cereal Germination under Low Oxygen: Molecular Processes
Previous Article in Journal
In Vitro Photoprotective, Anti-Inflammatory, Moisturizing, and Antimelanogenic Effects of a Methanolic Extract of Chrysophyllum lucentifolium Cronquist
Previous Article in Special Issue
Light-Mediated Reduction in Photosynthesis in Closed Greenhouses Can Be Compensated for by CO2 Enrichment in Tomato Production
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High-Temperature and Drought Stress Effects on Growth, Yield and Nutritional Quality with Transpiration Response to Vapor Pressure Deficit in Lentil

1
International Center for Agricultural Research in the Dry Areas (ICARDA), Rabat 10112, Morocco
2
Laboratoire de Biotechnologie et de Physiologie Végétales, Centre de Recherche BioBio, Faculté des Sciences, Mohammed V University Rabat, Rabat 10112, Morocco
3
AgroBioSciences (AgBS) Research Division, Mohammed VI Polytechnic University, Lot 660 Hay Moulay Rachid, Ben Guerir 43150, Morocco
4
Biotechnology Research Unit, Regional Center of Agricultural Research of Rabat, National Institute of Agricultural Research (INRA), Rabat 10090, Morocco
5
Vellore Institute of Technology (VIT), VIT School of Agricultural Innovations and Advanced Learning (VAIAL), Vellore 632014, Tamil Nadu, India
6
International Center for Agricultural Research in the Dry Areas (ICARDA), Beirut 1108 2010, Lebanon
*
Authors to whom correspondence should be addressed.
Plants 2022, 11(1), 95; https://doi.org/10.3390/plants11010095
Submission received: 29 November 2021 / Revised: 20 December 2021 / Accepted: 23 December 2021 / Published: 28 December 2021
(This article belongs to the Special Issue Plant Physiology and Environmental Stresses)

Abstract

:
High temperature and water deficit are among the major limitations reducing lentil (Lens culinaris Medik.) yield in many growing regions. In addition, increasing atmospheric vapor pressure deficit (VPD) due to global warming causes a severe challenge by influencing the water balance of the plants, thus also affecting growth and yield. In the present study, we evaluated 20 lentil genotypes under field conditions and controlled environments with the following objectives: (i) to investigate the impact of temperature stress and combined temperature-drought stress on traits related to phenology, grain yield, nutritional quality, and canopy temperature under field conditions, and (ii) to examine the genotypic variability for limited transpiration (TRlim) trait in response to increased VPD under controlled conditions. The field experiment results revealed that high-temperature stress significantly affected all parameters compared to normal conditions. The protein content ranged from 23.4 to 31.9%, while the range of grain zinc and iron content varied from 33.1 to 64.4 and 62.3 to 99.3 mg kg−1, respectively, under normal conditions. The grain protein content, zinc and iron decreased significantly by 15, 14 and 15% under high-temperature stress, respectively. However, the impact was more severe under combined temperature-drought stress with a reduction of 53% in protein content, 18% in zinc and 20% in iron. Grain yield declined significantly by 43% in temperature stress and by 49% in the combined temperature-drought stress. The results from the controlled conditions showed a wide variation in TR among studied lentil genotypes. Nine genotypes displayed TRlim at 2.76 to 3.51 kPa, with the genotypes ILL 7833 and ILL 7835 exhibiting the lowest breakpoint. Genotypes with low breakpoints had the ability to conserve water, allowing it to be used at later stages for increased yield. Our results identified promising genotypes including ILL 7835, ILL 7814 and ILL 4605 (Bakria) that could be of great interest in breeding for high yields, protein and micronutrient contents under high-temperature and drought stress. In addition, it was found that the TRlim trait has the potential to select for increased lentil yields under field water-deficit environments.

Graphical Abstract

1. Introduction

Lentil (Lens culinaris Medik.) is an important cool-season food legume that plays a significant role in human and animal nutrition as well as in maintaining and improving soil fertility [1,2]. Lentil grains are highly digestible and nutritious, representing a true staple crop in the dryland of many developing countries with affordable levels of dietary proteins (22–35%), vitamins, fiber, carbohydrates and essential micronutrients such as zinc and iron [3]. As lentil is highly nutritious and can be cooked more quickly than other pulses, it is the pulse most preferred by poor rural households worldwide, particularly those living in developing nations such as India, Bangladesh, Pakistan and Ethiopia [4]. However, micronutrient deficiencies are currently becoming more serious due to the high consumption of carbohydrate-rich cereal-based diets, which are low in micronutrients and vitamins [5,6]. Micronutrient malnutrition, commonly known as hidden hunger, affects more than 2 billion people in developing countries, while all over the world, >3 billion people suffer from micronutrient deficiencies [7,8,9,10,11]. In the last decade, efforts have been made to reduce the risks of malnutrition using several approaches including nutrient supplementation, diet diversification and food fortification. However, genetic biofortification using conventional and transgenic breeding methods has been recognized as the most effective and sustainable method for developing new crop varieties to combat hidden hunger and provide essential micronutrients and vitamins to the poor population through daily diet [6,12,13]. Highly nutritional genotypes have been developed through genetic biofortification in lentil [14], chickpea [15] and wheat [16], helping to minimize micronutrient deficiencies substantially in recent decades. However, the increase in the world population to around 10 billion by 2050 will further escalate malnutrition [17,18], and more efforts will be required to meet the human population’s nutritional needs for a healthy life. Enhancing the production of lentil would support food and nutritional security and improve the livelihoods of resource-poor farmers in developing nations. Still, lentil productivity in rainfed regions is expected to suffer from fluctuations in climate change, mainly due to the increased incidence of drought and higher temperatures [19].
High-temperature and drought stress have been reported as the major environmental factors that can markedly affect plant productivity and the quality of many cultivated crops [20]; however, their impact on cool-season food legumes appears to be more serious [21,22]. A temperature higher than 30 °C causes stress in most cool-season food legume crops. Both high-temperature and drought stress hamper plant growth by disturbing the normal physiology and morphology, thereby influencing an array of processes including growth, floral development, carbohydrates, protein content in grains, and micronutrient concentration (zinc and iron), which ultimately affect grain yield and nutritional quality [23,24]. Furthermore, the combined effect of high-temperature and drought stress on crops could be more severe than the individual stress impact. The reproductive stages of crops are more vulnerable to drought, high temperature and combined stress than the vegetative stages [14,25]. During anthesis and flowering periods, both stresses led to fertilization failures because of reduced pollen and ovule function and inhibited pollen development and sterility in legumes [26,27] and cereals [28,29]. While considering the current change in climate, it is expected that both the severity of high temperatures and the risk of drought will increase mainly in the subtropical region, impacting food security [30].
The impacts of high temperatures [31,32], water stress [33,34] and combined temperature-drought stress on lentil growth, yield and composition have been documented [20,25,35,36]. However, more information is needed on the physiological responses of lentil and their mechanisms to develop climate-resilient management strategies and new varieties. Since lentil is cultivated mainly in environments with high vapor pressure deficit (VPD) conditions (hot and dry areas) during the post-rainy season, where light rainfall is often noted, the study of the limited transpiration (TRlim) trait could be crucial in water-deficit conditions [37]. VPD is defined as the difference between the amount of moisture in the air and the moisture the air can hold when saturated. Temperature and relative humidity are the two major factors controlling VPD variation [38]. The high VPD, which usually occurs due to exacerbation of plant water stress in the middle of the day, affects plant photosynthesis, growth, yield and physiology [39]. The intensification of atmospheric water demand highly influences plant hydraulic status, and the key behind these impacts is plant stomatal behavior [40,41]. Regulation of stomatal conductance in high VPD conditions by limiting transpiration restricts the rate of water use and increases transpiration efficiency, helping to conserve water and support plant growth later in the season when drought develops [42,43]. TRlim in response to high VPD has been reported in several crop species including wheat (Triticum aestivum L.) [41,44], maize (Zea mays L.) [45,46], sorghum (sorghum bicolor L.) [47,48], soybean (Glycine max L.) [49,50], peanut (Arachis hypogaea L.) [51,52], cowpea (Vigna unguiculata L.) [53,54], chickpea (Cicer arietinum L.) [55,56] and lentil [57]. Based on the available literature, the study by Guiguitant et al. [57] is the only investigation that examined lentil genotype differences in transpiration response to high VPD conditions. In the study by Guiguitant et al. [57], 17 lentil genotypes were tested for the expression of TRlim and almost all tested lines showed a VPD breakpoint at approximately 3.4 kPa. In addition, Guiguitant et al. [57] reported a possible yield benefit of developing genotypes expressing the TRlim trait, especially in regions with low rainfall, and that the impact of this trait on lentil productivity differs with geography and environments. Furthermore, several studies have shown increases in crop yield by 25 to 75% in water-limited areas based on crop simulation models [45,58,59,60]. Thus, developing and/or identifying drought- and heat-tolerant varieties is crucial in improving lentil yield under water-limited and high-temperature environments, which would directly contribute to food and nutritional security in the near future.
Based on the results from preliminary field screening [35], 20 lentil genotypes were selected as primary candidates for temperature and drought stress studies. Two independent experiments were conducted (i) to investigate the impact of high temperature and drought stress on traits associated with phenology, grain yield, nutritional quality and canopy temperature under field conditions, and (ii) to examine the possible genotypic variation in lentil for transpiration response to high VPD conditions over controlled environments. In the second experiment, transpiration was measured for plants subjected to different levels of VPD ranging from 1.20 to 4.50 kPa. The generated information may help to inform selection strategies and decisions within breeding programs and identify options to mitigate the impacts of high temperatures and water deficit on lentil productivity.

2. Results

2.1. Weather Data

Accumulated rainfall for normal planting was about 225 mm between December and April, while the late planting received around 26 mm mostly before the anthesis period. Under normal conditions, the minimum and maximum temperatures differed from −2.40 to 11.7 °C and 9.18 to 30.6 °C, respectively, during the vegetative stage, while they varied from 0.35 to 8.58 °C and 13.3 to 30.2 °C during the reproductive stage, respectively. Under the late planting, minimum and maximum temperatures were from 0.35 to 15.3 °C and from 13.3 to 37.5 °C respectively, before the flowering stage, and from 5.61 to 14.9 °C and from 22.6 to 36.9 °C, respectively, during the reproductive stage. Under normal planting, the RH ranged from 52 to 98% in the vegetative stage and 41 to 91% in the reproductive stage. It varied in late planting from 32 to 96% in the vegetative stage and from 42 to 79% during the reproductive stage. The maximum VPD registered in normal planting was 2.11 kPa and 3.25 kPa in late planting (Figure 1).

2.2. Combined Analysis of Variance

The combined analysis of variance indicated a significant difference among studied genotypes for all traits (Table S1). ANOVA showed a highly significant (p < 0.001) effect of temperature stress and combined temperature-drought stress on all measured characteristics. However, genotype × treatment interaction was significant for all traits except NUPP and HSW. Furthermore, the coefficient of variation (CV) among treatments was less than 13% for most traits, but 35% for both GYP and NFPP, with 25 and 49% for NTPP and NUPP, respectively.

2.3. Temperature and Combined Heat-Drought Impact on Plant Phenology

PLH declined significantly, by 19.2%, in temperature stress and by 29.6% in the combined temperature-drought stress compared to its level in normal conditions. EGV ranged from 17.0 to 30.0 cm with an overall mean of 23.1 cm under normal conditions, whereas high-temperature and combined temperature-drought stress reduced EGV by 26.2 and 22.9%, respectively. DFF varied between 82 and 89 days in normal conditions, and a reduction of 28 days was observed for both stressed treatments. Likewise, days to 50% flowering (D50F) declined significantly, by 27 days, under both temperature and combined temperature-drought stress conditions (Table 1).
Temperature stress alone and in combination with drought severely reduced the D50P by 26 and 27 days, respectively. Under normal conditions, plants achieved physiological maturity in about 127 days, while days to maturity decreased significantly by 19 and 20 days under temperature stress and combined temperature-drought stress, respectively. NFPP varied from 40.8 to 234.2 in normal conditions and was significantly reduced by 44 to 69% in temperature stress and by 48 to 71% in combined temperature-drought stress, respectively. Similarly, NUPP decreased substantially by 21–62% and 3–43% in temperature stress and combined temperature-drought stress, respectively. NTPP ranged from 56 to 240 pods in normal conditions with an average of 122.8, while it decreased significantly by 44–72% under temperature stress and 46–75% under combined temperature-drought stress conditions. GYP ranged from 1.59 to 8.86 g with an overall mean of 4.13 g under normal conditions. The late planting treatments showed a significant reduction in grain yield, by 43% under temperature stress and 49% under the combined stress. HSW varied from 1.93 and 4.97 g in normal conditions, with a mean of 2.64 g, whereas it was decreased significantly by 8% and 19% under the temperature stress and combined temperature-drought stress, respectively (Table 1).

2.4. Influence of High-Temperature and Combined Temperature-Drought Stress on Fe, Zn and Protein Contents

Crude protein (CP) varied from 23.4 to 31.9% under normal conditions, with a mean of 28.9%. CP decreased significantly to 15% due to temperature stress and 53% under combined temperature-drought stress. Under normal conditions, the mean Zn concentration was 47.9 mg kg−1, ranging from 33.1 to 64.4 mg kg−1. High-temperature stress and combined temperature-drought stress significantly reduced Zn concentration by 14 and 18%, respectively. Similar results were observed for Fe with a reduction of 15% due to temperature stress and 20% under combined temperature-drought stress. Fe concentration varied from 62.3 to 99.3 mg kg−1 under normal conditions (Table 1).

2.5. Correlation Coefficient between Measured Traits

A significant positive correlation was observed between PLH and EGV (r = 0.38) under normal conditions. However, a highly significant negative correlation (p < 0.01) was observed among EGV, DFF and D50P. HSW had a negative correlation with DFF, D50F and D50P, and it had a positive correlation with EGV (r = 0.62, p < 0.01). GYP was positively correlated with NFPP and NTPP (r = 0.97 and r = 0.96 at p < 0.01, respectively). For quality traits, there were no correlations between CP and phenological and yield parameters under normal conditions. Positive correlations were observed between Fe and Zn (r = 0.68), D50P (r = 0.52) and DM (r = 0.38, p < 0.05), and negative correlations with EGV (r = −0.61, p < 0.01), HSW (r = −0.39, p < 0.05) and CP (r = −0.35, p < 0.05). In addition, a highly negative correlation was observed between Zn and CP (r = −0.47, p < 0.01) in normal conditions (Table S2).
A significant negative correlation was found between EGV and DFF, D50F, and D50P under both temperature stress and combined temperature-drought stress (Table S3). Similarly to normal conditions, GYP had a significant positive correlation with NTPP and NFPP (r = 0.90) under the two stressed conditions. Under temperature stress, Fe content showed a positive correlation with D50F (r = 0.37, p < 0.05), GYP (r = 0.34, p < 0.05) and Zn (r = 0.46, p < 0.01). A significant and positive correlation between Zn and Fe (r = 48, p < 0.01) was also observed in combined temperature-drought stress conditions. In addition, Zn was found to be positively correlated with HSW (r = 0.34, p < 0.05) and negatively correlated with CP under temperature stress (r = −0.36, p < 0.05) and combined temperature-drought stress (r = −0.45, p < 0.01). Likewise, a positive correlation was observed between CP and DM (r = 0.36, p < 0.05) under temperature stress.

2.6. Principal Component Analysis for Stress and Normal Conditions

Principal component analysis (PCA) was performed to study the relationships between traits in the studied lentil genotypes. Under normal conditions, the first three PCAs (PC1: 31.3%, PC2: 22.6% and PC3: 16.5%) explained 70.4% of the total variability (Table 2). The results of PCA showed that EGV, D50P, HSW and Fe content were the most important traits contributing to PC1 of distinct origin. In PC2, which described 22.6% of the total variance, NFPP, GYP and D50F demonstrated large contributions, while Zn content and CP accounted for much of the total variance in PC3 (Table S4). These 20 genotypes were clustered into three groups based on hierarchical cluster analysis (Figure S1A). The check Bakria (ILL 4605) and three genotypes (ILL 5919, ILL 7813 and ILL 3484) were identified in group 1 and had the highest GYP (5.55 g) and HSW (3.27 g), and moderate Zn, Fe and CP. In group 2, ten genotypes were classified and recorded the highest Fe (85.1 mg kg−1) and Zn (54.1 mg kg−1), with medium GYP and HSW. Six genotypes were grouped in cluster 3 and exhibited the highest CP (29.9%), medium Fe (72.4 mg kg−1) and HSW (2.5 g), and the lowest GYP (Table 3). The distribution of genotypes and measured traits along the first PC axes under normal conditions are shown in Figure 2.
For temperature stress, PC1 and PC2 (PC1: 34.0% and PC2: 24.3%) explained 58.3% of the total variation, while PC3 showed 16.5% of the total variation. PC1 was positively correlated with DFF, D50F, D50P, HSW and Zn, whereas PC2 was positively correlated with GYP, NTPP, NFPP and Fe content. Furthermore, negative correlations were observed between PC1 and EGV, NTPP and NFPP. PC3 had a significant positive correlation with Zn content and negative association with CP and DM (Table 2). Hierarchical cluster analysis was conducted, and the genotypes were grouped into three clusters (Figure S1B). Cluster 1 grouped only the check Bakria (ILL 4605), a moderate heat- and drought-tolerant variety, which had the highest HSW (4.95 g) and Zn and Fe contents with a mean of 56.6 and 73.8 mg kg−1, respectively. Five heat-tolerant genotypes (ILL 7814, ILL 6104, ILL 7835, ILL 7833 and ILL 6338) and one moderately heat- and drought-tolerant genotype (ILL 6363) were identified in cluster 2. This group had the highest GYP with a mean of 3.30 g and CP with an average of 24.7%, while Zn and Fe concentrations were 42.5 and 69.7 mg kg−1, respectively. Thirteen genotypes including two heat-tolerant lines (ILL 8029 and ILL 7286), three combined temperature-drought tolerant genotypes (ILL 6362, ILL 7804 and ILL 6075), five moderately tolerant lines and three susceptible lines were found in cluster 3 (Figure S1B). This group had moderate GYP with an average of 1.97 g, HSW with 2.35 g and CP with a mean of 24.4%, whereas Zn and Fe concentrations were 39.3 and 64.8 mg kg−1, respectively (Table 3). The biplot of PC1 and PC2 of the three clusters and traits under temperature stress conditions is shown in Figure 3.
The first three PCAs in the combined temperature-drought stress explained 66.2% of the entire dissimilarity (PC1:34.5%, PC2:19.3% and PC3:12.4%). PC1 showed a significant positive correlation with D50F, D50P, DM, GYP, NTPP, NFPP and CP, and a negative correlation with Zn, Fe and EGV. PC2 was positively correlated with GYP, NTPP and NFPP, and negatively correlated with D50F and HSW (Table 2). Based on hierarchical cluster analysis (Figure S1C), 11 genotypes including two heat-tolerant lines (ILL 6104 and ILL 8029), four heat-drought tolerant lines (ILL 6075, ILL 7804, ILL 7835 and ILL 7814), three moderately tolerant lines (ILL 7223, ILL 7819 and ILL 8025) and three susceptible genotypes (ILL 5919, ILL 7813 and ILL 7820) were found in group 1. This group was characterized by the highest Zn and Fe concentrations with an average of 41.3 and 65.0 mg kg−1, respectively, and medium GYP (1.84 g). Group 2 comprised two moderately tolerant genotypes (ILL 4605 and ILL 6359), and had the highest CP and HSW, with a mean of 14.2% and 3.56 g, respectively. Three heat-tolerant lines (ILL 7833, ILL 7286 and ILL 6338), two moderately tolerant (ILL 6363 and ILL 3484) and one heat-drought tolerant genotype (ILL 6362) were grouped in cluster 3. This group had the highest GYP with an average of 2.97 g, and moderate concentrations of Zn and Fe with a mean of 37.1 and 59.9 mg kg−1, respectively (Table 3). The biplot of PC1 and PC2 of the three clusters and measured traits under the combined temperature-drought stress is illustrated in Figure 4.

2.7. Canopy Temperature Variation under Normal and Stress Conditions

Analysis of variance for canopy temperature (CT) showed a significant difference between genotypes and treatments; however, there were no significant differences for genotype x treatment interaction (Table S5). Under normal conditions, the CT of plants increased significantly from an average of 29.1 °C at day 100 to 31.8 °C at day 115 after sowing. Under normal conditions, the minimum and maximum CT were 26.6–32.7 °C at day 100 and 29.7–33.4 °C at day 115 (Table S6). Under temperature stress, the plants recorded an average of 29.2 °C and 31.7 °C at day 65 and day 70, respectively. The average CT increased significantly from 31.5 °C at day 80 to 31.9 °C at day 90. The lowest maximum CT was recorded at day 65 with 31.3 °C, while the highest was achieved at day 90 with 33.8 °C. Under combined temperature-drought stress, the CT increased significantly from an average of 29.5 °C at day 65 to 31.9 at day 70. On day 80 and day 90, the CT was 31.6 and 31.7 °C, respectively. Plants at day 65 recorded the lowest maximum CT (31.3°C), while, the highest was recorded at day 90 with 34.0 °C.

2.8. Transpiration Response to VPD under Controlled Environments

The results revealed that nine lentil genotypes were well represented by the two-segment regression with a breakpoint (BP) (Table 4). However, 11 genotypes exhibited a linear response and failed to express the limited transpiration (TRlim) trait during the increase in VPD (Table 5). The R2 of the genotypes with BP ranged from 0.60 to 0.87 and varied from 0.66 to 0.95 for the genotypes without TRlim. The results displayed high variation between the assessed genotypes in TR to increasing VPD. The value of BP ± SE of the genotypes expressing TRlim varied from 2.76 ± 0.43 to 3.51 ± 0.54 kPa with an average of 3.14 ± 0.66 (Table 4).
High dissimilarity was distinguished for the slopes below and above BP increasing TR with further increases in VPD. The slope at VPD above BP ranged from −1.63 ± 5.34 to 28.1 ± 11.3 mg with an average of 14 mg H2O m−2 s−1 kPa−1. In contrast, the slope below BP varied from 18.1 ± 5.49 to 44.7 ± 7.37 mg with an average of 34.5 mg H2O m−2 s−1 kPa−1. Genotype ILL 8029 showed the highest BP at about 3.51 kPa, while the two genotypes ILL 7835 and ILL 7833 recorded the lowest BP at around 2.76 and 2.81 kPa, respectively (Figure 5). The genotype ILL 7833 expressed the lowest TR during the increase in VPD, with 18.1 and −1.63 mg H2O m−2 s−1 kPa−1 for left and right slopes, respectively. Similarly, ILL 7814 and ILL 7835 were characterized by a very low TR after expressing the BP, with slopes greater than the BP ranging from 5.59 to 10.8 mg H2O m−2 s−1 kPa−1. Among genotypes that did not express a BP, the check Bakria (ILL 4605) had the lowest slope with a TR of 16.5 ± 1.51 mg H2O m−2 s−1 kPa−1. Oppositely, ILL 8025 exhibited the highest slope with 36.5 ± 2.78 mg H2O m−2 s−1 kPa−1. The results also revealed that ILL 5919, ILL 7813 and ILL 8025 had the maximum TR at around 175 mg H2O m−2 s−1 kPa−1, which lost TRlim trait expression during the increase in VPD. ILL 6262 recorded the highest TR with 210 mg H2O m−2 s−1 kPa−1, followed by the genotypes ILL 6338 and ILL 3484, which had a TR of 190 mg H2O m−2 s−1 kPa−1 (Table 4).

3. Discussion

3.1. Phenology, Yield and Nutritional Quality under High-Temperature and Drought Stress

In this study, delayed sowing to expose the reproductive stage to heat stress significantly affected phenological traits, grain yield components, nutritional quality parameters and canopy temperature in lentil genotypes. However, combining drought with heat stress through limited irrigation had a more severe impact on plant growth, resulting in significantly lower grain yield and nutritional quality. These results are consistent with our previous studies that confirmed the influence of water limitation and high temperature on lentil [14,35]. Similar findings were also observed in chickpea [61], lentil [31,36], soybean [62] and groundnut [63] where combined temperature-drought stress was more detrimental than temperature stress alone.
Our findings revealed that the number of filled pods decreased because of high-temperature and combined temperature-drought stress. This is in agreement with earlier studies that showed temperature and/or drought stress negatively influenced the reproductive processes, mainly pollen and ovule fertility [64,65], resulting in the reduction of filled pods. Recently, Jiang et al. [66] described that temperature stress reduced the number of pollen grains per anther, induced smaller pollen grains and increased reactive oxygen species (ROS) production in pollen grains in pea. Still, it did not affect ROS accumulation in ovules and ovule number per ovary. Furthermore, high-temperature exposure when young floral buds were visible at the first formed reproductive node was more destructive to flower retention, seed set, pod development and seed yield compared to heat exposure started later, when flowers at the second reproductive node were fully open. Temperature and drought stress also impact pollen vacuolization, which plays a vital role in increasing the volume of pollen grains with the accumulation of cytoplasmic components [67]. A recent study by Fábián et al. [29] showed that the combination of high temperature and water stress altered the phenology of the plants, reduced pollen viability, shortened the duration of gametogenesis and grain filling, and modified the morphology and anatomy of the pistils. Functionality of female and male reproductive parts was reduced by 34 and 66%, respectively.
Our study showed a decrease in seed weight due to temperature stress, while combined temperature-drought stress magnified the impact. It is a fact that during the seed-filling stage, the decline in photosynthesis and leaf sucrose metabolism in response to temperature stress and drought stress led to lower carbohydrate availability for import into developing seeds [68]. Previous investigations suggested seed weight is reduced in response to high-temperature and/or drought stress during grain filling by altering endosperm cells [69] and reducing starch accumulation [70]. Starch is the central component in seeds for major global staple crops. A decrease in source strength and carbon availability for starch biosynthesis through the grain filling period will definitely reduce seed size and weight [71,72]. Recent studies revealed a decline in the activity of enzymes involved in starch biosynthesis under high-temperature or water stress. Eventually, their combination contributed to a decrease in the starch accumulation rate and duration and starch content [73,74]. Sehgal et al. [36] indicated a significant decrease in starch content under temperature stress and water stress in lentil seeds. However, the combined temperature-drought stress caused the most severe reduction due to the inhibition of starch synthesizing enzyme (starch synthase) and sucrose synthesizing activity (sucrose synthase), which were previously identified to have a correlated effect on seed size in chickpea [75].
A significant reduction in protein content was reported under temperature stress; however, the combined temperature-drought stress caused more decrease, which is in agreement with previous observations in lentil [14,36], and in chickpea [76] through reduced function of PSII, weakened nitrogen anabolism and strengthened protein catabolism. Previous findings by Zahedi et al. [77] and Triboï et al. [78] revealed a stable relation between protein content and the total quantity of nitrogen in wheat grain. They suggested that the synthesis of storage proteins is limited mainly by the availability of nitrogen. Similarly, Liu et al. [79] also reported a decrease in crude protein fractions in drought-stressed alfalfa (Medicago sativa L.) due to water limitation in addition to a reduction in N fixation. For that, improvement of N fixation could lead to higher biomass production and water-use efficiency and may increase pod yield under drought stress conditions.
High temperature induced a significant reduction in grain micronutrient concentrations. The combined temperature-drought stress further aggravated the decline, by 18% for Zn and 20% for Fe, compared to that under normal conditions. Similar findings were observed in other legume crops such as lentil, where iron and zinc contents were dramatically reduced in response to temperature and drought stress due to decreased root nutrient uptake, with reduced root biomass and metabolic rate [80] or by direct damage to roots [81]. In addition, the reduction in transpiration rate because of water deficit may also decrease nutrient absorption and the effectiveness of their use by the plants [36,82]. Choukri et al. [14] suggested that the decrease in iron and zinc was attributed to decreasing water availability under heat and drought stress conditions. Furthermore, Hummel et al. [83] reported that the variations in zinc, iron and protein under drought stress conditions are more affected by weather conditions than genotype. The present study also revealed a negative correlation between Zn and protein content in high-temperature stress and heat-drought stress conditions, and Fe and protein content in the combined temperature-drought stress, which are in disagreement with previous studies conducted by Ghanbari et al. [84] and Impa et al. [85]. They reported a positive correlation between protein content and micronutrients under temperature and drought stress conditions. However, Fe and Zn were positively correlated under the three treatments, which was also reported in earlier findings in lentil [86], chickpea [87,88] and pea [89]. Furthermore, our results showed a significant increase in canopy temperature of plants under high-temperature stress and combined temperature-drought stress. This effect might be explained by inhibition of stomatal conductance and transpiration reducing leaf water content, which resulted from an increase in leaf temperatures [90,91].

3.2. Transpiration Response to High VPD under Controlled Conditions

Our findings showed significant genetic variation among lentil genotypes in transpiration response to increasing VPD under controlled environments. Eleven genotypes consistently increased TR with increasing VPD, while nine exhibited a distinct response by limiting their TR when VPD reached about 3.1 kPa. Our outcomes indicated that most of the genotypes displayed lower BP compared to what was reported by Guiguitant et al. [57], who showed a BP at approximately 3.4 kPa in the majority of lentil genotypes, with the lowest BP at 3.31 kPa. In our study, two tolerant genotypes (ILL 7833 and ILL 7835) limited their TR at early VPD (2.8 kPa), representing the lowest BP over all the already published VPD experiments in lentil. Restriction of TR under high VPD was also reported in many other species such as chickpea [56], peanut [51], soybean [49], sorghum [92] and maize [93], in which the BP values mostly fluctuated from 1.1 to 2.7 kPa. However, Schoppach and Sadok [94] showed high BP ranging from 2.4 to 3.9 kPa in wheat, while Sinclair et al. [95] described a VPD threshold ranging from 2.6 to 3.38 kPa in peanut.
Over the whole range of tested VPD, almost all of the tolerant genotypes showed lower TRlim compared to the sensitive lentil lines, which exhibited a continuously increasing TR with increasing VPD. These results support our field findings [35], and they are in good agreement with the conclusion of Belko et al. [54] in cowpea. Genotypes with TRlim may have considerable potential for increased soil water conservation, indicating somehow more effective control of TR under limited-water conditions [43]. Several findings reported that genotypes with a conservative water behavior have the opportunity of using the conserved soil water to preserve physiological activity during the grain filling period and produce higher grain yield than genotypes without the TRlim trait in late-season water deficit conditions [37,42,43,57,60,96].
Interestingly, the check variety Bakria, a moderately tolerant line to drought and heat stress, exhibited the lowest TR among all tested lines with 16 mg H2O m−2 s−1 kPa−1 under increasing VPD, and probably had the lowest conductance among the studied genotypes. Among the tolerant lines, ILL 7833 and ILL 8029 had the lower TR presented in slope 1 with 18.1 and 22.7 H2O m−2 s−1 kPa−1, respectively. However, genotype ILL 8029 registered the maximum BP at 3.51 kPa, which was likely to be disadvantageous in expressing the TRlim trait under high levels of VPD [94,97]. Variation in TR shown in the two slopes, which normalized water use efficiency, can be explained by the difference in the hydraulic conductance to water flux in the plants. Several studies have indicated that the TRlim trait under high evaporative conditions of VPD is likely due to low hydraulic conductivity in leaves between the xylem and into the guard cells, which increased water use efficiency [93,98,99]. The data described here do not allow support for the claim of a hydraulic conductance limitation in leaves; however, the difference in VPD breakpoint values among lentil genotypes may indicate variation in stomatal conductance between studied genotypes.
Since the experiment was conducted under controlled conditions, the variations in temperatures and relative humidity were the main key factors to reach higher VPD levels; however, the high temperatures (>35 °C) may influence TR response of the plants. A recent study by Kar et al. [100] revealed that an increase in VPD had a powerful positive impact on plant water loss irrespective of the effect of temperature, light and relative humidity. Although previous studies have reported similar observations, the effects of temperature, relative humidity, radiation and water on TR had minor effects compared to VPD [94,101,102].
Genotypes with TRlim under high VPD offer an essential candidate for breeding programs to improve lentil yields in water deficit areas [57]. A low BP represents an imperative approach by the plants to maintain the greatest water conservation during critical drought periods when high VPD promotes TR. Our outcomes indicate that studied lentil genotypes could be used to improve yield potential over a range of drought-stressed environments. Hence, tolerant genotypes such as ILL 7833 and ILL 7835, which had the lowest BP, may be desirable for drought-prone conditions and improve lentil productivity in regions affected by drought stress. The TR profiles generated from the current investigation present an opportunity to understand the physiological behaviors of tolerant and sensitive lentil genotypes under controlled conditions. However, additional studies are needed to examine more physiological and genetic responses of lentil material under field and controlled conditions.

4. Materials and Methods

4.1. Plant Material

The material consisted of 20 lentil genotypes selected from our previous field screening study [35], where we evaluated a focused identification of germplasm strategy (FIGS) set against drought and temperature stresses. In the present study, we included moderately heat tolerant, moderately drought tolerant, heat tolerant, heat sensitive, drought tolerant and drought sensitive lines (Table S7).

4.2. Field Experiment

The germplasm was assessed under three independent experiments at the ICARDA experimental station Marchouch (33.56° N, 6.63° W, 392 m altitude) during the 2016–17 cropping season. These three experiments were deemed to represent three treatments, namely (i) normal date of planting (Treatment A), (ii) late planting with irrigation (Treatment B) and (iii) late planting without irrigation (Treatment C). All of the treatments were planted in an alpha lattice design with two replications. In the three treatments, each genotype was planted in a three-row plot of 1 m length, with a spacing of 0.30 m between rows. In each row, seeds were sown at 2 cm depth, maintaining 10 cm space between plants. Treatment A resulted in optimal growing conditions (150 mm well-distributed rainfall and temperature below 27 °C, without any heat or water stress to the plants). Treatment B (planted 65 days after normal planting date with irrigation at field capacity throughout crop duration) imposed high-temperature stress. In both late planting treatments, the plants were synchronized with temperatures above 32 °C during the reproductive stage. Regular irrigation at field capacity avoided any water stress to the plants during their growth and development. Treatment C (planted 65 days after normal planting date without irrigation during the reproductive stage) imposed combined high-temperature and drought stress. Treatment A was planted on 27 December 2016, without any supplementary irrigation during the cropping period, as the crop received enough well-distributed rainfall during the cropping season. Treatments B and C were planted on 1 March 2017. Irrigation was performed regularly to sustain water supply at field capacity using a sprinkler system throughout the crop duration in treatment B. In contrast, irrigation was stopped at the flower initiation stage onward in treatment C to impose water stress (<5 mm rainfall during the reproductive stage) combined with the temperature stress. All field management followed standard agricultural practices in lentil production [103].

4.3. Data Collection

Data were collected on early growth vigor (EGV) after 1 month of sowing, plant height (PLH), days to first flowering (DFF), days to 50% of flowering (D50F), days to 50% of podding (D50P) and days to physiological maturity (DM) on a plot basis. Five plants were randomly selected from each plot to measure number of total pods plant−1 (NTPP), number of filled pods plant−1 (NFPP) and number of unfilled pods plant−1 (NUPP), grain yield plant−1 (GYP) and hundred-seed weight (HSW). Canopy temperature (CT) was determined using a thermal infrared camera FLIR T4xx-series (Model T420-KIT-15). Thermal images were captured between 11:00 and 13:00 GMT time on a sunny day and analyzed using FLIR tools+ software (FLIR systems, Teledyne FLIR, Wilsonville, OR, USA). The seeds harvested from these experiments were analyzed for protein content and micronutrient contents (zinc and iron).

4.4. Crude Protein Content

A 0.3 g amount of ground lentil seeds from each treatment was digested with sulfuric acid and selenium mixture following the modified Kjeldahl procedure [104]. Based on the color reaction between ammonium and a weakly alkaline mixture of sodium salicylate and sodium hypochlorite, a UV visible spectrophotometer was used to determine the color development in the samples at 560 nm. Next, crude protein content (CP) was measured using nitrogen values multiplied by 6.25, and triplicate analyses were performed for each sample.

4.5. Iron and Zinc Determination

Seeds were ground by a Cyclone mill (Twister, 10 mm–250 um, Retsch). Iron (Fe) and zinc (Zn) concentrations were measured using a modified HNO3 and H2O2 method [105]. In the digestion block (QBlock series, Horiba), 0.5 g of each ground sample was placed in individual tubes and digested with 6 mL nitric acid (HNO3), followed by heat treatments at 90 °C for 1 h. To each tube, 3 mL of 30% hydrogen peroxide (H2O2) was added, and sample digestion was continued by heating for 15 min at 90 °C, and then 3 mL of 6 M hydrochloric acid (HCL) was added. Once samples were cooled, the solutions were filtered and diluted to 10 mL with distilled water. The mineral content analysis was carried out by inductivity coupled plasma-optical emission spectrometry (ICP-OES); (iCAP-7000 Duo, Thermo Fisher Scientific, Waltham, MA, USA) at the Cereals and Legumes Quality Laboratory, ICARDA, Rabat, Morocco.

4.6. Plant Growth Conditions in the Greenhouse

An independent experiment was conducted during summer 2018 in the greenhouse at ICARDA, Rabat, Morocco. Three seeds were sown at a depth of 2.5 cm in plastic pots (15 cm in diameter and 20 cm in height) filled with 1.5 kg of 50% sandy loam soil and 50% of compost garden soil (Floragard Vertriebs-GmbH product, Oldenburg) that included 18-10-20 N-P-K fertilizer. Here, the same set of 20 genotypes assessed under field conditions were studied under controlled conditions. Five replications were used for each genotype. The pots were positioned in a randomized design on steel grid platform installed at mid-height of 1.3 m. In each pot, a small hole was opened at the bottom of the end cap to facilitate drainage. Plants were grown under well-watered conditions, and the temperature in the greenhouse was maintained at 25 °C/18 °C (day/night) for 1 month. The photosynthetic photon flux density in the greenhouse was approximately 600 μmol m−2 s−1.
The day before the measurement, pots were watered to full capacity until their water flowed from the bottom of the pots, then the pots were allowed to drain overnight. After drainage, the pots were bagged with plastic bags and covered with polyethylene balls to keep soil evaporation to a minimum. A thermo-hygrograph sensor (Tinytag Ultra 2 TGU-4500 Gemini Datalogger Ltd., Chichester, UK) was positioned over the plant to measure the temperature and RH at 5 min intervals. Hourly atmospheric VPD was calculated on the basis of average temperature and RH. Plants were exposed to each humidity and temperature treatment for 1 h and then reweighed to measure the final weight. On an hourly basis, the transpiration rate (TR) was calculated as the weight difference between successive measurements using a balance with a resolution of 0.1 g. The measurement was conducted in the morning from 7:00 am (Morocco standard time) at low VPD until 7:00 pm, when the VPD decreased following the midday maximum. Measurements were first started with the low VPD (0–1.5 kPa) treatment, the medium VPD treatment (1.5–2.5 kPa), and finally, the high VPD (2.5–4.5 kPa) treatment. The experiment in the greenhouse allowed us to control the temperature and RH, exposing the plants to the expected ranges of VPD (1.18 to 4.5 kPa). The average temperature and RH during the experiment fluctuated from 19 to 40 °C and from 30 to 80%, respectively. A high VPD level was achieved by increasing temperature and decreasing RH inside the greenhouse. At the end of the experiment, the leaf area was measured using a leaf area meter, and the TR was expressed as water loss per unit of leaf area and time (g H2O2 cm−2 h−1).

4.7. Statistical Analysis

4.7.1. Field Data Analysis

Analysis of variance was carried out using the general linear model (GLM) in IBM SPSS statistics 23. Duncan’s post-hoc was applied to compare differences between the mean values at p < 0.05. Pearson’s correlation coefficient was determined for the three treatments. Principal component analysis (PCA) was performed using the Factoextra and FactoMineR packages in R version 4.1.0 and RStudio version 1. 3.1093 [106,107]. In addition, hierarchical cluster analysis was performed using Ward’s squared Euclidean distance method with the dendextend R package [108].

4.7.2. Vapor Pressure Deficit Analysis

The individual data for each VPD treatment of each genotype were used in the regression analysis of transpiration response. The data were first tested to fit the data to a two-segment linear regression using GraphPad Prism 7 (GraphPad Software Inc.). The results of a successful regression fit to the two-segment model were the coefficients defining two intersecting linear regressions:
If VPD < BP, TR = intercept 1 + slope 1 (VPD)
If VPD ≥ BP, TR = intercept 2 + slope 2 (VPD)
where BP is the breakpoint between the two linear segments; it is also an output of the GraphPad analysis and an estimate of the standard error for BP. The intercept represents the constant of the first and second line segments; the two slopes were statistically compared within GraphPad for a significant difference at p ≤ 0.05 level. In case there was a significant difference, the two-segment model characterized the outcomes for that genotype. All data for a genotype were assumed to fit a single linear regression model if two slopes were not identified to be significantly different.

5. Conclusions

The outcomes of the present study validate our previous field screening results. Almost all tolerant genotypes showed a high adaptation to high temperatures and combined temperature-drought stress under field conditions and exhibited a fairly clear TRlim at high VPD under controlled conditions. These genotypes would contribute to water saving in soil and improve the productivity under terminal heat or water-limited environments. Interestingly, two tolerant genotypes (ILL 7833 and ILL 7835) exhibited the lowest breakpoint found in lentil, which may be exclusively desirable for harsh regions affected by water deficit. Genotypes such as ILL 7835, ILL 7814 and Bakria (ILL 4605) combine good nutritional quality with moderate to high yield under both heat and combined heat-drought stress, which is of great interest for breeding programs. The parameters of the tested lentil genotypes in response to high temperatures, water deficit and higher VPD would facilitate formulating new experimental strategies to shed light on the molecular and genetic basis of heat and drought tolerance mechanisms in the future. Hence, combining crop physiology with molecular techniques could be a promising approach to improve lentils under high temperatures, water deficit and higher VPD.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/plants11010095/s1, Table S1: ANOVA analysis for 20 lines under normal and late planting trials at Marchouch, Table S2: Correlation between tested traits under normal conditions at Marchouch, Table S3: Correlation for late planting under combined temperature-drought stress (above) and temperature stress (below), Table S4: Contribution of traits to the first three PCA dimensions under normal conditions, temperature stress, and combined temperature-drought stress, Table S5: Analysis of variance for canopy temperature between normal conditions, temperature stress and combined temperature-drought stress at Marchouch, Table S6: Minimum, maximum, mean of canopy temperature under normal conditions, temperature stress and combined temperature-drought stress at Marchouch, Table S7: Description of origin, IG number and classification of lentil genotypes used in the current study. MHT: Moderately heat tolerant; MDT: Moderately drought tolerant, HT: heat tolerant, HS: heat sensitive, DT: drought tolerant, DS: drought sensitive, Figure S1. Dendrogram of normal (A), temperature stress (B) and combined temperature-drought stress (C) conditions.

Author Contributions

Conceptualization and methodology, N.E.H., H.C., M.E.G. and S.K.; Statistical analysis, N.E.H.; validation, S.K.; investigation, N.E.H., H.C. and K.H.; data curation, N.E.H.; writing—original draft preparation, N.E.H.; writing—review and editing, H.C., M.E.G., A.S., R.M., K.R., F.M. and S.K.; project administration and supervision, S.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was undertaken as part of and funded by the CGIAR Research Program on Grain Legumes and Dryland Cereals (GLDC).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article and Supplementary Materials.

Acknowledgments

The authors wish to recognize the intense work of the field technical staff.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sarker, A.; Kumar, S. Lentils in production and food systems in West Asia and Africa. Grain Legum. 2011, 57, 56–58. [Google Scholar]
  2. Erskine, W.; Sarker, A.; Kumar, S. Crops that feed the world Investing in lentil improvement toward a food secure world 3. Food Secur. 2011, 3, 127–139. [Google Scholar] [CrossRef]
  3. Kumar, S.; Rajendran, K.; Kumar, J.; Hamwieh, A.; Baum, M. Current knowledge in lentil genomics and its application for crop improvement. Front. Plant Sci. 2015, 6, 78. [Google Scholar] [CrossRef] [Green Version]
  4. Reda, A. Lentil (Lens Culinaris Medikus) Current Status and Future Prospect of Production in Ethiopia. Adv. Plants Agric. Res. 2015, 2, 00040. [Google Scholar] [CrossRef] [Green Version]
  5. Bouis, H.; Low, J.; McEwan, M.; Tanumihardjo, S. Biofortification: Evidence and Lessons Learned Linking Agriculture and Nutrition; FAO: Rome, Italy; WHO: Geneva, Switzerland, 2013. [Google Scholar]
  6. Jha, A.B.; Warkentin, T.D. Biofortification of Pulse Crops: Status and Future Perspectives. Plants 2020, 9, 73. [Google Scholar] [CrossRef] [Green Version]
  7. Cakmak, I.; Pfeiffer, W.H.; McClafferty, B. REVIEW: Biofortification of Durum Wheat with Zinc and Iron. Cereal Chem. J. 2010, 87, 10–20. [Google Scholar] [CrossRef] [Green Version]
  8. Depar, N.; Rajpar, I.; Memon, M.Y.; Imtiaz, M. Mineral nutrient densities in some domestic and exotic rice genotypes. Pak. J. Agric. Agric. Eng. Vet. Sci. 2011, 27, 134–142. [Google Scholar]
  9. Kumar, J.; Gupta, D.S.; Kumar, S.; Gupta, S.; Singh, N.P. Current Knowledge on Genetic Biofortification in Lentil. J. Agric. Food Chem. 2016, 64, 6383–6396. [Google Scholar] [CrossRef] [PubMed]
  10. Rehman, A.; Farooq, M.; Ozturk, L.; Asif, M.; Siddique, K.H.M. Zinc nutrition in wheat-based cropping systems. Plant Soil 2017, 422, 283–315. [Google Scholar] [CrossRef]
  11. Rehman, A.; Farooq, M.; Ullah, A.; Nadeem, F.; Im, S.Y.; Park, S.K.; Lee, D.-J. Agronomic Biofortification of Zinc in Pakistan: Status, Benefits, and Constraints. Front. Sustain. Food Syst. 2020, 4, 591722. [Google Scholar] [CrossRef]
  12. Saltzman, A.; Birol, E.; Bouis, H.; Boy, E.; De Moura, F.F.; Islam, Y.; Pfeiffer, W.H. Biofortification: Progress toward a more nourishing future. Glob. Food Secur. 2013, 2, 9–17. [Google Scholar] [CrossRef]
  13. Bouis, H.E.; Saltzman, A. Improving nutrition through biofortification: A review of evidence from HarvestPlus, 2003 through 2016. Glob. Food Secur. 2017, 12, 49–58. [Google Scholar] [CrossRef] [PubMed]
  14. Choukri, H.; Hejjaoui, K.; El-Baouchi, A.; El Haddad, N.; Smouni, A.; Maalouf, F.; Thavarajah, D.; Kumar, S. Heat and Drought Stress Impact on Phenology, Grain Yield, and Nutritional Quality of Lentil (Lens culinaris Medikus). Front. Nutr. 2020, 7, 596307. [Google Scholar] [CrossRef] [PubMed]
  15. Ullah, A.; Al-Sadi, A.; Al-Subhi, A.M.; Farooq, M. Characterization of chickpea genotypes of Pakistani origin for genetic diversity and zinc grain biofortification. J. Sci. Food Agric. 2020, 100, 4139–4149. [Google Scholar] [CrossRef]
  16. Singh, D.; Rajawat, M.V.S.; Kaushik, R.; Prasanna, R.; Saxena, A.K. Beneficial role of endophytes in biofortification of Zn in wheat genotypes varying in nutrient use efficiency grown in soils sufficient and deficient in Zn. Plant Soil 2017, 416, 107–116. [Google Scholar] [CrossRef]
  17. Shah, F.-U.; Sharif, M.K.; Bashir, S.; Ahsan, F. Role of healthy extruded snacks to mitigate malnutrition. Food Rev. Int. 2018, 35, 299–323. [Google Scholar] [CrossRef]
  18. Steiner, G.; Geissler, B.; Schernhammer, E.S. Hunger and Obesity as Symptoms of Non-Sustainable Food Systems and Malnutrition. Appl. Sci. 2019, 9, 1062. [Google Scholar] [CrossRef] [Green Version]
  19. Johnson, N.; Johnson, C.R.; Thavarajah, P.; Kumar, S.; Thavarajah, D. The roles and potential of lentil prebiotic carbohydrates in human and plant health. Plants People Planet 2020, 2, 310–319. [Google Scholar] [CrossRef] [Green Version]
  20. El Haddad, N.; Sanchez-Garcia, M.; Visioni, A.; Jilal, A.; El Amil, R.; Sall, A.T.; Lagesse, W.; Kumar, S.; Bassi, F.M. Crop Wild Relatives Crosses: Multi-Location Assessment in Durum Wheat, Barley, and Lentil. Agronomy 2021, 11, 2283. [Google Scholar] [CrossRef]
  21. Awasthi, R.; Gaur, P.; Turner, N.; Vadez, V.; Siddique, K.; Nayyar, H. Effects of individual and combined heat and drought stress during seed filling on the oxidative metabolism and yield of chickpea (Cicer arietinum) genotypes differing in heat and drought tolerance. Crop. Pasture Sci. 2017, 68, 823–841. [Google Scholar] [CrossRef]
  22. Sultana, R.; Choudhary, A.K.; Pal, A.K.; Saxena, K.B.; Prasad, B.D.; Singh, R. Abiotic Stresses in Major Pulses: Current Status and Strategies. In Approaches to Plant Stress and their Management; Springer: New Delhi, India, 2014; pp. 173–190. [Google Scholar] [CrossRef]
  23. Maqbool, A.; Abrar, M.; Bakhsh, A.; Çalışkan, S.; Khan, H.Z.; Aslam, M.; Aksoy, E. Biofortification Under Climate Change: The Fight Between Quality and Quantity. In Environment, Climate, Plant and Vegetation Growth; Springer: Cham, Switzerland, 2020; pp. 173–227. [Google Scholar] [CrossRef]
  24. Johansson, E.; Branlard, G.; Cuniberti, M.; Flagella, Z.; Hüsken, A.; Nurit, E.; Peña, R.J.; Sissons, M.; Vazquez, D. Genotypic and Environmental Effects on Wheat Technological and Nutritional Quality. In Wheat Quality for Improving Processing and Human Health; Igrejas, G., Ikeda, T.M., Guzmán, C., Eds.; Springer International Publishing: Cham, Germeny, 2020; pp. 171–204. ISBN 978-3-030-34163-3. [Google Scholar] [CrossRef]
  25. Sehgal, A.; Sita, K.; Kumar, J.; Kumar, S.; Singh, S.; Siddique, K.H.M.; Nayyar, H. Effects of Drought, Heat and Their Interaction on the Growth, Yield and Photosynthetic Function of Lentil (Lens culinaris Medikus) Genotypes Varying in Heat and Drought Sensitivity. Front. Plant Sci. 2017, 8, 1776. [Google Scholar] [CrossRef] [Green Version]
  26. Kumari, V.; Roy, A.; Vijayan, R.; Banerjee, P.; Verma, V.; Nalia, A.; Pramanik, M.; Mukherjee, B.; Ghosh, A.; Reja, H.; et al. Drought and Heat Stress in Cool-Season Food Legumes in Sub-Tropical Regions: Consequences, Adaptation, and Mitigation Strategies. Plants 2021, 10, 1038. [Google Scholar] [CrossRef]
  27. Sita, K.; Sehgal, A.; Kumar, J.; Kumar, S.; Singh, S.; Siddique, K.; Nayyar, H. Identification of High-Temperature Tolerant Lentil (Lens culinaris Medik.) Genotypes through Leaf and Pollen Traits. Front. Plant Sci. 2017, 8, 744. [Google Scholar] [CrossRef] [Green Version]
  28. Zahra, N.; Wahid, A.; Hafeez, M.B.; Ullah, A.; Siddique, K.H.; Farooq, M. Grain development in wheat under combined heat and drought stress: Plant responses and management. Environ. Exp. Bot. 2021, 188, 104517. [Google Scholar] [CrossRef]
  29. Fábián, A.; Sáfrán, E.; Szabó-Eitel, G.; Barnabás, B.; Jäger, K. Stigma Functionality and Fertility Are Reduced by Heat and Drought Co-stress in Wheat. Front. Plant Sci. 2019, 10, 244. [Google Scholar] [CrossRef]
  30. Raza, A.; Razzaq, A.; Mehmood, S.S.; Zou, X.; Zhang, X.; Lv, Y.; Xu, J. Impact of Climate Change on Crops Adaptation and Strategies to Tackle Its Outcome: A Review. Plants 2019, 8, 34. [Google Scholar] [CrossRef] [Green Version]
  31. Sita, K.; Sehgal, A.; Bhandari, K.; Kumar, J.; Kumar, S.; Singh, S.; Siddique, K.; Nayyar, H. Impact of heat stress during seed filling on seed quality and seed yield in lentil (Lens culinaris Medikus) genotypes. J. Sci. Food Agric. 2018, 98, 5134–5141. [Google Scholar] [CrossRef] [PubMed]
  32. Bhandari, K.; Siddique, K.; Turner, N.; Kaur, J.; Singh, S.; Agrawal, S.K.; Nayyar, H. Heat Stress at Reproductive Stage Disrupts Leaf Carbohydrate Metabolism, Impairs Reproductive Function, and Severely Reduces Seed Yield in Lentil. J. Crop. Improv. 2016, 30, 118–151. [Google Scholar] [CrossRef]
  33. Ahmed, B.; Mitu, A.; Sultana, M.; Tasnim, N.; Talukder, A. Yield and Physiological Indices of Lentil Genotypes against Drought Stress. Bangladesh Agron. J. 2021, 23, 97–102. [Google Scholar] [CrossRef]
  34. Salehi, M.; Haghnazari, A.; Shekari, F.; Faramarzi, A. The Study of Seed Yield and Seed Yield Components of Lentil (Lens culinaris Medik) under Normal and Drought Stress Conditions. Pak. J. Biol. Sci. 2008, 11, 758–762. [Google Scholar] [CrossRef] [PubMed]
  35. El Haddad, N.; Rajendran, K.; Smouni, A.; Es-Safi, N.E.; Benbrahim, N.; Mentag, R.; Nayyar, H.; Maalouf, F.; Kumar, S. Screening the FIGS Set of Lentil (Lens culinaris Medikus) Germplasm for Tolerance to Terminal Heat and Combined Drought-Heat Stress. Agronomy 2020, 10, 1036. [Google Scholar] [CrossRef]
  36. Sehgal, A.; Sita, K.; Bhandari, K.; Kumar, S.; Kumar, J.; Prasad, P.V.; Siddique, K.; Nayyar, H. Influence of drought and heat stress, applied independently or in combination during seed development, on qualitative and quantitative aspects of seeds of lentil (Lens culinaris Medikus) genotypes, differing in drought sensitivity. Plant Cell Environ. 2018, 42, 198–211. [Google Scholar] [CrossRef] [Green Version]
  37. Ghanem, M.E.; Kibbou, F.E.-Z.; Guiguitant, J.; Sinclair, T.R. Opportunities to improve the seasonal dynamics of water use in lentil (Lens culinaris Medik.) to enhance yield increase in water-limited environments. Chem. Biol. Technol. Agric. 2017, 4. [Google Scholar] [CrossRef]
  38. Devi, M.J.; Reddy, V.R. Transpiration Response of Cotton to Vapor Pressure Deficit and Its Relationship With Stomatal Traits. Front. Plant Sci. 2018, 9, 1572. [Google Scholar] [CrossRef] [Green Version]
  39. Grossiord, C.; Buckley, T.N.; Cernusak, L.A.; Novick, K.A.; Poulter, B.; Siegwolf, R.T.W.; Sperry, J.S.; McDowell, N.G. Plant responses to rising vapor pressure deficit. New Phytol. 2020, 226, 1550–1566. [Google Scholar] [CrossRef] [Green Version]
  40. Kulkarni, M.; Soolanayakanahally, R.; Ogawa, S.; Uga, Y.; Selvaraj, M.G.; Kagale, S. Drought Response in Wheat: Key Genes and Regulatory Mechanisms Controlling Root System Architecture and Transpiration Efficiency. Front. Chem. 2017, 5, 106. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Li, Y.; Li, H.; Li, Y.; Zhang, S. Improving water-use efficiency by decreasing stomatal conductance and transpiration rate to maintain higher ear photosynthetic rate in drought-resistant wheat. Crop. J. 2017, 5, 231–239. [Google Scholar] [CrossRef]
  42. Sinclair, T.R.; Devi, J.; Shekoofa, A.; Choudhary, S.; Sadok, W.; Vadez, V.; Riar, M.; Rufty, T. Limited-transpiration response to high vapor pressure deficit in crop species. Plant Sci. 2017, 260, 109–118. [Google Scholar] [CrossRef]
  43. Ghanem, M.E.; Sinclair, T.R. Limited-transpiration rate under high atmospheric vapour pressure deficit: A trait for developing genotypes for water-deficit conditions. CAB Rev. 2019, 14, 1–11. [Google Scholar] [CrossRef]
  44. Medina, S.; Vicente, R.; Nieto-Taladriz, M.T.; Aparicio, N.; Chairi, F.; Diaz, O.V.; Araus, J.L. The Plant-Transpiration Response to Vapor Pressure Deficit (VPD) in Durum Wheat Is Associated With Differential Yield Performance and Specific Expression of Genes Involved in Primary Metabolism and Water Transport. Front. Plant Sci. 2019, 9, 1994. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Shekoofa, A.; Sinclair, T.R.; Messina, C.D.; Cooper, M. Variation among Maize Hybrids in Response to High Vapor Pressure Deficit at High Temperatures. Crop. Sci. 2016, 56, 392–396. [Google Scholar] [CrossRef]
  46. Shrestha, R.K.; Lei, P.; Shi, D.; Hashimi, M.H.; Wang, S.; Xie, D.; Ni, J.; Ni, C. Response of maize (Zea mays L.) towards vapor pressure deficit. Environ. Exp. Bot. 2020, 181, 104293. [Google Scholar] [CrossRef]
  47. Karthika, G.; Kholova, J.; Alimagham, S.; Ganesan, M.; Chadalavada, K.; Kumari, R.; Vadez, V. Measurement of transpiration restriction under high vapor pressure deficit for sorghum mapping population parents. Plant Physiol. Rep. 2019, 24, 74–85. [Google Scholar] [CrossRef] [Green Version]
  48. Gholipoor, M.; Sinclair, T.R.; Prasad, P.V.V. Genotypic variation within sorghum for transpiration response to drying soil. Plant Soil 2012, 357, 35–40. [Google Scholar] [CrossRef]
  49. Gilbert, M.; Holbrook, N.; Zwieniecki, M.; Sadok, W.; Sinclair, T. Field confirmation of genetic variation in soybean transpiration response to vapor pressure deficit and photosynthetic compensation. Field Crop. Res. 2011, 124, 85–92. [Google Scholar] [CrossRef]
  50. Sadok, W.; Sinclair, T.R. Genetic variability of transpiration response to vapor pressure deficit among soybean (Glycine max [L.] Merr.) genotypes selected from a recombinant inbred line population. Field Crop. Res. 2009, 113, 156–160. [Google Scholar] [CrossRef]
  51. Shekoofa, A.; Rosas-Anderson, P.; Sinclair, T.R.; Balota, M.; Isleib, T.G. Measurement of Limited-Transpiration Trait under High Vapor Pressure Deficit for Peanut in Chambers and in Field. Agron. J. 2015, 107, 1019–1024. [Google Scholar] [CrossRef]
  52. Pradhan, D.; Shekoofa, A.; Sinclair, T.R. Temperature effect on peanut (Arachis hypogaea L.) transpiration response to vapor pressure deficit and its recovery. J. Crop. Improv. 2018, 33, 177–186. [Google Scholar] [CrossRef]
  53. Manandhar, A.; Sinclair, T.R.; Rufty, T.W.; Ghanem, M.E. Leaf Expansion and Transpiration Response to Soil Drying and Recovery among Cowpea Genotypes. Crop. Sci. 2017, 57, 2109–2116. [Google Scholar] [CrossRef]
  54. Belko, N.; Diop, N.; Cisse, N.; Ehlers, J.; Vadez, V.; Zaman-Allah, M.; Zombré, G. Restriction of transpiration rate under high vapour pressure deficit and non-limiting water conditions is important for terminal drought tolerance in cowpea. Plant Biol. 2012, 15, 304–316. [Google Scholar] [CrossRef]
  55. Sivasakthi, K.; Tharanya, M.; Kholová, J.; Muriuki, R.W.; Thirunalasundari, T.; Vadez, V. Chickpea Genotypes Contrasting for Vigor and Canopy Conductance Also Differ in Their Dependence on Different Water Transport Pathways. Front. Plant Sci. 2017, 8, 1663. [Google Scholar] [CrossRef] [Green Version]
  56. Zaman-Allah, M.; Jenkinson, D.M.; Vadez, V. Chickpea genotypes contrasting for seed yield under terminal drought stress in the field differ for traits related to the control of water use. Funct. Plant Biol. 2011, 38, 270–281. [Google Scholar] [CrossRef] [Green Version]
  57. Guiguitant, J.; Marrou, H.; Vadez, V.; Gupta, P.; Kumar, S.; Soltani, A.; Sinclair, T.R.; Ghanem, M.E. Relevance of limited-transpiration trait for lentil (Lens culinaris Medik.) in South Asia. Field Crop. Res. 2017, 209, 96–107. [Google Scholar] [CrossRef]
  58. Drake, P.L.; Froend, R.H.; Franks, P.J. Smaller, faster stomata: Scaling of stomatal size, rate of response, and stomatal conductance. J. Exp. Bot. 2013, 64, 495–505. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Ghanem, M.E.; Marrou, H.; Soltani, A.; Kumar, S.; Sinclair, T.R. Lentil Variation in Phenology and Yield Evaluated with a Model. Agron. J. 2015, 107, 1967–1977. [Google Scholar] [CrossRef]
  60. Messina, C.D.; Sinclair, T.R.; Hammer, G.L.; Curan, D.; Thompson, J.; Oler, Z.; Gho, C.; Cooper, M. Limited-Transpiration Trait May Increase Maize Drought Tolerance in the US Corn Belt. Agron. J. 2015, 107, 1978–1986. [Google Scholar] [CrossRef] [Green Version]
  61. Devasirvatham, V.; Tan, D.K.Y. Impact of High Temperature and Drought Stresses on Chickpea Production. Agronomy 2018, 8, 145. [Google Scholar] [CrossRef] [Green Version]
  62. Jumrani, K.; Bhatia, V.S. Impact of combined stress of high temperature and water deficit on growth and seed yield of soybean. Physiol. Mol. Biol. Plants 2017, 24, 37–50. [Google Scholar] [CrossRef]
  63. Hamidou, F.; Halilou, O.; Vadez, V. Assessment of Groundnut under Combined Heat and Drought Stress. J. Agron. Crop. Sci. 2012, 199, 1–11. [Google Scholar] [CrossRef] [Green Version]
  64. Mahrookashani, A.; Siebert, S.; Hüging, H.; Ewert, F. Independent and combined effects of high temperature and drought stress around anthesis on wheat. J. Agron. Crop. Sci. 2017, 203, 453–463. [Google Scholar] [CrossRef]
  65. Ruan, Y.-L.; Jin, Y.; Yang, Y.-J.; Li, G.-J.; Boyer, J.S. Sugar Input, Metabolism, and Signaling Mediated by Invertase: Roles in Development, Yield Potential, and Response to Drought and Heat. Mol. Plant 2010, 3, 942–955. [Google Scholar] [CrossRef]
  66. Jiang, Y.; Davis, A.R.; Vujanovic, V.; Bueckert, R.A. Reproductive development response to high daytime temperature in field pea. J. Agron. Crop. Sci. 2018, 205, 324–333. [Google Scholar] [CrossRef]
  67. Pacini, E.; Dolferus, R. The Trials and Tribulations of the Plant Male Gametophyte—Understanding Reproductive Stage Stress Tolerance. In Abiotic and Biotic Stress in Plants-Recent Advances and Future Perspectives; IntechOpen: London, UK, 2016. [Google Scholar]
  68. Cohen, I.; Zandalinas, S.I.; Huck, C.; Fritschi, F.B.; Mittler, R. Meta-analysis of drought and heat stress combination impact on crop yield and yield components. Physiol. Plant. 2020, 171, 66–76. [Google Scholar] [CrossRef]
  69. Nicolas, M.E.; Gleadow, R.; Dalling, M.J. Effect of Post-anthesis Drought on Cell Division and Starch Accumulation in Developing Wheat Grains. Ann. Bot. 1985, 55, 433–444. [Google Scholar] [CrossRef]
  70. Balla, K.; Rakszegi, M.; Li, Z.; Békés, F.; Bencze, S.; Veisz, O. Quality of winter wheat in relation to heat and drought shock after anthesis. Czech J. Food Sci. 2011, 29, 117–128. [Google Scholar] [CrossRef] [Green Version]
  71. Thitisaksakul, M.; Jiménez, R.C.; Arias, M.C.; Beckles, D. Effects of environmental factors on cereal starch biosynthesis and composition. J. Cereal Sci. 2012, 56, 67–80. [Google Scholar] [CrossRef]
  72. Liu, X.; Wang, X.; Wang, X.; Gao, J.; Luo, N.; Meng, Q.; Wang, P. Dissecting the critical stage in the response of maize kernel set to individual and combined drought and heat stress around flowering. Environ. Exp. Bot. 2020, 179, 104213. [Google Scholar] [CrossRef]
  73. Lu, H.; Hu, Y.; Wang, C.; Liu, W.; Ma, G.; Han, Q.; Ma, D. Effects of High Temperature and Drought Stress on the Expression of Gene Encoding Enzymes and the Activity of Key Enzymes Involved in Starch Biosynthesis in Wheat Grains. Front. Plant Sci. 2019, 10, 1414. [Google Scholar] [CrossRef] [PubMed]
  74. Ning, P.; Yang, L.; Li, C.; Fritschi, F.B. Post-silking carbon partitioning under nitrogen deficiency revealed sink limitation of grain yield in maize. J. Exp. Bot. 2018, 69, 1707–1719. [Google Scholar] [CrossRef] [Green Version]
  75. Kumar, A.; Turner, N.C. Growth and sucrose synthase activity of developing chickpea (Cicer arietinum L.) seeds under field conditions. Aust. J. Crop Sci. 2009, 3, 20–27. [Google Scholar]
  76. Awasthi, R.; Kaushal, N.; Vadez, V.; Turner, N.; Berger, J.; Siddique, K.; Nayyar, H. Individual and combined effects of transient drought and heat stress on carbon assimilation and seed filling in chickpea. Funct. Plant Biol. 2014, 41, 1148–1167. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Zahedi, M.; McDonald, G.; Jenner, C.F. Nitrogen supply to the grain modifies the effects of temperature on starch and protein accumulation during grain filling in wheat. Aust. J. Agric. Res. 2004, 55, 551–564. [Google Scholar] [CrossRef]
  78. Triboi, A.-M.; Martre, P.; Triboï-Blondel, A. Environmentally-induced changes in protein composition in developing grains of wheat are related to changes in total protein content. J. Exp. Bot. 2003, 54, 1731–1742. [Google Scholar] [CrossRef] [PubMed]
  79. Liu, Y.; Wu, Q.; Ge, G.; Han, G.; Jia, Y. Influence of drought stress on afalfa yields and nutritional composition. BMC Plant Biol. 2018, 18, 13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Heckathorn, S.A.; Giri, A.; Mishra, S.; Bista, D. Heat Stress and Roots. In Climate Change and Plant Abiotic Stress Tolerance; Wiley Online Library: Hoboken, NJ, USA, 2013; pp. 109–136. [Google Scholar]
  81. Huang, L.; Baumgartl, T.; Mulligan, D. Is rhizosphere remediation sufficient for sustainable revegetation of mine tailings? Ann. Bot. 2012, 110, 223–238. [Google Scholar] [CrossRef] [Green Version]
  82. Bista, D.R.; Heckathorn, S.A.; Jayawardena, D.M.; Mishra, S.; Boldt, J.K. Effects of Drought on Nutrient Uptake and the Levels of Nutrient-Uptake Proteins in Roots of Drought-Sensitive and -Tolerant Grasses. Plants 2018, 7, 28. [Google Scholar] [CrossRef] [Green Version]
  83. Hummel, M.; Hallahan, B.F.; Brychkova, G.; Ramirez-Villegas, J.; Guwela, V.; Chataika, B.; Curley, E.; McKeown, P.C.; Morrison, L.; Talsma, E.F.; et al. Reduction in nutritional quality and growing area suitability of common bean under climate change induced drought stress in Africa. Sci. Rep. 2018, 8, 16187. [Google Scholar] [CrossRef]
  84. Ghanbari, A.A.; Shakiba, M.R.; Toorchi, M.; Choukan, R. Nitrogen changes in the leaves and accumulation of some minerals in the seeds of red, white and chitti beans (’Phaseolus vulgaris’) under water deficit conditions. Aust. J. Crop Sci. 2013, 7, 706–712. [Google Scholar] [CrossRef]
  85. Impa, S.; Perumal, R.; Bean, S.; Sunoj, V.J.; Jagadish, S.K. Water deficit and heat stress induced alterations in grain physico-chemical characteristics and micronutrient composition in field grown grain sorghum. J. Cereal Sci. 2019, 86, 124–131. [Google Scholar] [CrossRef]
  86. Khazaei, H.; Podder, R.; Caron, C.T.; Kundu, S.S.; Diapari, M.; Vandenberg, A.; Bett, K.E. Marker–Trait Association Analysis of Iron and Zinc Concentration in Lentil (Lens culinaris Medik.) Seeds. Plant Genome 2017, 10, 1–8. [Google Scholar] [CrossRef] [Green Version]
  87. Tan, G.Z.H.; Das Bhowmik, S.S.; Hoang, T.; Karbaschi, M.R.; Long, H.; Cheng, A.; Bonneau, J.P.; Beasley, J.T.; Johnson, A.A.T.; Williams, B.; et al. Investigation of Baseline Iron Levels in Australian Chickpea and Evaluation of a Transgenic Biofortification Approach. Front. Plant Sci. 2018, 9, 788. [Google Scholar] [CrossRef]
  88. Samineni, S.; Mahendrakar, M.D.; Shankar, N.; Hotti, A.; Chand, U.; Rathore, A.; Gaur, P.M. Impact of heat and drought stresses on grain nutrient content in chickpea: Genome-wide marker-trait associations for protein, Fe and Zn. Environ. Exp. Bot. 2021, 194, 104688. [Google Scholar] [CrossRef]
  89. Ma, Y.; Zeng, X.-M.; Zhang, Y.; Wang, N.; Zheng, Y.; Wang, G.; Chen, C. Impact of the Choice of Land Surface Scheme on a Simulated Heatwave Event: The Case of Sichuan-Chongqing Area, China. Adv. Meteorol. 2017, 2017, 1–17. [Google Scholar] [CrossRef]
  90. Sehgal, A.; Sita, K.; Siddique, K.; Kumar, R.; Bhogireddy, S.; Varshney, R.K.; HanumanthaRao, B.; Nair, R.M.; Prasad, P.V.V.; Nayyar, H. Drought or/and Heat-Stress Effects on Seed Filling in Food Crops: Impacts on Functional Biochemistry, Seed Yields, and Nutritional Quality. Front. Plant Sci. 2018, 9, 1705. [Google Scholar] [CrossRef] [Green Version]
  91. Anjum, S.A.; Xie, X.; Wang, L.; Saleem, M.F.; Man, C.; Lei, W. Morphological, physiological and biochemical responses of plants to drought stress. Afr. J. Agric. Res. 2011, 6, 2026–2032. [Google Scholar]
  92. Shekoofa, A.; Balota, M.; Sinclair, T. Limited-transpiration trait evaluated in growth chamber and field for sorghum genotypes. Environ. Exp. Bot. 2014, 99, 175–179. [Google Scholar] [CrossRef]
  93. Sunita, C.; Sinclair, T.R.; Messina, C.D.; Cooper, M. Hydraulic Conductance of Maize Hybrids Differing in Transpiration Response to Vapor Pressure Deficit. Crop. Sci. 2014, 54, 1147–1152. [Google Scholar] [CrossRef]
  94. Schoppach, R.; Sadok, W. Differential sensitivities of transpiration to evaporative demand and soil water deficit among wheat elite cultivars indicate different strategies for drought tolerance. Environ. Exp. Bot. 2012, 84, 1–10. [Google Scholar] [CrossRef]
  95. Sinclair, T.R.; Pradhan, D.; Shekoofa, A. Inheritance of limited-transpiration trait in peanut: An update. J. Crop. Improv. 2018, 32, 281–286. [Google Scholar] [CrossRef]
  96. Sinclair, T.R. Theoretical Analysis of Soil and Plant Traits Influencing Daily Plant Water Flux on Drying Soils. Agron. J. 2005, 97, 1148–1152. [Google Scholar] [CrossRef]
  97. Vadez, V.; Kholová, J.; Yadav, R.S.; Hash, C.T. Small temporal differences in water uptake among varieties of pearl millet (Pennisetum glaucum (L.) R. Br.) are critical for grain yield under terminal drought. Plant Soil 2013, 371, 447–462. [Google Scholar] [CrossRef] [Green Version]
  98. Sinclair, T.R.; Zwieniecki, M.A.; Holbrook, N.M. Low leaf hydraulic conductance associated with drought tolerance in soybean. Physiol. Plant. 2007, 132, 446–451. [Google Scholar] [CrossRef] [PubMed]
  99. Fan, D.-Y.; Dang, Q.-L.; Xu, C.-Y.; Jiang, C.-D.; Zhang, W.-F.; Xu, X.-W.; Yang, X.-F.; Zhang, S.-R. Stomatal Sensitivity to Vapor Pressure Deficit and the Loss of Hydraulic Conductivity Are Coordinated in Populus euphratica, a Desert Phreatophyte Species. Front. Plant Sci. 2020, 11, 1248. [Google Scholar] [CrossRef] [PubMed]
  100. Kar, S.; Tanaka, R.; Korbu, L.B.; Kholová, J.; Iwata, H.; Durbha, S.S.; Adinarayana, J.; Vadez, V. Automated discretization of ‘transpiration restriction to increasing VPD’features from outdoors high-throughput phenotyping data. Plant Methods 2020, 16, 140. [Google Scholar] [CrossRef]
  101. Seversike, T.M.; Sermons, S.M.; Sinclair, T.R.; Carter, T.E., Jr.; Rufty, T.W. Temperature interactions with transpiration response to vapor pressure deficit among cultivated and wild soybean genotypes. Physiol. Plant. 2013, 148, 62–73. [Google Scholar] [CrossRef] [PubMed]
  102. Chenu, K.; Van Oosterom, E.J.; McLean, G.; Deifel, K.S.; Fletcher, A.; Geetika, G.; Tirfessa, A.; Mace, E.S.; Jordan, D.; Sulman, R.; et al. Integrating modelling and phenotyping approaches to identify and screen complex traits: Transpiration efficiency in cereals. J. Exp. Bot. 2018, 69, 3181–3194. [Google Scholar] [CrossRef] [PubMed]
  103. Kumar, S.; Rajendran, K. Lentil Ontology-Crop Ontology Curation Tool. Crop Ontology 547 for Agricultural Data. CGIAR Centers. Available online: https://cropontology.org/term/CO_339:ROOT (accessed on 20 December 2021).
  104. Baethgen, W.E.; Alley, M.M. A manual colorimetric procedure for measuring ammonium nitrogen in soil and plant Kjeldahl digests. Commun. Soil Sci. Plant Anal. 1989, 20, 961–969. [Google Scholar] [CrossRef]
  105. Gupta, D.S.; Thavarajah, D.; McGee, R.J.; Coyne, C.J.; Kumar, S.; Thavarajah, P.; Program, I.B. Genetic diversity among cultivated and wild lentils for iron, zinc, copper, calcium and magnesium concentrations. Aust. J. Crop. Sci. 2016, 10, 1381–1387. [Google Scholar] [CrossRef]
  106. Kassambara, A.; Mundt, F. Factoextra: Extract and Visualize the Results of Multivariate Data Analyses. 2020. Available online: https://cran.r-project.org/package=factoextra (accessed on 17 May 2020).
  107. Husson, F.; Josse, J.; Lê, S.; Mazet, J. FactoMineR: Multivariate Exploratory Data Analysis and Data Mining. 2019. Available online: https://cran.r-project.org/package=FactoMineR (accessed on 1 May 2020).
  108. Galili, T. dendextend: An R package for visualizing, adjusting and comparing trees of hierarchical clustering. Bioinformatics 2015, 31, 3718–3720. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Rainfall, maximum and minimum temperatures, and vapor pressure deficit (VPD) during normal and late planting experiments in Marchouch, Morocco during 2016/17 season.
Figure 1. Rainfall, maximum and minimum temperatures, and vapor pressure deficit (VPD) during normal and late planting experiments in Marchouch, Morocco during 2016/17 season.
Plants 11 00095 g001
Figure 2. Biplot of the first two axes of the PCA analysis summarizing the relationship between the measured variables and the classification of lentil genotypes under normal conditions.
Figure 2. Biplot of the first two axes of the PCA analysis summarizing the relationship between the measured variables and the classification of lentil genotypes under normal conditions.
Plants 11 00095 g002
Figure 3. Biplot of the first two axes of the PCA analysis summarizing the relationship between the measured variables and the classification of lentil genotypes under temperature stress conditions. HDS, heat stress and combined heat-drought susceptible; HDT, heat-drought tolerant; HT; heat tolerant and MT, moderately tolerant to heat stress and combined heat-drought stress.
Figure 3. Biplot of the first two axes of the PCA analysis summarizing the relationship between the measured variables and the classification of lentil genotypes under temperature stress conditions. HDS, heat stress and combined heat-drought susceptible; HDT, heat-drought tolerant; HT; heat tolerant and MT, moderately tolerant to heat stress and combined heat-drought stress.
Plants 11 00095 g003
Figure 4. Biplot of the first two axes of the PCA analysis summarizing the relationship between the measured variables and the classification of lentil genotypes under the combined temperature-drought stress conditions. HDS, heat-drought susceptible; HDT, heat-drought tolerant; HT; heat tolerant and MT, moderately tolerant to heat-drought stress.
Figure 4. Biplot of the first two axes of the PCA analysis summarizing the relationship between the measured variables and the classification of lentil genotypes under the combined temperature-drought stress conditions. HDS, heat-drought susceptible; HDT, heat-drought tolerant; HT; heat tolerant and MT, moderately tolerant to heat-drought stress.
Plants 11 00095 g004
Figure 5. Transpiration rate (mg H2O m−2 s−1) response to increasing vapor pressure deficit (VPD in kPa). Genotypes ILL 7833, ILL 7835 and ILL 7814 represent tolerant accessions with two segmental regression, while ILL 4605, ILL 8025 and ILL 7286 represent susceptible genotypes with linear transpiration response to increasing VPD. The breakpoint (BP), slopes and R2 are indicated.
Figure 5. Transpiration rate (mg H2O m−2 s−1) response to increasing vapor pressure deficit (VPD in kPa). Genotypes ILL 7833, ILL 7835 and ILL 7814 represent tolerant accessions with two segmental regression, while ILL 4605, ILL 8025 and ILL 7286 represent susceptible genotypes with linear transpiration response to increasing VPD. The breakpoint (BP), slopes and R2 are indicated.
Plants 11 00095 g005
Table 1. Minimum, maximum and mean morphological parameters of 20 lentil genotypes under normal, temperature stress and temperature-drought stress conditions.
Table 1. Minimum, maximum and mean morphological parameters of 20 lentil genotypes under normal, temperature stress and temperature-drought stress conditions.
TraitNormalTemperature StressTemperature-Drought Stress
MinMaxMean ± SEMinMaxMean ± SEMinMaxMean ± SE
PLH24.7534.2128.28 a ± 3.3919.7527.2522.86 b ± 0.3014.3325.0319.74 c ± 0.37
EGV17.0030.0023.12 a ± 0.5613.0021.0017.07 b ± 0.3012.0021.0017.82 c ± 0.35
DFF82.0089.0084.95 a ± 0.2260.0066.0061.55 b ± 0.19590.0064.0061.15 b ± 0.23
D50F86.00100.0093.08 a ± 0.4565.0074.0068.53 b ± 0.31630.0074.0068.32 b ± 0.39
D50P98.00106.00101.78 a ± 0.3671.0080.0076.20 b ± 0.44700.0080.0074.77 c ± 0.41
DM124.00131.00127.20 a ± 0.32101.00105.00103.53 b ± 0.17100.00105.00101.70 c ± 0.17
NFPP40.83234.25111.45 a ± 7.4312.75129.4057.75 b ± 4.6612.00119.5052.07 b ± 4.14
NUPP3.5025.2511.34 a ± 0.881.3320.006.15 b ± 0.512.0024.507.89 b ± 0.79
NTPP56.00240.00122.79 a ± 7.2315.75134.0063.90 b ± 4.7114.00130.0059.96 b ± 4.17
GYP1.598.864.13 a ± 0.280.434.562.35 b ± 0.170.664.322.12 b ± 0.14
HSW1.934.972.64 a ± 0.091.595.232.42 b ± 0.111.425.512.15 c ± 0.13
CP23.4031.9028.91 a ± 0.3522.0027.0024.36 b ± 0.2011.0014.9013.34 c ± 0.16
Zn33.1064.4047.86 a ± 1.5227.4059.1041.12 b ± 1.4130.4050.9039.28 c ± 0.80
Fe62.3099.3078.19 a ± 1.4451.3081.6066.71 b ± 1.4252.4073.5062.86 c ± 0.81
Means for each variable followed by the same letters are not significantly different (p < 0.05). PLH, plant height; EGV, early growth vigor; DFF, days to first flowering; D50F, days to 50% flowering; D50P, days to 50% podding, DM, days to 95% maturity; NFPP, number of filled pods plant−1; NUPP, number of unfilled pods plant−1; NTPP, number of total pods plant−1; GYP, grain yield plant−1; HSW, hundred-seed weight; CP, crude protein; Zn, zinc content and Fe, iron content.
Table 2. Eigenvalues and eigenvectors of the three PCA dimensions under normal conditions, temperature stress and the combined temperature-drought stress.
Table 2. Eigenvalues and eigenvectors of the three PCA dimensions under normal conditions, temperature stress and the combined temperature-drought stress.
TraitsNormalTemperature StressTemperature-Drought Stress
PC1PC2PC3PC1PC2PC3PC1PC2PC3
PLH−0.52 *0.21 ns−0.07 ns−0.31 ns−0.02 ns0.01 ns−0.27 ns−0.13 ns0.19 ns
EGV−0.78 **−0.47 *0.05 ns−0.49 *−0.42 ns0.38 ns−0.64 **−0.05 ns0.14 ns
DFF0.50 *0.40 ns−0.50 *0.86 **0.30 ns0.02 ns0.64 **−0.44 ns0.18 ns
D50F0.38 ns0.62 **−0.58 **0.86 **0.44 ns0.01 ns0.74 **−0.44 *0.44 ns
D50P0.88 **0.13 ns0.06 ns0.79 **0.42 ns−0.21 ns0.77 **−0.33 ns0.30 ns
DM0.46 *0.16 ns0.45 *0.22 ns0.41 ns−0.70 **0.89 **−0.11 ns0.19 ns
NFPP−0.50 *0.83 **0.14 ns−0.59 **0.76 **0.23 ns0.55 *0.81 **0.16 ns
NTPP−0.44 *0.84 **0.13 ns−0.57 **0.76 **0.23 ns0.59 **0.78 **0.09 ns
NUPP0.48 *−0.29 **−0.08 ns0.05 ns0.20 ns0.04 ns0.22 ns−0.20 ns−0.13 ns
GYP−0.49 *0.81 **0.11 ns−0.44 ns0.79 **0.33 ns0.49 *0.76 **0.19 ns
HSW−0.71 **−0.30 ns0.42 ns0.60 **−0.09 ns0.42 ns0.13 ns−0.50 *0.35 ns
CP−0.14 ns−0.07 ns−0.65 **−0.34 ns0.32 ns−0.60 **0.52 *−0.30 ns−0.39 ns
Zn0.35 ns0.18 ns0.77 **0.49 *−0.08 ns0.61 **−0.61 **0.14 ns0.35 ns
Fe0.74 *0.24 ns0.51 **0.36 ns0.47 *0.44 ns−0.64 **0.05 ns0.55 *
Eigenvalue4.383.172.314.082.921.984.832.701.73
Variance explained (%)31.2922.6516.4934.0024.3216.5234.5119.3012.42
Total variance (%)31.2953.2970.4334.0058.3274.8534.5153.8166.23
** 0.01, * 0.05, ns not significant. PLH, plant height; EGV, early growth vigor; DFF, days to first flowering; D50F, days to 50% flowering; D50P, days to 50% podding, DM, days to 95% maturity; NFPP, number of filled pods plant−1; NUPP, number of unfilled pods plant−1; NTPP, number of total pods plant−1; GYP, grain yield per plant; HSW, hundred-seed weight; CP, crude protein; Zn, zinc content; Fe, iron content.
Table 3. Grain yield, hundred-seed weight, crude protein, zinc content and iron content of the three clusters for normal, temperature stress and combined temperature-drought stress.
Table 3. Grain yield, hundred-seed weight, crude protein, zinc content and iron content of the three clusters for normal, temperature stress and combined temperature-drought stress.
ClusterExperimentGYPHSWCPZnFe
Cluster INormal5.553.2729.7844.1269.50
High-temperature1.654.9522.5256.6373.85
Temperature-drought1.842.0012.9241.2665.02
Cluster IINormal4.062.4427.9454.1385.15
High-temperature3.302.1724.6642.4969.73
Temperature-drought1.283.5614.2533.8558.65
Cluster IIINormal3.302.5429.9439.9472.39
High-temperature1.972.3524.3739.3064.77
Temperature-drought2.971.9613.8737.1259.99
GYP; grain yield, HSW; 100-seed weight, CP; crude protein, Zn; zinc content and Fe; iron content.
Table 4. Outputs from the analysis of the two-segment linear regression used for fitting the TR to the VPD variations.
Table 4. Outputs from the analysis of the two-segment linear regression used for fitting the TR to the VPD variations.
GenotypeBreakpointLeft SlopeRight SlopeR2
BP ± SECI (95%)S1CI (95%)S2CI (95%)
ILL 34843.27 ± 0.953.01 to 4.5044.4832.68 to 59.1924.953.45 to 39.940.87
ILL 60753.15 ± 1.142.73 to 4.5128.9121.74 to 45.8417.370.25 to 28.460.82
ILL 61043.00 ± 0.253.00 to 3.4536.4025.72 to 47.0612.26−2.57 to 26.430.83
ILL 63383.37 ± 0.613.00 to 3.4543.7331.41 to 60.6315.42−9.64 to 38.610.81
ILL 63623.01 ± 0.992.76 to 3.4644.7134.40 to 71.3928.142.88 to 41.590.87
ILL 78143.36 ± 0.562.79 to 4.1734.4225.90 to 77.815.59−26.00 to 28.450.85
ILL 78332.81 ± 0.471.67 to 3.5718.149.59 to 38.87−1.63−13.49 to 8.320.60
ILL 78352.76 ± 0.432.23 to 3.3037.3223.76 to 50.8810.79−3.32 to 24.910.79
ILL 80293.51 ± 0.543.07 to 4.0122.6914.10 to 31.2913.10−10.92 to 37.130.70
The BP is the breakpoint (kPa), S1 and S2 are the left slope and right slope, respectively (mg H2O m−2 s−1 kPa−1), CI is the confidence interval.
Table 5. Results of the 11 genotypes exhibiting single linear regression fits of TR response to vapor pressure deficit. The table includes slope with standard error (SE), X-intercept and R2 values.
Table 5. Results of the 11 genotypes exhibiting single linear regression fits of TR response to vapor pressure deficit. The table includes slope with standard error (SE), X-intercept and R2 values.
GenotypeSlope ± SESlope CI (95%)X-InterceptX-Intercept CI (95%)R2
ILL 591927.23 ± 2.5921.96 to 32.50−0.58−1.34 to −0.060.77
ILL 460516.53 ± 1.5113.46 to 19.61−0.39−1.13 to 0.130.79
ILL 635931.18 ± 2.8225.44 to 36.92−0.13−0.82 to 0.360.79
ILL 636327.07 ± 2.7221.55 to 32.58−0.75−1.65 to −0.150.74
ILL 722325.95 ± 1.0123.88 to 28.020.06−0.20 to 0.280.95
ILL 728622.04 ± 2.9815.94 to 28.13−1.32−2.84 to −0.440.66
ILL 780430.32 ± 2.1725.93 to 34.710.12−0.35 to 0.490.82
ILL 781329.08 ± 3.6721.54 to 36.62−0.15−1.19 to 0.480.70
ILL 781924.22 ± 2.2719.62 to 28.83−0.20−0.92 to 0.300.76
ILL 782023.01 ± 2.6117.60 to 28.43−0.61−1.62 to 0.030.78
ILL 802536.48 ± 2.7830.81 to 42.160.22−0.27 to 0.590.85
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

El Haddad, N.; Choukri, H.; Ghanem, M.E.; Smouni, A.; Mentag, R.; Rajendran, K.; Hejjaoui, K.; Maalouf, F.; Kumar, S. High-Temperature and Drought Stress Effects on Growth, Yield and Nutritional Quality with Transpiration Response to Vapor Pressure Deficit in Lentil. Plants 2022, 11, 95. https://doi.org/10.3390/plants11010095

AMA Style

El Haddad N, Choukri H, Ghanem ME, Smouni A, Mentag R, Rajendran K, Hejjaoui K, Maalouf F, Kumar S. High-Temperature and Drought Stress Effects on Growth, Yield and Nutritional Quality with Transpiration Response to Vapor Pressure Deficit in Lentil. Plants. 2022; 11(1):95. https://doi.org/10.3390/plants11010095

Chicago/Turabian Style

El Haddad, Noureddine, Hasnae Choukri, Michel Edmond Ghanem, Abdelaziz Smouni, Rachid Mentag, Karthika Rajendran, Kamal Hejjaoui, Fouad Maalouf, and Shiv Kumar. 2022. "High-Temperature and Drought Stress Effects on Growth, Yield and Nutritional Quality with Transpiration Response to Vapor Pressure Deficit in Lentil" Plants 11, no. 1: 95. https://doi.org/10.3390/plants11010095

APA Style

El Haddad, N., Choukri, H., Ghanem, M. E., Smouni, A., Mentag, R., Rajendran, K., Hejjaoui, K., Maalouf, F., & Kumar, S. (2022). High-Temperature and Drought Stress Effects on Growth, Yield and Nutritional Quality with Transpiration Response to Vapor Pressure Deficit in Lentil. Plants, 11(1), 95. https://doi.org/10.3390/plants11010095

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop