Next Article in Journal
Corrections of BDS Code Multipath Error in Geostationary Orbit Satellite and Their Application in Precise Data Processing
Next Article in Special Issue
Phenomic and Physiological Analysis of Salinity Effects on Lettuce
Previous Article in Journal
A Virtual Force Algorithm-Lévy-Embedded Grey Wolf Optimization Algorithm for Wireless Sensor Network Coverage Optimization
Previous Article in Special Issue
Time-Series Chlorophyll Fluorescence Imaging Reveals Dynamic Photosynthetic Fingerprints of sos Mutants to Drought Stress
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploration of Chlorophyll a Fluorescence and Plant Gas Exchange Parameters as Indicators of Drought Tolerance in Perennial Ryegrass

by
Piotr Dąbrowski
1,
Aneta H. Baczewska-Dąbrowska
2,*,
Hazem M. Kalaji
3,
Vasilij Goltsev
4,
Momchil Paunov
4,
Marcin Rapacz
5,
Magdalena Wójcik-Jagła
5,
Bogumiła Pawluśkiewicz
1,
Wojciech Bąba
6 and
Marian Brestic
7
1
Department of Environmental Improvement, Faculty of Civil and Environmental Engineering, Warsaw University of Life Sciences-SGGW, 02-776 Warsaw, Poland
2
Polish Academy of Sciences Botanical Garden-Center for Biological Diversity Conservation in Powsin, 02-973 Warsaw, Poland
3
Department of Plant Physiology, Faculty of Agriculture and Biology, Warsaw University of Life Sciences-SGGW, 02-776 Warsaw, Poland
4
Department of Biophysics and Radiobiology, Sofia University “St. Kl. Ohridski”, 1164 Sofia, Bulgaria
5
University of Agriculture in Kraków, Faculty of Agriculture and Economics, Department of Plant Physiology, 30-239 Krakow, Poland
6
Department of Plant Ecology, Institute of Botany, Faculty of Biology, Jagiellonian University, 30-387 Krakow, Poland
7
Department of Plant Physiology, Slovak University of Agriculture in Nitra, 949 76 Nitra, Slovakia
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(12), 2736; https://doi.org/10.3390/s19122736
Submission received: 20 April 2019 / Revised: 18 May 2019 / Accepted: 4 June 2019 / Published: 18 June 2019
(This article belongs to the Special Issue Chlorophyll Fluorescence Sensing in Plant Phenotyping)

Abstract

:
Perennial ryegrass (Lolium perenne L.) belongs to the common cultivated grass species in Central and Western Europe. Despite being considered to be susceptible to drought, it is frequently used for forming the turf in urban green areas. In such areas, the water deficit in soil is recognized as one of the most important environmental factors, which can limit plant growth. The basic aim of this work was to explore the mechanisms standing behind the changes in the photosynthetic apparatus performance of two perennial ryegrass turf varieties grown under drought stress using comprehensive in vivo chlorophyll fluorescence signal analyses and plant gas exchange measurements. Drought was applied after eight weeks of sowing by controlling the humidity of the roots ground medium at the levels of 30, 50, and 70% of the field water capacity. Measurements were carried out at four times: 0, 120, and 240 h after drought application and after recovery (refilling water to 70%). We found that the difference between the two tested varieties’ response resulted from a particular re-reduction of P700+ (reaction certer of PSI) that was caused by slower electron donation from P680. The difference in the rate of electron flow from Photosystem II (PSII) to PSI was also detected. The application of the combined tools (plants’ photosynthetic efficiency analysis and plant gas exchange measurements) allowed exploring and explaining the specific variety response to drought stress.

1. Introduction

Urban green areas fulfill numerous functions, the more important being: the reduction of pollutants by phytoremediation [1,2,3] and the regulation of the urban microclimate [4]. On that basis, the concept of green infrastructure was introduced. The basic purpose of that idea was to sustain the high complexity and dynamism of urban areas. Green infrastructure can be outlined as a network of natural and semi-natural areas, where combinations of different functions (ecological, social, and economic) coexist as a whole at different spatial scales, from urban centers to peri-urban areas. The turf, besides the trees and shrubs, is one of the most important parts of green infrastructure. When kept in good condition, it can improve the aesthetic value of the whole city.
Developing of urban areas is a main cause of the extensive transformation and alteration of land cover, biogeochemical cycles, air quality, and climate. Moreover, the environmental changes at the local scale may even have a stronger impact on vegetation than the global changes. The urban environment is characterized by changes in the surface energy balance, the high concentration of gases, pollutants, and particulate matter (PM) or rainfall patterns, and the promotion of the urban heat island effect [5,6]. These features of urban areas may enhance the stress pressure on the vegetation’s functionality in the form of abiotic and biotic factors, especially water stress [7].
Drought considerably affects plant productivity, which is connected mainly with the reduction of photosynthetic activity. That phenomena may be driven by both closing stomata and complex non-stomatal effects, which are not fully recognized nor understood [8]. The changes in Photosystem II (PSII) activity observed during drought may be the effects of both limitations and the effects dependent on the degrees and on the species. Under water deficit, inhibition of PSII activity has been proven for some agricultural species, e.g., wheat [9], barley [10], trees [11], and also Lolium/Festuca grasses [12,13,14]. On the other hand, the detailed information about the simultaneous drought effect on Photosystem I (PSI) and PSII activity in perennial ryegrass turf varieties is missing.
Perennial ryegrass (Lolium perenne L.) is the grass species of common use for turf in urban areas in all of Central and Western Europe; it is also used as a forage crop [15]. Since the advent of turf varieties in the 1960s, the development of perennial ryegrass varieties has increased steadily. In 2017, more than 400 turf varieties of that species were registered in the EU [16].
Many authors have stated that chlorophyll a fluorescence (ChlF) parameters are useful for the characterization of the responses of plants to a variety of environmental stress factors and for phenotyping [17,18,19,20]. On the basis of the ChlF technique, it is possible to measure the plant vigor under negative environmental conditions. The measured parameters of that phenomenon are strictly correlated with the PSII functioning [17,18,21]. Green plants under stress conditions change their photosynthetic metabolism. These changes may be direct or indirect, but they are usually correlated with the changes in the yield of ChlF. On the basis of the fluorescence quenching analyses, it is possible to deliver the information about the energy absorption, the utilization, the dissipation, and the electron transport in PSII.
The photosynthetic function of the plants under stress conditions can be estimated by the analysis of the polyphasic fast chlorophyll transient. That quick and nondestructive method (which is also called prompt fluorescence (PF)) was developed by Strasser and Strasser [22]. The photosynthetic apparatus’ photochemical efficiency is reflected by fluorescence kinetics, and it delivers valuable data about the functional and structural features of the mechanisms engaged in photosynthetic electron transport [23,24]. The fluorescence rises during the first second of illumination from the initial (Fo) to the maximal (Fm) value of fluorescence and is shown as a curve with several phases (labeled as O, K, J, I, P) [25,26,27].
The PSII antenna complex can also reemit the light absorbed for photochemical reactions, such as delayed fluorescence (DF). However, there are essential differences between PF and DF. The PF signal can be generated by a particular molecule, whereas DF depends on system relations [28,29,30,31,32]. All the redox reactions of the photosynthetic electron transport between PSII and PSI, as well as the total electron transfer reactions in the reaction centers (RC) of PSII (donor and acceptor side) are revocable. The agglomeration of electrons in the electron transport chain between PSI and PSII brings the electron back, and charge recombination in PSII RC occurs. Subsequently, the RC is re-excitated, and fast energy transfer from it repopulates the excited chlorophyll state of the PSII antenna [33].
The modulated light reflection (MR) signal measured at 820 nm delivers information about electron transport after the plastoquinone (PQ) and to the PSI acceptors [33,34]. Based on that, it is possible to detect the influence of stress factor on the redox state of the PSI RC and plastocyanin (PC).
In the opinion of Goltsev et al. [28], the parallel analyses of the PF, DF, and MR 820 are essential for obtaining the comprehensive relation of the impact of stress factors on electron transport chain (ETC). Because of this correlation, it is possible to reveal the importance of such simultaneous analyses, when adding substance donating or accepting electrons at characteristic locations in the ETC [30].
Regarding the hypothesis that the grasses from various places of origin have different adaptive responses to ecological conditions and drought stress may impair multiple sites of the photosynthetic electron chain, we wanted to know: (i) how the photosynthetic apparatus of two perennial ryegrass varieties responds to drought and which mechanisms of the photosynthetic processes are involved, (ii) if there are differences in response of PSII, PSI, and gas exchange between the tested varieties, and (iii) whether the recovery capacity of photosynthetic performance after re-watering can be estimated by chlorophyll fluorescence measurements. Therefore, the detailed in vivo analysis and the comparisons of the changes in PSI and PSII photochemistry induced by drought stress in two turf varieties of perennial ryegrass by means of the parameters derived from the prompt chlorophyll fluorescence and delayed chlorophyll fluorescence records and the MR 820 signal were the main aims of this study. We believe that our work expands the knowledge about the differences between these two varieties in the photosynthetic response to drought stress and provides clues to elucidate the possible photosynthetic mechanism by which perennial ryegrass tolerates drought stress. Broadening our knowledge about the impact of the altered response of plants to the fluctuating environment will help to improve predictions of the effect of future climate change on ecosystems and plant communities in urban areas.

2. Materials and Methods

2.1. Plants, Growth Conditions, and Experiment Design

