Salinity-Induced Palmella Formation Mechanism in Halotolerant Algae Dunaliella salina Revealed by Quantitative Proteomics and Phosphoproteomics

Abstract

Palmella stage is critical for some unicellular algae to survive in extreme environments. The halotolerant algae Dunaliella salina is a good single-cell model for studying plant adaptation to high salinity. To investigate the molecular adaptation mechanism in salinity shock-induced palmella formation, we performed a comprehensive physiological, proteomics and phosphoproteomics study upon palmella formation of D. salina using dimethyl labeling and Ti⁴⁺-immobilized metal ion affinity chromatography (IMAC) proteomic approaches. We found that 151 salinity-responsive proteins and 35 salinity-responsive phosphoproteins were involved in multiple signaling and metabolic pathways upon palmella formation. Taken together with photosynthetic parameters and enzyme activity analyses, the patterns of protein accumulation and phosphorylation level exhibited the mechanisms upon palmella formation, including dynamics of cytoskeleton and cell membrane curvature, accumulation and transport of exopolysaccharides, photosynthesis and energy supplying (i.e., photosystem II stability and activity, cyclic electron transport, and C4 pathway), nuclear/chloroplastic gene expression regulation and protein processing, reactive oxygen species homeostasis, and salt signaling transduction. The salinity-responsive protein–protein interaction (PPI) networks implied that signaling and protein synthesis and fate are crucial for modulation of these processes. Importantly, the 3D structure of phosphoprotein clearly indicated that the phosphorylation sites of eight proteins were localized in the region of function domain.

Introduction

The unicellular algae can develop a vegetative palmella in their life cycle, when exposed to various extreme environment conditions, such as salinity (Takouridis et al., 2015), copper (Sztrum et al., 2012), organic acids (Iwasa and Murakami, 1968), herbicide (Franqueira et al., 2000), oxidative stress (Wang et al., 2004a), and predators (Lurling and Beekman, 2006). The palmella formation in Chlamydomonas reinhardtii (Iwasa and Murakami, 1968; Lurling and Beekman, 2006; Sztrum et al., 2012) and Pediastrum tetras (Ellis, 1972) have been found to be regulated by free calcium concentration, starch accumulation, carotenoid/chlorophyll ratio, as well as the activities of catalase (CAT) and ascorbate peroxidase (APX). Similarly, for green motile flagellates of Haematococcus pluvialis, oxidative stress could induce its green motile cells transforming into enlarged red resting cysts (aplanospores) (Boussiba et al., 1999; Wang et al., 2004a; Han et al., 2012). In H. pluvialis, proteomics study revealed that the abundances of several proteins were induced, which were involved in reactive oxygen species (ROS) scavenging [e.g., superoxide dismutase (SOD), CAT and peroxidase (POD)], photosynthesis [e.g., ribulose-1,5-bisphosphate carboxylase/oxygenase large subunit (Rubisco LSU), phosphoglycerate kinase (PGK)], nitrogen assimilation [e.g., glutamine synthetase (GS)], and mitochondrial respiratory (e.g., mitochondrial ATPase β subunit) (Wang et al., 2004a). Moreover, the reduction of chlorophyll content, moderate declines in the maximal photosynthetic rate and the maximum quantum yield of photosystem (PS) II, as well as the significant increase in PS I activity were detected during the transformation of green vegetative cells to red aplanospores in H. pluvialis (Han et al., 2012). These indicates that the early stress response involves multiple enzymatic defense processes, which plays a critical role during the transition of green vegetative cells to red cysts (Wang et al., 2004a).

The halotolerant unicellular green algae Dunaliella salina is a model photosynthetic organism for studying plant adaptation to high salinity (Katz et al., 2007), which can adapt to a wide range of salinities ranging from about 0.05 M NaCl to saturation (around 5.5 M NaCl) (Liska et al., 2004). During their life cycle, D. salina may have a dominant palmella stage formed by a colony-like group of round non-motile cells to cope with extreme conditions (Borowitzka and Siva, 2007). When entering palmella stage, the cells usually lose their flagella and eyespot, become round, and excrete a slime exopolysaccharides (EPSs) layer outside the accumulation of green cells, but when returning to fresh medium at a “normal” salinity, algae cells usually reform their flagella and return to the motile and free-swimming state (Borowitzka and Siva, 2007).

It is reported that high salinity triggers palmella formation in D. salina (Montoya and Olivera, 1993), but the molecular mechanism of palmella formation is still unknown. Previous studies mainly focus on the mechanism of remarkable salinity and osmotic adaptability of D. salina. It is found that the de novo synthesis of compatible solute (i.e., glycerol and β-carotene) (Ben-Amotz and Avron, 1973), higher activity of Na⁺ extrusion systems (Pick, 2002), as well as active photosynthetic and energy metabolism (Liska et al., 2004; Katz et al., 2007; Alkayal et al., 2010) are crucial for salinity-response in D. salina. Importantly, proteomic research has provided new insights into the high salinity-responsive strategies in D. salina. The abundances of several proteins were salinity-regulated, which were involved in Calvin cycle, starch mobilization, energy production, protein synthesis/degradation, membrane structure stabilization, and signal transduction in D. salina (Liska et al., 2004; Katz et al., 2007). However, whether these fine-tuned molecular mechanism also happened in the process of salinity-induced palmella formation is poorly understood.

In this study, the salinity shock-induced palmella formation was analyzed in D. salina. By integrative analysis of morphology changes, photosynthesis, antioxidant enzyme activities, as well as protein abundance and phosphorylation level, we highlighted several important mechanisms upon palmella formation, such as photosynthetic modulation, ROS scavenging, gene expression regulation, and protein post-translational modification.

Experimental procedures

Cultivation, treatment, and biomass analysis

D. salina cells were cultivated in a modified medium containing 1 M NaCl under a 16/8 h light/dark cycle (light intensity 100 μmol photons × m⁻² × s⁻¹) at 26°C with shaking at 96 rpm (Katz and Avron, 1985). Cells in the stationary phase were transferred to fresh medium containing 1 M NaCl (control) and 3 M NaCl (salinity shock), respectively. After treatment, algae cells were used freshly or stored at −80°C for experiments.

The cell growth was evaluated by the absorbance of cultures at 630 nm in a spectrophotometer (Chen et al., 2009). The cell number was counted using a haemacytometer under light microscope. Three independent biological replicates for each sample were conducted for all the experiments. The morphology of D. salina cells were observed under Olympus BX53 Microscope (Olympus America Inc., Center Valley, PA, USA) equipped with Olympus DP72 digital camera system (Tokyo, Japan).

Assessment of cell viability

The cell viability was estimated according to the method of Mendes et al. (2013). 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) stock solution (60 μl; 5 mg/ml) was applied to 250 μl cells for 60 min at 26°C with shaking at 96 rpm, and then the samples were centrifuged at 8,000 × g for 5 min. The formazan crystals were dissolved in 150 μl of dimethylsulfoxide and vortexed for 10 min. After centrifugation at 12,000 × g for 5 min, the supernatants were transferred to 96-well cell culture plate. The absorbance at 490 nm was measured using an iMark™ Microplate Reader (Bio-Rad, Richmond, CA, USA).

Carbohydrate and glycerol content analysis

Carbohydrates were analyzed according to the method of Li et al. (2000). Ten milliliter culture was centrifuged at 12,000 × g for 10 min, and then 2 ml of supernatant was analyzed by the anthrone reagent method (Pistocchi et al., 1997). Glycerol content was determined according to the method of Chen et al. (2011).

Chlorophyll content, photosynthesis oxygen evolution, respiration rate, chlorophyll fluorescence and P700 analysis

The chlorophyll contents were determined using a method described by Lichtenthaler and Wellburn (1983). Photosynthesis oxygen evolution and respiration rate were recorded by Clark-type oxygen electrode (Hansatech, UK) according to the method of Ma et al. (2005). The chlorophyll fluorescence parameters and P700 were detected using a pulse-amplitudemodulated chlorophyll fluorometer (Dual-PAM-100; Walz, Effeltrich, Germany) and an emitter-detector-cuvette assembly with a unit 101ED (ED-101US) according to the method of Ma et al. (2008).

Enzyme activity assay

Antioxidant enzyme assays were performed essentially according to Suo et al. (2015). For the sample preparation of SOD, CAT, POD, APX, monodehydroascorbate reductase (MDHAR), dehydroascorbate reductase (DHAR), glutathione reductase (GR), glutathione S-transferase (GST), and glutathione peroxidase (GPX), 0.5 g of cells was homogenized on ice in 3 ml buffer [50 mM phosphate buffer (pH 7.8), 2% PVP-40, and 2 mM ascorbate (AsA)]. After centrifugation at 15,000 × g for 20 min at 4°C, the supernatants were used for enzyme activity assays according to the method of Zhao et al. (2016).

For glycerol-3-phosphate dehydrogenase (GPDH) activity assay, cells were ground to a fine powder in liquid nitrogen and resuspended in extraction buffer (100 mM Tris, 20 mM ascorbic acid, pH 6.9) (Chen et al., 2009). After being centrifuged at 34,000 × g for 20 min at 4°C, the supernatant was collected for GPDH activity assay. GPDH activity was determined using a plant GPDH assay kit according to manufacturer's instructions (IBL Hamburg, Germany).

Phosphoglucomutase (PGM) activity was determined using a PGM colorimetric assay kit according to manufacturer's instructions (BioVision, USA). For the H⁺-ATPase activity, 0.5 g of cells was used for the preparation of membranes. Plasmalemma vesicles were prepared as described by Kasamo (1988). The H⁺-ATPase activity was measured according to the method of Liang (1999).

Generation rate of and contents of H₂O₂, AsA, DHA, GSSG, and GSH

The generation rate of was measured by monitoring the nitrite formation from hydroxylamine in the presence of and the absorbance in the aqueous solution was read at 530 nm (Zhao et al., 2016). The content of H₂O₂ was measured by monitoring the A₄₁₀ of titanium-peroxide complex according to the method of Suo et al. (2015).

For the contents of reduced AsA, oxidized AsA (DHA), oxidized glutathione (GSSG), and reduced glutathione (GSH), 0.5 g of cells was homogenized on ice in 3 ml 5% trichloroacetic acid. After centrifugation at 15,000 × g for 15 min at 4°C, the supernatants were used for substrate content assays. Total AsA and reduced AsA were determined by recording the absorbance changes at 525 nm (Kampfenkel et al., 1995). DHA content was estimated from the difference between assays with and without dithiothreitol (DTT). The extracts were also used for GSH/GSSG assay according to the method of Baker et al. (1990). The concentration of GSH was calculated from the differences between total glutathione and GSSG assays (Haghjou et al., 2009).

Protein sample preparation

The proteins for differential stable isotope labeling experiments were extracted in the modified extraction buffer (900 mM sucrose, 100 mM Tris-HCl (pH 8.0), 65 mM DTT, 1 mM MgCl₂, 1% Triton X-100 (v/v), 1 mM phenylmethanesulfonyl fluoride, 1% protease inhibitor cocktail (v/v), 1 mM EDTA, 1 mM EDGA, 1 mM NaF, 1 mM Na₃VO₄, 1 mM Na₄O₇P₂, and 1 mM C₃H₇Na₂O₆P) according to the phenol-methanol methods of Suo et al. (2015). Protein concentration was determined using a Quant-kit according to manufacture instructions (GE Healthcare, USA).

