Effect of High Temperature Stress During the Reproductive Stage on Grain Yield and Nutritional Quality of Lentil (Lens culinaris Medikus)

Introduction

Lentil (Lens culinaris Medikus) is an ancient and early domesticated food legume crop originating from the Near East ~10,000 years ago (1). It has been part of the human diet since the early days of agriculture. Lentil belongs to the family Fabaceae and ranks fifth in global pulse production after common bean (Phaseolus vulgaris L.), dry pea (Pisum sativum L.), chickpea (Cicer arietnum L.), and cowpea (Vigna unguiculate L.) (2). The current global lentil production is 5.8 MT. Global lentil production is expected to increase due to the growing demand for plant-based protein and environment-friendly agriculture as it has a nitrogen-fixing capability to reduce energy costs in crop production (3).

Lentil is an integral part of the global food system with low inputs and provides a range of micronutrients (4). It is a rich protein source with 20–31% proteins, providing essential amino acids to the human body (5). Besides, lentil is also a good source of B-complex vitamins; it also contains a significant concentration of folate (6). Lentil seeds are relatively low in fat (<1%) and are an excellent supplement of the human diet in fiber and prebiotic carbohydrates (7, 8). Lentil provides a broad range of beneficial bioactive compounds such as polyphenols (9). Therefore, it is not surprising that lentil is the most desired legume in many lentil-producing regions worldwide because of the nutritional profiles. Besides their high nutritional value, lentils typically cook faster than other legumes because of their seed size and thin seed coat (10). Cooking time is an essential attribute in evaluating pulses, and the shorter cooking time is crucial to consumer preference as any fast cooking dish and end-use. Anti-nutritional factors (trypsin inhibitors, phytic acid, and tannins) reduce the digestibility of protein, carbohydrates, and micronutrient bioavailability in the human diet (11). Phytic acid is the primary storage form of phosphorus in legume seeds and plays a vital role in seed germination and seedling growth (12). Phytic acid is a potent chelator of positively charged cations including potassium (K), magnesium (Mg), calcium (Ca), iron (Fe), zinc (Zn), and mangenese (Mn) (13). Once consumed in human foods, phytic acid binds to these minerals in the intestinal tract to form mixed salts that are generally eliminated. This phenomenon poses a major risk of micronutrients deficiency (14).

Iron and Zinc deficiencies are the most widespread micronutrient deficiency globally. Micronutrient malnutrition is primarily due to the heavy reliance on staple food crop-based diets with low mineral bioavailability (15). Developing biofortified lentil varieties with high micronutrients bioavailability is a feasible approach to combat hidden hunger. However, several studies reported the genotypic and environmental factors affecting phytic acid and micronutrient concentrations in legume crops (16, 17). Furthermore, due to climate change, changes in temperature during the reproductive stage of legumes' development primarily affect their yield and nutritional value. Temperatures over 32°C during flowering and podding stages cause damage to reproductive organs, leading to significant losses in lentil grain yield (18). Outcomes from our previous study on the effect of heat stress on lentil seed quality revealed the adverse effect of rising temperature, at flowering and podding period, on seed quality components (19). Previous studies reported a significant reduction in Fe (18%), Zn (22%), and protein (14%) levels in lentils and observed an adequate genetic variation for Fe (0.5–10.9 mg/100 g), Zn (3.1–6.5 mg/100 g), and crude protein content (22.5–32.0 g/100 g) (19). In this regard, 36 lentil genotypes were selected based on their performance under heat stress conditions to test their response to grain yield and nutritional quality under controlled conditions. To achieve precision on the existing genetic variability of the nutritional quality of lentil with response to heat stress and the impact of rising temperature on the phytic acid and it effects on micronutrient bioavailability and cooking time in lentil seed, the present study was undertaken with the following objectives: (a) investigate high temperature effect on micronutrients, crude protein, phytic acid and cookability of selected genotypes, (b) to describe the variability in seed phytic acid concentration and estimate the bioavailability of Zn and Fe.

Materials and Methods

Plant Material

A set of 36 lentil genotypes was selected based on preliminary results of the quality traits and grain yield (19). Selected genotypes were also characterized by early flowering (<60 days after planting). Seeds were obtained from the International Center for Agricultural Research in the Dry Areas (ICARDA), Rabat, Morocco. Description of tested genotypes is presented in Supplementary Table S1. Plants were grown in plastic pots (15 cm in diameter and 20 cm in height) with 1.5 kg of a mixture of sandy loam soil and compost garden soil in a ratio of 2:2 (w/w). Six seeds were sown in each pot at 2 cm depth during the last week of December 2018 and reduced to four per pot after emergence. A randomized complete block design with two replicates was employed. To minimize the effects of variation in light profiles, the position of the pots was changed weekly during the experiment.

Greenhouse Experiment

The experiment was conducted under two temperature regimes maintained in separate chambers designated no stressed and heat stressed. Plants under no stressed were grown with a mean day/night temperature of 22/12°C, whereas plants with heat stress were grown with a mean temperature of 25/14°C. For both chambers, the day length was 14 h (500 μmol m⁻² s⁻¹ light intensity). Relative humidity varied between 70 and 80% throughout the experimental period in both chambers. Heat stress was initiated with the pre-flowering stage by the increasing temperature at an average rate of 2°C per day to be the final temperature of 32°C at flowering stage. The plants remained under these conditions till maturity. Besides the heat stress, plants were watered and fertilized as per protocol developed at ICARDA, Morocco. At physiological maturity, plants were hand-harvested, threshed, and cleaned seeds were used to determine the grain yield and nutritional quality traits.

Grain Yield Traits

Numbers of total pods, filled pods, unfilled pods, biological yield (g), grain yield (g), and hundred-seed weight (g) were determined for each treatment. All observations were taken on four plants per genotype for each replicate.

Seed Size and Seed Shape Parameters

Lentil seeds were randomly selected to determine the hundred-seed weight and seed geometry by image analysis using a high-speed seed counting device (OptoAgri2, Optomachine, France). The OptoAgri2 combines a camera with a high resolution depending on the grain size, a laboratory balance, and software for image processing with an algorithm developed by OPTOmachines. This device measures seed area, perimeter, length, width, circularity, diameter, thickness, and rugosity. In addition to seed shape parameters, the software gives information about the hundred-seed weight.

Mineral Concentration

