Skip to main content

Chlorophyll fluorescence emission can screen cold tolerance of cold acclimated Arabidopsis thaliana accessions



An easy and non-invasive method for measuring plant cold tolerance is highly valuable to instigate research targeting breeding of cold tolerant crops. Traditional methods are labor intensive, time-consuming and thereby of limited value for large scale screening. Here, we have tested the capacity of chlorophyll a fluorescence (ChlF) imaging based methods for the first time on intact whole plants and employed advanced statistical classifiers and feature selection rules for finding combinations of images able to discriminate cold tolerant and cold sensitive plants.


ChlF emission from intact whole plant rosettes of nine Arabidopsis thaliana accessions was measured for (1) non-acclimated (NAC, six week old plants grown at room temperature), (2) cold acclimated (AC, NAC plants acclimated at 4°C for two weeks), and (3) sub-zero temperature (ST) treated (STT, AC plants treated at -4°C for 8 h in dark) states. Cold acclimation broadened the slow phase of ChlF transients in cold sensitive (Co, C24, Can and Cvi) A. thaliana accessions. Similar broadening in the slow phase of ChlF transients was observed in cold tolerant (Col, Rsch, and Te) plants following ST treatments. ChlF parameters: maximum quantum yield of PSII photochemistry (FV/FM) and fluorescence decrease ratio (RFD) well categorized the cold sensitive and tolerant plants when measured in STT state. We trained a range of statistical classifiers with the sequence of captured ChlF images and selected a high performing quadratic discriminant classifier (QDC) in combination with sequential forward floating selection (SFFS) feature selection methods and found that linear combination of three images showed a reasonable contrast between cold sensitive and tolerant A. thaliana accessions for AC as well as for STT states.


ChlF transients measured for an intact whole plant is important for understanding the impact of cold acclimation on photosynthetic processes. Combinatorial imaging combined with statistical classifiers and feature selection methods worked well for the screening of cold tolerance without exposing plants to sub-zero temperatures. This opens up new possibilities for high-throughput monitoring of whole plants cold tolerance via easy and fully non-invasive means.


Cold tolerance is the ability of plants to withstand low temperatures and plays a crucial role in worldwide production of many important agricultural crops. In temperate climate zones, plants have developed strategies to adjust low temperature tolerance by employing a highly complex process of physiological re-arrangements, termed cold acclimation that is triggered by low but non-freezing temperatures. Numerous studies revealed that cold acclimation not only programs massive changes in transcriptome and metabolome [17] but also induces structural and compositional modifications of compatible solutes in various sub-cellular compartments [8]. Following acclimation, plants are more efficient in dealing with the impact of sustained cold or sudden temperature drops by activating a plethora of adjustments including accumulation of cryoprotective metabolites and proteins for their survival [9]. A correlation of cold acclimation capacities with habitat winter temperatures points to high metabolic or ecological costs of these adjustments [5, 10, 11]. Despite intensive research the molecular mechanisms of cold tolerance are still not fully understood and remain an area for intensive research, as understanding of mechanisms responsible for cold acclimation would allow breeding of cold tolerant crops.

The primary requirement for research associated with engineering of low temperature tolerant crops is to develop efficient and cost effective methods for measuring cold tolerance. Measurement of survival or re-growth scores following triggers of cold acclimation and de-acclimation treatments are quite lengthy and may not be very accurate [12, 13]. Most quantitative methods for assessing cold tolerance employ analysis of damage to plasma membranes by electrolyte leakage (EL) screens or the thylakoids via plastocyanin release after freeze-thaw cycles [14]. Generally, the temperature affording half of maximal damage (LT50), i.e. the temperature at which 50% of electrolytes or plastocyanin are released from their respective compartments, is being evaluated and used as a proxy for cold tolerance. However, since the measurement of purified plastocyanin is comparatively difficult [15], the most popular method to study plant cold tolerance is the evaluation of EL of detached leaves by using conductivity measurements [5, 16]. According to this method, detached leaves are placed in reaction tubes, which are progressively cooled down to certain temperatures. Low-temperature treated samples are then thawed at 4°C and electrical conductivity of a bathing solution is measured to calculate LT50 values that represent half-maximum damage to the plasma membrane [16]. EL is a widely accepted method for quantification of plant cold tolerance, although its results sometimes deviate from those obtained with alternative methods [17, 18].