Two turf varieties of perennial ryegrass (Lolium perenne L.), “Roadrunner” (Turf Seed, Inc., Gervais, OR, USA) and “Nira” (Małopolska Hodowla Roślin, Kraków, Poland), were the objects of the experiment. These twos varieties are one of the most commonly used for creating lawns in urban areas in Central Europe. The experiment was conducted in the greenhouse of Warsaw University of Life Sciences (SGGW), under natural lighting conditions (day/night, with a 16-hour day length). The average temperature was equal to 22 °C. The seeds (0.47 g/pot) of both varieties were sown on 11 March 2016 in 13 × 13 × 13 cm pots filled with a mixture of sand (70%) and clay (30%). The single dose of multi-component fertilizer Substral 100 (30 g·m−2) was applied. The drought stress was applied 8 weeks after sowing by smaller water dosing (compared to the control). The plants lost water through natural evapotranspiration. Water was dispensed to maintain a homogeneous field water content (FWC) of the substrate. Three treatments of water content were applied: 70% (control), 50%, and 30%. The soil moisture was controlled daily by use of ProCheck needle probe (Decagon Devices, Inc., Pullman, WA, USA). The experiment was set in a split-plot system, where the first treatment was varieties, and the second was level of drought. The whole experiment consisted of 36 pots.
All of the measurements were carried out in uniform conditions on flag leaves in 6 repetitions. The measurements were made at 4 time points. The first time point was performed before drying the soil. The next two time points were conducted 120 and 240 h after the stress application. The last was conducted 160 h after re-watering to 70%.

2.2. The Chlorophyll Fluorescence and the Gas Exchange Measurements

Chlorophyll fluorescence measurements were conducted on 8-week-old plants. Six leaves were dark adapted for 30 min before measurements started. The M-PEA Chlorophyll Fluorescence Measurements System (Multi-Function Plant Efficiency Analyzer, Hansatech, Pentney, UK) was used. The following measurement protocol was applied: measurement time 1.0 s, intensity of actinic light 3000 μmol·m−2·s−1, wavelength 635 ± 10 nm. The following chlorophyll fluorescence signals were measured: prompt chlorophyll a fluorescence, delayed chlorophyll a fluorescence, and 820-nm light reflection.
The JIP-test, described by Strasser et al. [33], was used to recalculate the characteristic points of photoinduced chlorophyll fluorescence transients to specific parameters of the light phase of photosynthesis. According to Goltsev et al. [35], the characteristic points (I1, I2, and D2) of the DF curve were also assessed. The I1 point is the first maximum of the curve, the I2 point the second maximum, and D2 the second minimum of the curve. Two ratios were calculated: (I1D2)/D2 and I1/I2. The MR signal was measured simultaneously with DF. The ratio MRt/MR0, where MRt is the modulated 820-nm reflection intensity at time t and MR0 is the value of the 820-nm reflection of the sample at the onset of the actinic illumination, was calculated. The minimum (MRmin) and maximum (MRmax) levels were also calculated, as well as ΔMRfast (difference between MR0 and MRmin) and ΔMRslow (difference between MRmax and MRmin).
Plants’ gas exchange in terms of net CO2 assimilation (A), transpiration rate (E), stomatal conductance (gs), and sub-stomatal CO2 concentration (Ci), were measured by a portable gas analyzer Lcpro+ (ADC, Bioscentific ltd, Hoddesdon, U.K.). That open-gas exchange system was operated during the measurement in differential mode at a 150-mmol·s−1 flow rate of ambient air. The air temperature during measurements was 22 °C. Saturating irradiance was about 1000 μmol·m−2·s−1. Air relative humidity was 60%. The measurements were made on 6 leaves in each treatment, always after the stabilization of conditions in the chamber. The descriptions of fluorescence parameters, DF, MR820, and gas exchange parameters are described in Table 1.

2.3. Statistical Analysis

Both chlorophyll fluorescence and gas exchange parameters were statistically analyzed by the two-way ANOVA model, where the first factor was the level of field water capacity and the second was the variety. The Fischer’s least significant differences test was used as a post hoc test at a 0.05 confidence level. The Statistica 10.0 program (Statsoft, Inc. Tulsa, OK, USA) was used to perform the statistical analysis. The relationship between PF, DF, MR 820 signal, and CO2 assimilation was assessed on the basis of the Pearson correlation coefficient. The values of the coefficients were considered only when significant at p = 0.05.

3. Results

3.1. Prompt Chlorophyll a Fluorescence

Water shortage significantly influenced the fluorescent transients in both varieties, especially at the I and P steps of the chlorophyll fluorescence induction curve (Figure 1). In the case of the “Roadrunner” variety grown under control conditions for more than 120 h, the fluorescence signal measured at the I step was at the level of 40 relative units (rel. u.), but after 240 h of 50% FWC stress application, it decreased significantly (by 25%). At this step, re-watering the plants to obtain 70% FWC caused an increase of the signal at this step to 35 rel. u. The signal at the I step after 240 h in plants growing on substrate with 30% FWC was at the level of 26 rel. u., but after this time, it was impossible to take measurements due a severe damage of the plants under that treatment. Similarly, we also were not able to make any measurements of plants’ gas exchange under these conditions. The P step measured in control plants was at the level 45 rel. u.; however, after 240 h of 50% FWC treatment, it was significantly lowered by 26.7%, and under 30% FWC, it was significantly reduced by 33.4%, as compared to control plants. Re-watering plants (raising FWC from 50% to 70%) caused a slight increase of the signals at this step. In the case of the “Nira: variety growing before the stress application, the I step was at the level of 40 rel. u., but after 240 h of 50% FWC treatment, it decreased significantly (by 25%). The I step after 240 h in plants growing on substrate with 30% was decreased by 35%. The P step measured in control plants was at the level 45 rel. u., but after 240 h of stress application, this value decreased significantly by 26% and 22.3% under 50% and 30% FWC conditions, respectively. The re-watering caused an increase in this value in both treatments, but not to the level of control plants.
To visualize the influence of drought stress on the dynamics of the chlorophyll fluorescence induction curve, the relative variable fluorescence, Vt = (Ft − FO)/(FM − FO), was calculated. In the next step, the calculations of the differences in relative variable fluorescence curves were made by subtracting the normalized fluorescence values (between O and P steps) recorded in control plants under drought stress (ΔV = Vt(Xh) − Vt(0h)) [33]. The character of the OJIP fluorescence transients recorded in plants under stress of both varieties differed from those recorded in non-stressed plants (Figure 2). The major changes in the curves’ course of stressed plants were detected in the O–J phase and I–P phase. In the O–J phase, the curve drawn for stressed plants had a higher course than for the control plants, but in the I–P phase, these curves had a lower course. There were also significant differences between the tested varieties.
For complete analysis of drought-induced changes in the OJIP curve, the differential curves were presented separately for main bands occurring during O-P transients: K-, L-, H-, and G-bands. The curves for these bands were constructed by subtracting the normalized fluorescence values (between O and K, O and J, J and P, or I and P, respectively) recorded in control plants from those recorded in plants under drought stress (Appendix A). The appearance of the K-band and L-peak in the fluorescence transients of stressed plants further pointed out the differences between tested varieties in terms of their response to the appearance of stress.
The case of L-band was quite different in comparison with the K-band for both varieties. The most visible changes occurred not 240 h, but 120 h after stress application. It should be noted that for that variety treated by 30% FWC, only during one treatment, the L-band was at a higher level than in the control plants.
The H-band was the parameter significantly decreased by the stress factor after 240 h. However, a significant decrease was also observed after 120 h.
The G-band parameter was significantly decreased by the stress application only after 240 h.
The OJIP transients were also transformed to the biophysical parameters [33]: the quantum yields (φPo, φEo, φRo, φDo, ΨEo, δRo), specific activities per reaction center (RC), and performance indices (PI). The values of the calculated parameters were normalized to those of the plants in the control (non-stress) conditions (0 h). The deviations of the action patterns of the stressed and control plants are shown as radar plots (Figure 3). The decrease in field water capacity caused significant differences in some parameters of quantum yields for both varieties. However, the changeable parameters depended on the variety. For the Roadrunner variety, the φDo and δRo parameters were changed under drought stress. Moreover, these changes were greater in plants treated with 30% FWC, than ones treated with 50% FWC. In the case of the Nira variety, only the φDo parameter was significantly increased by both drought treatments, but once again, the changes were greater in the plants treated with 30% FWC than in plants treated with 50% FWC. Regardless of the variety and the degree of drought, the re-watering caused practically full recovery of the parameters to the control values.
The performance index PIABS was used to quantify the PSII performance, and the performance index PItotal was used to measure the performance up to the reduction of PSI end electron acceptors. The performance indexes were also significantly modified by drought stress in both varieties; decreases in these parameters were observed. For the Roadrunner variety under 30% FWC, both parameters decreased significantly in comparison to others. Regardless of the variety and the degree of drought, the re-watering caused partial recovery of these parameters to the same values as in control plants.
Drought stress caused significant changes in the parameter classified as the specific activity of the reduction of end electron acceptors of PSI per reaction center (REo/RC). In both varieties, the values of REo/RC already increased significantly in the plants growing under 30% FWC. Moreover, the changes were also significant for the Roadrunner variety under 50% FWC.

3.2. Delayed Chlorophyll a Fluorescence

In this work, we focused on the first recorded delayed fluorescence (DF) decay component that appeared within a microsecond time range. In the initial part of the DF induction kinetics (during the first 1 s of illumination), that emission was the principal share. The changes of the delayed fluorescence induction curves depended on the field water capacity and the variety. The DF values diminished with increasing the time of exposure to stress for both varieties (Figure 4). For the Roadrunner variety growing under 50% FWC, a noticeable decrease in the intensity of curves after 240 h of stress application was observed. At that time point, the first maximum of the curve (I1) decreased from 154 rel. u. to 94 rel. u. For the re-watering plants, the curve had a similar course as after 240 h. In this variety treated with 30% FWC, the visible changes were also obtained 240 h after the application of stress, decreasing maximally from 142 rel. u. to 72 rel. u. For the Nira variety treated with 50% FWC, the first changes in the course of the curve were noted 120 h after stress application; at this time, the I1 maximum decreased from 161 rel. u. to 140 rel. u. Even greater changes were recorded after 240 h from the stress application. The value of the I1 parameter decreased to 101 rel. u. The re-watering caused an increase in this parameter to 111 rel. u. For the plants treated with 30% FWC, the changes were observed only 240 h after stress application, when the maximum of the curve decreased from 166 rel. u. to 97 rel. The re-watering did not cause any visible changes in the values of the maximum. It should be noted that there were no significant differences between the minimums of the curve (D2) measured in the control and drought-stressed plants.
To better demonstrate the changed shape of the DF, the ratios between the maxima and minima of the induction curves were determined (Figure 5). The ratio I1/I2 in plants treated with 50% and 30% FWC changed only for the Roadrunner variety after 120 and 240 h. The re-watering from 50%–70% FWC restored the I1/I2 to the value for control plants. The course of the changes of the ratio (I1 − D2)/D2 was similar to the changes of the ratio I1/I2. However, for this ratio, the changes were more visible and occurred also for the Nira variety.

3.3. The Modulated Light Reflection Signal Measured at 820 nm

The reflection at 820 nm was significantly affected by drought stress in both varieties (Figure 6). In both varieties, the values of the MRmin parameter were the lowest 240 h after the application of stress in plants treated with both field water capacities. However, the most significant decrease was noted for the Roadrunner variety under 30% FWC. The re-watering increased the course of the curve to a similar level as in control plants in all treatments. The MRmax parameter decreased significantly 240 h after applying drought stress in plants treated only with 30% FWC. For the Roadrunner variety, the decrease in that parameter was more pronounced. For the Nira variety, the re-watering caused alignment to the values similar to that of the control plants.
The rate of initial photooxidation of reaction center of PS I (P700) (Vox) for the Roadrunner variety treated with drought conditions increased in the comparison with control plants 240 h after the application of stress (Figure 7). However, the re-watering did not cause any changes: there were still significant differences between the treatments. In the case of the Nira variety, the most visible changes were recorded 120 h after the stress application, but they were less noticeable than for the Roadrunner variety. Moreover, the re-watering caused alignments to the values similar to those of the control plants: the rate of re-reduction of P700+ by electrons from PS II and the reduced PQ pool (Vred.). For the Roadrunner variety treated with 30% FWC, this was decreased in comparison to control plants 240 h after stress application. In the case of the Nira variety, the drought treatment did not cause any changes of this parameter.

3.4. Gas Exchange Analyses

Plants’ CO2 assimilation (A) was significantly influenced by drought in the case of both varieties (Figure 8). In the case of the Roadrunner variety, the significant differences between the values measured for all three treatments were detected 240 h after stress application. The CO2 assimilation measured in the control plants was equal to 3.2 µmol CO2 m−2·s−1 and decreased by 65% for the plants under 50% FWC and by 95% for the plants under 30% FWC. The re-watering from 50%–70% FWC caused the value of that parameter to increase to similar values as for the control plants. In the case of the Nira variety, significant differences were also observed 240 h after stress application. However, there were no significant differences between plants treated with 50 and 30% FWC. The CO2 assimilation measured in the control plants was equal to 4.2 µmol CO2 m−2·s−1 and decreased by 73% and 80% respectively in plants under drought stress. Re-watering from 50% and 30%–70% FWC increased the value of the parameter to a similar level as in the control plants.
Plants’ transpiration rate € in the case of the Roadrunner variety was affected significantly 240 h after the stress application. In the control plants, that parameter was equal to 1.5 mmol H2O m−2·s−1 and decreased to 0.4 mmol H2O m−2·s−1 for the plants treated with both levels of FWC. The re-watering from 50%–70% caused recovery of that parameter to similar values as for the control plants. In the case of the Nira variety, the transpiration was affected significantly also 240 h after the application of stress. The value measured on the control plants was at a level of 2.4 mmol H2O m−2·s−1. Drought stress caused a decrease in the E-value to 0.4 mmol H2O m−2·s−1 (83% in the comparison to the control) in the plants treated with both levels of FWC. The re-watering from 50% and 30%–70% FWC increased the values of that parameter, but they were still lower in comparison to control plants.
Stomatal conductance (gs) was also affected significantly by the drought stress. For the Roadrunner variety, a significant decrease was noted 240 h after stress application. The value measured for the control was equal to 0.30 mol H2O m−2·s−1. The value measured in the plants treated with the 50% FWC decreased by 67% and by 83% for the 30% FWC treatment when compared to the control. The values measured in treatments with both levels of FWC comprised one homogenous group. The re-watering caused an increase in the values of this parameter, but they were still lower than for the control plants. In the case of the Nira variety, the significant differences were noted 120 h after stress application. The value measured for control plants was equal to 0.36 m−2·s−1 and decreased by 34% in plants treated with 30% FWC. At the next time points (240 h), the value measured in the control plants was equal to 0.35 mol H2O m−2·s−1, decreasing by 57% for the plants treated with 50% FWC, and by 74% for the plants treated with 30% FWC. The significance of the differences between the 50% and 30% treatments was found. The re-watering from 50% and 30%–70% FWC caused an increase in the values of that parameter, but they were still lower in comparison to the control plants.
The internal CO2 concentration (Ci) for the Roadrunner variety was significantly influenced 240 h after the stress application. For the control plants, that parameter was at the level of 270 (µmol CO2 mol−1 air) and was higher only for the plants treated with 30% FWC (331 µmol CO2 mol−1 air). In the case of the Nira variety, the value of this parameter measured in control plants ranged from 250–300 (μmol CO2 mol−1 air), and there were no significant changes between them and the plants treated by drought.

3.5. Relationships among ChlF Parameters and Gas Exchange

On the basis of the statistical analysis, it can be noted that there were significant correlations between some PF parameters and the CO2 assimilation parameter (A) (Figure 9). First of all, the significant correlations were noted in both varieties only in stress conditions. Under drought stress in the Roadrunner variety, significant correlations were noted in the majority of the JIP test parameters. Only REo/RC and TRo/RC were not correlated with the A parameter. In the Nira variety also, only two parameters were not correlated with A: φRo and ETo/RC.
On the basis of the statistical analysis, it can be noted that there were significant correlations between some DF parameters and the CO2 assimilation parameter (Figure 10). For the Roadrunner variety in control conditions, there were no significant correlations between the A and all DF parameters, but in the plants treated with 30% FWC, A was significantly correlated with almost all DF parameters, and the correlation coefficient ranged from 0.84–0.94. Only the parameters I2 and D2 were not correlated with A. A similar phenomenon was noticed in the case of the Nira variety: A in the plants treated with drought conditions was significant correlated with all parameters (r ranged from 0.69–0.84), except I2 and D2. The statistical analysis showed also the significant correlations between the A parameter and the MR curve parameters in both varieties treated with drought conditions. There were no significant correlations in control plants. For the Roadrunner variety, the parameters MRmin, MRmax, Vred, and ΔMRfast were significant and positively correlated with the A parameter (r ranged from 0.84–0.90), while the parameters ΔMRslow and Vox were correlated negatively with the A parameter (r ranged from −0.90 to −0.78). In the Nira variety, the parameters MRmin and MRmax were correlated positively (r ranged from 0.59–0.80), while the parameters ΔMRslow and Vox were correlated negatively (r ranged from −0.59 to −0.57).

4. Discussion