Mineral concentration was measured using a previously published method (20). Five hundred milligrams from each sample were weighed into a glass tube and then 6 mL concentrated (70%) nitric acid (HNO₃) was added. The tubes were then placed into a digestion block (QBlock series, Ontario, Canada) and heated for 60 min at 90°C. Three milliliters of 30% hydrogen peroxide (H₂O₂) were then added, and the samples were heated at 90°C for another 15 min, after which most of the residue was digested. Then, 3 mL of 6 M hydrochloric acid (HCl) was added to each sample. After cooling to ambient temperature, the volume was adjusted to 10 mL and then filtered. Fe and Zn concentrations were measured using inductively coupled plasma-optical emission spectroscopy (ICP-OES); (ICAP-7000 Duo, Thermo Fisher Scientific, France). Calibration curves for Fe and Zn were made using serial dilution from 0.1 to 10 mg L⁻¹. Data validation was done using lab references and the NIST standard references.

Protein Concentration

The total nitrogen content of the samples was determined using the Kjeldahl procedure (21). Samples were digested with a mixture of sulfuric acid, selenium, and salicylic acid. The mixture was then digested at 300°C for 5 h using a digestion block (QBlock series, Ontario, Canada). The filtrate was treated with 5.5 ml of the buffer solution, 4.4 ml of sodium nitroprusside, and 2 ml of sodium hypochlorite. The mixture is placed in the dark at 37°C for 15 min. The sample's absorbance was measured at 650 nm. The percentage of protein content was calculated by multiplying the Kjeldahl nitrogen by the conversion factor of 6.25.

Phytic Acid Concentration

A Megazyme kit was used to measure the phytic acid concentration (22). One gram of ground lentil seed was digested with 20 mL of HCl (0.66 M) solution in 50 ml flacon tubes and placed in a shaker overnight (15 h) at room temperature. The day after, 1 mL of the extract was subjected to several enzymatic reactions to release inorganic phosphorus (Pi) from phytic acid. The inorganic phosphorus (Pi) was then reacted with ammonium molybdate to form 12-molybdophosphoric acid, which is subsequently reduced under acidic conditions to molybdenum blue. The amount of molybdenum blue formed in this reaction is proportional to the amount of Pi present in the original sample and hence to phytic acid. A UV-Visible spectrophotometer (T80 series, pg instruments, UK) was used to measure the absorbance of molybdenum blue at 655 nm. The phosphorus solutions were prepared as described by the Megazyme manual, using a standard phosphorus solution (24 mL, 50 μg/mL), and treated as samples for the colorimetric determination of phosphorus.

The inorganic phosphorus (Pi) is quantified as phosphorus from the generated calibration curve. At the same time, the calculation of phytic acid content assumes that the amount of phosphorus measured is exclusively released from phytic acid and comprises 28.2% of phytic acid.

$\begin{array}{l}\text{Phytic acid}=\frac{\text{Phosphorus}\left[\text{g}/100\text{g}\right]}{0.282}\end{array}$

Phytic Acid/Iron and Phytic Acid/Zinc Molar Ratios

The molar ratios phytic acid/iron and phytic acid/zinc were calculated according to the following equation:

$\begin{array}{l}\text{PA}/\text{MN ratio}=\frac{\text{PA}/\text{M}{\text{W}}_{\text{AP}}}{\text{MN}/\text{M}{\text{W}}_{\text{MN}}}\end{array}$

where: PA is phytate acid content; MW_PA is PA molecular weight (660.04 g mol⁻¹); MN is micronutrient content (zinc or iron); MW_MN is micronutrient molecular weight (Zn = 65.4 g mol⁻¹; Fe = 55.85 g mol⁻¹).

Cooking Time

The automated Mattson Cooker (Canadian Grain Commission, Canada) consisted of a cooking rack with 25 perforated depressions and 25 weighted plungers (80 g) of 2 mm tip diameter was used to measure lentil seed cooking time. Before cooking, the samples (2 g of lentil) were soaked in 50 mL of distilled water and maintained at room temperature (22 ± 2°C) for 12 h. Twenty-five soaked seeds were randomly selected and then placed in each of the 25 perforations of the rack. The tip of each plunger was rested on top of the seed. The loaded Mattson Cooker was then placed into a 2 L beaker containing 1.5 L of boiling distilled water, and the beaker was covered to reduce evaporation. The apparatus was maintained over direct heat using a hot plate. During testing, the hot plate was set at 390°C to maintain water boiling. Lentil seeds are considered cooked when the plungers penetrate the seeds and touch the sensor on the fourth rack (sensor assembly); the Easy CT program automatically recorded cooking time for each seed. Cooking time is defined as the number of minutes required for 80% of lentil seeds to be pierced.

Statistical Analysis

Summary data was reported as range and mean values ± standard deviation. The coefficient and the p-value of the correlations were calculated using Pearson's correlation. Analysis of variance (ANOVA) was used to test the effects of treatments, genotypes, and genotype x treatment interaction on the measured traits. Principal component analysis (PCA) was performed using the Factoextra (23) and FactoMineR (24) packages in R version 4.1.0 and RStudio version 1. 3.1093. Hierarchical cluster analysis was conducted using Ward's squared Euclidean distance method with the dendextend R package (25).

Results

The combined analysis of variance showed significant genotypic differences for all evaluated traits of 36 tested genotypes (Supplementary Table S2). The effect of high temperature was also significant (p < 0.001) for all measured traits, except for the number of unfilled pods, seed eccentricity, and Feret's diameter.

Similarly, genotype × treatment interaction was highly significant (p < 0.001) for all investigated traits, excluding seed width, area, perimeter, equivalent diameter, and rugosity.

Heat Stress Effect on Grain Yield and Yield-Related Traits

Grain yield decreased substantially (36.6%) by exposure to heat compared to the no stress plants, which had a mean of 2.3 g plant⁻¹. The biological yield was higher for the no stress plants; the range was 6.1–16.4 g, with a mean of 11.3 g. Whereas, biological yield in plants exposed to high temperature showed a marked reduction of 38.2% (Table 1). Under high temperature conditions, the total number of pods plant⁻¹ was particularly low compared with no stress conditions, with a mean of 61.1. Similarly, heat stress reduced the total number of filled pods plant⁻¹ by 30.9% and relatively increased the number of unfilled pods plant⁻¹.

TABLE 1

Table 1. Range and mean ± SD of grain yield, yield related traits, seed size and shape parameters under normal and heat stress conditions.

Heat Stress Effect on Seed Size and Shape Parameters

Heat application resulted in a slight increase in seed length, width, area, perimeter, and equivalent diameter; the mean was 5 mm, 4.7 mm, 18.5 mm², 18.8 mm, and 4.8 mm, respectively. However, the mean was 4.8 mm, 4.5 mm, 16.9 mm², 17.0 mm, and 4.6 mm in normal conditions, respectively. In contrast, heat stress has no impact on seed eccentricity, Feret's diameter, thickness, and rugosity. Compared to normal conditions, the hundred-seed weight decreased by 27.2% due to heat stress (Table 1).

