Skip to main content

Phospholipase D activation is required for 1-aminocyclopropane 1-carboxylic acid signaling during sexual reproduction in the marine red alga Neopyropia yezoensis (Rhodophyta)

Abstract

Background

1-aminocyclopropane 1-carboxylic acid (ACC) is the immediate precursor of the plant hormone ethylene. However, recent studies have suggested that ACC also acts as a signaling molecule to regulate development and growth independently from ethylene biosynthesis. In red algae, ACC stimulates the switch from a vegetative to a sexual reproductive phase. However, despite evidence that ACC signaling in plants and algae is widespread, the mechanistic basis of the ACC signaling pathway remains unknown.

Results

We demonstrate that exogenous ACC increased the activity of phospholipase D (PLD) and induced the accumulation of PLD transcripts in the marine red alga Neopyropia yezoensis. The product of PLD, the lipid second messenger phosphatidic acid (PA), also increased in response to ACC. Furthermore, the pharmacological inhibition of PLD by 1-butanol blocked ACC-induced spermatangia and carpospore production, but the inactive isomer t-butanol did not. In addition, 1-butanol prevented ACC-induced growth inhibition and inhibited transcript accumulation of genes upregulated by ACC, including extracellular matrix (ECM)-related genes, and alleviated the transcriptional decrease of genes downregulated by ACC, including photosynthesis-related genes.

Conclusions

These results indicate that PLD is a positive regulator of sexual cell differentiation and a negative regulator of growth. This study demonstrates that PLD and its product, PA, are components of ACC signaling during sexual reproduction in N. yezoensis.

Peer Review reports

Background

1-aminocylopropane-1-carboxylic acid (ACC) is known as the direct precursor of the plant hormone ethylene, which regulates a wide range of developmental processes and stress responses [1, 2]. Exogenous ACC application has been used as a proxy for ethylene in numerous experiments over decades of research on ethylene-associated signaling. However, there is emerging evidence from the model plant Arabidopsis thaliana that ACC itself may have a signaling role that is distinct from its role as a precursor of ethylene [3, 4]. For example, an acs octuple mutant with deficient ACC production exhibited a remarkable decrease in a seed set that was not observed in ethylene-insensitive mutants [5]. Other studies have uncovered that ACC signaling independent of the influence of ethylene can regulate vegetative processes such as early root and hypocotyl elongation in seedlings [6] and cell wall function [7]. In addition, recent studies suggest that ACC functions independently of ethylene in the basal land plant Marchantia polymorpha [8, 9]. Given the apparent significance of ACC signaling, elucidating the mechanistic basis of ACC signaling pathways is an important goal of contemporary research.

Phospholipase D (PLD) hydrolyzes membrane lipids, producing phosphatidic acid (PA) and free hydrophilic head groups, such as choline or ethanolamine [10]. PA is now regarded as a lipid signaling molecule that regulates numerous physiological processes in eukaryotes [11, 12]. In higher plants, PLD and its product PA mediate the signaling of various plant hormones, including abscisic acid (ABA), ethylene, jasmonic acid, and salicylic acid (SA) [10, 13]. To date, 1-butanol, an antagonist of PLD-dependent PA production, has been used to investigate the function of PLDs in plant hormone signaling. For example, one study of ABA and gibberellin (GA) signaling in barley aleurone tissue found that both pathways were inhibited by 1-butanol [14]. In addition, 1-butanol has also been found to block SA-induced stomatal movement in guard cells in A. thaliana [15].

The red alga Pyropia/Neopyropia (formerly Porphyra), which belongs to the Bangiales, is a significant marine crop that is harvested to produce nori, an edible seaweed commonly used to wrap sushi and onigiri. During the sexual life cycle of Pyropia/Neopyropia, the blade gametophytes bear non-flagellated male (spermatia) and female (carpogonia) gametes on the gametophytes. Fertilization occurs when the female gametes are retained on the gametophyte, and successive cell divisions produce clones of the zygote. It is the carpospores that then grow into filamentous sporophytes. Recent research from our group has demonstrated that the application of ACC induced gametogenesis and growth suppression in the monoecious species Neopyropia yezoensis and the dioecious species Pyropia pseudolinearis [16, 17]. Exogenous ACC analogs also promoted sexual reproduction in the same manner as ACC, whereas ethephon, an ethylene-releasing compound, did not stimulate sexual reproduction in N. yezoensis [18, 19]. In RNA-seq, transcripts associated with cell division, vesicular-trafficking, and extracellular matrix (ECM) were found to be up-regulated in gametophytes treated with ACC, while transcripts involved in translation, plastid transcription and photosynthesis were downregulated [16, 17]. In addition, the application of ACC generated reactive oxygen species (ROS) mediated by NADPH oxidase activity, increased ascorbate (AsA) synthesis, and decreased glutathione (GSH) synthesis [18]. Taken together, these results suggest that ACC may play a role as a signaling molecule independent from its role in ethylene signaling and that it may be involved in the regulation of sexual reproduction through changing the redox state of Pyropia/Neopyropia.

To date, several studies have suggested that PLD is involved in development processes and abiotic stress responses in Neopyropia species [20, 21]. In addition, the previous RNA-seq data [16] showed transcript accumulation of PLD gene upregulated by ACC. We therefore, investigated whether PLD and PA were required for the ACC response during sexual reproduction in N. yezoensis. We provide evidence that PLD and PA are required for the signal transduction events ultimately leading to ACC-induced sexual reproduction in N. yezoensis.

Methods

Algal materials and chemical treatments

