Abstract
Climate change affects the Mediterranean region stressing lentil crops during flowering and seed set. Early maturation and drought tolerance are desirable traits in these conditions. Phosphorus (P) is considered to enhance early flowering, maturity and thus yields. Four P rates (0, 30, 60, 90 kg P2O5 ha−1) were applied on four cultivars (Samos, Thessaly, Flip, Ikaria) during two seasons. Growing degree-days (GDD) were calculated for vegetative (V4–5, V7–8) and reproductive stages (R1, R2, R4, R6, R8). At R2 (full bloom) carbon isotope discrimination (Δ) was used to assess water-use efficiency. At R8 (full maturity), the seed weight (SW) was determined by harvest. Cultivars, P and the P × cultivar and P × growth season interactions affected the earliness in reproductive stages; P had no effect on GDD of vegetative stages. Phosphorus both induced earliness (Flip, Thessaly) and delayed maturity (Samos, Ikaria). GDD and SW were negatively correlated for the P × cultivar interaction at R1 (first bloom), R2, R4 (flat pod) and R6 (full pod filling) stages; being the strongest at R1. Negative correlations were evident for the P × growth season interaction at R2, R4 and R6 stages; being the strongest at R4. Cultivars and P did not affect Δ. A proper combination of cultivar and P rate can mitigate lentil yield losses under changing Mediterranean climate.
INTRODUCTION
Human activity contributes to increasing the concentration of greenhouse gases in the atmosphere thereby increasing the radiation absorbed, and subsequently affecting the air temperature. Climate change is a daunting environmental challenge because of its direct biophysical effects on plant physiological processes and consequently on crop production. In recent years, extreme climate events, including heat waves, heavy precipitation and drought incidents, have strongly affected Europe and led to crop failure (Reidsma et al. 2010).
Although higher temperatures in northern Europe are expected to increase crop yields, in the Mediterranean and south-western Balkans increasing temperatures, in combination with reduced rainfall, are expected to lead to large yield reductions (Olesen et al. 2011). However, it is reported that elevated levels of carbon dioxide (eCO2), due to climate change, can mitigate the effects of drought in grain legumes such as soybean (Li et al. 2013), as plants reduce stomatal and canopy conductance and therefore limit soil water absorption. While it is difficult to predict accurate scenarios regarding future climate change, crop productivity is a key factor affecting food security and economic stability (Hatfield et al. 2011).
Lentils (Lens culinaris Medik.) are grain legumes which have been grown around the Mediterranean basin since antiquity and still remain a structural block of the Mediterranean diet pyramid. Apart from their high nutritional value, legumes, and lentils in particular, constitute a food group which can effectively meet the demands imposed by climate change as the impact on the environment is minimal (Graham & Vance 2003). In the Mediterranean region, lentils are usually sown during the wet winter months, but decreasing soil moisture and increasing temperatures during the reproductive stages result in terminal drought and low seed yields (Shrestha et al. 2006). To cope with such phenomena, it is necessary to exploit genetic variation reported to reside among cultivated germplasm as well as within wild Lens spp. (Erskine et al. 2011).
Earliness in flowering is a desirable trait, especially under the semi-arid conditions prevailing in the Mediterranean region. Early sowing has been an effective strategy in many areas to avoid high temperatures and low rainfall during the reproductive stages, and thus to increase seed yields (Hajarpoor et al. 2014); however, earlier sowing was reported to increase susceptibility to disease due to longer favorable pre-winter weather conditions (Colbach et al. 1997). On the other hand, earliness owing to genotypic variation leading to drought avoidance could be exploitable. In addition, previous studies have reported the beneficial effect of increased phosphorus (P) rates on earliness in flowering and maturity. This effect is possibly caused by an alteration of the nitrogen/phosphorus (N/P) ratio in plant tissues and/or by the beneficial role of P in stimulation of the reproductive parts of plants (Sawan et al. 2008). Genotypic variation and environmental effects on water economy and productivity of C3 species can easily and reliably be assessed by carbon isotope discrimination (Δ values) in above-ground biomass. Δ values (a measure of the 13C/12C ratio in plant tissues compared to the air) are related to the ratio of intercellular to ambient CO2 concentrations (ci/ca) being a surrogate of the intrinsic water use efficiency [(WUEi, the ratio of CO2 assimilation rate to stomatal conductance (gs)] (Farquhar & Richards 1984). Moreover, Δ values have been proved a reliable, indirect and long-term indicator of water use efficiency at biomass level (WUE, the ratio of biomass produced to the water consumed to produce it) for many species; relationships between Δ and yield have also been established (Turner 1996).
In the Mediterranean region, soils are usually P-deficient, thus the use of P fertilizers remains relatively high (Harmsen & El Mahmoud 2004). Phosphorus is an essential macronutrient with a pivotal role in numerous plant processes. Since lentils meets their needs for N through N2-fixation, P is the most common limiting nutrient for both growth and grain yield (Yadav et al. 2009). Moreover, Jin et al. (2014) mentioned that eCO2 increased drought tolerance as well as water use efficiency (WUE) of field pea (Pisum sativum), especially when plants are grown under sufficient P supply.
Temperature is perhaps the most critical factor in all biological processes, therefore heat units are often used to estimate phenological development of plant species. Growing degree-days (GDD) provide a simple estimate of the accumulated heat available during the growth season or life cycle of an organism (Robertson et al. 2013). The concept of degree-days is based on three assumptions: (i) there is a basal temperature for all species below which no growth is observed, (ii) the growth of all kinds is proportional to the total accumulation of energy over a specific period, and (iii) species' maturity occurs only when a specific total of degree-days is achieved (Fealy & Fealy 2008).
Therefore, the aim of this work was to study the effect of four P rates (0, 30, 60 and 90 kg P2O5 ha−1) on the earliness of four lentil cultivars at seven growth stages, through the calculation of GDD. Furthermore, in order to elucidate whether P supply can affect WUE and seed yield under P-deficient Mediterranean soil, Δ values were determined in above-ground biomass at full bloom.
MATERIALS AND METHODS
Study area and plant material
A field experiment was conducted for two growth seasons (2013–2014, hereafter 2014 and 2014–2015, hereafter 2015) in the farm of the Aristotle University of Thessaloniki (AUTh; 40° 32′12 N, 22° 59′21 E, 6 m a.s.l.). The climate is typically Mediterranean. Mean monthly temperature and monthly precipitation during the growth seasons (December to May) are shown in Figure 1 while Table 1 presents soil characteristics of the experimental site.
Soil properties at 0–30 cm depth before the establishment of the experiments
. | 2014 . | 2015 . |
---|---|---|
Sand (g kg−1)a | 190 | 260 |
Silt (g kg−1)a | 360 | 300 |
Clay (g kg−1)a | 450 | 440 |
Texture | Clay | Clay |
pH (1:1 in H2O)b | 8.2 (0.1)c | 8.1 (0.1) |
Organic matter (g kg−1)d | 12.9 (1.8) | 11.6 (1.0) |
N-NO3 (mg kg−1)e | 37.5 (6.0) | 35.9 (3.3) |
P-Olsen (mg kg−1)f | 8.8 (1.4) | 8.4 (0.6) |
K (mg kg−1)g | 97 (0.9) | 112 (15.9) |
Na (mg kg−1)g | 65 (8.5) | 52 (8.1) |
Ca (mg kg−1)g | 2486 (39) | 2629 (60) |
Mg (mg kg−1)g | 798 (29) | 760 (18) |
CEC (cmol+ kg−1)h | 27 (1.0) | 27 (1.0) |
EC (dS m−1)i | 0.45 (0.02) | 0.46 (0.01) |
. | 2014 . | 2015 . |
---|---|---|
Sand (g kg−1)a | 190 | 260 |
Silt (g kg−1)a | 360 | 300 |
Clay (g kg−1)a | 450 | 440 |
Texture | Clay | Clay |
pH (1:1 in H2O)b | 8.2 (0.1)c | 8.1 (0.1) |
Organic matter (g kg−1)d | 12.9 (1.8) | 11.6 (1.0) |
N-NO3 (mg kg−1)e | 37.5 (6.0) | 35.9 (3.3) |
P-Olsen (mg kg−1)f | 8.8 (1.4) | 8.4 (0.6) |
K (mg kg−1)g | 97 (0.9) | 112 (15.9) |
Na (mg kg−1)g | 65 (8.5) | 52 (8.1) |
Ca (mg kg−1)g | 2486 (39) | 2629 (60) |
Mg (mg kg−1)g | 798 (29) | 760 (18) |
CEC (cmol+ kg−1)h | 27 (1.0) | 27 (1.0) |
EC (dS m−1)i | 0.45 (0.02) | 0.46 (0.01) |
The values presented are the average over three samples.
aHydrometer method.
bBy pH meter.
cStandard errors in parentheses.
dWet oxidation method.
eExtractable with KCl 2 M.
fOlsen method.
gExtractable with CH3COONH4.
hCH3COONa-CH3COONH4 method.
iIn saturated paste by conductance meter.
Monthly mean air temperature (°C) and precipitation (mm) during the two growth seasons.
Monthly mean air temperature (°C) and precipitation (mm) during the two growth seasons.
Under rainfed conditions, four lentil cultivars [Samos, Thessaly, Flip 2003-24 L (Flip), Ikaria], provided by the Hellenic Agriculture Organization (HAO)-DEMETER, Institute of Industrial and Fodder Plants, Larissa, Greece, were supplemented, before sowing, with four phosphorus (P) rates [0 (P0), 30 (P30), 60 (P60) and 90 kg P ha−1 (P90)] in a split-plot design with P rates in the main plots and cultivars in the subplots. The treatments were triplicated. Phosphorus was applied as triple superphosphate (460 g P2O5 kg−1). Seeding was conducted by hand on 6 December 2014 and 5 December 2015. Each subplot consisted of six rows, 4 m long and 0.25 m apart (6 m2), at a sowing density of 100–110 seeds m−2. Breeder's seeds were used for all the four cultivars.
For each subplot, the Julian date, by which specific growth stages were achieved, was recorded. Seven growth stages, following Erskine et al. (1990), were monitored: two vegetative, the stages of 4–5th and 7–8th leaves (V4–5 and V7–8, respectively), and five reproductive stages, namely first bloom (R1), full bloom (R2), flat pod stage (R4), full pod filling stage (R6) and full maturity (R8). The achievement of the growth stage was considered when at least 50% of the plants in each subplot had reached the defined stage.
Calculation of growing degree-days
Determination of carbon isotope discrimination and seed weight
At full maturity (R8 growth stage), a 0.125 m2 surface was harvested and threshed by hand in order to determine seed weight (SW) per subplot (in g m−2) after normalization at 10% seed moisture.
Statistical analyses
Data of GDD for each growth stage were subjected to analysis of variance (ANOVA) as an over-year, split-plot design with P rates in the main plots and cultivars in the subplots and with three replications. Means were compared by the least significant difference (LSD) test at P < 0.05 and analyses were run using the statistical software JMP 5.1 (SAS Institute Inc., Cary, NC, USA).
RESULTS AND DISCUSSION
In the present study, over the years, cv. Flip had the highest SW (406.34 g m−2) while it was the earliest among the three cultivars; it was followed by Thessaly (343.03 g m−2), Samos (299.36 g m−2) and finally Ikaria (283.59 g m−2), the latest maturing cultivar. The effects of climate change include high temperatures and prolonged drought periods in the Mediterranean region (Iglesias & Garrote 2015); thus earliness in flowering is a desirable trait under the prevailing hot and dry conditions. Moreover, lentil, as a cool season legume, is more sensitive to high temperatures, especially when they occur during the flowering and pod-filling stages resulting in high yield losses (Farooq et al. 2017). In most lentil producing areas, yield does not exceed half the potential yield, mostly due to scarcity of water combined with extreme temperatures, indicating the need to increase crop yield through genetic improvement and/or better crop management (Ghanem et al. 2015).
Breeding for stress tolerant cultivars is more intricate because of the complex nature of the mechanism that gives tolerance. Thus, a range of different crop management practices can be applied to mitigate the effects of drought stress including sowing date, early maturing cultivars and appropriate fertilizer use. In the present work, neither cultivars nor P rates had a significant effect on GDD of the vegetative stages but both factors affected the earliness of the reproductive stages; the interactions P rate × cultivar and P rate × growth season were significant as well. More specific, at the beginning of flowering (R1 growth stage), cv. Flip was the earliest showing no significant differences between P rates. On the contrary, in cvs. Samos, Thessaly and Ikaria, P additions induced earliness compared to the control (P0) treatments (Table 2). The highest GDD were recorded for P0 in cvs. Ikaria and Samos (1033.30 °C and 1042.15 °C, respectively). Similar results were evident at full bloom (R2 growth stage) when cv. Flip recorded the lowest GDD (998.89 °C) regardless of the P rate. At flat pod stage (R4), P additions began to stand out in comparison to P0. Treatments Ρ60 and Ρ90 induced earliness in Flip and Ikaria while in cv. Thessaly, all the P additions (P30, P60 and P90) scored lower GDD than the control (P0). In cv. Samos, only treatment P90 induced earliness (Table 2). At R6 stage, cv. Flip remained the earliest while in cv. Thessaly, all the P rates scored the same GDD. On the contrary, in cvs. Samos and Ikaria, P additions caused delayed maturity (Table 2). The stages and the duration required to mature a plant are determined by many factors such as climate, genetic material, agronomic practices and soil fertility (Miller et al. 2001). It is mentioned that P fertilization is more beneficial in legumes, as it regulates the N/P ratio in the soil resulting in vigorous growth and higher yields (Li et al. 2011). According to Rasheed et al. (2010), despite the fact that P application in lentils induces earliness in flowering, it can delay maturity due to increased nitrogenase activity of active nodules and due to modification of the ratio of P to other nutrients. For the P rate × cultivar interaction, significant negative correlations between GDD and SW were found at R1, R2, R4 and R6 stages, with that for R1 stage being the strongest (Figure 2).
Mean comparison of growing degree-days (GDD, °C) for the P rate × cultivar interaction for the four reproductive growth stages (R1, R2, R4 and R6) that showed significant differences
Interaction . | . | GDD (°C) . | |||
---|---|---|---|---|---|
. | R1 . | R2 . | R4 . | R6 . | |
Samos | Ρ0 | 1042.15a | 1168.50a | 1298.38a | 1443.52c |
Ρ30 | 1015.28b | 1162.53a | 1303.67a | 1466.52a | |
Ρ60 | 990.85c | 1139.11b | 1295.63a | 1466.52a | |
Ρ90 | 1004.08cb | 1139.11b | 1279.55b | 1466.52a | |
Thessaly | Ρ0 | 993.53c | 1114.49c | 1263.51c | 1429.17d |
Ρ30 | 971.37d | 1060.70d | 1235.07d | 1429.17d | |
Ρ60 | 969.21d | 1055.46d | 1235.07d | 1429.17d | |
Ρ90 | 969.21d | 1057.96d | 1235.07d | 1429.17d | |
Flip | Ρ0 | 921.70e | 998.89e | 1186.72e | 1335.85e |
Ρ30 | 914.80e | 998.89e | 1150.27f | 1335.85e | |
Ρ60 | 910.75e | 998.89e | 1140.17fg | 1327.53f | |
Ρ90 | 910.75e | 998.89e | 1129.85g | 1330.37f | |
Ikaria | Ρ0 | 1033.30a | 1147.78b | 1265.65c | 1443.52c |
Ρ30 | 993.60c | 1104.72c | 1265.65c | 1450.80b | |
Ρ60 | 990.85c | 1102.22c | 1242.45d | 1450.80b | |
Ρ90 | 990.85c | 1109.23c | 1242.45d | 1450.80b |
Interaction . | . | GDD (°C) . | |||
---|---|---|---|---|---|
. | R1 . | R2 . | R4 . | R6 . | |
Samos | Ρ0 | 1042.15a | 1168.50a | 1298.38a | 1443.52c |
Ρ30 | 1015.28b | 1162.53a | 1303.67a | 1466.52a | |
Ρ60 | 990.85c | 1139.11b | 1295.63a | 1466.52a | |
Ρ90 | 1004.08cb | 1139.11b | 1279.55b | 1466.52a | |
Thessaly | Ρ0 | 993.53c | 1114.49c | 1263.51c | 1429.17d |
Ρ30 | 971.37d | 1060.70d | 1235.07d | 1429.17d | |
Ρ60 | 969.21d | 1055.46d | 1235.07d | 1429.17d | |
Ρ90 | 969.21d | 1057.96d | 1235.07d | 1429.17d | |
Flip | Ρ0 | 921.70e | 998.89e | 1186.72e | 1335.85e |
Ρ30 | 914.80e | 998.89e | 1150.27f | 1335.85e | |
Ρ60 | 910.75e | 998.89e | 1140.17fg | 1327.53f | |
Ρ90 | 910.75e | 998.89e | 1129.85g | 1330.37f | |
Ikaria | Ρ0 | 1033.30a | 1147.78b | 1265.65c | 1443.52c |
Ρ30 | 993.60c | 1104.72c | 1265.65c | 1450.80b | |
Ρ60 | 990.85c | 1102.22c | 1242.45d | 1450.80b | |
Ρ90 | 990.85c | 1109.23c | 1242.45d | 1450.80b |
Within columns, means labeled with the same letter did not differ significantly at P < 0.05 using the LSD test.
Correlation between GDD and seed weight for the P rate × cultivar interaction at R1, R2, R4 and R6 growth stages, which showed the highest correlation coefficients. Where GDD, growing degree-days; R1, first bloom stage; R2, full bloom stage; R4, flat pod stage; R6, full pod filling stage.
Correlation between GDD and seed weight for the P rate × cultivar interaction at R1, R2, R4 and R6 growth stages, which showed the highest correlation coefficients. Where GDD, growing degree-days; R1, first bloom stage; R2, full bloom stage; R4, flat pod stage; R6, full pod filling stage.
For the P rate × growth season interaction, all reproductive stages but R8 showed significant differences in GDD (Table 3). In general, for R1, R2 and R4 stages, P additions caused earliness, especially in 2014. However, this trend was adverse in R6 growth stage. For the above interaction, early maturity affected SW positively as negative correlations between GDD and SW were found at R2, R4 and R6 stages with that at R4 being the strongest (Figure 3).
Mean comparison of growing degree-days (GDD, °C) for the P rate × growth season interaction for the four reproductive growth stages (R1, R2, R4 and R6) that showed significant differences
. | . | GDD (°C) . | |||
---|---|---|---|---|---|
Interaction . | . | R1 . | R2 . | R4 . | R6 . |
Ρ0 | 2014 | 1025.24a | 1161.98a | 1308.80a | 1463.19b |
2015 | 970.09bc | 1051.86d | 1198.33d | 1362.84c | |
Ρ30 | 2014 | 971.99bc | 1125.70b | 1287.55b | 1478.33a |
2015 | 975.53b | 1037.72e | 1189.78d | 1362.84c | |
Ρ60 | 2014 | 963.90c | 1122.63c | 1276.15c | 1478.33a |
2015 | 966.93bc | 1035.22e | 1180.50e | 1358.68d | |
Ρ90 | 2014 | 970.52bc | 1117.38bc | 1268.11c | 1478.33a |
2015 | 966.93bc | 1035.22e | 1175.34e | 1360.10cd |
. | . | GDD (°C) . | |||
---|---|---|---|---|---|
Interaction . | . | R1 . | R2 . | R4 . | R6 . |
Ρ0 | 2014 | 1025.24a | 1161.98a | 1308.80a | 1463.19b |
2015 | 970.09bc | 1051.86d | 1198.33d | 1362.84c | |
Ρ30 | 2014 | 971.99bc | 1125.70b | 1287.55b | 1478.33a |
2015 | 975.53b | 1037.72e | 1189.78d | 1362.84c | |
Ρ60 | 2014 | 963.90c | 1122.63c | 1276.15c | 1478.33a |
2015 | 966.93bc | 1035.22e | 1180.50e | 1358.68d | |
Ρ90 | 2014 | 970.52bc | 1117.38bc | 1268.11c | 1478.33a |
2015 | 966.93bc | 1035.22e | 1175.34e | 1360.10cd |
Within columns, means labeled with the same letter did not differ significantly at P < 0.05 using the LSD test.
Correlation between GDD and seed weight for the P rate × growth season interaction at R2, R4 and R6 growth stages. Where GDD, growing degree-days; R2, full bloom stage; R4, flat pod stage; R6, full pod filling stage.
Correlation between GDD and seed weight for the P rate × growth season interaction at R2, R4 and R6 growth stages. Where GDD, growing degree-days; R2, full bloom stage; R4, flat pod stage; R6, full pod filling stage.
Full maturity (R8) was influenced only by cultivars and the cultivar × growth season interaction, with cv. Flip recording the lowest GDD (1649.33 °C) and showing no significant differences for the two growth seasons (Table 4). On the contrary, cv. Ikaria showed the highest GDD, summed up to 1790.02 °C, a score very close to the upper limit (1806 °C) of GDD necessary for lentil maturity in Montana (Miller et al. 2001), an environment contrary to the Mediterranean of the present study. Early flowering and maturity genotypes have been proposed to avoid drought but because of their inability to react when there is available soil moisture they show low yields (Materne & Siddique 2009). However, no significant correlation was evident between GDD and SW for the cultivar × growth season interaction at R8 stage. It has been reported that P fertilization leads to earliness in many crops, however it was not always correlated with higher yields (Borges & Mallarino 2000).
Mean comparison of GDD for the cultivar × growth season interaction for the full maturity (R8) stage
Interaction . | . | GDD (°C) . |
---|---|---|
Samos | 2014 | 1748.23b |
2015 | 1731.01c | |
Thessaly | 2014 | 1748.23b |
2015 | 1704.49d | |
Flip | 2014 | 1645.03e |
2015 | 1653.63e | |
Ikaria | 2014 | 1748.23b |
2015 | 1831.82a |
Interaction . | . | GDD (°C) . |
---|---|---|
Samos | 2014 | 1748.23b |
2015 | 1731.01c | |
Thessaly | 2014 | 1748.23b |
2015 | 1704.49d | |
Flip | 2014 | 1645.03e |
2015 | 1653.63e | |
Ikaria | 2014 | 1748.23b |
2015 | 1831.82a |
Within columns, means labeled with the same letter did not differ significantly at P < 0.05 using the LSD test.
In general, P supply stimulates root growth by increasing root biomass and length especially in P-deficient soils (Li et al. 2011). As a result, P addition has been reported to improve drought tolerance by increasing WUE in leguminous crops (Jin et al. 2014). Carbon isotope discrimination (Δ) has been extensively used as an indirect, long-term assessment of WUE and it was found to correlate with yield in C3 species (Turner 1996; Lopes & Reynolds 2010). However, in the present study, neither cultivars nor P rates had a significant effect on Δ (Table 5) indicating that any yield differentiation among the cultivars or P treatments was not ascribed to a variation in water economy of the lentils. Previously, significant variation for Δ was found in lentil genotypes but no correlation with yield was reported (Matus et al. 1996), being in agreement with findings in other grain legumes like dry beans (Tsialtas et al. 2011).
Carbon isotope discrimination (Δ, ‰) measured on the above-ground biomass at R2 growth stage for cultivars and P rates
Cultivars . | Δ (‰) . |
---|---|
Samos | 22.78 |
Thessaly | 22.65 |
Flip | 22.85 |
Ikaria | 22.48 |
P rates | |
Ρ0 | 22.72 |
Ρ30 | 22.75 |
Ρ60 | 22.73 |
Ρ90 | 22.56 |
Cultivars . | Δ (‰) . |
---|---|
Samos | 22.78 |
Thessaly | 22.65 |
Flip | 22.85 |
Ikaria | 22.48 |
P rates | |
Ρ0 | 22.72 |
Ρ30 | 22.75 |
Ρ60 | 22.73 |
Ρ90 | 22.56 |
No significant differences were found at P < 0.05 using the LSD test.
CONCLUSIONS
Phosphorus effects on earliness were cultivar specific, inducing earliness in cvs. Flip and Thessaly and causing delayed maturity in Samos and Ikaria. The critical stage was found to be the beginning of flowering and earliness seems to be a factor contributing to higher yields. Neither cultivars nor P rates had a significant effect on water economy of the lentils. A proper combination of cultivar and P rate can be a means to increase lentil yields under changing Mediterranean conditions.
ACKNOWLEDGEMENTS
The authors are grateful to Dr D. Vlachostergios, HAO-DEMETER, Institute of Industrial and Fodder Plants, Larissa, Greece for providing the lentil seeds for the experimentation. We also thank Dr B. Harlow, Washington State University, School of Biological Sciences, Pullman, WA, USA for the carbon isotope determinations.