Abstract
Plants exhibit different physiological and molecular responses to adverse changes in their environment. One such molecular response is the sequestration of proteins, RNAs, and metabolites into cytoplasmic bodies called stress granules (cSGs). Here we report that, in addition to cSGs, heat stress also induces the formation of SG-like foci (cGs) in the chloroplasts of the model plant Arabidopsis thaliana. Similarly to the cSGs, (i) cpSG assemble rapidly in response to stress and disappear when the stress ceases, (ii) cpSG formation is inhibited by treatment with a translation inhibitor (lincomycin), and (iii) cpSG are composed of a stable core and a fluid outer shell. A previously published protocol for cSG extraction was successfully adapted to isolate cpSG, followed by protein, metabolite, and RNA analysis. Analogously to the cSGs, cpSG sequester proteins essential for SG formation, dynamics, and function, also including RNA-binding proteins with prion-like domain, ATPases and chaperones, and the amino acids proline and glutamic acid. However, the most intriguing observation relates to the cpSG localization of proteins, such as a complete magnesium chelatase complex, which is involved in photosynthetic acclimation to stress. These data suggest that cpSG have a role in plant stress tolerance.
Introduction
Cell compartmentalization is a way to separate and organize biochemical reactions in order to increase the efficiency of cellular processes. In addition to the membrane-based compartmentalization characteristic for eukaryotic cells, recent evidence has shown that cells can sequester and organize biomolecules, proteins, nucleic acids, and metabolites in non-membrane compartments using the process of liquid–liquid phase separation (LLPS). LLPS is driven by the collective protein–protein and protein–nucleic acid interactions, leading to the emergence of reversible biological condensates (; ; ). Stress granules (SGs) are an example of an evolutionarily conserved non-membrane compartment formed in response to stress and are composed of cytosolic proteins, mRNAs, and small molecules (; ; ).
The key to SG formation is self-assembly of intrinsically disordered proteins (IDPs). Prion-like domains provide structural flexibility and polymer-like behavior to IDPs, which enables the formation of biological scaffolds (; ). Key IDPs involved in SG assembly also contain RNA binding motifs that contribute to RNA sequestration within SGs (; ; ). Via a network of both direct and indirect interactions, SG protein–RNA scaffolds attract numerous additional proteins (; ) and, as recently shown, also small molecules (). SG-associated proteins can be classified into two classes: conserved and specific (; ; ). Conserved proteins are found across multiple SG isolations and represent functional classes required for SG assembly and dynamics, including the already mentioned IDPs, protein and RNA chaperones, ATPases, elongation initiation factors, and small ribosomal subunits. In contrast, sequestration of the specific proteins is dependent on organism, cell type, and developmental stage. Specific SG proteins include metabolic enzymes, e.g., Arabidopsis rhamnose synthase () and regulators, such as kinases (), transcription factors (), and receptors ().
The primary function postulated for SGs is the protection of mRNA and proteins from unfolding and degradation (; ). Sequestration of signaling and regulatory proteins also rapidly inhibit their activity, with consequences for the cellular processes that depend on them. To illustrate this point, in yeast, TORC1, a key regulator of the cell cycle, is sequestered into SGs during heat stress, where it is maintained in an inactive state. During post-stress recovery, SGs disassemble and TORC1 is released, contributing to efficient growth restoration ().
The most important characteristics of SGs are as follows (; ; ): (i) SGs form under stress conditions associated with stalled translation, (ii) SGs sequester both proteins and mRNAs, (iii) SGs are characterized by the presence of conserved proteins involved in SG formation and dynamics, such as IDPs, (iv) SGs rapidly assemble under stress conditions and disassemble when the stress ceases, (v) because of their liquid-like properties, SGs can fuse with each other and also with other biological condensates such as processing bodies and in the process exchange components (), and (vi) SGs are composed of a stable core and more dynamic, fluid shell.
Stress granules are historically associated with the cytosol; however, SG-like structures (chSGs) were also reported in the chloroplast of Chlamydomonas reinhardtii, a single-cell green alga (). Composed of both proteins and mRNAs, chSGs formed in response to oxidative stress and disassembled during stress recovery. Interestingly, the key protein for chSG formation and mRNA sequestration proved to be a large subunit of Rubisco. Moreover, the authors showed that mRNAs are not permanently trapped within chSGs but rather cycle between chSGs and polysomes, a mechanism contributing to the regulation of plastidial translation under stress conditions. Herein and, to our knowledge, for the first time, we report that, in response to heat, SG-like foci also form in the chloroplasts of higher plants. Plastidial SGs (cpSG) sequester proteins, mRNAs, and metabolites, including key photosynthetic regulators, indicative of the importance of cpSG for Arabidopsis stress response.
Materials and Methods
Preparation of 35S:SCO1-GFP Arabidopsis Lines
Arabidopsis thaliana Col-0 plants were transformed with 35S:SCO1-GFP (SCO1-corresponds to snowy cotyledon 1, nuclear-encoded protein localized in chloroplast, known to regulate translation in chloroplasts) provided by Dr. Chanhong Kim using a previously described () floral dip method. T0 seeds were screened for BASTA resistance. The green fluorescent protein (GFP) signal in T1 generation was confirmed using confocal microscopy and is in agreement with previously reported SCO1 cellular localization (). All the presented experiments were performed on the pool of seeds obtained from multiple transformation events.
Plant Growth Conditions and Heat Stress Treatment
Arabidopsis thaliana 35S:SCO1–GFP transgenic seedlings (1.5 mg seeds, c. 100 plants) were grown in sterile liquid cultures (100-ml Erlenmeyer glass flasks) in 35 ml half-strength Murashige and Skoog (1/2MS) liquid medium () with 1% sucrose. The seedlings were kept shaking for 7 days on orbital shakers with the speed of 130 rpm, in constant light (c. 80 μmol m–2 s–1) and temperature (22°C) conditions. After 1 week, the medium in all flasks was changed. After 10 days, heat stress treatment was performed by submitting the flasks to 42°C for 30 min in darkness, which aimed to induce cpSG formation. The control seedlings were kept at room temperature (RT) and in the light. cpSG formation was evaluated by confocal microscopy.
Plastidial Stress Granule Isolation Protocol
The plastidial stress granule extraction protocol was modified from . Extraction was performed in triplicate, starting from the same pool of plant material. Then, 8 g of pulverized material from control or stress (dark/heat) plant material was homogenized using a pre-cooled mortar and pestle with 8 ml of lysis buffer [50 mM Tris–HCl pH 7.4, 100 mM potassium acetate, 2 mM magnesium acetate, 0.5% NP-40, 0.5 mM DTT, 1 mM NaF, 1 mM Na3VO4, protease inhibitor cocktail (Sigma P9599) and 1 U μl–1 RNasin (Promega N2615)]. The resulting slurry was centrifuged at 4,000 g for 10 min at 4°C, the supernatant was removed, and the pellet was resuspended in 6 ml of lysis buffer. The obtained suspension was divided into three to four samples. These were centrifuged at 14,000 g for 10 min at 4°C. The pellets were resuspended in 1 ml lysis buffer, vortexed, and centrifuged at 14,000 g at 4°C for 10 min. As in the previous step, the supernatant was discarded, and the pellets were resuspended in 1 ml of lysis buffer. After a final centrifugation at 850 g for 10 min at 4°C, the supernatant was checked for the presence of cpSG by confocal microscopy and used in the further purification steps. A 1-ml volume (120 μl per sample) of Dynabeads protein A (Thermo Fisher 10002D, Hennigsdorf, Germany) was equilibrated for 30 min at RT with 1 ml of diethyl pyrocarbonate (DEPC)-treated phosphate-buffered saline (PBS) buffer (on a rotating wheel). The beads were separated with a magnet and washed once with DEPC-treated PBS buffer for 5 min at 4°C and then three times with the lysis buffer for 5 min at 4°C. A 370-μl volume of the cpSG fraction supplemented with RNasin (1:100 final dilution) was incubated (on a rotating wheel) with 60 μl of equilibrated Dynabeads for 15 min at room temperature. The beads were then separated, and the supernatant was incubated (on a rotating wheel) in a new tube with 12.5 μl anti-GFP rabbit IgG antibody (Life Technologies A11122, Hennigsdorf, Germany) for 60 min at RT. The excess anti-GFP antibody was separated by centrifugation at 14,000 g for 15 min at 4°C. The supernatant was removed, and the pellet was resuspended in 500 μl of lysis buffer supplemented with an additional 5 μl of RNasin. This resuspension was mixed (on a rotating wheel) with 60 μl of equilibrated Dynabeads for 15 min at RT, followed by several washing steps: three times for 5 min with lysis buffer at 4°C. In the final step, the supernatant was removed, and the SGs remained attached to the beads. Attached to the beads, protein and metabolites were extracted using a methyl tert-butyl ether/methanol/water (MTBE) solvent protocol described previously ().
GAPCP1-GFP Foci Extraction Protocol
GAPCP1-GFP (glyceraldehyde-3-phosphate dehydrogenase, enzyme of the glycolytic pathway) foci extraction was performed following the protocol for stress granules extraction described above. The only difference was the amount of starting material in which 3 g (instead of 8 g) of ground Arabidopsis (without stress) material was used. Proteins and metabolites were extracted according to .
Metabolite and Protein Extraction
Proteins, lipids, and polar compounds were extracted from the beads using the MTBE method () known to separate molecules into a pellet (proteins), an organic phase (lipid compounds), and an aqueous phase (polar compounds). Proteins and metabolites were dried for 4 h up to overnight, respectively, in a centrifugal evaporator and stored at −80°C until liquid chromatography–mass spectrometry analysis.
Protein Preparation
Protein pellets with beads were resuspended in 50 μl of denaturation buffer (6 M urea and 2 M thiourea in 40 mM ammonium bicarbonate). Cysteine residues were reduced by the addition of 2.5 μl of 100 μM DTT for 30 min at RT, followed by alkylation with 2.5 μl of 300 mM iodoacetamide for 20 min in the dark. The proteins were enzymatically digested using LysC/Trypsin Mix (Promega) according to the manufacturer’s instructions. Next, peptides were desalted on C18 SepPack columns and dried using a centrifugal evaporator to about 2 to 3 μl. Dried peptides were stored at −80°C until measurement.
Proteomics Measurements and Data Analysis of SCO1-GFP
Dried peptides were solubilized in a loading buffer (3% ACN, 0.1% FA) and measured on a Q Exactive HF coupled to an Easy nLC1000 UHPLC reverse-phase nano-liquid chromatography (Thermo Fisher Scientific, Hennigsdorf, Germany). The gradient ramped from 4 to 24% ACN over 20 min and then to 45% ACN over the next 15 min, followed by a 5-min washout with 80% ACN. The MS was run using a data-dependent MS/MS method with the following settings: full scans were acquired at a resolution of 120,000, AGC target of 3e6, maximum injection time of 50 ms, and an m/z ranging from 200 to 2,000. dd-MS2 scan was recorded at the 15,000 resolution with an AGC target of 1e5, maximum injection time of 150 ms, isolation window of 1.2 m/z, normalized collision energy 30, and dynamic exclusion of 20 s. Raw chromatograms were processed using the MaxQuant (Version 1.6.0.16, MPI of Biochemistry, Germany) software using the Arabidopsis TAIR database (Version 10, The Arabidopsis Information Resource, www.Arabidopsis.org).
Proteomics Measurements and Data Analysis of GAPCP-GFP
Dried peptides were solubilized in a loading buffer (3% ACN, 0.1% FA) and measured on a Q Exactive HF (Thermo Fisher Scientific, Hennigsdorf, Germany) coupled to a reverse-phase nano-liquid chromatography ACQUITY UPLC M-Class system (Waters). The gradient ramped from 5 to 24% ACN over 7 min and then to 36% ACN over the next 6 min, followed by a 2-min washout with 75% ACN. The MS was run using a data-dependent MS/MS method with the following settings: full scans were acquired at a resolution of 60,000, AGC target of 3e6, maximum injection time of 50 ms, and an m/z ranging from 300 to 1,600. dd-MS2 scan was recorded at 30,000 resolution with an AGC target of 1e5, maximum injection time of 100 ms, isolation window of 1.4 m/z, normalized collision energy 27, and dynamic exclusion of 30 s. Raw chromatograms were processed using the MaxQuant (Version 1.6.0.16, MPI of Biochemistry, Germany) software using the Arabidopsis TAIR database (Version 10, The Arabidopsis Information Resource, www.Arabidopsis.org).
Data Availability
The mass spectrometry proteomics data for SCO1-GFP and GAPCP-GFP were deposited to the ProteomeXchange Consortium via the PRIDE () partner repository with the table identifier PXD018348.
Metabolomic Analysis and Metabolite Annotation
Small molecules were separated by ultraperformance liquid chromatography and analyzed on an Exactive Orbitrap MS (Thermo Fisher Scientific) in positive and negative ionization modes, as described previously (). Data processing, including peak detection and integration and removal of isotopic peaks and chemical noise, was performed using Refiner MS 7.5 (GeneData). An in-house database of chemical compounds was used to annotate the obtained metabolic features (m/z at a given retention time), allowing 10 ppm and 0.1 min deviations from the reference compound mass and retention time, respectively.
Cloning of Plastidial Proteins and Transient Infiltration
Selected genes FNR, PORC, CYP-20, and CHLI2 were cloned into the vector pGWB554 () to overexpress RFP-tagged fusion protein, and Agrobacterium tumefaciens was transformed. Next, A. tumefaciens strain GV3101 harboring the construct of plastidial genes was used to transiently transform () leaves of 3-week-old rossettas expressing plastidial stress granule marker SCO1-GFP. All transient expression assays were performed on the independently grown plants. For each transformation event, three leaves from five independent plants were transformed. Then, 30 min of 42°C dark treatment was used to induce stress granule formation, which was then observed under a confocal microscope.
Confocal Microscopy
To validate the cellular localization of SCO1-GFP, images were acquired using a Leica DM6000B/SP5 confocal laser scanning microscope (Leica Microsystems, Wetzlar, Germany) on leaves from seedlings growing on the liquid medium as mentioned above. The SCO1-GFP signal was visualized using excitation of 488 nm laser and emission between 500 and 520 nm. For the transient infiltration experiment, co-localization images of SCO1-GFP and plastidial proteins were taken using sequential scan mode between lines, allowing separation for the excitation of two fluorophores at once. RFP signal was excited with 561-nm laser, and emission was detected between 595 and 620 nm.
1,6-Hexanediol Treatment
cpSG were extracted using centrifugation steps from the SG extraction protocol. Prepared cpSG-enriched lysate was treated with 10% of 1,6-hexanediol and immediately observed under a confocal microscope. A total of 15 t-stacks were collected using 20 s of intervals for 5 min of total time as suggested in . Measurement of fluorescent intensity per selected cpSG area was performed using Fiji () and region of interest option.
RNA Extraction From Purified cpSG
cpSG were purified using the SG isolation protocol. RNA was extracted from two technical replicates of cpSG isolated from control (RT) and heat-treated (42°C) samples as follows: to each tube of Dynabeads-purified cpSG, 100 μl of nuclease-free water was added, followed by an equal volume of acid phenol/chloroform/isoamyl alcohol (124:25:1, Invitrogen). The sample was vortexed for 20 s and centrifuged at 21,000 g for 5 min at room temperature. The top layer of this sample was added to a new tube, followed by 100 μl of chloroform. The sample was vortexed for 20 s and centrifuged at 21,000 g for 2 min at room temperature. The top layer (∼80 μl) of this sample was added to a new tube. RNA was precipitated by adding 1 μl of 20 μg/μl glycogen (Invitrogen), 10 μl of 3 M sodium acetate, pH 5.5, and 250 μl of 100% ethanol and incubating at −80°C for 4 h. The sample was centrifuged at 21,000 g for 30 min at 4°C, followed by two washes with 70% ethanol and centrifugation at 21,000 g for 10 min at 4°C. Pellet was resuspended in 7.5 μl of nuclease-free water and stored on ice until library preparation.
Library Preparation and RNA Sequencing of RNA From Purified cpSG
RNA sequencing libraries were prepared using the YourSeq Duet (FT + ’-DGE) RNAseq library kit (Amaryllis Nucleics). Input material was RNA pellets (7.5 μl) from cpSG, which were not quantified beforehand due to the low amount of RNA expected to be present. The full-transcript library preparation protocol was used as described by the manufacturer without modification, except that 19 cycles of PCR were used to amplify the libraries instead of the 14 cycles described in the protocol. Index #1 (ATGATTGA) and index #3 (GCTATTCT) were added to control cpSG RNA samples, and index #2 (GACTGCCT) and index #4 (AGCGTTAC) were added to heat-treated cpSG RNA samples. Prior to sequencing, the prepared libraries were quantified with a Qubit fluorometer (Thermo Fisher) and the Qubit dsDNA HS assay kit to ensure that a sufficient quantity of library (>10 μl of 2 ng/μl per sample) was present. Libraries were pooled and sequenced using the MiSeq Nano kit (Illumina), with 300 bp of paired-end reads (2 × 150 bp) generated. Adapter trimmed reads were mapped to the Arabidopsis genome (Araport11; ) using the RMTA workflow (v2.6.3; ) with HiSat2 and 0 FPKM filter options.
Results
Heat Stress Induces the Formation of Plastidial SGs
The protein composition of SG cores was first described in mammalian and yeast cells, using an SG isolation protocol combining differential centrifugation and affinity purification against SG marker proteins (). This protocol was recently adapted to Arabidopsis, yielding a list of 118 cytosolic proteins sequestered within SGs during heat stress (). Intriguingly and in addition to the cytosolic proteins, the SG isolates contained 28 proteins with a reported plastidial localization. Since the liquid phase condensates such as SGs were shown to fuse and, in the process, exchange components (), we decided to explore the possibility of SG-like foci forming in plastids. During lysate preparation, such hypothetical plastidial SGs (cpSG) would fuse with cytosolic SGs (cSGs). Supporting such a scenario, the list of 28 plastidial proteins contained RNA-binding proteins, ATPase, chaperones, and translation elongation factors functionally related to the known SG components (Supplementary Figure S1 and Supplementary Table S1).
To examine the possibility of SGs forming in chloroplasts of higher plants, we expressed one of the 28 plastidial proteins identified in the cSGs, snowy cotyledon 1 (SCO1), as a fusion with GFP under the control of the constitutive cauliflower mosaic virus (CaMV) 35S promoter. This is on par with previous studies, where transgenic lines expressing cSG marker proteins, under the control of the constitutive promoter, were used to characterize the composition of the cSGs (; ). SCO1 is a translation elongation factor (), and interestingly, it shares homology with the translation elongation factor ELF5A-2 found to be present in the cSGs (). In vitro-grown 5-day-old 35S:aSCO1-GFP Arabidopsis seedlings were subjected to control and heat stress conditions (42°C, darkness, 30 min). While under control conditions, SCO1-GFP signal was detected uniformly in the plastidial stroma, in agreement with the previously published SCO1 localization (); in response to heat, SCO1-GFP fluorescence turned into a dotted, star-like pattern reminiscent of the cSGs (Figure 1A). Importantly, in contrast to SCO1-GFP, the GFP protein expressed under the control of the 35S promoter showed no sign of aggregation under our experimental conditions ().
FIGURE 1
We next decided to investigate whether, in addition to the heat-induced assembly, SCO1-GFP foci share other properties of the cSGs. Firstly, as cSGs are characterized by rapid disassembly when the stress ceases, the fate of the SCO1-GFP foci was examined during heat recovery. After subjecting the seedlings to heat stress, the plants were moved back to the control conditions, and SCO1-GFP localization was followed in a time-course manner. SCO1-GFP1 foci formed within 30 min of heat stress and disassembled 40 min into stress recovery (Figure 1A). Secondly, the formation of cSGs is abolished by cycloheximide, a small-molecule inhibitor of protein translation. Analogously, we decided to evaluate whether ccpSG formation might be abolished by the application of lincomycin (Linc), a translation inhibitor specific for plastidial machinery, known to disturb the formation of the first peptide bond from newly synthesized polypeptides (; ). Treatment with Linc alone did not induce cpSG formation, suggesting that disruption of the translation machinery at the early stages of translation does not cause mRNA aggregation. However, incubation with Linc prior to the heat stress prevented the assembly of the SCO1-GFP foci (Figure 1B). This made us conclude that, even though the mechanism of action for Linc is different to cycloheximide, disruption of peptide extension at the initiation stage is enough to prevent mRNA sequestration into cpSG under heat stress conditions. Next, to determine whether, similarly to cSGs, SCO1-GFP foci are composed of a fluid shell surrounding a stable core (), we used a previously described treatment with 10% 1,6-hexanediol (), which interferes with the weak hydrophobic protein–protein and protein–RNA interactions characteristic for the SGs shell (; ; ). Directly after the addition of 1,6-hexanediol to the lysate prepared from the heat-treated SCO1-GFP seedlings, the SCO1-GFP foci were examined with a confocal microscope by a collection of time frames for 5 min, as described in . The 1,6-hexanediol treatment decreased the size of SCO1-GFP foci by approximately 30% (Figures 2A,B and Supplementary Table S2). These results demonstrated that SCO1-GFP foci are composed of a 1,6-hexanediol-sensitive shell and a 1,6-hexanediol-resistant and stable core. Finally, to determine what threshold temperature is required for the assembly of SCO1-GFP foci, SCO1-GFP seedlings were exposed to different temperature and time treatments. While we observed no SCO1-GFP foci formation at 27 and 30°C, 30 min of 33°C treatment was sufficient to induce the appearance of SCO-1-GFP foci (Figure 2C and Supplementary Figure S2). Similarly, cSGs were shown to form in temperatures above 34°C (). In summary, these results support the notion of SCO1-GFP foci having SG-like properties. From this point, we will refer to the SCO1-GFP foci as plastidial SGs (cpSG).
FIGURE 2
Protein and Metabolite Composition of the cpSG
In order to demarcate the protein and metabolite composition of the cpSG, we used SCO1-GFP seedlings grown in a liquid medium and subjected to either control conditions or a combination of heat and dark treatment (30 min, 42°C, and darkness). Following a previously published cSG isolation protocol (; Figure 3A), cpSG were enriched using a series of differential centrifugations, yielding a SG-enriched lysate that was used as input for affinity purification (AP) with anti-GFP antibodies and dynabeads. As observed before, under control conditions, SCO1-GFP was localized uniformly in the plastidial stroma, whereas under stress conditions, SCO1-GFP was sequestered into cpSG (Figure 3B). Importantly, no SCO1-GFP foci could be observed in the SG-enriched lysate from the control plants (Figure 3B). Purified cpSG, also referred to as cpSG isolates, were subjected to mass spectrometry (MS)-based proteomics and metabolomics analysis. A separate isolation was performed to search for cpSG-associated RNAs by extraction of all RNAs from purified cpSG followed by Illumina RNA-sequencing. As an additional negative control, we used seedlings over-expressing plastidial glyceraldehyde-3-phosphate dehydrogenase (35S:GAPCP1-GFP) (). Similarly to SCO1-GFP, GAPCP1-GFP assembles into dot-like foci (Supplementary Figure S3), but unlike the SCO1-GFP, GAPCP1-GFP foci are already present under optimal conditions. GAPCP1-GFP foci were isolated using the SG isolation protocol, with 35S:GFP lines used as a control to exclude non-specific interactors.
FIGURE 3
Multiple criteria were used to define cpSG-associated proteins and metabolites (Supplementary Tables S3–S5). As the cpSG were only present in the heat/dark-treated seedlings, we first selected for proteins and metabolites that were either (i) present in the cpSG isolate of stress-treated seedlings and absent in the cpSG isolate of control seedlings or (ii) significantly enriched (fold change > 2; t-test, p-value < 0.05) in the cpSG isolate from the stress-treated versus the control seedlings. We then excluded a total of 31 cytosolic, nuclear, mitochondrial, and membrane proteins, as these are unlikely to localize into plastidial SGs and can therefore be viewed as false positives. Subcellular localization was retrieved from the SUBA database (). Notably, the exclusion list contained 13 of the 118 proteins previously identified in the cSGs, e.g., polyadenylate-binding protein 8, eukaryotic translation initiation factor 5A-2, and ribonuclease TUDOR-1 (). Such overlap corroborates the previously raised possibility of cpSG and cSGs fusing during lysate preparation. The obtained list comprised 88 proteins and four metabolites (Figure 3C). Notably, of the 88 cpSG proteins, 14 were among the 28 plastidial proteins identified in the cSG isolates (). Considering that the measured plastidial proteome contains approximately 1,500 proteins (), such overlap is 8.5 times more than expected by chance and highly significant (Fisher exact test, p-value = 0.0006). Finally, the list of 88 cpSG proteins was compared with the list of 86 plastid proteins (Supplementary Table S4) associated with the 35S:GAPCP1-GFP condensate (negative list). The comparison revealed an overlap of eight proteins. Considering that the measured plastidial proteome contains approximately 1,500 proteins (), such overlap is again 1.6 times more than expected by chance, but the enrichment is not significant (Fisher exact test, p-value = 1), attesting to the specificity of the cpSG proteome presented here.
The sequence analysis () of the 88 ccpSG proteins revealed the presence of two proteins with prion-like intrinsically disorganized regions associated with SG assembly: CP29A is an RNA-binding protein, and RIP1 (MORF 8) is an RNA editing factor found in plastids and mitochondria. Five of the 88 ccpSG proteins including CP29A have RNA recognition motifs similar to those present in proteins involved in the nucleation of the cSG protein–RNA scaffolds. In fact, CP29 is a paralog of the cSG marker protein RBP47b. Additionally, 15 of the 88 proteins have nucleoside triphosphate hydrolase activity, three are chaperones, 11 constitute ribosomal subunits, and SCO1 and RABE1b are translational elongation factors (Figure 3C).
In addition to proteins which, based on what is known for cSGs, can be linked with cpSG assembly (e.g., CP29A), dynamics (e.g., ATPases), and function (e.g., chaperones) (), the cpSG interactome was characterized by the presence of numerous metabolic enzymes. Chlorophyll biosynthesis clearly stood out (Figure 3C). In addition to all three subunits of the magnesium chelatase complex (ALB1, GUN5, and CHLI1), cpSG sequestered additional enzymes and regulators involved in chlorophyll accumulation, including the magnesium chelatase-associated proteins GUN4 and CHLI2 and an enzyme downstream of the magnesium chelatase complex, PORB. The list of cpSG proteins also contained Rubisco activase and Rubisco accumulation factors required for Rubisco activity. Intriguingly, the enzymes responsible for the different steps of glucosinolate synthesis were found in both cSGs () and cpSG (here).
To confirm our proteomics results, we chose four proteins identified in cpSG: FNR2 (At1g20020), PORC (At1g03630), CYP20-3 (At3g62030), and CHLI2 (At5g45930). The C-terminal/RFP fusions of the four proteins expressed under control of the constitutive CaMV 35S promoter were transiently infiltrated into the leaves of the stable Arabidopsis 35S:SCO1-GFP marker line. To assess whether selected proteins localize into cpSG, infiltrated leaves were treated with 42°C heat stress for 30 min prior to the analysis with confocal microscopy. All four proteins were present in the punctate foci, reminiscent of the cpSG (Figure 4). It is important to note that not all the recorded foci co-localized with the SCO1-GFP. Based on what is known for cSGs (), this is likely reflective of the compositional heterogeneity present in the cpSG population.
FIGURE 4
In addition to proteins, the list of cpSG constituents included four metabolites: two fatty acids, stearic acid and palmitic acid, and two amino acids, glutamic acid and proline. Both proline and glutamic acid were also identified in cSGs (). Moreover, the cpSG localization of all four metabolites can be traced to the cpSG localization of their corresponding protein interactors, e.g., proline was shown to bind to chaperones ().
Finally, the RNA-seq analysis of cpSG revealed the presence of plastidial transcripts (Supplementary Table S6) with highest enrichment measured for ribosomal proteins and subunits of the ATP synthase complex.
Discussion
Stress granules are evolutionarily conserved aggregates of proteins, and untranslated mRNAs formed in response to stress. Here we demonstrate that SG-like foci also assemble in the chloroplasts of higher plants. Our supposition is based on multiple lines of evidence proving the broad similarity of the cpSG reported here with previously identified cSGs (; ; ): (i) cSG formation is a rapid, stress-dependent, and reversible process. Analogously, we demonstrated that cpSG form within 10 min of 42°C heat stress but are absent under control conditions. Moreover, cpSG disassemble in plants transferred to control temperature, 40 min into the heat stress recovery. We also showed that, similarly to cSGs, cpSG are induced by temperatures higher than 33°C (); (ii) Cycloheximide treatment is known to interfere with the formation of cSGs by trapping mRNAs required for the assembly of cpSG/cSGs scaffolds at polysomes (; ). Although Lync has a different mode of action in comparison to cycloheximide, its application prevented cpSG formation under heat stress conditions; (iii) Evaluation of cpSG properties by hexanediol treatment (; ) revealed that, similarly to cSGs, cpSG belong to the so-called liquid–solid demixing phase separation assemblies, which are composed of a stable core and a fluid shell. While a core is composed of a dense network of protein–protein-mRNA interactions, the shell is less concentrated and more dynamic in terms of its protein composition (; ; ).
As demarcated here, the cpSG proteome and metabolome bear a functional resemblance with those reported for cSGs. Remarkably, this resemblance extends not only to plant cSGs (; ; ) but also to SGs in yeast and mammalian cells (; ; ), attesting to the evolutionary conservation of proteins required for SG formation, dynamics, and function. RNA-binding proteins with prion-like domains, homologous to the CP29A reported here, are essential scaffold proteins because of their ability to form polymers and to sequester mRNAs (). Indeed CP29A was shown to bind and regulate the stability of multiple chloroplast mRNAs (). This makes CP29A a key candidate protein for cpSG assembly and function. Moreover, an important role in the regulation of cSG dynamics has been assigned to proteins with ATPase activity (), as ATP hydrolysis is required for stress granule assembly and dynamics. Enzymes with ATPase activity, e.g., DEAD box RNA helicase (RH3), were also found in cpSG. Interestingly, we also identified protoporphyrinogen IX oxidase, an enzyme from the tetrapyrrole biosynthetic pathway which was also reported to play an additional role in RNA editing of plastidial transcripts, leading to the modulation of NADH dehydrogenase-like complex activity (). A detailed analysis of RNA present in cpSG revealed the presence of two, ATCG00905 and ATCG01230, transcripts identified as edited targets by MORF2 (through induction by GUN1) (), which suggests that cpSG might be also formed as foci supporting RNA editing events in chloroplasts in the context of retrograde signaling. Finally, chaperones identified in cpSG, e.g., heat shock protein 90-5, likely function to prevent the misfolding of proteins concentrated in SGs. Along similar lines, proline that accumulates in cSGs and cpSG was shown to activate molecular chaperones and, by doing so, may prevent protein misfolding ().
Our proteomics analysis also revealed the presence of multiple stress-related proteins in Arabidopsis cpSG. The prime function of cSGs is to protect proteins from unfolding and aggregation. We speculate that cpSG may have a similar role. To exemplify this, we will use translation elongation factor Tu (RABE1b), which we detected in purified cpSG. RABE1b was shown to be prone to heat-induced irreversible aggregation both in vitro and in vivo, but while in vitro RABE1b aggregation could be observed even under slightly elevated temperatures, in vivo it required much higher temperatures, above 40°C (). The authors speculated about the existence of a cellular mechanism that would promote RABE1b stability and reduce its aggregation. Such a mechanism would be crucial because RABE1b is essential for plastid translation and heat stress tolerance. We speculate that reversible sequestration of RABE1b into cpSG, as reported here, would confer previously demonstrated in vivo resistance to unfolding and aggregation. Similarly, this could be true for Rubisco activase, also found in cpSG, which is extremely sensitive to aggregation by heat stress ().
In addition to its role in conferring protection, cSG sequestration was also shown to serve as a mechanism to temporarily deactivate regulatory proteins important, e.g., for plant growth (; ; ). We propose that such a mechanism could be relevant for chlorophyll biosynthetic enzymes such as the magnesium chelatase complex (; ; ), which were identified in cpSG. Because of their chemical reactivity, the levels of both chlorophyll and chlorophyll intermediates need to be tightly controlled. Rapid inactivation of chlorophyll synthesis by cpSG sequestration of the key biosynthetic enzymes offers an elegant mechanism to regulate chlorophyll levels in response to changing environmental conditions, such as heat and darkness imposed in our experiment. Once the stress ceases, cpSG disassemble, and chlorophyll synthesis could be quickly reactivated to meet photosynthetic demand.
Similarly to protein, mRNA sequestration within SGs can have a double role. On one hand, it protects mRNAs from degradation and, on the other hand, it can be used as means of translational regulation. For instance, under low-oxygen stress, cSG sequestration of the mRNAs encoding abundant, non-stress transcripts contributes to the preferential translation of stress-induced mRNAs. Upon reoxygenation, cSGs disassemble and the stored mRNAs return to polysomes ().
Here we present a compelling evidence for the existence of cpSG in higher plants. We demonstrated that cpSG sequester a plethora of proteins, metabolites, and mRNAs. By drawing a homology with cSGs, it is fair to speculate about the biological role of cpSG sequestration, on one hand as a protection mechanism and on the other hand as a means of rapid regulation. Future work will focus on (i) understanding the exact mechanism underlying cpSG assembly and disassembly, (ii) assessing compositional differences related to varying stress conditions, and, most importantly, (iii) validation of the biological significance of cpSG sequestration.
Statements
Data availability statement
The mass spectrometry proteomics data for SCO1-GFP and GAPCP-GFP were deposited to the ProteomeXchange Consortium via the PRIDE () partner repository with the table identifier PXD018348.
Author contributions
AS and MC conceived the project and wrote the manuscript. MC, AN-D, AN, JM, AM, and ES performed experimental work. AN-D and AN were involved in review and editing. AS provided resources. AS and MC agreed to serve as the authors responsible for contact and ensure communication. All authors contributed to the article and approved the submitted version.
Funding
The funding for the research included in the manuscript was fully provided by Max Planck Institute of Molecular Plant Physiology.
Acknowledgments
We would like to acknowledge Prof. Alisdair Fernie for critical discussion of our work, Änne Michaelis for technical support (metabolite measurement), and Stella Kössler for assistance with the affinity experiment.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2020.595792/full#supplementary-material
References
1
AlbertiS. (2017). The wisdom of crowds: regulating cell function through condensed states of living matter.J. Cell Sci.1302789–2796. 10.1242/jcs.200295
2
AlbrechtV.IngenfeldA.ApelK. (2006). Characterization of the snowy cotyledon 1 mutant of Arabidopsis thaliana: the impact of chloroplast elongation factor G on chloroplast development and plant vitality.Plant Mol. Biol.60507–518. 10.1007/s11103-005-4921-0
3
ApweilerR.BairochA.WuC. H.BarkerW. C.BoeckmannB.FerroS.et al (2004). UniProt: the Universal Protein knowledgebase.Nucleic Acids Res.32D115–D119. 10.1093/nar/gkh131
4
BhasinH.HulskampM. (2017). ANGUSTIFOLIA, a plant homolog of CtBP/BARS localizes to stress granules and regulates their formation.Front. Plant Sci.8:1004. 10.3389/fpls.2017.01004
5
BuchanJ. R.ParkerR. (2009). Eukaryotic stress granules: the ins and outs of translation.Mol. Cell.36932–941. 10.1016/j.molcel.2009.11.020
6
ChantarachotT.Bailey-SerresJ. (2018). Polysomes, stress granules, and processing bodies: a dynamic triumvirate controlling cytoplasmic mRNA fate and function.Plant Physiol.176254–269. 10.1104/pp.17.01468
7
ChengC. Y.KrishnakumarV.ChanA. P.Thibaud-NissenF.SchobelS.TownC. D. (2017). Araport11: a complete reannotation of the Arabidopsis thaliana reference genome.Plant J.89789–804. 10.1111/tpj.13415
8
ChotewutmontriP.BarkanA. (2018). Multilevel effects of light on ribosome dynamics in chloroplasts program genome-wide and psbA-specific changes in translation.PLoS Genet.14:e1007555. 10.1371/journal.pgen.1007555
9
CourchaineE. M.LuA.NeugebauerK. M. (2016). Droplet organelles?EMBO J.351603–1612. 10.15252/embj.201593517
10
DiamantS.EliahuN.RosenthalD.GoloubinoffP. (2001). Chemical chaperones regulate molecular chaperones in vitro and in cells under combined salt and heat stresses.J. Biol. Chem.27639586–39591. 10.1074/jbc.M103081200
11
GiavaliscoP.LiY.MatthesA.EckhardtA.HubbertenH. M.HesseH.et al (2011). Elemental formula annotation of polar and lipophilic metabolites using (13) C, (15) N and (34) S isotope labelling, in combination with high-resolution mass spectrometry.Plant J.68364–376. 10.1111/j.1365-313X.2011.04682.x
12
GibsonL. C.WillowsR. D.KannangaraC. G.von WettsteinD.HunterC. N. (1995). Magnesium-protoporphyrin chelatase of Rhodobacter sphaeroides: reconstitution of activity by combining the products of the bchH, -I, and -D genes expressed in Escherichia coli.Proc. Natl. Acad. Sci. U.S.A.921941–1944. 10.1073/pnas.92.6.1941
13
GilksN.KedershaN.AyodeleM.ShenL.StoecklinG.DemberL. M.et al (2004). Stress granule assembly is mediated by prion-like aggregation of TIA-1.Mol. Biol. Cell155383–5398. 10.1091/mbc.e04-08-0715
14
GomesE.ShorterJ. (2019). The molecular language of membraneless organelles.J. Biol. Chem.2947115–7127. 10.1074/jbc.TM118.001192
15
Gutierrez-BeltranE.MoschouP. N.SmertenkoA. P.BozhkovP. V. (2015). Tudor staphylococcal nuclease links formation of stress granules and processing bodies with mRNA catabolism in Arabidopsis.Plant Cell27926–943. 10.1105/tpc.114.134494
16
HamadaT.YakoM.MinegishiM.SatoM.KameiY.YanagawaY.et al (2018). Stress granule formation is induced by a threshold temperature rather than a temperature difference in Arabidopsis.J. Cell Sci.131:jcs216051. 10.1242/jcs.216051
17
HerrinD. L.NickelsenJ. (2004). Chloroplast RNA processing and stability.Photosynth. Res.82301–314. 10.1007/s11120-004-2741-8
18
HooperC. M.CastledenI. R.TanzS. K.AryamaneshN.MillarA. H. (2017). SUBA4: the interactive data analysis centre for Arabidopsis subcellular protein locations.Nucleic Acids Res.45D1064–D1074. 10.1093/nar/gkw1041
19
IkegamiA.YoshimuraN.MotohashiK.TakahashiS.RomanoP. G.HisaboriT.et al (2007). The CHLI1 subunit of Arabidopsis thaliana magnesium chelatase is a target protein of the chloroplast thioredoxin.J. Biol. Chem.28219282–19291. 10.1074/jbc.M703324200
20
JainS.WheelerJ. R.WaltersR. W.AgrawalA.BarsicA.ParkerR. (2016). ATPase-modulated stress granules contain a diverse proteome and substructure.Cell164487–498. 10.1016/j.cell.2015.12.038
21
KedershaN. L.GuptaM.LiW.MillerI.AndersonP. (1999). RNA-binding proteins TIA-1 and TIAR link the phosphorylation of eIF-2 alpha to the assembly of mammalian stress granules.J. Cell Biol.1471431–1441. 10.1083/jcb.147.7.1431
22
KosmaczM.GorkaM.SchmidtS.LuzarowskiM.MorenoJ. C.SzlachetkoJ.et al (2019). Protein and metabolite composition of Arabidopsis stress granules.New Phytol.2221420–1433. 10.1111/nph.15690
23
KroschwaldS.MaharanaS.MatejuD.MalinovskaL.NuskeE.PoserI.et al (2015). Promiscuous interactions and protein disaggregases determine the material state of stress-inducible RNP granules.eLife4:e06807. 10.7554/eLife.06807.060
24
KupschC.RuweH.GusewskiS.TillichM.SmallI.Schmitz-LinneweberC. (2012). Arabidopsis chloroplast RNA binding proteins CP31A and CP29A associate with large transcript pools and confer cold stress tolerance by influencing multiple chloroplast RNA processing steps.Plant Cell244266–4280. 10.1105/tpc.112.103002
25
LancasterA. K.Nutter-UphamA.LindquistS.KingO. D. (2014). PLAAC: a web and command-line application to identify proteins with prion-like amino acid composition.Bioinformatics302501–2502. 10.1093/bioinformatics/btu310
26
LeeM. W.YangY. (2006). Transient expression assay by agroinfiltration of leaves.Methods Mol. Biol.323225–229. 10.1385/1-59745-003-0:225
27
LiX.CaiC.WangZ.FanB.ZhuC.ChenZ. (2018). Plastid translation elongation factor Tu is prone to heat-induced aggregation despite its critical role in plant heat tolerance.Plant Physiol.1763027–3045. 10.1104/pp.17.01672
28
LokdarshiA.ConnerW. C.McClintockC.LiT.RobertsD. M. (2016). Arabidopsis CML38, a calcium sensor that localizes to ribonucleoprotein complexes under hypoxia stress.Plant Physiol.1701046–1059. 10.1104/pp.15.01407
29
LorkovicZ. J.BartaA. (2002). Genome analysis: RNA recognition motif (RRM) and K homology (KH) domain RNA-binding proteins from the flowering plant Arabidopsis thaliana.Nucleic Acids Res.30623–635. 10.1093/nar/30.3.623
30
MarkmillerS.SoltaniehS.ServerK. L.MakR.JinW.FangM. Y.et al (2018). Context-dependent and disease-specific diversity in protein interactions within stress granules.Cell172590.e13–604.e13. 10.1016/j.cell.2017.12.032
31
MolliexA.TemirovJ.LeeJ.CoughlinM.KanagarajA. P.KimH. J.et al (2015). Phase separation by low complexity domains promotes stress granule assembly and drives pathological fibrillization.Cell163123–133. 10.1016/j.cell.2015.09.015
32
Munoz-BertomeuJ.Cascales-MinanaB.MuletJ. M.Baroja-FernandezE.Pozueta-RomeroJ.KuhnJ. M.et al (2009). Plastidial glyceraldehyde-3-phosphate dehydrogenase deficiency leads to altered root development and affects the sugar and amino acid balance in Arabidopsis.Plant Physiol.151541–558. 10.1104/pp.109.143701
33
MurashigeT.SkoogF. (1962). A revised medium for rapid growth and bio assays with tobacco tissue cultures.Physiol. Plant.15473–497. 10.1111/j.1399-3054.1962.tb08052.x
34
NakagawaT.SuzukiT.MurataS.NakamuraS.HinoT.MaeoK.et al (2007). Improved Gateway binary vectors: high-performance vectors for creation of fusion constructs in transgenic analysis of plants.Biosci. Biotechnol. Biochem.712095–2100. 10.1271/bbb.70216
35
PapenbrockJ.MockH. P.TanakaR.KruseE.GrimmB. (2000). Role of magnesium chelatase activity in the early steps of the tetrapyrrole biosynthetic pathway.Plant Physiol.1221161–1169. 10.1104/pp.122.4.1161
36
PatelS. S.BelmontB. J.SanteJ. M.RexachM. F. (2007). Natively unfolded nucleoporins gate protein diffusion across the nuclear pore complex.Cell12983–96. 10.1016/j.cell.2007.01.044
37
PeriS.RobertsS.KrekoI. R.McHanL. B.NaronA.RamA.et al (2019). Read mapping and transcript assembly: a scalable and high-throughput workflow for the processing and analysis of ribonucleic acid sequencing data.Front. Genet.10:1361. 10.3389/fgene.2019.01361
38
PeskettT. R.RauF.O’DriscollJ.PataniR.LoweA. R.SaibilH. R. (2018). A liquid to solid phase transition underlying pathological huntingtin exon1 aggregation.Mol. Cell.70588.e6–601.e6. 10.1016/j.molcel.2018.04.007
39
ProtterD. S.ParkerR. (2016). Principles and properties of stress granules.Trends Cell. Biol.26668–679. 10.1016/j.tcb.2016.05.004
40
SalvucciM. E.Crafts-BrandnerS. J. (2004). Inhibition of photosynthesis by heat stress: the activation state of Rubisco as a limiting factor in photosynthesis.Physiol. Plant120179–186. 10.1111/j.0031-9317.2004.0173.x
41
SchindelinJ.Arganda-CarrerasI.FriseE.KaynigV.LongairM.PietzschT.et al (2012). Fiji: an open-source platform for biological-image analysis.Nat. Methods9676–682. 10.1038/nmeth.2019
42
ShannonP.MarkielA.OzierO.BaligaN. S.WangJ. T.RamageD.et al (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks.Genome Res.132498–2504. 10.1101/gr.1239303
43
SorensonR.Bailey-SerresJ. (2014). Selective mRNA sequestration by OLIGOURIDYLATE-BINDING PROTEIN 1 contributes to translational control during hypoxia in Arabidopsis.Proc. Natl. Acad. Sci. U.S.A.1112373–2378. 10.1073/pnas.1314851111
44
SzklarczykD.SantosA.von MeringC.JensenL. J.BorkP.KuhnM. (2016). STITCH 5: augmenting protein-chemical interaction networks with tissue and affinity data.Nucleic Acids Res.44D380–D384. 10.1093/nar/gkv1277
45
TakaharaT.MaedaT. (2012). Transient sequestration of TORC1 into stress granules during heat stress.Mol. Cell.47242–252. 10.1016/j.molcel.2012.05.019
46
UniackeJ.ZergesW. (2008). Stress induces the assembly of RNA granules in the chloroplast of Chlamydomonas reinhardtii.J. Cell Biol.182641–646. 10.1083/jcb.200805125
47
VizcainoJ. A.CsordasA.Del-ToroN.DianesJ. A.GrissJ.LavidasI.et al (2016). 2016 update of the PRIDE database and its related tools.Nucleic Acids Res.44:11033. 10.1093/nar/gkw880
48
WeberC.NoverL.FauthM. (2008). Plant stress granules and mRNA processing bodies are distinct from heat stress granules.Plant J.56517–530. 10.1111/j.1365-313X.2008.03623.x
49
WheelerJ. R.JainS.KhongA.ParkerR. (2017). Isolation of yeast and mammalian stress granule cores.Methods12612–17. 10.1016/j.ymeth.2017.04.020
50
WheelerJ. R.MathenyT.JainS.AbrischR.ParkerR. (2016). Distinct stages in stress granule assembly and disassembly.eLife5:e18413. 10.7554/eLife.18413.018
51
YounJ. Y.DunhamW. H.HongS. J.KnightJ. D. R.BashkurovM.ChenG. I.et al (2018). High-density proximity mapping reveals the subcellular organization of mRNA-associated granules and bodies.Mol. Cell.69517.e1–532.e1. 10.1016/j.molcel.2017.12.020
52
YounJ. Y.DyakovB. J. A.ZhangJ.KnightJ. D. R.VernonR. M.Forman-KayJ. D.et al (2019). Properties of stress granule and P-body proteomes.Mol. Cell.76286–294. 10.1016/j.molcel.2019.09.014
53
ZhangF.TangW.HedtkeB.ZhongL.LiuL.PengL.et al (2014). Tetrapyrrole biosynthetic enzyme protoporphyrinogen IX oxidase 1 is required for plastid RNA editing.Proc. Natl. Acad. Sci. U.S.A.1112023–2028. 10.1073/pnas.1316183111
54
ZhangX.HenriquesR.LinS. S.NiuQ. W.ChuaN. H. (2006). Agrobacterium-mediated transformation of Arabidopsis thaliana using the floral dip method.Nat. Protoc.1641–646. 10.1038/nprot.2006.97
55
ZhaoX.HuangJ.ChoryJ. (2019). GUN1 interacts with MORF2 to regulate plastid RNA editing during retrograde signaling.Proc. Natl. Acad. Sci. U.S.A.11610162–10167. 10.1073/pnas.1820426116
56
ZybailovB.RutschowH.FrisoG.RudellaA.EmanuelssonO.SunQ.et al (2008). Sorting signals, N-terminal modifications and abundance of the chloroplast proteome.PLoS One3:e1994. 10.1371/journal.pone.0001994
Summary
Keywords
stress granules, plastid, protome, chlorophyll biosynthesis, metabolites
Citation
Chodasiewicz M, Sokolowska EM, Nelson-Dittrich AC, Masiuk A, Beltran JCM, Nelson ADL and Skirycz A (2020) Identification and Characterization of the Heat-Induced Plastidial Stress Granules Reveal New Insight Into Arabidopsis Stress Response. Front. Plant Sci. 11:595792. doi: 10.3389/fpls.2020.595792
Received
17 August 2020
Accepted
07 October 2020
Published
30 October 2020
Volume
11 - 2020
Edited by
Adriano Nunes-Nesi, Universidade Federal de Viçosa, Brazil
Reviewed by
Hunseung Kang, Chonnam National University, South Korea; Robert Drant Willows, Macquarie University, Australia
Updates
Copyright
© 2020 Chodasiewicz, Sokolowska, Nelson-Dittrich, Masiuk, Moreno Beltran, Nelson and Skirycz.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Monika Chodasiewicz, kosmacz@mpimp-golm.mpg.de; Chodasiewicz@mpimp-golm.mpg.deAleksandra Skirycz, skirycz@mpimp-golm.mpg.de
This article was submitted to Plant Physiology, a section of the journal Frontiers in Plant Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.