Gametophytic blades of N. yezoensis strain TU-1 were cultured in a medium of sterile vitamin-free Provasoli’s enriched seawater (PES [22];) at 15 °C under cool-white fluorescent lamps at 40 μmol photons m− 2 s− 1 irradiance with a photoperiod regime of 10 h light:14 h dark. For PLD inhibitor experiments, we grew immature gametophytes (approximately 20 mm in blade length) that were microscopically determined to bear only vegetative cells in 90-mm-diameter petri dishes with 40 mL PES. Plates were then placed in a shaking incubator. Gametophytes were then treated with 0 or 50 μM ACC (Tokyo Chemical Industry, Tokyo, Japan) in the presence of 1-butanol (n-butyl alcohol) or t-butanol (tert-butyl alcohol) for 10 days. After 10 d, gametophytes that had formed spermatangia were observed under a Leica DM 5000 B microscope (Leica Microsystems, Tokyo, Japan), since the carpogonia of N. yezoensis are almost indistinguishable from surrounding vegetative cells. Formation of spermatangia in N. yezoensis generally begins in apical regions. Thus, in this study, gametophyte maturity level was evaluated using a four level scale as follows: “a”: an upper part without spermatangia, “b”: an upper part with spermatangia exhibiting slight discoloration, “c”: an upper part with spermatangia exhibiting clear discoloration, and “d”: an upper part with degraded cell walls and released spermata. In addition, after 14 d, the number of discharged carpospores attached to the bottom of the dishes was determined under a microscope as an index of the rate of formation of female gametes. Gametophyte blade lengths were also measured after 10 d, and the growth rate was calculated as the mean percentage of length increase per day using the following formula: Growth rate = [100(BLt − BL0)/BL0]/t, where BL0 = initial blade length, BLt = blade length as measured on day t of culture, and t = culture time in days.

PLD activity assay

Phospholipase D (PLD) activity was determined using an Amplex™ Red Phospholipase D Assay Kit (Thermo Fisher Scientific K.K., Tokyo, Japan), according to the manufacturer’s protocol with minor modifications. For the PLD activity assay, vegetative gametophytes from blades approximately 20 mm in length (0.05 g fresh weight; FW) were cultured in a 100 mL culture medium containing 50 μM ACC for 0, 3, and 7d. Samples were then ground in liquid nitrogen with a pestle and mortar. For each sample, the resulting homogenate was added to 0.5 mL of 50 mM Tris-HCl buffer at pH 8.0 and centrifuged at 4 °C for 5 min at 15,000×g. Enzyme activity was detected by adding 100 μL of the supernatant to the working solution of the Amplex Red reagent. The fluorescence was detected using a spectrofluorometer (FB-750, Jasco, Tokyo, Japan) at an excitation wavelength of 545 nm and an emission wavelength of 610 nm to decrease interference from autofluorescence caused by photosynthetic pigments. The resulting fluorescence was subtracted from the background fluorescence derived from samples that did not contain the Amplex Red reagent working solution. Relative PLD activity was calculated as a ratio of the measurement taken at 0 d to that of another taken after ACC treatment. All data are presented as mean ± standard deviation (SD) of five biological replicates.

PA measurement

Total PA content was measured by a coupled enzymatic reaction system using the Total Phosphatidic Acid Assay Kit (Cell Biolabs, Inc., San Diego, CA.). Experimental procedures were performed according to the manufacturer’s protocol with minor modifications. In brief, vegetative gametophytes from blades approximately 20 mm in length (0.05 g FW) were cultured in a 100 mL culture medium containing 50 μM ACC for 0, 3, and 7d. Samples were then ground in liquid nitrogen with a pestle and mortar. For each sample, the resulting homogenate was added to 0.75 mL methanol. Next, we added 1.15 mL 1 M NaCl and 1.25 mL chloroform to each sample and mixed the constituents thoroughly. After centrifugation at 4 °C for 10 min at 1500×g, the upper aqueous phase was discarded, and the lower chloroform phase was washed twice with a PEU solution, which was prepared by mixing 50 mL chloroform, 50 mL methanol, and 45 mL 1 M NaCl. Finally, the lower organic phase was dried under a gentle stream of nitrogen and dissolved in the provided Assay Buffer. Fluorescent signals were detected using a spectrofluorometer (FB-750, Jasco, Tokyo, Japan) at an excitation wavelength of 540 nm and an emission wavelength of 590 nm. Relative PLD production was quantified as a ratio of that PLD production at 0 d after ACC treatment. All data are presented as mean ± SD of five biological replicates.

Bioinformatics analysis

