Long-Distance Electrical and Calcium Signals Evoked by Hydrogen Peroxide in Physcomitrella

Abstract Electrical and calcium signals in plants are some of the basic carriers of information that are transmitted over a long distance. Together with reactive oxygen species (ROS) waves, electrical and calcium signals can participate in cell-to-cell signaling, conveying information about different stimuli, e.g. abiotic stress, pathogen infection or mechanical injury. There is no information on the ability of ROS to evoke systemic electrical or calcium signals in the model moss Physcomitrella nor on the relationships between these responses. Here, we show that the external application of hydrogen peroxide (H2O2) evokes electrical signals in the form of long-distance changes in the membrane potential, which transmit through the plant instantly after stimulation. The responses were calcium-dependent since their generation was inhibited by lanthanum, a calcium channel inhibitor (2 mM), and EDTA, a calcium chelator (0.5 mM). The electrical signals were partially dependent on glutamate receptor (GLR) ion channels since knocking-out the GLR genes only slightly reduced the amplitude of the responses. The basal part of the gametophyte, which is rich in protonema cells, was the most sensitive to H2O2. The measurements carried out on the protonema expressing fluorescent calcium biosensor GCaMP3 proved that calcium signals propagated slowly (>5 µm/s) and showed a decrement. We also demonstrate upregulation of a stress-related gene that appears in a distant section of the moss 8 min after the H2O2 treatment. The results help understand the importance of both types of signals in the transmission of information about the appearance of ROS in the plant cell apoplast.


