Intracellular Ser/Thr/Tyr phosphoproteome of the oral commensal Streptococcus gordonii DL1

Background To respond and adapt to environmental challenges, prokaryotes regulate cellular processes rapidly and reversibly through protein phosphorylation and dephosphorylation. This study investigates the intracellular proteome and Ser/Thr/Tyr phosphoproteome of the oral commensal Streptococcus gordonii. Intracellular proteins from planktonic cells of S. gordonii DL1 were extracted and subjected to 2D-gel electrophoresis. Proteins in general were visualized using Coomassie Brilliant Blue and T-Rex staining. Phosphorylated proteins were visualized with Pro-Q Diamond Phosphoprotein Gel Stain. Proteins were identified by LC-MS/MS and sequence analysis. Results In total, sixty-one intracellular proteins were identified in S. gordonii DL1, many of which occurred at multiple isoelectric points. Nineteen of these proteins were present as one or more Ser/Thr/Tyr phosphorylated form. The identified phosphoproteins turned out to be involved in a variety of cellular processes. Conclusion Nineteen phosphoproteins involved in various cellular functions were identified in S. gordonii. This is the first time the global intracellular Ser/Thr/Tyr phosphorylation profile has been analysed in an oral streptococcus. Comparison with phosphoproteomes of other species from previous studies showed many similarities. Proteins that are consistently found in a phosphorylated state across several species and growth conditions may represent a core phosphoproteome profile shared by many bacteria.


Background
Streptococcus gordonii is a pioneer oral colonizer, involved in the establishment of oral biofilms [1]. Due to the numerous different adhesins expressed on its surface, S. gordonii cells readily attach to and colonize dental surfaces [2]. Moreover, S. gordonii can passively migrate from small oral lesions through the blood stream and cause infective endocarditis by opportunistic infection of the heart valves [3]. This mainly saccharolytic species is considered a commensal with ubiquitous habitation in humans, and given the acid production and acid tolerance of S. gordonii, it prevails at the acidic conditions that periodically occur in oral biofilms [4][5][6]. When carbohydrate concentrations are low, S. gordonii like other oral streptococci utilizes a carbohydrate phosphotransferase transport system (PTS) with high glucose affinity [7]. Upon spikes in carbohydrate concentration, oral streptococci are at risk for "sugar killing" from damaging effects caused by accumulated glycolytic intermediates [8]. To evade such inflictions, oral streptococci switch to carbohydrate transport systems with lower affinity for glucose, e.g. the permease system. The glycolytic rate can also be regulated by switching to alternative pathways, e.g. through activation of lactate dehydrogenase (ldh) in pyruvate conversion for faster regeneration of NAD + [9]. In this way, cells reduce glucose uptake as well as drain themselves of glycolytic intermediates by producing large amounts of lactic acid. In addition, oral streptococci employ strategies colloquially termed as the "acid tolerance response" (ATR) that enhance cell survival in acidic environments [4,10,11]. The ATR is dependent on molecular chaperone activity [4,8] that sustains correct protein folding during biosynthesis even at low pH.
To keep up with environmental fluctuations, prokaryotes have been suggested to regulate the activity of proteins involved in the central carbon metabolism rapidly and reversibly through phosphorylation and dephosphorylation mechanisms [12,13]. The reversible regulatory phosphorylation events operate on a much faster time scale than changes in protein expression [14]. The earliest studies on protein phosphorylation and dephosphorylation as regulatory events in bacteria focused on phosphorylation of histidine and aspartate residues in relation to two-component systems [15]. However, phosphorylation events on serine, threonine, and tyrosine (Ser/Thr/Tyr) have also been found to play important roles in prokaryotic intracellular signalling. Although the most common group of Ser/Thr kinases (i.e. Hanks-type kinases) are often referred to as "eukaryotic kinases", many prokaryotic Ser/Thr kinases also belong to this group [16]. There are many examples of phosphorylation events regulating cellular activities in bacteria, relating to house-keeping functions as well as stress responses and virulence [12,17,18]. In Lactobacillus rhamnosus, phosphorylation of glycolytic enzymes is upregulated as a response to acidic stress [19]. In that study, the phosphorylation state of threonine and serine residues on glyceraldehyde 3-phosphate dehydrogenase (gapdh) differed between protein species depending on the growth conditions, while the overall abundance of the different proteins remained unchanged. This supports the idea that phosphorylation can regulate enzyme activity separately from regulation of gene expression. Different forms of the same protein may display increased or decreased catalytic activity, altered subcellular localization, or modified interaction with non-substrates [20,21]. Multisite phosphorylation can coordinate several such effects, determine the duration of a response, or mediate signal integration [21]. Protein phosphorylation may also completely alter the biological functions of proteins in the cell, a phenomenon referred to as protein moonlighting [22]. The regulatory effects of phosphorylation events on specific enzymes have been studied in some oral bacteria, mainly the oral streptococci [23][24][25], but to the best of our knowledge, global Ser/Thr/Tyr phosphorylation profiles of S. gordonii have not yet been detailed. The aim of the current study was to identify the intracellular protein expression profile, with special attention to Ser/Thr/Tyr phosphorylated proteins, in S. gordonii DL1.

