- Methodology article
- Open Access
Automated analysis of calcium spiking profiles with CaSA software: two case studies from root-microbe symbioses
BMC Plant Biology volume 13, Article number: 224 (2013)
Repeated oscillations in intracellular calcium (Ca2+) concentration, known as Ca2+ spiking signals, have been described in plants for a limited number of cellular responses to biotic or abiotic stimuli and most notably the common symbiotic signaling pathway (CSSP) which mediates the recognition by their plant hosts of two endosymbiotic microbes, arbuscular mycorrhizal (AM) fungi and nitrogen fixing rhizobia. The detailed analysis of the complexity and variability of the Ca2+ spiking patterns which have been revealed in recent studies requires both extensive datasets and sophisticated statistical tools.
As a contribution, we have developed automated Ca2+ spiking analysis (CaSA) software that performs i) automated peak detection, ii) statistical analyses based on the detected peaks, iii) autocorrelation analysis of peak-to-peak intervals to highlight major traits in the spiking pattern.
We have evaluated CaSA in two experimental studies. In the first, CaSA highlighted unpredicted differences in the spiking patterns induced in Medicago truncatula root epidermal cells by exudates of the AM fungus Gigaspora margarita as a function of the phosphate concentration in the growth medium of both host and fungus. In the second study we compared the spiking patterns triggered by either AM fungal or rhizobial symbiotic signals. CaSA revealed the existence of different patterns in signal periodicity, which are thought to contribute to the so-called Ca2+ signature.
We therefore propose CaSA as a useful tool for characterizing oscillatory biological phenomena such as Ca2+ spiking.
As a ubiquitous second messenger, calcium (Ca2+) mediates multiple signal transduction pathways in diverse types of plant responses to biotic and abiotic stimuli [1, 2]. Transient variations in calcium concentration occurring in the cytosol, nucleus and/or other compartments of the plant cell are believed to transduce extracellular signals into appropriate cellular responses . These range from changes in turgor pressure as in stomatal guard cells [4, 5], to the control of apical growth in pollen tubes  or the regulation of gene expression . One of the most intriguing and well-studied types of Ca2+ signaling is the generation of repeated peaks in Ca2+ concentration that can persist for relatively long periods of time (minutes to hours) and is commonly referred to as “Ca2+ spiking” [1, 4, 8].
Ca2+ spiking is characterized by distinctive features such as peak shape, amplitude, frequency and regularity, certain of which may confer a degree of specificity in the respective stimulus-response coupling. Amongst the best studied examples, the regulation of stomatal closure was demonstrated to be based on peak frequency , and in the case of host-endosymbiont signaling differences in Ca2+ spiking regularity may be important in the recognition of the bacterial or fungal partner [9, 10].
Unfortunately, the detailed characterization of such traits, which contribute to the so-called ‘calcium signature’ of the respective signal transduction pathway , is complicated by the intrinsic cell-to-cell variability of the oscillatory response [9, 10, 12, 13]. In addition, background noise often masks ‘true’ peaks with clusters of smaller random oscillations. For these reasons, Ca2+ spiking analyses always require large data-sets before general trends can be clearly established. As a result, manual large scale analyses and the application of statistical tools is often very time consuming and can lead to the introduction of human operator biases. Fully standardizing the analysis is a therefore a prerequisite for correctly analyzing Ca2+ spiking signatures.
Indeed, many analytical methods have been developed in other areas, such as neuroscience, to analyze and classify the repeated oscillations in a biological parameter such as membrane potential (e.g. [12–15]). Nevertheless, such analytical processes largely relate peak identification to their regular frequency. Since an intrinsic trait in symbiotic Ca2+ spiking signals in plants is a variable extent of irregularity [9, 12], such analyses based on peak prediction could not be successfully applied to our experimental data. We therefore chose to try a different approach and analytically characterize the phenomenon by identifying each spike based on its shape and how much this differs from background noise oscillations.
With this aim, we have prototyped original software that performs an automated Ca2+ spiking analysis (CaSA) based on our recent studies where a prototypal computational model was formalized to simulate Ca2+ spiking dynamics . The CaSA software that we present here can quickly and reliably analyze large data sets, independently of the instrument or method used for recording the biological signal, and perform automatic peak detection and several statistical analyses. A pipeline chart schematizing data preparation and CaSA usage can be found in Additional file 1. The CaSA software is coded in GNU OCTAVE , a high-level interpreted language, and can be run on GNU/Linux systems.
To evaluate the software reliability we have applied CaSA to nuclear Ca2+ spiking profiles obtained with root epidermal cells of the legume Medicago truncatula. We have chosen two case studies, both based on the induction of nuclear Ca2+ spiking by signals released from two symbiotic root-interacting microbes: the arbuscular mycorrhizal (AM) fungus Gigaspora margarita and the symbiotic nitrogen-fixing (SNF) bacterium Sinorhizobium meliloti. Legumes can establish mutualistic interactions with both types of root symbionts. The AM symbiosis develops inside the root tissues, where living cortical cells are colonized by specialized hyphal structures called arbuscules. In contrast, SNF rhizobia colonize specialized root-derived organs, called root nodules, where they convert atmospheric nitrogen into ammonia, that can be later assimilated into organic compounds that can be assimilated by the host .
For both AM and SNF, initial stages of the symbiotic association require reciprocal recognition between the host plant and the respective microbe. Rhizobia are identified by the host plant through the perception of secreted lipochito-oligosaccharides (LCOs) known as “Nod factors” . Recent evidence suggests that the molecules mediating the recognition of AM fungi may also be based on chito-oligosaccharides, either LCOs (so-called “Myc LCOs”)  or undecorated, short-chain chito-oligosaccharides such as chito-tetraose (CO4) . Downstream of their perception, these chitinaceous signals are transduced by a common symbiotic signalling pathway (CSSP) involving a common subset of plant proteins. One of these proteins, localized in the nucleoplasm, is a calcium- and calmodulin-dependent kinase (MtDMI3 in M. truncatula) . This key component of the CSSP is believed to interpret the respective nuclear Ca2+ spiking response which is activated following symbiotic signal perception [10, 22, 23].
AM fungi are beneficial to the plant since thay are able to extract inorganic phosphate (Pi) and other mineral nutrients from the soil with a higher efficiency compared to the non-symbiotic root system . Consistent with this, high Pi availability in the soil is well known to limit the overall level of root colonization by AM fungi [25–27], although the underlying molecular mechanisms are only starting to emerge . Furthermore, Pi availability has been shown to reduce the plant biosynthesis and secretion of strigolactones, key regulators of presymbiotic fungal development (for a review, ).
In our first case study we have investigated whether the Pi concentration in the growth medium also impacts the earliest plant responses to AM fungal signals by monitoring nuclear Ca2+ oscillations in the host epidermis. Our results using CaSA analysis have allowed us to detect a significant inhibition of presymbiotic signaling in response to high Pi levels.
How nuclear Ca2+ spiking specifically encodes AM- or SNF-related signals within the CSSP, and how these are then decoded by DMI3 to activate the appropriate downstream responses remains unclear. One possibility is that the pattern of the spiking profiles contains the encoded information in the form of a specific Ca2+ signature. For this reason the spiking profiles induced in legumes by both AM fungi and rhizobia have been the subject of several studies [9, 12]. In our second case study we have applied the CaSA software to analyze the spiking profiles induced by either Nod factor or the putative AM CO4 signals in the epidermis of M. truncatula roots. Our results underline highly significant differences in the regularity of peak-to-peak intervals, suggesting that this trait could be part of the signature used by the plant cell to discriminate between the two signals.
The CaSA software
CaSA software comprises a workflow of programming modules to analyze time-series data. It was prototyped in the GNU Octave language on the GNU/Linux platform. To run the GNU Octave software, the Gnuplot graphic utility and Zenity/GTK + tools are required.
The computation flow, presented in Figure 1, proceeds as follows:
The treatment and control data sets, loaded by the user in the form of text files, are filtered for noise using a polynomial SG filter. The same filtering also provides point by point first and second derivatives. The filtering is controlled by the parameters listed below. In order to fine tune the analysis to the specific characteristics of different experimental systems, all of the parameters can be adjusted by the user via dialogue windows (Additional file 2 and Additional file 3).
(w) Width of the sliding measurement window. This must be an odd value, corresponding to the total number of time points that will be included in the filter window. If n points are considered before and after the filtered point (see Methods), then w = 2 n + 1. In general, a window width covering the maximum peak duration should be selected.
(k) Degree of the approximating polynomial curve. A default value of 5 (odd) was chosen to capture peak asymmetry, which is a constant feature in our spiking records, but different values can be entered by the user to best fit any peak shape.
(τ) Threshold for peak detection in the function peak c (i). This value determines the threshold in peak height above which an oscillation is identified as a peak rather than background noise. A default threshold is proposed by the software, but again the user can enter a value that provides the best signal/noise discrimination for a specific data set.
The application of the derived peak function to the filtered data is used to obtain an instantaneous indicator of peak occurrence. In the case of the control data set, where only background noise is recorded, this function identifies a number of ‘false’ peaks that will be used in step 3 to validate the significance of peaks detected for the treatment data.
The amplitude of each peak in the treatment data set is compared to the amplitude of ‘false’ peaks identified in the control data set. In brief, a peak is validated whenever its amplitude is higher than that measured for a defined percentage of the false peaks from the control data set. Such a percentage is defined as a quantile (ranging between 0 and 1) and can be chosen by the user in the dialog window where the q parameter is requested. A default value of 0.9 (corresponding to 90%) is automatically proposed by the software.
The following quantitative information is then obtained for each data set: number of peaks detected, peak-to-peak intervals (waiting time). Actively-responding cells are defined as those with a minimum of three peaks during the 30 min period of observation. Basic statistical analyses are carried out on the output data, including the average waiting time and the waiting time autocorrelation for each sample, as well as the average peak number in the actively responding cells from the whole data set.
After the analysis, the CaSA software outputs a series of plots as PNG image files, where the original ratiometric data, the SG filtered curve, and each of the detected peaks (including the starting points, maximum and ending points) are indicated. Figure 2A and B show representative profiles recorded from cells that were treated with either the fungal exudate or sterile distilled water (as control). Figure 2C and D show the results of CaSA software analysis of this data. The original time series (in red) is overlayed onto the SG filtered curve (blue) and the starting points, maximum and ending points are indicated for each detected peak.
An output text file is also created, containing all the results of the computational analysis (an example is presented as Additional file 4), such as the peak count for each sample, the identification of “active” vs. “inactive” cells (based on the presence of at least 3 peaks) and the autocorrelation analysis of peak-to-peak intervals (waiting times). The file includes a legend briefly explaining the meaning of each quantitative parameter.
The CaSA software is available for download with the supplementary data (Additional file 5).
Case study 1: Ca2+spiking responses to AM fungal exudates as a function of the Pi concentration in the fungal and plant growth media
With the aim of studying whether Pi availability can influence the early pre-infection steps of the AM interaction, and in particular the activation of Ca2+ signaling within the CSSP, we recorded nuclear Ca2+ spiking profiles in epidermal tissues of M. truncatula roots in response to germinated AM fungal exudates when both roots and fungus were grown under different phosphate conditions, as described in the Methods section and reported in Table 1.
By applying the CaSA software to our data we were able to calculate different parameters of the spiking response. Firstly, the percentage of ‘actively responding cells (defined as those displaying at least three peaks over the 30 min measurement period) was calculated for each treatment. As reported in Figure 3, statistical tests showed that the percentage of responding cells when both the plant roots and the fungus were grown under low phosphate conditions (the standard growth conditions used for in vitro AM development) was significantly higher compared to all other experimental combinations (p < 0.05, using unpaired parametric Student t test and Wilcoxon- Mann–Whitney non parametric test). The average number of peaks in the active cells, directly derived from the CaSA output file, was also used for comparison. Figure 4 shows that total peak averages were statistically higher when both roots and fungus were grown under low Pi (exudate low phosphate - medium low phosphate; ELP-MLP) and statistically lower when both roots and fungus were grown under high Pi (exudate high phosphate - medium high phosphate; EHP-MHP), confirmed by both Kruskal-Wallis (p = 5.4 • 10-3) and post-hoc tests (p = 3.9 • 10–3). Since the high Pi conditions used in our experiments did not affect fungal viability or spore germination rate (not shown), our results indicate that an inhibitory effect of Pi is already present during the earliest stages of the interaction that involve the activation of host nuclear Ca2+ spiking in response to fungal symbiotic signals present in germinating spore exudates. This inhibition is highlighted by the lower intensity of the Ca2+ spiking response (in terms of both the percentage of responding cells and average number of peaks) and is most striking when both the plant and the fungus are grown in high Pi media.
Case study 2: Comparison of Ca2+spiking profiles in response to AM fungal and rhizobial signals
This case study aimed at identifying quantifiable features in the spiking patterns triggered by either Nod factors or CO4, the signal molecules that most effectively activate SNF- and AM-related nuclear Ca2+ signaling . Since early SNF and AM responses focus on different cell types - root hairs versus atrichoblasts respectively - in this study we have compared root hair responses to Nod factors with atrichoblast responses to CO4. Several examples of typical spiking patterns observed in M. truncatula root epidermal cells in response to S. meliloti Nod factor or CO4 are shown in Figure 5. We found that CaSA was able to quantify significant and consistent differences between the two populations of spiking data sets. It turned out that the most striking differences concerned the time intervals between adjacent peaks (waiting time). By applying autocorrelation analysis to these values, the software identified a different periodicity, thus providing a quantitative parameter that discriminates the two spiking responses. As presented in Figure 6, the distribution of autocorrelation values is significantly different for CO4 and Nod factor treatments (Pearson χ2 test, p = 2.15 • 10-6). This indicates that the sequence of peak-to-peak intervals in response to Nod factor is relatively regular: long waiting times are in most cases followed by long waiting times and short waiting times are followed by short waiting times. By contrast, the CO4 treatment elicits spiking which is generally characterised by negative autocorrelation values, indicating that most cells alternate short and long waiting times.
For comparative purposes, we have applied the same CaSA analysis to a data set of Ca2+ spiking responses induced in the M. truncatula epidermis by G. margarita exudates prepared in the absence of phosphate (see ELP condition in case study 1). Interestingly, the distribution of waiting time autocorrelation values is very similar to that for CO4-induced spiking, with no statistically significant difference (Pearson χ2 test, p = 0,01). This confirms the fact that short-length chitin oligomers can mimic the effect of AM fungal exudates in activating the CSSP and strengthens the proposal that these molecules are key signals in glomeromycete recognition by their hosts .
In conclusion, the CaSA software has allowed us to quantify significant differences in peak periodicity between the spiking patterns associated with AM fungal signaling (CO4 or fungal exudate) and rhizobial signaling (Nod factor).
The characterization of the molecular dialogues between symbiotic microbes and their plant hosts is of crucial importance in understanding the key recognition processes which precede successful associations. In the case of the endosymbiotic interactions between legumes and either glomeromycetes or rhizobia, research in many laboratories has to a large extent identified the genetic bases of the host mechanisms involved in the perception of fungal and bacterial signals, leading to the characterization of a common transduction pathway, the CSSP [30, 31]. It has been proposed that distinct signatures exist for the nuclear Ca2+ spiking which is at the core of the CSSP [30, 32]. Given the crucial role of Pi in the symbiotic nutrient trading between Glomeromycota and their host plants , growing interest is accumulating around the effects of available Pi on the establishment of AM [27, 34, 35]. The investigation of such processes has revealed a more complex picture than previously envisaged, suggesting the existence of a crosstalk in the plant perception of microbial signals and available nutrients.
Ca2+ spiking is the earliest known hallmark of symbiotic plant-microbe interactions and, as an increasing number of research groups focus on these topics, the variability and complexity of the observed spiking profiles is now becoming apparent. As a result, there is a crucial need for reliable tools that allow fast and reproducible analyses of spiking patterns. We propose the CaSA software as one such tool, and the examples we have chosen to illustrate its usefulness, targeted on two specific biological problems, and based on relatively large data sets, show that this approach can indeed identify and quantify distinct features of spiking profiles which can be used to define Ca2+ signatures.
CaSA software performance
Compared to manual peak count and data analysis, the use of CaSA analysis resulted to be advantageous in terms of both time and reliability. Although a few file conversion steps are also required when using CaSA, the software dramatically reduced the time required for data analysis and figure preparation. We estimated that, excluding image acquisition and confocal data import to a datasheet file, 30 minutes were required for manual peak counting in a set of 50 nuclei. CaSA reduced this time to less than 5 minutes, saving more than 80% of operator time. Moreover, CaSA also provides a few statistical analyses that would further increase the required time if calculated manually.
Even if CaSA can be considered as a prototypal software in its present release, a major effort has been made to simplify usage while at the same time providing direct access to all of the major computational functions implemented in the automated spiking analysis process. For this, the user is presented with a limited number of dialog windows where the key parameters of the analysis can be entered. In all cases default values are proposed which correspond to those that were the most suitable for our studies. Nevertheless, by adjusting each parameter, the analysis can be fine-tuned to very diverse types of data. For example, by changing the degree of the approximating polynomial from 5 to 2, the software can better detect peaks where asymmetry is less pronounced or totally absent.
In addition to the generation of a text file containing all the results of the analysis, the last step in the CaSA process includes the rendering of a series of image files where the data are plotted. This graphical representation of the analysis has many applications. First of all it is of great help when adjusting the software parameters to better fit the spiking records. The identification of errors in the procedure is extremely easy since each peak detected can be directly checked with the overlaid original data. Secondly, a rapid screening of the output image files allows the identification of individual records which may contain instrument-related mistakes or other anomalies. Such records can be easily excluded from the subsequent computation runs by entering their corresponding progressive numbers into the final dialog window entitled ‘Samples to be excluded’. This option is also of use when a subset of plots is needed, for example in view of data publication.
As for all software, CaSA analysis is not error free. In our tests, false positive and negative peaks did occasionally appear even when the analysis parameters were adjusted to give optimal performance. Nevertheless, the impact of such rare errors did not modify the final interpretation of any of the global trends that we considered in our case studies. In addition, repeated trials where we compared manual and automated peak counting revealed that manual peak detection varies significantly from one operator to another, and even when the same operator repeats the counting under different conditions. By contrast, CaSA provides the obvious advantage of full repeatability of the analysis. A typical case where manual peak identification could be ambiguous can be found in the second plot of Figure 5B, where an irregularly shaped peak (marked by an asterisk) displays an initial maximum followed by a second weaker oscillation. In our tests, different operators identified either one or two peaks in this shape. By contrast only the first maximum is consistently identified as a peak by CaSA.
In conclusion we propose the CaSA software as a versatile tool that can be of help to obtain rapid and reliable analyses of diverse types of oscillatory data, including but not limited to Ca2+ spiking.
Pi partially inhibits the presymbiotic dialogue in AM associations
A number of studies have analysed the physiological and molecular mechanisms that regulate the establishment of the AM symbiosis, highlighting the inhibitory effect of high phosphate availability on the extent of root colonization by AM fungi [26–28, 34, 36]. Here we have investigated whether the presence of high Pi in the growth media also affects earlier steps of the interaction, and in particular activation of the CSSP during the pre-symbiotic dialogue. Our analyses using CaSA suggest two distinct effects of high Pi concentration. Firstly, when both symbionts are independently grown under high Pi conditions, the root responsiveness to the fungal exudate is significantly reduced both in terms of the percentage of epidermal cells where the nuclear Ca2+ spiking signal is triggered and the total number of peaks recorded in responding cells. We interpret this as a partial inhibitory effect of Pi on the activation of the CSSP. This is a novel observation that suggests that at least part of the known inhibition of root colonization by AM fungi under comparable Pi concentrations [27, 34] may be due to a reduced pre-symbiotic host response. In further support of this, high Pi has also been shown to moderately lower the expression levels of Petunia genes acting upstream of Ca2+ spiking in the CSSP .
Secondly, when high Pi is present only during AM spore germination there is also a significantly weaker plant response in terms of average peak number compared to the reference low Pi conditions. This suggests that high Pi levels not only affect the plant host, but are also perceived by the fungus. The ability of the AM fungus to detect extracellular Pi availability is perhaps not surprising, if we consider that glomeromycetes are particularly effective in Pi scavenging from soils. The implication that the AM fungus could in some way limit its own signaling to the host plant under Pi-rich conditions is extremely intriguing and certainly deserves further investigation. In line with this,  recently demonstrated that Pi availability modulates the expression of several fungal genes, including the Pi transporter GintPT in the AM fungus Rhizophagus irregularis. It should be underlined that our reference ELP-MLP condition (35 μM Pi in the root culture medium and 0 μM in the spore germination medium) corresponds to the standard growth conditions that were found to be optimal for in vitro AM establishment [37, 38]. The identical condition, where spore germination was induced in sterile distilled water, was also recently used to characterize the potential signaling molecules (chito- and lipochito-oligosaccharides) produced by the presymbiotic mycelium of AM fungi [10, 13, 20].
A recent publication  also investigated the effects of Pi availability on CSSP activation with a significantly different setup. Part of the study was focused on a later stage of the interaction - after hyphopodium contact. Secondly, fungal exudates used to induce Ca2+ spiking were only prepared in the absence of Pi and from a different fungal species (Rhizophagus irregularis): we currently cannot exclude that the complex connection between AM development and Pi availability is also dependent on the plant/fungus coupling.
We therefore conclude that Pi has a limited but significant effect on host plant perception/transduction of AM signals present in fungal exudates from G. margarita. This attenuated Ca2+ spiking response in the presence of high Pi appears to be due both to fungal and host partners. Since the optimal synthesis and secretion of strigolactones, the root-secreted molecules which stimulate pre-symbiotic hyphal branching in AM fungi , also requires Pi deficiency , we can now consider that high Pi functions as a general negative regulator of the molecular pre-symbiotic dialogue.
Specific Ca2+signatures for AM and SNF?
Ever since the discovery that both AM fungal and rhizobial signals are transduced by the CSSP in legumes , the question immediately arose as to the mechanisms by which two different signals could be transduced by the same pathway and activate different downstream responses . Upstream of the CSSP, different receptors are probably involved in binding the respective AM fungal or rhizobial signals. Although a specific receptor for fungal signals has not yet been identified, this hypothesis is supported by the fact that putative Nod factor receptors (such as MtNFP) are not required for AM establishment . Similarly, downstream of the CSSP, different transcription factors are involved in gene regulation during AM (MtRAM1) and SNF (MtNSP1). It has also been proposed that a central and essential component of the CSSP, the nuclear kinase MtDMI3, is differentially activated when binding Ca2+ ions and/or calmodulin . This would indeed be in line with the generation of different calcium signals encoding AM- or SNF-specific information. The spiking profiles induced in SNF and AM have both been described as combining chaotic and stochastic elements [9, 12], while part of the response-specificity could be related by the type of epidermal cell responding to either AM fungi or rhizobia: atrichoblasts and root hairs respectively .
The study reported here is based on the hypothesis of  concerning the presence of a stochastic component in AM- and SNF-related Ca2+ spiking. Stochasticity can be interpreted as evidence for the transmission of a message through a communication channel . In this context, we have previously demonstrated the exponential distribution of waiting times in relation to the elapsed time from the initial stimulus , in agreement with the theoretically efficient use of a communication channel .
We therefore investigated the presence of a repetitive pattern in the sequence of peak-to-peak intervals, as a possible component of the calcium signature. Several techniques of pattern recognition can be found in literature which in most cases analyze the conformity of the experimental data to a predetermined pattern model . Since we had no indication for a precise pattern model that could fit our case, we introduced the calculation of waiting time autocorrelation. This allowed the characterization two distinct trends in the spiking patterns associated with either AM fungi or rhizobia.
Based on the results presented in this article, we now propose that the explicit mechanism encoding the Ca2+-mediated message can be found in the variability of the waiting time. By introducing the calculation of waiting time autocorrelation we have been able to characterize two distinct trends in the spiking patterns associated with either AM fungi or rhizobia. In line with the hypothesis of , these results provide evidence for a quantifiable trait in nuclear Ca2+ signaling that can play a role in the capacity of the legume root epidermis to discriminate between the two symbionts.
The CaSA software has been exploited in two different experimental contexts to further our understanding of the molecular/cellular signaling mechanisms underlying the establishment of symbiotic plant-microbe interactions. We therefore propose CaSA as an attractive platform to face so far unsolved questions concerning Ca2+ spiking signals and their modulation. Our findings suggest that AM fungal signal production could be limited in the presence of high inorganic phosphate levels, Since CO4 production is significantly enhanced by strigolactones , it would now be interesting to study strigolactone stimulation of the germinating AM spores as a function of different Pi concentrations. CaSA software has also revealed substantial, quantifiable differences between the Ca2+ spiking profiles triggered by either Nod factor or CO4, thus opening the question as to whether such differences remain when roots are treated with Myc-LCOs, whose chemical structures are more closely related to Nod factor .
Plant and fungal materials
Medicago truncatula genotype Jemalong A17 was used in this study. Agrobacterium rhizogenes-transformed root organ cultures (ROC) expressing the 35S:NupYC2.1 construct  were obtained according to . NupYC2.1 is a cameleon probe that exploits Förster resonance energy transfer (FRET) to highlight variations in the level of Ca2+ concentration in the nucleoplasm. Transformed roots with a good level of nuclear fluorescence were selected three weeks after inoculation, decontaminated and grown as ROCs on minimal (M) medium at 25°C in the dark . Since Ca2+ spiking responses to AM fungal signals are equivalent in composite plants and ROCs [10, 13], we have only used the latter in this work. As summarized in Table 1, ROCs were grown in the presence of 35 μm KH2PO4 (MLP), a standard condition for in vitro AM development , or 3.5 mM (MHP), a concentration that does not affect root growth but inhibits root colonization by AM fungi . M. truncatula ROCs were grown in vertically oriented Petri dishes to favour the development of a regular fishbone-shaped root system . The AM fungus used in this study was Gigaspora margarita isolate BEG 34 (International Bank for the Glomeromycota, University of Kent, UK).
The spores were harvested from cultures of Trifolium repens, as described in . In brief, spores were collected by repeated sieving under running water and stored at 4°C for one week to increase the germination rate.
To record nuclear Ca2+ levels we applied the protocol described in . Segments of primary roots of M. truncatula carrying one or two young laterals (2–4 cm in length) were placed in a 2-mm-thick microchamber on a microscope slide containing sterile distilled water. The water in the microchamber was then rapidly replaced by 100 μl of the treatment solution before starting image acquisition.
In our first case study, MHP and MLP roots were treated with fungal exudates. To produce the exudate (E), batches of 100 surface-sterilized G. margarita spores were placed in either 1 ml of sterile distilled water (ELP) or 3,5 mM KH2PO4 solution (EHP; see Table 1) and incubated for 7 d at 30°C in the dark to induce germination (germination rate > 90%). The germination medium was recovered by pipetting, concentrated 10-fold using a Lio5P lyophilizer (Cinquepascal, Milan, Italy) and stored at −20°C.
For our second case study, roots were treated with water solutions containing either 10-8 M CO4, 10-8 M Nod factor from S. meliloti or ELP. Nod factor treatments were performed on transformed roots from composite plants [13, 47], since Nod factor responsiveness is absent in excised ROCs . All records of nuclear Ca2+ spiking in M. truncatula root hairs in response to Nod factor were kindly provided by Mireille Chabaud and David Barker (LIPM, Toulouse, France). In both case studies, control treatments were included where the treatment solutions were replaced with sterile distilled water.
Confocal microscopy and measurement of changes in nuclear Ca2+levels
Confocal microscopy was used for all the FRET experiments. FRET-based detection and plotting of relative changes in nuclear Ca2+ levels were performed according to , by measuring the ratio of yellow fluorescent protein (YFP) to cyan fluorescent protein (CFP) signal intensity over time [50, 51]. A 40× water-immersion objective was used and the pinhole was set to 6 Airy Units so that the thickness of the optical sectioning would embrace the average diameter of epidermal cell nuclei. 512 by 512 pixel frames were collected every 5 s for 30 minutes. Since glomeromycota and rhizobia target different cell types, analyses involving fungal signals were done on ROC atrichoblasts, while Nod factors responses were recorded in root hairs from composite plants. The exact number of biological replicates for each condition is reported in Additional file 6.
The Leica Confocal Software was used to output the original values of YFP and CFP fluorescence intensity into text files (Additional file 7). Each file contained the values derived from five cells. Text files were then imported into a Microsoft Excel 2008 model file, where the ratio of YFP/CFP values was calculated for each time series (Additional file 8). The columns containing the ratiometric series were then exported to a new text file, to be used as an input file for the CaSA software. The CaSA software requires that both input files containing data series and control series are formatted as shown in Additional file 9. In brief, all time series must have the same length (in our case 360 time points); the series corresponding to each cell must be presented in tab-separated columns, where the first row contains the sample ID codes (e.g. “CELL_1”). As required for the Octave language, the first character of the text file must be #.
The analysis of time series displaying a spiking pattern (in our case a repetition of peaks in nuclear Ca2+ concentration) must address the following issues:
removing fluctuations due to measurement errors, to better discriminate signals from background noise;
identification of the instant in the timeline when each peak event occurs, in order to allow automated operations such as counting the number of peaks and measuring peak-to-peak intervals;
characterization of trends, to identify conserved traits in a population of records; automatically producing statistical analyses of the data obtained.
The first requirement implies the application of a smoothing filter to the original data. Such a smoothing process must level out background oscillations due to instrumental errors, so that only the major signal oscillations (peaks) stand out and can subsequently be identified and analyzed by the software. To this aim, the low-pass Savitzky-Golay filter family  is of common use and has demonstrated particularly efficient on noisy data sets . In brief, rather than having their properties defined in the Fourier domain and then converted to the time domain, the Savitzky-Golay filters are derived directly from a particular formulation of the data smoothing problem in the time domain. This class of filters also provides an estimate of point to point derivatives along a curve, allowing the assessment of analytical properties.
The operating assumptions for the use of the method underlying these filters are: 1) the data to be processed are recorded at fixed and uniform intervals; 2) the phenomenon analyzed can be represented by a continuous curve.
The method takes into account a sliding window of measurements embracing the point f0 on which the filtering is applied, and two finite sequences of points of length n respectively preceding and following f0.
A polynomial curve of degree k is calculated on this window of points to approximate the trend of the quantity with a continuous curve. Since in general 2 n + 1 > k +1 the method of least squares can be used. The filtered value of f0 is given by the value of the fitting curve in position 0 in relation to preceding values f− n, f− n + 1, · ··, f− 1 and following values f1, · ··, fn − 1fn.
From the polynomial approximation, an estimate of the first and second derivatives in position i can also determine the presence of a peak considering the conditions:
f'i = 0
f''i < 0
Other conditions, such as threshold exceeding and asymmetries of trends can then be integrated with the above conditions to form precise patterns of merit characterizing the peaks. For example, the expression
|·| = 1 if true, 0 if false
expresses the presence of a peak once the δ ξ γ parameters (which characterize value, trend and symmetry respectively) are fixed. In the practice, the above function has been implemented with a continuous function peak c (i) using the product of absolute values instead of the logic function. A threshold τ has then been introduced to discriminate peaks from background oscillations: peak (i) = 1 if and only if peak c (i) > τ.
By analyzing the trend in the neighborhood of f0, it is also possible to estimate the peak event start (t s ) and end (t e ), by scrolling the series of the derivatives of each peak from left to right. From this temporal information and values detected, the intensity of a phenomenon can be evaluated with the integration:
which gives the area under the curve f (t) from t s to t e This value allows the comparison of peaks from different time series. In our case, considering a series from the control data set, the above formulation allows the determination of a value of intensity below which a peak is not to be considered significant and must therefore be excluded from computation in the treatment data set.
Kudla J, Batistic O, Hashimoto K: Calcium signals: the lead currency of plant information processing. Plant Cell. 2010, 22: 541-563. 10.1105/tpc.109.072686.
Dodd AN, Kudla J, Sanders D: The language of calcium signaling. Annu Rev Plant Physiol Plant Mol Biol. 2010, 61: 593-620.
McAinsh MR, Pittman JK: Shaping the calcium signature. New Phytol. 2009, 181: 275-294. 10.1111/j.1469-8137.2008.02682.x.
McAinsh MR, Webb AAR, Taylor JE, Hetherington AM: Stimulus-induced oscillations in guard cell cytosolic free calcium. Plant Cell. 1995, 7: 1207-1219.
Allen GJ, Chu SP, Harrington CL, Schumacher K, Hoffmann T, Tang YY, Grill E, Schroeder JI: A defined range of guard cell calcium oscillation parameters encodes stomatal movements. Nature. 2001, 411: 1053-1057. 10.1038/35082575.
Hepler PK, Lovy-Wheeler A, McKenna ST, Kunkel JG: Ions and pollen tube growth. Plant Cell Mon. 2006, 3: 47-69. 10.1007/7089_043.
Whalley HJ, Knight MR: Calcium signatures are decoded by plants to give specific gene responses. New Phytol. 2013, 197: 690-693. 10.1111/nph.12087.
Meyer T, Stryer L: Calcium spiking. Ann Rev Biophys Bio. 1991, 20: 153-174. 10.1146/annurev.bb.20.060191.001101.
Kosuta S, Hazledine S, Sun J, Miwa H, Morris RJ, Downie JA, Oldroyd GED: Differential and chaotic calcium signatures in the symbiosis signaling pathway of legumes. Proc Natl Acad Sci USA. 2008, 105: 9823-9828. 10.1073/pnas.0803499105.
Chabaud M, Genre A, Sieberer BJ, Faccio A, Fournier J, Novero M, Barker DJ, Bonfante P: Arbuscular mycorrhizal hyphopodia and germinated spore exudates trigger Ca2+ spiking in the legume and nonlegume root epidermis. New Phytol. 2011, 189: 347-355. 10.1111/j.1469-8137.2010.03464.x.
McAinsh MR, Hetherington AM: Encoding specificity in Ca2+ signaling systems. Tr Plant Sci. 1998, 3: 32-36. 10.1016/S1360-1385(97)01150-3.
Hazledine S, Sun J, Wysham D, Downie JA, Oldroyd GED, Morris RJ: Nonlinear time series analysis of nodulation factor induced calcium oscillations: evidence for deterministic chaos?. PLoS One. 2009, 4: 6637-6637. 10.1371/journal.pone.0006637.
Genre A, Chabaud M, Balzergue C, Puech-Pagès V, Novero M, Rey T, Fournier J, Rochange S, Bécard G, Bonfante P, Barker DG: Short-chain chitin oligomers from arbuscular mycorrhizal fungi trigger nuclear Ca2+ spiking in Medicago truncatula roots and their production is enhanced by strigolactone. New Phytol. 2013, 198: 190-202. 10.1111/nph.12146.
Quian Quiroga R, Nadasdy Z, Ben-Shaul Y: Unsupervised spike sorting with wavelets and superparamagnetic clustering. Neural Comput. 2004, 16: 1661-1687. 10.1162/089976604774201631.
Wolf MT, Burdick JW: A Bayesian clustering method for tracking neural signals over successive Intervals. Trans Biomed Eng. 2009, 56: 2649-2659.
Sciacca E, Spinella S, Genre A, Calcagno C: Analysis of calcium spiking in plant root epidermis through cwc modeling. Electr Notes Theor Comp Sci. 2011, 277: 65-76.
GNU Octave. http://www.gnu.org/software/octave/.
Gage DJ: Infection and invasion of roots by symbiotic, nitrogen- fixing rhizobia during nodulation of temperate legumes. Microbiol Mol Biol Rev. 2004, 68: 280-300. 10.1128/MMBR.68.2.280-300.2004.
Dénarié J, Cullimore J: Lipo-oligosaccharide nodulation factors: a minireview new class of signaling molecules mediating recognition and morphogenesis. Cell. 1993, 74: 951-954. 10.1016/0092-8674(93)90717-5.
Maillet F, Poinsot V, André O, et al: Fungal lipochitooligosaccharide symbiotic signals in arbuscular mycorrhiza. Nature. 2011, 469: 58-64. 10.1038/nature09622.
Singh S, Parniske M: Activation of calcium- and calmodulin-dependent protein kinase (CCaMK), the central regulator of plant root endosymbiosis. Curr Op Plant Biol. 2012, 15: 444-453. 10.1016/j.pbi.2012.04.002.
Ehrhardt DW, Wais R, Long SR: Calcium spiking in plant root hairs responding to Rhizobium nodulation signals. Cell. 1996, 85: 673-681. 10.1016/S0092-8674(00)81234-9.
Sieberer BJ, Chabaud M, Timmers AC, Monin A, Fournier J, Barker DG: A nuclear-targeted cameleon demonstrates intranuclear Ca2+ spiking in Medicago truncatula root hairs in response to rhizobial nodulation factors. Plant Physiol. 2009, 151: 1197-1206. 10.1104/pp.109.142851.
Smith SE, Smith FA: Roles of arbuscular mycorrhizas in plant nutrition and growth: new paradigms from cellular to ecosystems scales. Annu Rev Plant Physiol Plant Mol Biol. 2011, 63: 227-250.
Smith SE, Read D: Mycorrhizal Symbiosis 3rd edition. London: Academic Press; 2008,
Balzergue C, Puech-Pagès V, Bécard G, Rochange SF: The regulation of arbuscular mycorrhizal symbiosis by phosphate in pea involves early and systemic signalling events. J Exp Bot. 2011, 62: 1049-1060. 10.1093/jxb/erq335.
Volpe V, Dell’aglio E, Bonfante P: The Lotus japonicus MAMI gene links root development, arbuscular mycorrhizal symbiosis and phosphate availability. Plant Sign Behav. 2013, doi:org/10.4161/psb.23414
Gu M, Chen A, Dai X, Liu W, Xu G: How does phosphate status influence the development of the arbuscular mycorrhizal symbiosis?. Plant Sign Behav. 2011, 6: 1300-1304. 10.4161/psb.6.9.16365.
Czarnecki O, Yang J, Weston DJ, Tuskan GA, Chen JG: A dual role of strigolactones in phosphate acquisition and utilization in plants. Int J Mol Sci. 2013, 14: 7681-7701. 10.3390/ijms14047681.
Oldroyd GED, Downie JA: Nuclear calcium changes at the core of symbiosis signalling. Curr Op Plant Biol. 2006, 9: 351-357. 10.1016/j.pbi.2006.05.003.
Parniske M: Arbuscular mycorrhiza: the mother of plant root endosymbiosis. Nature. 2008, 8: 763-775.
Shimoda Y, Han L, Yamazaki T, Suzuki R, Hayashi M, Imaizumi-Anraku H: Rhizobial and fungal symbioses show different requirements for calmodulin binding to calcium calmodulin-dependent protein kinase in Lotus japonicus. Plant Cell. 2012, 24: 304-321. 10.1105/tpc.111.092197.
Kiers ET, Duhamel M, Beesetty Y, et al: Reciprocal rewards stabilize cooperation in the mycorrhizal symbiosis. Science. 2011, 333: 880-882. 10.1126/science.1208473.
Breuillin F, Schramm J, Hajirezaei M, et al: Phosphate systemically inhibits development of arbuscular mycorrhiza in Petunia hybrida and represses genes involved in mycorrhizal functioning. Plant J. 2010, 64: 1002-1017. 10.1111/j.1365-313X.2010.04385.x.
Fiorilli V, Lanfranco L, Bonfante P: The expression of GintPT, the phosphate transporter of Rhizophagus irregularis, depends on the symbiotic status and phosphate availability. Planta. 2013, 237: 1267-1277. 10.1007/s00425-013-1842-z.
Sanders FE: The Effect of Foliar-Applied Phosphate on the Mycorrhizal Infection of Onion Roots. Endomycorrhizas. Edited by: Sanders FE, Mosse B, Tinker PB. London: Academic Press; 1975, 261-276.
Bécard G, Fortin JA: Early events of vesicular-arbuscular mycorrhiza formation on Ri T-DNA transformed roots. New Phytol. 1988, 108: 211-218. 10.1111/j.1469-8137.1988.tb03698.x.
Chabaud M, de Carvalho-Niebel F, Barker DG: Efficient transfor- mation of Medicago truncatula cv. Jemalong using the hypervirulent Agrobacterium tumefaciens strain AGL1. Plant Cell Rep. 2003, 22: 46-51. 10.1007/s00299-003-0649-y.
Balzergue C, Chabaud M, Barker DG, Bécard G, Rochange SF: High phosphate reduces host ability to develop arbuscular mycorrhizal symbiosis without affecting root calcium spiking responses to the fungus. Front Plant Sci. 2013, 4: 426.
Akiyama K, Matsuzaki K, Hayashi H: Plant sesquiterpenes induce hyphal branching in arbuscular mycorrhizal fungi. Nature. 2005, 435: 824-827. 10.1038/nature03608.
Kistner C, Winzer T, Pitzschke A, Mulder L, Sato S, Kaneko T, Tabata S, Sandal N, Stougaard J, Webb KJ, Szczyglowski K, Parniske M: Seven Lotus japonicus genes required for transcriptional reprogramming of the root during fungal and bacterial symbiosis. Plant Cell. 2005, 17: 2217-2229. 10.1105/tpc.105.032714.
Ben Amor B, Shaw SL, Oldroyd GE, Maillet F, Penmetsa RV, Cook D, Long SR, Dénarié J, Gough C: The NFP locus of Medicago truncatula controls an early step of Nod factor signal transduction upstream of a rapid calcium flux and root hair deformation. Plant J. 2003, 34: 495-506. 10.1046/j.1365-313X.2003.01743.x.
Gobbato E, Marsh JF, Vernié T, Wang E, Maillet F, Kim J, Miller JB, Sun J, Bano SA, Ratet P: A GRAS-type transcription factor with a specific function in mycorrhizal signalling. Curr Biol. 2012, 22: 2236-2241. 10.1016/j.cub.2012.09.044.
Liu W, Kohlen W, Lillo A, Op den Camp R, Ivanov S, Hartog M, Limpens E, Jamil M, Smaczniak C, Kaufmann K: Strigolactone biosynthesis in Medicago truncatula and rice requires the symbiotic GRAS-type transcription factors NSP1 and NSP2. Plant Cell. 2011, 23: 3853-3865. 10.1105/tpc.111.089771.
Cover TM, Thomas JA: Elements of Information Theory. Wiley-interscience: New York, USA; 1991,
Janicek R, Hotka M, Zahradníková A, Zahradníková A, Zahradník I: Quantitative analysis of calcium spikes in noisy fluorescent background. PLoS ONE. 2013, 8: 643-694.
Boisson-Dernier A, Chabaud M, Garcia F, Bécard G, Rosenberg C, Barker DG: Agrobacterium rhizogenes-transformed roots of Medicago truncatula for the study of nitrogen-fixing and endomycorrhizal symbiotic associations. Mol Plant-Microbe Interact. 2001, 14: 695-700. 10.1094/MPMI.2001.14.6.695.
Chabaud M, Venard C, Defaux-Petras A, Bécard G, Barker DG: Targeted inoculation of Medicago truncatula in vitro root cultures reveals MtENOD11 expression during early stages of infection by arbuscular mycorrhizal fungi. New Phytol. 2002, 156: 265-273. 10.1046/j.1469-8137.2002.00508.x.
Ramos J, Bisseling T: A method for the isolation of root hairs from the model legume Medicago truncatula. J Exp Bot. 2003, 54: 2245-2250. 10.1093/jxb/erg245.
Miyawaki A, Llopis J, Heim R, McCaffery JM, Adams JA, Ikura M, Tsien RY: Fluorescent indicators for Ca2+ based on green fluorescent proteins and calmodulin. Nature. 1997, 388: 882-888. 10.1038/42264.
Miyawaki A, Griesbeck O, Heim O, Tsien R: Dynamic and quantitative Ca2+ measurements using improved cameleons. Proc Natl Acad Sci USA. 1999, 9: 2135-2140.
Savitzky A, Golay MJE: Smoothing and differentiation of data by simplified least squares procedures. Anal Chem. 1964, 36: 1627-1639. 10.1021/ac60214a047.
Chen J, Jönsson P, Tamura M, Gu Z, Matsushita B, Eklundh L: A simple method for reconstructing a high-quality NDVI time-series data set based on the Savitzky-Golay filter. Remote Sens Environ. 2004, 91: 332-344. 10.1016/j.rse.2004.03.014.
We thank the anonymous reviewers for their constructive critics. We are grateful to Mireille Chabaud, Björn Sieberer and David Barker (LIPM, Toulouse, France) for sharing their records of Nod factor-induced Ca2+ spiking, and to David Barker for critical revision of the text. We also thank Mara Novero for her assistance in the production of the fungal exudate. Financial support for this research was granted to PB by Regione Piemonte BioBITs Project (Converging Technologies 2007, area: Biotechnology-ICT) and by the PRIN Project 2011 Pro-Root.
The authors declare that they have no competing interests.
GR carried out all the experimental work, contributed to developing the automated analysis and the CaSA software; SS and ES developed the automated analysis and the CaSA software; PB participated in conceiving and coordinating the study; AG participated in designing and coordinating the study and drafting the manuscript. All authors contributed to the final version of the manuscript, which they have read and approved.
Giulia Russo, Salvatore Spinella contributed equally to this work.
Electronic supplementary material
Additional file 1: Pipeline of CaSA usage. The flow chart schematizes file preparation and data flow through the CaSA software. The steps leading to input data file production are representative of our experiments, but CaSA only requires the input data files, independently of the upstream experimental setup. (PDF 71 KB)
Additional file 3: The figure shows how five different settings of the CaSA parameters (changes are highlighted in red) impact on peak identification (green triangles) in the resulting output plots. The optimal conditions for our experimental data set (which are proposed as a default by the CaSA software) are presented in the central box. For example, the width of the sliding measurement window (w) is set to 9 sampling points: with our sampling interval of 5 seconds, the window covers 45 (9 x 5) seconds, which is comparable with the average peak duration in our records. Longer or shorter windows can be chosen to best fit data sets with longer or shorter peak durations, respectively. Two examples are provided for each setting: one plot where the baseline shows a steep slope and peaks are more irregular (top), and one where peaks are very evident (bottom). While peaks are always identified correctly in the bottom plots, the fine tuning of CaSA parameters w, τ and q is crucial to optimize peak recognition in the top plots. For example, after increasing the w value to 13 (bottom box), changing τ from 1 to 2 or reducing q to 0.7, one or more peaks are left unmarked (asterisks). (PDF 363 KB)
Additional file 4: CaSA output file. This is an example of the output text file produced by CaSA. The file displays the results of automated analysis for each sample. Sample IDs are reported in the first row, followed by statistical analyses on waiting times, peak numbers, percentage of responsive samples in the dataset (displaying more than 2 peaks during the recorded period). (TXT 3 KB)
Additional file 5: CaSA software. This compressed (zip) file contains the CaSA software and associated files. All files should be downloaded to the same directory for the software to work. The software can be run from the terminal using the command line ./CaSA.m or by double-clicking on the CaSA.m icon from the file manager. In this case the file must be previously set as executable in the file properties. Input files should also be placed in the same directory as CaSA. (ZIP 150 KB)
Additional file 6: Table summarizing the number of biological replicates used to generate the data for this research. Since glomeromycota and rhizobia target different cell types, fungal signals effects were analysed in atrichoblasts from root organ cultures (ROC) and Nod factors (NF) effects in root hairs from composite plants. (PDF 49 KB)
Additional file 7: This text file is the standard output obtained from the LCS software by exporting 30 min recordings of the fluorescence intensity (stack profile function) in five epidermal nuclei from M. truncatula roots treated with the fungal exudate. (TXT 48 KB)
Additional file 8: This Microsoft Excel file was used to calculate the FRET values (sheet 1) from the original data exported by the Leica Confocal Software (sheet 2, Additional file 7). (XLSX 160 KB)
Additional file 9: This is the input file type required by the CaSA software. Two files are requested, one for the treated samples and one for the controls, where only background noise is recorded. Each column in the file contains the sample ID number and the FRET values as obtained by the Microsoft Excel spreadsheet (Additional file 8). (TXT 26 KB)
Authors’ original submitted files for images
About this article
Cite this article
Russo, G., Spinella, S., Sciacca, E. et al. Automated analysis of calcium spiking profiles with CaSA software: two case studies from root-microbe symbioses. BMC Plant Biol 13, 224 (2013). https://doi.org/10.1186/1471-2229-13-224
- Arbuscular mycorrhiza
- Calcium signaling
- Medicago truncatula
- Nitrogen fixation
- Plant-microbe interactions
- Automated data analysis