Low concentration of oxalic acid induces cell death in Arabidopsis roots
High concentration of oxalate is accumulated in the S. sclerotium-infected plant tissues (Marciano et al. 1983). We found that culturing S. sclerotiorum on solid medium resulted in acidification of environments with pH value dropping below 3.0 (Fig. 1). Therefore, we have to consider the acid nature (acidity) of OA in triggering cell death in plants. We first treated 6-d-old Arabidopsis seedlings with OA for 30 min at concentrations ranging from 0.05 to 1.0 mM (pH values ranging from 4.3 to 3.0). The death of Arabidopsis root cells was assayed by staining with Sytox green, a dead cell specific nucleic acid fluorescent dye. Low concentrations of OA induced cell death in Arabidopsis roots in a dose-dependent manner (Fig. 2a). Similar results were obtained when the OA-treated Arabidopsis roots were stained with fluorescein diacetate (FDA), which fluoresces in living cells but not in dead cells (Fig. 2b). Moreover, the dead cells in Arabidopsis roots treated with OA extended bidirectionally over time (Fig. S1).
It is known that propidium iodide (PI) can only be taken up by cells that have a compromised membrane (Autheman et al. 2013). So, it was widely used to stain the necrotic animal cells. However, for plant cells, PI also binds to pectin which forms the primary component of the cell wall (Rieger et al. 2011). So, the use of PI as a dead cell indicator was limited in plants. In our work, we used PI/FDA double staining method to check the OA-induced cell death. Arabidopsis roots were treated with 0.05 mM OA. After 30 min of OA treatment, the dead cells that could not be stained by FDA exhibited strong PI-staining inside the root cells, while in living regions, only the cell walls were stained (Fig. 2c). Consistently, the dead cell region extended bidirectionally in Arabidopsis roots in a time course of 45 min of OA treatment as assayed by PI/FDA double staining (Fig. S2). These results suggest that OA treatment results in emergence of compromised plasma membrane (PM) in Arabidopsis root cells.
Low concentration of OA induces cell death in Arabidopsis seedlings
Then we detected cell death in whole Arabidopsis seedlings upon OA treatment. We found that, after exposure to low concentrations of OA, both Arabidopsis roots and leaves could be stained with Evan's blue, a dead cell specific dye with a maximum absorbance at 600 nm wavelength (Fig. 3a, b). Meanwhile, the leaves became gradually wilting during cell death after OA treatment (Fig. 3a). So, we measured the rate of weight loss (RWL) of Arabidopsis seedlings after OA treatment and found that RWL is positively correlated with the results obtained from Evan's blue staining (Fig. 3b, c), suggesting that RWL could also be used as a readout to quantify the OA-induced cell death in Arabidopsis seedlings.
The acidity of OA induces cell death in Arabidopsis
Next, we sought to determine, between acidity and oxalate of OA, which induces cell death in Arabidopsis at a relatively low concentration (0.5 mM) for a short time of treatment. We found that just like OA (0.5 mM, pH 3.2), hydrochloric acid (HCl, pH 3.2) but not KOX (0.5 mM) induced cell death in Arabidopsis roots as assayed by Sytox green/FDA staining, or in seedlings as assayed by Evan's blue staining or determined by measuring RWL (Fig. 4a-e). These results suggest that the acidic nature of OA induced cell death in Arabidopsis.
The acidity of OA induces an unknown type of cell death in Arabidopsis
DNA-cleavage is a characteristic feature of PCD, which results in a ladder of internucleosomal fragments (Ryerson and Heath 1996). To further characterize the OA-induced cell death, we examined whether OA treatment results in DNA fragmentation. DNA extracted from Arabidopsis seedlings was separated with agarose gel. Consistent with a previous work (Errakhi et al. 2008), DNA isolated from 20 mM KOX-treated seedlings (24h) displayed DNA laddering. However, no DNA fragmentation was detected when seedlings were treated with OA (up to 1.0 mM) for 30 min (Fig. 4f). These results suggest that the type of OA-induced cell death is probably not PCD.
The global gene expression in Arabidopsis upon OA treatment
To further confirm the premise that the acidic nature of OA induces cell death and to further characterize its features, we performed RNA sequencing (RNA-seq) to analyze the genes differentially expressed upon treatment with OA (0.5 mM, pH 3.2), HCl (pH 3.2), KOX (0.5 mM), or H2O (the control) for 60 min. A total of 3074 and 3525 DEGs (differentially expressed genes, |log2FoldChange| > 1 and Padj < 0.05) were identified by comparing transcriptomes of OA- and HCl-treated Arabidopsis to the control, respectively (HCl- or OA-treated DEGs are provided in Dataset S1). However, there were no significant DEGs obtained by comparing KOX transcriptome with the control (Fig. 5a). The expression of a number of DEGs upon OA and HCl treatments determined by RNA-seq was confirmed by RT-qPCR (Fig. S3). Further hierarchical clustering analysis of DEGs upon treatment with OA, HCl, KOX, or H2O showed a high degree of gene coexpression between OA-treated versus HCl-treated transcriptomes, and between KOX-treated versus H2O-treated transcriptomes (Fig. 5b), indicating that OA and HCl have a comparable effect on global gene expression in Arabidopsis.
Furthermore, 2042 genes are overlapping between the 2462 genes (82.9%) induced by HCl and 2284 genes (89.4%) induced by OA, and 526 genes are overlapping between the 1063 genes (49.5%) downregulated by HCl and 790 genes (66.6%) by OA (Fig. 5c, d; Dataset S1). Furthermore, gene ontology (GO) analysis showed that genes up-regulated by both OA and HCl were primarily enriched in the defense and stress related biological process (BP) categories, such as immune system process, defense response to bacterium, defense response to fungus, regulation of defense response, as well as response to chitin, wounding, and jasmonic acid (JA) (Fig. 5e). Moreover, gene set enrichment analysis (GSEA) demonstrates that "defense response", "defense response to bacterium", "response to fungus", "calcium ion binding", and "cell death" gene sets are enriched in both HCl- and OA upregulated genes (Fig. S4).
Interestingly, of the 2042 overlapping genes induced by both HCl and OA treatment, 1468 ones are also up-regulated by S. sclerotiorum infection (Badet et al. 2017) (Fig. 5f). These results suggest that acid treatment, S. sclerotiorum infection, and other plant-biotic interactions trigger some overlapping transcriptomic responses, and OA likely plays a significant role in inducing reprogramming of gene expression in plants upon S. sclerotiorum infection.
The influx of calcium is likely required for the OA-induced cell death in Arabidopsis
Ca2+ influx is a prerequisite of HR-type cell death (Ma et al. 2008; Gao et al. 2013). Our transcriptome analysis also showed that the "calcium ion binding" gene set was enriched in both HCl- and OA-treated Arabidopsis transcriptomes (Fig. 6a, b). RT-qPCR analysis further confirmed the induction of "calcium ion binding" genes by both HCl and OA treatment (Fig. 6c, d). To examine whether calcium influx is required for OA-induced cell death in plants, we tested the effect of inhibiting calcium-channels by two calcium-channel inhibitors, LaCl3 and GdCl3 on OA-triggered cell death. The OA-induced cell death in Arabidopsis roots was obviously blocked by LaCl3 or GdCl3 pretreatment as assayed by staining with Sytox green or FDA (Fig. 7a, b). Similarly, although pretreatments with LaCl3 and GdCl3 resulted in a slight increase in cell death, the OA-induced cell death in Arabidopsis seedlings was significantly blocked by these two chemicals as assayed by Evan's blue staining or determined by measuring RWL (Fig. 7c-e). Moreover, the induction of cell death related genes, such as WRKY46, OXI1, and PAD4 by OA was significantly suppressed by LaCl3 or GdCl3 pretreatment (Fig. 8a-c). Together, these data suggest that the influx of calcium is likely required for the OA-induced cell death in plants.