Quantifying spatial heterogeneity of chlorophyll fluorescence during plant growth and in response to water stress
- Justine Bresson†1, 2, 3,
- François Vasseur†4,
- Myriam Dauzat1,
- Garance Koch1,
- Christine Granier1 and
- Denis Vile†1Email author
© Bresson et al.; licensee BioMed Central. 2015
Received: 11 December 2014
Accepted: 12 March 2015
Published: 26 March 2015
Effects of abiotic and biotic stresses on plant photosynthetic performance lead to fitness and yield decrease. The maximum quantum efficiency of photosystem II (F v/F m) is a parameter of chlorophyll fluorescence (ChlF) classically used to track changes in photosynthetic performance. Despite recent technical and methodological advances in ChlF imaging, the spatio-temporal heterogeneity of F v/F m still awaits for standardized and accurate quantification.
We developed a method to quantify the dynamics of spatial heterogeneity of photosynthetic efficiency through the distribution-based analysis of F v/F m values. The method was applied to Arabidopsis thaliana grown under well-watered and severe water deficit (survival rate of 40%). First, whole-plant F v/F m shifted from unimodal to bimodal distributions during plant development despite a constant mean F v/F m under well-watered conditions. The establishment of a bimodal distribution of F v/F m reflects the occurrence of two types of leaf regions with contrasted photosynthetic efficiency. The distance between the two modes (called S) quantified the whole-plant photosynthetic heterogeneity. The weighted contribution of the most efficient/healthiest leaf regions to whole-plant performance (called W max) quantified the spatial efficiency of a photosynthetically heterogeneous plant. Plant survival to water deficit was associated to high S values, as well as with strong and fast recovery of W max following soil rewatering. Hence, during stress surviving plants had higher, but more efficient photosynthetic heterogeneity compared to perishing plants. Importantly, S allowed the discrimination between surviving and perishing plants four days earlier than the mean F v/F m. A sensitivity analysis from simulated dynamics of F v/F m showed that parameters indicative of plant tolerance and/or stress intensity caused identifiable changes in S and W max. Finally, an independent comparison of six Arabidopsis accessions grown under well-watered conditions indicated that S and W max are related to the genetic variability of growth.
The distribution-based analysis of ChlF provides an efficient tool for quantifying photosynthetic heterogeneity and performance. S and W max are good indicators to estimate plant survival under water stress. Our results suggest that the dynamics of photosynthetic heterogeneity are key components of plant growth and tolerance to stress.
KeywordsArabidopsis thaliana Chlorophyll fluorescence imaging Heterogeneity of F v/F m values Modelling Photosynthetic performance Pixels distribution Plant growth Plant survival Sensitivity analysis
High-throughput phenotyping is increasingly used for dissecting the genetic and eco-physiological determinisms of plant performance and stress tolerance. Over the last decade, efficient automated imaging systems have been developed for the acquisition of visible, bioluminescence, fluorescence and multi-spectral images. A rising difficulty is now to extract valuable, i.e., biologically meaningful, preferably quantitative, information from the large collection of images generated by these systems .
Chlorophyll fluorescence (ChlF) imaging has become one of the most powerful and popular tools to track changes in the photosynthetic capacities of plants in response to abiotic and biotic factors [2-4]. Pulse-amplitude modulated ChlF techniques provide non-invasive assessment of the photosystem II (PSII) efficiency to supply electrons to the photosynthetic machinery. Light energy absorbed by chlorophyll molecules can undergo one of three competing fates: (i) driving photosynthesis (photochemistry); (ii) being dissipated as heat; or (iii) being re-emitted as ChlF. These three processes take place in a competitive manner, and under stress conditions, the photochemistry declines whereas heat dissipation and ChlF emission increase (for recent reviews, see [5,6]). ChlF is estimated by the quantification of the light re-emitted (in the red wavebands) after the application of a saturating flash (usually for a few seconds) to the photosynthetic organs . The saturating flash induces the transport of electrons through PSII centres, driving the reduction of QA, the primary stable electron acceptor of PSII. Once reduced, QA cannot accept new electrons before electrons are transferred to the next acceptor (the reaction centre is considered to be ‘closed’), and the excess of energy is dissipated through heat and fluorescence.
Amongst the different ChlF parameters, the F v/F m ratio is a useful and rapid parameter that reflects the maximum quantum efficiency of the PSII photochemistry . In dark-adapted leaves (in which all PSII reaction centres are in the ‘open’ state; QA fully oxidized), a measuring beam is applied to elicit the minimal value of ChlF, F 0 (i.e., basal fluorescence). F 0 represents the energy dissipation via light-harvesting antenna pigments when excitation energy is not being transferred to the PSII reaction centres. After reaching F 0, the application of a brief saturating pulse induces a maximum value of ChlF, F m (PSII reaction centres get ‘closed’ because of electron accumulation; QA fully reduced). The difference between F 0 and F m is the variable fluorescence, F v and F v/F m is given by (F m-F 0)/F m (for more details, see ). Low F v/F m indicate substantial photoinhibition or down-regulation of PSII that occurs when plants experience stress. It has been shown that F v/F m is a robust indicator of plant health. Healthy photosynthetic tissues of most plant species exhibit a mean F v/F m at ca. 0.83, while lower values are indicative of an impaired physiological status [8,9]. Rapid modifications of F v/F m are for instance reported in response to many environmental factors, such as water stress [8,10], temperature [11-13], wounding , photoinhibition [11,15], biotic interactions such as pathogenic as well as beneficial bacteria [16-19].
Soil water availability is one of the most important environmental factors for plant growth and development. The impact of water deficit on the photosynthetic performance of plants depends on the severity and duration of the stress. In the short-term, decrease in water supply usually induces stomata closure to maintain a favourable leaf water status, what in turn leads to a reduction of internal CO2 concentration . Hence, stomata closure under water stress promotes an imbalance between the photochemical activity of PSII and the electron requirement for carbon fixation, leading to over-excitations and subsequent photoinhibitory damages to PSII reaction centres . As a consequence, substantial decline in F v/F m in response to moderate water deficit is observed in various plant species (see references in ), and was closely related to decreased relative leaf water content . With increasing stress severity or duration, carbon starvation and hydraulic failure, which strongly alter F v/F m at the whole-plant level, lead to partial (or total) senescence or leaf abscission . Even though exacerbated leaf senescence can be lethal, sacrificing a few leaves might be a good strategy to ensure survival under severe resource limitation . Growth recovery following severe water stress is then associated with the (partial) re-establishment of the photosynthetic capacities of the senescing leaves, and/or with the development of new leaves with optimal photosynthetic performance .
ChlF imaging has revealed that photosynthetic performance is extremely heterogeneous at the leaf surface, as well as between leaves, when plants experience environmental stresses. For examples heterogeneity in ChlF is reported in response to changing CO2 concentration , light stimuli , ozone-induced perturbations , low growth temperature , chilling , pathogen attack , drought [10,30] or treatment with abscisic acid . Spatio-temporal heterogeneity across photosynthetic areas has been assessed by visual inspection of leaves [24,26,30], by measurements at spatially different small areas on the leaf surface [10,29,30], or by visual inspection of the shape of F v/F m distributions across leaves [26,28,29,32]. ChlF imaging of leaves of Arabidopsis grown under water stress for instance reveals a progressive decline of F v/F m beginning at the leaf’s tip . However we still lack an automatic and standardized method for the quantification of the spatial heterogeneity of F v/F m values, which is crucial to compare photosynthetic performance depending on the developmental stage, the genotype, or the environmental conditions.
Here, we analysed the distribution of F v/F m to estimate the spatial heterogeneity of photosynthetic efficiency (S) and the weighted contribution of the most efficient/healthiest leaf regions to whole-plant photosynthetic performance (W max). We first showed that the changes in S and W max were related to the survival of the Arabidopsis Col-0 accession to a severe water deficit (SWD). A sensitivity analysis of S and W max to simulated dynamics of F v/F m distributions showed to what extent S and W max can vary depending on plant tolerance and/or stress intensity. Finally, we found that a significant part of the variation in biomass accumulation in six contrasted Arabidopsis accessions is explained by the variation of W max in the course of plant development.
Analytical framework: severe water deficit strongly affects plant growth, photosynthetic efficiency and induces plant mortality
Computing and quantifying the heterogeneity of plant photosynthetic efficiency during growth and under severe water deficit
During SWD, F v/F m values at the plant surface became heterogeneous, as illustrated by the changes in the mean and distribution of F v/F m values (Figures 2A, B). We notably observed the establishment of multimodal distributions during SWD, reflecting the spatial variability of F v/F m in the rosette (Figure 2B). To explore the heterogeneity of F v/F m values during time course, we applied the Hartigan’s non-parametric significance test for unimodality [35-37]. As expected, the proportion of stressed plants showing multimodal distributions increased strongly after stress exposure (Figure 2C). Under WW conditions, the proportion of plants that displayed significant multimodal distributions also increased from < 10% to > 90% between 1 to 10 days after L4 stage. Stressed plants even displayed a slightly lower proportion of multimodal distributions compared to plants grown under WW conditions (Figure 2C).
For plants grown under WW conditions, each parameter was roughly constant during plant development (Figures 3B-G). In stressed plants, while μmax essentially followed the same variation of whole-rosette mean F v/F m (Figures 2A and 3C), μmin decreased to reach very low values (μmin = 0.24 ± 0.13 and 0.37 ± 0.17 for perishing and surviving plants, respectively; Figure 3B). Standard deviation σmax progressively increased during SWD establishment. However, while σmax of surviving plants recovered values similar to WW plants after rewatering, σmax continued to increase in perishing plants (Figure 3E). By contrast, standard deviation σmin increased more in surviving than in perishing plants, but recovered their initial value 13 days after rewatering (Figure 3D). In addition, the weight, i.e. the proportion, of the minimum mode ρmin increased to a greater extent in perishing plants (and the weight of the maximum mode ρmax decreased likewise) compared to surviving plants (Figures 3F, G).
Quantifying the effect of photosynthetic heterogeneity on whole-plant performance: description
Under SWD, S accurately represented the photosynthetic heterogeneity and allowed deciphering surviving and perishing plants. However, it failed to quantify the effect of photosynthetic heterogeneity on plant performance and stress tolerance, as shown by the overlap of S values between WW and perishing plants (Figure 4). This is because the deviation of both modes to the photosynthetic optimum is as important as the disparity between the two modes.
Quantifying the effect of photosynthetic heterogeneity on whole-plant performance: applications
In plants grown in WW conditions, W max increased progressively during development from 0 to ca. 0.85 (Figure 5B). This reflects the increase in the heterogeneity of whole-plant photosynthetic performance (i.e., a switch from unimodality to bimodality) with a very low and negligible effect of the minimum mode compared to the maximum mode. In stressed plants, the increase of W max was delayed and reduced (Figure 5B). In surviving plants, W max started to decrease at 15 days after L4 stage, and recovered shortly (2 days) after rewatering. At bolting, surviving plants exhibited a W max of ca. 0.65, i.e. 23% less than WW plants at the same developmental stage (Figure 5B). By contrast, in perishing plants, W max started to decrease after 14 days following the L4 stage and became negative ten days later.
We used simple mathematical functions to model the dynamics of the parameters of the bimodal distributions in various stressing conditions, and simulate the associated variations of S and W max (see Additional file 3). First, this simulation exercise was sufficient to reproduce what has been observed in this paper in plants grown under SWD that did not survive the stress. The parameters of these functions were then varied to simulate different scenarios of photosynthetic heterogeneity generated by different stress intensities. Our sensitivity analysis of W max showed that it becomes as negative as (1) the rate of decrease in means and (2) the increase in proportion of damaged leaf regions, are high (i.e., low stress tolerance, and/or diffuse stress effect, high stress intensity). Conversely, its decrease is delayed when the rates of decrease in means and rates of changes are low (i.e., high stress tolerance, stress effects with high patchiness, and/or low stress intensity; see Additional file 3).
Second, we investigated the effect of soil inoculation with Phyllobacterium brassicacearum STM196 strain, a plant growth-promoting rhizobacteria (PGPR) that improves plant tolerance to moderate water deficit  and also increases plant survival under SWD . Bresson et al., 2014  showed that STM196-inoculated surviving plants also exhibit a higher growth recovery after rewatering, leading to a higher plant biomass than non-inoculated plants . Here, we showed that STM196-inoculation induced a faster and higher increase in W max than non-inoculated plants from 2 days after rewatering (Additional file 5: Figure S4). This suggests that the positive effects of STM196 on growth recovery, biomass production and plant survival may be related to its effects on whole-plant photosynthetic heterogeneity.
Analysing the effects of environmental conditions on plant growth, survival and yield requires massive, rapid and non-invasive tools to track changes in plant performance. Non-invasive ChlF imaging has been developed to give insights into plant photosynthetic capacities and explore the ability of plants to tolerate various environmental stresses (e.g., [8,16,43]). Most often the mean values of various indices of ChlF, including the widely used F v/F m, of an organ or a plant is used to characterize the response to a stressor (e.g., [8,11]). However, a ChlF image is composed of a panel of pixels in a given range (F v/F m = [0; 1]). Hence, using mean values does not give a clear clue of the disparity of values that corresponds to contrasted physiology. Heterogeneity in the photosynthetic capacities of plants has been observed but rarely quantified in responses to a wide variety of external stimuli (e.g., [10,16,25-32]). For instance, the establishment of F v/F m heterogeneity in response to stress has been described by sampling small areas on the leaf surface [10,29,30], visual inspection of the shape of distributions [25-29,32], or by F v/F m clustering . However, this is prone to large variations depending on the species, experimenter and stress. ChlF heterogeneity is often admitted but its standardized, objective and reproducible quantification is still lacking. For instance, previous methods used threshold-based analysis to quantify the area and progression of senescence or damages [8,16]. Here, we proposed a method to quantify (1) the heterogeneity of F v/F m values at whole-plant level (S) and (2) the spatial efficiency of a photosynthetically heterogeneous plant (W max). Although we applied our method to measurements of F v/F m in the Arabidopsis rosette under a severe water deficit scenario, we argue that the approach can be used with other ChlF parameters (e.g., ΦII, NPQ) as well as in response to other stressing conditions that induce variations of the physiological status.
Photosynthetic heterogeneity is intrinsic to the development of plants
Our results showed that the distribution of F v/F m values shifted from unimodal to bimodal distributions both under WW and SWD conditions, and this despite a constant mean F v/F m in WW plants. This result indicates that heterogeneity in photosynthetic efficiency (i.e., the increase in the proportion of bimodal F v/F m distributions) does not appear only under stress but is intrinsic to the development of plants. Importantly, S and W max in WW plants also significantly increased during development. It therefore indicates that, even in the absence of visible senescence, (1) there were leaf regions exhibiting lower F v/F m, (2) low-efficiency leaf regions increased during development and, (3) the contribution of these latter was minor on whole-plant photosynthetic performance under WW conditions. There might be different sources of photosynthetic heterogeneity. First, at the whole-plant level, photosynthetic heterogeneity in plants might be caused by age-induced leaf senescence, i.e. by visible and non-visible cell death and nutrient remobilization, notably on the edges of the oldest leaves. In addition, the increase in the size of leaf veins with increasing leaf size can also induce a decrease in the mean F v/F m, as well as an increase in F v/F m heterogeneity. Second, at the sub-cellular level, some of PSII centres are inactive to linear electron transport. Functional PSII heterogeneity is for instance expected since 70-80% of PSII are located in the stacked grana region and the remaining PSII are located in the stroma-exposed region of the thylakoid membrane [44-47].
We also showed that the variation in photosynthetic heterogeneity might be a key trait related to plant growth, as suggested by the significant correlation between the increase in W max during development and biomass in six contrasted Arabidopsis accessions and despite no distinct differences in the mean F v/F m between genotypes. The analysis of the distributions of F v/F m values, as proposed with S and W max, allows the quantification of the whole-plant heterogeneity and may be more informative than the whole-plant mean value to investigate changes during plant development and genetic variation in plant performance.
The indicators of photosynthetic heterogeneity (S and W max) are linked to plant tolerance to severe water deficit
Our analysis revealed that SWD affected the establishment of the intrinsic heterogeneity in plants during development. The heterogeneity of F v/F m values (quantified by S) across the rosette increased differently depending on the state of the plants. Importantly, S was a more sensitive indicator of the plant physiological status than the mean F v/F m. Indeed the mean F v/F m was stable during the first 14 days in stressed plants, while a strong photosynthetic heterogeneity was already present (Figures 2A and 4). S allows the discrimination between surviving and perishing plants earlier, ca. four days, than the whole-rosette mean F v/F m. The lag time before recovery was also shorter in S values than the mean F v/F m.
Surprisingly, surviving plants displayed a higher increase of S than the others plants during stress establishment, and perishing plants exhibited S dynamics similar to plants grown under WW conditions. This did not reflect the lower absolute values of F v/F m in perishing plants. The higher photosynthetic heterogeneity in surviving plants can be explained by the establishment of a gradient of F v/F m values from the tip to the base in the oldest leaves, often observed under water stress  and with high F v/F m values in the youngest leaves (as suggested in this study, see Figure 2). The analysis of the different parameters of bimodal distributions shows that SWD did not induce a global decrease of F v/F m, but plants rather maintained leaf regions with near-optimum F v/F m and sacrificed other leaf regions. Moreover, plant survival to SWD was associated to a large variability in F v/F m of the most damaged/senescing leaves; but to a low variability in the healthiest leaves or leaf regions (Figure 3). After rewatering, we showed that surviving plants recovered optimal mean F v/F m values with decreasing S, by loss of senescing leaves and/or by development of new leaves with optimal F v/F m. This is in accordance with the survival strategy of plants aiming at recycling and reallocating resources from the oldest or senescing leaves to active growing organs . On the contrary, perishing plants displayed a decrease of F v/F m values in their oldest but also youngest leaves, resulting in a lower and constant value of bimodal separation S across the rosette. Perishing plants thus presented the same heterogeneity than well-watered plants.
However, the contribution of the most efficient leaf regions was more negatively impacted in the perishing plants. In surviving plants, W max decreased during stress establishment and recovered rapidly after rewatering. This indicated that the contribution of the healthiest leaf regions was sufficient to compensate the negative effects of senescence and to support plant survival. By contrast, in perishing plants, W max constantly decreased and reached negative values, reflecting that the contribution of low-efficiency leaf regions on whole-plant photosynthetic performance was too important and led to plant death. We also showed that increase in W max is related to the improvement of plant survival and biomass production upon inoculation with a PGPR that affects photosynthesis in plants . We argue that W max could be a good indicator to predict plant survival under water stress, and presumably in response to other stresses.
The results of the sensitivity analysis (Additional file 3) showed that W max is specifically sensitive to changes in the proportion of damaged leaf regions and in the lag time, i.e. the time before the appearance of the first symptoms. By contrast, S is independent of the proportion of leaf damaged but strongly impacted by the dynamics of the means. This simulation analysis suggests that the variation of W max during stress (specifically, the rate of decrease and the time to reach negative values) can be a good indicator of stress tolerance in response to many environmental stresses. We argue that the parameters of the mathematical functions used to model the dynamics of mixture distributions could be used as useful tolerance/sensitivity indices. Additional work is however needed to implement physiological hypotheses under the spatio-temporal dynamics of the mixture parameters.
Together, our results showed that the dynamics of S and W max allow quantifying photosynthetic heterogeneity and its relationship with plant performance, during plant development and under stress. Although the variation in mean F v/F m is informative about overall changes in plant performance, we showed that the modifications of F v/F m distributions are not accurately tractable by the modifications of mean F v/F m, or other threshold-based methods. For example, in comparison to Woo et al.  who used a 33% threshold for a mean F v/F m = 0.80 to discriminate surviving to perishing plants, here we showed that the bimodal separation S allowed discriminating plants four days earlier than mean F v/F m. Moreover, the quantitative measurement of photosynthetic heterogeneity can be followed, repeated and modelled in the course of plant development. For instance, negative values of W max indicated death for individual plants; this may represent a powerful tool to screening plants to water stress. Our study highlights that the management of the spatial photosynthetic heterogeneity may be key to plant survival, and that the F v/F m heterogeneity is a sensitive measure of plant responses to the environment. Further studies will undoubtedly improve our capacity to predict plant tolerance, including survival, to different stressing environmental conditions using the indices of photosynthetic heterogeneity.
Heterogeneity in photosynthetic performance has implications for overall plant performance. In this study, we characterized the heterogeneity of chlorophyll fluorescence in leaves impacted by severe water deficit. We used a quantitative measure of (1) the heterogeneity of photosynthetic efficiency S, and (2) the spatial efficiency of a photosynthetically heterogeneous plant, W max. S and W max gave a more accurate indication of the dynamics of leaf senescence or damages induced by water deficit than the whole-plant mean F v/F m. For instance, they allowed an earlier discrimination between surviving and perishing plants. S and W max were also strongly related to the genetic variability of growth between six contrasted accessions of Arabidopsis. Together our analyses suggest that S and W max could be useful indicators of plant responses to their abiotic and biotic environments. Other studies are now needed to explore further the physiological causes and implications of the quantitative variations of photosynthetic heterogeneity and then improve our ability to predict plant response to the environment.
Phenotyping platform PHENOPSIS
The PHENOPSIS platform  includes three automatons set up in growth chambers strictly controlled for environmental conditions. In each chamber, an automated arm is equipped with a balance (XB620C; Precisa, Dietikon, Zurich, Swiss) and a complete system of irrigation (A1645-6; Electro finish, Saint-Egrève, France) to accurately weigh and irrigate up to 504 Arabidopsis thaliana plants in individual pots. The arm is also equipped with multiple devices for non-destructive acquisition of plant phenotypic data such as Charge Coupled Device (CCD) cameras or fluorescence sensors . Light is provided by a bank of cool-white fluorescent tubes (36 W T8 Grolux, 120 cm; Sylvania) and quartz metal halide lamps (HPI-T Plus 400 W 645 E40; Philips). Light intensity is measured continuously at plant height, using a light sensor over the waveband of 400–700 nm (SKP215; Campbell Scientific, Logan, USA). Air temperature and relative humidity are measured every 20 s (HMP45C-L; Campbell Scientific). All measurements of temperature, light intensity and relative humidity are averaged and stored every 600 s in a data-logger (CR10X; Campbell Scientific) with data-logger support software (Loggernet V4; Campbell Scientific). The climatic regulation of the growth-chambers is controlled by Loggernet software allowing the control of the desired environment by employing an air drier or a water sprayer to modify air humidity, an air-cooler or a heater to modify air temperature.
Plant material, growth conditions and irrigation treatments
The experimentations performed in the PHENOPSIS automaton used A. thaliana (L.) Heynh, accession Columbia-0. Five seeds were sown at the soil surface in 250 cm3 cylindrical pots (10 cm high, 6 cm diameter) filled with a damped mixture (1:1, v:v) of loamy soil and organic compost Neuhaus N2 (see Additional file 6: Table S1 for soil chemical properties). Initial soil water content was controlled during pot filling by determining soil fresh weight (FWsoil) and soil dry weight (DWsoil, after 5 days at 80°C) every ten pots. Soil relative water content was calculated as RWCsoil = (FWsoil – DWsoil) × 100 × DWsoil −1. Subsequent changes in pot weight were attributed to a change in soil water status. The pots were kept in the dark for 2 days and were damped with sprayed deionised water three times a day until germination. Then, plants were cultivated under 12 h day length (180 μmol m−2 s−1 photosynthetic photon flux density, at plant height). During germination phase (7 days), air temperature was set to 20°C day and night, and air relative humidity was adjusted in order to maintain constant water vapour pressure deficit (VPD) at 0.6 kPa. Then, plants were grown at 20/17°C day/night and 0.8 kPa of VPD. Just before the beginning of water stress, seedlings of similar sizes and developmental stages were selected and were thinned to one to four plants per pot. Each pot was daily weighed and watered with a modified one-tenth-strength Hoagland solution  to reach the target RWCsoil. RWCsoil was maintained at 0.35 g H2O g−1 dry soil in the WW treatment (35%). SWD was started at L4 stage by stopping irrigation to decrease progressively RWCsoil to reach 0.06 g H2O g−1 dry soil (6%). After RWCsoil = 6% g H2O g−1 dry soil, irrigation was resumed by adding a daily constant volume of nutritive solution to reach the WW soil condition level, and was then maintained until final harvests at first flower open (stage 6.00; ). Soil water potential was determined by using a potentiometer (WP4-T dewpoint meter, Decagon Devices, Pullman, WA 99163, USA) during the soil drying.
An independent experiment was performed in the Max Planck Institute for Developmental Biology (Weigel lab, Tübingen, Germany) on six natural accessions from contrasted geographic origins: ICE107 (South Italia), ICE111 (South Italia), ICE228 (South Tyrol), ICE50 (Spain), Sha (Kazakhstan), Yeg-1 (Caucasus). Each accession was grown in four replicates. Five to ten seeds were sown at the soil surface of each pot and stratified during 2 days in the dark at 4°C. Plants were then grown at 16°C and under 8 h day length. At L4 stage, only one plant per pot was kept and grown until 48 DAS. Pots were manually watered 3 times a week to maintain good soil moisture. F v/F m was measured every 2–3 days from 17 to 48 DAS (using the same ChlF imaging system as in the PHENOPSIS automaton described below). At 48 DAS, rosettes were harvested, dried at 65°C for 4 days and weighed.
High resolution of chlorophyll fluorescence imaging
Acquisition of chlorophyll fluorescence images
ChlF measurements were performed using Imaging-PAM chlorophyll fluorometer and ImagingWin software application (ver. 2-45d, Heinz Walz GmbH) connected with PHENOPSIS automaton. ImagingWin software is driven by Optima PLC (ver. 2–14, build v323, by Optimalog SARL; Saint-Cyr-sur-Loire, France) that allows the automatic movement of Imaging-PAM implemented on the arm of robot (Figure 1A) and the ChlF measurement of each pot with identical settings. The measuring system consists of a 2/3″ Gigabit Ethernet IMAG K6-CCD camera (Manta, G-145B, ASG Allied Vision Technologies GmbH, Stadtroda, Germany), 1392 × 1040 pixel primary resolution with enhanced sensitivity by 4-pixel-binning resulting in 640 × 480 pixel images, coupled to an objective lens (F1.4/f = 12.5 mm; Cosmicar-Pentax, Hamburg, Germany) with a detector filter (RG665, 3 mm) and a short-pass interference filter (λ < 770 nm). Imaging-PAM is equipped with a powerful array of 44 high-power Luxeon LEDs for fluorescence excitation and actinic illumination with blue light (450 nm) as well as assessment of absorbed photosynthetically active radiation with the help of red light (650 nm) and near-infrared (NIR)-light (780 nm). The ChlF imaging system was equipped by a conic, black, metal shading hood of 21.5 cm height wherein the pot was lift up by the balance controlled by a hydraulic cylinder (Figure 1A). This system allows avoiding the illumination of neighbouring plants and achieving the optimal working of 18.5 cm distance from the LEDs resulting in an imaged area of approximately 9 x 12 cm with +/− 7% maximal deviation of intensity from the mean value. The image acquisition takes 30 s per plant, and different files are generated such as .PIM files (which are only usable by ImagingWin software) and .CSV files, which contain requested parameters such as F 0, F m and F v/F m averaged on a region of interest. Two images of ChlF acquisition are also generated: *.JPEG files which give a representation of the image acquired in false colour (from black pixel values (0) through red, yellow, green, blue to purple (ending at 1)) and *.TIFF files with contain different stacks of fluorescence parameters (see below).
Measurement of maximum efficiency of PSII (F v/F m)
The ChlF measurement was initiated by exposing dark-adapted leaf to measuring light pulses (1 Hz frequency, Intensity 2) for determination of F 0. The F m level of fluorescence is recorded during a saturating pulse (Si 9, width 800 ms). F v/F m, calculated as (F m-F 0)/F m, provides the maximum quantum yield of PSII (i.e., photosynthetic efficiency). In the study, F v/F m was daily performed, from the two first leaves to bolting, on dark-adapted plants (8–12 h of dark), under WW and SWD conditions.
Data extraction of photosynthetic efficiency and rosette expansion during time course
We developed an ImageJ macro “PHENOPSIS-Fluo” to semi-automatically extract whole-rosette F v/F m pixel values and total leaf area from the generated TIFF files (containing F 0 and F m pixel values, and NIR images). The analysis of ChlF images starts with image segmentation from NIR pictures, i.e. the automatic separation of the region of interest (here the rosette) from the background. Then, by subtracting and dividing F 0 and F m stacks, [(F m – F 0)/F m], the macro generates an image of F v/F m pixel values, which are represented in the 255 greyscale (0 and 255 corresponding to the minimum and maximum pixel values of the selection, respectively). The macro gives the whole-rosette F v/F m mean and the list of pixel values across the plant. Projected area of the rosette (RA proj ) was also determined from plant selection. Different plants in the same pot can be independently measured. Whole-rosette mean F v/F m extracted with the macro “PHENOPSIS-Fluo” was highly correlated with the mean of the rosette extracted from the commercial software ImagingWin (R 2 = 0.98; Additional file 7: Figure S5). The “PHENOPSIS-Fluo” macro is available on the PHENOPSIS website (http://bioweb.supagro.inra.fr/phenopsis/MacroImageJ.php).
Under SWD conditions, growth dynamics of surviving plants were modelled as a sigmoid curve fitted following RA proj = a / [1 + exp-[(d-a/2)/b]] where a is the maximum area, and d is the number of days after L4 stage. Under WW conditions, an exponential curve was fitted such as RA proj = exp (a – b × d).
Statistical analyses and modelling procedures
All analyses were performed using R 3.1. Comparisons of mean trait values between treatments were performed with Kruskal-Wallis non-parametric tests. Survival percentage was calculated as the proportion of surviving plants at the end of experiment compared to initial number of plants. Plant survival was analysed by Chi2 tests. Non-parametric significance test for unimodality, Hartigan’s dip test (R package ‘dip test’ [35-37]) was used to identity multimodal distribution, with p-value < 0.01 as the significance threshold for departure from unimodality. Then, the REBMIX algorithm for finite mixture models (; R package ‘rebmix’) as used to characterize each mode i of the bimodal distributions of F v/F m values (i.e., mean μ i , standard deviation σ i and weight ρ i ), using the Bayesian Information Criteria (BIC).
We developed a modelling procedure of the temporal dynamics of the means, standard deviations and proportions (weights) of the mixture distributions. We then simulated the variation of these parameters and, subsequently, the variation of S and W max (Additional file 3).
- QA :
Quinone A, the primary stable electron acceptor of PSII centres
- F v/F m :
Maximum quantum efficiency of PSII photochemistry (photosynthetic efficiency)
- F 0 :
Minimal fluorescence emission of a dark-adapted plant
- F m :
Maximum fluorescence emission after a short pulse of a saturating light
- F v :
Variable fluorescence from dark-adapted plant
Spatial heterogeneity of F v/F m
- W max :
Spatial efficiency of a photosynthetically heterogeneous plant
Severe water deficit
- RWCsoil :
Soil relative water content
Well watered condition
Days after stratification
Plant growth-promoting rhizobacteria
Charge Coupled Device
- FWsoil :
Soil fresh weight
- DWsoil :
Soil dry weight
Water vapour pressure deficit
- RA proj :
Projected area of the rosette
We are very grateful to G. Rolland, A. Bediee, F. Bouvery and C. Balsera for help during the experiments. We thank F. Fiorani for giving us advice on the design of the ChlF system. We thank D. Fourré, M. Desigaux, O. Meyerhoff and J. Kolbowski for helpful technical supports. We thank B. Touraine and F. Varoquaux for constructive remarks on plant-microbes interactions and for providing us the rhizobacteria (STM196). We thank B. Genty and S. Bieker for helpful comments and suggestions on this paper. We thank D. Weigel for his contribution and support at different stages of the study. We thank A. Christophe for the helpful discussions during the modelling procedure. We also thank two anonymous reviewers for helpful comments on this paper. JB was funded by French Ministry of Higher Education and Research and was supported by the Institutional Strategy of the University of Tuebingen (Deutsche Forschungsgemeinschaft, ZUK 63). FV was funded by a post-doctoral grant from the Max Planck Society. This work was made possible by a series of grants supporting the development of the plant phenotyping platform: the European Plant Phenotyping Network funded by the FP7 Research Infrastructures Program of the European Union (EPPN, grant agreement no. 284443) and the EIT Climate-KIC project AgWaterBreed.
- Granier C, Vile D. Phenotyping and beyond: modelling the relationships between traits. Curr Opin Plant Biol. 2014;18:96–102.View ArticlePubMedGoogle Scholar
- Baker NR, Rosenqvist E. Applications of chlorophyll fluorescence can improve crop production strategies: an examination of future possibilities. J Exp Bot. 2004;55:1607–21.View ArticlePubMedGoogle Scholar
- Porcar-Castell A, Tyystjärvi E, Atherton J, van der Tol C, Flexas J, Pfündel EE, et al. Linking chlorophyll a fluorescence to photosynthesis for remote sensing applications: mechanisms and challenges. J Exp Bot. 2014;doi:10.1093/jxb/eru191.
- Chaerle L, Leinonen I, Jones HG, Van Der Straeten D. Monitoring and screening plant populations with combined thermal and chlorophyll fluorescence imaging. J Exp Bot. 2007;58:773–84.View ArticlePubMedGoogle Scholar
- Murchie EH, Lawson T. Chlorophyll fluorescence analysis: a guide to good practice and understanding some new applications. J Exp Bot. 2013;doi:10.1093/jxb/ert208.
- Maxwell K, Johnson GN. Chlorophyll fluorescence - a practical guide. J Exp Bot. 2000;51:659–68.View ArticlePubMedGoogle Scholar
- Genty B, Briantais J-M, Baker NR. The relationship between the quantum yield of photosynthetic electron transport and quenching of chlorophyll fluorescence. Biochimica et Biophysica Acta (BBA)-General Subjects. 1989;990:87–92.View ArticleGoogle Scholar
- Woo NS, Badger MR, Pogson BJ. A rapid, non-invasive procedure for quantitative assessment of drought survival using chlorophyll fluorescence. Plant Methods. 2008;4:27.View ArticlePubMed CentralPubMedGoogle Scholar
- Hunt S. Measurements of photosynthesis and respiration in plants. Physiol Plant. 2003;117:314–25.View ArticlePubMedGoogle Scholar
- Sperdouli I, Moustakas M. Spatio-temporal heterogeneity in Arabidopsis thaliana leaves under drought stress. Plant Biol. 2012;14:118–28.PubMedGoogle Scholar
- Gray GR, Hope BJ, Qin XQ, Taylor BG, Whitehead CL. The characterization of photoinhibition and recovery during cold acclimation in Arabidopsis thaliana using chlorophyll fluorescence imaging. Physiol Plant. 2003;119:365–75.View ArticleGoogle Scholar
- Ehlert B, Hincha DK. Chlorophyll fluorescence imaging accurately quantifies freezing damage and cold acclimation responses in Arabidopsis leaves. Plant Methods. 2008;4:12.View ArticlePubMed CentralPubMedGoogle Scholar
- Janka E, Körner O, Rosenqvist E, Ottosen C-O. High temperature stress monitoring and detection using chlorophyll a fluorescence and infrared thermography in chrysanthemum (Dendranthema grandiflora). Plant Physiol Biochem. 2013;67:87–94.View ArticlePubMedGoogle Scholar
- Quilliam RS, Swarbrick PJ, Scholes JD, Rolfe SA. Imaging photosynthesis in wounded leaves of Arabidopsis thaliana. J Exp Bot. 2006;57:55–69.View ArticlePubMedGoogle Scholar
- Bolhàr-Nordenkampf HR, Öquist G. Chlorophyll fluorescence as a tool in photosynthesis research. In: Hall DO, Scurlock JMO, Bolhàr-Nordenkampf HR, Leegood RC, Long SP, editors. Photosynthesis and Production in a Changing Environment. Springer Netherlands; 1993. p. 193–206.
- Rousseau C, Belin E, Bove E, Rousseau D, Fabre F, Berruyer R, et al. High throughput quantitative phenotyping of plant resistance using chlorophyll fluorescence image analysis. Plant Methods. 2013;9:17.View ArticlePubMed CentralPubMedGoogle Scholar
- Rincon A, Valladares F, Gimeno TE, Pueyo JJ. Water stress responses of two Mediterranean tree species influenced by native soil microorganisms and inoculation with a plant growth promoting rhizobacterium. Tree Physiol. 2008;28:1693–701.View ArticlePubMedGoogle Scholar
- Heidari M, Golpayegani A. Effects of water stress and inoculation with plant growth promoting rhizobacteria (PGPR) on antioxidant status and photosynthetic pigments in basil (Ocimum basilicum L.). J Saudi Soc Agric Sci. 2012;11:57–61.Google Scholar
- Bresson J, Vasseur F, Dauzat M, Labadie M, Varoquaux F, Touraine B, et al. Interact to survive: Phyllobacterium brassicacearum improves Arabidopsis tolerance to severe water deficit and growth recovery. PLoS One. 2014;9:e107607.View ArticlePubMed CentralPubMedGoogle Scholar
- Chaves MM. Effects of water deficits on carbon assimilation. J Exp Bot. 1991;42:1–16.View ArticleGoogle Scholar
- Krause GH. Photoinhibition of photosynthesis. An evaluation of damaging and protective mechanisms. Physiol Plant. 1988;74:566–74.View ArticleGoogle Scholar
- McDowell NG. Mechanisms linking drought, hydraulics, carbon metabolism, and vegetation mortality. Plant Physiol. 2011;155:1051–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Munné-Bosch S, Alegre L. Die and let live: leaf senescence contributes to plant survival under drought stress. Funct Plant Biol. 2004;31:203–16.View ArticleGoogle Scholar
- Hensel LL, Grbić V, Baumgarten DA, Bleecker AB. Developmental and age-related processes that influence the longevity and senescence of photosynthetic tissues in Arabidopsis. Plant Cell. 1993;5:553–64.View ArticlePubMed CentralPubMedGoogle Scholar
- Genty B, Meyer S. Quantitative mapping of leaf photosynthesis using chlorophyll fluorescence imaging. Funct Plant Biol. 1995;22:277–84.Google Scholar
- Baker NR, Oxborough K, Lawson T, Morison JI. High resolution imaging of photosynthetic activities of tissues, cells and chloroplasts in leaves. J Exp Bot. 2001;52:615–21.View ArticlePubMedGoogle Scholar
- Leipner J, Oxborough K, Baker NR. Primary sites of ozone-induced perturbations of photosynthesis in leaves: identification and characterization in Phaseolus vulgaris using high resolution chlorophyll fluorescence imaging. J Exp Bot. 2001;52:1689–96.View ArticlePubMedGoogle Scholar
- Oxborough K, Baker N. An instrument capable of imaging chlorophyll a fluorescence from intact leaves at very low irradiance and at cellular and subcellular levels of organization. Plant Cell Environ. 1997;20:1473–83.View ArticleGoogle Scholar
- Hogewoning SW, Harbinson J. Insights on the development, kinetics, and variation of photoinhibition using chlorophyll fluorescence imaging of a chilled, variegated leaf. J Exp Bot. 2007;58:453–63.View ArticlePubMedGoogle Scholar
- Calatayud A, Roca D, Martínez P. Spatial-temporal variations in rose leaves under water stress conditions studied by chlorophyll fluorescence imaging. Plant Physiol Biochem. 2006;44:564–73.View ArticlePubMedGoogle Scholar
- Daley PF, Raschke K, Ball JT, Berry JA. Topography of photosynthetic activity of leaves obtained from video images of chlorophyll fluorescence. Plant Physiol. 1989;90:1233–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Lazár D, Sušila P, Nauš J. Early detection of plant stress from changes in distributions of chlorophyll a fluorescence parameters measured with fluorescence imaging. J Fluoresc. 2006;16:173–6.View ArticlePubMedGoogle Scholar
- Granier C, Aguirrezabal L, Chenu K, Cookson SJ, Dauzat M, Hamard P, et al. PHENOPSIS, an automated platform for reproducible phenotyping of plant responses to soil water deficit in Arabidopsis thaliana permitted the identification of an accession with low sensitivity to soil water deficit. New Phytol. 2006;169:623–35.View ArticlePubMedGoogle Scholar
- Boyes DC, Zayed AM, Ascenzi R, McCaskill AJ, Hoffman NE, Davis KR, et al. Growth stage-based phenotypic analysis of Arabidopsis: a model for high throughput functional genomics in plants. Plant Cell. 2001;13:1499–510.View ArticlePubMed CentralPubMedGoogle Scholar
- Hartigan JA, Hartigan P. The dip test of unimodality. The Annals of Statistics. 1985;13:70–84.View ArticleGoogle Scholar
- Hartigan PM. Algorithm AS 217: Computation of the Dip Statistic to Test for Unimodality. J R Stat Soc: Ser C: Appl Stat. 1985;34:320–5.Google Scholar
- Maechler M, Ringach D. diptest: Hartigan’s dip test statistic for unimodality–corrected code. R package version 075–74. 2012. Available online at: http://CRANR-projectorg/package=diptest [Retrieved: January 4, 2013].
- Nagode M. rebmix. An R Package for continuous and discrete finite mixture models. 2004.
- Zhang C, Mapes BE, Soden BJ. Bimodality in tropical water vapour. Q J Roy Meteorol Soc. 2003;129:2847–66.View ArticleGoogle Scholar
- Demmig-Adams B, Adams Iii WW, Barker DH, Logan BA, Bowling DR, Verhoeven AS. Using chlorophyll fluorescence to assess the fraction of absorbed light allocated to thermal dissipation of excess excitation. Physiol Plant. 1996;98:253–64.View ArticleGoogle Scholar
- Björkman O, Demmig B. Photon yield of O2 evolution and chlorophyll fluorescence characteristics at 77 K among vascular plants of diverse origins. Planta. 1987;170:489–504.View ArticlePubMedGoogle Scholar
- Bresson J, Varoquaux F, Bontpart T, Touraine B, Vile D. The PGPR strain Phyllobacterium brassicacearum STM196 induces a reproductive delay and physiological changes that result in improved drought tolerance in Arabidopsis. New Phytol. 2013;200:558–69.View ArticlePubMedGoogle Scholar
- 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:485–93.View ArticlePubMed CentralPubMedGoogle Scholar
- Chylla R, Garab G, Whitmarsh J. Evidence for slow turnover in a fraction of photosystem II complexes in thylakoid membranes. Biochimica et Biophysica Acta (BBA)-Bioenergetics. 1987;894:562–71.View ArticleGoogle Scholar
- Chylla RA, Whitmarsh J. Light saturation response of inactive photosystem II reaction centers in spinach. Photosynth Res. 1990;25:39–48.View ArticlePubMedGoogle Scholar
- Cao J. Chlorophyll a fluorescence transient as an indicator of active and inactive photosystem II in thylakoid membranes. Biochimica et Biophysica Acta (BBA)-Bioenergetics. 1990;1015:180–8.View ArticleGoogle Scholar
- Lavergne J, Leci E. Properties of inactive photosystem II centers. Photosynth Res. 1993;35:323–43.View ArticlePubMedGoogle Scholar
- Hoagland DR, Arnon DI. The water-culture method for growing plants without soil. Calif Agric Exper Stn Circ. 1950;347:1–32.Google Scholar
- R Development Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2009.Google Scholar
- Fabre J, Dauzat M, Negre V, Wuyts N, Tireau A, Gennari E, et al. PHENOPSIS DB: an Information System for Arabidopsis thaliana phenotypic data in an environmental context. BMC Plant Biol. 2011;11:77.View ArticlePubMed CentralPubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.