General protein expression profile
Intracellular proteomes from planktonic cells of S. gordonii DL1 were extracted and separated by 2DE. The total intracellular proteome was visualized with Coomassie Brilliant Blue stain (Sigma). In total, 222 protein spots were detected. Discrete spots were manually excised for identification with LC-MS/MS. Molecular weights (MWs) and isoelectric points (pIs) for identified proteins were estimated from the gels as well as gathered from the mass spectrometry data. MWs, pIs, MASCOT scores, number of matched peptides and % coverages are listed in Table 1. The sequenced peptides identified by LC-MS/MS can be found in supplemental material (Additional file 1). In total, 61 proteins were identified, many of which were present at multiple isoelectric points ( Fig. 1, Table 1).

Intracellular Ser/Thr/Tyr phosphorylation profile
Spots containing Ser/Thr/Tyr phosphorylated proteins were visualized with Pro-Q Diamond stain and phosphoimaging. The total intracellular proteome was visualized with T-Rex protein labelling for orientation of the phosphorylated spots. In total, 49 phosphorylated spots were detected. These were generally found on the acidic side of the 2DE gels (Fig. 2).
Phosphorylated protein spots where identified with LC-MS/MS ( Table 1). Six of the 49 phosphorylated spots remained unidentified, due to a lack of significant hits from the Mascot search (four spots) or absence from Coomassie gels, preventing excision for identification (two spots). In total, 19 putative Ser/Thr/Tyr phosphorylated proteins were identified, see Table 1.

Cellular processes associated with Ser/Thr/Tyr phosphorylated proteins
The proteins detected as Ser/Thr/Tyr phosphorylated are involved in various cellular processes (Fig. 3). All phosphoproteins involved in the carbon metabolism, except for HPr and tpi, were present as more than one phosphorylated species, occurring at different pIs on the gels. HPr and tpi were also present as one additional, non-phosphorylated species respectively (Figs. 1 and 2, Table 1). Pfk and pgm were both present as two species each, all phosphorylated. Enolase was identified in both unphosphorylated and phosphorylated forms. Four of the phosphorylated enolase species were found at the expected MW. Pyk was present as four adjacent species, two of which were phosphorylated. Eight spots were identified as gapdh, three of which were phosphorylated. Nine species of pgk were identified, of which three were phosphorylated. Three species of ldh were identified, two phosphorylated. Six phosphorylated proteins involved in biosynthesis were identified. Deob, pgam and             Two studies on the model organism Bacillus subtilis were also found and included (Table 2). From planktonic cells grown under varying conditions, the previous studies revealed in total 27 and 29 putative Ser/Thr/Tyr phosphorylated proteins in B. subtilis, 15 in L. rhamnosus, 51 in Neisseria meningitidis, 73 in Staphylococcus aureus, and between 10 and 26 in different species of mycoplasma. All phosphoproteins identified in the current study except for three (pfk, rmla, ppp), were also detected as Ser/Thr/Tyr-phosphorylated in those studies, with varying occurrences. The number of phosphorylated proteins in S. gordonii DL1 detected in this study make up approximately 1% of the total number of annotated proteins in the UniProt database [45] for this species. This corresponds well to the other species tested with similar methodology, whose Ser/Thr/Tyrphosphoproteomes as detected with Pro-Q Diamond stain varied between 0.7 to 2.6% of the total proteomes. Confirmed or putative Ser/Thr/Tyr phosphorylation sites in all proteins identified as phosphorylated in this study or homologues, except for rmla and ppp, were found in the literature ( Table 2).  (Table 2).

Discussion
Ser/Thr/Tyr phosphoproteome of S. gordonii DL1 Putative Ser/Thr/Tyr phosphorylated intracellular proteins in S. gordonii DL1 were identified from 2DE gels using Pro-Q Diamond stain and LC-MS/MS. From the 49 phosphorylated spots detected, 19 Ser/Thr/Tyr phosphoproteins were identified in S. gordonii. Determined or putative phosphorylation sites on Ser, Thr and/or Tyr in all phosphorylated proteins or homologues except for two, rmla and ppp, were found in the literature (Table 2).
A variety of cellular processes were associated with the phosphoproteins identified in S. gordonii DL1 (Fig. 3). This supports the notion that Ser/Thr/Tyr phosphorylation events are integrated in the pathways that regulate different cellular responses in this species. Many of the phosphorylated proteins were present at two or more pIs in S. gordonii. As is evident from Table 2, this is a common occurrence in other bacteria as well. Variably phosphorylated species of proteins are of interest because they may differ in function [21] and hence be of biological relevance. With 2DE, although laborious, the occurrence of physiologically distinct forms of one protein can be readily visualized in a more "hands on", comprehensible overview of the present proteomic profile, were differently modified protein species often occur at adjacent pIs. Studies on deletion mutants can be employed to investigate the role of specific kinases/phosphatases on phenotypic presentation or different  Table 1). Phosphorylated spots with no significant hits in the Mascot database search are labelled with an asterisk (*). Gels were produced in triplicates from three different cultures of S. gordonii DL1 bacterial functions [46]. The current material constitutes a first step towards future studies that may utilize such targeted approaches. Studies of this character can be carried out from any perspective of microbial physiology and biological interaction of interest. Future studies using clinical isolates as well as biofilm growth flow cell models with saliva as a nutritious substrate are also of interest to increase the resemblance to in vivo conditions.

Core phosphoproteome
The total number of identified Ser/Thr/Tyr phosphorylated proteins as detected with the Pro-Q Diamond stain in other bacteria grown in various growth conditions was between 10 and 73 [19,27,28,32,33,37,40,42], making up between 0.7 to 2.6% of the total proteomes as annotated in the UniProt database [45]. This is in accordance with the 19 phosphoproteins identified in S. gordonii DL1 (approximately 1% of the total proteome). All phosphoproteins identified in the current study except for three (pfk, rmla, ppp) were also previously detected as Ser/Thr/Tyr-phosphorylated with the Pro-Q Diamond stain in other bacteria ( Table 2). Proteins that are consistently found in their phosphorylated state across several species and growth conditions may represent a core phosphoproteome profile shared by many bacteria. The mapping of such a core phosphoproteome may facilitate the identification of phenotypic characteristics that deviate from this core pattern. Identification of distinct phenotype phosphorylation patterns that reflect microbial activities may be crucial for pursuing a further understanding of biofilm formation and colonization of the commensal microbiota. Based on our findings, pfk, rmla and ppp represent possible proteins of interest for further investigation. Two phosphorylated protein species of pfk were detected in S. gordonii. In Streptococcus pneumoniae, pfk was found to be phosphorylated on a tyrosine residue, and phosphorylation of this enzyme was suggested to be involved in regulation of the metabolic flux in the cell [31]. One phosphorylated protein species of ppp was detected in S. gordonii. The phosphatase stp1 in Streptococcus agalactiae was identified as a homologue to S. gordonii ppp (100% similarity) in the UniProt.org database [45]. Stp1 is involved in group B streptococcal virulence [47]. Two phosphorylated protein species of rmla was identified in S. gordonii. Activation and deactivation of rmla is an important regulating mechanism of peptidoglycan biosynthesis in gram positive bacteria, making it an interesting candidate in the search for new antibiotics [48]. Hardly anything can be found in the literature about potential regulation of catalytic activity by phosphorylation of rmla, but the current findings suggest that Ser/Thr/Tyr phosphorylation events are involved. These three proteins represent possible candidates for uniquely phosphorylated proteins in S. gordonii, however, because phosphorylation patterns are dependent on growth conditions, further experiments comparing the phosphoproteomes of different species grown in the same environments are needed.