The response of thylakoids membranes to sustained low-temperature treatments can be measured non-invasively by the analysis of ChlF emission. ChlF is re-emitted by chlorophyll a molecules following light absorption, and it is modulated by photochemical and non - photochemical events in the photosynthetic pigment–protein complexes PSII and PSI of the thylakoid membranes [19, 20]. The relative contribution of PSII and PSI in ChlF emission is also reported to change during cooling treatments [21]; thereby, this method can provide important insights into molecular processes of cold acclimation. The ChlF parameter FV/FM of detached leaves following freeze - thaw cycles has successfully been applied to quantify cold tolerance in Arabidopsis thaliana as well as in other plant species [reviewed in [8, 2224]. The analysis of polyphasic fluorescence rise of from initial low fluorescence (FO) to peak FP by JIP test [25] has also been reported to be useful for the selection of cold tolerance in wheat genotypes [26, 27]. Earlier, we measured the ChlF transients of detached leaves of differentially cold tolerant A. thaliana accessions during progressive cooling starting from room temperature to -15°C and found high correlation between LT50 measured by electrolyte leakage with ChlF parameters such as FO(-15°C)/FO(4°C) and [FS/FO]-15°C[10]. We demonstrated that application of advanced statistics-based combinatorial imaging methods to the sequence of time resolved ChlF images could be used to categorize cold tolerance levels by training of classifiers using fluorescence emission of detached leaves that were slowly cooled at mild sub-zero-temperatures of around -4°C [10]. Vaclavik et al. [7] reported that cold acclimation induced accumulation of Gluconapin and Flavon-3-ol glycosides, respectively, in cold tolerant vs. cold sensitive A. thaliana accessions, and thus metabolomic patterns, could be used for the screening of cold tolerance already in the cold acclimated state i.e. without exposing plants to freezing temperatures. Although metabolomic based methods can thus have better discriminating capacity of cold tolerant and cold sensitive accessions as compared to EL methods, they are not readily applicable in high-throughput screenings. Therefore, in the search for easy methods of sensing cold tolerance, we have tested the potential of ChlF transients measured on whole plant rosettes of the model species Arabidopsis thaliana. We involved NAC, AC, and STT plants of nine A. thaliana accessions that span the north–south range of the species [5]. A short light/dark protocol was applied for recording of ChlF transients of dark-adapted plants for about 202 s using a pulse amplitude modulation (PAM) based fluorometer [28]. We found that ChlF transients averaged over whole intact plant rosettes offered valuable information for assessing impacts of cold acclimation and cold induced alterations of the photosynthetic machinery. In addition to classical analysis of ChlF transients, we employed statistical classifiers and feature selection methods on captured ChlF images to search for highly contrasting features of cold tolerant vs. cold sensitive accessions. We have trained several classifiers and chose the best performing classifier, i.e. a QDC, for further use along with a high performing SFFS feature selection methods to identify features correlated with cold tolerance. The sets of images obtained for the STT state were then tested on AC state data sets, and we found that this method, i.e. combinatorial imaging, can be applied for assessing plant cold tolerance already in the cold acclimated state, i.e. without any sub-zero temperature treatments and by fully non-invasive means.


ChlF transients from whole plants are highly informative for categorization of cold tolerance in A. thaliana accessions

We measured the ChlF transients of whole plant rosettes in nine differentially cold tolerant Arabidopsis thaliana accessions in the NAC, AC, and STT state [Figure 1]. From the shape of ChlF transients it can be established that for the cold sensitive accessions (Co, C24, Can, and Cvi), the fluorescence intensity following peak FP quenched slowly in cold acclimated (AC) state, whereas the value of FP itself decreased significantly when these plants were treated with sub-zero temperature (ST, -4°C for 8 h in dark). In cold tolerant accessions (Te, Rsch, and Col-0), however, the shape of the ChlF transients remained similar for NAC and AC states; while after STT, quenching of fluorescence intensity following peak FP is significantly slowed down. Qualitative comparison of the ChlF transients of intermediate accessions (Ler and Nd) revealed a compound picture, and based on qualitative comparisons we tend to assign Ler as cold sensitive, while Nd can be classified cold tolerant.

Figure 1
figure 1

The chlorophyll fluorescence (ChlF) transients of nine differentially cold tolerant Arabidopsis thaliana accessions measured by short experimental light/dark protocol. The cold sensitive accessions are Co, C24, Can and Cvi (d); intermediate cold tolerant accession are Ler and Nd; and cold tolerant accessions are Col, Rsch and Te. The ChlF transient was integrated across whole plant rosettes and mean of three independent plants for each accession was presented.

ChlF parameters FV/FM and RFD of whole plant rosettes can measure cold tolerance in subzero-temperature (ST) treated plants

Fluorescence parameters FV/FM and RFD were evaluated from measured ChlF emission and presented as an average of three independent experiments, where values of each parameter were estimated by integrating across whole rosette leaves. We did not find correlations between the classical ChlF parameters and cold tolerance among the investigated A. thaliana accessions either for NAC or for the AC state. However, the ChlF parameters FV/FM and RFD, differed for cold sensitive, intermediate and cold tolerant accessions following STT (Figure 2). In STT state, FV/FM values for cold tolerant (Te, Col and Rsch), intermediate (Nd and Ler) and sensitive accessions (Co, C24, and Can) were around 0.85 ± 0.01, 0.83 ± 0.01 and 0.81 ± 0.01, respectively, while the averaged RFD values were calculated as 1.98 ± 0.19, 1.48 ± 0.07 and 1.40 ± 0.19, respectively. Non-paired t-test revealed FV/FM to differ significantly (p <0.01) for cold tolerant and cold sensitive accessions, while there was no significant difference between cold sensitive and intermediate or intermediate and tolerant accessions. In contrast, plant vitality index RFD showed significant differences between the cold sensitive and tolerant as well as between intermediate and tolerant accessions (p <0.05), but no significant difference was found between sensitive and intermediate tolerant accessions. For the accessions Cvi the values of ChlF parameters FV/FM and RFD are almost similar to that of intermediate accessions Ler and Nd, thereby its tolerance level of thylakoid may be categorized intermediate, while the plasma membrane behaved cold sensitive in EL measurements.

Figure 2
figure 2

Chlorophyll fluorescence (ChlF) parameters : maximum quantum yield of PSII photosystems (F V /F M ) [A] and fluorescence decrease ratio (R FD ) [B] of differentially cold tolerant accessions of Arabidopsis thaliana for non-acclimated (NAC), cold acclimated (AC), and sub-zero temperature treated (STT) states. The presented numeric values are mean of three independent plants with standard errors, and are integrated across whole plant rosettes.

Cold acclimation induced gain in photosynthetic performance of A. thaliana accessions

The ChlF parameter effective quantum efficiency of PS II photochemistryPSII) was significantly higher (p <0.01) in AC vs. NAC plants except for Cvi (Figure 3A). Interestingly, mild sub-zero temperature treatments (-4°C) for eight hours in AC plants led to a significant decline in the value of ФPSII in all accessions (Figure 3A) without ФPSII being correlated with the cold tolerance. Two weeks of cold acclimation caused a significant (p <0.01) increase of photochemical quenching (qP) in all accessions, while STT led to a decline (Figure 3B). As for ФPSII, no correlation of parameter qP with cold tolerance was observed.