Sequences of two NyPLD (NyPLD1 and NyPLD2) genes were retrieved from N. yezoensis genome sequence data [23]. Additionally, one new NyPLD gene (NyPLD3) was identified in the genome sequence data (Accession ASM982973v1 in the NCBI GenBank Database). The conserved domains of NyPLD proteins were confirmed by analysis using the “SMART” protein architecture research tool (http://smart.embl-heidelberg.de/).

Transcriptional analysis

RNA extraction and qRT-PCR were performed as described by [24]. Total RNA was extracted using an RNeasy Plant Mini Kit (Qiagen, Hilden, Germany) following the manufacturer’s instructions. Extracted RNA was purified using a TURBO DNA-free kit (Invitrogen/Life Technologies, Carlsbad, CA) to obtain DNA-free RNA. First strand cDNA was synthesized from 0.5 μg total RNA using the PrimeScript II First Strand cDNA Synthesis Kit (TaKaRa Bio, Shiga, Japan). cDNA was diluted 10-fold for qRT-PCR analysis; for each reaction, 1.0 μl of the diluted cDNA was used as a template in a 20 μL reaction volume containing KOD SYBR® qPCR Mix (TOYOBO, Osaka, Japan). qRT-PCR was then performed as per the manufacturer’s instructions using a LightCycler® 480 System (Roche Diagnostics, Basel, Switzerland) under the following conditions: 30 s at 95 °C followed by 40 cycles of 5 s at 95 °C and 31 s at 55 °C. mRNA levels were then calculated using the 2Ct method and were normalized to the observed expression level of 18S ribosomal RNA (Py18SrRNA) gene [16]. Relative expression levels over time were calculated as a ratio of the observed mRNA level at a given time point to the mRNA level present at 0 d after ACC treatment. qRT-PCR was performed in triplicate. Table S1 lists primers that were used for these analyses.

Statistical analysis

All data were expressed as the mean ± SD. For the PLD activity, PA content, and growth rate experiments, Mann–Whitney U tests were performed to determine whether differences in mean were statistically significant. For all analyses, p < 0.05 (significant) or p < 0.01 (highly significant) were used as thresholds of statistical significance.

Results

ACC stimulates PLD activation and PA production

We measured PLD activity in response to ACC to determine the role played by PLD in ACC-activated signal transduction in N. yezoensis. Gametophytes that were cultured after 7d of ACC treatment formed spermatangia exhibiting slight or clear discoloration on the upper parts, but gametophytes treated with ACC for 3d did not. Gametophytes treated with ACC displayed increased PLD activity at 3 d after treatment (1.70-fold) with peak activity at 7 d (2.03-fold) (Fig. 1). Next, we examined the effects of ACC treatment on PA production during sexual reproduction. The PA content of gametophytes then increased (1.22-fold) 3 d after ACC treatment, and rose again at 7 d (1.39-fold) (Fig. 2). Taken together, these results indicate that PA production caused by ACC treatment is dependent on PLD activity.

Fig. 1
figure 1

1-aminocylopropane-1-carboxylic acid (ACC)-stimulated PLD activity in N. yezoensis. Gametophytes were treated with ACC for 3, or 7d. Results are presented as measured relative activity compared to untreated gametophytes (0 d). Asterisks indicate a statistically significant difference at P < 0.05 between the control and the specified treatment. All data are presented as mean ± SD of five independent experiments

Fig. 2
figure 2

1-aminocylopropane-1-carboxylic acid (ACC) stimulated the production of phosphatidic acid (PA) in N. yezoensis. Gametophytes were treated with ACC for 3, or 7d. Results are presented as the relative amount compared to untreated gametophytes (0 d). Asterisks indicate a statistically significant difference at P < 0.05 between the control and the specified treatment. All data are presented as mean ± SD of five independent experiments

Characterization of PLD genes and their expression in response to ACC

In addition to our previous RNA-seq data, evidence of increased PLD activity and PA levels in gametophytes treated with ACC led us to characterize the PLD genes found in N. yezoensis. We identified three genes encoding PLDs in the N. yezoensis genome sequence (Table 1 and Fig. 3). NyPLD protein lengths ranged from 797 to 1201 amino acid residues, and BLAST analysis revealed that their protein sequences exhibited 79.5–90.4% identity with PLDs found in the red marine alga Porphyra umbilicalis. NyPLD proteins comprised highly conserved catalytic HKD motifs “HxKxxxxD” (HKD1 and HKD2) found at the C-terminal. Furthermore, we found C2 domains, which bind Ca2+ and other effectors including phospholipids and inositol phosphates [25], at the N-terminals of all NyPLDs. Given these characteristic domains, we identified these three NyPLDs as members of the C2-PLD family.

Table 1 The information of NyPLDs
Fig. 3
figure 3

Schematic representations of the conserved domain and motifs of NyPLDs. Domains and motifs were determined by searching the NyPLD protein sequences using the SMART protein architecture tool. C2; Ca2+/phospholipid binding domain, HKD; conserved catalytic region. Numbers correspond to amino acid positions from the first methionine residue

Next, we analyzed the expression of NyPLDs in response to exogenous ACC treatment using real-time PCR (Fig. 4). We found that the mRNA transcript levels of NyPLD1 had increased by 3 d after treatment (3.04-fold), and peaked at 7 d (7.79-fold). By contrast, the expression levels of NyPLD2 and NyPLD3 were only slightly upregulated and nearly unchanged, respectively, in response to ACC treatment. These results suggested that NyPLD1 is responsible for the increased PLD activity caused by ACC treatment.

Fig. 4
figure 4

1-aminocylopropane-1-carboxylic acid (ACC)-activated PLD expression. RNA samples were prepared from gametophytes that were treated with 50 μM ACC for 3, or 7d, and expression levels were assessed using the Ny18SrRNA gene for normalization. Results are presented as relative expression compared to untreated gametophytes (0 d). All data are presented as mean ± SD of three independent experiments

The PLD inhibitor 1-butanol blocks the ACC-induced PLD activity response

We next examined whether inhibition of PLD enzyme activity would block the induction of sexual reproduction caused by ACC treatment. PLDs can use primary alcohols as an acceptor of a phosphatidyl moiety instead of H2O, and in the presence of a primary alcohol such as 1-butanol, PLD catalyzes a transphosphatidylation reaction that leads to the formation of phosphatidylalcohol instead of PA [26]. Therefore, 1-butanol is generally used to inhibit PA formation by PLD. By contrast, tertiary alcohols such as t-butanol are not substrates of this transphosphatidylation reaction and can, therefore, be used as a control for any nonspecific butanol effect [26]. Moreover, we performed preliminary experiments that indicated that the application of 1-butanol in concentrations greater than 0.04% dramatically inhibited the growth of N. yezoensis. Thus, we treated vegetative gametophytes with 0.02% 1-butanol or t-butanol to assess whether their effects on ACC action differed.

We found that gametophytes treated with ACC produced colorless spermatangia on the upper parts of the thalli after 10 d of treatment; this was consistent with previous observations (Fig. 5). By contrast, the presence of 1-butanol inhibited ACC-induced spermatangia production, but the presence of the inactive t-butanol isomer did not (Fig. 5). In this study, we evaluated gametophyte maturity on a four-point scale described above (Fig. 6). When gametophytes were treated with a control or with butanol alone, all thalli (100%) were classified as “a” or “b” (Fig. 7A). By contrast, thalli treated with ACC alone (100%) or co-treated with ACC and t-butanol (93%) were classified as “c” or “d” (Fig. 7A). Moreover, when gametophytes were treated simultaneously with ACC and 1-butanol, most thalli (75%) were classified as “a” or “b” and showed suppression of the ACC-induced progression in maturity, likely caused by the 1-butanol treatment (Fig. 7A).

Fig. 5
figure 5

Effects of the PLD inhibitor 1-butanol in the presence of 1-aminocylopropane-1-carboxylic acid (ACC). Gametophytes were cultured in media containing 0 (Control) or 50 μM ACC (ACC). Co-treatments included a 10-day application of ACC with 0.02% 1-butanol (1-But + ACC), 0.02% t-butanol (t-But + ACC), 0.02% 1-butanol (1-But - ACC), and 0.02% t-butanol (t-But - ACC) alone for 10 d. We found that the thalli treated with ACC alone or by co-treatment with ACC and t-But formed many spermatangia in the upper regions and that these were often clear or discolored. Other thalli exposed to co-treatment with ACC and 1-But exhibited no discoloration in the upper regions. Scale bar = 10 mm

Fig. 6
figure 6

Mature levels of N. yezoensis for evaluation in the pharmacological experiments. a: the gametophyte without spermatangia, b: the gametophyte with spermatangia exhibited slight discoloration, c: the gametophyte with spermatangia exhibited clear discoloration, and d: the gametophyte degraded cell walls and released spermata. Arrowheads indicate spermatangia. Scale bar = 50 μm

Fig. 7
figure 7

A PLD inhibitor blocks 1-aminocylopropane-1-carboxylic acid (ACC)-induced sexual reproduction. A Suppression of progression of maturity in the presence of the PLD inhibitor. Gametophytes were cultured in media containing 0 (Control) and 50 μM ACC (ACC) or a combined treatment of ACC and 0.02% 1-butanol (1-But + ACC), ACC and 0.02% t-butanol (t-But + ACC), 0.02% 1-butanol without ACC (1-But), or 0.02% t-butanol (t-But) alone. Plant maturity levels are illustrated as Fig. 6. Data were calculated in four independent experiments with four thalli for each condition (n = 16). B Alleviation of ACC-induced growth inhibition in the presence of the PLD inhibitor. Data are expressed as mean ± SD of four independent experiments with four thalli for each condition (n = 16). Asterisks indicate significant differences at P < 0.01 between the controls and treatments. C We also observed suppression of ACC-induced carpospore formation in the presence of the PLD inhibitor. The number of carpospores released from gametophytes was counted under microscope after 14d of culture. Data are expressed as means ± SD of two independent experiments with four thalli for each condition (n = 8)

Next, we examined gametophyte growth rates. When gametophytes were cultured in media containing 50 μM ACC and without butanol, the growth rate was 6.0% d− 1, whereas untreated (control) gametophytes showed a growth rate of 15.5% d− 1 (Fig. 7B). However, we found that 1-butanol treatment prevented ACC-linked inhibition of growth (i.e., gametophytes treated showed a growth rate of 11.8% d− 1), whereas t-butanol did not show a significant effect on the inhibitory action of ACC (4.5% d− 1) (Fig. 7B). The mean number of carpospores released from the gametophytes treated with ACC was approximately 1412, whereas the mean of untreated control gametophytes was zero, and the mean of gametophytes co-treated with ACC and 1-butanol was approximately 1 (Fig. 7C). Taken together, these results indicate that the inactive t-butanol isomer had no inhibitory effect on ACC-induced carpospore production, in contrast to the significant effect of 1-butanol. Therefore, we conclude that treatment with 1-butanol significantly blocked the induction of sexual reproduction by exogenous ACC treatment in N. yezoensis.

Finally, we examined the inhibitory effects of 1-butanol on the expression of ACC-responsive genes using qRT-PCR. Based on previous RNA-seq data [16], we selected three highly upregulated and three highly downregulated genes as representative ACC-responsive genes (Table S2). As illustrated in Fig. 8, the expression of ACC-upregulated genes including ECM-related genes, such as NyGT14 and NyVWA1, strongly increased (31.96- to 34.73-fold, respectively) after 3 days of ACC treatment. However, the presence of 1-butanol strongly inhibited this increase in gene expression by treatment with ACC (we observed increases of only 1.88- to 2.31-fold, respectively). By contrast, the presence of t-butanol showed no effect on the expression levels of ACC-responsive genes. Moreover, the presence of 1-butanol alleviated the decrease in expression of genes that were downregulated by ACC, including photosynthesis-related genes, such as NySIG1 and NyPsbQ, whereas t-butanol had no effect on the ACC-induced downregulation of these genes.

Fig. 8
figure 8

The PLD inhibitor affects the expression of ACC-responsive genes in N. yezoensis. RNA samples were prepared from gametophytes that were cultured in media containing 0 (Control) and 50 μM ACC (ACC) and from those containing co-applications with 0.02% 1-butanol (1-But + ACC), 0.02% t-butanol (t-But + ACC), 0.02% 1-butanol (1-But), or 0.02% t-butanol (t-But) alone for 3 d. Expression levels were determined using the Ny18SrRNA gene for normalization. NyGT14, NyVWA1, and NyVWA2 were selected as representative genes that were upregulated by ACC, whereas NySIG1, NyPsbQ, and NyAly were selected as representative genes that were downregulated by ACC genes. These results are presented as relative expression level at the given time point relative to that present in untreated gametophytes (0 d). All data are presented as mean ± SD from three independent experiments

Discussion

Recent studies have proposed that in land plants and red algae, the ethylene precursor ACC may play a signaling role that is independent of the ethylene signaling pathway [4]. However, little is known regarding the mechanisms involved in regulating ACC signaling. In this study, we demonstrated that PLD activation is required for ACC signaling during sexual reproduction of N. yezoensis. This report is the first to describe the link between PLD and plant hormone signaling in red algae. In land plants, PLD has been thought to play a role in mediating a wide range of physiological processes, including hormone action, stress response, and cell differentiation [27, 28]. Plant PLDs comprise a large family in their genome; for example, there are 12 members of the PLD family in A. thaliana, and expression and functional studies indicate that plant PLDs have functional diversification [29]. In this study, ACC stimulated gene expression of only NyPLD1, suggesting that NyPLD1 is responsible for ACC signal transduction, but NyPLD2, 3 are associated with other cellular process. Previous studies suggest that NyPLD2, 3 may be involved in spore migration and heat stress responses in N. yezoensis [20, 21]. Application of a reverse genetic approach has the potential to answer the question why N. yezoensis possess three PLDs.

PLD-derived PA has been linked to vesicular trafficking processes including Golgi transport, endocytosis, and exocytosis [30]. This is important, since PA can promote membrane curvature at cisternal rims in the Golgi apparatus, which is required for vesicle formation and budding [31, 32]. For example, 1-butanol treatment decreased the accumulation of PA-mediated secretory vesicles in the apex, resulting in a loss of apical polarity in the pollen tube and inhibited tip growth [33]. Additionally, a previous data suggests that PLD and PA modulate auxin responses through the regulation of vesicle trafficking in A. thaliana [34]. In species related to N. yezoensis, the active formation of fibrous vesicles has been thought to play a role in protein turnover and cell wall formation during the differentiation of the spermata and carpospores [35]. In addition, our previous study demonstrated that the expression levels of vesicular trafficking-related genes, such as charged multivesicular body protein (CHMP), increased in response to ACC treatment during sexual reproduction in N. yezoensis and P. pseudolinearis [17]. CHMP proteins are components of the endosomal-sorting complex required for transport (ESCRT) [36], and PLD-produced PA can participate in the assembly of members of the ESCRT machinery in mammal cells [37]. These findings suggest that further research on the involvement of NyPLDs in vesicular trafficking processes may be important in understanding molecular mechanisms involved in red algal reproduction.

ECMs from macroalgae, which are commonly referred to as the cell wall, are complex assemblages of cellulose, various hemicelluloses, and unique sulfated polysaccharides [38]. In this study, 1-butanol treatment alleviated the ACC-downregulation of NyAly that was characterized as the gametophyte-specific expressed alginate lyase [39], which may be involved in the change of the composition of the cell wall during N. yezoensis life cycle. We also found that 1-butanol inhibited the transcription of ACC- upregulated genes, including ECM-related genes such as NyGT14. NyGT14 is a member of the putative glycosyltransferase family and encodes a homolog of N-acetyl glucosaminyl transferase that can transfer N-acetyl glucosamine (GlcNAc) to an acceptor substrate. Changes in the expression patterns of glycosyltransferase genes have been observed to lead to notable structural alteration of the N-glycans on the cell surface through modification of ECM proteins; this enzyme is thereby associated with various biological events including cell adhesion, migration, and cell differentiation [40]. In addition to NyGT14, 1-butanol inhibited the upregulation of two ACC-inducible genes that encode von Willebrand A (VWA) domain-containing proteins. The majority of well-characterized VWA domains are found in ECM proteins, and these domains are involved in protein–protein (e.g., receptor–ligand) interactions [41]. Although little is known about the role ECM proteins play in ACC signaling in N. yezoensis, we have also observed that the expression levels of putative ECM genes, including four NySPL-encoding spondin domain-containing proteins, were upregulated in response to ACC treatment [42].

Treatment with 1-butanol was found to alleviate transcriptional repression of genes that are downregulated by ACC treatment, including genes related to photosynthesis and NyPsbQ and NySIG. With respect to photosynthesis, photosystem II (PSII) is the key protein complex involved in light-energy conversion reactions [43]. Moreover, the components of PSII can be classified into core proteins, low-molecular-mass proteins, extrinsic oxygen-evolving complex (OEC) proteins, and the light-harvesting complex II protein [44]. In Arabidopsis, the loss of the PsbQ, one of the OEC proteins, induces significant changes in PS II function, particularly in plants grown in low light conditions; this indicates that the PsbQ protein is required for photoautotrophic growth [45]. Regarding the regulation of photosynthesis-related genes, PsbQ is encoded in the nuclear genome, whereas psbA and psbD encode PS II reaction center core proteins, D1 and D2 polypeptides, and are encoded by plastids [46]. The transcription of plastid genes, including psbA and psbD, is regulated by nuclear-encoded sigma factors; these bind directly to the core plastid RNA polymerase [47]. In one study in Arabidopsis, the authors observed that the transcription levels of psbA and psbD in SIG1 mutants were significantly lower than those of the wild type [48]. In the present study, the presence of 1-butanol alleviated ACC-induced downregulation of NySIG, a gene that encodes a homolog of plant sigma factors. Taken together, the alleviation of the downregulation of photosynthesis-related genes by 1-butanol may help mitigate ACC-induced growth inhibition.

Plant hormones generate ROS through the activation of NADPH oxidase, and this reaction is crucial event in the regulation of plant growth and development [49]. A previous study suggested that ROS generation plays an important role in ACC-induced sexual reproduction in N. yezoensis [18]. NADPH oxidases from higher plants possess an extended N-terminus that contain a Ca2+ binding EF-hand motif; the activation of these enzymes requires Ca2+ binding to the motif [50]. By contrast, red algal NADPH oxidases lack an EF-hand motif which serves as a calcium binding site [51,52,53], suggesting that red algal NADPH oxidases perceive the calcium signal through other than EF-hand motif [54]. There is evidence to show that plant NADPH oxidases are also directly activated by PLD-derived PA [55], whereas animal NADPH oxidases are indirectly activated through phosphorylation by PA-activated protein kinases [56]. Whether PA stimulates NADPH oxidase from red algae through direct or indirect means or by some more complicated interaction should be investigated in the future.

Plant PLDs are divided into two subclasses. PLDs containing the C2 domain belong to the C2-PLD subclass, whereas those with domains homologous to the domains of phox (PX) and/or pleckstrin (PH) are classified as PX/PH-PLDs [29]. By searching all domains of the protein sequences studied here, we noted that all NyPLDs contained a C2 domain. This makes intuitive sense: This domain is involved in the binding of Ca2+ and is known to be required for PLD activation. Previous studies in both plants and animals have demonstrated that ACC can elicit Ca2+ currents via glutamate receptors [57,58,59]. Therefore, ACC-induced Ca2+ signaling may play an important role in PLD activation in N. yezoensis.

Conclusion

We demonstrated that increases in PLD activity and PA production follow the accumulation of NyPLD1 transcripts in response to exogenous ACC treatment. In addition, the pharmacological inhibition of PLD activity alleviated ACC-induced growth repression and blocked ACC-induced formation of spermatangia and carpospores, indicating that PLD plays a negative role for plant growth and a positive role with respect to the differentiation of sex cells. These findings contribute significantly to our understanding of ACC signal transduction during sexual reproduction in red algae.

Availability of data and materials

All data supporting the findings in this study are presented within the manuscript.

References

  1. Lin ZF, Zhong SL, Grierson D. Recent advances in ethylene research. J Exp Bot. 2009;60:3311–36.

    Article  CAS  Google Scholar 

  2. Van De Poel B, Smet D, Van Der Straeten D. Ethylene and hormonal cross talk in vegetative growth and development. Plant Physiol. 2015;169:61–72. https://doi.org/10.1104/pp.15.00724.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Polko JK, Kieber JJ. 1-aminocyclopropane 1-carboxylic acid and its emerging role as an ethylene-independent growth regulator. Front Plant Sci. 2019;10:1602. https://doi.org/10.3389/fpls.2019.01602.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Van de Poel B. Ethylene's fraternal twin steals the spotlight. Nat Plants. 2020;6:1309–10. https://doi.org/10.1038/s41477-020-00796-8.

    Article  CAS  PubMed  Google Scholar 

  5. Tsuchisaka A, Yu G, Jin H, Alonso JM, Ecker JR, Zhang X, et al. A combinatorial interplay among the 1-aminocyclopropane-1-carboxylate isoforms regulates ethylene biosynthesis in Arabidopsis thaliana. Genetics. 2009;183:979–1003. https://doi.org/10.1534/genetics.109.107102.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Vanderstraeten L, Depaepe T, Bertrand S, Van Der Straeten D. The ethylene precursor ACC affects early vegetative development independently of ethylene signaling. Front Plant Sci. 2019;10:1591. https://doi.org/10.3389/fpls.2019.01591.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Tsang DL, Edmond C, Harrington JL, Nuehse TS. Cell wall integrity controls root elongation via a general 1-aminocyclopropane-1-carboxylic acid-dependent, ethylene-independent pathway. Plant Physiol. 2011;156:596–604. https://doi.org/10.1104/pp.111.175372.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Li DD, Flores-Sandoval E, Ahtesham U, Coleman A, Clay JM, Bowman JL, et al. Ethylene-independent functions of the ethylene precursor ACC in Marchantia polymorpha. Nat Plants. 2020;6:1335–44. https://doi.org/10.1038/s41477-020-00784-y.

    Article  CAS  PubMed  Google Scholar 

  9. Katayose A, Kanda A, Kubo Y, Takahashi T, Motose H. Distinct functions of ethylene and ACC in the basal land plant Marchantia polymorpha. Plant Cell Physiol. 2021;62:858–71. https://doi.org/10.1093/pcp/pcab042.

    Article  CAS  PubMed  Google Scholar 

  10. Munnik T. Phosphatidic acid: an emerging plant lipid second messenger. Trends Plant Sci. 2001;6:227–33. https://doi.org/10.1016/s1360-1385(01)01918-5.

    Article  CAS  PubMed  Google Scholar 

  11. Li M, Hong Y, Wang X. Phospholipase D- and phosphatidic acid-mediated signaling in plants. Biochim Biophys Acta. 2009;1791:927–35. https://doi.org/10.1016/j.bbalip.2009.02.017.

    Article  CAS  PubMed  Google Scholar 

  12. Wang XM, Devalah SP, Zhang WH, Welti R. Signaling functions of phosphatidic acid. Prog Lipid Res. 2006;45:250–78. https://doi.org/10.1016/j.plipres.2006.01.005.

    Article  CAS  PubMed  Google Scholar 

  13. Zhao J. Phospholipase D and phosphatidic acid in plant defence response: from protein-protein and lipid-protein interactions to hormone signalling. J Exp Bot. 2015;66:1721–36. https://doi.org/10.1093/jxb/eru540.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Ritchie S, Gilroy S. Abscisic acid signal transduction in the barley aleurone is mediated by phospholipase D activity. Proc Natl Acad Sci U S A. 1998;95:2697–702. https://doi.org/10.1073/pnas.95.5.2697.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Kalachova T, Lakovenko O, Kretinin S, Kravets V. Involvement of phospholipase D and NADPH-oxidase in salicylic acid signaling cascade. Plant Physiol Biochem. 2013;66:127–33. https://doi.org/10.1016/j.plaphy.2013.02.006.

    Article  CAS  PubMed  Google Scholar 

  16. Uji T, Matsuda R, Takechi K, Takano H, Mizuta H, Takio S. Ethylene regulation of sexual reproduction in the marine red alga Pyropia yezoensis (Rhodophyta). J Appl Phycol. 2016;28:3501–9. https://doi.org/10.1007/s10811-016-0904-6.

    Article  CAS  Google Scholar 

  17. Yanagisawa R, Sekine N, Mizuta H, Uji T. Transcriptomic analysis under ethylene precursor treatment uncovers the regulation of gene expression linked to sexual reproduction in the dioecious red alga Pyropia pseudolinearis. J Appl Phycol. 2019;31:3317–29.

    Article  CAS  Google Scholar 

  18. Uji T, Endo H, Mizuta H. Sexual reproduction via a 1-aminocyclopropane-1-carboxylic acid-dependent pathway through redox modulation in the marine red alga Pyropia yezoensis (Rhodophyta). Front Plant Sci. 2020;11:60. https://doi.org/10.3389/fpls.2020.00060.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Endo H, Mizuta H, Uji T. α-Aminoisobutyric acid mimics the effect of 1-aminocyclopropane-1-carboxylic acid to promote sexual reproduction in the marine red alga Pyropia yezoensis (Rhodophyta). J Appl Phycol. 2021;33:1081–7.

    Article  CAS  Google Scholar 

  20. Li L, Saga N, Mikami K. Ca2+ influx and phosphoinositide signalling are essential for the establishment and maintenance of cell polarity in monospores from the red alga Porphyra yezoensis. J Exp Bot. 2009;60:3477–89.

    Article  CAS  Google Scholar 

  21. Uji T, Sato R, Mizuta H, Saga N. Changes in membrane fluidity and phospholipase D activity are required for heat activation of PyMBF1 in Pyropia yezoensis (Rhodophyta). J Appl Phycol. 2013;25:1887–93.

    Article  CAS  Google Scholar 

  22. Provasoli L. Media and prospects for the cultivation of marine algae. In: Watanabe A, Hattori A, editors. Culture and collections of algae, Proc U S-Japan Conf, Hakone, September 1966. Jpn Soc Plant Physiol. Tokyo; 1968. p. 63–75.

    Google Scholar 

  23. Nakamura Y, Sasaki N, Kobayashi M, Ojima N, Yasuike M, Shigenobu Y, et al. The first symbiont-free genome sequence of marine red alga, susabi-nori (Pyropia yezoensis). PLoS One. 2013;8:e57122.

    Article  CAS  Google Scholar 

  24. Uji T, Gondaira Y, Fukuda S, Mizuta H, Saga N. Characterization and expression profiles of small heat shock proteins in the marine red alga Pyropia yezoensis. Cell Stress Chaperones. 2019;24:223–33.

    Article  CAS  Google Scholar 

  25. Corbalan-Garcia S, Gómez-Fernández JC. Signaling through C2 domains: more than one lipid target. Biochim Biophys Acta. 2014;1838:1536–47. https://doi.org/10.1016/j.bbamem.2014.01.008.

    Article  CAS  PubMed  Google Scholar 

  26. Munnik T, Arisz SA, deVrije T, Musgrave A. G protein activation stimulates phospholipase D signaling in plants. Plant Cell. 1995;7:2197–210. https://doi.org/10.2307/3870162.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Xue H, Chen X, Li G. Involvement of phospholipid signaling in plant growth and hormone effects. Curr Opin Plant Biol. 2007;10:483–9. https://doi.org/10.1016/j.pbi.2007.07.003.

    Article  CAS  PubMed  Google Scholar 

  28. Testerink C, Munnik T. Molecular, cellular, and physiological responses to phosphatidic acid formation in plants. J Exp Bot. 2011;62:2349–61. https://doi.org/10.1093/jxb/err079.

    Article  CAS  PubMed  Google Scholar 

  29. Qin C, Wang X. The Arabidopsis phospholipase D family. Characterization of a calcium-independent and phosphatidylcholine-selective PLD zeta 1 with distinct regulatory domains. Plant Physiol. 2002;128:1057–68.

    Article  CAS  Google Scholar 

  30. Jones D, Morgan C, Cockcroft S. Phospholipase D and membrane traffic. Potential roles in regulated exocytosis, membrane delivery and vesicle budding. Biochim Biophys Acta. 1999;1439:229–44. https://doi.org/10.1016/s1388-1981(99)00097-9.

    Article  CAS  PubMed  Google Scholar 

  31. Roth MG. Lipid regulators of membrane traffic through the Golgi complex. Trends Cell Biol. 1999;9:174–9. https://doi.org/10.1016/s0962-8924(99)01535-4.

    Article  CAS  PubMed  Google Scholar 

  32. Freyberg Z, Siddhanta A, Shields D. 'Slip, sliding away': phospholipase D and the Golgi apparatus. Trends Cell Biol. 2003;13:540–6. https://doi.org/10.1016/j.tcb.2003.08.004.

    Article  CAS  PubMed  Google Scholar 

  33. Monteiro D, Liu Q, Lisboa S, Scherer GE, Quader H, Malho R. Phosphoinositides and phosphatidic acid regulate pollen tube growth and reorientation through modulation of [Ca2+]c and membrane secretion. J Exp Bot. 2005;56:1665–74.

    Article  CAS  Google Scholar 

  34. Li G, Xue HW. Arabidopsis PLDzeta2 regulates vesicle trafficking and is required for auxin response. Plant Cell. 2007;19:281–95. https://doi.org/10.1105/tpc.106.041426.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Wang SJ, Xu ZD. Ultrastructural studies on the reproductive-organs of Porphyra haitanensis Chang,T.J. et Zheng,B.F. Hydrobiologia. 1984;116:213–7.

    Google Scholar 

  36. Adoro S, Park KH, Bettigole SE, Lis R, Shin HR, Seo H, et al. Post-translational control of T cell development by the ESCRT protein CHMP5. Nat Immunol. 2017;18:780.

    Article  CAS  Google Scholar 

  37. Egea-Jimenez AL, Zimmermann P. Phospholipase D and phosphatidic acid in the biogenesis and cargo loading of extracellular vesicles. J Lipid Res. 2018;59:1554–60. https://doi.org/10.1194/jlr.R083964.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Kloareg B, Badis Y, Cock JM, Michel G. Role and evolution of the extracellular matrix in the acquisition of complex multicellularity in eukaryotes: a macroalgal perspective. Genes. 2021;12:1059.

    Article  CAS  Google Scholar 

  39. Inoue A, Mashino C, Uji T, Saga N, Mikami K, Ojima T. Characterization of an eukaryotic PL-7 alginate lyase in the marine red alga Pyropia yezoensis. Curr Biotechnol. 2015;4:240–8. https://doi.org/10.2174/2211550104666150915210434.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Isaji T, Kariya Y, Xu QS, Fukuda T, Taniguchi N, Gu JG. Functional roles of the bisecting glcnac in integrin-mediated cell adhesion. Methods Enzymol Glycobiol. 2010;480:445–59. https://doi.org/10.1016/s0076-6879(10)80019-9.

    Article  CAS  Google Scholar 

  41. Whittaker CA, Hynes RO. Distribution and evolution of von Willebrand/integrin a domains: widely dispersed domains with roles in cell adhesion and elsewhere. Mol Biol Cell. 2002;13:3369–87. https://doi.org/10.1091/mbc.e02-05-0259.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Uji T, Ueda S, Mizuta H. Identification, characterization, and expression analysis of spondin-like and fasciclin-like genes in Neopyropia yezoensis, a marine red alga. Phycology. 2022;2:45–59.

    Article  Google Scholar 

  43. Nickelsen J, Rengstl B. Photosystem II assembly: from Cyanobacteria to plants. Annu Rev Plant Biol. 2013;64:609–35. https://doi.org/10.1146/annurev-arplant-050312-120124.

    Article  CAS  PubMed  Google Scholar 

  44. Lu Y. Identification and roles of photosystem II assembly, stability, and repair factors in Arabidopsis. Front Plant Sci. 2016;7. https://doi.org/10.3389/dpls.2016.00168.

  45. Yi XP, Hargett SR, Frankel LK, Bricker TM. The PsbQ protein is required in Arabidopsis for photosystem II assembly/stability and photoautotrophy under low light conditions. J Biol Chem. 2006;281:26260–7. https://doi.org/10.1074/jbc.M603582200.

    Article  CAS  PubMed  Google Scholar 

  46. Minagawa J, Takahashi Y. Structure, function and assembly of photosystem II and its light-harvesting proteins. Photosynth Res. 2004;82:241–63. https://doi.org/10.1007/s11120-004-2079-2.

    Article  CAS  PubMed  Google Scholar 

  47. Chi W, He BY, Mao J, Jiang JJ, Zhang LX. Plastid sigma factors: their individual functions and regulation in transcription. BBA-Bioenergetics. 2015;1847:770–8. https://doi.org/10.1016/j.bbabio.2015.01.001.

    Article  CAS  PubMed  Google Scholar 

  48. Macadlo LA, Ibrahim IM, Puthiyaveetil S. Sigma factor 1 in chloroplast gene transcription and photosynthetic light acclimation. J Exp Bot. 2020;71:1029–38. https://doi.org/10.1093/jxb/erz464.

    Article  CAS  PubMed  Google Scholar 

  49. Xia XJ, Zhou YH, Shi K, Zhou J, Foyer CH, Yu JQ. Interplay between reactive oxygen species and hormones in the control of plant development and stress tolerance. J Exp Bot. 2015;66:2839–56. https://doi.org/10.1093/jxb/erv089.

    Article  CAS  Google Scholar 

  50. Oda T, Hashimoto H, Kuwabara N, Akashi S, Hayashi K, Kojima C, et al. Structure of the N-terminal regulatory domain of a plant NADPH oxidase and its functional implications. J Biol Chem. 2010:1435–45. https://doi.org/10.1074/jbc.M109.058909.

  51. Hervé C, Tonon T, Collén J, Corre E, Boyen C. NADPH oxidases in eukaryotes: red algae provide new hints! Curr Genet. 2006;49:190–204. https://doi.org/10.1007/s00294-005-0044-z.

    Article  CAS  PubMed  Google Scholar 

  52. Luo Q, Zhu Z, Yang R, Qian F, Yan X, Chen H. Characterization of a respiratory burst oxidase homologue from Pyropia haitanensis with unique molecular phylogeny and rapid stress response. J Appl Phycol. 2015;27:945–55.

    Article  CAS  Google Scholar 

  53. Anderson A, Laohavisit A, Blaby IK, Bombelli P, Howe CJ, Merchant SS, et al. Exploiting algal NADPH oxidase for biophotovoltaic energy. Plant Biotechnol J. 2016;14:22–8.

    Article  CAS  Google Scholar 

  54. Shim E, Lee JW, Park H, Zuccarello GC, Kim GH. Hydrogen peroxide signaling mediates fertilization and post-fertilization development in the red alga Bostrychia moritziana. J Exp Bot. 2021:erab453. https://doi.org/10.1093/jxb/erab453.

  55. Zhang YY, Zhu HY, Zhang Q, Li MY, Yan M, Wang R, et al. Phospholipase D alpha 1 and phosphatidic acid regulate NADPH oxidase activity and production of reactive oxygen species in ABA-mediated stomatal closure in Arabidopsis. Plant Cell. 2009;21:2357–77. https://doi.org/10.1105/tpc.108.062992.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Waite KA, Wallin R, Qualliotinemann D, Mcphail LC. Phosphatidic acid-mediated phosphorylation of the NADPH oxidase component p47-phox- evidence that phosphatidic acid may activate a novel protein kinase. J Biol Chem. 1997;272:15569–78. https://doi.org/10.1074/jbc.272.24.15569.

    Article  CAS  PubMed  Google Scholar 

  57. Sheinin A, Nahum-Levy R, Shavit S, Benveniste M. Specificity of putative partial agonist, 1-aminocyclopropanecarboxylic acid, for rat N-methyl-D-aspartate receptor subunits. Neurosci Lett. 2002;317:77–80. https://doi.org/10.1016/s0304-3940(01)02430-2.

    Article  CAS  PubMed  Google Scholar 

  58. Inanobe A, Furukawa H, Gouaux E. Mechanism of partial agonist action at the NR1 subunit of NMDA receptors. Neuron. 2005;47:71–84. https://doi.org/10.1016/j.neuron.2005.05.022.

    Article  CAS  PubMed  Google Scholar 

  59. Mou W, Kao YT, Michard E, Simon A, Li D, Wudick M, et al. Ethylene-independent signaling by the ethylene precursor ACC in Arabidopsis ovular pollen tube attraction. Nat Commun. 2020;11:4082. https://doi.org/10.1038/s41467-020-17819-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgments

We are grateful to Drs. Katsutoshi Arai and Takafumi Fujimoto (Hokkaido University, Japan) for kindly providing access to the LightCycler 480 system.

Funding

This study was supported by a grant-in-aid for Young Scientists (19 K15907 to TU) from the Japan Society for the Promotion of Science.

Author information

Authors and Affiliations

Authors

Contributions

TU was responsible for the design of the experiments and interpretation of data. TU, TK, and SK performed the experiments. TU and HM wrote the manuscript. All authors have read and approved the final version of the manuscript.

Corresponding author

Correspondence to Toshiki Uji.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

The list of primers used for gene expression analysis by quantitative Real Time PCR. Table S2. The list of tested genes for ACC response in N. yezoensis.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Uji, T., Kandori, T., Konishi, S. et al. Phospholipase D activation is required for 1-aminocyclopropane 1-carboxylic acid signaling during sexual reproduction in the marine red alga Neopyropia yezoensis (Rhodophyta). BMC Plant Biol 22, 181 (2022). https://doi.org/10.1186/s12870-022-03575-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-022-03575-z

Keywords