Phosphoproteins in S. gordonii DL1 involved in carbohydrate transport and metabolism
The carbohydrate transporting protein HPr was detected as Ser/Thr/Tyr phosphorylated in S. gordonii. HPr is a known phosphoprotein with roles in the phosphotransferase system (PTS), for uptake and phosphorylation of certain carbohydrates in the central carbon metabolism.  The PTS has a high affinity for glucose, and a preference for glucose over alternative carbohydrates. HPr carries known phosphorylation sites on histidine, per se involved in the phosphotransfer catalytic activity, and serine, involved in the regulation of catalytic phosphotransfer activity [49]. Phosphorylation on HPr serine reduces the phosphotransfer activity, thereby reducing the uptake of PTS-carbohydrates, while simultaneously enhancing the uptake of alternative carbohydrates [7,26,49]. In this way, phosphorylation events on HPr can regulate both the quantity and types of carbohydrates that are processed by the cell. Seven glycolytic proteins (tpi, pfk, gapdh, pgk, pgm, eno, pyk), as well as ldh, involved in pyruvate conversion, were Ser/Thr/Tyr phosphorylated in S. gordonii. Prokaryotes have been suggested to regulate glycolytic activity rapidly and reversibly through protein phosphorylation and dephosphorylation in response to environmental changes [12,13]. Our findings indicate that these mechanisms are present in S. gordonii as well. Aside from controlling the metabolic rate and character of the metabolic end products of glycolysis, e.g. by activation of ldh, regulation of some glycolytic enzymes, e.g. tpi, controls the switch between different metabolic pathways [50,51]. Regulation of enzymes with reversible catalytic function, e.g. pgk, eno, may enable switching between catabolism and anabolism [45]. Rapid and reversible regulation of carbohydrate transport, control of rate of glycolysis, and alternative carbohydrate metabolizing pathways by phosphorylation, increase cell fitness by enabling adaptation of the cell metabolism in response to variations in carbohydrate concentration. Thereby, damage from processes such as sugar killing can be mitigated [8,9,12]. Our findings indicate that Ser/Thr/Tyr phosphorylation is involved in regulating the activities of these proteins. Control of glycolysis not only by the levels of metabolic substrates and products, but by transduction of a variety of other signals through enzyme phosphorylation events, suggests a high complexity of metabolic control in S. gordonii. This supports the idea that integration of metabolic activity and environmental factors in bacterial cells is more convoluted than previously believed.

Phosphoproteins in S. gordonii DL1 involved in acid tolerance response
To cope with the intermittent fluctuations in pH that naturally occur in oral biofilms, partly as a consequence of accumulation of lactic acid and other acids from the central carbon metabolism, oral streptococci employ the socalled acid tolerance response (ATR) [4,10,11]. Aside from responsive regulation of carbohydrate transport and metabolism, increased molecular chaperone activity comprises another important aspect of the ATR in oral streptococci [4,8]. The chaperones represent the group of proteins most frequently detected as Ser/Thr/Tyr  (Table 2). In the current study, S. gordonii displayed three phosphorylated molecular chaperones (DnaK, GroEL, tig). Ser/Thr/Tyr phosphorylation has been suggested to activate all these three chaperones in bacteria [27,37,39,52]. By modulating the character and acidity of metabolic end products that accumulate locally, while retaining effective metabolism and other cellular functions at lower environmental pH with the help of ATR responses, cells increase their competitiveness towards less aciduric species [53]. From an evolutionary point of view, these mechanisms may have contributed to the long-term survival of these species in the strenuous oral environment.

Phosphoproteins in S. gordonii DL1 involved in biosynthesis
Six phosphorylated proteins involved in biosynthesis were identified in S. gordonii DL1 (deob, imdh, ef-tu, rS1, pgam, rmla). These proteins are involved in amino acid, nucleotide and cell wall biosynthesis, as well as protein translation (Table 1). Except for the chaperones, ef-tu was the protein most commonly detected with Pro-Q Diamond stain on 2DE gels in the previously studied bacteria (Table 2). Phosphorylation on a threonine residue in Escherichia coli ef-tu seems to inactivate the protein by decoupling the tRNA from the ribosome A-site [54,55], triggering a rapid inhibition of protein translation. In contrast, phosphorylation on a threonine residue in the ribosomal protein rS1 in E. coli was found to activate protein translation [56]. Rapid and reversible regulation of biosynthesis through protein phosphorylation events enables the cell to be frugal regarding energy expenditure, and streamline cellular functions by modification of biomolecule production.

Moonlighting and multisite regulation of identified phosphoproteins
Several of the proteins identified as phosphorylated in S. gordonii DL1 have known moonlighting functions, often in adhesion, e.g. tpi [57], gapdh [58], eno [59], and GroEL [60], alternative metabolic pathways, e.g. gapdh in iron metabolism [61], or other functions. In prokaryotes, cell surface associated enolase is involved in adhesion to host components in connective tissue or saliva [59], and retention of exported enolase on the bacterial cell surface was found to increase in acidic environments [58]. Maintaining cell adhesion during acidic shifts that may inactivate other adhesins is a tentative strategy for increased competitiveness in oral biofilms. In E. coli enolase, apart from reducing its glycolytic activity, lysine phosphorylation also prevented the enzyme export related to its moonlighting functions [62]. These findings suggest that phosphorylation of enolase may affect cell adhesion during acidic shifts. Regulation of several of the phosphoproteins identified in S. gordonii DL1, such as chaperones DnaK and GroEL, also seem to involve multisite phosphorylation events [30,34,39,63]. This further supports the idea that regulation of cellular activities in prokaryotes is complex and involves integrated patterns of signal transduction.

Conclusion
This study clearly shows that Ser/Thr/Tyr phosphorylation is present in an array of cytoplasmic proteins from the oral commensal S. gordonii. In total, 61 intracellular proteins were identified in S. gordonii DL1, and 19 of these turned out to be present as phosphorylated. Most of the phosphorylated proteins are involved in the carbon metabolism, specifically related to glycolysis, which is no surprise for this saccharolytic oral streptococcus. Ser/Thr/Tyr phosphorylation presents a possible mechanism for regulation of multiple cellular processes in S. gordonii, and phosphorylated proteins species were often present at several pIs with potential variation in biological function. Many similarities were found between the identified phosphoproteome of S. gordonii DL1 and that of previously studied species, despite differences in basic cell physiology and the growth conditions applied. Identification of core regulatory pathways involved in the interaction between bacteria and their environment should provide useful insights regarding new strategies to manage biofilm-induced diseases. Studies in microbiology are often focused on investigating the role of bacteria in disease. However, the physiology of commensals that prevent shifts towards dysbiosis are just as relevant. To investigate mechanisms of biofilm formation and homeostasis in the oral cavity, it is essential to examine bacterial responses and regulation of responses in relation to specific environmental challenges.

Protein extraction
Harvested cells were subjected to three freeze-thaw cycles before washing and then resuspended in 700 μl lysis buffer containing 8 M urea, 2% (v/v) CHAPS, 64.8 mM DTT, 2% IPG buffer pH 4-7 (Pharmacia Amersham Biotech, Sweden). Ultrasonication of the samples was then performed with homogenizing 0.2 mm glass beads for 4 × 5 min (5 s pulses, amplitude 40, Vibra-Cell™ Ultrasonic Processor, SONICS), with alternate periods of cooling. Intact cells were sedimented by centrifugation at 17000×g for 10 min at 4°C and the supernatants (protein extracts) were stored at − 20°C. The protein concentration was determined using the 2-D Quant kit (GE Healthcare Life Sciences).
Gels for fluorescent staining were loaded with 50 μg protein, and proteins were pre-labelled with the T-Rex Labelling kit (NH DyeAGNOSTICS) before rehydration (see "Staining procedures" below). The gels were then fixed in the appropriate fix solution as recommended by the respective stain manufacturer (Coomassie stain 2% acetic acid and 40% ethanol and fluorescent stain 10% acetic acid and 50% ethanol). Gels for Coomassie staining were fixed for at least 1 h and gels for fluorescent staining were fixed overnight, protected from light. All gels were produced in triplicates from separate cell cultures.