Figure 3
figure 3

The effective quantum efficiency PS II PSII ) [A] and photochemical quenching (qP) [B] of cold sensitive (Co, C24, Cvi, Can), intermediate cold tolerant (Nd and Ler) and cold tolerant (Te, Rsch and Col) Arabidopsis thaliana accessions for non-acclimated (NAC), cold acclimated (ACC), and sub-zero temperature treated (STT) states. The presented numeric values are integrated across whole plant rosettes and mean of three independent plants with standard errors.

Combinatorial imaging

We applied statistical techniques of classifier and feature selection methods in order to identify features of cold tolerance from captured sequences of time-resolved ChlF images. This method is very powerful in identifying image sets from large sequences of time-resolved ChlF recordings that yield the highest contrast between groups to be compared. Earlier we demonstrated the usefulness of this approach for discriminating three species of family Lamiaceae at very early stages of growth [29] and for investigating features of cold tolerance at non-lethal temperatures [10]. The method depends on training and performance testing of randomly selected pixels of the most contrasting ChlF image sets, i.e. recordings for the highly cold tolerant accession Te and highly sensitive Co in the STT state [for technical details, see [10, 29, 30]. The performance, error rate and computational time to run the algorithms of five tested classifiers, linear discriminant classifier (LDC), quadratic discriminant classifier (QDC), k-nearest neighbors classifier (k-NNC), nearest neighbor classifier (NNC) and nearest mean classifier (NMC), are presented in Table 1. We found that algorithms of QDC are the best performing classifier (81% correct assignment of cold tolerance group among test images) in comparatively short time (<6 s, Table 1). Therefore, QDC was chosen and applied with the SFFS feature selection method to find the most contrasting image sets for highly cold tolerant (Te) vs. highly cold sensitive (Co) accessions. Figure 4 shows the performance curve for the most efficient feature selection method SFFS in combination with simulated classifier QDC, where x-axis represents number of images. Thus, in this experiment, the SFFS algorithm reduced the full data set of 218 images to three images identified as I21, I104 and I107 without compromising the classification performance (~80%, Figure 4). We obtained the linear combination: C = (0.3219)*I21 + (-0.5018)*I104+ (-0.2315)*I107, to perform best in discriminating cold tolerant and sensitive accessions. The coefficients of the linear combination were calculated according to Matouš et al. [30]. Figure 5 shows linear combinations of images for all nine accessions in the STT as well as AC state. The presented images are plotted on a false color scale where color represents the virtual fluorescence intensity of the pixel with reference to its corresponding value in the training data set. Thus, we can visualize a clear difference between cold sensitive, intermediate, and cold tolerant accessions in the STT state. When the combinatorial imaging was applied to AC plants it apparently failed to discriminate cold sensitive and intermediate accessions showing very similar patterns for leaf rosettes of the two groups of plants. However, the tolerant group of accessions could be clearly discriminated already in the AC state.

Table 1 The table presents performance, error rate and computational time of the tested classifiers
Figure 4
figure 4

The performance curve of sequential forward floating selection (SFFS) with images in the x-axis.

Figure 5
figure 5

Illustrates the combination of three most contrasting images for all nine accessions of Arabidopsis thaliana . After subzero temperature treatment the all nine accession has been visibly divided in three categories namely, most sensitive (Co, C24 and Can), intermediates (Ler and Nd) and the most tolerant accessions (Te, Rsch and Col).


Towards easy methods for assessing whole plant cold tolerance in natural accessions of A. thaliana

In an attempt to establish an easy and non-invasive method for screening plant cold tolerance, we compared different concepts for measuring cold tolerance based on chlorophyll fluorescence that have been approved for detached leaves of the model species Arabidopsis thaliana. As far as we know, for the first time we attempted to capture the ChlF transients of intact whole plant rosettes and found interesting modulation in it for the group of cold sensitive (Co, C24, Can, Cvi), intermediate tolerant (Ler, Nd) and cold tolerant (Te, Rsch, Col-0) A. thaliana accessions for NAC, AC and STT state (Figure 1). The ChlF transient measured by quenching analysis of modulated fluorescence is very useful methods that can monitor slight changes in the photochemical and non-photochemical process [31]. The slow quenching of ChlF transients following peak FP, i.e. broadening of ChlF transients, for cold sensitive group of accessions under cold acclimated state and for tolerant group of accessions under STT state (Figure 1) is complex and can be affected by series of events such as light induced intrathylakoid H+ concentration, non-photochemical quenching, inhibition of CO2 assimilation processes, ATP synthesis among others [reviewed in [32]. However, it reveals that sensitive and tolerant group of accessions possess two different strategies to utilize the absorbed irradiance following cold acclimation as well as when placed at mild subzero temperature (-4°C) for 8 h in darkness. Probably cold tolerant plants efficiently utilize acclimation temperature and keeping homeostasis of photosynthetic apparatus and yielded almost similar ChlF transients in NAC as well as in AC state that varied only when treated with sub-zero temperature. The reduction of FP for the group of cold sensitive accessions in STT state indicates perturbation in PSII functional properties. Significant inhibition of FV/FM for sensitive vs. tolerant accessions in STT states further support slight inactivation of PSII function. Ehlert and Hincha [23] reported that the easily recordable parameter FV/FM discriminates cold tolerant and sensitive detached leaves after a freeze-thaw cycle. Using a set of nine differentially cold tolerant natural accessions of Arabidopsis[5], we could verify also for whole plants that FV/FM differs for tolerant and sensitive accessions after a sub-zero temperature treatment, but plants with intermediate tolerance differed neither from sensitive nor from tolerant ones. It has been reported that the fluorescence decrease ratio, RFD, is more sensitive to a variety of stress factors than FV/FM, which mainly responds to extreme conditions [33, 34]. Indeed, RFD, which can be calculated from fluorescence peak FP and FS measured under prevailing actinic light, was not only able to discriminate sensitive and tolerant but also tolerant and intermediate accessions. Interestingly, both FV/FM and RFD data revealed that plastids of the accession Cvi behaved as intermediately tolerant, while the electrolyte leakage method classified this accession as sensitive based on damage to the plasma membrane during freeze-thaw treatments [5]. We have earlier reported that cold tolerance of plastids and the plasma membrane appear not strictly correlated in the accession Col-0 [8], and the differential tolerance of plastids and plasma membrane in Cvi is further evidence for independent principles underlying cold tolerance of the different types of membranes.

Considering plants not treated with sub-zero temperatures, which would be desirable for a cold tolerance screening program, neither FV/FM nor RFD were correlated with tolerance of the nine Arabidopsis accessions used in this study. Since both parameters are related to intactness of PSII, it is not unexpected that they are not affected by above-zero temperatures in a cold tolerant plant like Arabidopsis. However, it was noticed that parameters ФPSII, and qP that reflect effects of actinic light on chlorophyll fluorescence responded differentially to low above-zero and sub-zero temperature treatments. It has been reported that shifting plants from normal to low temperatures causes repression of genes related to photosynthesis [35], and among them are genes directly involved in light harvesting like e.g. proteins Lhca2*1 and Lhcb4*2 [36]. Down-regulation of photosynthesis genes is accompanied by a sudden suppression of photosynthesis at low temperatures [37]. However this might be a transient effect, since leaves that developed at low temperatures are reported to regain full photosynthetic activity [38]. Although a reduction in antenna size would not necessarily impact parameters ФPSII and qP, it is not self-evident why low temperatures should increase effective quantum yield. Savitch & co-workers [38] reported that low temperatures could increase photosynthetic capacity in species like e.g. wheat that maintain a high need for assimilates because of active growth at low temperatures. In the cold tolerant tomato species Lycopersicon peruvianum, a rise in qP during cold exposure was hypothesized to result from increased capacity of the Calvin-Benson cycle [39]. In this context it should be noted that temperature effects on light absorption and photosynthesis are often not easily separable from light effects, because in most studies light intensity is reduced at low temperature to avoid damage to the photosystems in the cold. At lower light intensity, qN is reduced and qP rises, and in turn also ФPSII is increased as long as no damage of the photosystems occurs. This effect may superimpose the temperature effect on light harvesting, and although the cold sensitive accession Cvi could be discriminated from all intermediate and tolerant accessions based on the lack of increase of ФPSII during acclimation, this parameter does not seem to be sufficient for screening cold tolerance in acclimated Arabidopsis plants.

Combinatorial imaging can discriminate cold acclimation induced cold tolerant plants

We have earlier reported that combinatorial imaging of ChlF transients combined with classifier and feature selection methods was able to discriminate detached leaves from cold sensitive, tolerant and intermediate accessions of Arabidopsis [10]. The main outcome of the current study is that this method can also be applied to whole plants, i.e. leaves to be compared need not be present on the same images, and thus this method is suitable for large scale screening of plant cold tolerance.

Combinatorial imaging combines sets of highly performing images from sequences of time-resolved ChlF images that provide strong discrimination capacity [10, 29, 30]. High performing images extract information from several thousands of pixels from hundreds of measured images (218 images per data for this experiment, and each image with resolution 512*512 pixels). This method does not trace underlying physiological phenomena rather its algorithms select image sets having optimal contrast between sensitive and tolerant features that ultimately allow discrimination on the basis of their cold tolerance levels [30]. Therefore, parameters identified for the classifiers are very likely species and treatment specific, and training of the classifiers needs to be performed for each application of the method. While this might appear disadvantageous, it also offers great flexibility, because successful discrimination is not dependent on any specific physiological process but can exploit every feature that allows discrimination [10, 29, 30, 40, 41]. In the current study, older leaves performed better than developing ones, probably because the higher fluorescence intensity provided a better signal-to-noise ratio. However, this does not implicate that the method is not suited to the screening of plant seedlings, because it analyzes relative fluorescence intensities distributed among the image pixels and for each sample type it is different [29]. Therefore, combinatorial imaging out-performs existing methods that depend on the application of freeze-thaw treatments.

The combinatorial imaging method could discriminate sensitive and tolerant accessions in the cold acclimated state but it failed to make any difference between cold sensitive and intermediate accessions. Because training was done for STT plants, this failure could result from differential responses of the photosynthetic apparatus to low above-zero and sub-zero temperatures as outlined above, or it could relate to insufficient resolution, which might be overcome by a comparison of larger sets of images. Since this would clearly result in higher computational load, a case specific decision for higher resolution vs. shorter analysis time might be necessary.

Even if only a two-class discrimination can be achieved with the combinatorial imaging method, it would be very useful for large scale screening of cold tolerance e.g. in recombinant inbred line populations or other plant sets consisting of large numbers of individuals representing different genotypes. The method is thus well suited for quantitative trait loci (QTL) mapping or mutant screenings to investigate genetic determinants of cold tolerance that may be used for plant breeding.


We have demonstrated that chlorophyll fluorescence emission from whole plant rosettes of Arabidopsis thaliana integrates information that can be used to discriminate cold sensitive and tolerant plants in the cold acclimated state when analyzed by the advanced statistics based combinatorial imaging approach [10, 29, 30]. This reveals the power of combinatorial imaging for identifying features of cold tolerant and sensitive accessions at cold acclimation state (Figure 5) where well known physiological parameters of ChlF emission (Figures 2 and 3) failed to provide any clue of discrimination. Moreover, capturing ChlF transients of whole plant rosettes following mild sub-zero temperature treatments (STT, -4°C for ~8 h in dark) is also very useful, because classical ChlF transients and the extracted parameters such as FV/FM and RFD can categorize cold sensitive, intermediate and tolerant accessions following STT. In addition, classical ChlF analysis, in contrast to combinatorial imaging methods, can yield physiologically relevant information that could be directly exploited for breeding efforts. The screening of whole plants cold tolerance via fully non-invasive means following cold acclimation, i.e. without any sub-zero temperature treatments, could be highly useful in high-throughput screening of cold tolerance where it is superior to data measured from single leaf or leaf discs by alternative methods such as metabolomics [7] or electrolyte leakage [5].


Plant material and growth conditions

Six plantlets of each of the nine accessions of Arabidopsis thaliana[5] were grown inside a cooling chamber in pots of 0.06 m for six weeks at day/night temperatures 20°C/18°C, light 90 μmol (photons) m-2 s-1 and a relative humidity of 70%. The accessions used in this study were: Cvi (Cap Verde Islands), Can (Canary Islands), C24 (genetically related to Co, Portugal), Co (Coimbra, Portugal), Col-0 (Columbia-0, genetically related to Gü, Germany), Nd (Niederzenz, Germany), Ler (Landsberg erecta, Poland), Rsch (Rschew, Russia), and Te (Tenela, Finland). The experiments were executed in three independent replica with three independent sets of plants used to measure ChlF in the non acclimated (NAC, six weeks grown plants), cold acclimated (AC, NAC plants acclimated at 4°C for another 2 weeks), and sub-zero temperature treated (STT, AC plants treated with mild sub-zero temperature of -4°C for 8 hours in dark) state.

Chlorophyll fluorescence measurement

Individual plants were used for ChlF measurements at room temperature using Handy FluorCam [; [23]. A short protocol of ~202 seconds modified according to Mishra et al. [10] was used for capturing time-resolved ChlF from plant rosettes. This protocol starts with the measurement of basal fluorescence (FO) and maximum fluorescence (FM) using measuring and saturating flashes. After a short dark period of ~20 seconds, actinic light of 40 μmol (photons) m-2 s-1 was applied to measure the fluorescence transients. Two strong flashes of saturating light were overlaid with actinic light to investigate activation of non-photochemical quenching followed by a third saturating flash 18 s after switching off the actinic light to see the relaxation of non-photochemical quenching mechanisms. Three replicate recordings were taken for NAC, AC and STT state, respectively. The images of the measured ChlF transients were averaged across the whole rosettes for quantitative evaluation of the fluorescence parameters or plotting of ChlF transients.

Combinatorial imaging using statistical classifier and feature selection method

The method of combinatorial imaging was applied to identify the discriminant between the accessions without using the whole data set consisting of 218 images in a time series. Each time series contains some repetitive images as well as images with low contrast, which were sorted out to reduce the size of the data sets. Application of statistical techniques of classifiers was chosen to avoid biasness. We randomly classified the data for Te and Co as training and testing sets. Using the training set, features discriminating the tolerant and sensitive accession were identified, and these features were next applied using the testing set for discriminating cold tolerant and sensitive accessions. Using this method we calculated the performance of several classifiers: LDC, QDC, k-NNC, NNC and NMC, and chose the best performing classifier for further analysis. The performance of each of the investigated classifiers was quantified by a number between 0–1: value‘0’ means random classification (1/2 of classifications into 2 equally represented classes correct, and another 1/2 is incorrect) and value ‘1’ meaning that the classifier was 100% successful [for details see [24]. In this experiments performance of QDC was 81% and computational time was less than 6 seconds. Therefore, QDC classifier was applied with sequential forward floating selection (SFFS) to reduce the number of images for effective classification [42]. The reduction is based on finding an image sub-set containing the most useful information for the visualization of highly contrasting features of cold tolerant vs. sensitive accessions and arranging the images in a form of descending order of their performance efficiency. The combinatorial images are developed by the combination of the three high performing images (x,y,z) multiplied with their coefficient (a,b,c) that are obtained by linear discriminant analysis (LDA): Combinatorial imaging (C) = (±a)*Ix + (±b)*Iy + (±c)*Iz.

Tool for the data analysis

The Matlab software package, version 6.5, with Pattern Reorganization toolbox (PRTools) was used for statistical analysis.



Chlorophyll a fluorescence




Cold acclimated


Sub-zero temperature (treated)


photosystem II


Photosystem I


Maximum quantum yield of PSII photochemistry, where FV is variable fluorescence (FV = FM-FO), FO is minimal fluorescence emission of a dark adapted plant with primary quinone acceptor (QA), oxidized and non-photochemical quenching inactive, and FM is maximum fluorescence emission of a dark adapted plant exposed to a short pulse of a strong light leading to a transient reduction of QA


Fluorescence decrease ratio, where FP is a fluorescence peak at the beginning of the transient with actinic light, and FS is steady state fluorescence


Quadratic discriminant classifier


Sequential forward floating selection


Electrolyte leakage


The half of maximal lethal temperature


Electrolyte leakage


Effective quantum efficiency of PSII photochemistry


Photochemical quenching


Non-photochemical quenching

k-NNC k:

- nearest neighbor classifier


Nearest neighbor classifier


Nearest mean classifier


Linear discriminant analysis.


  1. Jaglo-Otteson KR, Gilmour SJ, Zarka DG, Schabenberger O, Thomashow MF: Arabidopsis CBF1 overexpression induces COR genes and enhances freezing tolerance. Science. 1998, 280: 104-106. 10.1126/science.280.5360.104.

    Article  Google Scholar 

  2. Cook D, Fowler S, Fiehn O, Thomashow MF: A prominent role for the CBF cold response pathway in configuring the low-temperature metabolomes of Arabidopsis. Proc Natl Acad Sci U S A. 2004, 101 (45): 15243-15248.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Kaplan F, Kopka J, Haskell DW, Zhao W, Schiller KC, Gatzke N, Sung DY, Guy CL: Exploring the temperature-stress metabolome of Arabidopsis. Plant Physiol. 2004, 136: 4159-4168. 10.1104/pp.104.052142.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Guy C, Kaplan F, Kopka J, Selbig J, Hincha DK: Metabolomics of temperature stress. Physiol Plant. 2008, 132: 220-235.

    CAS  PubMed  Google Scholar 

  5. Hannah MA, Wiese D, Freund S, Fiehn O, Heyer AG, Hincha DK: Natural genetic variation of freezing tolerance in Arabidopsis. Plant Physiol. 2006, 142: 98-112. 10.1104/pp.106.081141.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Zuther E, Schulz E, Childs LH, Hinch DK: Clinal variation in the non-acclimated and cold-acclimated freezing tolerance of Arabidopsis thaliana accessions. Plant Cell Environ. 2012, 35: 1860-1878. 10.1111/j.1365-3040.2012.02522.x.

    Article  CAS  PubMed  Google Scholar 

  7. Vaclavik L, Mishra A, Mishra KB, Hajslova J: Mass spectrometry-based metabolomic fingerprinting for screening cold tolerance in Arabidopsis thaliana accessions. Anal Bioanal Chem. 2013, 405 (8): 2671-2683. 10.1007/s00216-012-6692-1.

    Article  CAS  PubMed  Google Scholar 

  8. Knaupp M, Mishra KB, Nedbal L, Heyer AG: Evidence for a role of raffinose in stabilizing photosystem II during freeze–thaw cycles. Planta. 2011, 234 (3): 477-486. 10.1007/s00425-011-1413-0.

    Article  CAS  PubMed  Google Scholar 

  9. Moellering ER, Bagyalakshmi M, Benning C: Freezing tolerance in plants requires lipid remodeling at the outer chloroplast membrane. Science. 2010, 330: 226-228. 10.1126/science.1191803.

    Article  CAS  PubMed  Google Scholar 

  10. Mishra A, Mishra KB, Höermiller II, Heyer AG, Nedbal L: Chlorophyll fluorescence emission as a reporter on cold tolerance in Arabidopsis thaliana accessions. Plant Signal Behav. 2011, 6 (2): 301-310. 10.4161/psb.6.2.15278.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Degenkolbe T, Giavalisco P, Zuther E, Seiwert B, Hinchey DK, Willmitzer L: Differential remodeling of the lipidome during cold acclimation in natural accessions of Arabidopsis thaliana. Plant J. 2012, 72: 972-982.

    CAS  PubMed  Google Scholar 

  12. Rapacz M, Gasior D, Zwierzykowski Z, Lesniewska-Bocianowska A, Humphreys MW, Gay AP: Changes in cold tolerance and the mechanisms of acclimation of photosystem II to cold hardening generated by anther culture of Festuca pratensis x Lolium multiflorum cultivars. New Phytol. 2004, 162 (1): 105-114. 10.1111/j.1469-8137.2004.01024.x.

    Article  CAS  Google Scholar 

  13. Novillo F, Alonso JM, Ecker JR, Salinas J: CBF2/ DREB1C is a negative regulator of CBF1/DREB1B and BF3/DREB1A expression and plays a central role in stress tolerance in Arabidopsis. Proc Natl Acad Sci U S A. 2004, 101: 3985-3990. 10.1073/pnas.0303029101.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Hincha DK: Cryoprotectin: a plant lipid-transfer protein homologue that stabilizes membranes during freezing. Philos T R Soc B. 2002, 357: 909-915. 10.1098/rstb.2002.1079.

    Article  CAS  Google Scholar 

  15. Hincha DK, Hofner R, Schwab KB, Heber U, Schmi JM: Membrane rupture is the common cause of damage to chloroplast membranes in leaves injured by freezing or excessive wilting. Plant Physiol. 1987, 83: 251-253. 10.1104/pp.83.2.251.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Steponkus PL, Lynch DV, Uemura M: The influence of cold-acclimation on the lipid-composition and cryobehavior of the plasma-membrane of isolated rye protoplasts. Philos T R Soc B. 1990, 326: 571-583. 10.1098/rstb.1990.0032.

    Article  CAS  Google Scholar 

  17. Azzarello E, Mugnai S, Pandolfi C, Masi E, Marone E, Mancuso S: Comparing image (fractal analysis) and electrochemical (impedance spectroscopy and electrolyte leakage) techniques for the assessment of the freezing tolerance in olive. Trees. 2009, 23: 159-167. 10.1007/s00468-008-0264-1.

    Article  Google Scholar 

  18. Yamori W, Noguchi K, Terashima I: Temperature acclimation of photosynthesis in spinach leaves: analyses of photosynthetic components and temperature dependencies of photosynthetic partial reactions. Plant Cell Environ. 2005, 28: 536-547. 10.1111/j.1365-3040.2004.01299.x.

    Article  CAS  Google Scholar 

  19. Govindjee: Sixty-three years since Kautsky: chlorophyll a fluorescence. Aust J Plant Physiol. 1995, 22: 131-160. 10.1071/PP9950131.

    Article  CAS  Google Scholar 

  20. Baker N: Chlorophyll fluorescence: a probe of photosynthesis in vivo. Annu Rev Plant Biol. 2008, 59: 89-113. 10.1146/annurev.arplant.59.032607.092759.

    Article  CAS  PubMed  Google Scholar 

  21. Agati G, Cerovic ZG, Moya I: The effect of decreasing temperature up to chilling values on the in vivo F685/F735 chlorophyll fluorescence ratio in Phaseolus vulgaris and Pisum sativum: The role of the photosystem I contribution to the 735 nm fluorescence band. Photochem Photobiol. 2000, 72: 75-84. 10.1562/0031-8655(2000)072<0075:TEODTU>2.0.CO;2.

    Article  CAS  PubMed  Google Scholar 

  22. Baker NR, Rosenqvist E: Applications of chlorophyll fluorescence can improve crop production strategies: an examination of future possibilities. J Exp Bot. 2004, 55 (403): 1607-1621. 10.1093/jxb/erh196.

    Article  CAS  PubMed  Google Scholar 

  23. Ehlert B, Hincha DK: Chlorophyll fluorescence imaging accurately quantifies freezing damage and cold acclimation responses in Arabidopsis leaves. Plant Methods. 2008, 4: 12-10.1186/1746-4811-4-12.

    Article  PubMed Central  PubMed  Google Scholar 

  24. Thalhammer A, Hincha DK, Zuther E: Measuring freezing tolerance: electrolyte leakage and chlorophyll fluorescence assays. Plant Cold Acclimation: Methods Mol Biol. 2014, 1166: 15-24. 10.1007/978-1-4939-0844-8_3.

    Article  CAS  Google Scholar 

  25. Strasser RJ, Srivastava A, Govindjee: Polyphasic chlorophyll-Alpha fluorescence transcient in plants and cyanobacteria. Photochem Photobiol. 1995, 61 (1): 32-42. 10.1111/j.1751-1097.1995.tb09240.x.

    Article  CAS  Google Scholar 

  26. Rapacz M, Gasior D, Koscielniak J, Kosmala A, Zwierzykowski Z, Humphreys MW: The role of the photosynthetic apparatus in cold acclimation of Lolium multiflorum. Characteristics of novel genotypes low-sensitive to PSII over-reduction. Acta Physiol Plant. 2007, 29 (4): 309-316. 10.1007/s11738-007-0040-7.

    Article  CAS  Google Scholar 

  27. Rapacz M, Wozniczka A: A selection tool for freezing tolerance in common wheat using the fast chlorophyll a fluorescence transient. Plant Breed. 2009, 128 (3): 227-234. 10.1111/j.1439-0523.2008.01566.x.

    Article  Google Scholar 

  28. Nedbal L, Soukupová J, Kaftan D, Whitmarsh J, Trtílek M: Kinetic imaging of chlorophyll fluorescence using modulated light. Photosynth Res. 2000, 66: 3-10.1023/A:1010729821876.

    Article  CAS  PubMed  Google Scholar 

  29. Mishra A, Matouš K, Mishra KB, Nedbal L: Towards discrimination of plant species by machine vision: advanced statistical analysis of chlorophyll fluorescence transients. J Fluoresc. 2009, 19: 905-913. 10.1007/s10895-009-0491-x.

    Article  CAS  PubMed  Google Scholar 

  30. Matouš K, Benediktyova Z, Berger S, Roitsch T, Nedbal L: Case study of combinatorial imaging: what protocol and what chlorophyll fluorescence image to use when visualizing infection of Arabidopsis thaliana by pseudomonas syringae?. Photosynth Res. 2006, 90: 243-253.

    Article  PubMed  Google Scholar 

  31. Maxwell K, Johnson GN: Chlorophyll fluorescence: a practical guide. J Exp Bot. 2000, 51: 659-668. 10.1093/jexbot/51.345.659.

    Article  CAS  PubMed  Google Scholar 

  32. Papageorgiou GC, Govindjee: Photosystem II fluorescence: slow changes – scaling from the past. J Photoch Photobio B. 2011, 104 (1–2): 258-270.

    Article  CAS  Google Scholar 

  33. Lichtenthaler HK, Langsdorf G, Lenk S, Buschmann C: Chlorophyll fluorescence imaging of photosynthetic activity with flash-lamp fluorescence imaging system. Photosynthetica. 2005, 43: 355-369. 10.1007/s11099-005-0060-8.

    Article  CAS  Google Scholar 

  34. Horgan DB, Zabkiewicz JA: Fluorescence decline ratio: comparison with quantum yield ration for plant physiological status and herbicide treatment responses. New Zealand Plant Protection. 2008, 61: 169-173.

    CAS  Google Scholar 

  35. Hannah MA, Heyer AG, Hincha DK: A global survey of gene regulation during cold acclimation in Arabidopsis thaliana. PLoS Genet. 2005, 1 (2): e26-10.1371/journal.pgen.0010026.

    Article  PubMed Central  PubMed  Google Scholar 

  36. Fowler F, Thomashow MF: Arabidopsis transcriptome profiling indicates that multiple regulatory pathways are activated during cold acclimation in addition to the CBF cold response pathway. Plant Cell. 2002, 14 (8): 1-

    Article  Google Scholar 

  37. Strand A, Hurry V, Gustafsson P, Gardeström P: Development of Arabidopsis thaliana leaves at low temperature releases the suppression of photosynthesis and photosynthetic gene expression despite the accumulation of soluble carbohydrates. Plant J. 1997, 12: 605-614. 10.1046/j.1365-313X.1997.00605.x.

    Article  CAS  PubMed  Google Scholar 

  38. Savitch LV, Leonardos ED, Krol M, Jansson S, Grodzinski B, Huner NPA, Oquist G: Two different strategies for light utilization in photosynthesis in relation to growth and cold acclimation. Plant Cell Environ. 2002, 25: 761-771. 10.1046/j.1365-3040.2002.00861.x.

    Article  CAS  Google Scholar 

  39. Linger P, Brüggemann W: Correlations between chlorophyll fluorescence quenching parameters and photosynthesis in a segregating Lycopersicon esculentum × L. peruvianum population as measured under constant conditions. Photosynth Res. 1999, 61: 145-156. 10.1023/A:1006241729948.

    Article  CAS  Google Scholar 

  40. Codrea MC, Hakala-Yatkin M, Kårlund-Marttila A, Nedbal L, Aittokallio T, Nevalainen OS, Tyystjärvi E: Mahalanobis distance screening of Arabidopsis mutants with chlorophyll fluorescence. Photosynth Res. 2010, 105 (3): 273-83. 10.1007/s11120-010-9578-0.

    Article  CAS  PubMed  Google Scholar 

  41. Mattila H, Valli P, Pahikkala T, Teuhola J, Nevalainen OS, Tyystjärvi E: Comparison of chlorophyll fluorescence curves and texture analysis for automatic plant identification. Precis Agric. 2013, 14: 621-636. 10.1007/s11119-013-9320-y.

    Article  Google Scholar 

  42. Jain AK, Duin RPW, Mao JC: Statistical pattern recognition: a review. IEEE Trans Pattern Anal Mach Intell. 2000, 22 (1): 4-37. 10.1109/34.824819.

    Article  Google Scholar 

Download references


Participation of AM and KBM was supported by the EfCOP - IPo project ENVIMET (CZ.1.07/2.3.00/20.0246) and by the EU (FP7) EPPN project (284443). We thank Ladislav Nedbal for his inspiration.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Kumud B Mishra.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

KBM, AGH and AM had planed and designed the experiments. KBM and AM performed the experiments. AM analyzed the data and drafted the manuscript. KBM and AGH contributed in the final version of the manuscript. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mishra, A., Heyer, A.G. & Mishra, K.B. Chlorophyll fluorescence emission can screen cold tolerance of cold acclimated Arabidopsis thaliana accessions. Plant Methods 10, 38 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: