==== Front Mol Cell Proteomics Mol Cell Proteomics Molecular & Cellular Proteomics : MCP 1535-9476 1535-9484 American Society for Biochemistry and Molecular Biology S1535-9476(23)00093-2 10.1016/j.mcpro.2023.100582 100582 Research Parallel Proteomic Comparison of Mutants With Altered Carbon Metabolism Reveals Hik8 Regulation of PII Phosphorylation and Glycogen Accumulation in a Cyanobacterium Huang Chengcheng 12‡ Duan Xiaoxiao 12‡ Ge Haitao 1 Xiao Zhen 12 Zheng Limin 12 Wang Gaojie 12 Dong Jinghui 12 Wang Yan 12 Zhang Yuanya 1 Huang Xiahe 1 An Hongyu 12 Xu Wu 3 Wang Yingchun ycwang@genetics.ac.cn 12∗ 1 State Key Laboratory of Molecular Developmental Biology, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing, China 2 University of Chinese Academy of Sciences, Beijing, China 3 Department of Chemistry, University of Louisiana at Lafayette, Lafayette, Louisiana, USA ∗ For correspondence: Yingchun Wang ycwang@genetics.ac.cn ‡ These authors contributed equally to this work. 22 5 2023 7 2023 22 5 2023 22 7 10058224 11 2022 18 5 2023 © 2023 The Authors 2023 https://creativecommons.org/licenses/by/4.0/ This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). Carbon metabolism is central to photosynthetic organisms and involves the coordinated operation and regulation of numerous proteins. In cyanobacteria, proteins involved in carbon metabolism are regulated by multiple regulators including the RNA polymerase sigma factor SigE, the histidine kinases Hik8, Hik31 and its plasmid-borne paralog Slr6041, and the response regulator Rre37. To understand the specificity and the cross-talk of such regulations, we simultaneously and quantitatively compared the proteomes of the gene knockout mutants for the regulators. A number of proteins showing differential expression in one or more mutants were identified, including four proteins that are unanimously upregulated or downregulated in all five mutants. These represent the important nodes of the intricate and elegant regulatory network for carbon metabolism. Moreover, serine phosphorylation of PII, a key signaling protein sensing and regulating in vivo carbon/nitrogen (C/N) homeostasis through reversible phosphorylation, is massively increased with a concomitant significant decrease in glycogen content only in the hik8-knockout mutant, which also displays impaired dark viability. An unphosphorylatable PII S49A substitution restored the glycogen content and rescued the dark viability of the mutant. Together, our study not only establishes the quantitative relationship between the targets and the corresponding regulators and elucidated their specificity and cross-talk but also unveils that Hik8 regulates glycogen accumulation through negative regulation of PII phosphorylation, providing the first line of evidence that links the two-component system with PII-mediated signal transduction and implicates them in the regulation of carbon metabolism. Graphical Abstract Highlights • Comparative proteomics for Synechocystis mutants defective in carbon metabolism. • Identification of proteins and protein networks involved in carbon metabolism. • PII phosphorylation by Hik8 negatively regulates glycogen accumulation. In Brief Comparative proteomics was simultaneously performed for five knockout mutants of Synechocystis each depleted of an important regulatory protein for carbon metabolism (RCM). Proteins uniquely or commonly regulated by the RCM proteins were identified. Particularly, dramatic hyperphosphorylation of PII on S49 was also identified uniquely in the hik8-knockout mutant that is concomitant with significantly reduced glycogen content. Further experimental evidence suggests that Hik8 regulates glycogen accumulation through negative regulation of PII phosphorylation. Keywords proteomics Synechocystis carbon metabolism cyanobacteria phosphorylation Abbreviations AGC automatic gain control AT Autotrophic Chl Chlorophyll C/N Carbon/nitrogen DCMU 3-(3,4-dichlorophenyl)-1,1-dimethyl urea DEPs Differentially expressed proteins FDR false discovery rate GgpS glucosylglycerolphosphate synthase HCD high-energy collisional dissociation LAHG Light-activated heterotrophic growth MS Mass spectrometry MT Mixotrophic OPPP Oxidative pentose phosphorylation pathway PfkA1 phosphofructokinase PHT Photoheterotrophic PRM Parallel reaction monitoring Pyk1 pyruvate kinase RCM Regulatory protein for carbon metabolism RP-HPLC Reversed phase-high performance liquid chromatography Synechococcus Synechococcus elongatus PCC 7942 Synechocystis Synechocystis sp. PCC 6803 TMT Tandem mass tag WT Wild type ==== Body pmcCyanobacteria are oxygenic photoautotrophic organisms inhabiting a wide range of environments (1). They are key players in shaping the current global ecosystem by converting CO2 into a utilizable carbon resource while evolving oxygen from the water through photosynthesis. Because of their resemblance to the chloroplast in the higher plants, a number of cyanobacteria strains have been effectively used as the model for photosynthesis research. More recently, cyanobacteria gained more interest because of their potential as the unique chassis (green Escherichia coli) for synthetic biology and the cell factory to produce clean and renewable biofuels (2, 3). Photosynthesis and carbon metabolism are two physically and functionally interconnected processes in cyanobacteria. In addition, nitrogen assimilation and metabolism are also interwoven and coordinated with photosynthesis and carbon metabolism. The light reaction of photosynthesis occurs on the thylakoid membrane, whereas CO2 fixation occurs in a specialized non-membranous microcompartment formed by multiple proteins, namely, carboxysome (4). In contrast, central carbon metabolism and nitrogen assimilation mainly occur in the cytoplasm. Many cyanobacteria strains can grow photoautotrophically (AT) using CO2 as the sole carbon source or photoheterotrophically (PHT) using the exogenously supplied organic form of carbon such as glucose or photomixotrophically (MT) using both as the carbon source. Under AT growth mode, CO2 is fixed via the Calvin cycle and hexose-phosphates are produced through gluconeogenesis. These sugar products are subsequently converted into storage carbohydrates, such as glycogen (5). Endogenously synthesized or exogenously supplied glucose can also be catabolized through the oxidative pentose phosphorylation pathway (OPPP) and/or glycolysis to supply carbon and energy for the growth of cyanobacteria (6, 7). In cyanobacteria, a number of proteins in sugar catabolic pathways are indispensable for heterotrophic growth in dark. For example, the mutants deficient with glucose-6-phosphate dehydrogenase (Zwf) or glyceraldehyde 3-phosphate dehydrogenase 1 (Gap1) are defective in viability in the dark and unable to carry out light-activated heterotrophic growth (LAHG) (8). As the Calvin cycle, OPPP, and glycolysis take place within a single cellular compartment, many reversible reactions, enzymes, regulators, and intermediate metabolites are shared by the three processes. Moreover, the processes are physically and functionally connected with the tricarboxylic acid cycle (TCA) cycle, nitrogen assimilation, amino acid and protein metabolism, lipid metabolism, and many other metabolic processes. A single perturbation of one process may generate significant consequences on the activities and outcomes of the others. For example, the addition of glucose to a cyanobacterium culture can enhance the activity of OPPP while partly repressing the activity of photosynthesis (9). In addition, the gluconeogenic activity of the upper part of glycolysis enhances the flux of the OPPP to supply NADPH depleted in the dark (10). The metabolic pathways, particularly the central carbon metabolism that interconnects both with photosynthesis and the other metabolic pathways such as nitrogen assimilation, must be intricately and precisely regulated. Multiple proteins have been reported as important regulators of carbon metabolism in Synechocystis sp. PCC 6803 (hereafter referred to as Synechocystis), including Hik8 (Sll0750), Hik31 (Sll0790), and its plasmid paralog Slr6041, Rre37 (Sll1330), and SigE (Sll1689) (hereafter collectively named as the regulatory proteins for carbon metabolism, RCM) (11, 12, 13, 14, 15, 16, 17). The histidine kinases Hik8 and Hik31/Slr6041 are sensors of the two-component systems (18, 19). The ortholog of Hik8 (SasA) and its cognate response regulator RpaA in Synechococcus elongatus PCC 7942 (hereafter referred to as Synechococcus) are output components of the circadian clock and regulate global gene expression with circadian rhythmicity (20, 21). Rre37 is a response regulator of a two-component system (11), and SigE is a group 2 sigma factor of bacterial RNA polymerase (11, 12, 13, 14, 15, 16, 17, 22, 23). A plethora of proteins involved in sugar catabolic pathways such as OPPP and glycolysis are downregulated in the RCM-deficient cyanobacterial mutants, indicating that they are positively regulated by the RCM proteins. Nevertheless, for each individual protein in these pathways, its expression could be either commonly regulated by more than one RCM protein with the same or different extent of regulations or uniquely regulated by only one RCM protein. Such regulatory specificity and cross-talk may warrant the plasticity and robustness of carbon metabolism, which is necessary for maintaining metabolic and energy homeostasis to fine-tune cyanobacteria growth in response to environmental and metabolic cues. For example, nitrogen starvation can induce the expression of a number of proteins involved in sugar catabolism through the operation of one or more of the RCM proteins (24, 25, 26). Unfortunately, since the studies regarding the RCM proteins were conducted by different groups, at different time, and with different experimental conditions (15, 17, 18, 27), it is difficult to quantitatively and qualitatively evaluate such regulatory specificity and the cross-talk. The scenario may become even more complicated when concomitant differential expression of proteins in other metabolic pathways is involved. It is also important to simultaneously identify such proteins and quantitatively determine the extent of their differential expression in the RCM-deficient mutants. Herein, we tried to address this problem by quantitatively and simultaneously identifying the proteins regulated by the five RCM proteins in the model organism Synechocystis. Synechocystis is one of the most extensively studied model cyanobacteria for photosynthesis (28, 29). It is the first cyanobacterium with a completely sequenced genome and is highly transformable (30). We generated knockout mutants for the five RCM protein-coding genes, and then quantitatively and simultaneously compared the proteomes across the five knockout mutants using a tandem-mass tag (TMT) labeling-based quantitative proteomics approach. Proteins differentially expressed in one or more of the mutants were identified and quantitatively compared. Novel insights into the regulation of carbon metabolism were provided from the result. Experimental Procedures Antibodies The primary antibodies for Gnd, Pgk, KaiC1, PntA, Zwf, and GltB were provided by PhytoAB. The antibodies for AtpB, ClpB2, Prk, PsaD, and SbtA were purchased from Agrisera. Cell Culture The wild type (WT) and mutant strains of Synechocystis were cultured in liquid BG-11 medium in medium light (50 μmol m−2 s−1 photons) with shaking at 30 °C. The cells were collected by centrifugation (4000g for 10 min) for biochemical and proteomic analyses when the culture reaches the optical density at 730 nm (OD730) of approximately 1.0. The optical density measures the turbidity of the cell culture (31), which positively correlates with cell density. The harvested cells were stored at −80 °C until use. For MT growth, 5 mM glucose was supplemented to the medium. For PHT growth, cells were cultured in the same condition as that used for MT growth except that 5 μM DCMU was added to the culture medium. For LAHG growth, the cells were grown in darkness in the presence of 5 mM glucose and a daily pulse of white light (at least 5 min). The starting concentrations of cells at OD730 of 0.05 and 0.2 were used for growth curve measurements under the phototrophic and the LAHG conditions, respectively. Mutants Generation The gene knockout mutants were generated through insertional mutation and homologous recombination using the procedures as we described (13, 32). To generate the strains carrying PII S49A mutation, a glnB (PII coding gene) knockout strain of Synechocystis was first generated as above described (ΔglnB). A DNA fragment containing wild-type glnB and its native promoter or a mutant copy carrying PII S49A substitution was inserted in the chromosome of the ΔglnB mutant at a site designated “neutral site” to generate the WT equivalent (WTe) strain or glnBS49A/WTe strain (33, 34). The insertional mutation of hik8 in the glnBS49A/WTe background was subsequently performed as above described to generate the mutant glnBS49A/WTe/Δhik8. All mutants were selected and confirmed as described earlier, and all primers used were included in supplemental Table S1. Pigment Analysis The Chl and total carotenoids in cells were extracted with N, N-dimethylformamide (DMF, Sigma-Aldrich, Saint Louis, MO), and the concentrations were determined using a previously described approach (35). The equations used for the calculation of pigment concentration are:Chla(μg/mL)=12.1×OD664−0.17×OD625. TotalCarotenoids(μg/mL)=(OD461−0.046×OD664)×4. 77K Fluorescence Measurement 77K fluorescence emission spectra were measured using an F-7000 Fluorescence Spectrophotometer (Hitachi). Cells were adjusted to a concentration of 15 μg Chl mL−1, and the fluorescence emission spectra (excitation at 435 nm, bandwidth of 5 nm) were recorded in the range of 650 to 800 nm and were normalized at 726 nm. Glycogen Content Determination Glycogen was determined as described (36). For each measurement, 20 mL of exponentially growing Synechocystis cells (OD730 = 1) were harvested, resuspended in 400 μL 30% (w/v) KOH, and incubated at 95 °C for 2 h. Glycogen was then precipitated with 75% (v/v) ice-cold ethanol and collected by centrifugation at 10,000g for 10 min. The pellet was sequentially washed with 70% and 98% (v/v) ethanol and dried at 60 °C for 10 min. The isolated glycogen was resuspended in 100 mM sodium acetate (pH 4.5) and enzymatically hydrolyzed to glucose with 2 mg/mL amyloglucosidase (Sigma-Aldrich) at 60 °C for 2 h. The glycogen content was determined with a glucose assay kit (Sigma-Aldrich) according to the manufacturer’s instructions. Protein Preparation Cell pellets were lysed in a buffer containing 0.4 M sucrose, 50 mM MOPS, 10 mM NaCl, 5 mM EDTA (pH 7.0), and 0.5 mM PMSF with a bead beater, and the insoluble debris was removed by centrifugation for 30 min at 5000g at 4 °C. After precipitation with ice-cold 10% trichloroacetic acid in acetone at −20 °C, total proteins were washed with acetone and resolubilized with 4% sodium dodecyl sulfate (SDS) in 0.1 M Tris-HCl, pH 7.6. A BCA protein assay kit (Thermo Scientific) was used to determine the protein concentration. Protein Digestion and TMT Labeling Proteins were digested using the filter-aided sample preparation (FASP) method according to a previously described method with slight modifications (32, 37). Briefly, the lysates (100 μg protein for each sample) were reduced with 10 mM DTT at 37 °C for 1 h and alkylated with 55 mM iodoacetamide (IAA, Sigma-Aldrich) at room temperature for 1 h in the dark. The alkylated lysates were transferred into the Microcon YM-30 centrifugal filter units (EMD Millipore Corporation), where the denaturing buffer was replaced by the 0.1 M triethylammonium bicarbonate (TEAB, Sigma-Aldrich, Saint Louis, MO), and then digested with sequencing grade trypsin (Promega) at 37 °C overnight. The resulting tryptic peptides were collected and labeled with 6-plex TMT reagents (Thermo Scientific) by incubating peptides with ethanol-dissolved TMT reagents for 2 h at room temperature in the dark. The labeling reaction was inactivated by the addition of 5% hydroxylamine, and the labeled samples were mixed together with equal ratios before being fractionated with reversed-phase (RP) high-performance liquid chromatography (HPLC). RP-HPLC and Desalting Offline basic RP-HPLC was performed using a Waters e2695 separations HPLC system coupled with a Phenomenex gemini-NX 5 μ C18 column (250 × 3.0 mm, 110 Å). The sample was separated with a 97 min basic RP-LC gradient as previously described (38). A flow rate of 0.4 mL/min was used for the entire LC separation. The separated samples were collected into 14 fractions, completely dried with a Speed-Vac concentrator, and stored at −20 °C. The dried peptides were resolubilized with 0.5% acetic acid and desalted using C18 Stage Tips (39). Desalted peptides were dried with a Speed-Vac concentrator, stored at −20 °C, and resuspended in 0.1% formic acid (FA) immediately before LC-MS/MS. MS Analysis For MS analysis, the peptides were resuspended in 0.1% FA and analyzed by a LTQ Orbitrap Elite mass spectrometer (Thermo Scientific) coupled online to an Easy-nLC 1000 in the data-dependent mode. Briefly, 2 μL of peptide sample (1 μg/μL) was injected into a 25-cm length, 75-μm inner diameter capillary analytic column packed with C18 particles of 5-μm diameter (SunChrom, Friedrichsdorf, Germany). The mobile phases for the LC include buffer A (2% acetonitrile, 0.1% FA) and buffer B (98% acetonitrile, 0.1% FA). The peptides were separated using a 90-min non-linear gradient consisting of 3% to 8% B for 10 min, 8% to 20% B for 60 min, 20% to 30% B for 8 min, 30% to 100% B for 2 min, and 100% B for 10 min at a flow rate of 300 nL/min. The source voltage and current were set at 2.5 KV and 100 μA, respectively. All MS measurements were performed in the positive ion mode. The precursors were measured by survey scans in the Orbitrap with a mass range of 300 to 1800 m/z at a resolution of 120,000 at m/z 400. The 15 most abundant precursor ions (top15) from each MS scan were isolated and fragmented by high-energy collisional dissociation (HCD) with 40% collision energy for MS/MS analysis using DDA mode. The automatic gain control (AGC) and maximum injection time were set at 1 × 106 and 200 ms, respectively. The MS/MS spectra were acquired at a resolution of 15,000 at m/z 400 in the Orbitrap. The duration of dynamic exclusion was set to 30 s to prevent repeat identification of peptide ions within the time duration. The PRM experiment was performed as described (40). PRM and DDA for the PRM spectral library generation were performed using Orbitrap Lumos (Thermo Fisher Scientific). The full scan mass range was 300 to 1500 m/z and the Orbitrap resolution was 120,000 at 200 m/z. The AGC target and maximum injection time were set at 1 × 106 and 30 ms, respectively. For DDA, the mass spectrometer was operated in data-dependent acquisition mode with a maximum duty cycle time of 3 s, the isolation window was 1.6 m/z, and the MS/MS spectra were acquired at a resolution of 30,000 at 200 m/z in the Orbitrap. The duration of dynamic exclusion was set to 30 s. For PRM, target precursor ions (target peptide list) (supplemental Table S6) were fragmented using HCD fragmentation with 32% collision energy and detected with Orbitrap at a mass resolution of 30,000 at 200 m/z, the isolation window was 1.4 m/z. Database Search The raw MS files were searched against the Synechocystis proteome sequence database using the software MaxQuant (version 1.5.4.1) (41). The database containing 3672 entries was downloaded from the CyanoBase (ftp://ftp.kazusa.or.jp/pub/CyanoBase/Synechocystis, released on 5/11/2009). The type of search was set to report ion MS2 and the 6-plex TMT was chosen, the minimum reporter parent ion interference was set to 0.75. Trypsin was chosen as the protease for protein digestion, and the maximum of 2 was set as the allowable miscleavages. N-terminal acetylation and methionine oxidation were chosen as the variable modifications and cysteine carbamidomethylation was chosen as the fixed modification. The mass tolerances were set to 4.5 ppm for the main search and 20 ppm for precursor and fragment ions. The minimum score for unmodified peptides and modified peptides were set to 15 and 40, respectively. Other parameters were set up using the default values. The false discovery rate (FDR) was set to 0.01 for both peptide and protein identifications. Experimental Design and Statistical Rationale Triplicated TMT labeling-based quantitative proteomic experiments were performed to analyze the differentially expressed proteins among the WT Synechocystis and the five RCM mutants, thereby covering three biological replicates of each sample. In each experiment, tryptic peptides from each sample were labeled with a distinct TMT reagent of the 6-plex TMT reagent in an alternating order as shown in Figure 1B to reduce quantitative bias (32).Fig. 1 Quantitative identification of Synechocystis proteomes among the WT and regulatory proteins for carbon metabolism (RCM) mutant strains under photoautotrophic (AT) condition.A, Schematic representation of the strategy for quantitative analysis of the Synechocystis proteomes using the TMT-labeling-based quantitative proteomics approach. Three biological replicates were included. B, The scheme of TMT labeling. Note that the orders of the TMT labeling are different among the replicates. C, Venn diagram shows the numbers of the total identified and quantified proteins. Bioinformatic and statistical analyses were mainly performed using the software Perseus (version 1.5.4.1) (42). Student’s t test was used to determine the significance of differential expression of proteins, and Fisher’s-exact test was used for the functional enrichment analysis. A p value < 0.05 was used as the cutoff for all statistical analyses. For the PRM experiment, four target peptides were selected (supplemental Table S6), which were the phosphopeptides bearing S49 phosphorylation of the protein PII and their non-phosphorylated counterparts (including all peptides resulting from miscleavages). The non-phosphorylated counterparts were used as the internal control because they are expected to show the opposite direction of change in abundances to that of the corresponding phosphopeptides, given that the abundance of PII is not changed among the samples to be compared (supplemental Figs. S9–S13). Three biological replicates were included to quantify the target peptides in AT/MT-grown WT and Δhik8. PRM data (Tier 2 level) were processed with Skyline (version 21.2.0.369) software (43). The DDA library was generated using msms.txt from MaxQuant (version 1.6.0.16) (44). A target peptide with a higher DOPT value (at least 0.8 in an experiment) was selected if the target peptide has two different charge states. The top 6 transitions of the target peptides (supplemental Table S7) were selected, peak picking was performed and manually inspected, and the peak area of the transitions was summed for quantification. Student’s t test was used to determine the significance of the abundance difference of the target peptides. Standard deviation was used to confirm the precision of target peptides. Results Generation and the Growth Experiment of the Knockout Mutants of Genes Encoding the Five RCM Proteins The knockout mutants were generated individually by an insertional mutation of each of the RCM-coding genes. The full segregation of the knockout mutants was confirmed by PCR detection of the wild-type (WT) DNA fragments of the corresponding genes, which should be completely undetectable in the mutants (supplemental Fig. S1) (13). The WT and the mutants were cultured in the AT, MT, PHT, and LAHG conditions and the growth curves were plotted (supplemental Fig. S2). Under the AT and MT conditions, no significant growth phenotype was observed for all the mutants. Under the PHT conditions, the growth of ΔsigE was severely inhibited, and the growth of Δhik8 was slightly affected as the culture became more yellowish. For all the other mutants, no significant difference of growth from that of the WT was observed. Under the LAHG conditions, the growth of ΔsigE and Δhik8 was nearly completely inhibited, and the growth of Δrre37 was also severely impaired but to a lesser extent compared to that of ΔsigE and Δhik8. In contrast, Δhik31 and Δslr6041 grew as nearly well as the WT. The growth phenotypes of all the mutants are consistent with existing reports (11, 18, 27, 45). Note that individual inactivation of hik31 or slr6041 did not cause an obvious growth defect, but inactivation of both genes was reported to impair the growth with D-glucose (46), suggesting that the two genes may be functionally redundant, at least partly. For the cells cultured in the AT mode, pigment measurement demonstrates that the chlorophyll (Chl) and carotenoid content were not significantly changed in all the mutants compared with the WT (supplemental Fig. S2E). Measurement of 77K fluorescence shows only minor differences for the total amount and the stoichiometry of photosystem I and II among the WT and the mutants (supplemental Fig. S2F). Together, the results of the phenotypic analyses indicate that all the mutants grew nearly equally well under the AT conditions. Since differential protein expression in mutants could be induced directly by mutation or indirectly by secondary factors such as growth defect. All the cells used for proteomic analysis were subsequently cultured under the AT conditions to minimize the potential growth-defect induced differential protein expression. Quantitative Proteomic Identification of Differentially Expressed Proteins Across the RCM-Deficient Mutants The proteomes of the WT and the five mutant cells were analyzed using a 6-plex TMT-labeling-based quantitative proteomics approach, which allows simultaneous quantification of the relative abundances of proteins in up to six independent samples (47). The overall strategy of the proteomic analysis is illustrated (Fig. 1A), reshuffling the order of TMT labeling between replicates was employed to avoid quantitative bias from each TMT reagent (Fig. 1B) (13). In total, 2543 proteins were identified with 2189 proteins containing quantitative TMT information in at least two replicates (Fig. 1C), and the latter represents the selected subset of proteins for further statistical and bioinformatics analyses (supplemental Table S2). The overall reproducibility of the TMT quantitation is high, as indicated by the high correlation efficiency between any two of the biological replicates (R2 > 0.91) (supplemental Fig. S3). The fold-changes of each protein in all the mutants relative to the WT were calculated and Z-scored, and used for a clustering analysis. All the replicates from each of the samples were correctly clustered (Fig. 2A), further demonstrating high reproducibility of the quantification. Student’s t test with a threshold p < 0.05 and a fold-change threshold were used to sequentially filter the DEPs (32, 37, 40, 48). To determine the fold-change threshold, the TMT ratios (mutant/WT) for each mutant were firstly calculated for all proteins in all three replicates. The mean of the TMT ratios for each replicate was then calculated, which is expected to be close to 1 because samples in TMT channels were mixed with an equal ratio. The TMT ratios in each replicate were then normalized by the corresponding mean and the distribution of the resulting normalized ratios was calculated and plotted for all replicates of all mutants (supplemental Fig. S4). The result shows that at least 95.5% of the normalized ratios are smaller than 1.3, suggesting that a fold change of 1.3 is a reasonable threshold for quantitation with an estimated FDR smaller than 5%. In addition, due to the well-known precursor ion interference that usually causes a compressed ratio in isobaric labeling-based quantitative proteomics (49, 50), a protein with a measured 1.3-fold change in abundance may have a greater extent of differential expression. By applying both Student's t test p < 0.05 and a fold change of 1.3, 40, 116, 49, 78, and 129 proteins were identified as differentially expressed proteins (DEPs)in Δhik31, Δhik8, Δrre37, ΔsigE, and Δslr6041, respectively (Fig. 2B). A few DEPs were further validated through immunoblotting (Fig. 2C), and all these proteins were previously implicated in carbon metabolism except KaiC1 (7, 15, 27, 51, 52, 53, 54). The circadian clock protein KaiC1 and its ortholog in Synechococcus interact directly with Hik8 and its ortholog SasA, respectively (20, 26, 55), and the interaction is important for the regulation of circadian rhythmicity (20). Meanwhile, immunoblotting for other proteins without significant TMT-measured changes consistently show only slight variations in abundance among all the RCM mutants (supplemental Fig. S5). For Δhik31 and ΔsigE, DNA microarray-based transcriptomic data from similar photoautotrophic conditions is available in reported studies and was compared with our proteomics data, and the majority of overlapping genes show a similar trend of differential expression at both mRNA and protein levels (supplemental Fig. S6). Only a few genes show opposite directions of differential expression at the two levels, this is not surprising because it is well accepted that mRNA and protein levels are not always positively correlated (56). Similar comparisons were also performed for Δhik8 and Δrre37. The two mutants do not have either a full list of transcriptomics data (Δhik8) or the transcriptomics data obtained under the photoautotrophic growth conditions (Δrre37). Therefore, we manually collected the transcription data of individual genes from the literature and performed the comparisons (11, 17, 27, 57). The results are similar to those of Δhik31 and ΔsigE (supplemental Fig. S6).Fig. 2 Determination of differentially expressed proteins (DEPs) among the RCM mutants.A, Clustering analysis of the z-scored TMT reporting ion intensities to evaluate the reproducibility of the quantitative analysis. B, Distribution of fold changes of DEPs in the RCM mutants. Student’s t test p value (p < 0.05) and fold change threshold (FC>1.3) were used as the thresholds to filter for the DEPs. C, Western blotting confirmation of indicated DEPs, the log2 transformed fold-changes measured by TMT were shown to the right side of each panel. Ponceau staining was used as the loading control. The number of DEPs overlapping among the RCM mutants or specific to a particular mutant is shown by the Venn diagram (Fig. 3A). In total, only 4 DEPs were overlapping across all five mutants. These include Sll7085, Slr6074 and Slr1161 that were downregulated and glucosylglycerolphosphate synthase (GgpS) that was upregulated in all mutants (Fig. 3B). For the mutant-specific DEPs, the numbers are 11, 60, 14, 25, and 97 respectively for Δhik31, Δhik8, Δrre37, ΔsigE and Δslr6041, representing the list of proteins uniquely regulated by the corresponding RCM protein (supplemental Table S3). For better visualization of proteins regulated by at least two of the RCM proteins, a regulatory network for 80 proteins with a significant change in at least two mutant strains was constructed using the five RCM proteins as the hubs (Fig. 3C and supplemental Table S4). The network clearly demonstrates the cross-talk of regulation that might be important for Synechocystis to grow and cope with metabolic and environmental changes.Fig. 3 The differentially expressed proteins (DEPs) commonly and specifically regulated by RCM proteins.A, Venn diagram showing the DEPs commonly and specifically regulated by the RCM proteins (FC>1.3, p < 0.05). B, DEPs commonly regulated by all the five RCM proteins. C, the regulatory network for 80 DEPs showing differential expression in at least two RCM mutants. The RCM proteins were included as the hub and highlighted in yellow. Commonly and Differentially Regulated Functions in the Mutants To gain insight into the functions that are commonly or differentially regulated by the five RCM proteins, the enrichment of the Gene Ontology (GO) terms was assayed by Fisher’s-exact test for the DEPs in each RCM mutant. Among the upregulated proteins, two GO terms, the trehalose biosynthesis process, and the glucosylglycerol biosynthetic process, are enriched in all mutants (Fig. 4A) suggesting that all the five RCM proteins are involved in the negative regulation of the two processes. A few other GO terms are enriched in three, two, or only one mutant, indicative of the common or unique function regulated by the corresponding RCM proteins. Similarly, among the downregulated proteins, the GO term response to the biotic stimulus is enriched in four out of the five mutants, with the exception of the mutant Δslr6041. The other GO terms are enriched in three, two or only one mutant (Fig. 4B).Fig. 4 Enriched functions in the differentially expressed proteins (DEPs) across the five RCM mutants. Fisher’s-exact test for the functional enrichment in the upregulated (A) and downregulated (B) proteins across the RCM mutants. The gene ontology term biological process (GOBP) was used for the enrichment assay, p value < 0.05 and enrichment factor >1.5 were used as the thresholds. It is expected that the functions related to carbon metabolism are impacted in the RCM mutants. Indeed, a number of GO terms associated with carbon metabolism were enriched in the DEPs. For the upregulated proteins, the nonoxidative branch of the pentose phosphate shunt is enriched in Δhik31, Δrre37, and ΔsigE. For the downregulated proteins, carbohydrate metabolic process, glucose metabolic process, pentose phosphate shunt, fructose-phosphate metabolic process, glycogen catabolic process, mannose metabolic process, and glycolysis are enriched in up to two mutants. The result further illustrates the specificity and cross-talk of the regulations on carbon metabolism by the RCM proteins. Remarkably, the deletion of hik31 or its plasmid-borne copy slr6041 differentially impacts a number of functions (Fig. 4), though the two proteins share more than 99% sequence similarity. The result is consistent with a previous report (18), and suggests that the two proteins are not completely redundant, though they can partially complement each other. For better visualization and functional categorization of the DEPs, all DEPs with known functional categories of the CyanoBase were grouped and displayed in a heat map that shows the fold-changes of all DEPs in all mutants (Fig. 5 and supplemental Table S5). Again, common and differential regulation of the functions by the five RCM proteins are evident. The overall expression patterns of the DEPs are more similar between ΔsigE and Δhik8 and between Δhik31 and Δslr6041. In Δrre37, the expression patterns of some DEPs such as ribosomal proteins are more similar to that in ΔsigE and Δhik8, whereas some others are more similar to that in Δhik31 and Δslr6041 such as lipoproteins and porins. In contrast, significant downregulation of lipoproteins and porins (Slr1908, Slr1841, Slr1272, and Sll0772) occurs in Δrre37, Δhik31, and Δslr6041, but not in the other two mutants, although they display the trend of downregulation more or less in all mutants. It is evident that the numbers of DEPs in Δhik8, ΔsigE, and Δslr6041 are greater than those in Δrre37 and Δhik31. Notably, several ribosomal proteins were downregulated significantly in Δhik8, including Rpl29, Rps14, Rpl13, Rpl17, Rpl23, Rpl10, Rps10, and Rps7. The result is in line with a previous microarray result (27).Fig. 5 Heatmap visualization of the differentially expressed proteins (DEPs) assigned with known functional categories. Only DEPs in the known functional categories according to CyanoBase were included, and the DEPs in the categories hypothetical and unknown were not included. The color scale shows the abundance of each protein in each RCM mutant relative to that in the WT: blue indicates low level, red indicates high level, and white indicates no change. ∗p value < 0.05. For DEPs involved in carbon metabolism, 6 DEPs are involved in the pentose phosphate pathway, and their overall downregulation was apparent in ΔsigE and Δhik8, but not in the other three mutants. Particularly, transaldolase (Tal) and 6-phosphogluconate dehydrogenase (Gnd) were downregulated in both ΔsigE and Δhik8; Zwf, glucose 6-phosphate dehydrogenase assembly protein (OpcA), and pentose-5-phosphate-3-epimerase were significantly downregulated in ΔsigE, and they were also downregulated more or less in Δhik8 but not statistically significant. The 6-phosphogluconolactonase was downregulated only in Δhik8. Four DEPs are involved in glycolysis, and their overall downregulation in ΔsigE, Δhik8 and Δrre37 is apparent except for phosphoglycerate kinase, which was upregulated in Δhik8. Particularly, Gap1 and phosphofructokinase (PfkA1) were downregulated significantly in all three mutants and show the trend of weak downregulation in the other two mutants, pyruvate kinase (Pyk1) was downregulated significantly in ΔsigE and Δhik8 but not in Δrre37. Overall, the differential expression pattern of the proteins involved in sugar catabolism is consistent with previous reports (15, 27). It is noticeable that no significant downregulation of the aforementioned DEPs occurs in Δhik31 or Δslr6041, possibly due to the functional redundancy of Hik31 and Slr6041 (18). Of all DEPs involved in sugar catabolism, Gap1 is important for glycolysis and Zwf catalyzes the first and the rate-limiting step of the OPPP (7, 54). Depletion of Gap1 and Zwf was reported to impair the viability of Synechocystis under the LAHG condition (58, 59), and the extents of downregulation of the two proteins in Δhik8, ΔsigE, and Δrre37 correlate well with the extents of the reduced viability of the mutants in LAHG condition (supplemental Fig. S2D). Together, these results suggest that the decreased viability of Δhik8, ΔsigE, and Δrre37 could mainly be attributed to the downregulation of Zwf and Gap1. Decreased Glycogen Content and Upregulated PII Phosphorylation in AT-Grown Δhik8 Glycogen is a major carbon and energy storage compound in cyanobacteria (60). Glycogen catabolism and cellular respiration are critical for cell viability in darkness (26, 61, 62, 63). Depletion of the RCM proteins could affect glycogen accumulation. Measurement of the glycogen content reveals that it was significantly decreased in Δhik8 as previously reported (27), but not in the other mutants (Fig. 6A). Remarkably, the proteins involved in glycogen synthesis were not significantly downregulated in Δhik8 (supplemental Table S5), suggesting there might exist a novel mechanism underlying Hik8-dependent glycogen accumulation. The accumulation of glycogen was reported to be regulated by the intracellular C/N balance, which is in turn sensed and regulated by the phosphorylation status of the protein PII in cyanobacteria (64, 65, 66). PII acts as a multitasking signal-integrating regulator whose phosphorylation state mirrors the cellular ATP and 2-oxoglutarate (2-OG) levels (67). The canonical PII S49 phosphorylation is one of the most dominant phosphorylation events in Synechocystis and can be easily detected by liquid chromatography-tandem mass spectrometry (LC-MS) (68, 69). To reveal the PII S49 phosphorylation status in the RCM mutants, we re-searched the raw mass spectrometry files against the Synechocystis proteome database including phosphorylation as the variable modification, and quantitatively identified with high confidence a tryptic peptide bearing the PII S49 phosphosite in all replicates (Fig. 6, B and C). Importantly, the PII S49 phosphorylation level, as quantified by TMT, was upregulated more than 6-fold in Δhik8 in comparison with that in the WT, and also upregulated in Δslr6041 but to a much less extent (Fig. 6B). In all other mutants, PII S49 phosphorylation level was not significantly different from that in the WT (Fig. 6B). To confirm the upregulation of PII S49 phosphorylation and to ask whether the high level of phosphorylation is constitutive or regulable in Δhik8, we specifically quantified the abundances of the phosphopeptide in Figure 6C in AT/MT-grown WT and Δhik8 cells using a parallel reaction monitoring (PRM)-based target proteomics approach and calculated the S49 phosphorylation occupancy (40). The PRM result confirms the significant upregulation of PII S49 phosphorylation in AT-grown Δhik8, which is a nearly 26-fold increase in abundance compared with that in the WT (Fig. 6D and supplemental Fig. S7). Notably, the fold increase quantified by PRM is much greater than quantified by TMT due to the ratio compression effect of the isobaric labeling-based approach (70, 71). Nevertheless, both methods corroborate the strong increase of PII S49 phosphorylation in AT-grown cells. In the MT-grown Δhik8, the S49 phosphorylation level significantly decreased compared with that in AT-grown Δhik8 (Fig. 6D and supplemental Fig. S8), whereas in the MT-grown WT cells, the PII S49 phosphorylation was significantly upregulated to a level nearly commensurate with that in the AT-grown Δhik8, and the upregulation is consistent with our previous measurement using a TMT-labeling-based proteomics approach (32). PII S49 phosphorylation occupancy was also calculated for the cells grown under both conditions, the result revealed that nearly 70% of PII molecules are phosphorylated in Δhik8, whereas less than 4% are phosphorylated in the WT under the AT conditions (Fig. 6E, supplemental Table S8). Together, the PRM results suggest that the PII S49 phosphorylation in Δhik8 is regulable, at least through the change of nutrition mode.Fig. 6 Concomitant upregulation of PIIphosphorylation and decrease of glycogen content in Δhik8.A, Determination of glycogen contents among the WT and RCM mutants under AT condition. The cells were cultured under medium light intensity. ∗∗∗: p value < 0.001. B, quantitative identification of PII S49 phosphorylation across the RCM mutants. Bars show the median abundances of a phosphopeptide from PII bearing the S49 phosphosite measured from all replicates of the corresponding RCM mutants. Error bars represent standard deviations. ∗∗p value < 0.01. C, a representative mass spectrum of the phosphopeptide indicated in (B). The region containing TMT report ions (boxed by dashed lines) is horizontally zoomed in and displayed in the inset to show the relative abundance of the phosphopeptide in each sample. D, Bar graph shows the peak area of the phosphopeptides in (C) in the WT and Δhik8 cultured under AT and mixotrophic (MT) conditions. The measurement was performed using parallel reaction monitoring (PRM)-based target proteomics approach. ∗∗∗∗p value < 0.0001, ∗∗p value < 0.01. E, Calculation of PII S49 phosphosite occupancy in the indicated Synechocystis strains and growth conditions (left panel). The calculation was performed using a previously described method (106). The phosphopeptide YRGSpEYTVEFLQK and its non-phosphorylated counterparts (GSEYTVEFLQK and YRGSEYTVEFLQK) were used for quantitative PRM measurement (supplemental Tables S6–S8), the total PII level was detected by Western blotting and the quantitation was performed using the software ImageJ (right panel). ∗∗p value < 0.01, ∗∗∗∗p value < 0.0001. PII S49A Mutation Restores Glycogen Accumulation and Dark Viability of Δhik8 To investigate whether there is a causal relationship between the upregulation of PII S49 phosphorylation and the decrease of the glycogen content in AT-grown Δhik8, we introduced a site-specific mutation on glnB, the coding gene of PII, to introduce a non-phosphorylatable S49A mutation in PII in the background of both Δhik8 and the WT equivalent strains (WTe) (supplemental Fig. S14 and Experimental procedures). Glycogen determination shows that there is no significant difference in glycogen content between the two strains, and the glycogen content in glnBS49A/WTe/Δhik8 was dramatically increased to a level comparable to that in the WT, and the level is at least two-fold higher than that in Δhik8 (Fig. 7A). These data indicate PII S49A mutation can successfully recover the glycogen content in Δhik8 but with little effect on the glycogen accumulation in the WT. Together, these results strongly suggest that hyperphosphorylation of PII in Δhik8 accounts for the reduced glycogen content.Fig. 7 PIIS49A substitution restored the glycogen content and rescued the dark viability of Δhik8.A, Bars show the glycogen content measured for the indicated strains. All strains were cultured under AT conditions. Error bars: standard deviations. NS: not significant. B, Left: growth curves of the indicated strains under the light-activated heterotrophic growth (LAHG) condition. Right: The images of the cultures photographed at the indicated time points to show the color changes of the cells during the time course. C, After 96 h LAHG growth, the indicated strains (also shown in the red box in (B)) were transferred to continuous light and illuminated for another 48 h. D, A working model for Hik8 regulation of glycogen accumulation and dark viability through negative regulation of PII phosphorylation. To further ask whether PII S49A mutation can rescue Δhik8 in LAHG, a lethal condition for Δhik8 as observed by us and others (supplemental Fig. S2D) (27), the growth experiment was performed for the WT and the mutants in the LAHG condition. The growth curve shows that glnBS49A/WTe/Δhik8 can grow in LAHG (Fig. 7B). The concentration of the cells continued to increase during the time course of culturing, though the growth rate was slower than that of the WT and the color of the cells were yellowish. In contrast, the growth of Δhik8 was completely stopped after 24 h incubation in LAHG and the color of the culture became pale white. To further confirm the viability of the mutants, Δhik8 and glnBS49A/WTe/Δhik8 were re-illuminated with light after 96 h incubation in the dark. The culture of glnBS49A/WTe/Δhik8 became fresh green after 48 h incubation in light (Fig. 7C), indicating that the mutant cells survived the 96 h dark incubation. In contrast, the color of the Δhik8 culture did not change, indicative of cell death. Noticeably, the growth of glnBS49A/WTe was also significantly repressed as indicated by the slower growth and yellowish color of the culture, though the repression was not as severe as observed for the other two mutants (Fig. 7B). The observation suggests that PII S49A mutation cause pleiotropic effects that are detrimental for the WT but beneficial for Δhik8. This also explains that the PII S49A mutation did not restore the growth of Δhik8 to the comparable level of that of the WT. Based on these results, we proposed a model to depict the Hik8-regulated PII phosphorylation, glycogen accumulation, and dark viability (Fig. 7D). Shown in the model, regulation of the circadian clock (KaiC1-Hik8-RpaA) maintains PII phosphorylation at a low level in the photoautotrophically growing WT Synechocystis. Depletion of Hik8 leads to hyperphosphorylation of PII, resulting in the decrease of the glycogen content and the dark viability of Synechocystis. In addition to Hik8 and RpaA, KaiC1 is the only other circadian clock protein shown in the model, because KaiC1 directly interacts with Hik8 and is the only upregulated circadian oscillator protein in Δhik8 (supplemental Fig. S17) (20, 26, 55). KaiC1 is not a typical serine/threonine kinase but was reported to have autokinase and autophosphatase activities (72). Therefore, it is also possible that upregulated KaiC1 directly hyperphosphorylates PII in Δhik8 (Fig. 7D). In addition to RpaA, other unknown factors could also be involved to regulate PII phosphorylation in downstream of Hik8 (Fig. 7D). Discussion Cyanobacteria coordinately regulate diverse types of cellular metabolism in response to environmental and internal changes. The regulatory schemes that control carbon metabolism are very complex due to the large number of proteins, metabolites, and pathways involved. For more robust and precise regulation in response to the changing environmental and internal conditions, it is desirable that a subset of proteins in a pathway are commonly regulated by multiple proteins or a subset of proteins in different pathways are exclusively regulated by a single protein. For carbon metabolism, a few proteins have been implicated in such an intricate and elegant regulation, including SigE, Hik8, Rre37, Hik31, and Slr6041(11–17). However, the identities of the common and separate targets regulated by these RCM proteins have not been comprehensively investigated. In this study, a 6-plex TMT-labeling-based quantitative proteomics strategy was explored to simultaneously compare the proteomes of the mutants each depleted of one of the RCM proteins. Common and separate targets of the RCM proteins were identified (Fig. 3C and supplemental Table S3), functions that are commonly regulated by multiple RCM proteins or specifically regulated by a single RCM protein were uncovered (Fig. 4). The results provide a holistic perspective of proteome changes to understand the regulatory mechanism of carbon metabolism in Synechocystis. Remarkably, four proteins were found to be commonly regulated by all the five RCM proteins (Fig. 3B). Of these, transcription of ggpS was reported be upregulated in response to osmotic stress and salt stress (73, 74, 75), and the transcription of slr1161 was reported to be repressed by the stresses of high light, 3-(3,4-dichlorophenyl)-1,1-dimethyl urea (DCMU), and 2,5-dibromo-3-methyl-6-isopropyl-p-benzoquinone (DBMIB) (76, 77, 78). The unanimous upregulation and downregulation of GgpS and Slr1161, respectively, in all five RCM mutants suggest that the mutants need to cope with the stresses generated more or less by the depletion of the RCM proteins. Sll7085 is a member of the CRISPR-associated proteins superfamily (48, 79), which is downregulated in all mutants but to greater extents in ΔsigE and Δhik8. The coding gene of Sll7085 is part of an operon on the endogenous plasmid pSYSA (80), and the operon also encodes the other components of the same CRISPR complex. It remains elusive why Sll7085 but not the other components of the CRISPR complex were dramatically downregulated in the two mutants. Slr6074 is not previously associated with any stress responses, indicating it could be a bona fide carbon metabolism-related target of RCM proteins. Indeed, the knockout mutant of slr6074 displayed severe defects in trophic growth, and the results for this part of the work will be reported separately. The dark viability of ΔsigE, Δhik8, and Δrre37 in the presence of glucose was all impaired as also shown by previous reports (supplemental Fig. S2D) (15, 62, 81). Impaired sugar catabolism has been proposed to account for this phenotype because of the inability to provide the energy necessary for dark survival through catabolizing glucose (8, 27). Moreover, the accumulation of glucose to a certain level could also be toxic for cyanobacteria (13, 46, 82, 83). In Synechocystis, and many other cyanobacteria strains as well, OPPP is the major route for catabolizing glucose in the dark (7). This, together with the downstream part of glycolysis that starts from the reaction catalyzed by Gap1, channels the vast majority of carbon flux for the production of energy, reducing equivalent, and building blocks for biomass (84). In Synechocystis, Nostoc sp. ATCC 29133, and Synechococcus, mutants depleted of Zwf displayed impaired dark viability in the presence of glucose (10, 58, 85). In Δhik8 and ΔsigE, four proteins in OPPP including Zwf, OpcA, Gnd, and Tal were unanimously downregulated, indicative of defective OPPP. Thus, the impaired dark viability of the two mutants could be mainly attributable to the downregulation of these proteins. In addition, three proteins in glycolysis including Pyk1, PfkA1, and Gap1 were downregulated in both mutants, and this may also contribute more or less to the impaired dark growth (84). In Δrre37, however, only two of the aforementioned proteins in sugar catabolism were downregulated (PfkA1 and Gap1), and the extent of downregulation is not as great as that in ΔsigE and Δhik8. This result may correlate with the relatively less serve repression of the growth of glucose in the dark as observed (supplemental Fig. S2D). Nevertheless, the mechanism of impaired dark viability in the presence of glucose is complex, alteration of other proteins in carbon metabolism may also be attributable to this phenotype. Dramatically decreased glycogen content and massive upregulation of PII S49 phosphorylation in Δhik8 but not in the other RCM mutants are the most prominent findings in the present study (Fig. 6). Decreased glycogen content in an AT-grown hik8-deficient mutant was previously reported (27), but it is a bit surprising that this occurs only in Δhik8 but not in other mutants such as ΔsigE, considering that expression patterns of many proteins in carbon metabolism are very similar between in Δhik8 and ΔsigE. The observation underscores the specificity of Hik8 in regulating carbon metabolism, particularly, glycogen accumulation. The massive increase of PII phosphorylation in Δhik8 but not in the other mutants perfectly correlates, though negatively, with the decreased glycogen content (Fig. 6), and this prompted us to investigate their causal relationship. The successful restoration of glycogen content to a level comparable with that in the WT by PII S49A mutation suggests that hyperphosphorylation of PII is more probably attributable to the reduced glycogen content in Δhik8. It is highly remarkable that the glycogen contents in different Synechocystis strains in the current study are at a basal (∼1 μg/OD730) or sub-basal level (Fig. 6A), which is up to two orders of magnitude lower than that induced by nitrogen starvation (86, 87). It has been well established that nitrogen starvation-induced PII phosphorylation positively correlates with massive accumulation of glycogen in WT Synechocystis (88, 89), but this does not necessarily contradict our result that hyperphosphorylation of PII negatively correlates with basal level glycogen content in Δhik8. It is expectable that basal and stimulated levels of glycogen accumulation might be differentially regulated at both molecular and metabolic levels. Recently, Orthwein et al. reported under nitrogen starvation a small protein PirC (Sll0944) specifically binds to and inhibits 2,3-phosphoglycerate-independent phosphoglycerate mutase (PGAM), an enzyme important for the lower part of the glycolysis pathway. Inhibition of PGAM would re-direct the carbon flux from the photosynthetic Calvin cycle to glycogen synthesis, and result in massive accumulation of glycogen in nitrogen-starved cells. Under nitrogen repletion, the binding of PirC by PII sequesters PirC from and activates PGAM, which directs carbon flux to the downstream TCA cycle and represses glycogen synthesis (88, 89). In our case, PII S49A mutation restored glycogen from a sub-basal level to the basal level in Δhik8 but did not significantly change the basal level glycogen content in the WT (Fig. 7A). If PII S49A mutation inhibits PII-PirC interaction and promoters PirC binding to and inhibiting PGAM, and consequently increases glycogen content in Δhik8, then it is expectable that S49A mutation would also significantly increase the glycogen content in the WT, which is obviously not the case (Fig. 7A). Thus, the mechanism of PirC-regulated glycogen synthesis could be more appropriate to explain the induced massive accumulation of glycogen that requires a dramatic directional change in metabolic carbon flux. However, whether the mechanism is equally applicable to explain the regulation of glycogen accumulation at a basal or sub-basal level is still questionable. Together, although PII S49 phosphorylation was regarded as a well-known signal for the intracellular C/N balance (90, 91, 92), its regulation on glycogen accumulation, particularly at the basal level under the control of Hik8, is just revealed. To our knowledge, this is the first line of evidence that links PII phosphorylation to a histidine kinase of the two-component system in cyanobacteria. Notwithstanding the effective restoration of the glycogen content and the prevention of cell death in Δhik8 by PII S49A mutation (Fig. 7), the dark growth of glnBS49A/WTe/Δhik8 with glucose was not restored to a level commensurate with that of the glnBS49A/WTe, and the later was also partly repressed in growth compared with the WT (Fig. 7B). The observation suggests that PII S49A substitution favors dark survival of Δhik8 in the presence of glucose. The increased glycogen resulting from S49A substitution in glnBS49A/WTe/Δhik8 and its catabolism might provide energy, through respiratory electron transport, necessary for Δhik8 to survive the prolonged darkness (27). Nevertheless, the PII S49A substitution is also harmful to cell growth as displayed by the growth phenotype of glnBS49A/WTe (Fig. 7B). This is not surprising because PII S49A is unmodifiable and hence loses the flexibility to regulate or fine-tune C/N metabolism through reversible phosphorylation, which is harmful for the cyanobacterium to cope with the adverse environment, such as prolonged darkness. Nevertheless, the network of carbon metabolism is highly complex and interconnected with nitrogen metabolism (64, 65). Multiple pathways including a number of reactions (uni-direction or bi-direction) are confined within the same cellular compartment. Coordination of these pathways and reactions requires tight and precise spatiotemporal regulation that occurs at multiple levels on mRNAs, proteins, and metabolites, and different parts of a pathway may also be differentially regulated to precisely divert the carbon flux to different destinations in response to environmental changes. Therefore, a detailed mechanism of Hik8-regulated glycogen accumulation, PII phosphorylation, and dark survival with glucose is far from elucidated. Reversible PII S49 phosphorylation is a well-established mechanism for sensing the state of C/N homeostasis in many cyanobacteria strains and has been extensively discussed (67). Nevertheless, the functions of phosphorylated and unphosphorylated forms of PII are not well studied. PII is necessary for ammonium-responsive inhibition of nitrate assimilation (93). It was reported that neither S49A nor S49D substitution affects nitrate assimilation and ammonium-responsive inhibition of nitrate assimilation in Synechocystis and Synechococcus (94, 95), and S49D mutation dramatically reduced PII-interaction with Amt1, NrtD, and UrtE, which are involved in transporting ammonium, nitrate, and urea, respectively (96). Other than this scarce information is available regarding the phenotype of cyanobacteria strains carrying PII S49A substitution. Our results strongly suggest that PII phosphorylation is involved in the regulation of the tropic growth of Synechocystis. Specifically, the capability of PII being phosphorylated is important for dark growth (Fig. 7B and supplemental Fig. S15), but not as important for AT and MT growth (supplemental Fig. S16). Hik8 is a well-known two-component sensory histidine kinase involved in the regulation of circadian rhythmicity. The Hik8 ortholog in Synechococcus, SasA, interacts directly with the circadian clock protein KaiC (20), and the latter has both autokinase and autophosphatase activities towards phosphorylation and dephosphorylation of S431 and T432 on KaiC, respectively. The dual phosphorylation of KaiC was reported to be the basis of circadian rhythm in cyanobacteria (72). Prominently, the Synechocystis ortholog of KaiC, namely KaiC1, is the only circadian clock protein upregulated in Δhik8 (Fig. 2C and supplemental Fig. S17). KaiC1 also interacts directly with Hik8 (20, 26, 55). KaiC1 could regulate the circadian rhythm of the phosphotransfer activity from Hik8/SasA to its cognate response regulator RpaA/Rre31 (20, 97, 98, 99), and subsequently regulate rhythmic oscillation of gene expression and metabolic activities (21). RpaA depletion induces a major metabolic shift of the mutant that largely phenocopies Δhik8 (100). It is conceivable that the inactivation of Hik8 changed the phosphorylation status of RpaA and its activity as a transcription factor and affected the expression of one or more genes that regulate the phosphorylation of PII. The observed upregulation of KaiC1 in Δhik8 could be a feedback response of Hik8 inactivation (supplemental Fig. S17). Indeed, PII phosphorylation is also modulated in response to light and darkness (32), indicating that PII phosphorylation could also be under circadian regulation. Unfortunately, to date the identity of the kinases responsible for PII phosphorylation remains elusive, and it has been suggested that multiple protein kinases that complement with each other are involved (101). Overexpression RpaA in Synechocystis was reported to slow down glycogen degradation in dark (102). It is conceivable that the inactivation of RpaA or Hik8 may accelerate glycogen degradation, resulting in reduced glycogen content. In Synechococcus, the glycogen content displayed a typical circadian oscillation peaking at the subjective dusks and reaching troughs at the subjective dawns (103, 104), such a circadian rhythm of glycogen content is more likely regulated by the output signals from the circadian clock mediated by Hik8-RpaA. The abolished circadian rhythm in Δhik8 is expected to disrupt the circadian rhythm of the glycogen accumulation and arrest the mutant in a state with constitutively low glycogen content as also suggested by Kawasaki et al (104), and this process is probably mediated through hyperphosphorylation of PII. Together, our findings could help to make connections among Hik8, KaiC1, and PII, and to place PII-controlled C/N homeostasis under the regulation of the circadian clock with Hik8-RpaA as the major output pathway. Data Availability The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD035864 and PXD035898 (105). Supplemental data This article contains supplemental data (106). Conflict of interest The authors declare no competing interests. Supplementary Data Supplemental Figures S1–S18 Supplemental Table S1 Supplemental Table S2 Supplemental Table S3 Supplemental Table S4 Supplemental Table S5 Supplemental Table S6 Supplemental Table S7 Supplemental Table S8 Acknowledgments We thank Dr Jinlong Wang, Dr Weiyang Chen and Dr Longfa Fang for the suggestions in data analysis. Funding and additional information The work was supported by a grant 2019YFA0802203 from the Ministry of Science and Technology of China, a grant from Strategic Priority Research Program (XDA24040202) of 10.13039/501100002367 Chinese Academy of Sciences , and a grant from 10.13039/501100001809 National Natural Science Foundation of China (32170253). Author contributions Yan Wang supervision and methodology; C. H., H. G., X. D., Z. X. Yan Wang., J. D., G. W., and L. Z. investigation; C. H., H. G., and X. D. formal analysis; X. H. and Y. Z. project administration; C. H. and Yingchun Wang writing – original draft; W. X. and Yan Wang writing – review and editing. ==== Refs References 1 Stanier R.Y. Cohen-Bazire G. Phototrophic prokaryotes: the cyanobacteria Annu. Rev. Microbiol. 31 1977 225 274 410354 2 Oliver N.J. Rabinovitch-Deere C.A. Carroll A.L. Nozzi N.E. Case A.E. Atsumi S. Cyanobacterial metabolic engineering for biofuel and chemical production Curr. Opin. Chem. Biol. 35 2016 43 50 27614173 3 Berla B.M. Saha R. Immethun C.M. Maranas C.D. Moon T.S. Pakrasi H.B. Synthetic biology of cyanobacteria: unique challenges and opportunities Front. Microbiol. 4 2013 246 24009604 4 Rae B.D. Long B.M. Badger M.R. Price G.D. Functions, compositions, and evolution of the two types of carboxysomes: polyhedral microcompartments that facilitate CO2 fixation in cyanobacteria and some proteobacteria Microbiol. Mol. Biol. Rev. 77 2013 357 379 24006469 5 Ball S.G. Morell M.K. From bacterial glycogen to starch: understanding the biogenesis of the plant starch granule Annu. Rev. Plant Biol. 54 2003 207 233 14502990 6 Grundel M. Scheunemann R. Lockau W. Zilliges Y. Impaired glycogen synthesis causes metabolic overflow reactions and affects stress responses in the cyanobacterium Synechocystis sp. PCC 6803 Microbiology 158 2012 3032 3043 23038809 7 Yang C. Hua Q. Shimizu K. Metabolic flux analysis in Synechocystis using isotope distribution from 13C-labeled glucose Metab. Eng. 4 2002 202 216 12616690 8 Osanai T. Azuma M. Tanaka K. Sugar catabolism regulated by light- and nitrogen-status in the cyanobacterium Synechocystis sp. PCC 6803 Photochem. Photobiol. Sci. 6 2007 508 514 17487300 9 Kurian D. Jansen T. Maenpaa P. Proteomic analysis of heterotrophy in Synechocystis sp. PCC 6803 Proteomics 6 2006 1483 1494 16421936 10 Summers M.L. Wallis J.G. Campbell E.L. Meeks J.C. Genetic evidence of a major role for glucose-6-phosphate dehydrogenase in nitrogen fixation and dark growth of the cyanobacterium Nostoc sp. strain ATCC 29133 J. Bacteriol. 177 1995 6184 6194 7592384 11 Azuma M. Osanai T. Hirai M.Y. Tanaka K. A response regulator Rre37 and an RNA polymerase sigma factor SigE represent two parallel pathways to activate sugar catabolism in a cyanobacterium Synechocystis sp. PCC 6803 Plant Cell Physiol. 52 2011 404 412 21216748 12 Dorrich A.K. Mitschke J. Siadat O. Wilde A. Deletion of the Synechocystis sp. PCC 6803 kaiAB1C1 gene cluster causes impaired cell growth under light-dark conditions Microbiology 160 2014 2538 2550 25139948 13 Gao L. Shen C. Liao L. Huang X. Liu K. Wang W. Functional proteomic discovery of Slr0110 as a central regulator of carbohydrate metabolism in Synechocystis species PCC 6803 Mol. Cell. Proteomics 13 2014 204 219 24169622 14 Iijima H. Watanabe A. Takanobu J. Hirai M.Y. Osanai T. rre37 overexpression alters gene expression related to the tricarboxylic acid cycle and pyruvate metabolism in Synechocystis sp. PCC 6803 ScientificWorldJournal 2014 2014 921976 15 Osanai T. Kanesaki Y. Nakano T. Takahashi H. Asayama M. Shirai M. Positive regulation of sugar catabolic pathways in the cyanobacterium Synechocystis sp. PCC 6803 by the group 2 sigma factor sigE J. Biol. Chem. 280 2005 30653 30659 15944148 16 Osanai T. Oikawa A. Azuma M. Tanaka K. Saito K. Hirai M.Y. Genetic engineering of group 2 sigma factor SigE widely activates expressions of sugar catabolic genes in Synechocystis species PCC 6803 J. Biol. Chem. 286 2011 30962 30971 21757761 17 Tabei Y. Okada K. Tsuzuki M. Sll1330 controls the expression of glycolytic genes in Synechocystis sp. PCC 6803 Biochem. Biophys. Res. Commun. 355 2007 1045 1050 17331473 18 Nagarajan S. Sherman D.M. Shaw I. Sherman L.A. Functions of the duplicated hik31 operons in central metabolism and responses to light, dark, and carbon sources in Synechocystis sp. strain PCC 6803 J. Bacteriol. 194 2012 448 459 22081400 19 Nagarajan S. Srivastava S. Sherman L.A. Essential role of the plasmid hik31 operon in regulating central metabolism in the dark in Synechocystis sp. PCC 6803 Mol. Microbiol. 91 2014 79 97 24237382 20 Iwasaki H. Williams S.B. Kitayama Y. Ishiura M. Golden S.S. Kondo T. A kaiC-interacting sensory histidine kinase, SasA, necessary to sustain robust circadian oscillation in cyanobacteria Cell 101 2000 223 233 10786837 21 Markson J.S. Piechura J.R. Puszynska A.M. O'Shea E.K. Circadian control of global gene expression by the cyanobacterial master regulator RpaA Cell 155 2013 1396 1408 24315105 22 Osanai T. Kuwahara A. Iijima H. Toyooka K. Sato M. Tanaka K. Pleiotropic effect of sigE over-expression on cell morphology, photosynthesis and hydrogen production in Synechocystis sp. PCC 6803 Plant J. 76 2013 456 465 23941239 23 Tokumaru Y. Uebayashi K. Toyoshima M. Osanai T. Matsuda F. Shimizu H. Comparative targeted proteomics of the central metabolism and photosystems in SigE mutant strains of Synechocystis sp. PCC 6803 Molecules 23 2018 1051 29723969 24 Osanai T. Imamura S. Asayama M. Shirai M. Suzuki I. Murata N. Nitrogen induction of sugar catabolic gene expression in Synechocystis sp. PCC 6803 DNA Res. 13 2006 185 195 17046957 25 Osanai T. Numata K. Oikawa A. Kuwahara A. Iijima H. Doi Y. Increased bioplastic production with an RNA polymerase sigma factor SigE during nitrogen starvation in Synechocystis sp. PCC 6803 DNA Res. 20 2013 525 535 23861321 26 Osanai T. Shirai T. Iijima H. Kuwahara A. Suzuki I. Kondo A. Alteration of cyanobacterial sugar and amino acid metabolism by overexpression hik8, encoding a KaiC-associated histidine kinase Environ. Microbiol. 17 2015 2430 2440 25403325 27 Singh A.K. Sherman L.A. Pleiotropic effect of a histidine kinase on carbohydrate metabolism in Synechocystis sp. strain PCC 6803 and its requirement for heterotrophic growth J. Bacteriol. 187 2005 2368 2376 15774880 28 Nakamura Y. Kaneko T. Tabata S. CyanoBase, the genome database for Synechocystis sp. strain PCC6803: status for the year 2000 Nucleic Acids Res. 28 2000 72 10592184 29 Vermaas W. Molecular genetics of the cyanobacterium Synechocystis sp. PCC 6803: principles and possible biotechnology applications J. Appl. Phycol. 8 1996 263 273 30 Kaneko T. Sato S. Kotani H. Tanaka A. Asamizu E. Nakamura Y. Sequence analysis of the genome of the unicellular cyanobacterium Synechocystis sp. strain PCC 6803. II. Sequence determination of the entire genome and assignment of potential protein-coding regions DNA Res. 3 1996 109 136 8905231 31 Toennies G. Gallant D.L. The relation between photometric turbidity and bacterial concentration Growth 13 1949 7 20 18150958 32 Fang L. Ge H. Huang X. Liu Y. Lu M. Wang J. Trophic mode-dependent proteomic analysis reveals functional significance of light-independent chlorophyll synthesis in Synechocystis sp. PCC 6803 Mol. Plant 10 2017 73 85 27585879 33 Bustos S.A. Golden S.S. Light-regulated expression of the psbD gene family in Synechococcus sp. strain PCC 7942: evidence for the role of duplicated psbD genes in cyanobacteria Mol. Gen. Genet. 232 1992 221 230 1372952 34 Hauf W. Schmid K. Gerhardt E.C. Huergo L.F. Forchhammer K. Interaction of the nitrogen regulatory protein GlnB (PII) with biotin carboxyl carrier protein (BCCP) controls Acetyl-CoA levels in the Cyanobacterium Synechocystis sp. PCC 6803 Front. Microbiol. 7 2016 1700 27833596 35 Moran R. Formulae for determination of chlorophyllous pigments extracted with n,n-dimethylformamide Plant Physiol. 69 1982 1376 1381 16662407 36 Xu Y. Guerra L.T. Li Z. Ludwig M. Dismukes G.C. Bryant D.A. Altered carbohydrate metabolism in glycogen synthase mutants of Synechococcus sp. strain PCC 7002: cell factories for soluble sugars Metab. Eng. 16 2013 56 67 23262095 37 Chen W. Zheng L. Dong J. Ge H. Huang X. Wang G. A systematic survey of the light/dark-dependent protein degradation events in a model cyanobacterium Mol. Cell. Proteomics 20 2021 100162 38 Udeshi N.D. Svinkina T. Mertins P. Kuhn E. Mani D.R. Qiao J.W. Refined preparation and use of anti-diglycine remnant (K-epsilon-GG) antibody enables routine quantification of 10,000s of ubiquitination sites in single proteomics experiments Mol. Cell. Proteomics 12 2013 825 831 23266961 39 Rappsilber J. Ishihama Y. Mann M. Stop and go extraction tips for matrix-assisted laser desorption/ionization, nanoelectrospray, and LC/MS sample pretreatment in proteomics Anal. Chem. 75 2003 663 670 12585499 40 Peterson A.C. Russell J.D. Bailey D.J. Westphall M.S. Coon J.J. Parallel reaction monitoring for high resolution and high mass accuracy quantitative, targeted proteomics Mol. Cell. Proteomics 11 2012 1475 1488 22865924 41 Cox J. Mann M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification Nat. Biotechnol. 26 2008 1367 1372 19029910 42 Cox J. Mann M. 1D and 2D annotation enrichment: a statistical method integrating quantitative proteomics with complementary high-throughput data BMC Bioinformatics 13 2012 S12 43 MacLean B. Tomazela D.M. Shulman N. Chambers M. Finney G.L. Frewen B. Skyline: an open source document editor for creating and analyzing targeted proteomics experiments Bioinformatics 26 2010 966 968 20147306 44 Tyanova S. Temu T. Cox J. The MaxQuant computational platform for mass spectrometry-based shotgun proteomics Nat. Protoc. 11 2016 2301 2319 27809316 45 Osanai T. Sato S. Tabata S. Tanaka K. Identification of PamA as a PII-binding membrane protein important in nitrogen-related and sugar-catabolic gene expression in Synechocystis sp. PCC 6803 J. Biol. Chem. 280 2005 34684 34690 16109709 46 Kahlon S. Beeri K. Ohkawa H. Hihara Y. Murik O. Suzuki I. A putative sensor kinase, Hik31, is involved in the response of Synechocystis sp. strain PCC 6803 to the presence of glucose Microbiology 152 2006 647 655 16514145 47 Dayon L. Hainard A. Licker V. Turck N. Kuhn K. Hochstrasser D.F. Relative quantification of proteins in human cerebrospinal fluids by MS/MS using 6-plex isobaric tags Anal. Chem. 80 2008 2921 2931 18312001 48 Ge H. Fang L. Huang X. Wang J. Chen W. Liu Y. Translating divergent environmental stresses into a common proteome response through the histidine kinase 33 (Hik33) in a model cyanobacterium Mol. Cell. Proteomics 16 2017 1258 1274 28668777 49 Ting L. Rad R. Gygi S.P. Haas W. MS3 eliminates ratio distortion in isobaric multiplexed quantitative proteomics Nat. Methods 8 2011 937 940 21963607 50 Wang J. Zhang Y. Huang X. Lu D. Wang Y. Evaluation of the potential risk of advanced peak determination in distorting isobaric labeling-based single-shot proteome quantitation Proteomics 20 2020 e1900255 51 Muro-Pastor M.I. Reyes J.C. Florencio F.J. Ammonium assimilation in cyanobacteria Photosynth. Res. 83 2005 135 150 16143848 52 Shibata M. Katoh H. Sonoda M. Ohkawa H. Shimoyama M. Fukuzawa H. Genes essential to sodium-dependent bicarbonate transport in cyanobacteria: function and phylogenetic analysis J. Biol. Chem. 277 2002 18658 18664 11904298 53 Sauer U. Canonaco F. Heri S. Perrenoud A. Fischer E. The soluble and membrane-bound transhydrogenases UdhA and PntAB have divergent functions in NADPH metabolism of Escherichia coli J. Biol. Chem. 279 2004 6613 6619 14660605 54 Hagen K.D. Meeks J.C. The unique cyanobacterial protein OpcA is an allosteric effector of glucose-6-phosphate dehydrogenase in Nostoc punctiforme ATCC 29133 J. Biol. Chem. 276 2001 11477 11486 11152472 55 Kobler C. Schultz S.J. Kopp D. Voigt K. Wilde A. The role of the Synechocystis sp. PCC 6803 homolog of the circadian clock output regulator RpaA in day-night transitions Mol. Microbiol. 110 2018 847 861 30216574 56 Griffin T.J. Gygi S.P. Ideker T. Rist B. Eng J. Hood L. Complementary profiling of gene expression at the transcriptome and proteome levels in Saccharomyces cerevisiae Mol. Cell. Proteomics 1 2002 323 333 12096114 57 Okada K. Horii E. Nagashima Y. Mitsui M. Matsuura H. Fujiwara S. Genes for a series of proteins that are involved in glucose catabolism are upregulated by the Hik8-cascade in Synechocystis sp. PCC 6803 Planta 241 2015 1453 1462 25732003 58 Scanlan D.J. Sundaram S. Newman J. Mann N.H. Carr N.G. Characterization of a zwf mutant of Synechococcus sp. strain PCC 7942 J. Bacteriol. 177 1995 2550 2553 7730289 59 Jansén T. Kurian D. Raksajit W. York S. Mäenpää S. Characterization of trophic changes and a functional oxidative pentose phosphate pathway in Synechocystis sp. PCC 6803 Acta Physiol. Plant. 32 2010 511 518 60 Klotz A. Forchhammer K. Glycogen, a major player for bacterial survival and awakening from dormancy Future Microbiol. 12 2017 101 104 28106464 61 Kuwahara A. Arisaka S. Takeya M. Iijima H. Hirai M.Y. Osanai T. Modification of photosynthetic electron transport and amino acid levels by overexpression of a circadian-related histidine kinase hik8 in Synechocystis sp. PCC 6803 Front. Microbiol. 6 2015 1150 26539179 62 Tabei Y. Okada K. Horii E. Mitsui M. Nagashima Y. Sakai T. Two regulatory networks mediated by light and glucose involved in glycolytic gene expression in cyanobacteria Plant Cell Physiol. 53 2012 1720 1727 22915573 63 Krasikov V. Aguirre von Wobeser E. Dekker H.L. Huisman J. Matthijs H.C. Time-series resolution of gradual nitrogen starvation and its impact on photosynthesis in the cyanobacterium Synechocystis PCC 6803 Physiol. Plant. 145 2012 426 439 22289076 64 Schwarz R. Forchhammer K. Acclimation of unicellular cyanobacteria to macronutrient deficiency: emergence of a complex network of cellular responses Microbiology 151 2005 2503 2514 16079330 65 Klotz A. Georg J. Bucinska L. Watanabe S. Reimann V. Januszewski W. Awakening of a dormant cyanobacterium from nitrogen chlorosis reveals a genetically determined program Curr. Biol. 26 2016 2862 2872 27720620 66 Forchhammer K. Schwarz R. Nitrogen chlorosis in unicellular cyanobacteria - a developmental program for surviving nitrogen deprivation Environ. Microbiol. 21 2019 1173 1184 30318768 67 Forchhammer K. Selim K.A. Carbon/nitrogen homeostasis control in cyanobacteria FEMS Microbiol. Rev. 44 2020 33 53 31617886 68 Yang M.K. Qiao Z.X. Zhang W.Y. Xiong Q. Zhang J. Li T. Global phosphoproteomic analysis reveals diverse functions of serine/threonine/tyrosine phosphorylation in the model cyanobacterium Synechococcus sp. strain PCC 7002 J. Proteome Res. 12 2013 1909 1923 23461524 69 Wang Y. Ding S.J. Wang W. Jacobs J.M. Qian W.J. Moore R.J. Profiling signaling polarity in chemotactic cells Proc. Natl. Acad. Sci. U. S. A. 104 2007 8328 8333 17494752 70 Elschenbroich S. Kislinger T. Targeted proteomics by selected reaction monitoring mass spectrometry: applications to systems biology and biomarker discovery Mol. Biosyst. 7 2011 292 303 20976349 71 Pan S. Aebersold R. Chen R. Rush J. Goodlett D.R. McIntosh M.W. Mass spectrometry based targeted protein quantification: methods and applications J. Proteome Res. 8 2009 787 797 19105742 72 Nishiwaki T. Satomi Y. Kitayama Y. Terauchi K. Kiyohara R. Takao T. A sequential program of dual phosphorylation of KaiC as a basis for circadian rhythm in cyanobacteria EMBO J. 26 2007 4029 4037 17717528 73 Rübsam H. Kirsch F. Reimann V. Erban A. Kopka J. Hagemann M. The iron-stress activated RNA 1 (IsaR1) coordinates osmotic acclimation and iron starvation responses in the cyanobacterium Synechocystis sp. PCC 6803 Environ. Microbiol. 20 2018 2757 2768 29468839 74 Kirsch F. Pade N. Klähn S. Hess W.R. Hagemann M. The glucosylglycerol-degrading enzyme GghA is involved in acclimation to fluctuating salinities by the cyanobacterium Synechocystis sp. strain PCC 6803 Microbiology 163 2017 1319 1328 28857041 75 Chen L. Wu L. Zhu Y. Song Z. Wang J. Zhang W. An orphan two-component response regulator Slr1588 involves salt tolerance by directly regulating synthesis of compatible solutes in photosynthetic Synechocystis sp. PCC 6803 Mol. Biosyst. 10 2014 1765 1774 24718497 76 Kobayashi M. Ishizuka T. Katayama M. Kanehisa M. Bhattacharyya-Pakrasi M. Pakrasi H.B. Response to oxidative stress involves a novel peroxiredoxin gene in the unicellular cyanobacterium Synechocystis sp. PCC 6803 Plant Cell Physiol. 45 2004 290 299 15047877 77 Hihara Y. Sonoike K. Kanehisa M. Ikeuchi M. DNA microarray analysis of redox-responsive genes in the genome of the cyanobacterium Synechocystis sp. strain PCC 6803 J. Bacteriol. 185 2003 1719 1725 12591891 78 Hihara Y. Kamei A. Kanehisa M. Kaplan A. Ikeuchi M. DNA microarray analysis of cyanobacterial gene expression during acclimation to high light Plant Cell 13 2001 793 806 11283337 79 Scholz I. Lange S.J. Hein S. Hess W.R. Backofen R. CRISPR-Cas systems in the cyanobacterium Synechocystis sp. PCC6803 exhibit distinct processing pathways involving at least two Cas6 and a Cmr2 protein PLoS One 8 2013 e56470 80 Kaneko T. Nakamura Y. Sasamoto S. Watanabe A. Kohara M. Matsumoto M. Structural analysis of four large plasmids harboring in a unicellular cyanobacterium, Synechocystis sp. PCC 6803 DNA Res. 10 2003 221 228 14686584 81 Osanai T. Oikawa A. Numata K. Kuwahara A. Iijima H. Doi Y. Pathway-level acceleration of glycogen catabolism by a response regulator in the cyanobacterium Synechocystis species PCC 6803 Plant Physiol. 164 2014 1831 1841 24521880 82 Hihara Y. Ikeuchi M.J.P.R. Mutation in a novel gene required for photomixotrophic growth leads to enhanced photoautotrophic growth of Synechocystis sp. PCC 6803 Photosyn. Res. 53 1997 243 252 83 Beuf L. Bédu S. Durand M.C. Joset F. A protein involved in co-ordinated regulation of inorganic carbon and glucose metabolism in the facultative photoautotrophic cyanobacterium Synechocystis PCC 6803 Plant Mol. Biol. 25 1994 855 864 8075401 84 Koksharova O. Schubert M. Shestakov S. Cerff R. Genetic and biochemical evidence for distinct key functions of two highly divergent GAPDH genes in catabolic and anabolic carbon flow of the cyanobacterium Synechocystis sp. PCC 6803 Plant Mol. Biol. 36 1998 183 194 9484473 85 Janse´n T. Kurian D. Raksajit W. York S. Summers M.L. Ma¨enpa¨a P. Characterization of trophic changes and a functional oxidative pentose phosphate pathway in Synechocystis sp. PCC 6803 Acta Physiol. Plant. 32 2010 511 518 86 Koch M. Doello S. Gutekunst K. Forchhammer K. PHB is produced from glycogen turn-over during nitrogen starvation in Synechocystis sp. PCC 6803 Int. J. Mol. Sci. 20 2019 1942 31010017 87 Muro-Pastor M.I. Cutillas-Farray A. Perez-Rodriguez L. Perez-Saavedra J. Vega-de Armas A. Paredes A. CfrA, a novel carbon flow regulator, adapts carbon metabolism to nitrogen deficiency in cyanobacteria Plant Physiol. 184 2020 1792 1810 32900980 88 Orthwein T. Scholl J. Spat P. Lucius S. Koch M. Macek B. The novel PII-interactor PirC identifies phosphoglycerate mutase as key control point of carbon storage metabolism in cyanobacteria Proc. Natl. Acad. Sci. U. S. A. 118 2021 e2019988118 89 Forchhammer K. Selim K.A. Huergo L.F. New views on PII signaling: from nitrogen sensing to global metabolic control Trends Microbiol. 30 2022 722 735 35067429 90 Kloft N. Rasch G. Forchhammer K. Protein phosphatase PphA from Synechocystis sp. PCC 6803: the physiological framework of PII-P dephosphorylation Microbiology 151 2005 1275 1283 15817794 91 Forchhammer K. Global carbon/nitrogen control by PII signal transduction in cyanobacteria: from signals to targets FEMS Microbiol. Rev. 28 2004 319 333 15449606 92 Forchhammer K. Tandeau de Marsac N. The PII protein in the cyanobacterium Synechococcus sp. strain PCC 7942 is modified by serine phosphorylation and signals the cellular N-status J. Bacteriol. 176 1994 84 91 8282715 93 Lee H.M. Flores E. Forchhammer K. Herrero A. Tandeau De Marsac N. Phosphorylation of the signal transducer PII protein and an additional effector are required for the PII-mediated regulation of nitrate and nitrite uptake in the cyanobacterium Synechococcus sp. PCC 7942 Eur. J. Biochem. 267 2000 591 600 10632730 94 Kobayashi M. Takatani N. Tanigawa M. Omata T. Posttranslational regulation of nitrate assimilation in the cyanobacterium Synechocystis sp. strain PCC 6803 J. Bacteriol. 187 2005 498 506 15629921 95 Takatani N. Kobayashi M. Maeda S. Omata T. Regulation of nitrate reductase by non-modifiable derivatives of PII in the cells of Synechococcus elongatus strain PCC 7942 Plant Cell Physiol. 47 2006 1182 1186 16854940 96 Watzer B. Spat P. Neumann N. Koch M. Sobotka R. Macek B. The signal transduction protein PII controls ammonium, nitrate and urea uptake in cyanobacteria Front. Microbiol. 10 2019 1428 31293555 97 Takai N. Nakajima M. Oyama T. Kito R. Sugita C. Sugita M. A KaiC-associating SasA-RpaA two-component regulatory system as a major circadian timing mediator in cyanobacteria Proc. Natl. Acad. Sci. U. S. A. 103 2006 12109 12114 16882723 98 Kageyama H. Nishiwaki T. Nakajima M. Iwasaki H. Oyama T. Kondo T. Cyanobacterial circadian pacemaker: kai protein complex dynamics in the KaiC phosphorylation cycle in vitro Mol. Cell 23 2006 161 171 16857583 99 Gutu A. O'Shea E.K. Two antagonistic clock-regulated histidine kinases time the activation of circadian gene expression Mol. Cell 50 2013 288 294 23541768 100 Scheurer N.M. Rajarathinam Y. Timm S. Kobler C. Kopka J. Hagemann M. Homologs of circadian clock proteins impact the metabolic switch between light and dark growth in the cyanobacterium Synechocystis sp. PCC 6803 Front. Plant Sci. 12 2021 675227 101 Spat P. Macek B. Forchhammer K. Phosphoproteome of the cyanobacterium Synechocystis sp. PCC 6803 and its dynamics during nitrogen starvation Front. Microbiol. 6 2015 248 25873915 102 Iijima H. Shirai T. Okamoto M. Kondo A. Hirai M.Y. Osanai T. Changes in primary metabolism under light and dark conditions in response to overproduction of a response regulator RpaA in the unicellular cyanobacterium Synechocystis sp. PCC 6803 Front. Microbiol. 6 2015 888 26379657 103 Pattanayak G.K. Phong C. Rust M.J. Rhythms in energy storage control the ability of the cyanobacterial circadian clock to reset Curr. Biol. 24 2014 1934 1938 25127221 104 Kawasaki K. Iwasaki H. Involvement of glycogen metabolism in circadian control of UV resistance in cyanobacteria Plos Genet. 16 2020 e1009230 105 Perez-Riverol Y. Bai J. Bandla C. Garcia-Seisdedos D. Hewapathirana S. Kamatchinathan S. The PRIDE database resources in 2022: a hub for mass spectrometry-based proteomics evidences Nucleic Acids Res. 50 2022 D543 D552 34723319 106 Olsen J.V. Vermeulen M. Santamaria A. Kumar C. Miller M.L. Jensen L.J. Quantitative phosphoproteomics reveals widespread full phosphorylation site occupancy during mitosis Sci. Signal. 3 2010 ra3