Heat Stress Effect on Nutritonal Quality

The mean Fe concentration declined in genotypes exposed to heat stress compared to no stress conditions. The range varied from 6.0 to 8.8 mg/100 g, with a mean of 7.2 mg/100 g. Whereas, under no stress conditions, the mean was 8.6 mg/100 g. Comparatively, Zn concentration was negatively influenced in response to high temperature; the values range from 5.0 to 6.6 mg/100 g with a mean of 5.4 mg/100 g. As compared to no stress conditions, high temperature caused a reduction in protein concentration. The maximum value for protein under no stress conditions was 27.4 g/100 g, while in heat stress, it reduced to 24.3 g/100 g. High temperature induced an increase in total phytic acid, which varied from 0.4 to 1.4 g/100 g; the mean was 1.0 g/100 g. While, under no stress conditions, total phytic acid varied from 0.5 to 1.2 g/100 g, with a mean of 0.9 mg/g. As for cooking time, the values ranged from 5 to 14 min; the mean cooking time was 8.1 min under no stress conditions. In contrast, cooking time showed a significant reduction (28.3%) under high temperature treatment, with a mean of 5.8 min (Table 2).

TABLE 2

Table 2. Range and mean ± SD of seed quality traits and cooking time under normal and heat stress conditions.

Heat Stress Effect on Micronutrient Bioavailability

The mean of both phytic acid/micronutrients (Fe and Zn) molar ratios showed reduced value compared to heat stress conditions. The high temperature increased PA/Fe and Zn molar ratios the mean was 11.3. and 16.8 (Table 3). Under no stress conditions, all studied lentil genotypes had a phytic acid/iron molar ratio above the critical value of 14, that estimates a good Fe bioavailability. Whereas, ILL 224, ILL950, ILL 4881, ILL 5968, ILL 5416, ILL 4804, ILL 597, ILL 6281, Zaaria, ILL 5261, and ILL 6644 had a phytic acid/zinc molar ratio above the critical values of 15. On the other hand, under heat stress conditions, the phytic acid/iron molar ratio was above the threshold of 14 for ILL 5505 ILL 2230 and ILL 956. As for phytic acid/zinc ratio, there was a shift of 13 genotypes from moderately to low bioavailable Zn and two genotypes from low to moderately bioavailable Zn, while the rest maintain the same category as in the normal conditions (Table 3).

TABLE 3

Table 3. Iron and zinc bioavailability of tested genotypes under no stress and heat stress conditions.

Association Between Seed Yield, Seed Shape and Nutritional Quality Traits

The correlation of grain yield with total number of pods and filled pods was positive and highly significant in both environments. Grain yield showed a highly significant positive correlation with biological yield under no stress conditions (r = 0.62; p < 0.001), while under heat stress treatment, this correlation was positive but not significant (Supplementary Table S3). Another positive and highly significant correlation, in both no stress and heat stress treatments, was observed between total number of filled pods and grain yield (r = 0.79; p < 0.001) and (r = 0.66; p < 0.001), respectively. The correlation between seed shape features and both grain yield and 100-seed weight of was statistically non-significant under both no stress and heat stress conditions. However, seed geometric features showed highly significant correlations among them, these correlations followed the same pattern under both environments with the exception of few parameters. Seed circularity showed a significant negative correlation with seed length (r = −0.37; p < 0.05), width (r = −0. 43; p < 0.01), area (r = −0. 37; p < 0.05), perimeter (r = −0. 40; p < 0.05) and equivalent diameter (r = −0. 41; p < 0.05) under high temperature treatment, whereas, under normal conditions these correlations were negative but remained non-significant. Under no stress conditions, seed thickness was negatively correlated with seed circularity (r = −0.64; p < 0.001) and Feret's diameter (r = −0. 36; p < 0.05), while under heat stress conditions seed thickness showed a significant negative correlation only with seed circularity (r = −0. 43; p < 0.01). On the other hand, all seed shape parameters showed no correlation with hundred-seed weight (Supplementary Table S3).

The correlations between grain yield and hundred-seed weight with different quality traits were also computed under both no stress and heat stress conditions (Figure 1). Interestingly, almost all the correlations were found to be statistically non-significant. There was no significant correlation between seed geometric parameters and seed quality traits under high temperature treatment, while under normal conditions, cooking time was found to be negatively correlated with seed circularity (r = −0.33; p < 0.05) and positively correlated with seed thickness (r = 0.35; p < 0.05). Among the nutritional quality traits, a significant positive association was identified between zinc and iron content under both no stress and heat stress conditions (r = 0.56 and r = 0.69; p < 0.001, respectively). In contrast, a significant negative correlation was detected between zinc and protein under both conditions. The correlation between iron and protein was negative and significant under heat stress treatment (r = −0.45; p < 0.01), whereas iron content did not appear to be significantly correlated with protein content under normal conditions. A highly significant negative association was obtained between phytic acid and cooking time under no stress conditions (r = 0.82, p < 0.001), while under heat stress conditions the association was negative but not significant. Cooking time and phytic acid showed no significant correlation with grain yield and other quality traits under normal and high temperature conditions. When looking at phytic acid/micronutrients ratio, PA/Zn showed a highly positive and significant correlation between cooking time (r = 0.78; p < 0.001), phytic acid (r = 0.91; p < 0.001) and PA/Fe ratio (r = 0.95; p < 0.001) under no stress conditions. On the other hand, under heat stress conditions, this ratio showed a negative and significant correlation with cooking time (r = −0.52; p < 0.01) and a highly positive significant correlation with phytic acid (r = 0.90; p < 0.001) and PA/Fe (r = 0.94; p < 0.001).

FIGURE 1

Figure 1. Pearson correlations between grain yield and seed qaulity traits. The diagonal panels show histograms for each trait. The lower and upper triangular panels, respectively, show scatter plots and Pearson correlation coefficients between grain yield and seed quality traits under no stress (A) and heat stress (B) conditions. GY, Grain yield (g); HSW, hundred-seed weight (g); ST, Seed thickness; CT, Cooking time (min); CP, Crude protein (g/100 g); PA, Phytic acid (g/100 g); Zn, Zinc content (mg/100 g); PA/Zn, Phytic acid/Zinc ratio; Fe, Iron content (mg/100 g); PA/Fe, Phytic acid/Iron ratio. *, **, and *** indicate significance at 0.05, 0.01, and 0.001 probability levels, respectively.