Drought is one of the main factors that inhibits plant functioning [36]. Excessive drought brings increased osmotic stress, promptly [37]. Despite numerous studies conducted to understand the mechanism by which drought stress affects plant photosynthesis, this topic is still under investigation [38,39]. Generally, factors that limit photosynthesis under drought stress can be divided into stomatal or nonstomatal ones. However, in most study cases, both effects seem to act simultaneously [40]. Reduction of the A parameter is attributed to stomatal limitation if Ci decreases and gs increases with the decrease in soil FWC. If not, the reduction of A is attributed to nonstomatal limitation [41]. However, such reactions are dependent on the studied genotype, the growth stage of the plant, and the level of drought stress.
Ohashi et al. [42] showed that, under drought stress, stomatal closure contributes to maintaining the high leaf water content; however, it leads to a decrease in leaf photosynthesis. Our results proved that under drought stress, the photosynthesis (A) and the transpiration (E) decreased. Similar phenomena were noticed for stomatal conductance (gs). The photosynthesis intensity of stressed plants recovered to a value similar to that of plants grown under control conditions after re-watering the soil. This proves that both perennial ryegrass varieties were able to repair metabolic processes within a short time, with the particular case of electron transfer between the photosystems. Chlorophyll a fluorescence parameters measured after re-watering reached values similar to those measured for control plants (nearly complete reactivation of photosynthetic processes occurred). These results are in agreement with those obtained by Pukacki and Kamińska-Rożek [43].
The mechanism limiting photosynthesis under drought stress is complex, especially when drought reaches a serious level of severity. Our study showed that the analysis of the leaf gas exchange parameters alone was insufficient to expose the extent of injury in the photosynthetic apparatus caused by drought completely. Additional analyses related to nonstomatal limitation of photosynthesis (based on chlorophyll fluorescence measurements), which can contribute to the reduced activity in the PSII and PSI, helped us to obtain better understanding on the mechanisms of plant tolerance to drought [39].
Measurements of chlorophyll a fluorescence kinetics deliver exhaustive information on the structure and function of the photosynthetic apparatus, in particular PSII [25,44,45]. Thus, drought can directly or indirectly change the kinetics of chlorophyll a fluorescence [46,47]. Changes in chlorophyll a fluorescence parameters connected with drought response were compared before in drought-tolerant and -susceptible genotypes of Lolium perenne, Festuca arundinacea, and their hybrids [12,48]. Accompanying proteomic and lipidomic studies revealed differences in thylakoid connected proteins (e.g., ATP-ase, oxygen-evolving enhancer protein), in enzymes linked to light-independent photosynthetic processes (e.g., RuBisCO activase, chloroplast fructose-bisphosphate aldolase), as well as in proteins involved in plastid gene expression [11,12]. Furthermore, considerable changes in both the lipid composition of chloroplast membranes and the abundance of enzymes involved in their metabolism were observed [12,13].
The JIP-test analyses in the cases of Augea capensis Thunb. and Zygophyllum prismatocarpum E. Meyer ex Sond., from the Namibian Desert, showed that drought caused an increase in the fluorescence of both species under all experimental conditions in comparison to the control plants. The regulation of photochemical activity is carried out primarily through the deactivation of PSII reaction centers. In the opinion of those authors, the PIABS parameter calculated from the JIP-test was a precise and sensitive pointer of the physiological condition of plants, as well as in the field and under laboratory conditions. There was a high correlation between the decreased CO2 assimilation capacity and the decreased PIABS values, which is a multiparametric expression of three independent parameters: the density of reaction centers, which corresponds to the absorption flux, the quantum yield of trapping, and the probability that a trapped exciton will move an e into the ETC beyond QA [49]. That finding provided the evidence of the connection of OJIP fluorescence changes during drought in field or laboratory conditions with changes in the general ability of photosynthesis. The integrated regulation of the entire photosynthesis process takes place in such a way that the internal balance between the efficiency of the photosynthetic light phase reaction and the efficiency of reactions leading to the CO2 assimilation is maintained [36]. Our study confirmed that PIABS is one of the most sensitive parameters amongst JIP-tests for drought stress.
The decrease in the field water capacity caused significant differences in some parameters for quantum yields in both varieties. However, the parameter that would be changed depends on the variety. For the Roadrunner variety, the φDo and δRo parameters changed heavily under drought stress. In the case of the Nira variety, only the φDo parameter was significantly reduced by both drought treatments. The quantum yield parameters are known to be accurate as the indicators of abiotic stress factors. That is because they reflect the lower efficiency of the regulation of the excitation energy utilization and dissipation by the photosynthetic membrane. It is also known that minimizing the reduction of quantum yield under stress should be advantageous for the yield [10]. In the opinion of Oukarroum et al. [50], a lower δRo level indicated a decreased outflow of electrons at the acceptor side of PSI caused by an inactivation of ferredoxin NADP+-reductase. Interestingly, the drought tolerance-dependent differences in ferredoxin NADP+ reductase were observed before in Festuca arundinacea, a drought-tolerant grass genetically similar to and easy to produce hybrids with Lolium perenne [12].
Drought application caused the induction curves for both varieties to be at significantly lower levels in comparison to the control. Similar results were observed by Oukarroum et al. [25] for barley plants subjected to drought stress, albeit mainly for the I–P phase (reduction of electron acceptors), and the variation in fluorescence kinetics resulting from various degrees of tolerance in individual varieties. Drought stress may also change characteristic points of the OJIP curve and reduce fluorescence intensity at the J, I, and P steps [51]. In our study, the L-peak was influenced by increasing FWC, which might be evidence that the stability and structure of PSII RCs, as well as antennae, were affected. At the same time, the occurrence of the K-peak was noted in the OJIP transient, and ΔWK rose with extended treatment duration. This phenomenon can be associated with oxygen evolving complex (OEC) inactivation and/or inhibition of electron transport on the donor or acceptor side of PSII [27]. It was also noted that the K-peak increased significantly under drought stress. Because of that, the L- and K-peaks may be advised as potential indicators for physiological instability under drought stress before the appearance of their visible signs. The oxygen-evolving enhancer protein 2 was reported to be differentially accumulated in drought-tolerant and -susceptible genotypes of Festuca arundinacea [12]. In our work, differences in the activity of the oxygen evolving complex between the tested cultivars can also be proven by the differences in the H- and G-bands. A higher activity can be expected in the case of the Nira variety. Therefore, this is consistent with the results of the study cited above. In this work, we proved that the DF signal may provide important information about the processes in PSII and the photosynthetic membrane in drought-stressed plants. Nevertheless, such an elaborate model as the JIP-test has not yet been proposed for the DF induction curves.
Based on our work, it can be concluded that the shape of the DF induction curve was influenced by this stress. This phenomenon was confirmed by Goltsev et al. [35]. From the analysis of DF, further insights can be derived by Strasser et al. [33] and Goltsev et al. [35]. Goltsev et al. [35] proved that during the initial stages of the desiccation of bean plants, the overall intensities of DF were not affected significantly until the RWC reached about 20%. When the water content was decreased below 20%, the intensity of DF decreased. The authors concluded that in a clear contrast to the overall intensity, the shape of the DF induction curve was sensitive to water stress in the initial stages of the drought stress. In the range of RWC between 100 and 50%, the characteristic peak in the DF induction curve (I2) disappeared. The authors of that work concluded that from the DF parameters, the I1/I2 ratio was most sensitive to mild drought stress. A well expressed linear dependence on RWC between 100% and 40% RWC was observed. Interestingly, again, an almost linear dependence, but not with the same slope, was observed at lower RWC: I1/I2 sharply decreased between 20% and 5% and then, at even lower RWC, again slightly increased.
In particular, I1, I2, and D2 parameters were negatively affected by drought in both tested varieties, indicating an enhanced electron transfer capacity at both the donor side and the acceptor side of PSII; also, re-oxidation of PQH2, as well as oxidation of P700, became slower [52]. The reduction of I2 and D2 probably means that the re-oxidation capacity of PQH2, the formation of the proton gradient, and the reducing activity of the PSII complex were reduced. We also demonstrated that in the case of the Roadrunner variety, that water loss from 70%–30% of FWC caused a significant increase in the ratio between the two maximums in the DF induction curves. This is evidence of the high sensitivity of this technique to show the photosynthetic responses that are not detectable by other methods.
The differences in the drought response between the Roadrunner and the Nira varieties were also demonstrated on the basis of the analysis of the MR parameters. For the Roadrunner variety, a decrease in the Vred parameters was observed. A similar trend, but not so pronounced, was observed for the Nira variety. According to Oukarroum et al. [53] and Gao et al. [52], the decrease in that parameter under drought stress was evidence of a slower re-reduction of P700+, which is the result of slower electron donation from P680. These authors showed also a simultaneous reduced rate of photo-induced electron transfer through PSI (decrease in Vox), which was the result of the reduction of the PSI photochemical activity. Our results did not confirm the decrease of the Vox parameter under drought stress. In the opinion of many authors [34,54], the MR signal is mostly influenced by the rate of electron flow from PSII to PSI. In addition, there is also a possibility that a limitation on the acceptor side of PS I was the reason for the changes in the MR signal. On the basis of our results, it can be concluded that the ΔMRslow and ΔMRfast parameters measured for the drought-stressed plants were significantly changed in comparison to the control. In addition, our results did not confirm the correlation between those parameters and the gas exchange parameters. These results are in opposition to the findings of some other studies. Schansker et al. [34] suggested that these parameters could be affected by stress. These authors concluded that the fast phase of the signal corresponds to the kinetics of the photoinduced changes in the P700 redox state, being significantly changed only by strong stress. At the same time, the slow phase of the MR signal, reflecting a P700+ re-reduction, decreased gradually with stress intensity. In the opinion of Goltsev et al. [55] and Salvatori et al. [56], the deactivation of this kinetic phase was demonstrated as the reduction of the electron transport rate through plastoquinone to P700+, which is the reason for the more stable fast phase of the MR signal.
Altogether, the simultaneous analyses of the JIP-test, DF, MR, and gas exchange showed that drought stress led to deactivation of PSII RCs, decreasing in the connectivity between independent PSII units, and the limitation of electron transport beyond QA. Moreover, the damage of PSI acceptor-side electron transporters and the split between the activities of PSI and PSII were noted. The results obtained from the four types of simultaneously-measured parameters provided much more detailed information about the drought-induced changes in the photosynthetic electron chain than in the case of the application of any of them, separately. This conclusion was confirmed by Zhou et al. [57] for maize. Moreover, in comparison with the less resistant variety, Roadrunner, the destructive influence of drought stress on the electron transport chain of the more resistant variety Nira showed a lesser extent and appeared earlier. In the case of the resistant Nira variety, the drought stress showed a less obvious negative effect on all sections of the OJIP curve, as well as on the energetic connectivity between PSII and PSI oxidation and reduction. At all those points, changes were pronounced faster in the case of the Roadrunner variety. These results proved that the resistant variety may preserve higher photosynthetic activity under drought stress. Thus, it can be said that the simultaneous measurement of OJIP, DF, and MR may be used to distinguish the tolerant perennial ryegrass varieties from the sensitive ones.

5. Conclusions