Staining procedures
Coomassie gels were stained overnight in 300 ml Coomassie brilliant blue staining solution containing 17% Coomassie Brilliant Blue G -Colloidal Concentrate (Sigma) and 21% ethanol. After staining, Coomassie gels were destained with 25% ethanol for approximately 1 h and stored in Ultra-High Quality water (UHQ) at 4°C until scanning and excision of spots. For fluorescent staining, lysine residues of the general proteome were pre-labelled with T-Rex Labelling kit (NH DyeAGNOS-TICS) before rehydration, to facilitate orientation of the phosphorylation profiles. Pre-labelling was performed according to the manufacturer's protocol. Briefly, 2 μl T-Rex solvent and 50 μg of the extracted proteins were transferred into a T-Rex labelling vial and incubated on ice for 30 min. Thereafter the gels were prepared as described above. After fixation, each T-rex stained gel was washed with UHQ for 3 × 15 min and pre-scanned once at 647 nm with a photomultiplier tube (PMT) setting of 800, pixel size 100 μm, before staining with Pro-Q Diamond Phosphoprotein Gel Stain (Invitrogen™, Thermo-Fisher Scientific) according to the manufacturer's instruction. Phosphoprotein staining was performed protected from light, and with minor modifications based on results from Agrawal and Thelen [65]. Each gel was stained with Pro-Q Diamond Phosphoprotein Gel Stain diluted 1:2 in UHQ to a final volume of 250 ml for 2 h. Gels were destained in 250 ml Pro-Q Diamond Phosphoprotein Gel Destaining Solution (Invitrogen™, Thermo-Fisher Scientific) per gel for 3 × 30 min on a shaker and then washed in UHQ for 3 × 5 min.

Phosphoimage analysis
Phosphoimaging was performed with a Fujifilm FLA-9000 (Science Imaging Scandinavia AB). T-Rex labelled proteins were visualized at 647 nm, PMT 800, 100 μm (displayed as red). Pro-Q Diamond stained proteins were visualized at 532 nm, PMT 600, 100 μm (displayed as green). Phosphorylation profiles were analysed by three operators and then compared and discussed for calibration, with good concordance. Quantification of Pro-Q Diamond Phosphoprotein Gel Stain signals from phosphorylated spots was not possible due to variations between triplicates.

Identification of streptococcal proteins by LC-MS/MS
Protein spots were manually excised from Coomassie brilliant blue-stained gels and subjected to in-gel digestion with trypsin as previously described [64]. The resulting protein fragments were separated with liquid chromatography (LC) and characterized using tandem mass spectrometry (MS/MS) (Aberdeen Proteomics, University of Aberdeen). In short, proteins were treated with DTT for reduction (60°C, 20 min), iodacetamide for S-alkylation (25°C, 10 min) and trypsin for digestion (37°C, 8 h). Following drying by rotary evaporation (SC110 Speedvac, Savant Instruments), the peptide extract was dissolved in 0.1% formic acid. Analysis of peptide solutions was carried out using an HTCultraPTM Discovery system (Bruker Daltonics) coupled to an Ul-tiMate 3000 LC system (Dionex). Separation of tryptic peptides was performed on a monolithic capillary column (200 μm internal diameter × 5 cm, Dionex). The gradient consisted of 5% acetonitrile in UHQ containing 0.5% formic acid, to 37% acetonitrile in UHQ containing 0.45% formic acid over 12 min at a flow rate of 2.5 μl min − 1 . Data-dependent mode was employed to acquire peptide fragment mass spectra, AutoMS (2). The scan range was 300-1500 m/z, with three averages, and up to three precursor ions that were selected from the MS scan (100-2200 m/z). Active exclusion of precursors was performed within a 1.0 min window, as well as exclusion of all singly charged ions. Detection and deconvolution of peptide peaks were carried out automatically using Data Analysis software (Bruker). Mass lists in the form of Mascot Generic files were created automatically and used as the input for Mascot MS/MS ion searches in the NCBInr database for Firmicutes (Gram Positive Bacteria) using the Matrix Science web server [66]. Parameters were set to 0.5 Da peptide mass tolerance, methionine oxidation and carboxyamidomethylation of cysteine. The maximum number of missed cleavages was set to 1. Proteins were identified by at least two peptides with a Mascot score of 33 or higher, the latter as suggested by Koenig et al., 2008 [67].

Literature search
The PubMed, Web of Science and Cochrane databases were searched for studies that examined global Ser/Thr/ Tyr phosphoproteomes in human commensal and pathogenic bacteria using Pro-Q Diamond Phosphoprotein Gel Stain on 2DE gels. Studies on the model organism Bacillus subtilis were also included.
Additional file 1. The sequenced peptides identified by LC-MS/MS can be found in supplemental material (Additional file 1). Sequenced peptides identified by LC-MS/MS; Ser/Thr/Tyr phosphorylated protein species (Table A), and non-phosphorylated protein species (Table B).