- Methodology article
- Open Access
A procedure for maize genotypes discrimination to drought by chlorophyll fluorescence imaging rapid light curves
Plant Methodsvolume 13, Article number: 61 (2017)
Photosynthesis can be roughly separated into biochemical and photochemical processes. Both are affected by drought and can be assessed by non-invasive standard methods. Gas exchange, which mainly assesses the first process, has well-defined protocols. It is considered a standard method for evaluation of plant responses to drought. Under such stress, assessment of photochemical apparatus by chlorophyll fluorescence needs improvement to become faster and reproducible, especially in growing plants under field conditions. For this, we developed a protocol based on chlorophyll fluorescence imaging, using a rapid light curve approach.
Almost all parameters obtained by rapid light curves have shown statistical differences between control and drought stressed maize plants. However, most of them were affected by induction processes, relaxation rate, and/or differences in chlorophyll content; while they all were influenced by actinic light intensity on each light step of light curve. Only the normalized parameters related to photochemical and non-photochemical quenching were strongly correlated with data obtained by gas exchange, but only from the light step in which the linear electron flow reached saturation.
The procedure developed in this study for discrimination of plant responses to water deficit stress proved to be as fast, efficient and reliable as the standard technique of gas exchange in order to discriminate the responses of maize genotypes to drought. However, unlike that, there is no need to perform daily and time consuming calibration routines. Moreover, plant acclimation to the dark is not required. The protocol can be applied to plants growing in both controlled conditions and full sunlight in the field. In addition, it generates parameters in a fast and accurate measurement process, which enables evaluating several plants in a short period of time.
In recent years, it has been reported an increase in the frequency and intensity of extreme weather events such as cold, heat, flooding and drought especially in those regions with intensive plant cultivation and large production of food in the world [1, 2]. Such adverse weather events have potential to greatly reduce the crops productivity and, consequently, jeopardize the food supply worldwide . Whenever the environmental variables undergo changes, moving away from a range considered ideal for plant cultivation towards an extreme limit, they creates the conditions for the occurrence of what is known as abiotic stresses . Abiotic stresses can be defined as environmental conditions that reduce growth and yield below optimum levels . Currently, the water availability in the soil is the variable that contributes the most to stress by drought, among others environmental variables with the greatest potential to cause stress in plants .
In order to deal with any stressful situation, plants trigger several mechanisms that work in all the levels of organization . Under drought, it is widely accepted that one of the primary effects on plants is a reduction on the stomatal opening, with a consequent decrease in the water vapor and CO2 conductance [8, 9]. This mechanism reduces the water vapor loss, but, as a side effect, it also restricts CO2 entry in the leaf mesophyll . Therefore, the measurement of gas exchange is considered a standard technique for studies related to this kind of stress. Consequently, the variables derived from the measurement of gas exchange, especially CO2 assimilation rate, stomata conductance and transpiration have been the most used to discriminate the plant responses to stress by water deficit .
Despite the fact that its initial manifestation is mainly perceived on gas exchange, stress by water deficit is one of the most comprehensive stresses regarding the effects on plant metabolism. It results from the fact that, in general, water shortage itself affects several processes, and it is usually associated with an increase in temperature and light intensity , enhancing the effects of drought mainly in tropical environments. Consequently, parallel to the effect on the gas exchange and the photosynthetic metabolism, damages to the structures and consequently to the photochemical processes occurring inside the chloroplasts are enhanced .
Drought not only causes degradation of photosynthetic pigments, but also disorganization of the thylakoid membranes [13, 14]. In such membranes are anchored the components of the primary photochemistry of photosynthesis, represented mainly by the light-harvesting complexes and their respective reaction centers of photosystems I and II . The latter is responsible for using the energy of the absorbed light to drive the electron flow through the chloroplast membrane system [16,17,18].
The loss of pigments reduces the absorption of light energy and, in the first instance, may even attenuate the potential damages caused by the drought on the photochemical apparatus . However, the disorganization of the thylakoid membranes causes damage to the reaction centers of photosystem II , which itself impairs the primary reactions of light energy conversion to chemical energy . Due to the restriction of CO2 entry in chloroplasts during drought, there is a decrease in primary photochemical processes . In order to avoid damage to the photochemical apparatus, the captured light energy is partially diverted to different processes [21,22,23].
The light energy absorbed by chlorophylls associated with PSII can be used to drive photochemistry, lost from PSII as heat or emitted as chlorophyll fluorescence [13, 14, 17, 24, 25]. The chlorophyll fluorescence parameters that can be used to estimate the flux of excitation energy into three competing pathways are Y(II), Y(NPQ) and Y(NO) . In the first two cases the emission of fluorescence does not occur, that is, it is quenched, and the resulting processes are called photochemical and non-photochemical quenching of chlorophyll fluorescence, respectively. Separation of fluorescence quenching into photochemical and non-photochemical components can be reached by chlorophyll fluorescence technique, using the saturation pulse method . This method basically measures the fluorescence signal emitted under dark (Fo, Fm) or actinic illumination (Fs, Fm′) .
From the saturation pulse method, and considering the model of antenna pigment organization in the thylakoids as a puddle or a lake , several parameters were developed which correlate with photochemical and non-photochemical quenching [26, 27]. Thus, Y(II), qP, and Fv/Fm are considered as photochemical quenching parameters, while NPQ is non-photochemical . qL and qN represent respectively the photochemical and non-photochemical quenching in the lake model, while qP represents photochemical quenching in the puddle model . For all this, chlorophyll fluorescence is considered as a probe of the photochemical apparatus . In addition, it has been proven to be an efficient and reproducible method for evaluating plants under stress .
The measurements obtained by the chlorophyll fluorescence technique can be rapid, highly sensitive and non-invasive [17, 25, 30]. They can be performed on intact leaves still attached to plants. This technique has been used to discriminate the responses to drought in controlled conditions mainly from induction curves [31,32,33] since the simple measurement of Fv/Fm seems to work only in cases of severe drought [34, 35]. However, to evaluate a large number of plants under real environmental conditions in the field, induction curves (IC) would not be the most recommended approach. This is because IC is a procedure in which plant photosynthetic apparatus, after a previous period of dark adaptation, sufficient to completely oxidize the PSII reaction center, is re-submitted to illumination . The protocols currently employed for IC’s running allows the correct measurement of Fo and Fm, which are fundamental for the calculation of Fv/Fm and all parameters related to photochemical and non-photochemical fluorescence quenching [17, 25]. However, IC’s require a prior period of darkness of at least 15 min and the measurement process itself takes longer, i.e., about 5 min per sample. To overcome these challenges, especially when using chlorophyll fluorescence technology for generating images of increasingly larger areas in high throughput systems, an increasing trend nowadays [37,38,39], researchers have used new approaches that work very well on model plants under highly controlled conditions [40, 41]. Nonetheless, such approaches do not work for large plants grown in the field. In this case, to get a fast image of chlorophyll fluorescence parameters, rapid light curves (RLC’s) can be considered as an option.
RLC’s capture the instantaneous responses of plants in a considerable range of light intensity by the simultaneous evaluation of a number of parameters. The biggest advantage of this approach is that it does not demand acclimatization to dark. Moreover, it can be applied to plants grown in full sunlight which saves time and increases the measuring range. However, if the RLC’s are performed without pre-dark acclimation, the initial fluorescence measurements obtained do not correspond to the theoretical values of Fo and Fm. Consequently, the calculated fluorescence quenching does not match to actual values. Therefore, RLC’s are usually used to determine the maximum electron transport rate as well as other parameters [42, 43]. Nevertheless, if the aim of the study does not depend on an accurate measure of Fo, Fm and photochemical and non-photochemical quenching, such parameters obtained from RLC’s without prior dark-acclimation can be used . Previous studies using RLC’s approach for discrimination of plant responses to stresses only compared chlorophyll fluorescence parameters obtained along the light curve [13, 42, 43].
In this study, we used adult maize plants of different genotypes grown under field conditions which were subjected to drought in the pre-flowering stage. Since the cultivation conditions restricted the approaches options for evaluation by chlorophyll fluorescence, we applied upward RLC’s on plants without pre-dark-acclimation to discriminate the genotype responses to drought. The RLC’s were dissected in its components light steps and the measured chlorophyll fluorescence parameters in each of them were correlated with gas exchange. Based on the induction effects of the light curve itself and the effects of drought on the photochemical apparatus of plants, we argued that not all measured chlorophyll fluorescence parameters nor all light steps can be used for correct discrimination of maize genotypes to drought. Only normalized photochemical and non-photochemical parameters, measured in the light step in which the linear electron flow reached saturation, can be indicated for this. All of them were as efficient and reliable as the standard technique of gas exchange in order to discriminate the responses of maize genotypes to stress by water deficit. From this study, it was proposed a procedure to be used for this aim, which standardizes the type and portion of the leaf to be evaluated, the RLC’s parameters configuration and the ideal light step in which the most appropriate parameters potentiate the differences.
Changes in chlorophyll fluorescence parameters over RLC’s and chlorophyll content in control and drought stressed maize plants
All chlorophyll fluorescence parameters obtained throughout the RLC’s showed statistically significant differences between control and drought stressed maize plants. For most parameters, the biggest differences between stressed and not stressed plants occurred from the middle to the end portion of the RLC’s, except for the fluorescence yield (F) and for the quantum yield of nonregulated energy dissipation [Y(NO)] (Fig. 1; Additional files 1, 2).
The RLC’s relative to F and Y(NO) showed similar shape for stressed and non-stressed plants (Fig. 1a, b). The means of these two parameters, measured without actinic illumination (light step zero), were higher in plants subjected to drought in comparison to control. After starting actinic illumination and a gradual increase in light levels, both parameters increased for both groups of plants. However, the rises occurred only in the first step of actinic light in drought stressed plants and up to the fourth light step in the control ones. Then, F and Y(NO) began to decline in a faster and more intense manner in plants under stress. At the end of the RLC’s, the means of both parameters in stressed plants become smaller than those of control plants.
The RLC’s representing parameters known as maximum fluorescence yield in the light (Fm′) and effective PSII quantum yield [Y(II)] also showed a very similar pattern to each other (Fig. 1c, d). Both parameters showed an expected trend of decreasing the extent to which actinic light has increased. However, the control plants showed higher Fm′ and Y(II) means throughout all the RLC’s steps. Similarly to what had happened to F and Y(NO), the decline in Fm′ and Y(II) means was much higher and faster in plants subjected to drought stress compared to the control ones. The major differences in Fm′ and Y(II) means between stressed and control plants occurred at the intermediate portion of the RLC’s, whereas in extreme ends there was a tendency of shortening of these differences.
The parameters representing regulated dissipation of light energy as heat [qN, Y(NPQ) and NPQ] (Fig. 1e–g) increased along the RLC’s for both groups of maize plants, stressed and non-stressed. Nevertheless, the increase was faster and more robust on drought stressed plants, so that the means of all those parameters remained higher in this group of plants compared to controls. On the other hand, the parameters related to photochemical dissipation of light energy [Y(II), qP and qL] (Fig. 1d, h, i) decreased. For these parameters, the decrease was more pronounced in stressed plants.
Instead of ETR, we measured a similar parameter called the apparent rate of photosynthesis (PS). While for ETR a PAR-absorptivity of 0.84 is assumed, for PS the PAR-absorptivity is determined using an IMAGING-PAM routine. Calculated in this way, the PS values reached around 50. Then, to allow them to be visualized by the color bar, whose range varies from 0 to 1, the PS values were divided by 50, hence PS/50. As expected, the apparent rate of photosynthesis (PS/50) went up with the increase in the intensity of the actinic light up to saturation level. The treatments by which the plants were subjected have influenced the intensity and speed of the PS/50 increase as well as the actinic light intensity in which saturation occurred (Fig. 1j). It was observed that, in stressed plants, PS/50 increased to a less extent degree and saturation occurred at a lower light level, around 185 µmol m−2 s−1, while on the control plants it increased in a more pronounced way and saturation occurred at 280 µmol m−2 s−1.
Since the non-normalized chlorophyll fluorescence parameters, such as F and Fm′, are affected by leaf chlorophyll content, the chlorophyll content index (CCI) was measured in the leaves of control and drought stressed maize plants throughout the experiment. This parameter did not change in the leaves of control plants over the experimental period, while it was observed a fall in drought stressed plants (Fig. 2; Additional file 3).
Relationship among RLC’s chlorophyll fluorescence parameters and gas exchange
Most of the parameters obtained from every RLC’s light step showed statistically significant correlation with those of gas exchange (Table 1). However, in the initial steps of actinic lighting, the correlation coefficient for most parameters lies in the range considered as from moderate to weak. As the intensity of the actinic light increased, it was seen an increase in the correlation coefficient until 185 µmol m−2 s−1 for most parameters, except for F and Y(NO). For these parameters, the correlation decreased from the initial to the middle portion of the light curve and increased again afterward.
Despite the level of statistical significance in all light steps, the low coefficient correlations for F and Y(NO), preponderantly, are considered as moderate to weak. On the other hand, Fm′, Y(II), PS/50, qP and qL as well as qN, NPQ and Y(NPQ) showed strong correlations with the parameters of gas exchange, except Ci. On average, the correlation coefficient for Y(II) and its derived PS/50 were slightly higher than for the other parameters.
Since the correlation strength between the most gas exchange and chlorophyll fluorescence parameters was similar, with exception of Ci, we chose the rate of CO2 assimilation (A) to graph it. Thus, a scatter plot showing the strength of the correlation between A and the fluorescence parameters in the light step in which there was the greatest difference between stressed and control plants, i.e. at 280 µmol m−2 s−1, is shown (Fig. 3).
Discrimination of maize genotypes responses to drought by chlorophyll fluorescence RLC’s
It is known that the changes caused by drought in the parameters of leaf gas exchange occur in parallel and interrelated to changes in chlorophyll fluorescence. Therefore, we hypothesized that these changes could be followed by the chlorophyll fluorescence parameters obtained from rapid light curves.
The responses of the maize plants submitted to drought occurred as expected, based on the main gas exchange parameters chosen. There was a sharp drop in A, gs, and E, and an increase in Ci over 7 days after withholding (Fig. 4, Additional files 4, 5). To show that the chlorophyll fluorescence parameters follow a similar dynamic, we chose the value of each of them obtained in the light step at 280 µmol m−2 s−1 to plot over time (Fig. 5; Additional files 6, 7). Such a choice was made because it was one of the light steps less affected by induction processes. In this light step, PS/50 saturation had already been achieved on both control and stressed plants. Moreover, it was the light step in which the most parameters showed the greatest differences between stressed and non-stressed plants in the shortest time from the beginning of the RLC’s, around 60 s (Fig. 1).
Similar to what had already occurred with the gas exchange parameters, it was observed reasonably unchanged in all chlorophyll fluorescence parameters in control plants over seven days. Also, as expected, there was a decrease in parameters related to photochemical dissipation of light energy [Y(II), qP, qL and PS/50] (Fig. 5d, h–j) consistent the extent to which the plants were entering the drought stress. The same applies to the non-normalized parameters related to fluorescence emission (F and Fm′) (Fig. 5a, c). On the other hand, the parameters related to the heat dissipation of light energy [qN, Y(NPQ) and NPQ] (Fig. 5e–g) simultaneously increased. Meanwhile, the parameter related to unregulated dissipation of light energy [Y(NO)] (Fig. 5b) virtually unchanged due to drought. At the end of the 7th day after water withholding, most parameters showed statistical differences between control and drought stressed plants which can be viewed through the images (Fig. 6; Additional files 8, 9, 10).
Changes in chlorophyll fluorescence parameters over RLC’s in control and drought stressed maize plants
Clearly, there was a pattern of responses for each chlorophyll fluorescence parameter gotten along the RLC’s. This pattern repeated itself in control plants as well as in those subjected to stress by water deficit in all maize genotypes (Fig. 1; Additional file 1). It is noteworthy the initial increase followed by a reduction in F and Y(NO), both in stressed and non-stressed plants. Probably the kinetics of these parameters was influenced by the induction processes that act on the photochemical apparatus, which are affected by the procedures adopted to carry the RLC’s out. Such induction processes tend to occur, even if it takes a period of time as short as just 10 s of darkness before starting the RLC’s .
These results demonstrate that both stressed and control plants suffered induction effects; however, two predominantly different and competitive mechanisms were involved. Thus, the kinetic of F and Y(NO) along the RLC’s was possibly determined by the main mechanism of fluorescence quenching that each plant group have preferably used. Thus, the stressed plants, due to the limitations imposed on photosynthesis by drought, diverted a considerable part of the electrons flow to the non-photochemical processes, where induction is faster so that F and Y(NO) decreased sooner in this group of plants. On the other hand, control plants drove the electrons flow preferably for the photochemical process whose induction is slower which delayed the decline in both parameters. At the end of the RLC’s, when both processes probably were fully activated, the differences between F and Y(NO) from stressed and control plants were reduced. Even so, based on the values of these parameters, it was possible to detect statistically significant differences between stressed and control plants on the initial as well as on the end portion of the RLC’s.
The maximum fluorescence emission in the light (Fm′) is an indicative of energy dissipation by non-photochemical quenching processes, which are highly regulated in plant leaves . In this study, Fm′ arbitrary values from stressed plants were lower than those of control throughout the RLC’s. However, the simple comparison of arbitrary Fm′ values of the two groups of plants would not be a good indication of the differences between them. The reason is that Fm′ is influenced by differences in chlorophyll content  and, at least in the beginning of the RLC’s, it also is affected by NPQ relaxation . In this study, both traits were affected by drought (Figs. 1, 2). It would be more appropriate, therefore, to perform an analysis disregarding the arbitrary values of Fm′ and taking into account only the kinetics of it along the curve, which has no relation to the chlorophyll content. By using this approach, the fall in the Fm′ was much bigger and faster in plants subjected to water deficit. It is probably related to the faster start of the NPQ, which has shown much higher rates than in control plants. Similar to what occurred to F, the difference in Fm′ values between stressed and control plants was lower at the extreme ends of the RLC’s. For that reason, the intermediate portion of the RLC’s would be the most recommended part to discriminate the responses of plants to water deficit using that parameter.
The sharp reduction in the parameters related to photochemical quenching [qL, qP and Y(II)], and the concurrent and proportional increase in those related to non-photochemical quenching [qN, NPQ and Y(NPQ)], suggest that maize plants submitted to drought maintained their capacity to regulate dissipation of light energy, despite the stressful condition. This fact, along with the rapid reduction of F in this group of plants, constitutes evidence that the xanthophyll’s cycle was efficient in giving vent to excess electrons flow generated for the increase in light intensity in the chloroplasts. In fact, Saccardy et al.  had shown that when subjected to drought, maize plants responded by an increase in the xanthophyll’s cycle pool size as well as the in the proportion of de-epoxidized (antheraxanthin and zeaxanthin) to epoxidized (violoxanthin) xanthophylls.
Higher pressure on the photochemical apparatus, represented by an increase in light intensity along the RLC’s, had only a limited effect on the enlargement of the difference in Y(II) or PS/50 between control and drought plants. PS/50 reached saturation around 280 µmol m2 s−1 in control plants, while in those stressed saturation occurred in a lower light intensity, around 185 µmol m2 s−1. The increase in the difference between these two groups of plants could, in this case, only occur if Y(II) and, consequently, PS/50 fell to zero in stressed plants. This possibility would depend mainly on the increase in stress intensity and on the genotype sensitivity, but not by an increase in light intensity. In fact, after reaching saturation, an increase in light intensity would have an opposite effect, reducing the difference between these groups of plants, as a consequence of the drop in both Y(II) and PS/50, mainly in control plants. Based on this found, there is no need to move beyond 280 µmol m2 s−1 for the purpose of discriminating the maize genotypes responses to drought by using these two parameters. On the other hand, smaller light intensity values do not express the maximum potential difference that can exist between the two groups of plants.
RLC’s chlorophyll fluorescence parameters enabling discrimination of maize genotypes responses to drought
Based on the results of this study, we found that we must be careful in interpretation and using parameters obtained by RLC’s to discriminate plant responses to drought. The use of such parameters for that purpose only because they have shown statistically significant correlation with those of gas exchange (Table 1; Fig. 3), or due to the fact they have presented differences between control and stressed plants (Figs. 1, 5, 6), can lead to errors. It is important that their use for this purpose be supported by the knowledge of what each of them actually represents, and to what trait in the plant it is linked. In addition, one should take into account that secondary factors, which are associated with drought, may interfere with chlorophyll fluorescence parameters obtained by RLC’s.
F, Fm′ and Y(NO) well illustrate this point. The low correlation coefficients for F and Y(NO), both related to light energy dissipation in a non-regulated manner [16, 26], for instance, indicate that this process does not maintain a strong link with the leaf gas exchange. Furthermore, the analysis of the results from non-normalized parameters, such as F and Fm′, becomes complicated. They were affected by the loss of chlorophyll from the leaves throughout the drought period (Fig. 2), with a probable and consequent reduction in the amount of absorbed light . It must also be considered a remarkable effect of induction over these parameters, observed throughout the actinic lighting steps from RLC’s. As a result, depending on the level of light used for comparison, stressed plants can have a smaller or larger values of F and Y(NO) in comparison to controls. This makes these parameters somewhat ambiguous to discriminate plant responses to drought.
The parameters related to the photochemical [Y(II), qP and qL] and non-photochemical [qN, NPQ, and Y(NPQ)] quenching showed, differently of the parameters linked to fluorescence emission, strong correlation with those derived from gas exchange. According to the theory of partitioning of absorbed light energy among photochemistry, fluorescence and heat, Y(II) measures the proportion of the light energy absorbed by PSII which is used in photochemistry [13, 14, 17, 24, 25]. As such, it can give a measure of the rate of linear electron transport from PSII to PSI. Consequently, there is a strong linear relationship between this parameter and quantum efficiency of carbon assimilation .
By limiting photochemistry, as occurred in drought maize plants, Y(II) would decrease in direct proportion of CO2 assimilation, whereas heat dissipation would increase in the inverse proportion. For the first premise to be true, there must be no change in fluorescence yield [Y(NO)], which was actually confirmed (Fig. 5b). In addition, other electron drains [45,46,47] cannot be functional. Indeed, there is no evidence in some studies that other processes such as photorespiration or Mehler reaction may be involved in fluorescence extinction in maize plants under drought stress [44, 46]. And, even though if there may have been other electrons drains, as suggested by others studies on maize under stress , including drought [48, 49], in this study they were not strong enough to prevent the correlation among gas exchange and chlorophyll fluorescence parameters.
It was found that any of these normalized photochemical and non-photochemical parameters obtained by RLC’s were able to discriminate maize plant responses to drought. Therefore, any of them can be used for such purpose. This result is independent on whether or not the measurement was performed using the actual values of chlorophyll fluorescence parameters. However, not all light steps can be recommended for this because the correlation coefficients from all parameters varied along the RLC’s, depending on the intensity of actinic light. They were lower in the initial steps of RLC’s, but increased along with actinic light, reaching the highest values as soon as the light reached saturation, around 280 µmol m2 s−1. From there, any increase in light intensity virtually had no effect on correlation coefficients.
It is known there is a spatial heterogeneity in chlorophyll fluorescence parameters along the leaf blade of plants, regardless of the previous conditions under which they were submitted [16, 25, 36, 50, 51]. This phenomenon was also observed along the maize leaf blade (Fig. 6), which was more pronounced in certain parameters, such as Y(II) (Additional file 9) and PS/50 (Additional file 10). Taking this into account, in addition to the unfeasibility of measuring the whole leaf section, we were careful to choose the area of interest (AOI) of the leaves to be evaluated (see “Methods”). In this same AOI, dozens of data obtained from several replicates were used. Therefore, irrespective of heterogeneity, we showed there was a strong correlation between parameters of gas exchanges with those of chlorophyll fluorescence. Finally, we emphasize that the central point of the manuscript was to show that it is possible to discriminate responses of maize genotypes to drought based on RLC’s chlorophyll fluorescence measurements.
Almost all chlorophyll fluorescence parameters generated by RLC’s showed differences between control and drought maize genotypes. However, not all parameters and nor all light steps could be used for correct discrimination to drought. Most parameters were affected by induction processes, relaxation rate, and/or differences in chlorophyll content; while they all were influenced by actinic light intensity on each light step of light curve. Only normalized parameters relative to photochemical and non-photochemical quenching, measured in the light step in which the linear electron flow reached saturation, were indicated for this. Such measurements were as fast, efficient and reliable as the standard technique of gas exchange in order to discriminate the responses of maize genotypes to stress by water deficit. These findings are important as theses parameters can be used for discrimination of plants subjected to drought both in controlled and field conditions. This may represent an advance in the current trend of large-scale plant phenotyping by imaging systems.
Plant material and growing conditions
The experiments were carried out in a greenhouse at National Agroenergy Research Center (https://www.embrapa.br/en/agroenergia), in Brasilia, Brazil (S—15.732°, W—47.900°) from November 2015 to February 2016. Three weather variables (light intensity, temperature, and humidity) fluctuated accordingly to environmental conditions. Only water supply to the plants and soil moisture content were controlled. Four maize genotypes showing contrasting responses to drought based only on grain yield were used: susceptible (BRS 1010), intermediate (BRS 1030) and tolerant (DKB 390 and 2B 707). They were chosen based on previous results from experiments carried out at National Maize and Sorghum Research Center (https://www.embrapa.br/en/milho-e-sorgo) under field conditions. The seeds were sown in plastic pots (20 kg), filled with typical soil used for maize cultivation in the region (dystrophic Red Latosol according to Brazilian Soil Classification), limed and fertilized based on results from a physical–chemical analysis. After germination, two plants were held per pot, on a daily replenishment of water at field capacity (100% of the available water), based on the weight.
Drought stress and assessed leaf
Stress treatment by drought started about 60 days after sowing, at the pre-flowering stage, when all the plants reached the vegetative growth stage V16 (16 leaf collars emitted). A group of five plants from each genotype was subjected to drought stress by withholding water until the water potential decreased to about −1.5 MPa. They were kept in these conditions for 12 consecutive days, while another group of five plants of the same genotype remained fully irrigated (control). All chlorophyll fluorescence, gas exchange and chlorophyll content index measurements were performed on the adaxial side of a healthy and fully expanded leaf 16, in a previously marked area of approximately 40 cm2 (Additional file 11). This area was located in the middle third of the leaf, at a distance of 25 cm from the apex. It had been chosen in pilot experiments because: (1) it presented highest photosynthetic rates based on gas exchange measurements; (2) under control conditions, it presented lower variation in photosynthetic rates among individuals of the same genotype and even among different genotypes, as can be seen in the data; (3) It is located at a certain distance from the apex, the region most affected by drought stress. Such characteristics are important for studying long-term stresses in which the same attached leaf is evaluated over time.
Chlorophyll fluorescence parameters
Plants submitted to drought stress were assessed by the chlorophyll fluorescence technique (Saturation Pulse Method). For this, chlorophyll a fluorescence was measured using an imaging fluorimeter Walz model IMAGING-PAM Maxi version (Heinz Walz GmbH, Effeltrich, Bayern, Germany) driven by the ImaginWin version 2.40b software. The measuring head consisting of LED-Array Illumination Unit IMAG-MAX/L and a CCD camera IMAG-MAX/K4 was mounted on IMAG-MAX/GS stand. The setup was set on a tripod to make easier the adjustment of the measuring head height for matching the height of the plant to be imaged (Additional file 12). The following settings were used: measuring light = 1; saturation pulse = 10 (2800 µmol m−2 s−1); gain = 1; dumping = 2; red gain = 25; red intensity = 4; NIR intensity = 7; Fm factor = 1.055; F factor = 0.999. For measuring the chlorophyll fluorescence parameters, upward light curves were used. The measurements were performed on individual attached leaves placed under the measuring head and kept in the dark for about 10 s. After that, a light curve with increasing actinic light steps (0, 20, 55, 110, 185, 280, 335, 395, 460) was initiated, using blue light provided by LEDs in 10 s intervals, at the end of which it was applied a saturation pulse. During RLC’s running, a black fabric was used to cover the measuring head in order to avoid external light input at the sample stage. The chlorophyll fluorescence measurements were performed between 9:00 and 12:00 a.m. in a leaf circular area (approximately 12 cm2). This area of interest (AOI) represents the center of the previously marked area in which all measurements were performed, as described in the previous section. It is a circular area that encompasses just half width of the leaf to avoid midrib. The initial values of F and Fm′ obtained in the RLC’s at the light step zero were considered to be Fo and Fm just in order to calculate parameters related to the photochemical and non-photochemical quenching. However, they do not represent the actual values of those parameters, which are obtained when the PSII reaction center is completely oxidized. Also, for the calculation of parameters that depend on it, such as Y(II), the F values were considered as Fs. From the measured parameters, all derived parameters were calculated according to the equations described in the equipment manual, from the ImaginWin version 2.40b software. The Y(II) was determined using the equation: Y(II) = (Fm′ − Fs)/Fm′, as originally described by Genty et al. . Other equations were basically the same described in Van Kooten and Snel  and Maxwell and Johnson  as follow: qP = (Fm′ − Fs)/(Fm′ − Fo′), qN = 1 − (Fm′ − Fo′)/(Fm − Fo), NPQ = (Fm − Fm′)/Fm′. The qL, Y(NPQ) and Y(NO) were determined according to the equations described by Kramer et al. : qL = qP × Fo′/Fs, Y(NPQ) = 1 − Y(II) − 1/(NPQ + 1 + qL(Fm/Fo − 1)), Y(NO) = 1/(NPQ + 1 + qL(Fm/Fo − 1)). Instead of electron transport rate (ETR), the ImaginWin software generates images of a similar parameter denominated PS which is calculated according to the equation: [PS = 0.5 × Y(II) × PAR × Abs.]. In order to display images of this parameter on a false color scale ranging from 0 to 1, the PS value was divided by a number, which correspond to the expected limit of maximal PS. For this study, it was used the standard setting which is 50 and the resulting parameter is called PS/50. This means that the pixel value 1 is reached when PS/50 = 1. The amount of light absorption determined by routine absorptivity was used for the calculation of PS/50. The absorptivity (Abs.) value represents a measure of the fraction of the incident red light which is absorbed by the leaf. It was calculated by the following equation: Abs. = 1 − R/NIR. The values of calculated and measured chlorophyll fluorescence parameters represent the mean values of all pixels within an AOI.
Gas exchange parameters
An LI-COR 6400XT (LI-COR, Lincoln, NE, USA) infra-red gas analyzer, equipped with a size measuring head with 2 × 3 cm and a lighting system artificial LED model 6400-02B was used for measuring of gas exchange. The equipment was configured to maintain the relative humidity within the measuring chamber between 50 and 60% at 30 °C, light intensity in 2000 μmol m−2 s−1 and flow rate at 500 μmol s−1 while collecting data. CO2 concentration was maintained at 400 ppm using a CO2 mixer model 6400-01 from cylinder (12 g) CO2. Each gas exchange measurement in each maize leaf was carried out before chlorophyll fluorescence measurement. By using equations described in the LI-COR 6400XT user manual, the following parameters were obtained by means of software OPEN version 6.3: A = net CO2 assimilation rate (µmol CO2 m−2 s−1), gs = stomatal conductance to water vapor (mol H2O m−2 s−1), E = transpiration rate (mmol H2O m−2 s−1) and C i = intercellular CO2 concentration (µmol CO2 mol air−1).
Chlorophyll content index
Chlorophyll content index (CCI) was measured using a chlorophyll meter Opti-science model CCM-200 Plus (Opti-Sciences Inc., Hudson, NH, USA). In the selected area of the leaf (see Additional file 10) 5 measures were performed. The average of these measures represented the CCI of the leaf.
Experimental design and statistical analysis
A completely randomized design was used to assign two treatments (control × drought) in four maize genotypes (BRS 1030, BRS1010, DKB 390 and 2B 707) and five replicates. Data were analyzed using a one-way analysis of variance (ANOVA). When the treatments or their interactions were significant, the means were compared by Tukey’s test (p ≤ 0.05). Pearson’s product moment correlation coefficients were used to estimate relationships between gas exchange and chlorophyll fluorescence parameters. All the statistical and correlations analyses were performed by using the statistical program STATISTICA version 12 (www.statsift.com, Tulsa, OK, USA).
- Fo, Fo′:
minimum fluorescence yield on dark- and light-adapted leaf, respectively
fluorescence yield (not necessarily in the steady-state) before application a saturate pulse
fluorescence yield on light-adapted leaf (steady-state) before application a saturate pulse
- Fm, Fm′:
maximum fluorescence yield on dark- and light-adapted leaf, respectively
maximum PSII quantum yield
effective PSII quantum yield
quantum yield of regulated energy dissipation
quantum yield of nonregulated energy dissipation
measuring of absorption light by the leaf
coefficient of non-photochemical quenching (lake model)
coefficient of photochemical quenching (puddle model)
coefficient of photochemical quenching (lake model)
rapid light curve
- PSI, PSII:
photosystem I and photosystem II, respectively
electron transport rate
apparent rate of photosynthesis
- A :
net CO2 assimilation rate
- gs :
stomatal conductance to water vapor
- E :
- C i :
intercellular CO2 concentration
chlorophyll content index
Fischer EM, Knutti R. Anthropogenic contribution to global occurrence of heavy-precipitation and high-temperature extremes. Nat Clim Change. 2015;5:560–4.
Mancosu N, Snyder RL, Kyriakakis G, Spano D. Water scarcity and future challenges for food production. Water. 2015;7:975–92.
Food and Agriculture Organization of the United Nations. The state of the world’s land and water resources for food and agriculture (SOLAW)—managing systems at risk. London: Rome and Earthscan; 2011.
Vinebrooke RD, Cottingham KL, Norberg J, Scheffer M, Dodson SI, Maberly SC, Sommer U. Impacts of multiple stressors on biodiversity and ecosystem functioning: the role of species co-tolerance. Oikos. 2004;104:451–7.
Cramer GR, Urano K, Delrot S, Pezzotti M, Shinozaki K. Effects of abiotic stress on plants: a systems biology perspective. BMB Plant Biol. 2011;11:1–14.
Fita A, Rodriguez-Burruezo A, Boscaiu M, Prohens J, Vicente O. Breeding and domesticating crops adapted to drought and salinity: a new paradigm for increasing food production. Front Plant Sci. 2015;6:1–14.
Fraire-Velázquez S, Balderas-Hernández VE. Abiotic stress in plants and metabolic responses. In: Vahdati K, Leslie C, editors. Abiotic stress—plant responses and applications in agriculture. Rijeka: InTech; 2013. p. 25–48.
Chaves MM, Pereira JS, Maroco J, Rodrigues ML, Ricardo CPP, Osorio ML, Carvalho I, Faria T, Pinheiro C. How plants cope with water stress in the field. Photosynthesis and growth. Ann Bot. 2002;89:907–16.
Xu Z, Zhou G. Responses of leaf stomatal density to water status and its relationship with photosynthesis in a grass. J Exp Bot. 2008;59:3317–25.
Lawson T, Blatt MR. Stomatal size, speed, and responsiveness impact on photosynthesis and water use efficiency. Plant Physiol. 2014;164:1556–70.
Osakabe Y, Osakabe K, Shonozaki K, Tran L-SP. Response of plants to water stress. Front Plant Sci. 2014;5:1–8.
Peñuelas J, Munné-Bosch S, Llusià J, Filella I. Leaf reflectance and photo- and antioxidant protection in field-grown summer-stressed Phillyrea angustifolia. Optical signals of oxidative stress? New Phytol. 2004;162:115–24.
Batra NG, Sharma V, Kumari N. Drought-induced changes in chlorophyll fluorescence, photosynthetic pigments, and thylakoid membrane proteins of Vigna radiata. J Plant Interact. 2014;1:712–21.
Ladjal M, Epron D, Ducrey M. Effects of drought preconditioning on thermo tolerance of photosystem II and susceptibility of photosynthesis to heat stress in cedar seedlings. Tree Physiol. 2000;20:1235–41.
Staehelin LA, Arntzen CJ. Regulation of chloroplast membrane function: protein phosphorylation changes the spatial organization of membrane components. J Cell Biol. 1983;97:1327–37.
Murchie EH, Lawson T. Chlorophyll fluorescence analysis: a guide to good practice and understanding some new applications. J Exp Bot. 2013;13:3983–98.
Kalaji HM, Goltsev V, Bosa K, Allakhverdiev SI, Strasser RJ, Govindjee. Experimental in vivo measurements of light emission in plants: a perspective dedicated to David Walker. Photosynth Res. 2012;114:69–96.
Amarnath K, Bennett DIG, Schneider AR, Fleming GR. Multiscale model of light harvesting by photosystem II in plants. Proc Natl Acad Sci. 2016;113:1156–61.
Mafakheri A, Siosemardeh A, Bahramnejad B, Struik PC, Sohrabi Y. Effect of drought stress on yield, proline and chlorophyll contents in three chickpea cultivars. Aust J Crop Sci. 2010;4:580–5.
Goltsev V, Zaharieva I, Chernev P, Kouzmanova M, Kalaji HM, Yordanov I, Krasteva V, Alexandrov V, Stefanov D, Allakhverdiev SI, Strasser RJ. 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–8.
Cousins AB, Adam NR, Wall GW, Kimbal BA, Pinter JR, Ottman MJ, Leavitt SW, Webber AN. Photosystem II energy use, non-photochemical quenching and the xanthophyll cycle in Sorghum bicolor grown under drought and free-air CO2 enrichment (FACE) conditions. Plant Cell Environ. 2002;25:1551–9.
Massacci A, Nabiev SM, Pietrosanti L, Nematov SK, Chernikova TN, Thor K, Leipner J. Response of the photosynthetic apparatus of cotton (Gossypium hirsutum) to the onset of drought stress under field conditions studied by gas-exchange analysis and chlorophyll fluorescence imaging. Plant Physiol Biochem. 2008;46:189–95.
Wei L, Wang L, Yang Y, Wang P, Guo T, Kang G. Abscisic acid enhances tolerance of wheat seedlings to drought and regulates transcript levels of genes encoding ascorbate-glutathione biosynthesis. Front Plant Sci. 2015;6:458.
Müller P, Li X-P, Niyogi KK. Non-photochemical quenching. A response to excess light energy. Plant Physiol. 2001;125:1558–66.
Baker NR. Chlorophyll fluorescence: a probe of photosynthesis in vivo. Annu Rev Plant Biol. 2008;59:89–113.
Kramer DM, Johnson G, Kiirats O, Edwards GE. New fluorescence parameters for the determination of QA redox state and excitation energy fluxes. Photosynth Res. 2004;79:209–18.
Van Kooten O, Snel JFH. The use of chlorophyll fluorescence nomenclature in plant stress physiology. Photosyn Res. 1990;25:147–50.
Maxwell K, Johnson GN. Chlorophyll fluorescence—a practical guide. J Exp Bot. 2000;51:659–68.
Li X, Liu F. Drought stress memory and drought stress tolerance in plants: biochemical and molecular basis. In: Hossain MA, Wani SH, Bhattacharjee S, Burritt DJ, Tran LSP, editors. Drought stress tolerance in plants, vol. 1. Cham: Springer; 2016. p. 17–44.
Barbagallo RP, Oxborough K, Pallett KE, Baker NR. Rapid, noninvasive screening for perturbations of metabolism and plant growth using chlorophyll fluorescence imaging. Plant Physiol. 2003;132:486–93.
Frachebould Y, Haldimann P, Leipner J, Stamp P. Chlorophyll fluorescence as a selection tool for cold tolerance of photosynthesis in maize (Zea mays L.). J Exp Bot. 1999;50:1533–40.
Lichtenthaler H, Babani F. Detection of photosynthetic activity and water stress by imaging the red chlorophyll fluorescence. Plant Physiol Biochem. 2000;38:889–95.
Zivcak M, Brestic M, Balatova Z, Drevenakova P, Olsovska K, Kalaji HM, Yang X, Allakhverdiev SI. Photosynthetic electron transport and specific photoprotective responses in wheat leaves under drought stress. Photosynth Res. 2013;117:529–46.
Woo NS, Badger MR, Pogson BJ. A rapid, non-invasive procedure for quantitative assessment of drought survival using chlorophyll fluorescence. Plant Methods. 2008;4:1–14.
Jansen M, Gilmer F, Biskup B, Nagel KA, Rascher U, Fischbach A, Briem S, Dreissen G, Tittmann S, Braun S, Jaeger ID, Metzlaff M, Schurr U, Scharr H, Walter A. Simultaneous phenotyping of leaf growth and chlorophyll fluorescence via GROWSCREEN FLUORO allows detection of stress tolerance in Arabidopsis thaliana and other rosette plants. Funct Plant Biol. 2009;36:902–14.
Mishra A, Heyer AG, Mishra HB. Chlorophyll fluorescence emission can screen cold tolerance of cold acclimated Arabidopsis thaliana accessions. Plant Methods. 2014;38:2–10.
Klukas C, Chen D, Pape J-M. IAP: an open-source information system for high-throughput plant phenotyping. Plant Physiol. 2014;165:506–18.
Fahlgren N, Feldman M, Gehan MA, Wilson MS, Shyu C, Bryant DW, Hill ST, McEntee CJ, Warnasooriya SN, Kumar I, Ficor T, Turnipseed S, Gilbert KB, Brutnell TP, Carrington JC, Mockler TC, Baxter I. A Versatile phenotyping system and analytics platform reveals diverse temporal responses to water availability in Setaria. Mol Plant. 2015;8:1520–35.
Flood PJ, Kruijer W, Schnabel SK, van der Schoor R, Jalink H, Snel JFH, Harbinson J, Aarts MGM. Phenomics for photosynthesis, growth and reflectance in Arabidopsis thaliana reveals circadian and long-term fluctuations in heritability. Plant Methods. 2016;12:1–14.
van Rooijen R, Aarts MGM, Harbinson J. Natural genetic variation for acclimation of photosynthetic light use efficiency to growth irradiance in Arabidopsis. Plant Physiol. 2015;167:1412–29.
Cruz JA, Savage LJ, Zegarac R, Kovac WK, Chen J, Kramer DM. Dynamic environmental photosynthetic imaging reveals emergent phenotypes. Cell Syst. 2016;2:365–77.
White AJ, Critchley C. Rapid light curves: a new fluorescence method to assess the state of the photosynthetic apparatus. Photosyn Res. 1999;59:63–72.
Ralph PJ, Gademann R. Rapid light curves: a powerful tool to assess photosynthetic activity. Aquat Bot. 2005;82:222–37.
Saccardy K, Pineau B, Roche O, Cornic G. Photochemical efficiency of photosystem II and xanthophyll cycle components in Zea mays leaves exposed to water stress and high light. Photosyn Res. 1998;56:57–66.
Genty B, Briantais J-M, Baker NR. The relationship between the quantum yield of photosynthetic electron transport and quenching of chlorophyll fluorescence. Biochim Biophys Acta. 1989;990:87–92.
Furbank RT, Badger MR. Photosynthetic oxygen exchange in attached leaves of C4 monocotyledons. Aust J Plant Physiol. 1982;9:553–8.
Fryer MJ, Andrews JR, Oxborough K, Blowers DA, Baker NR. Relationship between CO2 assimilation, photosynthetic electron transport, and active O2 metabolism in leaves of maize in the field during periods of low temperature. Plant Physiol. 1998;116:571–80.
Benesová M, Holá D, Fischer L, Jedelsky PL, Hnilicka F, Wilhelmová N, Rathová O, Kocová M, Procházková D, Honnerová J, Fridrichová L, Hnilicková H. The physiology and proteomics of drought tolerance in maize: early stomatal closure as a cause of lower tolerance to short-term dehydration? PLoS ONE. 2012;7:e38017.
Lavinski AO, Magalhães PC, Ávila R, Gomes-Jr PC, Carneiro NP. Analysis of maize photosynthesis parameters and whole plant oxidative damage under long-term drought. Adv Crop Sci Technol. 2015;S1:007.
Lichtenthaler HK, Langsdorf G, Lenk S, Buschmann C. Chlorophyll fluorescence imaging of photosynthetic activity with the flash-lamp fluorescence imaging system. Photosynthetica. 2005;2005(43):355–69.
Lichtenthaler HK, Ac A, Marek MV, Kalina J, Urban O. Differences in pigment composition, photosynthetic rates and chlorophyll fluorescence images of sun and shade leaves of four tree species. Plant Physiol Biochem. 2007;45:577–88.
CAF Sousa, PC Magalhães and RL Gomide conceived and designed the experiments; CAF Sousa, DS Paiva, RACN Casari and NG Oliveira performed the experiments and collected the data; NG Oliveira performed all the statistical analyses; CAF Sousa, AK Kobayashi, HBC Molinari and MT Sousa Junior analyzed the data and wrote the paper. All authors read and approved the final manuscript.
The authors would like to thank Maize Breeding Program from National Maize and Sorghum Research Center, Brazil, especially to Dr. Lauro José Moreira Guimarães, for supplying the seeds of maize genotypes. N. G. de Oliveira would like to thank the Coordination for the Improvement of Higher Education Personnel (CAPES), a public foundation subordinated to the Ministry of Education of Brazil, for the fellowship grant.
CAF Sousa received his DSc degree in plant physiology from University of Campinas in 2001. He is currently responsible for the research area for plant phenotyping by imaging at National Agroenergy Research Center from Embrapa, Brazil.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated or analyzed during this study are included in this published article [and its supplementary information files].
This research was supported by the Agency of Funding for Studies and Projects linked to the Brazilian Ministry of Science, Technology and Innovation (Grant MCTI/FINEP Number 0111008900) and also by the Brazilian Agricultural Research Corporation (Embrapa/Macroprograma5/Project PHENOCORN: 05.12.12.001.00.00).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.