With respect to PA/Fe ratio, a highly positive and significant correlation was obtained with cooking time (r = 0.77; p < 0.001) and phytic acid (r = 0.91; p < 0.001) under no stress conditions. Whereas, under high temperature conditions, PA/Fe ratio revealed a highly significant positive association with phytic acid (r = 0.88; p < 0.001) and a significant negative correlation with cooking time (r = −0.51; p < 0.01). In contrast, PA/Fe ratio showed significant negative correlation with iron content under both conditions.

Principal Component Analysis

Under no stress conditions, the first two principal components explained 66.8% of the total variability (Table 4). PC1 and PC2 accounted for 42.4 and 24.5% of the total variation, while PC3 explained 16.2% of the variation. Cooking time, phytic acid and phytic acid/iron and zinc ratios showed the highest contribution to the first principal component (PC1). Crude protein, iron and zinc contents contributed the most for PC2, whereas grain yield and hundred-seed weight were the most contributing traits to PC3 (Figure 2). On the other hand, under heat stress conditions, the first two axes of the PC explained 56.9% of the total variation. PC1 accounted for 35.4% of the total variation, and PC2 and PC3 explained 21.5 and 56.9% of the total variation. The PC1 variation was mainly due to phytic acid and phytic acid/ iron and zinc ratios. The variation in PC2 was associated with grain yield, hundred-seed weight, and zinc content, while the variation in PC3 was due to iron and crude protein contents.

TABLE 4

Table 4. Percentage of contribution of different traits to the three major principal components with percentage variation under the different treatments.

FIGURE 2

Figure 2. Principal component analysis, hierarchical clustering and various traits contributing to the variability under no stress (A,B) and heat stress (C,D) conditions. GY, Grain yield (g); HSW, hundred-seed weight (g); CT, Cooking time (min); CP, Crude protein (g/100 g); PA, Phytic acid (g/100 g); Zn, Zinc content (mg/100 g); PA/Zn, Phytic acid/Zinc ratio; Fe, Iron content (mg/100 g); PA/Fe, Phytic acid/Iron ratio.

Hierarchical cluster analysis clearly classified tested genotypes into three distinct groups; the mean values of each group are listed in Table 5. Under no stress conditions, group_1 was characterized by high grain yield, hundred-seed weight and iron and zinc contents; the mean was 2.6 g, 1.5 g and 9.2 mg/100 g, respectively. Group_1 was also characterized by a short cooking time and a low phytic acid. The second group contained genotypes with high grain yield, hundred-seed weight and protein content, the mean was 2.5 g, 1.3 g and 25.6 g/100 g, respectively. Genotypes belonging to group_3 had a high level of phytic acid, long cooking time and a low bioavailable zinc. Under heat stress conditions, cluster analysis also revealed the existance of three groups. Genotypes falling into group_1 showed high grain yield and phytic acid and long cooking time,While group_2 was characteriezed by high protein content and low bioavailable zinc content. Group_3 was distinguished by a very high phytic acid level and a very low bioavailable iron and zinc content (Table 5).

TABLE 5

Table 5. Mean ± SD of seed quality traits, cooking time and phytic acid/ micronutrient ratio of the three clusters under no stress and heat stress conditions.

Discussion

Our results indicated that temperature beyond 30°C has a detrimental impact on seed yield and nutritional quality of investigated lentil genotypes. All genotypes were to some extent adversely affected by heat stress regarding biological yield, grain yield, and hundred-seed weight, which was in accordance with our previous study where heat stress decreased biological yield, grain yield and hundred-seed weight by 42.2, 53.5, and 26.8%, respectively (19). Our observations in this regard were similar to those in many studies. Sehgal et al. (26) revealed that application of heat stress during the reproductive stage reduced seed weight by 65–68% in two heat-tolerant genotypes. Supporting results were recently published, where terminal heat at the reproductive stage adversely affected both biological yield and grain yield but not hundred-seed weight (18). Similarly, heat stress during the grain filling period significantly reduced seed yield of chickpea genotypes (27) and common bean (28). Reduction in grain yield can be explained by the fact that high temperature during the reproductive stage negatively affects flower production (29). Failure of pod set can be explained by the negative impact of high temperature on pollen germination, pollen viability, pollen tube growth through the style (30, 31). In addition, our results indicated that the number of filled pod was adversely reduced under high temperature treatment, which is in good agreement with the above mentioned findings.

Seed size determination in lentils has historically relied on measuring the hundred or thousand seed weight (32). However, this method cannot give an indication of the actual seed shape parameters. Therefore, we used a computer-assisted image analysis method to determine both seed weight and seed shape. Our results showed that high temperature treatment induced slight changes in seed shape parameters. Seed length, width, perimeter and equivalent diameter showed a slight increase, while seed area was noted as moderately influenced character under heat stress. On the other hand, seed eccentricity, Feret's diameter, circularity and rugosity were not affected by heat stress. Little information is available in this regard, however, studies on other crops have revealed that seed shape related traits, such as seed width, seed length, and seed weight, are complex quantitative traits governed by polygenes and highly influenced by environmental conditions (33). Eşref and Nazmi (34) investigated the effect of different levels of moisture content on certain physical properties of the yellow lentil, results from this study showed that an increase in moisture content level led to an increase in seed area, diameter, thickness and diameter. On the other hand, correlation analysis in our study revealed no significant association between hundred-seed weight and seed shape parameters under both stress and no stress conditions. In contrast to our findings, Saha et al. (35) demonstrated that the hundred-seed weight of lentil had a higher positive correlation with seed diameter (r = 0.88; p < 0.01). Wang et al. (36) showed that the hundred-seed weight showed a high significant correlation with seed length, width and thickness in soybean, whereas, no correlation was detected between the hundred-seed weight and length/width ratio, length/thickness ratio, and width/thickness ratio.

The effect of high temperature on seed quality during reproductive stage was observed in many crops. Our previous findings have also highlighted the impact of heat stress on iron, zinc and crude protein (19). Similarly, in the present study, high temperature resulted in a significant reduction in both micronutrients. It has been reported that a moderately high temperature reduced the rate of nutrient uptake by roots in many crops (37). Therefore, increasing temperature due to global warming will likely have overall negative effects on plant nutrient relations, which will surely contribute to decreases in micronutrients level of lentil seed. Supporting results reported that high temperature resulted in a substantial decline in iron and zinc contents in lentil (38). High temperature also had a clear effect on protein content of tested genotypes. Recently, El Haddad et al. (39) reported that high temperature significantly decreased protein content, iron and zinc concentration in lentil by 15, 14, and 15%, respectively. Reduction in the rate of protein content accumulation and storage in legume seeds can be explained by inhibition of protein synthesis processes due to heat stress (29). Supporting studies have shown similar results in lentil (38), chickpea (40), soybean (41) and bean (42).