Analysis of obtained data showed that the reaction of the photosynthetic machinery of perennial ryegrass to drought stress is a very complicated process. On the basis of the simultaneous measurements of PF, DF, MR 820 signals, and plant gas exchange parameters, we were able to observe specific effects of water shortage on the photosynthesis processes. That include the reduction of electron flow rate through PSII and inhibition of electron transfer from the reduced PQ pool to the PSI reaction center (P700+) and reoxidation of QA. Moreover, an inhibition of the quantum yields of photoinduced electron transport in PSII reaction centers to QA together with a reduction of the fast phase of photoinduced kinetics of the MR signal were also detectedunder this stress. Additionally, we found that the combination of such applied in this research methods is potentially useful for the non-destructive estimation of relative water content in leaves. However, due significant differences found between the two tested varieties, it is recommended to extend research to a larger population of perennial ryegrass accessions.

Author Contributions

P.D., H.M.K. and B.P., conceptualization. P.D. and A.H.B.-D., investigation and data curation. M.R., H.M.K., M.W.-J. and M.B. interpreted the basic results and drew the more important conclusions. V.G., M.P. and W.B. were responsible for the data processing, bioinformatics, and statistical analysis. All the authors contributed to the writing of the final version and approved the manuscript. Contribution percentage of authors: Piotr Dąbrowski (5%), Aneta H. Baczewska-Dąbrowska (50%), Hazem M. Kalaji (5%), Vasilij Goltsev (5%), Momchil Paunov (5%), Marcin Rapacz (5%), Magdalena Wójcik-Jagła (5%), Bogumiła Pawluśkiewicz (5%), Wojciech Bąba (5%), and Marian Brestic (10%).

Funding

This research was funded by the Slovak Research and Development Agency under the project APVV-15-0721.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Differential curves of ΔWK (obtained by subtracting the control curve from the first sample in the time from 0.02 to 3 ms) of perennial ryegrass varieties (Roadrunner and Nira) under different field water capacity of soil (50 and 30% FWC) and time (0, 120, 240 and 400 h after stress application): (A) Roadrunner under 50% FWC, (B) Roadrunner under 30% FWC, (C) Nira under 50% FWC, (D) Nira under 30% FWC. n = 6. Relative units.
Figure A1. Differential curves of ΔWK (obtained by subtracting the control curve from the first sample in the time from 0.02 to 3 ms) of perennial ryegrass varieties (Roadrunner and Nira) under different field water capacity of soil (50 and 30% FWC) and time (0, 120, 240 and 400 h after stress application): (A) Roadrunner under 50% FWC, (B) Roadrunner under 30% FWC, (C) Nira under 50% FWC, (D) Nira under 30% FWC. n = 6. Relative units.
Sensors 19 02736 g0a1
Figure A2. Differential curves of ΔWL (obtained by subtracting the control curve from the first sample in time from 0.02 to 0.3 ms) of perennial ryegrass varieties (Roadrunner and Nira) under different field water capacity of soil (50 and 30% FWC) and time (0, 120, 240 and 400 h after stress application): (A) Roadrunner under 50% FWC, (B) Roadrunner under 30% FWC, (C) Nira under 50% FWC, (D) Nira under 30% FWC. n = 6. Relative units.
Figure A2. Differential curves of ΔWL (obtained by subtracting the control curve from the first sample in time from 0.02 to 0.3 ms) of perennial ryegrass varieties (Roadrunner and Nira) under different field water capacity of soil (50 and 30% FWC) and time (0, 120, 240 and 400 h after stress application): (A) Roadrunner under 50% FWC, (B) Roadrunner under 30% FWC, (C) Nira under 50% FWC, (D) Nira under 30% FWC. n = 6. Relative units.
Sensors 19 02736 g0a2
Figure A3. Differential curves of ΔWH (obtained by subtracting the control curve from the first sample time from 1.9 to 30 ms) of perennial ryegrass varieties (Roadrunner and Nira) under different field water capacity of soil (50 and 30% FWC) and time (0, 120, 240 and 400 h after stress application): (A) Roadrunner under 50% FWC, (B) Roadrunner under 30% FWC, (C) Nira under 50% FWC, (D) Nira under 30% FWC. n = 6. Relative units.
Figure A3. Differential curves of ΔWH (obtained by subtracting the control curve from the first sample time from 1.9 to 30 ms) of perennial ryegrass varieties (Roadrunner and Nira) under different field water capacity of soil (50 and 30% FWC) and time (0, 120, 240 and 400 h after stress application): (A) Roadrunner under 50% FWC, (B) Roadrunner under 30% FWC, (C) Nira under 50% FWC, (D) Nira under 30% FWC. n = 6. Relative units.
Sensors 19 02736 g0a3
Figure A4. Differential curves of ΔWG (obtained by subtracting the control curve from the first sample in time from 30 to 1000 ms) of perennial ryegrass varieties (Roadrunner and Nira) under different field water capacity of soil (50 and 30% FWC) and time (0, 120, 240 and 400 h after stress application): (A) Roadrunner under 50% FWC, (B) Roadrunner under 30% FWC, (C) Nira under 50% FWC, (D) Nira under 30% FWC. n = 6. Relative units.
Figure A4. Differential curves of ΔWG (obtained by subtracting the control curve from the first sample in time from 30 to 1000 ms) of perennial ryegrass varieties (Roadrunner and Nira) under different field water capacity of soil (50 and 30% FWC) and time (0, 120, 240 and 400 h after stress application): (A) Roadrunner under 50% FWC, (B) Roadrunner under 30% FWC, (C) Nira under 50% FWC, (D) Nira under 30% FWC. n = 6. Relative units.
Sensors 19 02736 g0a4

