Visualizing and quantifying ABA-triggered [Ca2+] dynamics at single-cell resolution
We postulated that an ultrasensitive GECI might reliably detect ABA-induced Ca2+ signals in various plant cell types. To support this hypothesis, we performed time-lapse recordings of ABA responses in the roots of transgenic GCaMP6s plants10,27. ABA triggered rapid and specific Ca2+ waves with dual peaks at 100–200 s in the root tip never previously observed, which were distinct from ABA-activated Ca2+ oscillations in guard cells23,25,26 (Extended Data Fig. 1). To simplify and improve the quantification of GCaMP6s and maintain its Ca2+ sensitivity and high dynamic range, we created a synthetic gene that expressed a single protein with two fluorescent proteins—the green fluorescent GCaMP6s and red fluorescent dTomato—fused by a short SSGS linker (Fig. 1a,b). In this simple Ca2+ ratiometric sensor (CRS), the C-terminal dTomato served as a stable internal control insensitive to Ca2+ and did not interfere with the Ca2+ binding of the N-terminal GCaMP6s (Fig. 1c). While the Ca2+-induced green fluorescence of CRS was inhibited by EGTA chelation of Ca2+, the red fluorescence of CRS was unaffected (Fig. 1c). Thus, the ratio of Ca2+-sensitive green fluorescence to Ca2+-insensitive red fluorescence of CRS with excitation at 488 nm and 554 nm, respectively, indicates the relative Ca2+ signal changes independent of protein expression levels. The Kd of CRS for Ca2+ was 165 nM by in vitro binding assays and was comparable to that of GCaMP6s (Fig. 1d)27,30. As a positive control, we first measured [Ca2+] dynamics after nitrate stimulation in the mesophyll protoplasts transiently expressing CRS. The ratio of green to red fluorescence reported a similar dynamic Ca2+ profile as previously recorded using GCaMP6s (Fig. 1e)10. Remarkably, we also detected reliable ABA-induced [Ca2+] dynamics with a distinct cytosolic pattern and lower amplitude in the transfected mesophyll protoplasts, but not in the mock control without ABA (Fig. 1f,g)10.
To enable analyses of ABA-induced Ca2+ responses in different plant organs and cell types, we generated CRS transgenic plants. Because ABA-induced Ca2+ oscillations were previously observed in guard cells, we first validated the use of CRS in transgenic plants by analysing ABA-induced Ca2+ response in guard cells. Ca2+ oscillations in the cytosol of guard cells were observed only after ABA stimulation (Fig. 2a and Extended Data Fig. 2). We then monitored specific ABA-triggered Ca2+ responses in different organs, tissues, and cell types, not previously feasible. After ABA treatment, the stimulated mesophyll cells initiated a slower Ca2+ wave featuring a different shape and kinetic compared with the faster and more sustained Ca2+ wave in the root tip or root differentiation zone, or the rapid and sharper Ca2+ wave in the lateral root primordium (Fig. 2b-e). CRS-based live imaging in plants revealed for the first time that ABA could trigger diverse Ca2+ signatures distinct from Ca2+ oscillations in guard cells and might activate different cellular Ca2+ sensors to relay ABA signalling responses in various organs and cell types, as implicated by ABA reporters and next-generation ABACUS biosensors1–7,11−23,31,32.
Dynamic spatiotemporal patterns of ABA-triggered Ca2+ signatures
Nitrate and ABA appeared to activate different subcellular patterns of [Ca2+] changes and dynamics in mesophyll protoplasts (Fig. 1e-g)10. Unlike GCaMP6s detecting nanomolar [Ca2+] changes in both the cytosol and the nucleus10, CRS was restricted in detecting cytosolic [Ca2+]cyt. To overcome the subcellular limit of CRS, we created new versions of CRS localized at the plasma membrane or specific organelles to further characterize the high spatiotemporal resolution of Ca2+ signatures in different subcellular compartments (Fig. 3a). Importantly, transgenic plants expressing CRS fusion proteins showed no obvious growth defect phenotype (Extended Data Fig. 3). We analyzed ABA-triggered Ca2+ signatures in transgenic plants stably expressing the CRS fusion protein variants targeted to the plasma membrane (CRS-PM), cytosol (CRS), or nucleus (CRS-NLS). As the root tip exhibited prominent ABA-responsive reporter gene expression31 and Ca2+ dynamics (Extended Data Fig. 1 and Fig. 2c), we monitored time-resolved ABA-triggered [Ca2+] change patterns in subcellular compartments in the uniform-shaped epidermal cells in the root meristematic zone (Fig. 2c, Fig. 3b-f and Extended Data Fig. 4, 5). When treated with ABA, cells expressing CRS-PM displayed the most rapid elevation of Ca2+ signals (< 50 s), which appeared to represent the first phase of two-phase cytosolic Ca2+ increases in timing and amplitude (Extended Data Fig. 1). The second phase of the elevated cytosolic Ca2+ signal was sustained for more than 5 min when significant decline occurred near the plasma membrane after 2 min. Unexpectedly, slower (approximately 100 s) but striking Ca2+ oscillation was observed with CRS-NLS (Fig. 3b-d and Extended Data Figs. 4 and 5). Imaging cytosolic [Ca2+] changes in single cells expressing CRS variants suggests that the initial Ca2+ signal increase due to ABA treatment might be originated at the discrete nanodomains of the plasma membrane, and then spread throughout the plasma membrane before migrating to the cytosol and into the nucleus (Fig. 1f and Fig. 3b-d and Extended Data Figs. 4 and 5). The findings are consistent with the detection of ABA receptors and PP2C/ABI1-CPK21 signalling complexes on the plasma membrane nanodomains33–35
To directly visualize high spatiotemporal ABA-Ca2+ signalling processes at the single-cell resolution, we performed time-lapse imaging analysis along the plasma membrane of single CRS-expressing mesophyll cells exhibiting slower ABA-activated [Ca2+] dynamics (Fig. 1f, g and Fig. 2b). We found that ABA-triggered Ca2+ signal elevations could start from discrete sites and were not evenly distributed along the plasma membrane before migrating into the cytosol (Extended Data Fig. 6a). In the root tip with faster ABA-Ca2+ signalling (Extended Data Fig. 1 and Fig. 2c), the CRS seedlings were treated with extracellular Ca2+ chelator, EGTA or BAPTA, or a plasma membrane Ca2+ channels blocker, GdCl3, the ABA-induced cytosolic Ca2+ waves were abolished (Extended Data Fig. 6b). To further demonstrate the ABA-triggered Ca2+ signal was first initiated near the plasma membrane and distinct from the internal nitrate-Ca2+ signalling in GCaMP6s expressing mesophyll protoplasts10, we generated a transgenic line harbouring dual CRS-PM and CRS-NLS sensors (CRS-PM/CRS-NLS). When CRS-PM/CRS-NLS seedlings were treated with ABA, Ca2+ signals first appeared near the plasma membrane at 10–20 s and then in the nucleus at 50 s in the single epidermal cells of the primary root meristem (Fig. 3e). We further conducted a higher-resolution kymography analysis of CRS-PM in response to ABA by dividing the perimeter of a single root epidermal cell into 440 adjacent regions (Fig. 3f). Combining live imaging analyses in single leaf and root cells with CRS variants provided strong evidence that elevated Ca2+ signals originated along the plasma membrane by ABA was not homogeneously distributed. CRS enabled the ultrasensitive detection of distinct extracellular and intracellular Ca2+ sources and dynamics in ABA and nitrate signalling, respectively10.
Calibration of CRS for quantifying intracellular [Ca2+]
To quantify the precise in vivo nanomolar cytosolic [Ca2+]cyt in response to ABA, we recorded [Ca2+] changes in the root tip of CRS transgenic plant in response to ABA using two references generated by extracellular ATP (eATP)-activated micromolar increase in [Ca2+] and Ca2+ chelators to deplete free Ca2+ in the root tip cells (Fig. 4a)24–26. To calibrate and calculate the Ca2+ concentration of CRS in the root tip cells, we first treated CRS seedlings with eATP to saturate cellular Ca2+ and obtain the maximum fluorescence intensity and ratio. We then incubated the same plants in EGTA and BAPTA-AM to acquire the minimum fluorescence intensity and ratio. Since the in vitro Kd of CRS for Ca2+ is 165 nM, the maximum and minimum fluorescence intensity ratios were set at 330 nM and 0 nM, respectively, to calculate the intracellular [Ca2+] activated by ABA in the root tip (Fig. 4a,b). ABA increased [Ca2+]cyt to ~ 164.3 nM in the root tip cells expressing CRS.
ABA-triggered nanomolar [Ca2+] wave via ABA receptors and SNRK2.2/2.3/2.6
Despite decades of intensive research11–26, how ABA activates [Ca2+] increase remained unknown in most cell types except guard cells14,23. Arabidopsis PYR1 and PYL1/2/4/5/8 are functionally redundant ABA receptors that regulate ABA inhibition of primary root formation, while SNRK2.2/2.3/2.6 are key activators of ABA signalling2,3,34,35. To determine whether ABA-triggered [Ca2+] increase requires ABA receptors and SNRK2.2/2.3/2.6, we generated stable transgenic lines expressing CRS in the sextuple pyrpyl112458 and triple snrk2.2,2.3,2.6 mutants34,35. The ABA-induced [Ca2+] increase was eliminated in the pyrpyl112458-CRS and snrk2.2,2.3,2.6-CRS mutants, (Fig. 4c,d). To show that the abolished [Ca2+] increase was not due to CRS failure, we treated the same seedlings with eATP after recording the ABA responses. Striking [Ca2+] rise was induced by eATP in pyrpyl112458-CRS and snrk2.2,2.3,2.6-CRS, indicating that the ABA-triggered [Ca2+] increase is specifically regulated downstream of ABA receptors and SNRK2.2/2.3/2.6 (Fig. 4e,f).
The delayed root cap differentiation cycle by ABA-Ca-CPK signalling
CRS detected conspicuous ABA-triggered Ca2+ signalling in the root tip (Extended Data Fig. 1, Fig. 2c, and Fig. 4) where the root cap is located and plays essential functions in stem cell niche protection, gravitropism, thigmotropism, water/nutrient sensing and uptake, root system architecture, as well as protection from biotic and abiotic stress36–40. To explore ABA-Ca2+ signalling functions and regulatory mechanisms in the root cap, we examined the effect of ABA on the unique root cap cycle encompassing the differentiation programs of the columella and lateral root cap cells from stem cell regulation, differentiation, and maturation, to detachment (Fig. 5a) 36–40. We first defined the periodicity of the root cap cycle by live imaging of the root cap morphology changes at 3-h intervals for 48 h starting with 5-day-old seedlings using an ECHO Revolve microscope (ECHO) (Fig. 5b). The average cycle period of root cap initiation and sloughing is approximately 30.4 ± 0.5 h (Fig. 5b). ABA treatment delayed the root cap cycle to approximately 34 ± 0.4 h. To further confirm the involvement of ABA in regulating the root cap cycle, we examined ABA regulation in the sextuple ABA receptor pyrpyl112458 mutant (Fig. 5b). The pyrpyl112458 mutant was insensitive to ABA treatment in delaying root cap cycle. Since ABA triggered a maximum [Ca2+]cyt of 164.3 nM in the root tip of the CRS seedlings (Fig. 4a), previous studies suggested that CPK10/30/32 with a higher Ca2+ sensitivity than other CPK subgroups may participate in ABA signalling10,15,18–20,41. To test if CPK10/30/32 may serve as the Ca2+ sensors to decode ABA-Ca2+ signalling in modulating the root cap cycle, we assessed the ABA response in the chemical-inducible triple icpk mutant (Fig. 5c)10. In the icpk seedling without 3MBiP
(1-isopropyl-3-(3-methylbenzyl)-1H-pyrazolo[3,4-d]pyrimidin-4-amine) induction10, the average period of root cap cycle was similar to that in WT at approximately 31.5 ± 0.7 h. ABA treatment delayed the root cap cycle to 37.5 ± 0.7 h, but was reduced to approximately 34.3 ± 0.4 h in icpk (Fig. 5c).
To elucidate the molecular function of ABA in regulating the root cap cycle, we investigated the expression of key marker genes involved in regulating the differentiation processes. Transcription factor FEZ controls the columella stem cells and the epidermal/lateral root cap stem cells responsible for the initiation of the root cap. Transcription factor SOMBERO (SMB) represses FEZ to promote root cap differentiation, whereas BEARSKIN 2 (BRN2) is a central transcription factor for the maturation and detachment processes by activating target genes encoding enzymes involved in cell wall degradation, e.g., CELLULASE 5 (CEL5), PECTINESTERASE (PE11), and XYLOGLUCAN ENDOTRANSGLUCOSYLASE 5 (XTH5), and programmed cell death, including RIBONUCLEASE 3 (RNS3), and METACASPASE 9 (MC9) (Fig. 5a)36–40. RT-qPCR analysis was performed on RNA isolated from the root tip of WT and pyrpyl112458 mutants without or with ABA. The results revealed that ABA activated a typical marker gene KIN114 but repressed FEZ, SMB, BRN2, CEL5, XTH5, PE11, RNS3,and MC9 expression (Fig. 5d and Extended Data Fig. 7a). The ABA repression of these genes was abolished in the pyrpyl112458 mutant (Fig. 5d, Extended Data Fig. 7a). Similarly, ABA repression of SMB, BRN2, CLE5, PE11, RNS3, and MC9 was significantly reduced in the icpk mutant. The ABA repression of FEZ, CEL5, and XTH5 required ABA receptors but was CPK10/30/32 independent (Fig. 5e, Extended Data Fig. 7b). These data further support the hypothesis that ABA-Ca2+-CPK signalling is involved in regulating the root cap cycle by suppressing the expression of key genes in the developmental program of the entire root cap cycle.
Altering the root cap cycle plasticity by unique salt-ABA-Ca2+ signalling
Current knowledge supports that high salt stress triggers a transient and strong [Ca2+] increase and the classical salt-overly-sensitive (SOS) signalling pathway mediated by the Ca2+ sensor SOS3 (CBL4), the protein kinase SOS2 (CIPK24) and the Na+/H+ antiporter SOS1 in roots28,42,43. Little was known about salt stress effects on root cap cells. In CRS transgenic plants, we confirmed that 200 mM NaCl activated a sharp 50-s cytosolic Ca2+ spike in the root elongation zone (Fig. 6a)28. Unexpectedly, in contrast to the prevailing NaCl-activated micromolar [Ca2+] spike, salt stress activated sustained nanomolar [Ca2+] wave in root cap cells (Fig. 6b) and significantly delayed the root cap cycle (Fig. 6c). This novel salt-stimulated [Ca2+] increase was predominantly dependent on ABA signalling as it was significantly attenuated in the pyrpyl112458-CRS and snrk2.2,2.3,2.6-CRS mutants (Fig. 6b). Our result suggests the involvement of two salt signalling pathways, ABA-independent and ABA-dependent, and the later orchestrates the salt stress-induced changes in [Ca2+] in the root cap. To assess the impact of salt stress on the regulation of the root cap cycle, we established a long-term observation method to investigate salt stress-induced phenotypes. Five-day-old seedlings were treated with 50 mM NaCl and the differentiation, maturation, and detachment of the root cap were monitored under a microscope every 3 h for 60 h. The average period of the root cap cycle from initiation to shedding was delayed from 31.4 ± 0.7 h (in mock treatment) to approximately 41.1 ± 0.6 h with NaCl treatment (Fig. 6c). The delayed root cap cycle showed insensitivity in the pyrpyl112458 and snrk2.2,2.3,2.6 mutants with 50 mM NaCl treatment (Fig. 6c). Additionally, RT-qPCR analysis was performed with RNA isolated from the root tip of WT seedlings with or without 200 mM NaCl treatment. Salt stress activated the expression of KIN1, a salt-stress- and ABA-responsive gene21, while repressed FEZ, SMB, BRN2, PE11, RNS3, MC9, XTH5, and CEL5 expression (Fig. 6d, Extended Data Fig. 8) similar to the ABA repression (Fig. 5d,e and Extended Data Fig. 7). The data presented here connect salt stress to the previously unknown endogenous ABA and Ca2+ signalling in controlling the root cap cycle (Fig. 6e).