Due to the presence of a large genetic variability coupled with high heritability estimates (43, 44), lentil holds a great potential toward genetic biofortification for increased Fe and Zn content. However, many studies reported that iron and zinc absorption from diets based on legumes was found to be low (45). A major constraint to iron and zinc intake is the presence of antinutrients like phytic acid (46). In this study, tested genotypes were evaluated for their level of phytic acid under normal and high temperature conditions. Our results in this regard showed significant variation for phytic acid among studied genotypes and between treatments. Under no stress conditions, the variation for phytic acid ranged from 0.5 to 1.2 g/100 g with a mean of 0.9 g/100 g. Thavarajah et al. (20) reported a smaller range of variation (0.3–0.4 g/100 g) for phytic acid in lentil genotypes from three different market classes grown at two locations in Saskatchewan, Canada. Whereas, the Spanish brown cultivar Pardina grown in Spain, had higher level of total seed phytic acid concentration (1.2/100 g), that exceeded the mean value of phytic acid in our study. Similarly, Paul and Shang (47) found more similar mean value in 17 lentil genotypes grown in Bangladesh (0.7/100 g). On the other hand, high temperature treatment resulted in an increase of phytic acid level of tested genotypes. In agreement with our findings, Thavarajah et al. (48) indicated that phytic acid concentrations in lentil seeds were significantly higher in the rising temperature regime. Similar results were also found in common bean (49). Accumulation of phytic acid under stress conditions is likely related to its function in limiting oxidative stress in legumes (50). Accordingly, phytic acid has been reported to accumulate in chickpea seeds in response to drought stress (51).

To estimate the potential bioavailability of iron and zinc in lentil seeds, phytic acid/micronutrients molar ratios were calculated. It is considered that a higher molar ratio can decline mineral bioavailability and vice versa. According to Gibson (52), phytic acid/zinc molar ratios <5, between 5 and 15, and >15 have been associated with high, moderate and low zinc bioavailability. Whereas, phytic acid/iron molar ratio >14 can be a cause of iron deficiency (53, 54). In the present study, the calculated PA/Fe mean molar ratio under normal conditions was 8.6 with all genotypes below 14. This suggests that iron absorption might not be significantly impaired by phytic acid present in these genotypes. However, the mean PA/Fe molar ratio increased under high temperature treatment. Genotypes ILL 2230, ILL 5505, and ILL 956 showed a ratio above the critical value (13), and likely have a low bioavailable iron in their seeds. Regarding phytic acid/zinc ratio, twenty-five of the tested genotypes have a ratio below 15 and fell into moderate zinc bioavailability category, while other genotypes were low in bioavailable zinc. On the other hand, the accumulation of phytic acid due to rising temperature resulted in a shift of the PA/Zn molar ratio of the majority of genotypes from moderate to low bioavailability category. Collectively, these results showed that zinc content is the most affected micronutrient in terms of bioavailability, which can be explained by the fact that zinc forms the most stable complex with phytate contributing to a low bioavailability of zinc (55). Phytic acid/micronutients ratios in our study were less than those previously recorded levels in chickpea (56), faba bean (57), dry pea and kidney bean as well as in common bean, black gram and mung bean (58).

Our study showed a high genotypic variability for cooking time under both conditions. The mean cooking time under no stress conditions was 8 min, which was less than those previously reported by other authors (59, 60), while, the mean cooking time of genotypes under high temperature treatment was lower than the normal conditions. There is no available information related to the impact of heat stress on cooking time of lentil. However, it has been documented that cooking time is affected by many factors including climatic conditions, harvesting processes and seed storage conditions (61, 62). Furthermore, cooking time is mainly influenced by mineral nutrition, an adequate soil supply of phosphorus is essential to enhance lentil cooking quality (63). Correlation analysis in our study showed that cooking time was positively correlated with phytic acid under normal conditions. Supporting results were previously published, where high seed phytic acid content was related to high cookability in lentil (64, 65). It is noteworthy to mention that a decrease in phytic acid resulted in hard-to-cook phenomenon of many legume crops, including lentil. Galiotou-Panayotou et al. (66) stated that a decrease in phytic acid from 10.4 to 7.2 mg g⁻¹ induced a 42% increase in hardness. This suggests that low cookability of lentil seed grown under high temperature was due to the high level of phytic acid.Cooking time showed no significant correlation with 100-seed weight under both conditions, which was in accordance with finding reported by Ninou et al. (67). However, significant positive correlation has been reported in lentil (68) and common bean (69). Correlation analysis also showed that cooking time was positively correlated with seed thickness under normal conditions. In contrast, a non-significant correlation between cooking time and seed thickness has been demonstrated in lentil (67). Moreover, this seed shape parameter also showed non-significant correlation with cooking time under high temperature treatment. It should be pointed out that seed thickness has been reported to be largely influenced by environmental conditions during seed development (70). However, in our study seed thickness showed no reduction under high temperature treatment. From the very few studies reported here on factors influencing cooking time in lentil, our results suggests that a decrease in cooking time is connected with an increase in seed phytic acid level rather than its shape parameters.

In this study, Principal Component Analysis was used to define important traits contributing to the total variation in 36 lentil genotypes. The variation was mainly explained by the nutritional quality traits and grain yield, which were in accordance with our previous findings (19). In addition, Hierarchical Cluster Analysis was used to classify these genotypes based on similar nutritional quality traits and grain yield under both normal and stress conditions. Among the three identified groups, group_1 had the highest grain yield, hundred seed weight and micronutrients level. This group was also characterized by a highly bioavailable iron and zinc, a short cooking time and a low phytic acid. Genotypes belonging to this group could be useful for biofortification programs. Under heat stress conditions, genotypes belonging to group_1 showed high grain yield and phytic acid and long cooking time. These genotypes were able to defend themselves from heat stress by increasing the phytic acid level in their seeds. The use of multivariate analysis coupled with hierarchical cluster analysis provides an effective means for identifying the components that determine the variation when considering several traits simultaneously, and for grouping genotypes with similar traits. However, these findings can be further verified in a genotype x environment study to assess their nutritional quality and yield performance across different environmental conditions.

Conclusion