References

  1. Radziemska, M.; Gusiatin, Z.M.; Bilgin, A.A. Potential of using immobilizing agents in aided phytostabilization on simulated contamination of soil with lead. Ecol. Eng. 2017, 102, 490–500. [Google Scholar] [CrossRef]
  2. Radziemska, M.; Bilgin, A.; Vaverková, M.D. Application of mineral–based amendments for enhancing phytostabilization in Lolium perenne L. cultivation. Clean-Soil Air Water 2018, 46. [Google Scholar] [CrossRef]
  3. Bęś, A.; Warmiński, K.; Adomas, B. Long-term responses of Scots pine (Pinus sylvestris L.) and European beech (Fagus sylvatica L.) to the contamination of light soils with diesel oil. Environ. Sci. Pollut. Res. 2019, in press. [Google Scholar] [CrossRef]
  4. Sung, C.Y. Mitigating surface urban heat island by a tree protection policy: A case study of The Woodland, Texas, USA. Urban For. Urban Green. 2013, 12, 474–480. [Google Scholar] [CrossRef]
  5. Hidalgo, J.; Masson, V.; Baklanov, A.; Pigeon, G.; Gimeno, L. Advances in urbanclimate modeling. Ann. N. Y. Acad. Sci. 2008, 1146, 354–374. [Google Scholar] [CrossRef]
  6. Fusaro, L.; Salvatori, E.; Mereu, S.; Marando, F.; Scassellati, E.; Abbate, G.; Manes, F. Urban and peri-urban forests in the metropolitan area of Rome: Ecophysiological response of Quercus ilex L. in two greeninfrastructures in an ecosystem services perspective. Urban For. Urban Green. 2015, 14, 1147–1156. [Google Scholar] [CrossRef]
  7. Nosalewicz, A.; Siecińska, J.; Kondracka, K.; Nosalewicz, M. The functioning of Festuca arundinacea and Lolium perenne under drought is improved to a different extend by the previous exposure to water deficit. Environ. Exp. Bot. 2018, 156, 271–278. [Google Scholar] [CrossRef]
  8. Flexas, J.; Medrano, H. Drought-inhibition of photosynthesis in C3 plants: Stomatal and non-stomatal limitations revisited. Ann. Bot. 2002. [CrossRef]
  9. Zivcak, M.; Kalaji, H.M.; Shao, H.-B.; Olsovska, K.; Brestic, M. Photosynthetic proton and electron transport in wheat leaves under prolonged moderate drought stress. J. Photoch. Photobiol. B 2014, 137, 107–115. [Google Scholar] [CrossRef]
  10. Kalaji, H.M.; Rastogi, A.; Živčák, M.; Brestic, M.; Daszkowska-Golec, A.; Sitko, K.; Alsharafa, K.Y.; Lotfi, R.; Stypiński, P.; Samborska, I.A.; et al. Prompt chlorophyll fluorescence as a tool for crop phenotyping: An example of barley landraces exposed to various abiotic stress factors. Photosynthetica 2018, 56, 953–961. [Google Scholar] [CrossRef]
  11. Kalaji, H.M.; Račková, L.; Paganová, V.; Swoczyna, T.; Rusinkowski, S.; Sitko, K. Can chlorophyll-a fluorescence parameters be used as bio-indicators to distinguish between drought and salinity stress in Tilia cordata Mill? Environ. Exp. Bot. 2018, 152, 149–157. [Google Scholar] [CrossRef]
  12. Kosmala, A.; Perlikowski, D.; Pawłowicz, I.; Rapacz, M. Changes in the chloroplast proteome following water deficit and subsequent watering in a high and a low drought tolerant genotype of Festuca arundinacea. J. Exp. Bot. 2012, 63, 6161–6172. [Google Scholar] [CrossRef]
  13. Perlikowski, D.; Kiersznowska, S.; Sawikowska, A.; Krajewski, P.; Rapacz, M.; Eckhardt, Ä.; Kosmala, A. Remodeling of leaf cellular glycerolipid composition under drought and re-hydration conditions in grasses from the Lolium-Festuca complex. Front. Plant Sci. 2016, 7, 1027. [Google Scholar] [CrossRef]
  14. Digrado, A.; Bachy, A.; Mozaffar, A.; Schoon, N.; Bussotti, F.; Amelynck, C.; Dalcq, A.-C.; Fauconnier, M.-L.; Aubinet, M.; Heinesch, B.; et al. Long-term measurements of chlorophyll a fluorescence using the JIP-test show that combined abiotic stresses influence the photosynthetic performance of the perennial ryegrass (Lolium perenne) in a managed temperate grassland. Physiol. Plant 2017, 161, 355–371. [Google Scholar] [CrossRef]
  15. Dąbrowski, P.; Pawluśkiewicz, B.; Kalaji, H.M.; Baczewska, A.H. The effect of light availability on leaf area index, biomass production and plant species composition of park grasslands in Warsaw. Plant Soil Environ. 2013, 59, 543–548. [Google Scholar] [CrossRef]
  16. Common Catalogue of Varieties of Agricultural Plant Species. 36th complete edition. (2017/C 433/01). Available online: https://eur-lex.europa.eu/legal-content/EN/ALL/?uri=CELEX%3AC2017%2F433%2F01 (accessed on 6 June 2019).
  17. Maxwell, K.; Johnson, G.N. Chlorophyll fluorescence: A practical guide. J. Exp. Bot. 2000, 51, 659–668. [Google Scholar] [CrossRef]
  18. Baker, N.R.; Rosenqvist, E. Applications of chlorophyll fluorescence can improve crop production strategies: An examination of future possibilities. J. Exp. Bot. 2004, 55, 1607–1621. [Google Scholar] [CrossRef]
  19. Ouni, Y.; Mateos-Naranjo, E.; Abdelly, C.; Lakhdar, A. Interactive effect of salinity and zinc stress on growth and photosynthetic responses of the perennial grass, Polypogon monspeliensis. Ecol. Eng. 2016, 95, 171–179. [Google Scholar] [CrossRef]
  20. Pollastrini, M.; Holland, V.; Brüggemann, W.; Bruelheide, H.; Dânilâ, I.; Jaroszewicz, B.; Valladares, F.; Bussotti, F. Taxonomic and ecological relevance of the chlorophyll a fluorescence signature of tree species in mixed European forests. New Phytol. 2016, 212, 51–65. [Google Scholar] [CrossRef]
  21. Dąbrowski, P.; Pawluśkiewicz, B.; Baczewska, A.H.; Kalaji, H.M.; Oglęcki, P. Chlorophyll a fluorescence of perennial ryegrass (Lolium perenne L.) varieties under long term exposure to shade. Zemdirbyste-Agriculture 2015, 102, 305–312. [Google Scholar] [CrossRef]
  22. Strasser, B.; Strasser, R. Measuring fast fluorescence transients to address environmental questions: The JIP-test. In Photosynthesis: From Light to Biosphere; Garab, G., Ed.; Kluwer: Dordrecht, Netherlands, 1995; Volume 5, pp. 977–980. [Google Scholar]
  23. Kalaji, H.M.; Bosa, K.; Kościelniak, J.; Żuk-Gołaszewska, K. Effects of salt stress on photosystem II efficiency and CO2 assimilation of two Syrian barley landraces. Environ. Exp. Bot. 2011, 73, 64–72. [Google Scholar] [CrossRef]
  24. Stirbet, A. Govindjee On the relation between the Kautsky effect (chlorophyll a fluorescence induction) and photosystem II: Basics and applications of the OJIP fluorescence transient. J. Photochem. Photobiol. B Biol. 2011, 104, 236–257. [Google Scholar] [CrossRef] [PubMed]
  25. Oukarroum, A.; Madidi, S.E.; Schansker, G.; Strasser, R.J. Probing the responses of barley cultivars (Hordeum vulgare L.) by chlorophyll a fluorescence OLKJIP under drought stress and re-watering. Eviron. Exp. Bot. 2007, 60, 438–446. [Google Scholar] [CrossRef]
  26. Dąbrowski, P.; Baczewska, A.H.; Pawluśkiewicz, B.; Paunov, M.; Alexantrov, V.; Goltsev, V.; Kalaji, H.M. Prompt chlorophyll a fluorescence as a rapid tool for diagnostic changes in PSII structure inhibited by salt stress in Perennial ryegrass. J. Photochem. Photobiol. B Biol. 2016, 157, 22–31. [Google Scholar] [CrossRef] [PubMed]
  27. Meng, L.-L.; Song, J.-F.; Wen, J.; Zhang, J.; Wei, J.-H. Effects of drought stress on fluorescence characteristics of photosystem II in leaves of Plectranthus scutellarioides. Photosynthetica 2016, 54, 414–421. [Google Scholar] [CrossRef]
  28. Goltsev, V.; Zaharieva, I.; Lambrev, P.; Yordanov, I.; Strasser, R. Simultaneous analysis of prompt and delayed chlorophyll a fluorescence in leaves during the induction period of dark to light adaptation. J. Theor. Biol. 2003, 225, 171–183. [Google Scholar] [CrossRef]
  29. Zaharieva, I.; Goltsev, V. Advances on Photosystem II investigation by measurement of delayed chlorophyll fluorescence by a phosphoroscopic method. Photochem. Photobiol. 2003, 77, 292–298. [Google Scholar] [CrossRef]
  30. Guo, Y.; Tan, J. A kinetic model structure for delayed fluorescence from plants. BioSystems 2009, 95, 98–103. [Google Scholar] [CrossRef]
  31. Oukarroum, A.; Goltsev, V.; Strasser, R.J. Temperature Effects on Pea Plants Probed by Simultaneous Measurements of the Kinetics of Prompt Fluorescence, Delayed Fluorescence and Modulated 820 nm Reflection. PLoS ONE 2013, 8, e59433. [Google Scholar] [CrossRef]
  32. Dąbrowski, P.; Kalaji, H.M.; Baczewska, A.H.; Pawluśkiewicz, B.; Mastalerczuk, G.; Borawska-Jarmułowicz, B.; Paunov, M.; Goltsev, V. Delayed chlorophyll a fluorescence, MR 820, and gas exchange changes in perennial ryegrass under salt stress. J. Lumin. 2017, 183, 322–333. [Google Scholar] [CrossRef]
  33. Strasser, R.J.; Tsimilli-Michael, M.; Qiang, S.; Goltsev, V. Simultaneous in vivo recording of prompt and delayed fluorescence and 820-nm reflection changes during drying and after rehydration of the resurrection plant Haberlea rhodopensis. Biochim. Biophys. Acta Bioenerget. 2010, 1797, 1313–1326. [Google Scholar] [CrossRef]
  34. Schansker, G.; Srivastava, A.; Govindjee; Strasser, R.J. Characterization of the 820-nm transmission signal paralleling the chlorophyll a fluorescence rise (OJIP) in pea leaves. Funct. Plant Biol. 2003, 30, 785–796. [Google Scholar] [CrossRef] [Green Version]
  35. Goltsev, V.; Zaharieva, I.; Chernev, P.; Strasser, R. Delayed fluorescence in photosynthesis. Photosynth. Res. 2009, 101, 217–232. [Google Scholar] [CrossRef] [PubMed]
  36. Van Heerden, P.D.R.; Swanepoel, J.W.; Kruger, G.H.J. Modulation of photosynthesis by drought in two desert scrub species exhibiting C3-mode CO2 assimilation. Environ. Exp. Bot. 2009, 61, 124–136. [Google Scholar] [CrossRef]
  37. Kocheva, K.V.; Busheva, M.C.; Georgiev, G.I.; Lambrev, P.H.; Goltsev, V.N. Influence of short-term osmotic stress on the photosynthetic activity of barley seedlings. Biol. Plant. 2005, 49, 145–148. [Google Scholar] [CrossRef]
  38. Lawson, T.; Oxborough, K.; Morison, J.L.; Baker, N.R. The responses of guard and mesophyll cell photosynthesis to CO2, O2, light, and water stress in a range of species are similar. J. Exp. Bot. 2003, 54, 1743–1752. [Google Scholar] [CrossRef] [PubMed]
  39. Liao, J.X.; Wang, G.X. Effects of Drought Stress on Leaf Gas Exchange and Chlorophyll Fluorescence of Glycyrrhiza uralensis. Rus. J. Ecol. 2014, 46, 532–538. [Google Scholar] [CrossRef]
  40. Zhang, S.Y.; Zhang, G.C.; Gu, S.Y.; Xia, J.B.; Zhao, J.K. Critical responses of photosynthetic efficiency of goldspur apple tree to soil water variation in semiarid loess hilly area. Photosynthetica 2010, 48, 589–595. [Google Scholar] [CrossRef]
  41. Lang, Y.; Wang, M.; Xia, J.; Zhao, Q. Effects of soil drought stress on photosynthetic gas exchange traits and chlorophyll fluorescence in Forsythia suspensa. For. Res. 2018, 29, 45–53. [Google Scholar] [CrossRef]
  42. Ohashi, Y.; Nakayama, N.; Saneokal, H.; Fujita, K. Effects of drought stress on photosynthetic gas exchange, chlorophyll fluorescence and stem diameter of soybean plants. Biol. Plant. 2006, 50, 138–141. [Google Scholar] [CrossRef]
  43. Pukacki, P.M.; Kamińska-Rożek, E. Effect of drought stress on chlorophyll a fluorescence and electrical admittance of shoots in Norway spruce seedlings. Trees 2005, 19, 539–544. [Google Scholar] [CrossRef]
  44. Guo, P.; Baum, M.; Varshney, R.K.; Graner, A.; Grando, S.; Ceccarelli, S. QTLs for chlorophyll and chlorophyll fluorescence parameters in barley under post-flowering drought. Euphytica 2008, 163, 203–214. [Google Scholar] [CrossRef]
  45. Longenberger, P.S.; Smith, C.W.; Duke, S.E.; McMichael, B.L. Evaluation of chlorophyll fluorescence as a tool for the identification of drought tolerance in upland cotton. Euphytica 2009, 66, 25–33. [Google Scholar] [CrossRef]
  46. Hura, T.; Grzesiak, S.; Hura, K.; Thiemt, E.; Tokarz, K.; Wędzony, M. Physiological and biochemical tools useful in drought-tolerance detection in genotypes of winter triticale: Accumulation of ferulic acid correlates with drought tolerance. Ann. Bot. 2007, 100, 767–775. [Google Scholar] [CrossRef] [PubMed]
  47. Razavi, F.; Pollet, B.; Steppe, K.; Labeke, M.C. Chlorophyll fluorescence as a tool for evaluation of drought stress in strawberry. Photosynthetica 2008, 46, 631–633. [Google Scholar] [CrossRef]
  48. Perlikowski, D.; Kosmala, A.; Rapacz, M.; Kościelniak, J.; Pawłowicz, I.; Zwierzykowski, Z. Influence of short-term drought conditions and subsequent re-watering on the physiology and proteome of Lolium multiflorum/Festuca arundinacea introgression forms with contrasting levels of tolerance to long-term drought. Plant Biol. 2014, 16, 385–394. [Google Scholar] [CrossRef] [PubMed]
  49. Slabbert, R.M.; Krüger, G.H.J. Assessment of changes in photosystem II structure and function as affected by water deficit in Amaranthus hypochondriacus L. and Amaranthus hybridus L. Plant Physiol. Biochem. 2011, 49, 978–984. [Google Scholar] [CrossRef]
  50. Oukarroum, A.; Bussotti, F.; Goltsev, V.; Kalaji, H.M. Correlation between reactive oxygen species production and photochemistry of photosystems I and II in Lemna gibba L. plants under salt stress. Environ. Exp. Bot. 2015, 109, 80–88. [Google Scholar] [CrossRef]
  51. Baker, N.R. Chlorophyll fluorescence: A probe of photosynthesis in vivo. Annu. Rev. Plant. Biol. 2008, 59, 89–113. [Google Scholar] [CrossRef]
  52. Gao, J.; Li, P.; Ma, F.; Goltsev, V. Photosynthetic performance during leaf expansion in Malus micromalus probed by chlorophyll a fluorescence and modulated 820 nm reflection. J. Photochem. Photobiol. B Biol. 2014, 137, 144–150. [Google Scholar] [CrossRef]
  53. Oukarroum, A.; Strasser, R.J.; Schansker, G. Heat stress and the photosynthetic electron transport chain of the lichen Parmelina tiliacea (Hoffm.) Ach. in the dry and the wet state: Differences and similarities with the heat stress response of higher plants. Photosynth. Res. 2012, 111, 303–314. [Google Scholar] [CrossRef] [PubMed]
  54. Hamdani, S.; Qu, M.; Xin, C.-P.; Li, M.; Chu, C.; Govindjee; Zhu, X.-G. Variations between the photosynthetic properties of elite and landrace Chinese rice cultivars revealed by simultaneous measurements of 820 nm transmission signal and chlorophyll a fluorescence induction. J. Plant Physiol. 2015, 177, 128–138. [Google Scholar] [CrossRef] [PubMed]
  55. Goltsev, V.; Zaharieva, I.; Chernev, P.; Kouzmanova, M.; Kalaji, H.M.; Yordanov, I.; Krasteva, V.; Alexandrov, V.; Stefanov, D.; Allakhverdiev, S.I.; et al. Drought-induced modifications of photosynthetic electron transport in intact leaves: Analysis and use of neural networks as a tool for a rapid non-invasive estimation. Biochim. Biophys. Acta 2012, 1817, 1490–1498. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Salvatori, E.; Fusaro, L.; Gottardini, E.; Pollastrini, M.; Goltsev, V.; Strasser, R.J.; Bussotti, F. Plant stress analysis: Application of prompt, delayed chlorophyll fluorescence and 820 nm modulated reflectance. Insights from independent experiments. Plant Physiol. Biochem. 2014, 85, 105–113. [Google Scholar] [CrossRef]
  57. Zhou, R.; Kan, X.; Chen, J.; Hua, H.; Li, Y.; Ren, J.; Feng, K.; Liu, H.; Deng, D.; Yin, Z. Drought-induced changes in photosynthetic electron transport in maize probed by prompt fluorescence, delayed fluorescence, P700 and cyclic electron flow signals. Environ. Exp. Bot. 2019, 158, 51–62. [Google Scholar] [CrossRef]