Introduction
As an early terrestrial plant, the moss Physcomitrella (new botanical name Physcomitrium patens) had to cope with exposure to UV radiation, ozone, wounding and many other abiotic and biotic stressors that enhance the synthesis of reactive oxygen species (ROS) (Tucker et al. 2005, Rensing et al. 2008, 2020. Physcomitrella is an emerging model plant with a fully sequenced genome (Rensing et al. 2008, Lang et al. 2018. It is widely used to study defense mechanisms against environmental stress factors (Frank et al. 2005, Saidi et al. 2009, Koselski et al. 2019, Hoernstein et al. 2023. It was previously demonstrated in laboratory of Plant Physiology and Biophysics, Maria Curie-Sklodowska University, Lublin, Poland that Physcomitrella is capable of generating action potentials (APs) after illumination with light of sufficient (over-threshold) intensity, cooling, glutamate (Glu) treatment, etc. (Koselski et al. 2008(Koselski et al. , 2019(Koselski et al. , 2020. Electrophysiological experiments with the application of ion channel inhibitors and the manipulation of ion gradients across the plasma membrane indicated that Ca 2+ fluxes from external and internal stores were involved in the generation of APs in Physcomitrella. They interacted with K + and, to a lesser extent, with Cl − fluxes. Our recent study on GCaMP mutants allowing the monitoring of changes in the cytoplasmic Ca 2+ concentration ([Ca 2+ ] cyt ) demonstrated that the local Glu application caused an increase in [Ca 2+ ] cyt confined to the site of the Glu treatment, whereas AP was transmitted to distant cells (Koselski et al. 2020). This was rather an unexpected result because, in other plant species examined before, APs or variation potentials (known also as slow wave potentials, SWPs) spread together with Ca 2+ waves and seem mutually dependent , Gilroy et al. 2016, Scherzer et al. 2022). In the model vascular plant Arabidopsis thaliana, two glutamate receptors (GLRs) AtGLR3.3 and AtGLR3.6 were identified as key players in wound-induced SWP generation and transmission (Mousavi et al. 2013, Farmer et al. 2020. Both genes are predominantly expressed in vascular tissues: phloem and xylem contact cells, respectively. Glu (and other amino acids) released from wounded tissues binds to the Ligand Binding Domain subunit of the GLR-like channels facilitating the opening of the channel pore and in consequence initiating a Ca 2+ wave and other downstream responses, like jasmonate (JA) key genes and other stress-related gene expressions (Grenzi et al. 2022).
According to the recent findings in vascular plants, in addition to the ion fluxes facilitating regenerative [Ca 2+ ] cyt wave formation and transmission, a ROS-based component supplements the long-distance signaling system (Baxter et al. 2014, Kurusu et al. 2015, Gilroy et al. 2016, Marcec et al. 2019. NADPH oxidases [respiratory burst oxidase homologs (RBOHs)] localized in the plasma membrane seem to be good candidates that fit with this scheme. RBOHs cause the reduction of oxygen to superoxide anion (O 2 − ), which is quickly converted to hydrogen peroxide (H 2 O 2 ). They possess EF-hand motifs that are activated upon binding of Ca 2+ to enhance ROS production (Wojtaszek 1997, Marcec et al. 2019. Before the machinery of ROS scavenging reduces the ROS level back to normal, its temporary excess affects many metabolic and signaling pathways and their components. H 2 O 2 produced in the apoplast can be quickly transported to the cytosol via aquaporins (Tian et al. 2016). Recently, a specific apoplastic H 2 O 2 receptor, H 2 O 2 -Induced Ca 2+ Increases 1 (HPCA1), was discovered . It is a membrane-spanning protein composed of an extracellular H 2 O 2 sensor and a cytoplasmic kinase component, which is postulated to activate Ca 2+ -permeable channels after H 2 O 2 binding to the sensor (Foyer 2020. Thus, ion channels can be affected from both sides of the plasma membrane.
It was previously demonstrated that ROS can activate ion channels in different plant species classified in different branches of the phylogenetic tree of characean algae (Demidchik et al. 1997) to dicots and monocots (Wu et al. 2015). In vascular plants, H 2 O 2 and hydroxyl radicals affect Ca 2+ -permeable channels in the plasma membrane and in the internal store-the vacuole. In A. thaliana, two classes of non-selective calcium-permeable channels are postulated to be involved: GLR channels activated by Glu channel and cyclic nucleotide-gated channels (CNGCs; Finka et al. 2012, Mousavi et al. 2013. Additionally, annexin1 has been reported to be involved in ROS-induced [Ca 2+ ] cyt elevation (Richards et al. 2014). It has been demonstrated that homologs of Glu channels and CNGCs were involved in signaling processes in Physcomitrella (Finka and Goloubinoff 2014, Ortiz-Ramírez et al. 2017, Koselski et al. 2020.
In addition to calcium-permeable channels, K + channels in the plasma membrane have been found to respond to an increase in the ROS level (Demidchik et al. 2010). Massive K + efflux in ROS-treated plants was one of the first observations of plant response to this type of stress factors (Nassery 1979, Demidchik et al. 2010. In Arabidopsis, guard cell outward rectifying K + and stellar potassium outward rectifier channels have been identified to be responsible for that effect (Demidchik 2018). These channels are also regarded as good candidates to pass an outward current during the AP repolarization phase in excitable plants (Cuin et al. 2018, Scherzer et al. 2022. The goal of our study was to examine for the first time the effect of the H 2 O 2 application on long-distance electrical and calcium signaling systems in the model moss Physcomitrella. It was selected as a relatively simple plant, possessing no typical vascular tissues (Reski 1998), to address the question of evolution of signaling systems in early terrestrial plants. Electrical potential changes and [Ca 2+ ] cyt transients were measured in wild-type (WT) plants and glr1 KO mutants treated with H 2 O 2 . We demonstrated that in the WT plants, H 2 O 2 evoked the electrical signals recorded both in the gametophyte leaves and protonema and the calcium signals predominantly in protonema. The susceptibility of the protonema to the treatment was much higher than in the leaf cells (0.5 mM versus 5 mM). These signals were blocked by La 3+ -a calcium channel inhibitor-and by EDTA-a Ca 2+ chelator. In the glr1 KO mutants, the electrical signals were reduced but not totally blocked, which reveals that GLR channels are not crucial in H 2 O 2 -induced membrane potential changes and Ca 2+ fluxes in Physcomitrella. The involvement of other factors affecting [Ca 2+ ] cyt is discussed. Looking for the downstream effects of the signals, we demonstrate an enhanced expression of stressrelated genes.

Electrical signals in leaf cells
The microelectrode measurements carried out on leaf (phylloid) cells from the Physcomitrella gametophyte proved that the application of H 2 O 2 evoked systemic electrical signals in the form of membrane potential changes transmitted along the plant. The long-distance responses transmitted from the basal part of the plant (from the rhizoid side) and recorded in apical leaf cells are presented in Fig. 1A. The different H 2 O 2 concentrations influenced both the amplitude and the membrane potential characteristic of the response. Stimulation of the basal part of the plant with 0.05 mM and 0.1 mM H 2 O 2 evoked irregular membrane potential changes with similar amplitude, duration and rate of depolarization (Fig. 1A, Table 1). After an increase of the H 2 O 2 concentration to 0.5 mM, significant changes in basic parameters describing the membrane potential changes were obtained (Fig. 1A, Table 1). The cell membrane depolarized twice as fast as in the lower concentrations (14.5 ± 1.6 mV/s, n = 18) and reached more positive values. The amplitude of the depolarization (A) amounted to 68 ± 3 mV (n = 19) and was ∼30 mV higher than after treatment with 0.05 and 0.1 mM H 2 O 2 , respectively. Another difference occurred in the shape of the responses, where after the peak of depolarization evoked by 0.5 mM H 2 O 2 , a several-minute plateau of the membrane potential was recorded.
Letters a-m indicate statistically significant differences. depolarization (R dep = 2.3 ± 0.4 mV/s, n = 8), lower amplitude (22 ± 5 mV, n = 9) and lower voltage at the peaks of the signals (V (max) = −142 ± 8 mV, n = 9). The additional feature was an absence of a characteristic plateau of the membrane potential. The shape of the membrane potential changes resembled the responses obtained after direct stimulation of a cell in the apical part of the plant (Fig. 1C) or removal of the basal part ( Fig. 1D), where a 10-fold increase in the H 2 O 2 concentration (to 5 mM) evoked responses similar to these transmitted from the basal part after 0.5 mM H 2 O 2 administration. These results indicate that the highest susceptibility to H 2 O 2 occurs in the basal part of the plant having protonema cells and rhizoids, which act as an initial place for the generation of long-distance electrical signals propagated toward the plant apex. The results also showed that an increase in the H 2 O 2 concentration facilitates the generation of fully developed responses resembling APs. APs belong to the basic long-distance electrical signals recorded in plant cells whose generation is dependent on an influx of calcium ions into the cytoplasm. We decided to study the dependence of the recorded responses on the presence of an inhibitor of calcium channels (2 mM lanthanum) or a calcium chelator (0.5 mM EDTA), respectively ( Fig. 2A, B, Table 2). In addition to the effects of lanthanum and EDTA, we also examined the possibility to reverse the evoked effects. Each plant was stimulated twice-after initial immersion for 3-4 h in lanthanum or EDTA and then after the exchange of the solution back to the standard solution. Immersion of the plants in lanthanum caused a total blockage of the response in 10 of the 19 tested cells. The responses in the other plants had a significantly reduced amplitude (to 9 ± 2 mV, n = 9).
Lanthanum also shifted the membrane potential (V (0) ) to positive values (to −130 ± 3 mV, n = 19). The depolarization after the lanthanum application was slower (R dep amounted to 0.5 ± 0.2 mV/s, n = 8) and reached more negative values (V (max) = −128 ± 5, n = 9) than in the control plants. After a washout of lanthanum from the measuring chamber, the amplitude of the responses increased to 27 ± 4 mV (n = 4), which indicated that the inhibition of the responses by lanthanum is partially reversible. EDTA totally blocked the responses in 9 of the 18 tested cells. Similar to lanthanum, EDTA evoked the reduction of the amplitude (to 13 ± 4 mV, n = 9), the shift of the membrane potential to more positive values (to −122 ± 5 mV, n = 18) and the reduction of the rate of depolarization (to 1.4 ± 0.7 mV/s, n = 7). The shift of the membrane potential after EDTA was partially reversible since a washout resulted in a restoration of the membrane potential close to the values measured in standard solution (−144 ± 4 mV, n = 18).

Electrical signals in protonema cells
The results of the measurements carried out on the leaf cells indicated that it is hard to evoke electrical signals in such cells, and the responses can start mainly from the protonema cells and/or rhizoids-probably the target for the action of H 2 O 2 .
Letters a-f indicate statistically significant differences.
To study this hypothesis, we decided to examine the effect of H 2 O 2 on membrane potential changes in protonema cells. In these measurements, stimulation was carried out by microinjection of H 2 O 2 in three regions: initially into a chain of protonema cells adjacent to the tested cell with the inserted microelectrode, then directly into the tested cell and at the end into the gametophyte base.
The results confirmed that, although electrical signals in protonema cells are transmitted from cell to cell, the responses recorded in the same cell differed depending on the region of stimulation (Fig. 4, Table 3, Supplementary Video S1). The stimulation of the gametophyte base was the most effective, as it evoked electrical signals in each tested cell immediately upon stimulation. In comparison to the electrical signals recorded in the leaf cells, the signals from the protonema reached a higher amplitude (81 ± 4 mV, n = 8) and the duration measured in half of the amplitude (t 1/2 = 486 ± 100 s, n = 8). Surprisingly, weaker effects were achieved by the direct stimulation of the tested cell or the stimulation of adjacent cells located close to the tested cell. This method of stimulation evoked electrical signals with a smaller amplitude (66 ± 6 mV, n = 8) than when H 2 O 2 was applied in the gametophyte base. The responses recorded in the protonema after the direct or indirect stimulation also exhibited a low depolarization rate (1.6 ± 0.6 mV/s, n = 8), compared to the responses in the basal part of the gametophyte (9.3 ± 2.5 mV/s, n = 8) or those recorded in the leaf cells. Cell-to-cell transmission of the electrical signal evoked by stimulation of the single protonema cells was rarely recorded and occurred only in 2 out of 18 tested plants ( Supplementary  Fig. S1, Supplementary Video S2).
One of the candidates responsible for cell-to-cell communication in plants is the GLR, which can act as a non-selective calcium-permeable channel. The participation of the GLR in the transmission of electrical signals recorded in our experiments was tested with the use of a glr1 mutant of Physcomitrella (Ppglr KO ). In the Ppglr1 KO mutants, as in the WT, electrical signals propagated from the base of the gametophyte to the leaf and protonema cells (Figs. 3, 4, Supplementary Video S3). In comparison to the WT, electrical signals recorded in the leaf cells of Ppglr1 KO had a reduced amplitude (to 46 ± 3 mV, n = 26) and a lower depolarization rate (9.6 ± 1.2 mV/s, n = 24). The amplitude and rate of depolarization of responses in the protonema cells from Ppglr1 KO recorded after the stimulation of the basal gametophyte part reached similar values to those recorded in the WT protonema cells (78 ± 3 mV, n = 10 and 8.9 ± 3 mV/s, n = 10, respectively). As in the WT protonema cells, the direct or indirect stimulation of the cells from Ppglr1 KO evoked a smaller amplitude (63 ± 4 mV, n = 10) and a lower depolarization rate (1.1 ± 0.3 mV/s, n = 10) than after the stimulation of the basal part of the Ppglr1 KO gametophyte. All these changes in the parameters of membrane potential changes in the protonema cells are presented in Table 3.

Calcium signals
Calcium signals were recorded in plants expressing GCaMP3, i.e. a fluorescent calcium indicator. Fluorescence measurements of calcium signals recorded after the stimulation of the basal part of the plant indicated that the application of 0.5 mM H 2 O 2 evoked calcium waves which propagated from cell to cell in thread-like protonema cells. In contrast to electrical signals, the calcium signals were slower and appeared a few minutes after the application of the stimulus (Fig. 5, Supplementary  Video S4). In some plants, the calcium signals did not start at the stimulation site but appeared at some distance from the stimulation, acting as a new source for calcium waves. In such a situation, calcium wave propagation was observed in both directions-into and from the stimulation site. The rate of calcium signal propagation along the protonema toward the site of the stimulation was similar to that in the opposite direction, reaching 5.5 ± 0.7 μm/s (n = 4) and 5.2 ± 0.3 μm/s (n = 7), respectively. A characteristic trait of some recorded calcium signals was its decrement with distance ( Supplementary Fig. S2, Supplementary Video S5). As in the case of electrical signals, the observation of calcium signals in the stimulated leaves was possible after the increase of the H 2 O 2 concentration to 5 mM (Supplementary Video S6). In turn, in contrast to electrical signals, calcium signals in leaf cells after the application of 0.5 mM H 2 O 2 to the basal part of the plant were observed occasionally only in several cells from single leaves but never (in none of 30 tested plants) in all leaves (Supplementary Video S7).

Expression of genes in distant regions of the moss
To analyze whether the electrical signal that was triggered upon local treatment would be accompanied by changes in gene expressions in distant tissues, we selected candidate genes and performed quantitative real-time PCR (qPCR). The candidates were selected based on homology to Arabidopsis genes differentially expressed after local stress stimulus (Zandalinas et al. 2019). In that study, a high light stimulus was locally imposed on selected Arabidopsis rosette leaves, and changes in gene expression in the locally treated and distant, nontreated leaves (local and systemic responses, respectively) were −155 ± 2 (n = 10) −77 ± 3 (n = 10) 78 ± 3 (n = 10) 294 ± 25 (n = 9) 8.9 ± 3 c,d (n = 10) Direct stimulation and indirect stimulation correspond to the measuring variants presented in Fig. 4. Letters a-d indicate statistically significant differences. Explanations are provided in Table 1. Letters a-h indicate statistically significant differences.
investigated by RNAseq analysis. From this published list, we selected candidate genes that were upregulated in the systemic response if their expression was also inducible via H 2 O 2 treatment and depended on the function of the RBOH D (RBOHD). Here, the gene encoding a galacturonosyltransferaselike 10 protein (AT3G28340) was 5-fold upregulated 5 min after the light stress treatment (Zandalinas et al. 2019  one gene (Supplementary Fig. S3, Pp3c16_25090V3.1) exhibited a significant increase in gene expression (P = 0.0305) and hence was selected for further analysis. For Pp3c5_28420V3.1, an increase in gene expression was detectable ( Supplementary  Fig. S3), but the difference was not significant. Nevertheless, this gene was included in further experiments. For these two selected candidates, we tested whether both genes are regulated in gametophore sections distant from a local treatment with H 2 O 2 . Apices of gametophores were harvested after treatment with 0.5 mM H 2 O 2 at the base for 8 min (Supplementary Fig. S4), as well as corresponding untreated apices. A significant upregulation of Pp3c16_25090V3.1 (P = 0.00184) was observed (Fig. 6), whereas an upregulation of Pp3c5_28420V3.1 was indicated but was not statistically significant (P = 0.06075). Thus, a local stress stimulus (here H 2 O 2 at the gametophore base) triggers alteration of gene expression within 8 min in distant sections (here apex) in Physcomitrella.

Discussion
Plant cell signaling is based mainly on the generation and transmission of different types of signals including electrical, Fig. 6 Gene expression analysis by qPCR of the selected gene candidates. Apices were used from untreated gametophores and gametophores treated at the base with 0.5 mM H 2 O 2 for 8 min. Dots represent biological replicates (mean values from three technical replicates). Bars are the mean values from the four biological replicates with standard deviation. Relative expression (2 (−∆CT) ) of Pp3c5_28420V3.1 and Pp3c16_25090V3.1 is calculated against the reference genes L21 (Pp3c13_2360V3.1) and LWD (Pp3c22_18860V3.1) according to Livak and Schmittgen (2001). Significance levels are based on a one-way ANOVA with subsequent Tukey's HSD post hoc test (**P < 0.01); n.s. = not significant. calcium and ROS. Investigations of the relationship between these signals allow unraveling signaling pathways triggered in response to biotic and abiotic stresses. In this study, we tried to answer the question of whether the external application of ROS (H 2 O 2 ) evokes a systemic response in the form of electrical and calcium signals to analyze the interdependence of these signals and their impact on gene expression in distant regions.
The importance of ROS in the plant defense system has been evidenced by many experiments (Huang et al. 2019). One of the first studies focused on the involvement of ROS in cell-tocell communication in Arabidopsis (Miller et al. 2009) proved that local stimuli produced ROS waves propagating through the plant at a rate of 8.4 cm/min and the response was dependent on the RBOHD gene encoding plant NADPH oxidase involved in the production of ROS. The rate of propagation of such RBOHDrelated ROS signals recorded in different tissues ranged from ∼400 to 1400 μm/s and was dependent on the type of stress (Choi et al. 2017). The accumulation of extracellular ROS was observed along the path of propagation (Miller et al. 2009), indicating that each cell along the path is able to activate RBOHD and release ROS in adjacent cells. In such an auto-propagation process, named ROS-induced ROS release, H 2 O 2 generated by RBOHD can be regarded as a long-distance signal. However, in the study by Miller et al. (2009), the external application of H 2 O 2 did not evoke a ROS wave indicating the involvement of some other signal molecules in the propagation of the ROS wave.
Calcium is a candidate for such a signal molecule and, together with ROS, can cooperate in long-distance transmission of information about stimuli. For example, calciumdependent protein kinase CPK5 playing a role in plant immunity is important for cell-to-cell communication based on ROS waves (Dubiella et al. 2013). There are also other ways of activation of RBOH proteins by calcium, including direct binding of this ion to the EF-hand motif on the RBOH protein (Kimura et al. 2012) or binding of phosphatidic acid to the same protein whose accumulation in the cell can be induced by calcium (Zhang et al. 2009). In addition to the calcium-induced ROS release process, ROS-induced calcium release is possible. For example, it has been shown that ROS can activate different calcium-permeable channels, e.g. hyperpolarization-activated Ca 2+ channels in root cells (Demidchik et al. 2007), Ca 2+ influx channels in guard cells (Pei et al. 2000) or Ca 2+ -permeable channels regulated by annexin1 (Richards et al. 2014). Recent studies carried out on Arabidopsis allowed the discovery of a plasma membrane leucine-rich repeat receptor kinase, i.e. HPCA1, which links the perception of apoplastic H 2 O 2 with Ca 2+ signaling . One of the most probable mechanisms for the cooperation between calcium and ROS in long-distance transmission of signals is that an increase in cytoplasmic calcium evoked by a local stimulus can induce the production of ROS via Ca 2+ -RBOH interactions and the accumulation of ROS in the apoplast. The transport of ROS from the apoplast to the cell, probably carried out by aquaporins, would evoke ROS-induced calcium release, which in turn could induce ROS production .
As shown by the results of our study, H 2 O 2 -evoked calcium signals are not systemic responses, such as those observed earlier in P. patens under osmotic stress or salt stimulation (Storti et al. 2018). Calcium waves with a velocity of ∼400 μm/s in response to local salt stress were also recorded in Arabidopsis , where the blockage of calcium channels by lanthanum inhibited not only the calcium waves but also the ROS-regulated transcriptional marker (ZAT12), indicating that calcium and ROS waves are closely linked. Assuming the interconnection of calcium signals and ROS waves, it is unlikely that the extracellular application of H 2 O 2 in our experiments evoked the ROS wave since the H 2 O 2 -evoked calcium signals admittedly appeared at some distance from the stimulation site but propagated with a decrement (Fig. 5, Supplementary  Videos S4 and S5, Supplementary Fig. S2). The other feature of the H 2 O 2 -evoked calcium signals in Physcomitrella, i.e. the slower rate of propagation than in the case of the ROS wave (>5 μm/s), also indicates a low probability of the appearance of a calcium-associated ROS wave. The ability to evoke local calcium signals but not self-propagating calcium waves by the extracellular H 2 O 2 application implies the engagement of some other factors taking part in the transmission of information about such stimuli. It seems that, in P. patens, electrical signals, which appear at the moment of stimulation and propagate along the whole plant, are a proper carrier of the information about H 2 O 2 enhancement. There still remains the question of the relationships between electrical, calcium and ROS signals.
The relationship between ROS and electrical signals was previously confirmed in Arabidopsis mutants lacking RBOHD, where propagating electrical signals after heat stress or high light had a significantly reduced amplitude or were totally blocked (Suzuki et al. 2013). Taking into account the ability of extracellular ROS (mainly H 2 O 2 ) to activate different ion channels (Demidchik 2018), including calcium-permeable channels that can participate in long-distance electrical signals, it seems that ROS can promote electrical signals along the ROS wave path. This hypothesis can also be supported by the similar velocity of different types of plant self-propagating electrical signals to the velocity of the ROS wave [like APs (20-400 cm/min) and system potentials (5-10 cm/min); Zimmermann et al. 2009]. In Physcomitrella, the H 2 O 2 -evoked electrical signals propagated immediately from the basal part along the whole plant (Figs. 1, 4, Supplementary Video S1), indicating that the signals appeared before the slowly propagating calcium waves. In cooperation with ROS and calcium waves, electrical signals, which are probably the fastest carriers of information, can initiate the whole array of plant defense responses. It is also probable that, by reaching cells distant from the stimulation site, electrical signals initiate a calcium wave response only in some cells, in which the threshold of calcium signal generation is lower than in other cells. This assumption is proved by the observation of calcium waves in the protonema cells, where, at a distance from the stimulation site in the basal part of the plant, only a few protonema cells in the thread-like chain generated calcium signals propagating in two directions-into and out of the site of stimulation (Fig. 5, Supplementary Fig. S2, Supplementary  Videos S4 and S5).
One of the candidates for ion channels participating in long-distance electrical signals in plants is GLR-like channels, which was confirmed in experiments carried out on Arabidopsis (Mousavi et al. 2013, Salvador-Recatalà et al. 2014, Salvador-Recatalà 2016. Wound-induced electrical signals recorded in this species were dependent on clade III GLRs (GLR 3.3 and GLR 3.6) and played a crucial role in the distal production of JAs taking part in plant defense responses (Mousavi et al. 2013). It is a matter of discussion if those channels are directly responsible for Ca 2+ fluxes since they are predominantly expressed in endomembranes (Nguyen et al. 2018, Farmer et al. 2020. In Physcomitrella, only two GLR genes have been identified (Verret et al. 2010). Both genes (PpGLR1 and PpGLR2) are paralogs to GLR clade III from Arabidopsis (De Bortoli et al. 2016) and encode channels participating in chemotaxis and reproduction of Physcomitrella (Ortiz-Ramírez et al. 2017). A patch-clamp study indicated that GLR1 is a calcium-permeable ion channel localized in the cell membrane and partially inhibited by GLR antagonists (Ortiz-Ramírez et al. 2017). These data indicate that GLR1 in Physcomitrella, similar to GLR 3.3 and GLR 3.6 from Arabidopsis, can be important in the transmission of long-distance electrical signals. This hypothesis was not fully confirmed in our study since H 2 O 2 evoked long-distance electrical signals in the Physcomitrella glr1 KO mutants (Figs. 3, 4,  Supplementary Video S3); however, the response differed in the amplitude ( Table 2). Such a 'weak' effect of the GLR1 knockout suggests that, in addition to GLR channels, other channels must be engaged in the long-distance propagation of ROSinduced electrical signals. The most important channels taking part in the generation of electrical signals in Physcomitrella are probably calcium permeable, given the blockage of H 2 O 2evoked responses by the calcium channel inhibitor (lanthanum) or the calcium chelator (EDTA) (Fig. 2). It is also probable that, in addition to ion channels in the plasma membrane, a significant role in the transmission of electrical signals is assigned to intracellular calcium channels. One of the bestknown intracellular channels permeable, e.g., to calcium is the two-pore channel 1 (TPC1) located in the vacuolar membrane, the tonoplast (Dadacz-Narloch et al. 2011). A patchclamp study carried out in our laboratory demonstrated that TPC channels in Physcomitrella vacuoles conduct Ca 2+ currents (Koselski et al. 2013). Involvement of the channel in the transmission of long-distance calcium waves has been demonstrated , but there is no information about the role of TPC1 in the transmission of electrical signals. As calciumpermeable channels, TPC1 may be involved in the release of calcium from intracellular compartments (Qudeimat et al. 2008) leading to calcium-induced calcium release (CICR), a desired phenomenon for long-distance signal transmission, but this assumption arouses controversy (Ward andSchroeder 1994, Pottosin et al. 1997). The role of H 2 O 2 in CICR is questionable since microelectrode ion flux estimation and patch-clamp studies carried out on vacuoles from Beta vulgaris demonstrated that H 2 O 2 suppressed the Ca 2+ efflux from the vacuole and slow vacuolar currents carried by TPC1 (Pottosin et al. 2009).
Effectively, the treatment with 0.5 mM H 2 O 2 at the base of gametophores was sufficient to trigger an increase in gene expression in the apex of a component of the homogalacturonan biosynthesis (Fig. 6, Pp3c16_25090V3.1). Here, the homologous candidate Pp3c5_28420V3.1 was not significantly upregulated although an increasing trend was detectable. This agrees with publicly available data, indicating that both genes are not co-expressed. Galacturonosyltransferase-like proteins such as Pp3c16_25090V3.1 act in the pectin assembly (homogalacturonan biosynthesis) of the primary cell wall (reviewed by Loix et al. 2017). Under oxidative stress, cell wall pectins also represent a source for the biosynthesis of ascorbic acid (Valpuesta andBotella 2004, García-Caparró s et al. 2021), which in turn is used to detoxify ROS such as H 2 O 2 . Both selected candidates are homologs of an Arabidopsis isoform (AT3G28340) whose gene expression is regulated via ROS waves (Zandalinas et al. 2019). However, it should be noted that Galacturonosyltransferase-like proteins comprise a large gene family with at least 21 members in Arabidopsis and 17 genes in Physcomitrella (Van Bel et al. 2022) and clear ortholog relations are not yet resolved. The expression of the Arabidopsis homolog (AT3G28340) was 5-fold increased only in systemic leaves distant from a local high light stress impulse (Zandalinas et al. 2019), and those data further indicate that the increase of expression was a response to a ROS wave. In contrast, the two selected homologs in Physcomitrella were not regulated by high light stress (Supplementary Fig. S5), but their expression increased after heat stress. Heat stress in plants is accompanied by the elevated production of ROS such as H 2 O 2 (reviewed by Mittler et al. 2022). In summary, these data show that the gene expression of at least one of the two Physcomitrella candidates (Pp3c16_25090V3.1) is responsive to ROS. Consequently, the increase of expression in the untreated apex (Fig. 6) was likely based on a propagating ROS wave.
Taken together, our study demonstrates differences in the generation and propagation of H 2 O 2 -evoked electrical and calcium signals in the model moss Physcomitrella. Many of the applied variants of measurements indicated that the basal part of the gametophyte is the most excitable region, probably because a large number of protonema cells are juvenile cells from an early stage of new gametophyte development. In comparison to the leaf cells, the responses in the protonema exhibited higher amplitudes and lasted longer, which may indicate higher susceptibility of such cells to ROS (Tables 2 and 3,  Figs. 1, 4). The main difference between the electrical and calcium signals was the velocity of propagation, which was higher in the electrical signals. The other difference was the ability to propagate without a decrement; in contrast to the electrical signals, calcium signals diminished with distance, even if they appeared at some distance from stimuli ( Supplementary  Fig. S2, Supplementary Video S5). Given these differences, we propose that H 2 O 2 -evoked long-distance electrical signals are the first to reach distant regions of the plant and activate calcium signals, but not in every cell. The protonema cells were the most susceptible to calcium signals, indicating that the signals play a key role in young and developing cells. A similar observation was reported in our previous work focused on Glu-evoked calcium signals (Koselski et al. 2020). Electrical and calcium signals were not the only effect of H 2 O 2 application. Gene expression analysis proved that apart from the signals, an increase of stress-related gene expression is observed in not stimulated distant regions of plants (Fig. 6). The increase in the gene expression appeared after 8 min-the time close to the duration of electrical signals. The similar time scale of electrical signals and gene expression raises the question about interdependence of both phenomena.

Cultivation of plant material used for electrophysiological and calcium imaging analysis
Physcomitrella gametophytes were grown on Knop solid agar medium (Reski and Abel 1985) in 160-mm diameter Petri dishes. The WT plants and knockout glr1 mutants (Ppglr1 KO ) were used for the analysis of membrane potential changes in leaf and protonema cells with the use of the microelectrode technique. A Physcomitrella mutant expressing GCaMP3 was used for the fluorescence imaging of changes in the calcium concentration. The plants were grown in a growing chamber (Conviron Adaptis A1000, Conviron, Winnipeg, Canada) at a photoperiod of 16-/8-h light/dark, with light intensity 50 μmol/m 2 s and at a temperature set to 23 ∘ C.  Reski and Abel (1985), and 10 mL per liter of an ME stock solution (309 mg/l H 3 BO 3 , 845 mg/l MnSO 4 × 1 H 2 O, 431 mg/l ZnSO 4 × 7 H 2 O, 41.5 mg/l KI, 12.1 mg/l Na 2 MoO 4 × 2 H 2 O, 1.25 mg/l CoSO 4 × 5 H 2 O and 1.46 mg/l Co(NO 3 ) 2 × 6 H 2 O) was also prepared as described by Egener et al. (2002) and Schween et al. (2003). The suspension culture was dispersed weekly with an ULTRA-TURRAX (IKA, Staufen, Germany) at 18,000 rpm for 90 s.

Cultivation of plant material used for gene expression analysis
Hydroponic cultures of Physcomitrella gametophores were cultivated as described by Hoernstein et al. (2023). Glass rings covered with mesh (PP, 250-m mesh, 215-m thread, Zitt-Thoma GmbH, Freiburg, Germany) were prepared as described by Erxleben et al. (2012). Protonema suspension was adjusted to a final density of 440 mg/l (dry weight per volume) as described by Decker et al. (2017) and evenly distributed in equal volumes on the mesh surface. Glass rings covered with protonema were placed in Magenta ® vessels (Sigma-Aldrich, St. Louis, MO, USA). KnopME medium was supplemented until it touched the bottom of the mesh. The medium was changed every 4 weeks.
Suspension cultures and hydroponic cultures were cultivated in standard light conditions (55 μmol photons/m 2 s) at 22 ∘ C in a 16-/8-h light/dark cycle.

Measurements of membrane potential in leaf cells
The method of membrane potential measurements with microelectrodes was similar to that employed by Koselski et al. (2020). Plastic Petri dishes used in the microelectrode measurements were divided into two chambers with a barrier. The barrier had a small (1-mm width) gap sealed with Vaseline. Before the experiments, the plants were incubated for 3-6 h at a light intensity of 50 μmol/m 2 in a bath solution containing (in mM) 1 KCl, 1 CaCl 2 , 50 sorbitol and 2 HEPES, pH 7.5 (buffered by Tris). The stimuli were introduced by the application of 500 μl of a bath solution supplemented with 0.5 mM H 2 O 2 . The lanthanum chloride or EDTA influence on the H 2 O 2 -evoked responses described in this study was assessed by the application of one of these substances into one of the Petri dish compartments containing the basal part of the plant. Borosilicate glass capillaries (1B150F-6, World Precision Instruments, Sarasota, FL, USA) were used to make micropipettes with the use of a P-30 micropipette puller (Sutter Instrument Co., Novato, CA, USA), filled with 100 mM KCl and connected to the FD223 electrometer (World Precision Instruments, Sarasota, FL, USA). A Sensapex SMX (SensApex, Oulu, Finland) electronic micromanipulator was used for positioning and insertion of the microelectrode. The reference electrode was composed of an Ag/AgCl 2 wire inside a plastic tube filled with 100 mM KCl and ending with a porous tip. The measured data were acquired by a Lab-Trax-4 device (World Precision Instruments, Sarasota, FL, USA) working with LabScribe2 software, which also allowed analyses of the data. The analysis of statistical differences was performed in SigmaStat 4.0 (Systat Software Inc., Palo Alto, CA, USA). Recordings of changes in the membrane potential were recorded with a 2-Hz data collection frequency. Figures were prepared with the use of Sigma Plot 9.0 (Systat Software Inc, California, USA) and CorelDraw 12 (Corel Corporation, Ottawa, Canada) software.

Measurements of membrane changes in protonema cells
The plants were prepared in the same way as for the measurements carried out on the leaves. The ME was positioned and inserted using a PatchStar (Scientifica, East Sussex, UK) micromanipulator and observed under an Olympus IX71 (Olympus, Tokyo, Japan) microscope with a camera (Artcam-500MI, Tokyo, Japan) working with QuickPHOTO Camera software (version 2.3, Promicra, Prague, Czech Republic). A CellTram Vario microinjector (Eppendorf, Hamburg, Germany) with a borosilicate glass micropipette (with a diameter of ∼3 μm) was used for the application of H 2 O 2 onto the cell surface. As in our previous paper (Koselski et al. 2020), 1 mM methyl blue was used for staining of the stimulating solution, which allows the observation of its dispersion. Live recording of membrane potentials visible in LabScribe 3 software and microscope camera images were recorded by Open Broadcaster Software (ver. 23.2.1; Open Broadcaster Software, Massachusetts, USA).

Fluorescence calcium imaging
Films and images of calcium concentration changes in the GCaMP3 P. patens mutants were recorded with the use of NIS-Elements AR software (ver. 5.20.00; Nikon, Tokyo, Japan) working with a Nikon Eclipse Ti fluorescence microscope equipped with a Nikon Plan UW 2× WD:7.5 objective and a Nikon DS-Ri2 camera (Nikon, Tokyo, Japan). The excitation of fluorescence was provided by a Prior lumen 200 metal arc lamp (Prior Scientific Instruments Ltd, Cambridge, England) and a 495-nm dichroic mirror with a standard GFP excitation filter 470 ± 20 nm. Images of fluorescence emission were recorded with 1-s exposure and a standard GFP fluorescence emission barrier filter 525 ± 50 nm (Nikon, Tokyo, Japan). Stimulation of the plants was carried out by microinjection using a CellTram Vario microinjector (Eppendorf, Hamburg, Germany). The standard solution supplemented with H 2 O 2 was visible due to tinting by 0.025 μM fluorescein.

Treatment of plants used for gene expression analysis
Gametophores from hydroponic cultures were used for the treatment with H 2 O 2 . Half of the gametophores from one glass ring were cut approximately in the middle, and only the upper half of the gametophores (∼100 mg) was used as a control sample. Cut gametophore apices were gently dried by dabbing with filter paper. The glass ring with the remaining uncut gametophores was transferred into a new Magenta ® vessel containing KnopME with 0.5 mM H 2 O 2 . After 8 min, the upper half of the remaining uncut gametophores (treated sample) was harvested as described earlier. In total, gametophores from four independent hydroponic cultures (four biological replicates) were sampled.

RNA extraction and qPCR
Extraction of RNA was done with the innuPREP Plant RNA Kit (Analytic Jena, Jena, Germany) using the extraction buffer 'PL' . Of the total RNA, 5 μg was treated with DNAseI (Thermo Fisher Scientific, Waltham, MA, USA) at 37 ∘ C for 1 h, and integrity of the RNA was checked on agarose gels. Two-microgram DNAseI digested RNA was used for reverse transcription using the TaqMan TM Reverse Transcription kit (N8080234, Thermo Fisher Scientific) with random hexamer primers. Reverse transcription was performed at 42 ∘ C for 1 h. A non-transcribed control without the addition of MultiScribe TM Reverse Transcriptase enzyme was included. Primers for the qPCR were designed using the Primer3Plus software (https://www.primer3plus.com/; Untergasser et al. 2012) with qPCR settings, and an efficiency of 2 was confirmed using a 1:2 dilution series of cDNA. Melting curve analysis was performed to exclude the presence of off-targets. qPCR was performed in 96-well plates using the SensiFast TM SYBR No-ROX Kit (Bioline, Memphis, Tennessee, USA) in technical triplicates for each biological replicate. Fifty-nanogram cDNA was used for each technical triplicate, and the PCR reaction was performed in a LightCycler ® 480 (Roche, Basel, Switzerland). −RT and water controls for each primer pair were included. The PCR reaction was performed in 45 cycles with a melting temperature of 60 ∘ C. The expression analysis of the genes Pp3c5_28420V3 and Pp3c16_25090V3 [genes of interest (GOIs)] was performed as described by Bohlender et al. (2020) in relation to the housekeeping (reference) genes L21 coding for the ribosomal protein L21 (Pp3c13_2360V3, Beike et al. 2015) and LWD coding for the transducin/WD40 repeat-like superfamily protein (Pp3c22_18860V3, Schuessele et al. 2016). All primers are listed in Supplementary Table S1.
The expression levels were calculated relative to the reference genes according to Livak and Schmittgen (2001) using the software for the LightCycler® 480 (V1.5.0, Roche, Basel, Switzerland). Relative expression is represented as 2 (−∆CT) with ∆CT = CT [GOI] − CT [reference] . Fig. 6 is created and statistics were calculated in R (R Core Team 2022). Statistical significance was tested via one-way ANOVA with subsequent Tukey's honest significant difference (HSD) post hoc test. Significance was accepted at P < 0.05.

Supplementary Data
Supplementary data are available at PCP online.

Data Availability
The data underlying this article are available in the article and in its online supplementary material. The other data underlying this article will be shared on request to the corresponding author.