The current study has shown differences in micronutrients (Fe and Zn), protein content, phytic acid and cooking time, along with the evaluation of their responses to high temperature treatment. Our findings revealed that high temeprature will not only decrease the nutritional quality of lentil seed, but will also incease antinutrients like phytic acid, leading to a decrease in cooking time of lentil. In this regard, we have focused on the effect of phytic acid because of its major influence on iron and zinc bioavailability. Our results showed that the phytic acid/micronutrients molar ratios were somehow smaller under no stress conditions, while these ratios were higher under heat stress especially for phytic acid/zinc ratio, indicating that the estimated bioavailability of these micronutrients is low and is mainly affected by the phytic acid. However, a total of 11 promising genotypes have been identifid under high temperature conditions, namely ILL 5968, ILL 4841, Chakkouf, ILL 5261, Bichette, ILL 4881, Zaaria, ILL 494, ILL 2181, ILL 6493 and ILL 597. These genotypes were characteriased by high yielding short cooking time, high iron and moderate zinc bioavailability. Our findings could be useful for generating varieties with appropriate levels of phytic acid and micronutrients, which can contribute in mitigating micronutrients malnutrition.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Author Contributions

HC performed the research work and prepared first draft of the manuscript. HC, KA, and KH raised experiments in the greenhouse and recorded phenotyping data. HC and AE-B analyzed the seed samples for grain micronutrients in cereal and legumes quality laboratory. HC and NE performed statistical analysis of data. SK planned and supervised research activity. 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.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2022.857469/full#supplementary-material

References

1. Cubero JI, Perez De La Vega M, Fratini R. Origin, phylogeny, domestication and spread. Lentil Bot Prod Uses. (2009) 13−33. doi: 10.1079/9781845934873.0013

CrossRef Full Text | Google Scholar

2. FAOSTAT. Food Agric. Organ. Corp. Stat. Database. Available online at: http://www.fao.org/ (accessed January 18, 2021).

Google Scholar

3. Akibode CS, Maredia MK. Global and Regional Trends in Production, Trade and Consumption of Food Legume Crops (2012). 89. doi: 10.22004/ag.econ.136293

CrossRef Full Text | Google Scholar

4. Thavarajah D, Thavarajah P, Wejesuriya A, Rutzke M, Glahn RP, Combs GF, et al. The potential of lentil (Lens culinaris L.) as a whole food for increased selenium, iron, and zinc intake: preliminary results from a 3 year study. Euphytica. (2011) 180:123–8. doi: 10.1007/s10681-011-0365-6

CrossRef Full Text | Google Scholar

5. Urbano G, Porres JM, Frias J, Vidal-Valverde C. Nutritional value. In: Yadav, SS Mc Neil D, Stevenson PC, editors. Lentil An Anc. Crop Mod. Times. Springer, Netherlands (2007). p. 47–93.

Google Scholar

6. Sen Gupta D, Thavarajah D, Knutson P, Thavarajah P, McGee RJ, Coyne CJ, et al. Lentils (Lens culinaris L.), a rich source of folates. J Agric Food Chem. (2013) 61:7794–9. doi: 10.1021/jf401891p

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Graham C, Thavarajah D, Beck R. Carbohydrate concentration in lentils (Lens culinaris Medikus.): genotypic and environmental effects. Commun Soil Sci Plant Anal. (2017) 48:2447–54. doi: 10.1080/00103624.2017.1411937

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Johnson N, Johnson CR, 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–9. doi: 10.1002/ppp3.10103

CrossRef Full Text | Google Scholar

9. Ganesan K, Xu B. Polyphenol-rich lentils and their health promoting effects. Int J Mol Sci. (2017) 18:2390. doi: 10.3390/ijms18112390

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Tang J, Sokhansanj S, Sosulski FW. Moisture-absorption characteristics of laird lentils and hardshell seeds. Cereal Chem. (1994) 71:423–7.

Google Scholar

11. Wang N, Daun JK. Effects of variety and crude protein content on nutrients and anti-nutrients in lentils (Lens culinaris). Food Chem. (2006) 95:493–502. doi: 10.1016/j.foodchem.2005.02.001

CrossRef Full Text | Google Scholar

12. Urbano G, Lopez-Jurado M, Aranda P, Vidal-Valverde C, Tenorio E, Porres J. The role of phytic acid in legumes: antinutrient or beneficial function? J Physiol Biochem. (2000) 56:283–94. doi: 10.1007/BF03179796

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Majumder S, Datta K, Datta SK. Rice biofortification: high Iron, Zinc, and vitamin-A to fight against “hidden hunger.” Agronomy. (2019) 9:803. doi: 10.3390/agronomy9120803

CrossRef Full Text | Google Scholar

14. Hunt JR. Bioavailability of iron, zinc, and other trace minerals from vegetarian diets. Am J Clin Nutr. (2003) 78:633S−9S. doi: 10.1093/ajcn/78.3.633S

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Wakeel A, Farooq M, Bashir K, Ozturk L. Micronutrient malnutrition and biofortification: recent advances and future perspectives. Plant Micronutr. Use Effic. (2018) 225–43. doi: 10.1016/B978-0-12-812104-7.00017-4

CrossRef Full Text | Google Scholar

16. Vasić M, Tepić A, Mihailović V, Mikić A, Gvozdanović-Varga J, Šumić Z, et al. Phytic acid content in different dry bean and faba bean landraces and cultivars. Rom Agric Res. (2012) 29:79–85.

Google Scholar

17. Singh J, Kanaujia R, Srivastava AK, Dixit GP, Singh NP. Genetic variability for iron and zinc as well as antinutrients affecting bioavailability in black gram (Vigna mungo (L.) Hepper). J Food Sci Technol. (2017) 54:1035–42. doi: 10.1007/s13197-017-2548-1

PubMed Abstract | CrossRef Full Text | Google Scholar

18. El Haddad N, Rajendran K, Abdelaziz S, Es-Safi NE, Nadia B, Mentag R, et al. Screening the FIGS set of lentil (Lens culinaris Medikus) germplasm for tolerance to terminal heat and combined drought-heat stress. Agronomy. (2020) 10:1036. doi: 10.3390/agronomy10071036

CrossRef Full Text | Google Scholar

19. Choukri H, Hejjaoui K, El-Baouchi A. Heat and drought stress impact on phenology, grain yield, and nutritional quality of lentil (Lens culinaris Medikus). Front Nutr. (2020) 7:596307. doi: 10.3389/fnut.2020.596307

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Thavarajah P, Thavarajah D, Vandenberg A. Low phytic acid lentils (Lens culinaris L.): a potential solution for increased micronutrient bioavailability. J Agric Food Chem. (2009) 57:9044–9. doi: 10.1021/jf901636p

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Baethgen WE, Alley MM. A manual colorimetric procedure for measuring ammonium nitrogen in soil and plant Kjeldahl digests. Commun Soil Sci Plant Anal. (1989) 20:961–9. doi: 10.1080/00103628909368129

CrossRef Full Text | Google Scholar