Figure 1. Induction curves of chlorophyll a fluorescence of perennial ryegrass varieties (Roadrunner and Nira) grown under different field water capacities of soil (50 and 30% field water content (FWC)) at 0, 120, 240, and 400 h after stress application. (A) Roadrunner under 50% FWC, (B) Roadrunner under 30% FWC, (C) Nira under 50% FWC, (D) Nira under 30% FWC. n = 6 relative units. rel. u., relative units.
Figure 1. Induction curves of chlorophyll a fluorescence of perennial ryegrass varieties (Roadrunner and Nira) grown under different field water capacities of soil (50 and 30% field water content (FWC)) at 0, 120, 240, and 400 h after stress application. (A) Roadrunner under 50% FWC, (B) Roadrunner under 30% FWC, (C) Nira under 50% FWC, (D) Nira under 30% FWC. n = 6 relative units. rel. u., relative units.
Sensors 19 02736 g001
Figure 2. Differential curves of ΔVt (obtained by subtracting the control curve from the first sample) of perennial ryegrass varieties (Roadrunner and Nira) under different field water capacities of soil (50 and 30% FWC) and times (0, 120, 240, and 400 h after stress application): (A) Roadrunner under 50% FWC, (B) Roadrunner under 30% FWC, (C) Nira under 50% FWC, (D) Nira under 30% FWC. n = 6 relative units.
Figure 2. Differential curves of ΔVt (obtained by subtracting the control curve from the first sample) of perennial ryegrass varieties (Roadrunner and Nira) under different field water capacities of soil (50 and 30% FWC) and times (0, 120, 240, and 400 h after stress application): (A) Roadrunner under 50% FWC, (B) Roadrunner under 30% FWC, (C) Nira under 50% FWC, (D) Nira under 30% FWC. n = 6 relative units.
Sensors 19 02736 g002
Figure 3. JIP-test parameters normalized to the values before stress application (0 h) as radar plots under different field water capacities of soil (50 and 30% FWC). (A) Roadrunner under 50% FWC, (B) Roadrunner under 30% FWC, (C) Nira under 50% FWC, (D) Nira under 30% FWC. Relative units. Means of one parameter marked by an asterisk differ significantly (p < 0.05, n = 6).
Figure 3. JIP-test parameters normalized to the values before stress application (0 h) as radar plots under different field water capacities of soil (50 and 30% FWC). (A) Roadrunner under 50% FWC, (B) Roadrunner under 30% FWC, (C) Nira under 50% FWC, (D) Nira under 30% FWC. Relative units. Means of one parameter marked by an asterisk differ significantly (p < 0.05, n = 6).
Sensors 19 02736 g003
Figure 4. Delayed fluorescence induction curves of perennial ryegrass varieties (Roadrunner and Nira) under different field water capacities of soil (50 and 30% FWC) and times (0, 120, 240, and 400 h after stress application): (A) Roadrunner under 50% FWC, (B) Roadrunner under 30% FWC, (C) Nira under 50% FWC, (D) Nira under 30% FWC. n = 6. Relative units.
Figure 4. Delayed fluorescence induction curves of perennial ryegrass varieties (Roadrunner and Nira) under different field water capacities of soil (50 and 30% FWC) and times (0, 120, 240, and 400 h after stress application): (A) Roadrunner under 50% FWC, (B) Roadrunner under 30% FWC, (C) Nira under 50% FWC, (D) Nira under 30% FWC. n = 6. Relative units.
Sensors 19 02736 g004
Figure 5. (I1-D2)/D2 and I1/I2 ratios ± S.D. of perennial ryegrass varieties (Roadrunner and Nira) under different field water capacities of soil (50 and 30% FWC) and times (0, 120, 240, and 400 h after stress application): (A) Roadrunner under 50% FWC, (B) Roadrunner under 30% FWC, (C) Nira under 50% FWC, (D) Nira under 30% FWC. n = 6. Relative units. Mean values within a term marked by the same letter did not show significant differences (p < 0.05, n = 6).
Figure 5. (I1-D2)/D2 and I1/I2 ratios ± S.D. of perennial ryegrass varieties (Roadrunner and Nira) under different field water capacities of soil (50 and 30% FWC) and times (0, 120, 240, and 400 h after stress application): (A) Roadrunner under 50% FWC, (B) Roadrunner under 30% FWC, (C) Nira under 50% FWC, (D) Nira under 30% FWC. n = 6. Relative units. Mean values within a term marked by the same letter did not show significant differences (p < 0.05, n = 6).
Sensors 19 02736 g005
Figure 6. Kinetics of modulated light reflection at 820 nm (all data are normalized to the initial measured value of the signal) of perennial ryegrass varieties (Roadrunner and Nira) under different field water capacities of soil (50 and 30% FWC) and times (0, 120, 240, and 400 h after stress application): (A) Roadrunner under 50% FWC, (B) Roadrunner under 30% FWC, (C) Nira under 50% FWC, (D) Nira under 30% FWC. n = 6 relative units.
Figure 6. Kinetics of modulated light reflection at 820 nm (all data are normalized to the initial measured value of the signal) of perennial ryegrass varieties (Roadrunner and Nira) under different field water capacities of soil (50 and 30% FWC) and times (0, 120, 240, and 400 h after stress application): (A) Roadrunner under 50% FWC, (B) Roadrunner under 30% FWC, (C) Nira under 50% FWC, (D) Nira under 30% FWC. n = 6 relative units.
Sensors 19 02736 g006
Figure 7. Relative rates of reaction centers of PS I (P700) oxidation (A,C) and P700+ re-reduction (B,D) ± S.D. calculated from the MR820 signal of Roadrunner (A,B) and Nira (C,D) under different field water capacities of soil (50 and 30% FWC) and times (0, 120, 240, and 400 h after stress application). Relative units. Mean values within a term marked by the same letter did not show significant differences (p < 0.05, n = 6).
Figure 7. Relative rates of reaction centers of PS I (P700) oxidation (A,C) and P700+ re-reduction (B,D) ± S.D. calculated from the MR820 signal of Roadrunner (A,B) and Nira (C,D) under different field water capacities of soil (50 and 30% FWC) and times (0, 120, 240, and 400 h after stress application). Relative units. Mean values within a term marked by the same letter did not show significant differences (p < 0.05, n = 6).
Sensors 19 02736 g007
Figure 8. Gas exchange parameters ± S.D. of Roadrunner and Nira under different field water capacities (70, 50, and 30% FWC) and times after the application of stress (0, 120, 240, and 400 h). (A) CO2 assimilation (A) in Roadrunner; (B) CO2 assimilation (A) in Nira; (C) H2O transpiration (E) in Roadrunner; (D) H2O transpiration (E) in Nira; (E) stomatal conductance (gs) in Roadrunner; (F) stomatal conductance (Gs) in Nira; (G) internal CO2 concentration (Ci) in Roadrunner; (H) internal CO2 concentration (Ci) in Nira. Mean values within a term marked by the same letter did not show significant differences (p < 0.05, n = 6).
Figure 8. Gas exchange parameters ± S.D. of Roadrunner and Nira under different field water capacities (70, 50, and 30% FWC) and times after the application of stress (0, 120, 240, and 400 h). (A) CO2 assimilation (A) in Roadrunner; (B) CO2 assimilation (A) in Nira; (C) H2O transpiration (E) in Roadrunner; (D) H2O transpiration (E) in Nira; (E) stomatal conductance (gs) in Roadrunner; (F) stomatal conductance (Gs) in Nira; (G) internal CO2 concentration (Ci) in Roadrunner; (H) internal CO2 concentration (Ci) in Nira. Mean values within a term marked by the same letter did not show significant differences (p < 0.05, n = 6).
Sensors 19 02736 g008
Figure 9. Pearson correlation coefficients (r) between prompt fluorescence parameters and CO2 assimilation (A) (correlations significant at p < 0.05, n = 6). Relations between parameters marked by a solid line were calculated for plants in control conditions, and those marked by a dotted line were calculated for plants under 30% FWC.
Figure 9. Pearson correlation coefficients (r) between prompt fluorescence parameters and CO2 assimilation (A) (correlations significant at p < 0.05, n = 6). Relations between parameters marked by a solid line were calculated for plants in control conditions, and those marked by a dotted line were calculated for plants under 30% FWC.
Sensors 19 02736 g009
Figure 10. Pearson correlation coefficients (r) between delayed fluorescence parameters and CO2 assimilation (A) and between modulated reflection at 820 nm and CO2 assimilation (A) (correlations significant at p < 0.05, n = 6). Relations between parameters marked by a solid line were calculated for plants in control conditions, and those marked by a dotted line were calculated for plants under 30% FWC.
Figure 10. Pearson correlation coefficients (r) between delayed fluorescence parameters and CO2 assimilation (A) and between modulated reflection at 820 nm and CO2 assimilation (A) (correlations significant at p < 0.05, n = 6). Relations between parameters marked by a solid line were calculated for plants in control conditions, and those marked by a dotted line were calculated for plants under 30% FWC.
Sensors 19 02736 g010
Table 1. The description of fluorescence parameters (after Strasser et al. 2010 [33]), delayed fluorescence (DF), MR820, and gas exchange parameters. All fluorescence parameters are in relative units. RC, reaction center.
Table 1. The description of fluorescence parameters (after Strasser et al. 2010 [33]), delayed fluorescence (DF), MR820, and gas exchange parameters. All fluorescence parameters are in relative units. RC, reaction center.
VJ = (FJ − F0)/(FM − F0)Relative variable fluorescence at the J-step
φPo = 1 − F0/FM Maximum quantum yield of primary photochemistry (at t = 0)
φEo = (1 − F0/FM)(1 − VJ) Quantum yield of electron transport (at t = 0)
φRo = (1 − F0/FM)(1 − VI) Quantum yield of reduction of end electron acceptors at the PSI acceptor side (reaction center-RE)
φDo = F0/FM Quantum yield (at t = 0) of energy dissipation
ΨEo = 1 − VJ Probability (at t = 0) that a trapped exciton moves an electron into the electron transport chain beyond Q A
δRo = (1 − VI)/(1 − VJ) Efficiency/probability with which an electron from the intersystem electron carriers moves to reduce end electron acceptors at the PSI acceptor side (RE)
t(FM) Time (in m/s) to reach the maximal fluorescence intensity FM
PIABS = γRC/(1 − γRC) × φPo/(1 − φPo) × ΨEo/(1 − ΨEo) Performance index (potential) for energy conservation from exciton to the reduction of intersystem electron acceptors
PItotal = PIABS × δRo/(1 − δRo) Performance index (potential) for energy conservation from exciton to the reduction of PSI end acceptors
ABS/RC = (1 − γRC)/γRC Absorption flux (of antenna Chls) per RC
M0 Approximated initial slope (in ms−1) of the fluorescence transient V = f(t)
TR0/RC = M0(1/VJ) Trapping flux (leading to QA reduction) per RC
ET0/RC = M0(1/VJ0 Electron transport flux (further than Q A) per RC
RE0/RC = M0(1/VJ)(1 − VI) Electron flux reducing end electron acceptors at the PSI acceptor side per RC
DI0/RC = (ABS/RC − TR0/RC) Dissipated energy flux per RC (at t = 0)
RC/CS0 = φPo (VJ/M0) F0 Density of RCs (QA reducing PSII reaction centers)
I1 and I2Maximum of the DF induction curve
D2Minimum of the DF induction curve
MR0 Modulated 820-nm reflection intensity at Time “0”
MRmin Minimum of modulated 820-nm reflection intensity
MRmax Maximum of modulated 820-nm reflection intensity
ΔMRfastFast phase (oxidation) of reflection intensity = MR0 − MRmin
ΔMRslowSlow phase (reduction) of reflection intensity = MRmax − MRmin
ACO2 assimilation (µmol CO2 m−2·s−1)
EH20 transpiration rate (mmol H2O m−2·s−1)
GsStomatal conductance (mol H2O m−2·s−1)
CiSub-stomatal CO2 concentration (µmol CO2 mol1 air)

Share and Cite

MDPI and ACS Style

Dąbrowski, P.; Baczewska-Dąbrowska, A.H.; Kalaji, H.M.; Goltsev, V.; Paunov, M.; Rapacz, M.; Wójcik-Jagła, M.; Pawluśkiewicz, B.; Bąba, W.; Brestic, M. Exploration of Chlorophyll a Fluorescence and Plant Gas Exchange Parameters as Indicators of Drought Tolerance in Perennial Ryegrass. Sensors 2019, 19, 2736. https://doi.org/10.3390/s19122736

AMA Style

Dąbrowski P, Baczewska-Dąbrowska AH, Kalaji HM, Goltsev V, Paunov M, Rapacz M, Wójcik-Jagła M, Pawluśkiewicz B, Bąba W, Brestic M. Exploration of Chlorophyll a Fluorescence and Plant Gas Exchange Parameters as Indicators of Drought Tolerance in Perennial Ryegrass. Sensors. 2019; 19(12):2736. https://doi.org/10.3390/s19122736

Chicago/Turabian Style

Dąbrowski, Piotr, Aneta H. Baczewska-Dąbrowska, Hazem M. Kalaji, Vasilij Goltsev, Momchil Paunov, Marcin Rapacz, Magdalena Wójcik-Jagła, Bogumiła Pawluśkiewicz, Wojciech Bąba, and Marian Brestic. 2019. "Exploration of Chlorophyll a Fluorescence and Plant Gas Exchange Parameters as Indicators of Drought Tolerance in Perennial Ryegrass" Sensors 19, no. 12: 2736. https://doi.org/10.3390/s19122736

APA Style

Dąbrowski, P., Baczewska-Dąbrowska, A. H., Kalaji, H. M., Goltsev, V., Paunov, M., Rapacz, M., Wójcik-Jagła, M., Pawluśkiewicz, B., Bąba, W., & Brestic, M. (2019). Exploration of Chlorophyll a Fluorescence and Plant Gas Exchange Parameters as Indicators of Drought Tolerance in Perennial Ryegrass. Sensors, 19(12), 2736. https://doi.org/10.3390/s19122736

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