Protein in-solution digestion, dimethyl labeling and phosphopeptide enrichment

The protein was resuspended in the denaturing buffer containing 8 M urea and 100 mM triethylammonium bicarbonate (pH 8.0) (Bian et al., 2012; Song et al., 2014), and then was reduced by DTT and alkylated by iodoacetamide. After that, the solution was diluted with 100 mM triethylammonium bicarbonate and trypsin was added with an enzyme-to-protein ratio of 1/25 (w/w) and incubated at 37°C overnight.

For the light and heavy dimethyl labeling, 500 μl of CH₂O (4%, v/v) and CD₂O (4%, v/v) was added to the tryptic digests from D. salina cells cultured in 1 M NaCl (control) or 3 M NaCl (salinity shock for 6 h) medium, respectively, then 500 μl of freshly prepared NaBH₃CN (0.6 M) was subsequently added. The detailed labeling and desalted procedure were the same as described in the method of Song et al. (2014).

The phosphopeptides from 2.0 mg of labeled protein digestion mixture were enriched by Ti⁴⁺-immobilized metal ion affinity chromatography (IMAC) microspheres (Zhou et al., 2008, 2013) and then were resuspended in 30 μl 5% formic acid.

LC-MS/MS analysis

For the reversed-phase (RP) liquid chromatography (LC)-MS/MS analysis, a 12 cm × 75 μm i.d. capillary column packed with C18 AQ particles (3 μm, 120 Å) was used as the separation column, and 0.1% formic acid in water and in acetonitrile was used as mobile phases A and B, respectively (Song et al., 2012).

For the quantitative proteomic analysis, 20 μg tryptic digests were loaded onto the strong cation exchange (SCX) monolithic column and then were eluted by 11 salt stepwise elutions to perform the 2-D RPLC-MS/MS analysis (Wang et al., 2007). Each RPLC-MS/MS elution step was performed in a 72 min gradient from 5 to 35% acetonitrile. For the two-dimensional RPLC-MS/MS phosphoproteomic analysis, a RP-strong cation exchange chromatography (SCX) biphasic column combined with C18 analytical column was used (Bian et al., 2012, 2016). Briefly, the phosphopeptides were manually loaded onto the RP-SCX column and then eluted from the RP segment to SCX segment by a 160 min RP gradient LC-MS/MS (0 mM). Then, a series of stepwise elutions with salt concentrations of 24, 40, 56, 72, 100, 200, and 1,000 mM NH₄AC were used to elute phosphopeptides from the SCX column to the second dimensional C18 separation column. Finally, the RPLC-MS/MS separation was performed with a 117 min gradient from 5 to 25% acetonitrile.

The LTQ-Orbitrap Velos mass spectrometer (Thermo Fisher Scientific) was operated in data-dependent MS/MS acquisition mode. Full mass spectra were acquired in the Orbitrap at a resolution of 60,000 (m/z 400). The 20 most intense precursors were selected for fragmentation via collision induced dissociation in the LTQ. For the phosphopeptide analysis, the multistage activation was enabled. The dynamic exclusion function was set as follows: repeat count 2, repeat duration 30 s, and an exclusion duration of 60 s.

Protein identification and relative quantification

Protein identification and quantification are based on the MaxQuant software (version 1.1.1.36) according to the standard workflow (Cox and Mann, 2008). Because the m/z of the light and heavy labeled peptides were different at MS1 level, this mass difference can be distinguished by the mass spectrometer with resolution (R > ~10,000). Quantification is achieved by comparing the intensities of the labeled peptides at MS1 level (Boersema et al., 2009). At least 1.5 fold change is consider as differential proteins in response to salinity.

The database search was performed on the Andromeda search engine against the Chlorophyta protein database (from NCBI, containing 156,988 entries) concatenated with reversed sequences for evaluating of FDR (Cox et al., 2011). Carbamidomethylation on cysteine was set as a fixed modification, whereas oxidation on methionine was set as a variable modification. For phosphopeptides, Phospho (STY) was also set as variable modifications. FDR were set to 1% at phosphorylation site, peptide, and protein group level. For protein quantitation, only the proteins with a BenjaminiHochberg corrected p-value < 0.05 based on significance B were reported by the MaxQuant software (Cox and Mann, 2008). A more detailed parameters setting for database search was illustrated in our previous studies (Song et al., 2012). The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE (Vizcaino et al., 2014) partner repository with the dataset identifier PXD005443 and PXD005501.

Phosphorylation site localization

Only the phosphorylation sites with a reported localization probability of at least 0.75 were used for further analysis (Sharma et al., 2014). For quantification analysis, we used the intensity-weighted ratio significance values reported by MaxQuant to determine significantly changed phosphorylation sites.

Protein classification, subcellular location and protein–protein interaction prediction