22. Megazyme Phytic Acid Assay Kit. Available online at: https://www.megazyme.com/phytic-acid-assay-kit (accessed January 1, 2021).

Google Scholar

23. Extrac Visualize the Results of Multivariate Data Analyses. R Package Version 1.0.3. Available online at: https://CRAN.R-project.org/package=factoextra (accessed January 25, 2021).

Google Scholar

24. Lê S, Josse J, Husson F. FactoMineR: an R package for multivariate analysis. J Stat Softw. (2008) 25:1–18. doi: 10.18637/jss.v025.i01

CrossRef Full Text | Google Scholar

25. Galili T. dendextend: an R package for visualizing, adjusting and comparing trees of hierarchical clustering. Bioinformatics. (2015) 31:3718–20. doi: 10.1093/bioinformatics/btv428

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Sehgal A, Sita K, Kumar J, Kumar S, Singh S, Siddique KHM, et al. 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. doi: 10.3389/fpls.2017.01776

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Awasthi R, Kaushal N, Vadez V, Turner NC, Berger J, Siddique KHM, et al. Individual and combined effects of transient drought and heat stress on carbon assimilation and seed filling in chickpea. Funct Plant Biol. (2014) 41:1148–67. doi: 10.1071/FP13340

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Suárez JC, Polanía JA, Contreras AT, Rodríguez L, Machado L, Ordoñez C, et al. Adaptation of common bean lines to high temperature conditions: genotypic differences in phenological and agronomic performance. Euphytica. (2020) 216:1–20. doi: 10.1007/s10681-020-2565-4

CrossRef Full Text | Google Scholar

29. Bhandari K, Sharma KD, Rao BH, Siddique KHM, Gaur P, Agrawal SK, et al. Temperature sensitivity of food legumes: a physiological insight. Acta Physiol Plant. (2017) 39:68. doi: 10.1007/s11738-017-2361-5

CrossRef Full Text | Google Scholar

30. Kumar S, Thakur P, Kaushal N, Malik JA, Gaur P, Nayyar H. Effect of varying high temperatures during reproductive growth on reproductive function, oxidative stress and seed yield in chickpea genotypes differing in heat sensitivity. Arch Agron Soil Sci. (2013) 59:823–43. doi: 10.1080/03650340.2012.683424

CrossRef Full Text | Google Scholar

31. Sita K, Sehgal A, Kumar J, Kumar S, Singh S, Siddique KHM, et al. Identification of high-temperature tolerant lentil (Lens culinaris Medik.) genotypes through leaf and pollen traits. Front. Plant Sci. (2017) 8:744. doi: 10.3389/fpls.2017.00744

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Tullu A, Kusmenoglu I, McPhee KE, Muehlbauer FJ. Characterization of core collection of lentil germplasm for phenology, morphology, seed and straw yields. Genet Resour Crop Evol. (2001) 48:143–52. doi: 10.1023/A:1011254629628

CrossRef Full Text | Google Scholar

33. Breseghello F, Sorrells ME. QTL analysis of kernel size and shape in two hexaploid wheat mapping populations. Field Crop Res. (2007) 101:172–9. doi: 10.1016/j.fcr.2006.11.008

CrossRef Full Text | Google Scholar

34. Eşref I, Nazmi I. Effects of moisture content on some physical properties of the yellow lentil. J Agric Sci. (2016) 22:307–16. doi: 10.1501/Tarimbil_0000001389

CrossRef Full Text | Google Scholar

35. Saha GC, Sarker A, Chen W, Vandemark GJ, Muehlbauer FJ. Inheritance and linkage map positions of genes conferring resistance to stemphylium blight in lentil. Crop Sci. (2010) 50:1831–9. doi: 10.2135/cropsci2009.12.0709

CrossRef Full Text | Google Scholar

36. Wang X, Fu G, Yuan S, Zhang H. Genetic diversity analysis of seed appearance quality of Chinese soybean mini core collection. Procedia Eng. (2011) 18:392–7. doi: 10.1016/j.proeng.2011.11.063

CrossRef Full Text | Google Scholar

37. Heckathorn SA, Giri A, Mishra S, Bista D. Heat stress and roots. Clim. Chang. Plant Abiotic Stress Toler. (2013) 1:109–36. doi: 10.1002/9783527675265.ch05

CrossRef Full Text | Google Scholar

38. Sita K, Sehgal A, Bhandari K, Kumar J, Kumar S, Singh S, et al. 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–41. doi: 10.1002/jsfa.9054

PubMed Abstract | CrossRef Full Text | Google Scholar

39. El Haddad N, Choukri H, Ghanem ME, Smouni A, Mentag R, Rajendran K, et al. 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. doi: 10.3390/plants11010095

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Behboudian MH, Ma Q, Turner NC, Palta JA. Reactions of chickpea to water stress: yield and seed composition. J Sci Food Agric. (2001) 81:1288–91. doi: 10.1002/jsfa.939

CrossRef Full Text | Google Scholar

41. Krishnan HB, Kim W-S, Oehrle NW, Smith JR, Gillman JD. Effect of heat stress on seed protein composition and ultrastructure of protein storage vacuoles in the cotyledonary parenchyma cells of soybean genotypes that are either tolerant or sensitive to elevated temperatures. Int J Mol Sci. (2020) 21:4775. doi: 10.3390/ijms21134775

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Ghanbari AA, Shakiba MR, 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–12.

Google Scholar

43. Sen Gupta D, Thavarajah D, McGee RJ, Coyne CJ, Kumar S, Thavarajah P. Genetic diversity among cultivated and wild lentils for iron, zinc, copper, calcium and magnesium concentrations. Aust J Crop Sci. (2016) 10:1381–7. doi: 10.21475/ajcs.2016.10.10.pne6

CrossRef Full Text | Google Scholar

44. Kumar J, Thavarajah D, Kumar S, Sarker A, Singh NP. Analysis of genetic variability and genotype × environment interactions for iron and zinc content among diverse genotypes of lentil. J Food Sci Technol. (2018) 55:3592–605. doi: 10.1007/s13197-018-3285-9

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Sandberg A-S. Bioavailability of minerals in legumes. Br J Nutr. (2002) 88:281–5. doi: 10.1079/BJN/2002718

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Raboy V. Progress in breeding low phytate crops. J Nutr. (2002) 132:503S−5S. doi: 10.1093/jn/132.3.503S

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Paul DK, Shang RK. Nutritive values and phytate content of some new varieties of wheat and legumes. Pak J Biol Sci. (2004) 7:1006–9. doi: 10.3923/pjbs.2004.1006.1009

CrossRef Full Text | Google Scholar

