Previous Article | Next Article ![]()
Eukaryotic Cell, March 2006, p. 457-468, Vol. 5, No. 3
1535-9778/06/$08.00+0 doi:10.1128/EC.5.3.457-468.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Institut für Allgemeine Botanik, Friedrich-Schiller-Universität Jena, Am Planetarium 1, 07743 Jena,1 Thermo Electron GmbH, Im Steingrund 4-6, 63303 Dreieich, Germany2
Received 6 November 2005/ Accepted 20 December 2005
|
|
|---|
|
|
|---|
Our current approach is to apply proteomics to the investigation of posttranslational modifications in C. reinhardtii. One of the key modifications of proteins, which is crucial in the control of many regulatory pathways and can effect protein function, activity, stability, localization, and interactions, is phosphorylation. Therefore, information about the phosphoproteome of C. reinhardtii is extremely useful for understanding a variety of cellular processes in this alga, which may serve as a basis for examining these processes in other organisms. Several phosphoproteins of C. reinhardtii, such as the
heavy chain of outer arm dynein (16) and IC138, a WD repeat dynein intermediate chain (10, 15), have already been identified. Also, several radial spoke proteins of the flagella have been characterized as phosphoproteins via in vivo 32P labeling (26). Metabolic enzymes such as cytosolic glutamine synthetase (27) and proteins of the photosynthetic machinery (e.g., Cp29-like protein [12]) of C. reinhardtii have also been identified as phosphoproteins.
Proteome analysis of phosphoproteins is a challenging task (20, 28). Phosphoproteins can possess more than one phosphorylation site, and the phosphorylation status of these sites can fluctuate depending on the physiological conditions under which the cells are kept. This leads to a great variety of phosphoproteins. In addition, the ratio of the phosphorylated to nonphosphorylated form of a protein can be very low. Although proteins can be identified down to the femtomole, and even attomole, level with modern mass spectrometry (MS), many phosphoproteins within a crude extract (especially those of cell signaling pathways) are not abundant enough to be unambiguously identified by MS. For this reason, enrichment of such proteins is often a prerequisite for efficient phosphoproteome analysis.
Different methods that can be used for this purpose have been described in the literature (reviewed in references 11 and 28). One of them, immobilized metal-ion affinity chromatography (IMAC), is based on the presence of negatively charged phosphate groups and enriches for phosphorylated Ser, Thr, and Tyr. This method has been applied, for example, to the analysis of phosphoproteins from yeast (4, 7) and from a lymphoma cell line (30). The IMAC method relies on direct identification of phosphopeptides by MS, in contrast to other methods that chemically substitute for the phosphate residue (11, 28). However, in tandem MS (MS2), phosphopeptide precursor ions can exhibit a neutral loss of phosphoric acid (98 Da). The reason for this loss is that phosphopeptides (phosphoserine and phosphothreonine) can undergo gas-phase ß elimination when subjected to collision-induced fragmentation (20). Because m/z (mass/charge) values, and not absolute masses, are measured in a mass spectrometer, doubly and triply charged peptide ions show an apparent loss of 49 and 32.66, respectively. In the MS2 spectrum, the presence of a neutral loss therefore indicates phosphorylation and can be used as a selection parameter for phosphopeptides.
Recently, a new type of linear-ion-trap electrospray ionization (ESI)-MS was developed, permitting the acquisition of data-dependent neutral-loss experiments. With this instrument (Finnigan LTQ; Thermo Electron Corp., San Jose, CA), neutral-loss analysis can be executed during MS measurements (MS2 scan). If a pair of peaks of the most prominent ion of the MS2 spectrum, versus the full-scan MS spectrum, is found with a mass difference of 98, 49, or 32.66 (depending on the charge of the peptide ion), a phosphorylation-specific neutral loss is indicated. This prominent ion of the MS2 spectrum (representing the dephosphorylated peptide) is then selected for an automatically triggered MS3 scan. Thus, sufficient fragment ion information with a significant cross-correlation (Xcorr) factor (2) can be obtained. The Xcorr factor describes the cross-correlation between the experimentally measured MS2/MS3 spectrum and an in silico-generated MS2/MS3 spectrum of candidate peptides in the databases. The data-dependent neutral-loss-triggered MS3 mode of operation has been recently applied successfully to the study of phosphoproteins from the yeast signaling pathway (7).
So far, relatively little information about phosphoproteomes in plants is known. Among the few examples, phosphoproteins of Arabidopsis plasma membranes (23) and phosphoproteins in the moss Physcomitrella patens (9) have been characterized. Here, we have enriched phosphoproteins from C. reinhardtii by using IMAC in combination with Ga3+ as the metal component. For this purpose, cells were grown in presence of a strong inhibitor of Ser/Thr phosphatases (okadaic acid) to increase the phosphorylation status of the proteins. It was previously shown that the growth of cells in the presence of such an inhibitor can robustly increase phosphorylation on specific sites (21), and such inhibitors have already been used successfully for phosphoproteome analysis (30). Using the information from MS2 and neutral-loss-triggered MS3 spectra, we identified 328 phosphoproteins in C. reinhardtii with 360 phosphopeptides. These results provide insight not only into phosphorylation events within various metabolic and signaling pathways but also into tactic movements, based on comparative proteome analysis with proteins from the centriole and the flagella.
|
|
|---|
Preparation of crude extracts of C. reinhardtii. Cells from a 100-ml culture were harvested by centrifugation and stored in liquid nitrogen. Extract preparation was based on the protocols of Shu et al. (30) and Mittag (22) with some modifications. Cells were resuspended in TriPure isolation reagent (Roche Applied Science), and glass beads (diameter, 0.25 to 0.30 mm) that had been prewashed once with TriPure isolation reagent were added. Cells were disrupted by vortexing this mixture five times for 1 min each, and the tubes were placed on ice for 1 min between each vortexing step. After a centrifugation step (200 x g, 2 s), the supernatant was subjected to the TriPure procedure for removal of RNA and DNA, according to the Roche manual. Proteins were precipitated with isopropanol, and the pellet was washed three times with 0.3 M guanidine hydrochloride in 95% ethanol and finally with 100% ethanol. A 6 M guanidine hydrochloride solution (100 µl) was added to dissolve the protein pellet. For tryptic digestion, the solubilized proteins were diluted 10-fold in 100 mM ammonium bicarbonate and the protein concentration was determined with the Bio-Rad protein assay. The digest was performed with 20 µg of trypsin (Promega) and 1 mg of solubilized protein overnight at 37°C. Then, acetonitrile and formic acid were added to the sample to reach final concentrations of 2% (vol/vol) acetonitrile and 0.1% (vol/vol) formic acid. The tryptic peptide solution was centrifuged (10,000 x g, 4°C, 5 min), and the supernatant was put on a fast protein liquid chromatography (LC) system (Amersham Biosciences)-driven 1-ml column (SOURCE RPC; Amersham Biosciences [now part of GE Healthcare]) for desalting. Peptides were eluted with 500 µl of 80% (vol/vol) acetonitrile-0.1% (vol/vol) formic acid. The desalted sample was subjected to the IMAC procedure.
Enrichment of proteins by IMAC. The IMAC procedure was performed according to Shu et al. (30) with some modifications. For IMAC, we used homemade microcolumns prepared according to the procedure of Erdjument-Bromage et al. (3). Fifty microliters of POROS 20 MC (catalog no. 1-5428-02; Applied Biosystems) metal-chelating resin (66% [wt/wt] slurry) was pipetted into a microtip column (Eppendorf gel loading tip). Then, two 50-µl volumes of 0.1% (vol/vol) acetic acid were passed through the column. The charging of the metal-chelating resin was carried out by applying 150 µl of 100 mM GaCl3. To remove excess GaCl3, the column was washed with 50 µl of 0.1% (vol/vol) acetic acid. Afterward, the tryptically digested and desalted peptide solution (see above) was loaded onto the activated IMAC column. The column was washed successively with 50 µl of 0.1% (vol/vol) acetic acid, 50 µl of 50% (vol/vol) acetonitrile-0.1% (vol/vol) acetic acid, 50 µl of 50% (vol/vol) acetonitrile-0.1% (vol/vol) acetic acid-100 mM sodium chloride, and finally 50 µl of 0.1% (vol/vol) acetic acid. Then, phosphopeptides were eluted with three 20-µl volumes of 200 mM Na2HPO4 and desalted with a ZipTip microtip (Perseptive Biosystems) according to the method of Stauber et al. (31). The ZipTip contained C18 reversed-phase chromatographic medium. In addition, a 50-µl mixture of POROS R2 (Applied Biosystems) and methanol (1:5) was loaded, and the column was washed two times with 5% (vol/vol) methanol-5% (vol/vol) formic acid. To improve the binding of the phosphopeptides to the ZipTip microtip, formic acid was added to the sample to a final concentration of 5% (vol/vol) prior to loading. Then, the column was again washed two times with 50 µl of 5% (vol/vol) methanol-5% (vol/vol) formic acid. Peptides were eluted with two 50-µl volumes of 60% (vol/vol) methanol-5% (vol/vol) formic acid. The resulting 100-µl sample solution was dried in a SpeedVac concentrator (about 2 to 3 h), and the pellet was then stored at 70°C.
Peptide identification by nano-LC-ESI-MS (MS2 and neutral-loss-triggered MS3). After the dried pellet was dissolved in 5 µl of 5% (vol/vol) acetonitrile-0.1% (vol/vol) formic acid, the phosphopeptides were subjected to LC-ESI-MS using a nanoscale C18 column (flow rate, 600 nl/min) coupled online with a linear-ion-trap mass spectrometer (Finnigan LTQ; Thermo Electron Corp., San Jose, CA). The LC system consisted of a Surveyor high-performance liquid chromatograph (Thermo Electron), including a Surveyor autosampler and MS pump. A gradient was used to elute peptides from the reversed-phase column (Picofrit) (length, 10 cm; inner diameter, 100 µm; Magic-C18 particle size, 5 µm [Spectronex]). The successive steps of the applied gradient were as follows: 10 min, 100% (vol/vol) A; 10 min, gradually shifting to 70% (vol/vol) A-30% (vol/vol) B; 20 min, gradually shifting to 40% (vol/vol) A-60% (vol/vol) B; 10 min, gradually shifting to 20% (vol/vol) A-80% (vol/vol) B; 10 min, 20% (vol/vol) A-80% (vol/vol) B; and 1 min, gradually shifting to 100% (vol/vol) A, where A consists of 0.1% (vol/vol) formic acid in water and B consists of 0.1% (vol/vol) formic acid in acetonitrile. The instrument was run by the data-dependent neutral-loss method, cycling between one full MS scan and MS2 scans of the four most-abundant ions. The detection of a neutral-loss fragment (98, 49, or 32.66 Da) in the MS2 scans triggered immediately an MS3 scan of the precursor ion representing the dephosphorylated peptide. The MS2 and MS3 data were used to search the available Chlamydomonas reinhardtii databases (see below) using Bioworks software (version 3.1; Thermo Electron Corp., San Jose, CA) including the SEQUEST algorithm (19). The software parameters were set to detect a modification of 79.96 Da in Ser, Thr, or Tyr in MS2 and MS3 spectra. When phosphoserine and phosphothreonine undergo gas-phase ß elimination, dehydroalanine (Dha) and 2-aminodehydrobuturic acid (methyldehydroalanine [MeDha]), respectively, are produced. Thus, modifications of 18.00 Da in Ser and Thr residues were additionally used for database searches with MS3 data. Searches were done for tryptic peptides, allowing two missed cleavages. Mass tolerance was set to 1.5 Da for the peptide precursor ion in MS mode. For fragment ions (MS2 and MS3 modes), mass tolerance was set to 0.5 Da. It should be noted that Bioworks software (version 3.1) includes masses of fragment ions up to 1 Da when this setting is applied. Scores for the Xcorr factor (2) were set to the following limits: Xcorr of >1.5 if the charge of the peptide was 1, Xcorr of >2 if the charge of the peptide was 2, and Xcorr of >2.5 if the charge of the peptide was 3. Peptide searches were done with the following databases: genomic databases, including their predicted proteins (version 2), available from the Joint Genome Institute (JGI) facility (http://genome.jgi-psf.org/chlre2/chlre2.home.html); the Kazusa EST database (http://www.kazusa.or.jp/en/plant/chlamy/EST); the Chlamy Center EST database (http://www.chlamy.org) (31); the NCBI mitochondrion database (gi|11467088|ref|NC_001638.1|); and the chloroplast database (http://www.chlamy.org/chloro/default.html). As Bioworks software (version 3.1) ignores stop codons from translated DNA sequences of genome and EST databases, an in-house program was developed to automatically delete such false positives.
Phosphothreonine immunoblotting. Lysates from control and okadaic acid-treated cells were immunoblotted with an antiphosphothreonine antibody. C. reinhardtii cells were either left untreated or treated with 1.5 µM okadaic acid for 29 h under dim light. Crude extracts were prepared as described previously (22), and proteins were denatured in sodium dodecyl sulfate (SDS) sample buffer (30). The proteins (50 and 100 µg, respectively) were resolved by 9% SDS-polyacrylamide gel electrophoresis, transferred to a nitrocellulose membrane, and immunoblotted with antiphosphothreonine antibody, according to the manufacturer's protocol (catalog no. 9381; Cell Signaling Technologies) with one modification: the commercially available soya-based powder Slimfast was used instead of nonfat dried milk powder in the blocking buffer.
|
|
|---|
To analyze whether treatment with okadaic acid for 29 h would significantly increase the levels of phosphoproteins in C. reinhardtii, we prepared crude extracts from okadaic acid-treated and untreated cells and analyzed proteins by immunoblotting with antiphosphothreonine antibodies (for details, see Materials and Methods). Immunoblotting showed that okadaic acid significantly increases Thr phosphorylation in C. reinhardtii (Fig. 1). An increase in the level of phosphoproteins in okadaic acid-treated versus untreated cells was also observed when the yield of phosphoproteins within these two phosphoproteomes was compared (see "Identification of phosphopeptides from C. reinhardtii by nano-LC-ESI-MS," below).
![]() View larger version (90K): [in a new window] |
FIG. 1. Enrichment of phosphorylated proteins by okadaic acid treatment of C. reinhardtii. Cells were either left untreated () or treated with 1.5 µM okadaic acid (+), as described in Materials and Methods, harvested, and used for preparation of a crude extract. Proteins (50 and 100 µg) were denatured in SDS sample buffer, resolved on a 9% SDS gel, transferred to a nitrocellulose membrane, and immunoblotted with antiphosphothreonine polyclonal antibody (Cell Signaling Technologies).
|
![]() View larger version (22K): [in a new window] |
FIG. 2. Schema for the enrichment and analysis of phosphopeptides from C. reinhardtii. The details are described in Materials and Methods. Nano-LC-ESI-MS (MS2 and MS3) analysis was carried out in a mass spectrometer with a linear ion trap, permitting the acquisition of data-dependent neutral loss (Finnigan LTQ; Thermo Electron Corp., San Jose, CA).
|
|
View this table: [in a new window] |
TABLE 1. Identification of phosphopeptides from beta-casein, which were added to the total phosphopeptide sample from C. reinhardtii, by LC-ESI-MS2 and neutral-loss-triggered MS3
|
In four independent experiments, cells were grown in the presence of 1.5 µM okadaic acid for 29 h and harvested, and the phosphopeptides were obtained by the procedure described in Fig. 2. To increase the confidence of the phosphopeptide assignment, we chose only phosphopeptides that fulfilled three criteria for further analysis. (i) They appeared at least in two of the four independent experiments. (ii) They showed a significant Xcorr factor (see Materials and Methods). (iii) They could be identified (a) within the protein sequence of a predicted gene model, (b) within the ORF of an EST assembly, or (c) within a potential part of an ORF (potential exon) predicted after translation of the genome sequence itself. The latter (c) was included only if this potential exon had significant homology (a BLAST E score of
1 x 105) to an already known protein from any other organism. Thus, the yield of the phosphopeptides correlates to a high degree with the current number of gene models predicted by JGI from the available genome sequence (version 2) and the available EST data.
Using these criteria, we were able to identify 186 phosphoproteins of known and putative function (see Table S1 in the supplemental material) and 142 novel phosphoproteins of as-yet-unknown function (see Table S2 in the supplemental material). Certain phosphopeptides could be identified solely by their MS2 spectra with a significant Xcorr factor. These include almost all phosphotyrosine-containing peptides but also Ser/Thr phosphorylated peptides. Some phosphopeptides were identified by MS2 spectra and additionally by the corresponding MS3 spectra of their neutral-loss peptide (indicated in bold in Tables S1 and S2 in the supplemental material). An example is depicted in Fig. 3. A detailed view of the full MS scan with the precursor ion m/z 914.6 is shown in Fig. 3A. The MS2 spectrum of the peptide ion showing the neutral loss led to identification of the phosphopeptide RWPEDGINTpEDMVWLAEEANKR with an Xcorr factor of 5.54 (Fig. 3B). Analysis of the MS3 fragmentation spectrum from the peptide ion m/z 881.6 revealed the same peptide with the 2-aminodehydrobuturic acid residue (MeDha) as a result of the neutral-loss event with an Xcorr factor of 8.0, showing nearly a complete set of the corresponding y- and b-ion series (Fig. 3C).
![]() View larger version (28K): [in a new window] |
FIG. 3. Identification of a phosphopeptide from a MORN (for "membrane occupation and recognition nexus") repeat protein (gene model C_190173) by nano-LC-ESI-MS2 and by neutral-loss-triggered MS3. (A) Full MS scan in the m/z range of 870 to 960 with the prominent peptide ion 914.6. (B) Identification of the phosphopeptide RWPEDGINTpEDMVWLAEEANKR by the MS2 fragmentation pattern of peptide ion 914.6. "p" indicates the phosphorylation site on Thr. (C) Identification of the peptide RWPEDGIN(MeDha)EDMVWLAEEANKR by the MS3 fragmentation pattern of the detected neutral-loss fragment 881.6. "MeDha" indicates the site of the neutral loss of phosphoric acid from the phosphothreonine. Only prominent y- and b-fragment ions have been labeled.
|
![]() View larger version (21K): [in a new window] |
FIG.4. Identification of a phosphopeptide from a novel protein of unknown function (gene model C_50124) by neutral-loss-triggered nano-LC-ESI-MS3. (A) Full MS scan in the m/z range of 500 to 600 detects the prominent peptide ion 546.7. (B) The MS2 fragmentation spectrum of peptide ion 546.7 reveals the neutral loss (fragment ion 497.4) but no fragment information sufficient for peptide identification. (C) Identification of the peptide KAAGFFG(MeDha)GR by neutral-loss-triggered MS3 of peptide ion 479.4. "MeDha" indicates the site of the neutral loss of phosphoric acid from the phosphothreonine. Only prominent y- and b-fragment ions have been labeled.
|
Functional analysis of the phosphoproteome of C. reinhardtii.
We searched for the functions of all phosphoproteins by using the information on the JGI website, including annotation notes as well as information about conserved domains. Further, we did an NCBI homology search (NCBI protein BLAST). An overview of the output of the functional analysis is shown in Table 2, and a detailed view is presented in Tables S1 and S2 in the supplemental material. Thereby, proteins were labeled as novel proteins of as-yet-unknown function (see Table S2 in the supplemental material) if (i) the protein had not yet been characterized for C. reinhardtii, (ii) it had either no or weak homology (a BLAST E score of
1 x 105) to a protein from any other organism, or (iii) it had significant homology to a protein from any other organism but the function of the protein in the other organism was unknown and no conserved domain could be found indicating any putative function.
|
View this table: [in a new window] |
TABLE 2. Functional analysis of identified phosphoproteins
|
We also performed a comparative proteome analysis with already known subproteomes from the chloroplast 70S ribosome (33, 34), thioredoxin interaction partners (18), the flagella (24), and the centriole (13). In this way, we were able to identify several proteins of these subproteomes that are subject to phosphorylation (Table 3). A few of them even appeared in more than one subproteome; these are indicated in Table 3. In summary, 1 protein of the chloroplast 70S ribosome, 7 from the thioredoxin targets, 4 from the centriole, and 32 from the flagella were detected. The relatively high yield of phosphoproteins from the flagella is reflected by higher phosphopeptide coverage (41 phosphopeptides in all). Some of the phosphoproteins of the flagella are related to certain human diseases, which are discussed below.
|
View this table: [in a new window] |
TABLE 3. Comparative proteome analysis with the phosphoproteome of C. reinhardtii
|
|
|
|---|
Addition of okadaic acid, an inhibitor of Ser/Thr phosphatases, was provided to increase the phosphorylation status of proteins that contain phosphoserine or phosphothreonine. Treatment of cells with the Ser/Thr phosphatase inhibitor calyculin had been already successfully applied for determining the phosphoproteome of a lymphoma cell line (30). In this case, it was shown that incubation of cells with the inhibitor led to a robust increase in the levels of phosphothreonine-containing proteins in comparison to untreated cells. By use of the inhibitor, the authors were able to identify a total of 107 phosphoproteins in a WEHI-231 B lymphoma cell line. Incubation of cells with okadaic acid had also been already studied, though with regard to a specific protein (human p53). It was shown thereby that the ratio of nonphosphorylated to specific phosphorylated amino acids from p53 was significantly shifted toward the phosphorylated form, permitting efficient determination of phosphorylation sites by LC-MS (21). For example, Ser315 of p53 had a ratio of 0.95 (nonphosphorylated form) to 0.05 (phosphorylated form) in wild-type cells and a ratio of 0 (nonphosphorylated form) to 1 (phosphorylated form) in okadaic acid-treated cells. This shows clearly the potential impact of inhibitor treatment in specifically increasing phosphorylation sites. Additionally, in the case of C. reinhardtii, the yield of phosphoproteins was significantly enhanced in the presence of okadaic acid in comparison to untreated cells. This was tested (i) by immunoblotting with antiphosphothreonine antibody and (ii) by comparing the yield of the phosphoproteins from untreated versus treated cells, in which it increased by 38%. Although okadaic acid enhanced the yield, phosphatase inhibition is not without biological consequences. In our case, treatment for 5 days caused defects in motility for about 75 to 80% of the cells and presumable cell death for about 10% of the cells. For this reason, we used an incubation time of 29 h, producing cell viability comparable to that of untreated controls.
To enrich the phosphopeptides in a complex trypsin-digested peptide mixture from C. reinhardtii, we used IMAC, which had already been successfully used to investigate, e.g., phosphoproteins from yeast (4, 7) and a lymphoma cell line (30). Since it had been reported (28, 30) that a change in the metal component from Fe3+ to Ga3+ could significantly increase the yield of phosphopeptides versus nonphosphorylated peptides, we used IMAC in combination with Ga3+ as the metal component. The internal phosphopeptide standard from beta-casein showed that MS2 spectra are sufficient to identify phosphopeptides, as was the case for the monophosphopeptide of beta-casein. However, it also became clear that in some cases the neutral-loss-triggered MS3 event can be crucial for the identification of a phosphorylation site, as was found for the tetraphosphopeptide of beta-casein. A similar situation was found for the identification of phosphopeptides from C. reinhardtii. Many of them could be identified based on the MS2 spectra alone or in combination with MS3 spectra. Notably, some phosphopeptides could be identified based only on the MS3 spectra, showing the advantage of the data-dependent neutral-loss method that automatically triggers MS3 scan events.
In all, 328 phosphoproteins, along with 360 phosphopeptides, from C. reinhardtii were identified. The yield of phosphoproteins was correlated directly with the current number of correctly predicted gene models (based on the genome sequence, version 2) and available EST data. Gene model predictions for Arabidopsis thaliana were reported by Peltier et al. (25) to contain significant errors for about 30% of the analyzed models, varying from incorrect N and C terminus predictions to errors in intron/exon boundary predictions and missing exons. In our case, 31% of the analyzed models were supported by EST data covering at least part of the exons of each gene model. However, 13% of the positive hits were based purely on EST data, since no gene models were available in these cases. This is due possibly to incorrect gene model predictions and certainly to the fact that parts of the genome of C. reinhardtii still contain regions with undetermined nucleotides (due to the high GC content of its DNA, which affects sequencing). It can be predicted that the yield of positive hits will increase after the final version of the C. reinhardtii genome has been released.
Our selection model provided rigorous criteria. We selected only phosphopeptides that occurred in at least two of the four independent experiments and had significant Xcorr factors (see Materials and Methods). Further hits that covered parts of potential exons within the genomic sequence of C. reinhardtii but showed no or only weak homology to currently known proteins and were not supported by gene models or EST data were neglected. It is possible that they have been omitted from the gene models due to errors in gene model prediction, but they may also represent false positives.
It is estimated that about 30% of all proteins of a eukaryotic cell (with regard to vertebrates) are subject to phosphorylation (20). Even if the possibility that the genome still has nondetermined nucleotides and that some gene models may be not correct is taken into account, one can predict that the current data do not cover all phosphoproteins from C. reinhardtii. Currently, there are about 19,000 predicted gene models (5), accounting for the same number of proteins in C. reinhardtii. If one assumes that about 30% of all proteins from C. reinhardtii are phosphorylated, a potential yield of 5,700 phosphoproteins would have to be considered. Taking into account the possibility that about 30% of the gene models might be wrong, as found previously (25), one would still expect 3,990 phosphoproteins. Thus, the yield of phosphoproteins identified in this study would cover about 8.2% of the total. Of course, the level of phosphoproteins in C. reinhardtii could be less than 30%, and indeed there is some indication, for example, that it differs in comparison to lymphoma cells (30). When the levels of proteins labeled with the antiphosphothreonine antibody (Cell Signaling Technologies) in the lymphoma cell line (20 µg of protein per lane [see Fig. 2 in reference 30]) are compared with those in C. reinhardtii (50 µg of protein per lane [Fig. 1]), one can detect more labeled phosphoproteins in the lymphoma cell line even though the amount of protein loaded per lane was less.
There are several reasons why the yield of phosphoproteins can be decreased: (i) very low abundance phosphopeptides are not enriched sufficiently by IMAC and thus have been missed; (ii) multiply phosphorylated peptides are more enriched by IMAC, as reported previously (20), and thus monophosphopeptides may be more easily missed; and (iii) phosphorylation of many proteins requires specific stimuli, and therefore any particular growth condition for the organism reflects only some of these stimuli. Thus, it was shown previously that calcium can repress the phosphorylation of an 85-kDa axonemal protein of C. reinhardtii at a concentration of
106 M, while the phosphorylation of a 95-kDa protein (designated b4) was increased at a similar concentration (29). These proteins seem not to be present in the current phosphoproteome. In this context, it could be of interest that a flagellar, putative plasma membrane calcium-transporting ATPase (C_700061) was identified among the phosphoproteins. Another example of stimulus-dependent phosphorylation is the light-harvesting complex II protein CP29. It contains four phosphorylated amino acids when the cells are exposed to so-called State 2 conditions, while in State 1 condition-exposed cells only two phosphorylated amino acids could be found (12). The State 1 to State 2 transition in the photosynthetic membranes involves the functional coupling of phosphorylated light-harvesting complexes of photosystem II to photosystem I. Thus, a set of different phosphoproteins will be present in the cell depending on the growth conditions (e.g., light, temperature, nutritional supplements, time of harvesting). Under our conditions, two phosphopeptides from CP29 (gene model C_10030) that cover the same amino acids were found. In the first phosphopeptide, one threonine (Thr6) was phosphorylated, while in the other, two threonines were phosphorylated within the same peptide. Only Thr6, not Thr10, of CP29 was reported to be phosphorylated under specific State conditions (12).
Our phosphoproteome data show that not only rather abundant phosphoproteins (e.g., metabolic enzymes) but also proteins that might be assumed to be present at low concentrations in the cell were detected. These include proteins from signaling pathways, such as several transporters, receptors, kinases, and phosphatases. Some of them, such as protein kinases, which can be subjected to positive and negative regulation by phosphorylation (17), could be expected to appear in the phosphoproteome. Further, putative transcription factors (e.g., plant homeodomain zinc fingers or SBP domain-containing proteins, which function as transcription factors in early flower development) and putative RNA-binding proteins containing, e.g., RRM domains, as well as DNA/RNA modifying proteins were detected. Also, proteins having domains that communicate protein-protein interactions (e.g., PAS, ankyrin repeats, TPR, Leu-rich repeats) were found.
Some of the proteins that were identified (e.g., CP29) had been already known to represent phosphoproteins of C. reinhardtii. Another example, the
dynein heavy chain (DHC) from the outer arm of the flagellum, was shown to have multiple sites of phosphorylation (16). From quantitative analysis it was concluded that DHC is phosphorylated at a minimum of six sites. In our study, we were able to detect four phosphopeptides from DHC with a total of eight phosphorylation sites. A further example is cytosolic glutamine synthetase, which was described previously as a phosphoprotein (27). One phosphopeptide with two phosphorylation sites was found in this case.
Comparative proteome analysis (Table 3) revealed that 32 proteins (with a total of 41 phosphopeptides) of the flagella and 4 proteins of the centriole are phosphorylated. A large number of phosphoproteins within the flagella could be expected, since several kinases and phosphatases have been shown to be located directly in the flagella (24, 35, 36), suggesting that dephosphorylation/phosphorylation plays a key role within the eukaryotic cilia. In addition to the above-mentioned DHC, which has multiple phosphorylation sites, a protein (gene model C_410060) from C. reinhardtii that is similar to the mouse hydrocephaly protein hydin (HY3) also contains numerous phosphorylation sites (eight sites situated in three phosphopeptides). HY3 is an example of a protein that is conserved in mammals and is involved in mammalian diseases. Thus, a mutation in mouse HY3 results in lethal communicating hydrocephalus with perinatal onset (1). Other phosphoproteins from the flagella (Table 3) are also related to diseases, such as, for example, predicted proteins from gene models C_290071 (similar to cold autoinflammatory syndrome 1 protein), C_420012 (polycystic kidney and hepatic disease 1-like protein), and C_590099 (a novel protein similar to human polycystic kidney disease 2-like protein 1). These findings open the possibility of specifically studying the mechanisms of posttranslational control for these proteins. Furthermore, some phosphoproteins from the flagella having a putative or unknown function can now be analyzed in a more focused way. However, the phosphoproteins identified within the flagella also reveal the limits of our analysis. For example, the phosphoprotein IC138, a WD repeat dynein intermediate chain that is hyperphosphorylated in paralyzed flagellar mutants (10, 15), was not found in the current phosphoproteome. Also, a significant number of radial spoke proteins of the flagella, which were identified as phosphoproteins by in vivo pulse labeling with 32P (26), are missing in the analyzed phosphoproteome.
Comparison of the proteome from the 70S ribosome (33, 34) with the phosphoproteome (Table 3) revealed only one hit, namely, RAP41. RAP41 represents one of two proteins that are specific for the complete 70S ribosome and could not be found in its 50S and 30S subunits (33). No components of the 50S and 30S subunits of the chloroplast ribosome could be found in the phosphoproteome. In contrast, proteins with significant homology to those from the cytosolic 80S ribosome (60S and 40S subunits) were found in the phosphoproteome.
Seven thioredoxin target proteins in C. reinhardtii (18) turned out to be phosphoproteins. This finding suggests that they are controlled in a complex manner by phosphorylation and redox regulation.
The present analysis has depicted a significant number of phosphoproteins of C. reinhardtii within the whole cell. It aimed to get a first insight into the variety of phosphoproteins within C. reinhardtii. Nevertheless, many phosphoproteins from this alga are still missing, based on prediction data for the phosphoprotein content within a eukaryotic cell (20). An efficient approach for the future to increase the yield of phosphoproteins may be to concentrate on the phosphoproteome of subcellular fractions. For example, flagella that bear many phosphoproteins could be isolated and used specifically for analysis. However, one must also point out that certain phosphoproteins may be ultimately missed by subjecting only subcellular fractions to analysis. Another approach may be to grow and harvest cells in the presence of different stimuli that might up-regulate phosphorylation of certain proteins. Thus, knowledge of stimulus-dependent phosphorylation would become available at the same time.
Our work was supported by grants from the Deutsche Forschungsgemeinschaft.
Supplemental material for this article may be found at http://ec.asm.org/. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»