Protein functional domains were predicted using the position-specific iterated BLAST and pattern-hit initiated BLAST programs (http://www.ncbi.nlm.nih.gov/BLAST/). Combined BLAST alignments with Gene Ontology and knowledge from literatures, proteins were classified into different categories.

The subcellular locations of proteins were determined using five internet tools: Yloc (http://abi.inf.uni-tuebingen.de/Services/YLoc/webloc.cgi), LocTree3 (https://rostlab.org/services/loctree3/), Plant-mPLoc (http://www.csbio.sjtu.edu.cn/bioinf/plant-multi/.), ngLOC (http://genome.unmc.edu/ngLOC/index.html), and TargetP (http://www.cbs.dtu.dk/services/TargetP/). The confident results have the consistent predictions from at least two tools. For the inconsistent prediction results among five tools, subcellular localizations for corresponding proteins were predicted based on literatures.

The protein–protein interactions (PPI) were predicted using the web-tool STRING 10 (http://string-db.org). The proteins homologs in Arabidopsis were analyzed by sequence BLASTing in TAIR database (http://www.arabidopsis.org/Blast/index.jsp), and then the homologs were subjected to the molecular interaction tool of STRING 10 for creating the proteome-scale interaction network (Suo et al., 2015).

Protein 3D structure analysis and multiple sequence alignment

The protein three-dimensional (3D) structure homology model was pre-calculated in SWISS-MODEL workspace (https://swissmodel.expasy.org) (Arnold et al., 2006). The 3D structures and the phosphorylation site were displayed using the SPDBV (version 4.1) software (Arnold et al., 2006). Multiple sequence alignment was performed with sequences from multiple species using Clustal Omega software (http://www.ebi.ac.uk/Tools/msa/clustalo/).

Statistical analysis

All physiological results were presented as means ± standard deviation of at least three replicates. Salinity-responsive proteins (SRPs) were quantified in at least two of the three replicates. The data were analyzed by Student's t-test using the statistical software SPSS 18.0 (SPSS Inc., Chicago, IL). A p-value less than 0.05 was considered statistically significant.

Results

Cell growth and viability upon salinity-induced palmella formation

The D. salina cells were cultured under medium containing 1 M NaCl (Katz and Avron, 1985). The cell number was increased gradually, and cell growth entered stationary phase after 10 days of culturing (Figure 1A), and then the cells were directly transferred to medium containing 3 M NaCl (Figure 1B and Supplemental Figure S1). The number of palmella cells was about 0.63 ± 0.16 × 10⁶ per milliliter (about 5% of the total cell number) at 2 h after salinity shock (HAS), and reached to 2.08 ± 0.25 × 10⁶ per milliliter (17%) at 6 HAS. The maximum of palmella cell number was detected as 6.09 ± 0.38 × 10⁶ per milliliter (49%) at 1 days after salinity shock (DAS), indicating the D. salina cells formed a palmella stage under the salinity induction at 1 DAS. The palmella stage continued for 3 days and reverted to free-swimming cell stage at 4 DAS (Figure 1B). The change of free-swimming cell number was also reflected from the medium color during this course (Supplemental Figure S1). After salinity shock, the medium color turned from green at 0 HAS to light green gradually at 1 DAS, and then reverted to dark green at 6 DAS (Supplemental Figure S1).

Figure 1

Besides, the change of cell viability evaluated by MTT assay was consistent with the number of free-swimming cells (Figure 1C). The cellular viability was gradually reduced at 6 HAS, reached to the lowest value at 1 DAS, and gradually increased from 4 DAS. In addition, the cell morphology has obvious changes during salinity shock-induced palmella formation. The free-swimming cells under 1 M NaCl were green, biflagellate, and rod to ovoid shaped, as well the cell length was about 12 μm (Figure 1D). The cell morphology didn't show obvious changes at 0 HAS (Figure 1E), but some cells started to lose their flagella and changed to round and immotile at 2 HAS (Figure 1F). More than 17% of free-swimming cells formed immotile palmella at 6 HAS, and the accumulated palmella cells lost their flagella and eyespot, appearing round (Figure 1G). All these indicate that 3 M NaCl shock induces the formation of a short palmella stage, but the palmella cells reverted to free-swimming cells at 4 DAS. During this course, there were three morphotypes at 6 HAS, which were predominantly motile cells, motile form with a prominent palmelloid cells (non-motile, mucilage rich), and palmelloid form with a weakly motile cells (Buchheim et al., 2010). Thus, 6 HAS is a critical time point for the transformation from free-swimming cells to palmella cells.

Osmotic regulation upon salinity-induced palmella formation

The glycerol content was significantly increased, and the activity of GPDH for glycerol synthesis was also obviously increased at 6 HAS (Figure 2A). Besides, the content of EPSs was significantly induced, but the cellular carbohydrate content was slightly decreased in cells at 6 HAS (Figure 2B). In addition, the activity of a carbohydrate metabolic enzyme PGM, which catalyzes the bidirectional interconversion of glucose-1-phosphate and glucose-6-phosphate, was significantly reduced at 6 HAS (Figure 2C). Importantly, the activity of plasma membrane (PM) H⁺-ATPase was significantly induced at 6 HAS, indicating a proton-motive force generator for sucrose transporter was enhanced (Figure 2D).

Figure 2

Chlorophyll content and fluorescence parameters

The contents of chlorophylls in control cells and salinity-shocked cells at 6 HAS were detected (Figure 3). The contents of chlorophyll a (Figure 3A), chlorophyll b (Figure 3A), and total chlorophyll (Figure 3B) were all declined, but the chlorophyll a/b (Figure 3B) didn't show obvious changes at 6 HAS.

Figure 3

Photosynthetic oxygen evolution and respiration rates were all decreased in the salinity- shocked cells when compared with those in the control cells (Figure 3C). In addition, chlorophyll fluorescence parameters were monitored to evaluate the photosynthetic performance. The PSII related parameters, including PSII maximum quantum yield (Fv/Fm) (Figure 3D), the excitation capture efficiency of open centers (Fv'/Fm') (Figure 3D), the effective PSII quantum yield (Y(II)) (Figure 3E), non-photochemical quenching (NPQ) (Figure 3E), and electron transport rate (ETR(II)) (Figure 3F), were declined in salinity-shocked cells. Furthermore, PSI related parameters, such as maximum fluorescence yield (Pm') (Figure 3G), non-photochemical quantum yields of PSI caused by donor-side limitation (Y(ND)) (Figure 3H) and by acceptor-side limitation (Y(NA)) (Figure 3H), were declined, whereas photochemical quantum yield of PSI (Y(I)) (Figure 3G) and electron transport rate (ETR(I)) were increased in salinity-shocked cells (Figure 3F).

Antioxidant enzyme activities and metabolite contents upon salinity shock-induced palmella formation

To evaluate the ROS level and the dynamics of ROS scavenging system in cells during palmella formation, the generation rate, H₂O₂ content, the activities of nine antioxidant enzymes, as well as the contents of metabolites in ROS scavenging system were analyzed (Figure 4). The generation rate was declined whereas the H₂O₂ content and GPX activity were increased at 6 HAS (Figures 4A,E). The activity of SOD was initially increased to catalyze the dismutation of into oxygen and H₂O₂ at 6 HAS (Figure 4B). However, the activities of POD and CAT in charge of the conversion of H₂O₂ to H₂O were decreased at 6 HAS (Figures 4B,C). Moreover, the activities of enzymes (i.e., APX, MDHAR, DHAR, and GR) as well as the contents of AsA and DHA in GSH-AsA cycle were also decreased at 6 HAS (Figures 4C–E,G), so did the contents of GSSG and GSH (Figure 4H) and the activity of GST (Figure 4F).

Figure 4

Proteomic analysis upon salinity shock-induced palmella formation

By using differential stable isotope labeling coupled with mass spectrometry approaches, 809 proteins were identified in D. salina cells (Figures 5A,B, Supplemental Table S1, and Supplemental Figure S2). This composes the largest protein database of D. salina so far (Supplemental Table S1). Among them, 509 proteins were quantified in three biological replicates (Figures 5A,B). On the basis of the Gene Ontology, BLAST alignment, KEGG database and information from the literature, the 809 proteins were classified into 14 functional categories (Figure 5C and Supplemental Table S1). Among them, other metabolism (20%), protein synthesis (13%), photosynthesis (13%), carbohydrate and energy metabolism (11%) were over-represented (Figure 5C).

Figure 5

Among these proteins, 151 SRPs were quantified in at least two of the three replicates, including 63 salinity-increased proteins, and 88 salinity-decreased proteins (Figure 5D, Tables 1, 2, and Supplemental Table S2). They were classified into 13 functional categories: membrane and transport, cell structure, carbohydrate and energy metabolism, other metabolisms, photosynthesis, transcription, protein synthesis, protein folding and transporting, protein degradation, ROS scavenging, signaling, cell division, differentiation and fate, and function unknown (Figure 5E and Supplemental Table S2). Among them, photosynthesis (24%), other metabolism (19%) and protein synthesis (13%) were predominant (Figure 5E). Most proteins involved in other metabolism, ROS scavenging, transcription, protein synthesis, and signaling were salinity-decreased (Figure 5E and Tables 1, 2). Besides, we found three SRPs involved in membrane and transport, and six involved in cell structure were changed upon palmella formation (Tables 1, 2).

In addition, mitochondria ATP synthase, NADH:ubiquinone oxidoreductase, and glycine cleavage system T protein were accumulated in cells at 6 HAS (Table 1). However, nine chloroplast- and/or cytoplasm- located proteins involved in pentose phosphate pathway, glyco-metabolism, pyruvate metabolism, and starch/sucrose metabolism were decreased at 6 HAS (Table 2). Among 36 photosynthetic SRPs, those in charge of light harvesting, PS I stability, electron transport, and ATP synthesis were increased, but the majority of proteins in carbon fixation and chlorophyll synthesis were decreased (Tables 1, 2), which were in line with the changes of photosynthesis and chlorophyll content upon palmella formation. Furthermore, 28 other metabolism-related SRPs were mainly involved in fatty acid biosynthesis, amino acid metabolism, terpenoid backbone biosynthesis, vitamins metabolism, purine and pyrimidine metabolism, and sulfur metabolism (Tables 1, 2).

Additionally, the protein profiles revealed the transcription, translation, as well as protein processing and fate were altered upon palmella formation. Four decreased histones and four RNA processing-related proteins indicate the transcription was inhibited upon palmella formation (Table 2). Besides, 20 SRPs involved in chloroplastic or cytoplastic protein synthesis (Tables 1, 2). In addition, the changes of 16 protein processing-related proteins imply that the protein folding, processing and transport were altered by salinity shock (Tables 1, 2). Furthermore, five protein degradation-related enzymes were increased, which were localized in cytoplasm, chloroplast, or mitochondria, respectively (Tables 1, 2). Besides, five signaling proteins and one cell division-related protein were decreased upon palmella formation (Table 2).

Phosphoproteomic changes during palmella formation

To investigate the changes of phosphoproteins upon palmella formation, stable isotope dimethyl labeling-based MS approaches were applied to analyze phosphoproteins, 137 phosphorylation sites were identified (Figure 6A), representing 100 phosphoproteins in 12 functional categories (Figure 6C). Among them, 54 phosphorylation sites were quantified (Figure 6B), representing 40 unique phosphoproteins (Figure 6D and Supplemental Tables S3, S4). The amino acid sequences of the 137 identified and 54 quantified phosphorylation sites were shown in Supplemental Table S4. The number of phosphopeptides or phosphorylation sites significantly changed at least a 1.5-fold in three replicates were shown in Supplemental Table S5. Among quantified phosphoproteins, photosynthesis (13%), signaling (13%) and cell structure (13%) were over-represented (Figure 6D). Finally, 35 salinity-responsive phosphoproteins (SRPPs) including 14 salinity-increased phosphoproteins and 22 salinity-decreased phosphoproteins were examined in D. salina during palmella formation (Figure 6E, Table 3, Supplemental Table S6, and Supplemental Figure S3). They were involved in photosynthesis, carbohydrate and energy metabolism, other metabolisms, transcription, protein synthesis, protein folding, protein degradation, signaling, and cell structure (Table 3 and Supplemental Table S6).

Figure 6

Phosphoprotein 3D structure homology model analysis

In this study, all the 35 SRPPs were subjected to SWISS-MODEL database (https://swissmodel.expasy.org) for building the 3D structure homology model. Among them, nine homologs of SRPPs were aligned to their 3D structure homology models (Figure 7 and Supplemental Table S7).

Figure 7

Importantly, the 3D structure clearly indicated that the phosphorylation sites of eight proteins were localized in the region of function domain (Figure 7), presenting useful information for underlying their specific functions. For pyruvate phosphate dikinase (PPDK), T506 was localized in beta turn, and was part of the phosphoenolpyruvate (PEP)-utilizers enzyme mobile domain (IPR008279), which was a “swiveling” beta/beta/alpha domain for catalyzing the transfer of a phosphoryl group from PEP to a histidine residue (Figure 7A). In RNA polymerase subunit 8 (RPB8), T72 was localized in OB-fold nucleic acid binding domain (IPR012340), which was a part of the five-stranded beta-barrel structure for nucleic acid recognition (Figure 7B). Besides, salinity-decreased T147 was located in DEAD-box helicases (DBH), which was part of DEAD/DEAH box helicase domain involved in various aspects of RNA metabolism (Figure 7C). For ribosomal proteins, the salinity-induced S69 of ribosomal protein S3a (RPS3a) was localized in beta sheet (Figure 7D), while T127 of ribosomal protein S6 (RPS6) was located in beta turn, which were all the conserved site of ribosomal protein S3Ae/S6e family (IPR018281/IPR001377) (Figure 7E). The salinity-induced T304 of membrane AAA-metalloprotease (FtsH) was located in beta turn of AAA⁺ domain (IPR003959), which has conserved alpha-beta-alpha core structure and walker A and B motifs of P-loop NTPase (Figure 7G). The salinity-induced Y251 of glycogen synthase kinase 3 (GSK3) was localized in protein kinase domain (IPR003959), which catalyzes the transfer of the gamma phosphate from ATP to one or more amino acid residues in a protein substrate, resulting in a conformational change affecting protein function (Figure 7H). In addition, salinity-decreased T349 phosphorylation of α-tubulin localized in the carboxy-terminal region which was the alpha domain interface region for motor protein binding (Figure 7I and Supplemental Figure S4). However, the salinity-induced S6 phosphorylation of ubiquitin-conjugating enzyme E2 (UBC) localized outside of its conserved function domain (Figure 7F). Whether the increase of S6 phosphorylation of UBC has an effect on conjugation of ubiquitin to the target protein needs to be confirmed.

Prediction of salinity-responsive protein–protein interaction upon palmella formation

To discover the relationship of the 151 SRPs, the PPI networks were generated using the web-tool STRING10 (http://www.string-db.org). After BLASTing in TAIR database (http://www.arabidopsis.org/Blast/index.jsp), 135 homologs in Arabidopsis of the 151 SRPs were analyzed, and then subjected to the molecular interaction tool of STRING 10 for creation of proteome-scale interaction network (Supplemental Table S8). Among them, 118 proteins were depicted in the STRING database, and illuminated in nine functional modules with tightly-connected clusters (stronger associations represented by thicker lines) in the network (Figure 8A). Besides, 23 SRPPs were depicted in the STRING database (Figure 8B and Supplemental Table S9). The relationship of proteins in various modules indicates that signaling and protein synthesis/processing/turnover are crucial for the modulation of light harvesting, carbon assimilation and energy apply, as well as the cytoskeleton upon palmella formation. Importantly, the reverse phosphorylation of proteins plays important roles in regulation of proteins interactions.

Figure 8

Discussion

Salinity shock-induces cell morphology changes upon palmella formation

D. salina has the ability to grow over an extremely wide range of salinity from less than sea water to NaCl saturation (Montoya and Olivera, 1993). The normal cells of D. salina are ellipsoid to fusiform, without cell wall, but covered with mucilaginous glycocalyx. The biflagellate motile cells contain a single large cup-shaped posterior chloroplast (Borowitzka and Siva, 2007). The osmotic shock, such as rapid decrease/increase of salinity, would induce palmella formation of D. salina (Preetha et al., 2012). The cells in palmella stage lose their flagella and eyespot, become round, and excrete a slime layer to format a multicellular aggregated colonies within a common mucilage (Montoya and Olivera, 1993).

Upon palmella formation, a salt-induced Bin/Amphiphysin/Rvs domain-containing protein would facilitate to bind membrane for membrane curvature modulation in D. salina (Figure 9A and Table 1; Rao and Haucke, 2011). Osmotic and salt stresses cause microtubule disassembly and reorganization (Komis et al., 2002; Shoji et al., 2006). In this study, the salt-induced β-tubulin implies that the tubulin accumulation in cytoplasmic pool of palmella cells may prepare for the flagella reformation of free-swimming cells (Figure 9A), while the salinity-decreased phosphorylation of T349 would facilitate the binding of α-tubulin with motor proteins to regulate the microtubule stability upon palmella formation (Figures 7I, 9A, and Supplemental Figure S4; Ban et al., 2013; Wang et al., 2013).

Figure 9

We found the flagella formation was modulated in response to salinity (Figure 9A). Flagellar associated proteins function as a molecular chaperon, and their declined abundance and phosphorylation level would contribute to the flagella lost during the palmella formation. Consistently, the abundance and phosphorylation at S581 of kinesin-like calmodulin binding protein (KCBP) were also decreased, but the normalized phosphorylation was unchanged upon palmella formation (Figure 9A). KCBP was tightly connected with tubulin in PPI network (Figure 8). The highly conserved KCBP is mainly localized near the basal body in D. salina (Shi et al., 2013) and C. reinhardtii (Dymek et al., 2006), which has a myosin tail homology-4 region in the N-terminal tail and a calmodulin-binding region following the motor domain (Dymek et al., 2006). KCBP plays an important role in microtubule-based intracellular motility. Therefore, the decrease of the abundance and phosphorylation of KCBP implies its motor activity may be inhibited upon palmella formation (Figure 9A and Table 3), but we cannot confirm its activity was regulated by phosphorylation. In addition, the increased abundance and phosphorylation level of striated fiber-assemblin would facilitate to regulate the flagellar root microtubule stability upon palmella formation (Figure 9A; Lechtreck et al., 2002).

Accumulation of exopolysaccharides and ROS homeostasis are crucial upon palmella formation

The Dunaliella cells were enclosed in a thin elastic plasma membrane surrounded by mucus “surface coat”, but lack rigid polysaccharide cell wall (Ben-Amotz and Avron, 1990). It is known that accumulation of glycerol in cells was crucial for osmotic homeostasis in D. salina in response to salinity (Pick, 2002; Liska et al., 2004). In this study, salinity shock-induced accumulation of glycerol and increase of the activity of glycerol metabolic enzyme are important for palmella formation (Figure 2A).

Similarly, total EPS contents were all increased in halotolerant cyanobacterium Microcoleus vaginatus (Chen et al., 2006), halotolerant bacterium Rhodopseudomonas acidophila (Sheng et al., 2006), and medicinal mushroom Phellinus linteus (Zou et al., 2006). In this study, EPS slime layer surrounded D. salina cells was enhanced (Figure 1G), due to the content of EPSs was increased upon palmella formation (Figure 2B). The salinity-induced EPSs would facilitate the accumulation of water and the reduction of ion influx, protecting the membrane system upon palmella formation. The EPSs around D. salina cells were excreted from cells through plasma membrane sugar transporters, which were energy-dependent H⁺-symporters (Büttner and Sauer, 2000). Its energization was via the proton-motive force generated by the PM H⁺-ATPase, and the modulation of H⁺-ATPase activity would immediately affect the sugar transport kinetics (Doidy et al., 2012). In this study, the abundance and activity of PM H⁺-ATPase were all increased to activate the efflux of sugars upon palmella formation (Figures 2D, 9B, and Table 1; Carpaneto et al., 2010). We also found ten carbohydrate metabolism-related proteins were all salinity-decreased upon palmella formation, which were involved in glycolysis, pentose phosphate pathway, as well as starch mobilization and glucose metabolism, respectively (Figure 9B and Tables 2, 3). The normalized phosphorylation level and activity of PGM were all decreased upon salinity-induced palmella formation (Table 3 and Supplemental Table S1), which implies that the phosphorylation-dependent activity of PGM is inhibited in D. salina under salinity-induced oxidative stress. Moreover, the phosphorylation level of glucose-6-phosphate isomerase and phosphomannose isomerase were also reduced upon palmella formation. It has been reported that glycolytic enzymes were oxidized and inactivated when cells were subjected to oxidative stress, leading to a metabolic reshuffling of glucose equivalents through the pentose phosphate pathway for providing a necessary reducing power (NADPH) for antioxidant defense mechanism in cells (Shanmuganathan et al., 2004). These indicate that carbohydrate metabolism may be reduced in cells upon salinity-induced palmella formation, similar to the salinity-decreased total carbohydrates in halotolerant cyanobacterium M. vaginatus (Chen et al., 2006).

In addition, ROS homeostasis is crucial for rapid metabolism transition upon salinity-shock palmella formation. In this study, the salinity shock-decreased activities of most ROS scavenging enzymes (i.e., POD, APX, MDHAR, DHAR, GR, and GST), the abundances of 2-Cys peroxiredoxin (PrxR) and APX (Figures 4B–F, 9J and Table 2), as well as the contents of GSSG and GSH (Figure 4H) imply that GSH-AsA cycle, PrxR/thioredoxin pathway, and APX pathway might be inhibited under salinity shock. However, the salinity-increased activities of SOD and GPX indicate that the dismutation of superoxide into oxygen and H₂O₂ and subsequently reduction in GPX pathway are enhanced for ROS scavenging upon palmella formation (Figures 4B,E, 9J). Importantly, the abundance of SOD was oxidative stress-increased during the early transition of green vegetative cells to red cysts in H. pluvialis (Wang et al., 2004a). Thus, the enhanced SOD and GPX pathways may be the key strategy for ROS homeostasis upon salinity-induced palmella formation.

Photosynthesis modulation upon palmella formation

Upon salinity shock-induced palmella formation, the photosynthesis of D. salina was significantly inhibited, being reflected from the obvious decline of chlorophyll content and photosynthesis oxygen evolution. A similar phenomenon (Figures 1D–G) also occurred upon copper-induced palmella formation of C. reinhardtii (Sztrum et al., 2012), and the transformation process of green vegetative cells into red aplanospores of H. pluvialis (Scibilia et al., 2015).

The excess light energy absorbance was reduced. The decline of the abundances of chlorophyll biosynthesis-related enzymes (i.e., Mg-protoporyphyrin IX chelatase, and gln-glu-trna synthetase) indicated that the chlorophyll biosynthesis may be inhibited upon palmella formation (Figure 9C and Table 2), which was also observed in H. pluvialis in response to nitrogen starvation (Scibilia et al., 2015). Besides, the decrease of chloroplast oxygen-evolving protein and PSII 47 kDa protein would attribute to the salt-inhibited photosynthesis oxygen evolution activity, which was supposed to facilitate the reduction of excess light energy absorbance in palmella cells (Figure 9C and Table 2).

The stability and activity of PSII were modulated by salinity shock. The PSII photochemistry activity [e.g., Fv/Fm, Fv'/Fm', and Y(II)] and linear electron flux were all obviously reduced upon palmella formation, which was similar with what happened in high light-stressed D. salina (Gu et al., 2014). However, the abundances of two PSII assembly-related proteins [PSII protein H (PsbH) and PSII protein VI] were increased upon palmella formation (Figure 9C and Table 1). PsbH is a 10-kDa phosphoprotein associated with the inner antenna PSII 47 kDa protein, which is conserved and essential for assembling PSII in algae and higher plants (Summer et al., 1997). The psbH mutants appeared a PSII deficient phenotype and lack of a functional PSII complex (Summer et al., 1997). Interestingly, we also found T8 phosphorylation of PsbH was decreased upon palmella formation (Figure 9C and Table 3). However, a T2 knockdown mutant of PsbH in C. reinhardtii has similar phonotype of wild-type strains, indicating that T2 phosphorylation of PsbH probably doesn't affect its function (O'Connor et al., 1998). Besides, PSII protein VI is the beta subunit of cytochrome b559 essential for PSII assembly, whose increase would enhance the cyclic electron transport (CET) or in a side path of electron flow for protecting PSII from photoinhibition (Figure 9C and Table 1; Burda et al., 2003).

Enhancement of CET would facilitate ATP synthesis. The salinity-suppressed ETR(II) and salinity-increased ETR(I) abundance implied that CET was tend to enhanced upon palmella formation in D. salina. The increased abundances of electron carrier protein cytochrome b₆/f complex and various PSI complex proteins (i.e., P700 chlorophyll a apoprotein A1, P700 chlorophyll a apoprotein A2, PsaG, PSI reaction center subunit II, PSI subunit VII, assembly protein Ycf3, and light-harvesting chlorophyll a/b protein 3) would contribute to the induced CET (Figure 9C and Table 1). This would facilitate funneling excess electrons to generate ATP without increasing oxygen evolution (Zhang et al., 2010; Gu et al., 2014). Interestingly, the normalized phosphorylation level of light harvesting chlorophyll a/b binding proteins (LHCII) of PSII (LHCb) was enhanced upon palmella formation (Figure 9C and Table 3). The reversible phosphorylation of LHCII is generally considered as an adaptation mechanism to balance energy distribution between PSII and PSI for regulating redox homeostasis in chloroplasts (Grieco et al., 2012). Our proteomic results revealed that the abundances of several enzymes in Calvin cycle were decreased upon salinity shock-induced palmella formation, such as ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO), phosphoglycerate kinase, and fructose-1, 6-bisphosphatase (Figure 9C and Table 2). The photosynthetic CO₂ fixation was inhibited due to the down-regulation of transcript levels of Calvin cycle enzymes in D. salina under high salt stress (Liska et al., 2004). Similarly, the decreased RuBisCO was also examined in H. pluvialis during the early transition of green vegetative cells to red cysts under oxidative stress (Wang et al., 2004a). Additionally, the abundances of chloroplastic ribulose phosphate-3-epimerase and transketolase also imply that the Calvin cycle may be inhibited in D. salina under salinity shock (Figure 9B and Table 2).

We also found the phosphorylation level of chloroplast-localized PPDK was decreased upon palmella formation (Figure 9C and Table 3), and PPDK was tightly connected with sedoheptulose-1,7-bisphosphatase in PPI network (Figure 8B). PPDK catalyzes the formation of PEP that is the initial acceptor of CO₂ in C4 pathway. In U. linza, the elevated PPDK transcription and enzyme activity enhanced C4 carbon metabolism under high salinity stress (Xu et al., 2013). In this study, the phosphorylation site T506 of PPDK is localized in beta turn of PEP-utilizers enzyme mobile domain (IPR008279) region (Figure 7A). It is reported that the reverse phosphorylation of PPDK (T456 in maize) would regulate the switch of its active (dephosphorylation) and inactive (phosphorylation) states (Chastain et al., 2002). Thus, the decreased phosphorylation at T506 of PPDK (homologous site with maize) would induce its activity, resulting in enhanced C4 pathway activity upon the salinity-induced palmella formation.

Several basic metabolisms are reduced in dormant palmella cells

The photosynthetic oxygen evolution, PSII activity, and photosynthetic CO₂ fixation were all reduced after free-swimming cells lost their flagella to be immotile (Figures 1, 3). Additionally, the abundances of several enzymes involved in fatty acid metabolism, amino acids metabolism, terpenoid backbone biosynthesis, and purine/pyrimidine metabolism were decreased (Figure 9D and Tables 2, 3), implying these pathways may be inhibited upon palmella formation. Similarly, proteins involved in hydrophobic fatty acid biosynthesis and amino acid synthesis were decreased in C. reinhardtii in response to salt stress for less than 24 h (Mastrobuoni et al., 2012) and nitrogen starvation for 6 h (Longworth et al., 2012). Besides, amino acid metabolism-related genes were down-regulated in brown alga Ectocarpus siliculosus under hypersaline conditions (Dittami et al., 2011). In addition, several nitrogen assimilation and pyruvate kinase metabolism-related proteins were decreased during the oxidative stress-induced early transition of green vegetative cells to red cysts in H. pluvialis (Wang et al., 2004a). All these indicate that the palmelloid cells become dormant under salinity shock condition. However, the energy production was induced upon palmella formation. This is consistent with what happened during oxidative stress-induced early transition of green vegetative cells to red cysts in H. pluvialis (Wang et al., 2004a), as well as in hypersaline-stressed brown alga E. siliculosus (Dittami et al., 2011), which is crucial for energy supply in dormant or stressed algae cells.

Nuclear and chloroplastic gene expression regulation upon palmella formation

The metabolic changes were regulated by gene expression pattern upon palmella formation. We found the abundances of ten transcription-related proteins and phosphorylation level of RPB8 were changed (Figure 9E and Tables 1–3). Three decreased RNA processing-related proteins (i.e., Sm protein B, nucleolar protein snoRNP, and exon junction complex) indicates the pre-mRNA splicing and localization are probably reduced (Figure 9F and Table 2; Tange et al., 2004). However, G-patch domain-containing protein, which functions in RNA recognition, RNA binding or splicing, was increased upon palmella formation (Figure 9F and Table 1). In addition, DBH is involved in various aspects of RNA metabolism (e.g., nuclear transcription, pre-mRNA splicing, and RNA decay) was salinity-altered, In our results, the salinity-reduced DBH is tightly connected with other protein synthesis-related proteins (Figure 8), and its phosphorylation site T147 located in helices (Figure 7C), implying the phosphorylation is perhaps involved in the regulation of DBH function (Figure 9F).

Besides, the phosphorylation at T72 of RPB8 was salinity-induced upon palmella formation, and T72 is localized in OB-fold nucleic acid binding domain (IPR012340) for nucleic acid recognition (Figures 7B, 9E and Table 3). Similarly, the S2 phosphorylation of RNA polymerases (RNA pol) II was increased in wild-type fission yeast under nitrogen starvation (Sukegawa et al., 2011), and RNA pol α subunit were also stress-increased in other algae, such as copper-stressed marine alga Scytosiphon gracilis (Contreras et al., 2010), salt-treated Bifidobacterium longum NCIMB 8809 (Sánchez et al., 2005), and acid (low pH)-stressed Streptococcus mutans (Len et al., 2004). Although the function of RPB8 phosphorylation remains unknown, its increase would affect RNA pol III assembly and catalytic function upon palmella formation of D. salina (Voutsina et al., 1999). Additionally, RPB8 was tightly connected with PRP1 splicing factor (PRP1) in PPI network (Figure 8B). Interestingly, the phosphorylation at T255 of PRP1 was decreased upon palmella formation (Figure 9F and Table 3). The N-terminus highly-conserved site of PRP1 can be phosphorylated by PRP4 during the spliceosome activation (Lützelberger et al., 2010). The decline of PRP1 phosphorylation implies that spliceosome activity is probably salinity reduced upon palmella formation.

De novo protein synthesis plays an important role in abiotic stress adaptation in plants. We found that eukaryotic initiation factor (eIF), eIF 4A-like protein, elongation factor (EF) 1 alpha 2 and EF2 were all reduced (Figure 9G and Table 2), indicating that translation initiation and peptide elongation are decreased upon palmella formation. Importantly, protein phosphorylation is involved in the initiation of protein translation (Jackson et al., 2010). We found that the phosphorylation levels of RPS6 (T127), RPS3a (S69), and ribosomal protein L12 (RPL12) (S27), and other ribosomal proteins (RPS3 and RPL23) were increased, but the abundances of ribosomal proteins S14 and ribosomal proteins L7 were decreased during palmella formation (Figure 9G and Tables 1–3). Among them, RPS6 phosphorylation in plants leads to the selective recruitment of ribosomal mRNAs to polysomes, regulating in the growth pattern of plants in response to environment changes (Reinbothe et al., 2010). In addition, the phosphorylation sites of RPS6 and RPS3a were all the conserved site of ribosomal protein family (Figures 7D,E). Therefore, the induced phosphorylation of T127 in RPS6 could be responsible for the elevated level of protein synthesis upon palmella formation. Besides, RPS3 is crucial in translation initiation via participation in rearrangements of the 40S subunit structure and promotion of subsequent recognition of the start codon and interaction with mRNA (Valášek, 2012). In our studies, the increased RPS3 abundance and the phosphorylation at S69 of RPS3a would contribute to translation initiation upon palmella formation. In addition, the RPL12 is a member of ribosomal L7/12 stalk, which functions in restoring the biological activity of “core” ribosomal particles (Ilag et al., 2005). Although the phosphorylation of RPL12 was reported, its exact function is still not clear.

Chloroplast and cytoplasm protein processing and turnover upon palmella formation

Under stress conditions, the misfolded proteins accumulated in algae cells should be refolded or removed. We found that 23 proteins and four phosphoproteins were involved in protein processing and degradation upon palmella formation (Figure 9H and Tables 1–3). The phosphorylation at S768 and S773 of prefoldin β subunit (PFDb) was decreased upon palmella formation (Figure 9H and Table 3). PFD is a cofactor of the group II chaperonins for capturing and transferring the unfolded actin and tubulin for correct folding (Millán-Zambrano and Chávez, 2014). Furthermore, the decrease of cytosolic chaperonin T-complex protein 1 β subunit (TCP) was examined upon palmella formation (Figure 9H and Table 2), which assists in the folding of newly synthesized actin, tubulin, and other polypeptides in vivo (Sternlicht et al., 1993). Previous proteomics results have showed that the cold-decreased PFD in cold-tolerant Zoysia japonica (Xuan et al., 2013) and salt-decreased TCP in Thellungiella halophila (Pang et al., 2010), indicating the protein folding is inhibited under stress. In our studies, the decrease of PFDb phosphorylation and the abundance of TCP might be related with microtubule organization upon palmella formation.

Besides, three endoplasmic reticulum-located chaperones, calreticulin (CALR), protein disulfide isomerase 1, and luminal binding protein Bip1, were all obviously increased upon palmella formation in D. salina (Figure 9H and Table 1). CALR is involved in the folding and quality control of newly synthesized proteins and glycoproteins, which is highly conserved and crucial for plant development and stress response (Garg et al., 2015). In addition, protein disulfide isomerase aids in the formation of proper disulfide bonds during protein folding in the endoplasmic reticulum (Appenzeller-Herzog and Ellgaard, 2008), which were salinity-increased in salt-tolerant Medicago sativa (Rahman et al., 2015), barley (Mostek et al., 2015), rice (Ghaffari et al., 2014), and halotolerant yeast (Rhodotorula mucilaginosa) (Lahav et al., 2004). Additionally, luminal binding protein Bip1 functions in precursor protein import and translocation (Wang et al., 2004b), which was high salt-induced in halotolerant yeast (Lahav et al., 2004). All these imply that salinity shock might enhance the protein folding and importation for maintaining the correct conformation and subcellular location of the oxidative proteins upon palmella formation.

In addition, two cytoplasmic heat shock protein 70s (HSP70) were decreased, but one chloroplast HSP70 was increased upon palmella formation (Figure 9H and Tables 1, 2). HSP70s assists the folding of de novo synthesized polypeptides and the translocation of precursor proteins in response to osmotic stress (Wang et al., 2004b), which was reported to increase in salinity-treated D. salina (Liska et al., 2004; Katz et al., 2007), and upon oxidative stress-induced early transition from green vegetative cells to red cysts in H. pluvialis (Wang et al., 2004a). This implies that the transportation and processing of certain newly synthesized peptide into chloroplasts are enhanced upon palmella formation of D. salina. Moreover, the chloroplast-located preprotein translocase SecY subunit is a component of the SecYEG translocon, being driven by the chloroplast ATPase SecA. Both SecY and SecA are essential for protein transportation from cytoplasm to chloroplasts (du Plessis et al., 2011). Their abundance changes may imply that the nuclear encoding proteins are necessary for chloroplasts in salinity shock-stressed D. salina (Figure 9H and Tables 1, 2). The increase of cyclophilin in D. salina would enhance the protein folding upon palmella formation (Figure 9H and Table 1; Zhang et al., 2012). The decrease of a mitochondrial GroEL and three chloroplastic chaperonin 60s would inhibit protein folding of RuBisCO, and the abundance of chaperonin 60s correlate positively with RuBisCO in plant cells (Figure 9H and Table 2; Holland et al., 1998).

Ubiquitin-dependent selective degradation of proteins is crucial upon palmella formation in D. salina (Ciechanover, 1998). In this study, salinity-induced abundance of 26S proteasome regulatory complex, as well as phosphorylation at S6 of UBC indicate that ubiquitin-dependent protein degradation in cytoplasm is enhanced upon palmella formation (Figure 9I and Tables 1, 3). Pervious proteomic study reported the increased 26S proteasome regulatory complex in D. salina under long-term salt stress (Liska et al., 2004). The increased phosphorylation at S6 of UBC located in the kinase motifs of CK2 implies that the conjugation of ubiquitin to the target protein may be enhanced (Supplemental Table S6).

The damaged proteins in chloroplasts need to be degraded through an efficient proteolytic system during chloroplast biogenesis, maintenance, and stress response (Ramundo et al., 2014). Upon palmella formation, the abundance and phosphorylation level of FtsH, as well as the abundance of ATP-dependent Clp protease proteolytic subunit were salinity-increased in chloroplasts in D. salina (Figure 9I and Tables 1, 3). FtsH can efficiently degrade proteins with low thermodynamic stability, and was increased in short-term salinity-stressed D. salina, osmotic shock-treated Oenococcus oeni (Bourdineaud et al., 2003), and heat/light-stressed cyanobacterium Synechocystis PCC 6803 (Silva et al., 2003). Additionally, FtsH was tightly connected with photosynthesis-related proteins in PPI network (Figure 8B). Therefore, the increase of T304 phosphorylation in AAA⁺ domain of FtsH would facilitate its assembly/disassembly with photosynthetic proteins to perform its peptidase function (Figures 7G, 9I and Table 3). Besides, the chloroplast gene encoding Clp protease is essential for cell viability (Ramundo et al., 2014), which was increased in D. salina under long-term salt stress (Katz et al., 2007). All these would enhance the removal of abnormal, modified, and mistargeted proteins in chloroplasts.

The mitochondria peptidase M16 was increased in palmella cells (Figure 9I and Table 1) and in D. salina cells under long-term salt stress (Liska et al., 2004). The peptidase M16 was also increased in suspension cells of A. thaliana (Ndimba et al., 2005), gametophore of Physcomitrella patens (Wang et al., 2008), and roots of O. sativa (Li et al., 2010) under salt stress. These results suggest the proteolytic reaction in mitochondrial processing may be important for stress-induced palmella formation.

Salinity-responsive signaling pathways are reduced upon palmella formation

The palmella formation is regulated by a sophisticate signaling network in D. salina. The salinity-decreased protein abundance and/or phosphorylation level imply that several crucial salinity-responsive signaling pathways (e.g., mitogen-activated protein kinase (MAPK) signaling, brassinosteroid (BR) signaling, and Snf1-like protein kinase (SnRK) signaling) tend to be inhibited upon salinity shock-induced palmella formation (Figure 9K and Tables 2, 3). Among them, MAPKs are highly conserved serine/threonine kinases in combination with their upstream activators, which convey osmotic stress signals to appropriate effectors and contribute to adaptation to the high salt stress (Moustafa et al., 2014). Additionally, SH2 domains are crucial for protein docking to phosphorylated tyrosine residues on other proteins, which are common in adapter proteins that aid in the signal transduction of receptor tyrosine kinase pathways (Koytiger et al., 2013). The decrease of MAPK and SH2 domain containing protein indicated that the active salt-responsive MAPK signaling might be inhibited in the dormant palmelloid cells. Moreover, in the MAPK signaling pathway, GSK3 can phosphorylate MAP1B to regulate flagellar length and assembly (Wilson and Lefebvre, 2004). The phosphorylation of Y218 in C. reinhardtii GSK3 regulated its active states (Wilson and Lefebvre, 2004). The GSK3 in D. salina shared the conserved phosphorylation site with that in C. reinhardtii. This conserved phosphorylation site was localized in protein kinase domain (Figure 7H). Thus, the decreased phosphorylation of Y251 in GSK3 of D. salina would inhibit its activity to reduce the flagellar stability upon palmella formation.

The salinity-decreased calmodulin-like protein and KCBP would contribute to inhibition of microtubule-binding activity of motor domain during palmella formation of D. salina (Shi et al., 2013). The decreased phosphorylation level at S173 and S174 of SnRK indicates that SnRK-related signaling pathway for regulation of ion homeostasis and ROS production may be reduced upon palmella formation (Gong et al., 2002; Diedhiou et al., 2008; Mao et al., 2010; Kulik et al., 2012). The decreased phosphrylation level at S497 of Bsu 1 phosphatase would reduce its activity, leading to the inhibition of nuclear transcription factors to repress BR-responsive gene expression upon palmella formation (Mora-García et al., 2004).

Conclusions

In the life cycle of unicellular algae, palmella stage is critical for cell surviving in various stress conditions. Salinity-induced palmella formation is a fine-tuned cellular process. By integrating analysis of physiological, quantitative proteomics, and phosphoproteomics data, we revealed the specific molecular mechanisms upon palmella formation (Figure 9). They mainly include (1) cell membrane curvature and cytoskeleton dynamics are modulated for cell morphological changes, (2) accumulations of glycerol and EPSs are enhanced for protection of membrane system, (3) SOD and GPX pathways are specific for ROS scavenging, (4) the activities of photosynthesis oxygen evolution and PSII are inhibited, but the CET and C4 pathway are enhanced, (5) nuclear and chloroplastic gene expression are regulated in response to salinity, and (6) chloroplast and cytoplasm protein processing and turnover are enhanced. All these provide novel insights into the underlying salinity shock-induced palmella formation.

References

• 1

  AlkayalF.AlbionR. L.TillettR. L.HathwaikL. T.LemosM. S.CushmanJ. C. (2010). Expressed sequence tag (EST) profiling in hyper saline shocked Dunaliella salina reveals high expression of protein synthetic apparatus components. Plant Sci.179, 437–449. 10.1016/j.plantsci.2010.07.001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 2

  Appenzeller-HerzogC.EllgaardL. (2008). The human PDI family: versatility packed into a single fold. Biochim. Biophys. Acta1783, 535–548. 10.1016/j.bbamcr.2007.11.010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3

  ArnoldK.BordoliL.KoppJ.SchwedeT. (2006). The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling. Bioinformatics22, 195–201. 10.1093/bioinformatics/bti770

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4

  BakerM. A.CernigliaG. J.ZamanA. (1990). Microtiter plate assay for the measurement of glutathione and glutathione disulfide in large numbers of biological samples. Anal. Biochem.190, 360–365. 10.1016/0003-2697(90)90208-Q

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5

  BanY.KobayashiY.HaraT.HamadaT.HashimotoT.TakedaS. (2013). Alpha-tubulin is rapidly phosphorylated in response to hyperosmotic stress in rice and Arabidopsis. Plant Cell Physiol.54, 848–858. 10.1093/pcp/pct065

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 6

  Ben-AmotzA.AvronM. (1973). The role of glycerol in the osmotic regulation of the halophilic alga Dunaliella parva. Plant Physiol.51, 875–878. 10.1104/pp.51.5.875

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7

  Ben-AmotzA.AvronM. (1990). The biotechnology of cultivating the halotolerant algae Dunaliella. Trends Biotechnol.8, 121–126. 10.1016/0167-7799(90)90152-N

  • CrossRef
  • Google Scholar

• 8

  BianY.LiL.DongM.LiuX.KanekoT.ChengK.et al. (2016). Ultra-deep tyrosine phosphoproteomics enabled by a phosphotyrosine superbinder. Nat. Chem. Biol.12, 959–966. 10.1038/nchembio.2178

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 9

  BianY.YeM.SongC.ChengK.WangC.WeiX.et al. (2012). Improve the coverage for the analysis of phosphoproteome of HeLa cells by a tandem digestion approach. J. Proteome Res.11, 2828–2837. 10.1021/pr300242w

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10

  BoersemaP. J.RaijmakersR.LemeerS.MohammedS.HeckA. J. (2009). Multiplex peptide stable isotope dimethyl labeling for quantitative proteomics. Nat. Protoc.4, 484–494. 10.1038/nprot.2009.21

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11

  BorowitzkaM. A.SivaC. J. (2007). The taxonomy of the genus Dunaliella (Chlorophyta, Dunaliellales) with emphasis on the marine and halophilic species. J. Appl. Phycol.19, 567–590. 10.1007/s10811-007-9171-x

  • CrossRef
  • Google Scholar

• 12

  BourdineaudJ. P.NehmeB.TesseS.Lonvaud-FunelA. (2003). The ftsH gene of the wine bacterium Oenococcus oeni is involved in protection against environmental stress. Appl. Environ. Microbiol.69, 2512–2520. 10.1128/AEM.69.5.2512-2520.2003

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13

  BoussibaS.BingW.YuanJ. P.ZarkaA.ChenF. (1999). Changes in pigments profile in the green alga Haeamtococcus pluvialis exposed to environmental stresses. Biotechnol. Lett.21, 601–604. 10.1023/A:1005507514694

  • CrossRef
  • Google Scholar

• 14

  BuchheimM. A.KirkwoodA. E.BuchheimJ. A.VergheseB.HenleyW. J. (2010). Hypersaline soil supports a diverse community of dunaliella (Chlorophyceae)1. J. Phycol.46, 1038–1047. 10.1111/j.1529-8817.2010.00886.x

  • CrossRef
  • Google Scholar

• 15

  BurdaK.KrukJ.BorgstädtR.StanekJ.StrzalkaK.SchmidG. H. (2003). Mössbauer studies of the non-heme iron and cytochrome b559 in a Chlamydomonas reinhardtii PSI- mutant and their interactions with alpha-tocopherol quinone. FEBS Lett.535, 159–165. 10.1016/S0014-5793(02)03895-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 16

  BüttnerM.SauerN. (2000). Monosaccharide transporters in plants: structure, function and physiology. Biochim. Biophys. Acta1465, 263–274. 10.1016/S0005-2736(00)00143-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17

  CarpanetoA.KoepsellH.BambergE.HedrichR.GeigerD. (2010). Sucrose-and H⁺-dependent charge movements associated with the gating of sucrose transporter ZmSUT1. PLoS ONE5:e12605. 10.1371/journal.pone.0012605

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18

  ChastainC. J.FriesJ. P.VogelJ. A.RandklevC. L.VossenA. P.DittmerS. K.et al. (2002). Pyruvate, orthophosphate dikinase in leaves and chloroplasts of C(3) plants undergoes light-/dark-induced reversible phosphorylation. Plant Physiol.128, 1368–1378. 10.1104/pp.010806

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 19

  ChenH.JiangJ. G.WuG. H. (2009). Effects of salinity changes on the growth of Dunaliella salina and its isozyme activities of glycerol-3-phosphate dehydrogenase. J. Agric. Food Chem.57, 6178–6182. 10.1021/jf900447r

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20

  ChenH.LaoY. M.JiangJ. G. (2011). Effects of salinities on the gene expression of a (NAD⁺)-dependent glycerol-3-phosphate dehydrogenase in Dunaliella salina. Sci. Total Environ.409, 1291–1297. 10.1016/j.scitotenv.2010.12.038

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 21

  ChenL. Z.LiD. H.SongL. R.HuC. X.WangG. H.LiuY. D. (2006). Effects of salt stress on carbohydrate metabolism in desert soil alga Microcoleus vaginatus Gom. J. Integr. Plant Biol.48, 914–919. 10.1111/j.1744-7909.2006.00291.x

  • CrossRef
  • Google Scholar

• 22

  CiechanoverA. (1998). The ubiquitin–proteasome pathway: on protein death and cell life. EMBO J.17, 7151–7160. 10.1093/emboj/17.24.7151

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 23

  ContrerasL.MoenneA.GaillardF.PotinP.CorreaJ. A. (2010). Proteomic analysis and identification of copper stress-regulated proteins in the marine alga Scytosiphon gracilis (Phaeophyceae). Aquat. Toxicol.96, 85–89. 10.1016/j.aquatox.2009.10.007

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 24

  CoxJ.MannM. (2008). MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat. Biotechnol.26, 1367–1372. 10.1038/nbt.1511

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25

  CoxJ.NeuhauserN.MichalskiA.ScheltemaR. A.OlsenJ. V.MannM. (2011). Andromeda: a peptide search engine integrated into the MaxQuant environment. J. Proteome Res.10, 1794–1805. 10.1021/pr101065j

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 26

  DiedhiouC. J.PopovaO. V.DietzK. J.GolldackD. (2008). The SNF1-type serine-threonine protein kinase SAPK4 regulates stress-responsive gene expression in rice. BMC Plant Biol.8:49. 10.1186/1471-2229-8-49

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 27

  DittamiS. M.GravotA.RenaultD.GoulitquerS.EggertA.BouchereauA.et al. (2011). Integrative analysis of metabolite and transcript abundance during the short-term response to saline and oxidative stress in the brown alga Ectocarpus siliculosus. Plant Cell Environ.34, 629–642. 10.1111/j.1365-3040.2010.02268.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28

  DoidyJ.GraceE.KühnC.Simon-PlasF.CasieriL.WipfD. (2012). Sugar transporters in plants and in their interactions with fungi. Trends Plant Sci.17, 413–422. 10.1016/j.tplants.2012.03.009

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29

  du PlessisD. J.NouwenN.DriessenA. J. (2011). The Sec translocase. Biochim. Biophys. Acta1808, 851–865. 10.1016/j.bbamem.2010.08.016

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30

  DymekE. E.GodutiD.KramerT.SmithE. F. (2006). A kinesin-like calmodulin-binding protein in Chlamydomonas: evidence for a role in cell division and flagellar functions. J. Cell Sci.119, 3107–3116. 10.1242/jcs.03028

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31

  EllisR. J. (1972). Control of palmelloid formation in the green alga Pediastrum. Plant Cell Physiol.13, 663–672. 10.1093/oxfordjournals.pcp.a074778

  • CrossRef
  • Google Scholar

• 32

  FranqueiraD.OrosaM.TorresE.HerreroC.CidA. (2000). Potential use of flow cytometry in toxicity studies with microalgae. Sci. Total Environ.247, 119–126. 10.1016/S0048-9697(99)00483-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 33

  GargG.YadavS.YadavG. (2015). Key roles of calreticulin and calnexin proteins in plant perception under stress conditions: A review. Adv. Life Sci.5, 18–26. 10.5923/j.als.20150501.03

  • CrossRef
  • Google Scholar

• 34

  GhaffariA.GharechahiJ.NakhodaB.SalekdehG. H. (2014). Physiology and proteome responses of two contrasting rice mutants and their wild type parent under salt stress conditions at the vegetative stage. J. Plant Physiol.171, 31–44. 10.1016/j.jplph.2013.07.014

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35

  GongD.GongZ.GuoY.ChenX.ZhuJ. K. (2002). Biochemical and functional characterization of PKS11, a novel Arabidopsis protein kinase. J. Biol. Chem.277, 28340–28350. 10.1074/jbc.M107719200

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36

  GriecoM.TikkanenM.PaakkarinenV.KangasjarviS.AroE. M. (2012). Steady-state phosphorylation of light-harvesting complex II proteins preserves photosystem I under fluctuating white light. Plant Physiol.160, 1896–1910. 10.1104/pp.112.206466

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37

  GuW.LiH.ZhaoP.YuR.PanG.GaoS.et al. (2014). Quantitative proteomic analysis of thylakoid from two microalgae (Haematococcus pluvialis and Dunaliella salina) reveals two different high light-responsive strategies. Sci. Rep.4:6661. 10.1038/srep06661

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38

  HaghjouM. M.ShariatiM.SmirnoffN. (2009). The effect of acute high light and low temperature stresses on the ascorbate–glutathione cycle and superoxide dismutase activity in two Dunaliella salina strains. Physiol. Plant.135, 272–280. 10.1111/j.1399-3054.2008.01193.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 39

  HanD.WangJ.SommerfeldM.HuQ. (2012). Susceptibility and protective mechanisms of motile and non motile cells of Haematococcus pluvialis (chlorophyceae) to photooxidative stress1. J. Phycol.48, 693–705. 10.1111/j.1529-8817.2012.01147.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 40

  HollandN.BelkindA.HollandD.PickU.EdelmanM. (1998). Stress-responsive accumulation of plastid chaperonin 60 during seedling development. Plant J.13, 311–316. 10.1016/j.bbrc.2012.04.137

  • CrossRef
  • Google Scholar

• 41

  IlagL. L.VidelerH.McKayA. R.SobottF.FuciniP.NierhausK. H.et al. (2005). Heptameric (L12) 6/L10 rather than canonical pentameric complexes are found by tandem MS of intact ribosomes from thermophilic bacteria. Proc. Natl. Acad. Sci. U.S.A.102, 8192–8197. 10.1073/pnas.0502193102

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 42

  IwasaK.MurakamiS. (1968). Palmelloid formation of Chlamydomonas. Physiol. Plant.21, 1224–1233. 10.1111/j.1399-3054.1968.tb07353.x

  • CrossRef
  • Google Scholar

• 43

  JacksonR. J.HellenC. U.PestovaT. V. (2010). The mechanism of eukaryotic translation initiation and principles of its regulation. Nat. Rev. Mol. Cell Biol.11, 113–127. 10.1038/nrm2838

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 44

  KampfenkelK.Van MontaguM.InzeD. (1995). Extraction and determination of ascorbate and dehydroascorbate from plant tissue. Anal. Biochem.225, 165–167. 10.1006/abio.1995.1127

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 45

  KasamoK. (1988). Inhibition of Tonoplast and Plasma Membrane H⁺-ATPase Activity in Rice (Oryza sativa L.) Culture Cells by Local Anesthetics. Plant Cell Physiol.29, 215–221. 10.1093/oxfordjournals.pcp.a077483

  • CrossRef
  • Google Scholar

• 46

  KatzA.AvronM. (1985). Determination of intracellular osmotic volume and sodium concentration in dunaliella. Plant Physiol.78, 817–820. 10.1104/pp.78.4.817

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 47

  KatzA.WaridelP.ShevchenkoA.PickU. (2007). Salt-induced changes in the plasma membrane proteome of the halotolerant alga Dunaliella salina as revealed by blue native gel electrophoresis and nano-LC-MS/MS analysis. Mol. Cell. Proteomics6, 1459–1472. 10.1074/mcp.m700002-mcp200

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 48

  KomisG.ApostolakosP.GalatisB. (2002). Hyperosmotic stress induces formation of tubulin macrotubules in root-tip cells of Triticum turgidum: their probable involvement in protoplast volume control. Plant Cell Physiol.43, 911–922. 10.1093/pcp/pcf114

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 49

  KoytigerG.KaushanskyA.GordusA.RushJ.SorgerP.MacBeathG. (2013). Phosphotyrosine Signaling Proteins that Drive Oncogenesis Tend to be Highly Interconnected. Mol. Cell. Proteomics12, 1204–1213. 10.1074/mcp.M112.025858

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 50

  KulikA.Anielska-MazurA.BucholcM.KoenE.SzymanskaK.ZmienkoA.et al. (2012). SNF1-related protein kinases type 2 are involved in plant responses to cadmium stress. Plant Physiol.160, 868–883. 10.1104/pp.112.194472

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 51

  LahavR.NejidatA.AbeliovichA. (2004). Alterations in protein synthesis and levels of heat shock 70 proteins in response to salt stress of the halotolerant yeast Rhodotorula mucilaginosa. Antonie Van Leeuwenhoek85, 259–269. 10.1023/b:anto.0000020361.81006.2b

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 52

  LechtreckK. F.RostmannJ.GrunowA. (2002). Analysis of Chlamydomonas SF-assemblin by GFP tagging and expression of antisense constructs. J. Cell Sci.115, 1511–1522. Available online at: https://apps.crossref.org/titleList/

  • Pubmed Abstract
  • Google Scholar

• 53

  LenA. C.HartyD. W.JacquesN. A. (2004). Stress-responsive proteins are upregulated in Streptococcus mutans during acid tolerance. Microbiology150, 1339–1351. 10.1099/mic.0.27008-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 54

  LiH. S.SunQ.ZhaoS. J.ZhangW. H. (2000). Principles and Techniques of Plant Physiological Biochemical Experiment.Beijing: Higher Education Press.

  • Google Scholar

• 55

  LiX. J.YangM. F.ChenH.QuL. Q.ChenF.ShenS. H. (2010). Abscisic acid pretreatment enhances salt tolerance of rice seedlings: proteomic evidence. Biochim. Biophys. Acta1804, 929–940. 10.1016/j.bbapap.2010.01.004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 56

  LiangY. (1999). Effects of silicon on enzyme activity and sodium, potassium and calcium concentration in barley under salt stress. Plant Soil209, 217–224. 10.1023/A:1004526604913

  • CrossRef
  • Google Scholar

• 57

  LichtenthalerH. K.WellburnA. R. (1983). Determinations of total carotenoids and chlorophylls a and b of leaf extracts in different solvents. Biochem. Soc. Trans.11, 591–592. 10.1042/bst0110591

  • CrossRef
  • Google Scholar

• 58

  LiskaA. J.ShevchenkoA.PickU.KatzA. (2004). Enhanced photosynthesis and redox energy production contribute to salinity tolerance in Dunaliella as revealed by homology-based proteomics. Plant Physiol.136, 2806–2817. 10.1104/pp.104.039438

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 59

  LongworthJ.NoirelJ.PandhalJ.WrightP. C.VaidyanathanS. (2012). HILIC-and SCX-based quantitative proteomics of Chlamydomonas reinhardtii during nitrogen starvation induced lipid and carbohydrate accumulation. J. Proteome Res.11, 5959–5971. 10.1021/pr300692t

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 60

  LurlingM.BeekmanW. (2006). Palmelloids formation in Chlamydomonas reinhardtii: defence against rotifer predators?. Ann. Limnol. Int. J. Lim.42, 65–72. 10.1051/limn/2006010

  • CrossRef
  • Google Scholar

• 61

  LützelbergerM.BottnerC. A.SchwelnusW.Zock-EmmenthalS.RazanauA.KauferN. F. (2010). The N-terminus of Prp1 (Prp6/U5-102 K) is essential for spliceosome activation in vivo. Nucleic Acids Res.38, 1610–1622. 10.1093/nar/gkp1155

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 62

  MaW.ShiD.WangQ.WeiL.ChenH. (2005). Exogenous expression of the wheat chloroplastic fructose-1, 6-bisphosphatase gene enhances photosynthesis in the transgenic cyanobacterium, Anabaena PCC7120. J. Appl. Phycol.17, 273–280. 10.1007/s10811-005-4850-y

  • CrossRef
  • Google Scholar

• 63

  MaW.WeiL.WangQ. (2008). The response of electron transport mediated by active NADPH dehydrogenase complexes to heat stress in the cyanobacterium Synechocystis 6803. Sci. China Life Sci.51, 1082–1087. 10.1007/s11427-008-0139-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 64

  MaoX.ZhangH.TianS.ChangX.JingR. (2010). TaSnRK2.4, an SNF1-type serine/threonine protein kinase of wheat (Triticum aestivum L.), confers enhanced multistress tolerance in Arabidopsis. J. Exp. Bot.61, 683–696. 10.1093/jxb/erp331

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 65

  MastrobuoniG.IrgangS.PietzkeM.AßmusH. E.WenzelM.SchulzeW. X.et al. (2012). Proteome dynamics and early salt stress response of the photosynthetic organism Chlamydomonas reinhardtii. BMC Genomics13:215. 10.1186/1471-2164-13-215

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 66

  MendesL. F.Zambotti-VillelaL.ColepicoloP.Marinho-SorianoE.StevaniC. V.YokoyaN. S. (2013). Metal cation toxicity in the alga Gracilaria domingensis as evaluated by the daily growth rates in synthetic seawater. J. Appl. Phycol.25, 1939–1947. 10.1007/s10811-013-0036-1

  • CrossRef
  • Google Scholar

• 67

  Millán-ZambranoG.ChávezS. (2014). Nuclear functions of prefoldin. Open Biol.4:140085. 10.1098/rsob.140085

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 68

  MontoyaH. T.OliveraA. G. (1993). Dunaliella salina from saline environments of the central coast of Peru. Hydrobiologia155–161. 10.1007/978-94-011-2076-0_11

  • CrossRef
  • Google Scholar

• 69

  Mora-GarcíaS.VertG.YinY.Caño-DelgadoA.CheongH.ChoryJ. (2004). Nuclear protein phosphatases with Kelch-repeat domains modulate the response to brassinosteroids in Arabidopsis. Genes Dev.18, 448–460. 10.1101/gad.1174204

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 70

  MostekA.BornerA.BadowiecA.WeidnerS. (2015). Alterations in root proteome of salt-sensitive and tolerant barley lines under salt stress conditions. J. Plant Physiol.174, 166–176. 10.1016/j.jplph.2014.08.020

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 71

  MoustafaK.AbuQamarS.JarrarM.Al-RajabA. J.Tremouillaux-GuillerJ. (2014). MAPK cascades and major abiotic stresses. Plant Cell Rep.33, 1217–1225. 10.1007/s00299-014-1629-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 72

  NdimbaB. K.ChivasaS.SimonW. J.SlabasA. R. (2005). Identification of Arabidopsis salt and osmotic stress responsive proteins using two-dimensional difference gel electrophoresis and mass spectrometry. Proteomics5, 4185–4196. 10.1002/pmic.200401282

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 73

  O'ConnorH. E.RuffleS. V.CainA. J.DeakZ.VassI.NugentJ. H.et al. (1998). The 9-kDa phosphoprotein of photosystem II. Generation and characterisation of Chlamydomonas mutants lacking PSII-H and a site-directed mutant lacking the phosphorylation site. Biochim. Biophys. Acta1364, 63–72. 10.1016/S0005-2728(98)00013-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 74

  PangQ.ChenS.DaiS.ChenY.WangY.YanX. (2010). Comparative proteomics of salt tolerance in Arabidopsis thaliana and Thellungiella halophila. J. Proteome Res.9, 2584–2599. 10.1021/pr100034f

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 75

  PickU. (2002). Adaptation of the halotolerant alga Dunaliella to high salinity, in Salinity: Environment-Plants-Molecules, eds LäuchliA.LüettgeU. (Dordrecht: Kluwer Academic Publishers), 97–112. 10.1007/0-306-48155-3_5

  • CrossRef
  • Google Scholar

• 76

  PistocchiR.GuerriniF.BalboniV.BoniL. (1997). Copper toxicity and carbohydrate production in the microalgae Cylindrotheca fusiformis and Gymnodinium sp. Eur. J. Phycol.32, 125–132. 10.1080/09670269710001737049

  • CrossRef
  • Google Scholar

• 77

  PreethaK.JohnL.SubinC. S.VijayanK. K. (2012). Phenotypic and genetic characterization of Dunaliella (Chlorophyta) from Indian salinas and their diversity. Aquat. Biosyst.8:1. 10.1186/2046-9063-8-27

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 78

  RahmanM. A.AlamI.KimY. G.AhnN. Y.HeoS. H.LeeD. G.et al. (2015). Screening for salt-responsive proteins in two contrasting alfalfa cultivars using a comparative proteome approach. Plant Physiol. Biochem.89, 112–122. 10.1016/j.plaphy.2015.02.015

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 79

  RamundoS.CaseroD.MühlhausT.HemmeD.SommerF.CrèvecoeurM.et al. (2014). Conditional depletion of the Chlamydomonas chloroplast ClpP protease activates nuclear genes involved in autophagy and plastid protein quality control. The Plant Cell26, 2201–2222. 10.1105/tpc.114.124842

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 80

  RaoY.HauckeV. (2011). Membrane shaping by the Bin/amphiphysin/Rvs (BAR) domain protein superfamily. Cell. Mol. Life Sci.68, 3983–3993. 10.1007/s00018-011-0768-5

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 81

  ReinbotheC.PollmannS.ReinbotheS. (2010). Singlet oxygen signaling links photosynthesis to translation and plant growth. Trends Plant Sci.15, 499–506. 10.1016/j.tplants.2010.05.011

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 82

  SánchezB.Champomier-VergèsM. C.AngladeP.BaraigeF.de Los Reyes-GavilánC. G.MargollesA.et al. (2005). Proteomic analysis of global changes in protein expression during bile salt exposure of Bifidobacterium longum NCIMB 8809. J. Bacteriol.187, 5799–5808. 10.1128/JB.187.16.5799-5808.2005

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 83

  ScibiliaL.GirolomoniL.BerteottiS.AlboresiA.BallottariM. (2015). Photosynthetic response to nitrogen starvation and high light in Haematococcus pluvialis. Algal Res.12, 170–181. 10.1016/j.algal.2015.08.024

  • CrossRef
  • Google Scholar

• 84

  ShanmuganathanA.AveryS. V.WillettsS. A.HoughtonJ. E. (2004). Copper-induced oxidative stress in Saccharomyces cerevisiae targets enzymes of the glycolytic pathway. FEBS Lett.556, 253–259. 10.1016/S0014-5793(03)01428-5

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 85

  SharmaK.D'SouzaR. C.TyanovaS.SchaabC.WiśniewskiJ. R.CoxJ.et al. (2014). Ultradeep human phosphoproteome reveals a distinct regulatory nature of Tyr and Ser/Tr-based signaling. Cell Rep.8, 1583–1594. 10.1016/j.celrep.2014.07.036

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 86

  ShengG.YuH.YueZ. (2006). Factors influencing the production of extracellular polymeric substances by Rhodopseudomonas acidophila. Int. Biodeterior. Biodegrad.58, 89–93. 10.1016/j.ibiod.2006.07.005

  • CrossRef
  • Google Scholar

• 87

  ShiK.CuiL.JiangH.YangL.XueL. (2013). Characterization of the microtubule-binding activity of kinesin-like calmodulin binding protein from Dunaliella salina. Res. Microbiol.164, 1028–1034. 10.1016/j.resmic.2013.08.009

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 88

  ShojiT.SuzukiK.AbeT.KanekoY.ShiH.ZhuJ. K.et al. (2006). Salt stress affects cortical microtubule organization and helical growth in Arabidopsis. Plant Cell Physiol.47, 1158–1168. 10.1093/pcp/pcj090

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 89

  SilvaP.ThompsonE.BaileyS.KruseO.MullineauxC. W.RobinsonC.et al. (2003). FtsH is involved in the early stages of repair of photosystem II in Synechocystis sp PCC 6803. The Plant Cell15, 2152–2164. 10.1105/tpc.012609

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 90

  SongC.WangF.ChengK.WeiX.BianY.WangK.et al. (2014). Large-scale quantification of single amino-acid variations by a variation-associated database search strategy. J. Proteome Res.13, 241–248. 10.1021/pr400544j

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 91

  SongC.YeM.LiuZ.ChengH.JiangX.HanG.et al. (2012). Systematic analysis of protein phosphorylation networks from phosphoproteomic data. Mol. Cell. Proteomics11, 1070–1083. 10.1074/mcp.M111.012625

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 92

  SternlichtH.FarrG. W.SternlichtM. L.DriscollJ. K.WillisonK.YaffeM. B. (1993). The t-complex polypeptide 1 complex is a chaperonin for tubulin and actin in vivo. Proc. Natl. Acad. Sci. U.S.A.90, 9422–9426.

  • Pubmed Abstract
  • Google Scholar

• 93

  SukegawaY.YamashitaA.YamamotoM. (2011). The fission yeast stress-responsive MAPK pathway promotes meiosis via the phosphorylation of Pol II CTD in response to environmental and feedback cues. PLoS Genet.7:e1002387. 10.1371/journal.pgen.1002387

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 94

  SummerE. J.SchmidV. H.BrunsB. U.SchmidtG. W. (1997). Requirement for the H phosphoprotein in photosystem II of Chlamydomonas reinhardtii. Plant Physiol.113, 1359–1368.

  • Pubmed Abstract
  • Google Scholar

• 95

  SuoJ.ZhaoQ.ZhangZ.ChenS.CaoJ.LiuG.et al. (2015). Cytological and proteomic analyses of Osmunda cinnamomea germinating spores reveal characteristics of fern spore germination and rhizoid tip-growth. Mol. Cell. Proteomics14, 2510–2534. 10.1074/mcp.M114.047225

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 96

  SztrumA. A.SabatiniS. E.RodríguezM. C. (2012). Isocitrate lyase activity and antioxidant responses in copper-stressed cultures of Chlamydomonas reinhardtii (Volvocales, Chlorophyceae). Phycologia51, 135–143. 10.2216/10-56.1

  • CrossRef
  • Google Scholar

• 97

  TakouridisS. J.TribeD. E.GrasS. L.MartinG. J. (2015). The selective breeding of the freshwater microalga Chlamydomonas reinhardtii for growth in salinity. Bioresour. Technol.184, 18–22. 10.1016/j.biortech.2014.10.120

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 98

  TangeT. Ø.NottA.MooreM. J. (2004). The ever-increasing complexities of the exon junction complex. Curr. Opin. Cell Biol.16, 279–284. 10.1016/j.ceb.2004.03.012

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 99

  ValášekL. S. (2012). ‘Ribozoomin’–translation initiation from the perspective of the ribosome-bound eukaryotic initiation factors (eIFs). Curr. Protein Peptide Sci.13, 305–330. 10.2174/138920312801619385

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 100

  VizcainoJ. A.DeutschE. W.WangR.CsordasA.ReisingerF.RiosD.et al. (2014). ProteomeXchange provides globally coordinated proteomics data submission and dissemination. Nat. Biotechnol.32, 223–226. 10.1038/nbt.2839

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 101

  VoutsinaA.RivaM.CarlesC.AlexandrakiD. (1999). Sequence divergence of the RNA polymerase shared subunit ABC14. 5 (Rpb8) selectively affects RNA polymerase III assembly in Saccharomyces cerevisiae. Nucleic Acids Res.27, 1047–1055. 10.1093/nar/27.4.1047

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 102

  WangF.DongJ.JiangX.YeM.ZouH. (2007). Capillary trap column with strong cation-exchange monolith for automated shotgun proteome analysis. Anal. Chem.79, 6599–6606. 10.1021/ac070736f

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 103

  WangL.PiaoT.CaoM.QinT.HuangL.DengH.et al. (2013). Flagellar regeneration requires cytoplasmic microtubule depolymerization and kinesin-13. J. Cell Sci.126, 1531–1540. 10.1242/jcs.124255

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 104

  WangM. C.PengZ. Y.LiC. L.LiF.LiuC.XiaG. M. (2008). Proteomic analysis on a high salt tolerance introgression strain of Triticum aestivum/Thinopyrum ponticum. Proteomics8, 1470–1489. 10.1002/pmic.200700569

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 105

  WangS. B.ChenF.SommerfeldM.HuQ. (2004a). Proteomic analysis of molecular response to oxidative stress by the green alga Haematococcus pluvialis (Chlorophyceae). Planta220, 17–29. 10.1007/s00425-004-1323-5

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 106

  WangW.VinocurB.ShoseyovO.AltmanA. (2004b). Role of plant heat-shock proteins and molecular chaperones in the abiotic stress response. Trends Plant Sci.9, 244–252. 10.1016/j.tplants.2004.03.006

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 107

  WilsonN. F.LefebvreP. A. (2004). Regulation of flagellar assembly by glycogen synthase kinase 3 in Chlamydomonas reinhardtii. Eukaryot. Cell3, 1307–1319. 10.1128/EC.3.5.1307-1319.2004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 108

  XuJ.ZhangX.YeN.ZhengZ.MouS.DongM.et al. (2013). Activities of principal photosynthetic enzymes in green macroalga Ulva linza: functional implication of C(4) pathway in CO(2) assimilation. Sci. China Life Sci.56, 571–580. 10.1007/s11427-013-4489-x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 109

  XuanJ.SongY.ZhangH.LiuJ.GuoZ.HuaY. (2013). Comparative proteomic analysis of the stolon cold stress response between the C4 perennial grass species Zoysia japonica and Zoysia metrella. PLoS ONE8:e75705. 10.1371/journal.pone.0075705

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 110

  ZhangH.HanB.WangT.ChenS.LiH.ZhangY.et al. (2012). Mechanisms of plant salt response: insights from proteomics. J. Proteome Res.11, 49–67. 10.1021/pr200861w

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 111

  ZhangT.GongH.WenX.LuC. (2010). Salt stress induces a decrease in excitation energy transfer from phycobilisomes to photosystem II but an increase to photosystem I in the cyanobacterium Spirulina platensis. J. Plant Physiol.167, 951–958. 10.1016/j.jplph.2009.12.020

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 112

  ZhaoQ.SuoJ.ChenS.JinY.MaX.YinZ.et al. (2016). Na₂CO₃-responsive mechanisms in halophyte Puccinellia tenuiflora roots revealed by physiological and proteomic analyses. Sci. Rep.6:32717. 10.1038/srep32717

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 113

  ZhouH.YeM.DongJ.CorradiniE.CristobalA.HeckA. J.et al. (2013). Robust phosphoproteome enrichment using monodisperse microsphere-based immobilized titanium (IV) ion affinity chromatography. Nat. Protoc.8, 461–480. 10.1038/nprot.2013.010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 114

  ZhouH.YeM.DongJ.HanG.JiangX.WuR.et al. (2008). Specific phosphopeptide enrichment with immobilized titanium ion affinity chromatography adsorbent for phosphoproteome analysis. J. Proteome Res.7, 3957–3967. 10.1021/pr800223m

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 115

  ZouX.SunM.GuoX. (2006). Quantitative response of cell growth and polysaccharide biosynthesis by the medicinal mushroom Phellinus linteus to NaCl in the medium. World J. Microbiol. Biotechnol.22, 1129–1133. 10.1007/s11274-006-9153-1

  • CrossRef
  • Google Scholar