48. Thavarajah P, See C-T, Vandenberg A. Phytic acid and Fe and Zn concentration in lentil (Lens culinaris L.) seeds is influenced by temperature during seed filling period. Food Chem. (2010) 122:254–9. doi: 10.1016/j.foodchem.2010.02.073

CrossRef Full Text | Google Scholar

49. Hummel M, Hallahan BF, Brychkova G, Ramirez-Villegas J, Guwela V, Chataika B, 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:1–11. doi: 10.1038/s41598-018-33952-4

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Kido EA, Neto JR, Silva RL, Belarmino LC, Neto JP, Soares-Cavalcanti NM, et al. Expression dynamics and genome distribution of osmoprotectants in soybean: identifying important components to face abiotic stress. BMC Bioinform. (2013) 14:1–11. doi: 10.1186/1471-2105-14-S1-S7

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Joshi-Saha A, Reddy KS. Repeat length variation in the 5′UTR of myo-inositol monophosphatase gene is related to phytic acid content and contributes to drought tolerance in chickpea (Cicer arietinum L.). J Exp Bot. (2015) 66:5683–90. doi: 10.1093/jxb/erv156

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Gibson RS. Zinc: the missing link in combating micronutrient malnutrition in developing countries. Proc Nutr Soc. (2006) 65:51–60. doi: 10.1079/PNS2005474

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Lestienne I, Mouquet-Rivier C, Icard-Vernière C, Rochette I, Treche S. The effects of soaking of whole, dehulled and ground millet and soybean seeds on phytate degradation and Phy/Fe and Phy/Zn molar ratios. Int J Food Sci Technol. (2005) 40:391–9. doi: 10.1111/j.1365-2621.2004.00941.x

CrossRef Full Text | Google Scholar

54. Luo Y, Gu Z, Han Y, Chen Z. The impact of processing on phytic acid, in vitro soluble iron and Phy/Fe molar ratio of faba bean (Vicia faba L.). J Sci Food Agric. (2009) 89:861–6. doi: 10.1002/jsfa.3525

CrossRef Full Text | Google Scholar

55. Weaver CM, Kannan S. Phytate and mineral bioavailability. In Rukma Reddy N, Sathe SK, editors. Food Phytates. Boca Raton, FL: CRC Press (2001). p. 227–40.

Google Scholar

56. Kaya M, Küçükyumuk Z, Erdal I. Phytase activity, phytic acid, zinc, phosphorus protein contents in different chickpea genotypes in relation to nitrogen zinc fertilization. Afr J Biotechnol. (2009) 8.

Google Scholar

57. Lazarte CE, Carlsson N.-G., Almgren A, Sandberg A.-S., et al. Phytate, zinc, iron and calcium content of common Bolivian food, and implications for mineral bioavailability. J. Food Compos. Anal. (2015) 39:111–9. doi: 10.1016/j.jfca.2014.11.015

CrossRef Full Text | Google Scholar

58. Seneviratne M, Gunasinghe WKR, Gama-Arachchige NS, Dissanayake N, Karunaratne A. Mineral bioavailability in three locally consumed pulses processed using popular methods: interpreted using molar ratios with phytic acid. Ceylon J Sci. (2012) 41:19–26. doi: 10.4038/cjsbs.v41i1.4534

CrossRef Full Text | Google Scholar

59. Wang N, Daun JK. Determination of cooking times of pulses using an automated Mattson cooker apparatus. J Sci Food Agric. (2005) 85:1631–5. doi: 10.1002/jsfa.2134

CrossRef Full Text | Google Scholar

60. Pirhayati M, Soltanizadeh N, Kadivar M. Chemical and microstructural evaluation of ‘hard-to-cook'phenomenon in legumes (pinto bean and small-type lentil). Int J Food Sci Technol. (2011) 46:1884–90. doi: 10.1111/j.1365-2621.2011.02697.x

CrossRef Full Text | Google Scholar

61. Iliadis C. Effects of harvesting procedure, storage time and climatic conditions on cooking time of lentils (Lens culinaris Medikus). J Sci Food Agric. (2001) 81:590–3. doi: 10.1002/jsfa.848

CrossRef Full Text | Google Scholar

62. Yousif AM, Kato J, Deeth HC. Effect of storage on the biochemical structure and processing quality of adzuki bean (Vigna angularis). Food Rev Int. (2007) 23:1–33. doi: 10.1080/87559120600865172

CrossRef Full Text | Google Scholar

63. Hussain M, Shah SH, Nazir MS. Qualitative response of three lentil (Lens culinaris Medic.) cultivars to phosphorus application. Pak J Agri Sci. (2003) 40:25–7.

Google Scholar

64. Bhatty RS. Comparisons of good-and poor-cooking lentils. J Sci Food Agric. (1995) 68:489–96. doi: 10.1002/jsfa.2740680413

CrossRef Full Text | Google Scholar

65. Liu K, Bourne MC. Cellular, biological, and physicochemical basis for the hard-to-cook defect in legume seeds. Crit Rev Food Sci Nutr. (1995) 35:263–98. doi: 10.1080/10408399509527702

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Galiotou-Panayotou M, Kyriakidis NB, Margaris I. Phytase–phytate–pectin hypothesis and quality of legumes cooked in calcium solutions. J Sci Food Agric. (2008) 88:355–61. doi: 10.1002/jsfa.2973

CrossRef Full Text | Google Scholar

67. Ninou E, Papathanasiou F, Vlachostergios DN, Mylonas I, Kargiotidou A, Pankou C, et al. Intense breeding within lentil landraces for high-yielding pure lines sustained the seed quality characteristics. Agriculture. (2019) 9:175. doi: 10.3390/agriculture9080175

CrossRef Full Text | Google Scholar

68. Erskine W, Williams PC, Nakkoul H. Genetic and environmental variation in the seed size, protein, yield, and cooking quality of lentils. Field Crop Res. (1985) 12:153–61. doi: 10.1016/0378-4290(85)90061-9

CrossRef Full Text | Google Scholar

69. Santos CS, Carbas B, Castanho A, Bronze MR, Serrano C, Vasconcelos MW, et al. Relationship between seed traits and pasting and cooking behaviour in a pulse germplasm collection. Crop Pasture Sci. (2018) 69:892–903. doi: 10.1071/CP18205

CrossRef Full Text | Google Scholar

70. Domoney C, Duc G, Ellis THN, Ferrándiz C, Firnhaber C, Gallardo K, et al. Genetic and genomic analysis of legume flowers and seeds. Curr Opin Plant Biol. (2006) 9:133–41. doi: 10.1016/j.pbi.2006.01.014

PubMed Abstract | CrossRef Full Text | Google Scholar