|
Advertisement | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Molecular & Cellular Proteomics 7:1609-1638, 2008.
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|---|
|
|
|---|
Differential protein accumulation and activities between M and the BS chloroplasts of maize and sorghum (both NADP-malic enzyme (ME) type C4) have been studied using various low throughput techniques. This has shown that BS and M cells each accumulate a distinct set of C4 photosynthetic enzymes and has established the general pathway for C4 photosynthesis (5–7). A quantitative analysis of purified M and BS chloroplast stromal proteomes determined BS/M accumulation ratios for 125 proteins, covering a wide range of plastid functions and allowing integration of information from past studies (8).
M thylakoids have a complete linear electron transport chain, containing the Photosystem II (PSII), cytochrome b6f (cytb6f), and PSI complexes, similar to C3 plants. In contrast, fully differentiated BS thylakoids have little functional PSII and normal or increased PSI levels, whereas cytb6f and ATP synthase complexes are quite evenly distributed between BS and M thylakoids. As a consequence, BS thylakoids mostly carry out cyclic electron flow and have low rates of linear electron flow (9–11). After several conflicting reports on protein levels of PSII and light-harvesting complex (LHC) II (11–13), a subsequent study partially clarified both activity and protein accumulation levels of PSI and PSII complexes in BS and M thylakoids (14). However, a more comprehensive overview of accumulation levels and assembly state of the thylakoid-bound photosynthetic apparatus remains to be determined. Thylakoids also contain NAD(P)H dehydrogenase (NDH) complexes of which several chloroplast-encoded and nucleus-encoded subunits were identified in maize (15) with a clear enrichment of NDH complexes in the BS membranes (16–18). However, a better understanding of the composition and assembly state of the NDH complex, as well as interactions with other electron transport components, will be important to understand its contribution to cyclic electron flow.
Proteome analysis of thylakoid membranes (19–22) and associated lipoprotein particles (plastoglobules) (23, 24) as well as the plastid chromosome (the nucleoid) (25) from Arabidopsis chloroplasts also identified many proteins with non-photosynthetic functions (26). In addition forward and reverse genetics studies in Arabidopsis (27, 28) as well as maize (29) identified several nucleus-encoded thylakoid proteins affecting biogenesis and chloroplast functions (30). Finally a number of stromal metabolic enzymes may functionally interact with the thylakoid, e.g. to transfer reducing equivalents to metabolic pathways, such as GcpE/IspG in the methylerythritol phosphate pathways (31). Thus, it is highly likely that numerous thylakoid (-associated) proteins with non-photosynthetic functions accumulate to different levels in the two cell types and may contribute to BS/M chloroplast differentiation; their quantification will help to better understand cell specific differentiation.
Chloroplasts are surrounded by a double envelope membrane containing the protein import apparatus (32, 33) and transporters for exchange of ions and metabolites and also participate in metabolic pathways and plastid-nucleus signaling (34–37). Hundreds of proteins have been experimentally identified in the inner and/or outer envelope proteome of Arabidopsis (38–40). However, no systematic comparison of envelope proteomes of differentiated BS and M maize chloroplasts has been carried out, but major differences can be expected in e.g. translocators of carbohydrates (34, 41).
Gel-based and, in particular, MS-based techniques to determine quantitative differences between proteomes have greatly improved in recent years (42, 43). Some of these techniques (two-dimensional (2D) IEF gels, ICAT, and label-free MS-based quantification using MS ion chromatograms) were applied in our previous comparative analysis of differentiated BS and M stromal proteomes (8). However, quantitative comparison of M and BS chloroplast membrane proteomes represents a significant challenge because of their hydrophobicity (19, 21, 44, 45). Moreover the presence of about 100 abundant photosynthetic proteins in thylakoids representing some 98% of the protein mass provides an additional challenge to quantify the low abundance proteins. In addition, although most maize genes are probably represented in the collection of maize EST and unigenes, no complete and assembled maize genome is yet available, thus making protein identification and quantification more difficult.
The prime objectives of this study were to determine maize BS and M cell-specific differences in (i) assembly state and composition of the four photosynthetic complexes and the NDH complex, (ii) auxiliary functions of the thylakoid proteome, and (iii) metabolic and transport functions of M and BS thylakoids and envelopes. Therefore, we carried out a quantitative proteome analysis of differentiated BS and M membranes using techniques compatible with membrane proteomes and also taking advantage of a new, fast, and highly accurate mass spectrometer, the LTQ-Orbitrap (46–48). This clarified the organization of photosynthetic machineries of BS and M thylakoids, revealed large NDH complexes and associated proteins with unknown function, and determined functional differentiation of transporters and biosynthetic pathways. These results complement our recent comparative analysis of soluble stromal BS and M chloroplast proteomes (8) and will provide new entry points for future studies to unravel cell-specific chloroplast differentiation in maize. Proteomics data and functional annotation are available via the Plant Proteome Database.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
2D Blue Native Gel Analysis and MS-based Protein Identification—
After purification of BS and M membranes, identical amounts of BS and M proteins were solubilized in β-n-dodecyl maltoside (β-DM) (Sigma) and separated on blue native gels as described previously (50). The gel lanes were then cut out, and proteins were denatured, reduced, and alkylated (8) and separated by Tricine-SDS-PAGE (51) as described previously. The resulting focused 2D gels were stained with Coomassie Brilliant Blue R-250 (USB Corp.) and scanned (Epson Perfection 4490). Seven independent M and BS chloroplast preparations were resolved on seven pairs of 2D BN-PAGE gels. Image analysis was performed using Phoretix software (Nonlinear Dynamics Ltd.). About 180 spots/gel were detected with spot volumes spanning 4 orders of magnitude. Spot matching was manually verified with a large time investment in spot matching and verification. Ambiguous matches were resolved by protein identification from MS. After background removal and normalization of each spot "volume" to the total gel spot volume, virtual average M and BS gels were created where each spot volume represented the average of volumes of matched spots in M or in BS gels. We included on the average gels only those matched spots that were present at least three times on M or BS gels. BS/M accumulation ratios were obtained for 85 spots. Spot quantification data are shown in supplemental Table 1.
Stained gel spots were excised manually or using a ProPic robot (Genomic Solutions, Ann Arbor, MI). The spots were washed, digested with modified trypsin (Promega), and extracted using a ProGest robot (Genomic Solutions). The extracted peptides were dried and resuspended in 15 µl of 5% formic acid (FA). Protein identification was performed by peptide mass fingerprinting using MALDI-TOF MS in reflectron mode (Voyager DE-STR from Perseptive Biosystems) and on-line LC-ESI-MS/MS (Q-TOF, Micromass). All protein identification data are shown in supplemental Table 1.
The mass spectra were obtained automatically by MALDI-TOF MS in reflectron mode followed by automatic internal calibration using tryptic peptides from autodigestion. Peak lists (.mgf files) from the MALDI data were generated using MoverZ software m/z (freeware edition, Proteometrics, Inc.) using a minimum signal to noise ratio of 3.0 and peak resolution of 5000. A large portion of the spectral annotations (in particular assignments of monoisotopic masses) were manually verified. The resulting peptide mass lists were used to search the Maize EST assembly from TIGR (ZmGI v16.0, 56,364 entries) by Mascot (v2.2) in automated mode using the following search parameter criteria: significant protein Mowse score p < 0.05, no missed cleavages allowed, variable methionine oxidation, fixed carbamidomethylation of cysteines, and minimum mass accuracy of 100 ppm. These search result pages were extracted and further filtered to pass the following criteria for positive identification: (i) at least five or more matching peptides with mass error distribution clustered within 25 ppm and (ii) at least 15% sequence coverage. Only peptides without missed cleavages (by trypsin) were considered with methionine oxidation as a variable modification and carbamidomethylation as a fixed modification. In exceptional cases (i.e. proteins less than 20 kDa and matching gel coordinates) four matching peptides were considered as positive identification. The criteria for positive protein identification used here result in less than 1% false positives as determined by searching the data from these gel spots using the target-decoy database consisting of the ZmGI v16.0 database (56,364 entries) concatenated with a decoy where all the sequences were shuffled. FPR was calculated as follows: 2 x decoy_hits/total_hits. PDF files of annotated spectra with assigned peptide masses are provided in the supplemental data. Data for protein identification by MALDI-TOF peptide mass fingerprinting are listed in supplemental Table 1, supplemental Fig. 1, and associated additional files showing spectral annotations.
The Q-TOF instrument was operated in positive ion mode with a sample cone voltage of 35 kV, capillary voltage of 3.3 kV, and source temperature of 90 °C. The samples were run in data-dependent mode where each full MS scan was followed by three consecutive MS/MS scans. The MS survey scans (m/z 350–1550) had a scan time of 1 s and an interscan time of 0.08 s. MS/MS spectra were automatically acquired when the peak intensity rose above a threshold of 10 counts·s–1. Normalized collision energies for peptide fragmentation were set using the charge state recognition files for 1+, 2+, 3+, and 4+ peptide ions provided by MassLynx (Waters). For tandem MS acquisition we used a scan range from m/z 50 to 1550 with a scan time of 1.92 s, an interscan time of 0.08 s, and a dynamic exclusion window of 240 s. Argon was used as the collision gas. Peak lists (.mgf files) from the Q-TOF data were generated using Mascot Distiller (v2.0).
The peak lists from Q-TOF were searched against the Maize EST assembly from TIGR (ZmGI v16.0, 56,364 entries) by Mascot (v2.2) with p < 0.05, with a maximum mass error of 0.8 and 2 Da for product and precursor ions, respectively, and with a minimum ion score of 20. Mascot results were subsequently filtered to increase the stringency for positive identification by MS/MS as follows. If only one or two matching peptides were found, clear partial Y-ion series and partial complementary B-ion series needed to be present (as determined by manual inspection) with a minimum peptide ion score of 31 for one-peptide identifications and minimum peptide ion scores of 22 for two-peptide identifications. If three or more sequence tags were found, manual inspection was not a requirement, but the protein Mowse score must be over 70 with only peptides with ion scores of 20 or more contributing. These criteria resulted in a protein false positive rate below 1% as determined from target-decoy database searches using the ZmGI v16.0 database (56,364 entries) concatenated with a decoy where all the sequences were shuffled. Protein FPR was calculated as follows: 2 x decoy_hits/total_hits. Protein identification data are displayed in supplemental Table 1, and additional Mascot protein report pages for single peptide identifications are provided as supplemental files.
Seven independent biological replicates (pairs of BS and M membranes) were used for this 2D BN analysis. The minimum requirement for spot quantification was spot presence in at least three biological replicates for either M or BS membranes. The average BS/M ratio across the replicates and their coefficients of variation (cvs) for each quantified accession, normalized average spot volume with their cvs for each BS and M sample, and the histograms with spot volume for each of the seven replicates are listed in supplemental Table 1. We report only one representative accession for each group in Table I; redundant accessions are reported in supplemental Table 1.
|
The recalibrated peak lists were searched against the ZmGI v16.0 database (56,364 entries) concatenated with a decoy where all the sequences were shuffled. Each of the peak lists was searched using Mascot v2.2 (maximum p value of 0.01) for full tryptic peptides using a precursor ion tolerance window set at ±6 ppm, variable methionine oxidation, fixed cysteine carbamidomethylation, and a minimal ion score threshold of 42 that yielded a peptide false discovery rate below 1%. Peptide FPR was calculated as follows: 2 x decoy_hits/total_hits. The false protein identification rate of protein identified with two or more peptides was 0. To reduce the false protein identification rate of proteins identified by one peptide, the Mascot search results were further filtered as follows. The ion score threshold was increased to 45, and mass accuracy on the precursor ion was required to be within ±3 ppm. All filtered results were uploaded into the PPDB. Protein identification data are also provided in supplemental Table 2.
To assess the reproducibility between technical replicates, we used the G-test of independence (53). Recently this test has been used to assess significance of expression changes in label-free comparative proteomics experiments for eukaryotes and prokaryotes (54, 55). We applied the G-test to three technical replicates of each of the MS and BS preparations. In both of the cases, about 2% of all proteins were determined to be significantly changed between the three replicates, showing that the technical replicates were highly reproducible.
To determine protein accumulation ratios between BS and M, we calculated two parameters for quantification: (i) normalized number of matched spectra for each protein across all 27 gel slices per replicate (
Counts) and (ii) normalized total Mowse protein score for each protein across all 27 gel slices per replicate (
Mowse). Only those proteins for which the sum of
Counts across the three replicates was higher than 10 in either cell type were used for quantification. G-tests were used to help determine this cutoff as explained above. The BS/M protein ratios were then calculated based on averaged (across the three replicates)
Mowse score or
Counts for each protein per cell type (Table I and supplemental Table 7). We report only one representative accession for each group in Table I; redundant accessions and all peptide sequences are reported in supplemental Table 2.
iTRAQ Analysis—
Two independent biological pairwise replicates of BS and M cells were isolated, and 30 µg protein from each of the four samples were precipitated and digested as described previously (23). The peptides were dissolved with 30 µl of 0.5 M triethylammonium bicarbonate and mixed with the four iTRAQ reagents (Applied Biosystems). The samples BS1, M1, BS2, and M2 were labeled for 90 min with the four iTRAQ reagents 114, 115, 116, and 117 and then pooled. The sample was then diluted 10-fold to a final concentration of 25% ACN and 1% FA, and half of the final volume was separated using strong cation exchange chromatography on a PolySULFOETHYL ATM column (200 x 2.1 mm, 5 µm, 300 Å) from PolyLC Inc. connected to an Agilent 1100 HPLC system. Solution A was 25% ACN, and solution B was 25% ACN, 0.5 M ammonium formate (pH 3.0, set with FA). 500 µl of the sample were injected and eluted at a flow rate of 200 µl/min as follows: 3% B, 0–10 min; 3–10% B, 10–15 min; 10–60% B, 15–25 min; 60–100% B, 25–35 min; 100% B, 35–45 min; 100–3% B, 45–47 min; and 3% B, 47–55 min. The fractions were collected in microtiter plates, lyophilized, resuspended with 15 µl of 5% FA, and analyzed by MALDI-TOF MS. The 10 fractions containing the majority of the labeled peptides were selected for analysis by Q-TOF MS. The data were processed with Mascot Distiller v2.1 (Matrix Science) resulting in .mgf files, and the areas of the reporter ions were extracted from the .mgf files and searched against ZmGI v16.0 (56,364 entries) using Mascot (v2.2) with p < 0.05 using the following criteria: one missed cleavage allowed, variable methionine oxidation, cysteine carbamidomethylation and iTRAQ label as fixed modifications, and a minimum peptide ion score threshold of 20.
Search results were extracted from the XLM pages and further filtered for overlap with other accessions using an in-house program written by Dr. Qi Sun. All peptides where the difference in score between the highest and the second highest interpretations of the same query was less than 10 points were manually checked, and if it was not possible to distinguish which was the best interpretation, both were removed. Areas from spectra derived from the same sequence (e.g. repetitive fragmentation of the same precursor, different charge states, or variable modifications) were pooled. The areas were normalized by adding the areas from all peptides with all four reporter ions. Spectra lacking one or two reporter ions were kept if both ions from the same biological replicate (e.g. BS1 and M1 or BS2 and M2) were present and used to calculate the protein average. Peptides were grouped by accession. In the case of seven proteins quantified with six to 20 peptides, one or two peptides showing a ratio clearly different from the average of the majority of the peptides were removed from the calculation of the average protein accumulation.
In total, 183 accessions were identified of which 116 were assigned an average BS/M accumulation ratio based on quantification of at least two peptides, and the corresponding S.D. and cv are reported (Table I and supplemental Table 7). Proteins exclusively identified in M or BS proteomes were removed from the quantification analysis. The average cv across the quantified peptides per protein for all protein ratios was 0.26. In a number of cases, the MS data matched two or more closely related proteins, particular in the LHC family and ribosomal proteins. We report only one representative accession for each group in Table I; redundant accessions and all peptide sequences used for iTRAQ quantification are reported in supplemental Table 5. All sequences used for iTRAQ quantification are listed in supplemental Table 6.
The Plant Proteome Database and Functional Assignment of Identified Proteins—
The Mascot output files were automatically processed by in-house software, 3 and a number of output parameters were uploaded into the PPDB. Because many of the maize ZmGI accessions lack functional annotation, we functionally annotated all identified ZmGI accessions using a combination of best BLAST hits in the predicted rice proteome (Oryza sativa Gene Index) and the predicted Arabidopsis thaliana proteome, ATHv6 (from The Arabidopsis Information Resource). Pairwise BLAST search results between ATHv6, O. sativa Gene Index v4, and ZmGI v16.0 are available via PPDB. Each identified protein was assigned to a molecular function using the hierarchical, non-redundant classification system developed for MapMan (56), adjusted after manual verification and information from the literature, and incorporated into the PPDB. In some cases ZmGI accessions numbers were identified as a set of accessions sharing entirely the set of identified peptides; this was essentially due to a lack of good assembly of ZmGI accessions within small gene families (e.g. LHCs). In these situations we cite only one accession number in quantification tables; the subset accessions are listed in tables containing the corresponding identification information (supplemental Tables 1, 2, and 5).
| RESULTS |
|---|
|
|
|---|
Overview of the Comparative Proteome Analysis—
To overcome the challenges of comparative membrane proteomics in maize, we used a combination of three comparative proteomics methods that are compatible with membrane proteins and that are based on different principles for quantification as follows: (i) quantification based on spot image analysis of 2D gels using blue native gels (50) as the first dimension and SDS-PAGE as the second dimension (further referred to as "2D BN"), (ii) MS-based quantification using label-free analysis by nano-reverse phase LC-LTQ-Orbitrap MS (further referred to as the "label-free" method), and (iii) MS-based quantification using differential peptide labeling with an amine-reactive isobaric tagging reagent (iTRAQ).
The mass spectral data were searched against the assembled transcript sequences of maize TIGR ZmGI v16.0 supplemented with complete chloroplast and mitochondria genome sequences obtained from the National Center for Biotechnology Information (NCBI). The collective data set was functionally annotated based on manual annotation for about 400 accessions and BLAST alignments to rice and Arabidopsis and the MapMan functional classification system (56). Mass spectral information extracted from the Mascot search pages, functional annotation, and the highest scoring BLAST results against Arabidopsis and rice are available through the PPDB. Matched peptides projected on translated reading frames can also be viewed via PPDB. A number of ZmGI unigenes were identified that appeared to be assembled incorrectly; this is indicated in the annotations.
2D BN Analysis—
2D BN has been used to analyze protein composition and oligomeric states of chloroplast envelope and thylakoid membranes (57, 58). Following similar procedures, the purified BS and M membranes were solubilized in β-DM and run on BN gels followed by separation of each gel lane by SDS-PAGE. Seven independent biological replicates (pairs of BS and M membranes) were used. The minimum requirement for spot quantification was spot presence in at least three biological replicates for either M or BS membranes. Representative 2D BN gels of M and BS chloroplast membranes with annotation of the major complexes and proteins are shown in Fig. 1A. The differential BS/M accumulation ratio was determined for 85 spots (supplemental Table 1), and examples of quantifications are shown as bar diagrams (Fig. 1B).
|
270 kDa, and a partial core at
200 kDa) were observed. Accumulation of the PSII monomer and a partially assembled PSII core was respectively 3- and 5-fold reduced in BS membranes, whereas the high molecular weight PSII complexes observed in M thylakoid were below detection in BS membranes. Interestingly the accumulation of LHCII trimers (
140-kDa complex X) was only 1.6-fold reduced in the BS membranes suggesting that a subset of LHCII proteins do function in the BS consistent with previous findings (13) (Fig. 1B). Furthermore BS thylakoids were characterized by a 40% higher accumulation of PSI and its antennae LHCI subunits in particular within the high molecular mass forms of PSI-LHCI oligomeric assemblies (>700 kDa). A 30% increased BS accumulation of the ATP synthase complex (300 kDa) was also observed. NDH complex (complex V) was detected at
550–600 kDa and was 400% enriched in the BS membranes (Fig. 1B). Interestingly thylakoid lumenal protein isomerase TLP20 was identified with complex V, suggesting a possible association with the NDH complex. Finally an abundant and novel complex of
800 kDa (complex III) was found in BS membranes; it includes two homologous proteins (TC299048 and TC299050) of unknown function that co-migrate with NdhF on 2D BN (spot 200; Fig. 1B) and show a strong co-migration pattern with the NdhD and -F subunits (see below for 1D BN profile analysis). FNR1 showed a strong M accumulation (BS/M = 0.2), the FtsH thylakoid protease accumulated with a BS/M ratio of 0.9, and a strong M expression was determined for the Tic110 subunit of the inner envelope translocon.
1D BN-PAGE and Label-free Quantitative Analysis by LTQ-Orbitrap MS—
The 2D BN shows distinct differences between BS and M membranes, but proteins of lower abundance are not visible on the gels and thus cannot be quantified. Moreover, spots containing more than one protein do not provide reliable quantification. Therefore, we separated the BS and M proteomes by BN-PAGE only and then directly analyzed and quantified the proteins by MS using the fast and accurate LTQ-Orbitrap. Each BN-PAGE lane was cut in 27 bands, proteins were cleaved by in-gel digestion with trypsin, and extracted peptides were analyzed by MS/MS (Fig. 2, A and B). The on-line chromatography and MS acquisition were optimized for label-free quantification (59). MS data were analyzed against ZmGI v16.0 followed by additional filtering (see "Experimental Procedures"). Each sample was analyzed in triplicate to reduce variation due to chromatography and MS acquisition. The chromatography was very reproducible (not shown), which was also reflected by similar total numbers of matched spectra (counts) and total Mowse scores between the replicate runs (Fig. 2B).
|
Counts) and (ii) normalized total Mowse protein score for each protein across all 27 gel slices per replicate (
Mowse). For information on the average cv, see "Experimental Procedures." To increase the "robustness" of the data set for quantification, we discarded those proteins for which the sum of
Counts across the three replicates was lower than 10 in either cell type (corresponding to a minimum
Mowse of
420). G-tests were used to help determine this cutoff (see "Experimental Procedures"). This reduced the number of proteins only found in BS or M cells to 11 and 51, respectively; and 289 proteins were found in both BS and M. The BS/M protein ratios were then calculated based on averaged (across the three replicates)
Mowse score or
Counts for each protein per cell type (Table I and supplemental Table 7). Cross-correlation of the BS/M protein ratios calculated by the two parameters showed a strong linear correlation (r2 = 0.97), indicating that the two parameters yielded comparable results (Fig. 3A).
|
The use of the 1D BN as a protein separation method prior to the MS analysis also provided information about oligomeric state and was complementary to the 2D BN analysis. We calculated the normalized
Mowse score for each protein per gel slice (supplemental Table 3) and created profiles of their oligomeric assembly state in BS and M membranes. Examples of profiles for subunits of photosynthetic complexes are shown in Fig. 4; profiles corresponded well to these observed from gel image analysis on 2D BN. These native mass profiles will be used further below to characterize BS/M differentiation.
|
In total, 183 accessions were identified of which 116 were assigned an average BS/M accumulation ratio based on quantification of at least two peptides, and the corresponding S.D. and cv are reported (Table I and supplemental Table 7). Proteins exclusively identified in M or BS proteomes were removed from the quantification analysis. In a number of cases, the MS data matched two or more closely related proteins, particularly in the LHC family and ribosomal proteins. We report only one representative accession for each group of homologues in Table I; redundant accessions and all peptide sequences used for iTRAQ quantification are reported in supplemental Table 6.
The Identified and Quantified Proteomes and Consistency—
In total 610 proteins were identified with the highest number (543 proteins) obtained by the 1D BN analysis. Proteins were functionally classified using the MapMan functional classification system developed for Arabidopsis (56) (supplemental Table 7). 398 proteins passed the minimal criteria for quantification (described above) by at least one of the three quantification methods and were assigned a BS/M ratio (summarized in Table I). An additional 192 proteins that did not pass the quantification quality filters are only listed in supplemental Table 7. Generally the three methods gave similar BS/M ratios with an average cv across the three methods of 0.28, but we did note that the BS/M ratios determined by iTRAQ appeared to be generally closer to unity than the ratios determined by the other two methods (Table I). This suggests that the label-free quantification (also named spectral counting) has a better dynamic resolution than iTRAQ. Furthermore the quantification from the 2D BN gel image analysis of the more abundant proteins showed similar or slightly more extreme BS/M ratios than the label-free quantification.
A direct comparison of the BS/M ratios of 23 proteins with quantifications from our previous stromal analysis (8) showed good consistency and also confirmed that the BS/M cross-contamination was well below 15% (Fig. 3B and supplemental Table 8). Examples of abundant markers of the C4 carbon fixation pathway (average ratios across three methods) were: PPDK, BS/M = 0.4; malate dehydrogenase (MDH), BS/M = 0.24; glyceraldehyde dehydrogenase subunit B (GAPB), BS/M = 0.39; Rubisco large subunit (RBCL), BS/M = 2.71; Rubisco small subunit (RBCS), BS/M = 7.69; and Rubisco activase (RCA), BS/M = 5.42. Furthermore only a handful of non-chloroplast proteins were identified, and none passed the criteria for quantification (see supplemental Table 7), indicating that the chloroplasts were of high purity.
Differential Accumulation of Photosystem II Subunits and Subcomplexes—
PSII can be viewed as the assembly of the PSII reaction center surrounded by the core subunits, the oxygen-evolving complex (OEC16,23,33), and the major and minor LHCII antennae proteins. Fifteen proteins of the reaction center and core are encoded by the plastid genome, whereas the other subunits are encoded by the nuclear genome; in the context of BS/M differentiation this might be important because transcriptional and translational controls of these two groups are different.
The PSII Core and OEC—
The chloroplast-encoded PSII core subunits D1 (NP_043004), CP47 (NP_043049; TC283413; TC225511), CP43 (TC241707), and cytb559
(NP_043041; TC242401) showed average BS/M ratios between 0.25 and 0.72 (Table I). The nucleus-encoded core subunit PsbR (TC279588) showed a BS/M ratio of 0.41 (Table I), whereas nucleus-encoded subunits PsbTn-1 (TC305359), PsbW (TC282454), and PsbH phosphoprotein (NP_043052) were only identified in M thylakoids (supplemental Table 7). OEC16, -23, and -33 proteins represented by six ZmGI accessions showed a BS/M accumulation ratio between 0.22 and 0.54, which was on average lower than the rest of the PSII core (Table I). In contrast, two distantly related OEC16-like proteins with unknown function (AI001247 and TC307806) each showed much higher BS/M accumulation ratios of 2.39 and 2.81, respectively, strongly suggesting that they do not belong to PSII. In fact AI001247 showed a weak similarity to the NdhL subunit of Synechocystis elongatus PCC6803, and both proteins TC307806 and AI001247 co-migrated with the NDH complexes (see below). We also identified three distant homologues of OEC23, namely TC307123, TC295379, and TC289659, with BS/M ratios of 1.17, 0.61, and 2.82. Interestingly Arabidopsis homologues named PPL1 (AT3G553330) and PPL2 (AT2G39470) for TC295379 and TC289659, respectively, were recently characterized in Arabidopsis (62). PPL1 was suggested to be required for PSII stability and repair, whereas PPL2 was suggested to be required for accumulation of the NDH complex. This does fit very well with the BS/M ratios for the two maize proteins. Moreover bona fide maize OEC subunits were identified essentially in the low molecular weight fraction in 1D BN analysis probably because they were released from the membrane upon detergent treatment. In contrast, AI001247 (distant OEC16 and putative NdhL), TC289659 (PPL2 homologue and putative NDH subunit), and TC307806 (distant OEC23) were all found in high molecular mass complexes between 250 and 750 kDa (supplemental Table 3). The PPL2 and "distant OEC16" homologues showed a clear co-migration pattern and similar BS/M expression ratios with different oligomeric assemblies of the NDH complex, suggesting that they might be functionally related with the NDH complex.
LHCII Antennae—
The major and minor LHCII antennae were represented by 18 different accessions with an average BS/M ratio of 0.68 (cv = 0.33) (Table I). Six ZmGI accessions matched best to LHCII1.5; within this cluster two accessions appeared to have higher BS/M ratios than the other LHCII proteins, which may indicate that these proteins function in both compartments. Indeed although Schuster et al. (11) observed the absence of LHCII polypeptides in BS chloroplasts in 2–3-week-old maize seedlings, others determined that LHCII gene family members showed different cell-specific expression patterns (measured at the mRNA level) in 14-day-old seedlings (63). These results were supported by the observation of different ultrastructures of the LHCII in freeze fracture experiments between maize M and BS thylakoids (13), suggesting that complexes of LHCII proteins in the M and BS are of different composition.
Light Stress Proteins with Chlorophyll Binding Domains—
PsbS (TC300425) was identified with a very high
Mowse score (20,612) in the label-free experiment and a BS/M ratio of 0.5, similar to the PSII core, underscoring the involvement of PsbS in dissipation of the excess energy from PSII (64). PsbS 1D BN profiles show large peaks (at
650,
350, and
200 kDa) encompassing PSII assemblies (supplemental Table 3). However, these interactions are relatively weak because most PsbS was found in low molecular weight fraction 25 consistent with observations that PsbS links LHCII antennae complexes with the PSII core rather than being a stable core subunit (65). We also identified and quantified maize homologues of small light stress proteins with chlorophyll binding domains, Ohp1 (TC310009), Ohp2 (BM268267 and TC290705), Lil3.2 (TC282388), and Lil4/Sep1 (TC313218) that was only found in M thylakoids. The two Ohp2 accessions were strongly enriched in M thylakoids (BS/M ratio, 0.18 and 0.35), whereas the Ohp1 BS/M ratio was 0.69. 1D BN protein profiles for Ohp2 (BM268267) showed a co-migration with the PSI and PSII dimers suggesting possible interaction and its contribution in PSII and/or PSI stress response (supplemental Table 3). In contrast, Lil3.2 was equally distributed over the two cell types (BS/M = 1), suggesting a role in PSI stress response.
PSI, LHCI, and FNR—
The PSI and LHCI proteins, represented by 25 and 13 quantified accessions, respectively, showed preferential BS expression with average BS/M ratios of 1.57 (cv = 0.28) and 1.72 (cv = 0.43), respectively (Table I and supplemental Table 7). These observations are consistent with previously suggested higher accumulation of PSI in BS chloroplasts of C4 plants based on an immunoblot against the PsaA protein (12). The exception was LHCI-680B (BQ619330), which showed higher M accumulation. PSI-associated PsaP (TC306448) protein was expressed equally in both compartments. Strong differential BS/M expression was observed by all thee methods for FNR1 in agreement with our previous observation for the soluble stromal proteomes (8). FNR1, represented by two accessions, was either strongly expressed in the M (TC310667) or equally distributed over both compartments (TC281599).
Cytb6f and ATP Synthase Complexes—
The three quantified subunits of cytb6f (cytf, Rieske Fe-S, and subunit IV), represented by six accessions, showed very little BS/M differentiation with an average BS/M ratio of 1.03 (cv = 0.26). This is logical because cytb6f is needed for both linear and cyclic electron flow. Some variability was observed within two ZmGI accessions corresponding to the Rieske subunit (TC286498 and TC286511), suggesting that two different isoforms are differentially expressed between M and BS (Table I and supplemental Table 7). All nine subunits of the ATP synthase complex (represented by 15 ZmGI accessions) were quantified with an average BS/M ratio of 1.33 (cv = 0.37), indicating a fairly equal distribution of the ATP synthase over BS and M thylakoids. The somewhat higher cv is most likely due to redundant and incorrect ZmGI unigene assemblies (Table I).
The NDH Complex, New Interaction Partners, and IMMUTANS—
Thylakoid NDH complex is involved in one of the two pathways for cyclic electron flow around PSI and also in chlororespiration. In C3 plants, the NDH complex is of much lower abundance than the four major photosynthetic complexes (estimated at 0.2% of total thylakoid protein (66)). The plastid genome encodes for 11 putative NDH subunits (A–K), and several have been identified in both C3 (67) and C4 (15, 18, 68) plants. Recently five nucleus-encoded Ndh subunits (M, N, and O and CRR3 and -7) were identified in Arabidopsis (69–71).
Here we quantified seven chloroplast-encoded (A, D, F, H, J, I, and K) and three nucleus-encoded maize NDH subunits (M, N, and O), all showing strong BS expression with an average BS/M ratio of 3.0 (cv = 0.32) (Table I). In addition we identified NdhB (CF039487) with low scores in complex V (spot 213; data not shown) in 2D BN and quantified the maize homologue of PPL2 (BS/M = 2.8) as mentioned above. A misassembled accession (TC294644) containing fragments of NdhA and -G and an accession (TC297967) containing fragment sequences similar to NDH-J and -K annotated here as NdhJorK were both quantified with similar BS/M ratios (2.67 and 2.70). Preferential accumulation of the NDH complex was described in BS chloroplasts of sorghum (68) and maize with BS/M ratios between 2.5 and 3 (16, 17), which is consistent with our observations. We did not detect any of the CRR proteins, involved in NDH biogenesis or NDH subunits (67, 70, 71), but we note that the maize homologue of soluble CRR1 (TC242930) was found in our previous stromal analysis (8).
Analysis of the 1D BN label-free profiles showed that all identified NDH subunits assemble in high molecular weight complexes with the exception of NdhA, which was only found as a monomer (band 24). Detachment of NdhA was also observed previously using either sucrose gradient fractionation (15) or affinity purification (69), suggesting that it is located at the periphery of the complex. Based on peaks detected for each of the NDH subunits in the 1D BN profiles (and by 2D BN analysis for the
800- and
550-kDa complexes), we conclude that the membrane-embedded and -connecting NDH subunits co-migrate in two assemblies at
1000 and
550 kDa (Figs. 4B and 5 and supplemental Table 4). Although we did not identify all membrane-embedded NDH subunits (because of their small size and hydrophobicity), it is likely that these
550- and
1000-kDa complexes represent a monomer and a dimer, respectively, of a fully assembled NDH complex (Fig. 5, A and C). The 1D BN profiles showed additional NDH assemblies with a highly abundant and well defined peak at
800 kDa (containing NdhD and -F) and at
650 kDa (containing NdhJorK, -M, -N, and -O) and a less abundant complex at
320 kDa (NdhD and -F), at
220 kDa (NdhH, -J, -JorK, -M, and -N), and at
170 kDa (NdhH, -I, -J, -JorK, -M, -N, and -O). These are likely the result of destabilization due to the detergent, and they do suggest a modular organization of the NDH complex. Consistently the NDH complex in tobacco and maize was purified as an assembly of 550 kDa (15, 72). Using 2D BN-PAGE and four NDH-directed antibodies, the Ndh complex in maize was found at 550 and
1000 kDa, and subcomplexes were found at 300 and 250 kDa (18).
|
550 kDa (Figs. 3 and 4 and supplemental Table 4). The 2D BN gel analysis showed TLP20 co-migrating with the
550-kDa NDH complex, whereas FKBP was found in the same spot (spot 192) as NdhJ and PPL2.
Interestingly a second set of four unknown proteins (TC291961, TC299050, TC299048, and TC305219) and a ferredoxin (Fd) domain protein (TC307918) showed 1D BN profiles very similar to those of NdhD and -F with peaks at
1000,
800, and
320 kDa with a predominant and well defined peak at
800 kDa (band 4) (Fig. 4, B and C, and supplemental Tables 3,4). Surprisingly abundant, this
800-kDa complex was also observed on 2D BN gels (Fig. 1, complex III). An additional ZmGI accession (TC234236) similar to TC299050 was identified on 2D BN at
1000 kDa (Fig. 1A, spot 207) also preferentially expressed in the BS (Table I). TC299050 and TC299048 are similar to Arabidopsis proteins AT1G15980 and AT1G64770. These Arabidopsis proteins were identified in the chloroplast (73) and thylakoids (59) but have no known function and no clear predicted functional domains. The strong clustering of BS/M ratios between these proteins and the NDH subunits and their co-migration in three different native assemblies suggest that they constitute interacting partners of the NDH complex (Fig. 5, B and C).
A recent study in Arabidopsis proposed that the PIFI protein (AT3G15840) (postillumination fluorescence increase protein) was a new and essential component of Ndh (74). We identified a maize homologue of PIFI (TC279329) with a BS/M ratio of 1.88; however, it was not associated with NDH complex assemblies. A maize homologue (TC307307) of the Arabidopsis IMMUTANS (IM) protein was identified only once in this study and only in BS membranes. IM is a chloroplast alternative oxidase and is likely to be the elusive terminal oxidase in chlororespiration. IM appears to be a versatile electron sink especially early in chloroplast development, and its inactivation in Arabidopsis leads to leaf variegation (75). Our data suggest surprisingly low accumulation levels as compared with NDH.
Regulation of Light Harvesting Capacity—
Reversible phosphorylation of light-harvesting antennae LHCII and subunits of PSII regulate the state transitions (76, 77) and the PSII repair cycle (57), respectively. State transitions involve reversible interactions of LHCII with PSII and PSI, thus changing the balance of the light harvesting capacity of the two photosystems, and involve lateral movement of LHCII antennae proteins. Two Arabidopsis thylakoid serine-threonine kinases, STN8 and STN7 (78–80), have been identified. Whereas STT7 and STN7 phosphorylate LHCII and are required for state transitions, STN8 specifically phosphorylates the N-terminal threonine residues of PSII core subunits D1 and D2 and CP43 proteins and Thr-4 in the PsbH protein (80).
We identified maize homologues of both STN7 (TC305194) and STN8 (TC311729) kinases. Both proteins were only identified in the M proteome suggesting enrichment in the M thylakoid (we note that STN8 was identified in all three replicates with good scores; STN7 was only identified in one replicate; see supplemental Table 7). In addition, 1D BN profiles indicated that STN7 migrates as a dimer (
110 kDa) or associated as a monomer with its target, the LHCII trimer (supplemental Table 3). STN8 migrated in two complexes (
180–200 and
100 kDa).
A soluble thylakoid phosphoprotein TSP9 was shown to interact with PSII antennae as well as the PSI and PSII cores, and it was suggested to play a role in regulation of light harvesting capacity (81). We identified a maize homologue of TSP9 (TC312266) that was more highly expressed in M thylakoids (BS/M = 0.45), and its 1D BN profile showed peaks in bands 10 (
500 kDa) and 13 (
300 kDa) with the majority found as monomers (in band 24). We also identified homologues of Arabidopsis PGRL1 protein (82) involved in NDH-independent cyclic electron flow, one accumulating with a very pronounced ratio in the M thylakoid (TC312147; BS/M = 0.26) and the other in the BS membranes (TC305270; BS/M = 6.58).
The Chloroplast Expression Apparatus—
The plastid chromosome binds to the inner envelope and thylakoid membranes (83). We quantified four maize homologues of Arabidopsis proteins found to interact with the chloroplast chromosome. The Arabidopsis matrix attachment filament protein 1 (MFP1; AT3G16000) is a coiled coil protein that is believed to anchor the plastid chromosome to thylakoids (84). We quantified MFP1-like proteins (TC304489 and TC296233) with BS/M ratios of 1.9 and 2.4, respectively. These proteins are present in high molecular mass complexes of
700,
500, and
250 kDa. We also quantified maize homologues of TCP34 (TC312665) and pTAC16 (TC281749) with BS/M ratios of 0.43 and 0.54, respectively. TCP34 is a tetratricopeptide repeat (TPR) protein found in association with a transcriptionally active protein-DNA complex (TAC) from chloroplasts; a regulatory function of TCP34 in plastid gene expression was proposed (85). pTAC16 was identified in TACs from A. thaliana and mustard chloroplasts and was shown to be involved in plastid gene expression (25).
We identified 27 chloroplast ribosomal subunits (mostly 50 S subunits) either only identified in M membranes or identified with a low BS/M accumulation ratio (average 0.5). Also plastid ribosome-associated proteins (PSRPs) showed M preferential accumulation, including PSRP-7 (TC295920 and TC280502), PSRP-2 (TC312649), and elongation factor Tu (TC298927 and TC298925). Two accessions (TC281611 and TC281610) corresponding to the ribosome recycling factor yielded opposite expression ratios (3.0 and M only, respectively), suggesting specialization according to cell type. Chloroplast-encoded membrane proteins are synthesized on ribosomes docked onto the thylakoid membrane (86). It is the 50 S rather than the 30 S ribosomal particle that makes the interaction to the thylakoid membrane (via the Sec translocon); this is consistent with the higher number of identified 50 S subunits.
Chloroplast Protein Import, Membrane Biogenesis, and Protein Targeting—
Nuclear-encoded chloroplast proteins are imported via the Toc-Tic complex, consisting of several outer and inner membrane proteins (33). Components of the Tic complex consistently showed preferential M accumulation (average BS/M of 0.46 or only in M), namely Tic110 (TC300617 and TC300618), Tic21 (TC307215 and TC293198), and Tic40 (TC310831 and TC310825) (Table I). Tic110 was identified in a 170-kDa complex with its Tic40 and Tic21 partners in both 2D BN and label-free native profiles. An additional complex of lower intensity also was observed for these three Tic subunits around 500 kDa. The major outer envelope preprotein receptor Toc159 (TC300201) was also identified, but there was not enough information for quantification (supplemental Table 7). Assuming a similar envelope/thylakoid ratio, the enrichment for Tic components in the M chloroplast suggests more active protein import as compared with BS chloroplasts. It is important to note that other inner envelope proteins, not involved in protein import, are clearly enriched in the BS samples, e.g. chloroplast envelope quinone oxidoreductase (TC311227) (87) with a BS/M ratio of 3.7, HP45 unknown function (TC314108), and the 2-oxoglutarate/malate translocator (DIT2; TC282161), indicating that the low BS/M ratio for Tic components cannot be explained by the amount of envelopes in the purified membranes.
We identified two maize proteins that are homologues to Arabidopsis proteins TF1 (thylakoid formation 1) (88) and VIPP1 (vesicle-inducing protein in plastids) (89), both involved in thylakoid membrane formation. Loss of VIPP1 expression is deleterious to thylakoid membrane formation (89), and VIPP1 appears to be needed for formation of the thylakoid bilayer rather than the functional assembly of thylakoid protein complexes (90). Maize TF1 (TC285636) was only identified with a low BS/M ratio (0.61), whereas maize VIPP1 (TC281010) protein was identified with a high BS/M ratio (2.67) in very high molecular mass complexes (band 1) and as monomers.
Several pathways for insertion and translocation of thylakoid proteins have been discovered in C3 and C4 plants (91, 92). The Sec pathway involving cpSecA,Y,E is used for insertion of both chloroplast- and nucleus-encoded proteins (92, 93). Two accessions for cpSecY (TC288462 and NP004391) were quantified with quite different BS/M ratios of 0.82 and 2.51, respectively (Table I). Other intrachloroplast targeting components such as TatC (TC300579) and two accessions for cpSecA were identified but could not be quantified (supplemental Table 7), whereas cpSRP54 (TC306667) was identified only in the M. The expression of the thylakoid signal peptidase TPP2 (TC279895) involved in cleavage of lumenal transit peptides showed a BS/M of 1.3.
DnaJ Domain Proteins and Immunophilins—
DnaJ proteins are best known as co-chaperones and interacting partners of Hsp70 (DnaK) and GrpE chaperones. The Arabidopsis genome encodes about 90 DnaJ domain proteins (94), and several have been identified in Arabidopsis thylakoid membranes (19). The function of some DnaJ domain proteins in chloroplasts have been determined and include ARC6 involved in plastid division (95), maize BSD2 involved in biogenesis of the Rubisco complex (2), and the Orange 1 protein (Or1) in Brassica oleracea required for carotenoid accumulation in chromoplasts (96).
We identified and quantified the BS/M accumulation patterns for several maize DnaJ domain proteins. Accessions TC292740, TC300370, and TC300372 were enriched in M membranes, whereas TC296329 and TC297356 were strongly enriched in BS membranes. 1D BN profiles indicated that all of these DnaJ proteins are present in high molecular mass complexes at
580 kDa (band 8) and
650 kDa (band 6) situated just above and below the abundant
600-kDa (band 7) band containing PSI-LHCI and PSII dimers. A maize homologue of Or1 (TC288173) was found preferentially expressed in the BS (BS/M = 2.65) and was only found as a monomer (or small complex).
We identified nine thylakoid-localized peptidyl-prolyl isomerases, corresponding to five known Arabidopsis thylakoid lumen proteins. Two of these Arabidopsis proteins, TLP20 and TLP40, are believed to hold most of the isomerase activity (97). In particular the TLP20 homologues appear to be relatively abundant in the maize thylakoid proteome based on their high cumulative Mowse scores (up to 14,000) and spectral counts (up to 600). Two maize homologues (TC295339 and TC295338) of the isomerase protein TLP20 in Arabidopsis had clear preferential BS accumulation (BS/M = 2.3 and 1.9), whereas the maize homologues of TLP40 (TC298326 and TC289507) showed preferential accumulation in the M thylakoid. This M enrichment for TLP40 is in line with its proposed role in light-mediated protein phosphorylation of PSII subunits (98). Also FKBP-type protein (BM499069) showed preferential accumulation in M (BS/M = 0.78), but its role is unknown. Peptidyl-prolyl cis-trans isomerases CO454223, BG321212, and TC296787 homologous to FKBP13 involved in the biogenesis of the Rieske cytb6f protein (99) all showed high BS/M ratios (1.7, 3.5, and 3.1, respectively).
PSII and PSI Assembly Factors—
We identified a maize homologue of the Arabidopsis "low PSII accumulation 1 protein" (LPA1) (TC312066) with good scores in all three replicates of M thylakoids but not in BS thylakoids. LPA1 has two TPR repeats and is implicated in the synthesis of the chloroplast-encoded D1 protein and assembly of PSII in A. thaliana (100). The observation of LPA1 in M thylakoids and not in BS thylakoids is consistent with its proposed role as a specific factor in PSII biogenesis. We also identified the maize homologue of Arabidopsis protein HCF136, shown to be involved in assembly of the PSII reaction center assembly factor (101) and also in maize (4). Maize HCF136 (TC296744) distributed equally between M and BS. The equal distribution of HCF136 suggests that HCF136 is not specific for PSII or that it is a remnant of the PSII synthesis and assembly machinery that was not (yet) lost during BS differentiation and PSII down-regulation. It is interesting to note that HCF136 in M membranes showed a predominant accumulation in a complex of
180 kDa (band 18), whereas HCF136 predominantly accumulated at
100 kDa (band 24) in BS membranes. This difference in oligomeric organization might indicate a partial loss of function in the absence of high levels of active PSII in the BS.
We identified two PSI assembly factors, YCF3 (TC293155) and YCF4 (NP_043035). YCF3 was identified in all three replicates in M (albeit with low scores and thus did not pass the filter for Table I) but not in BS, whereas YCF4 accumulated preferentially in the BS (BS/M = 2.7) and was identified with good scores. Both proteins are essential for the stability and accumulation, but not synthesis, of PSI in tobacco and Chlamydomonas reinhardtii (Refs. 102–104; for a review, see Ref. 30). Sucrose gradient fractionation and immunoblots indicate that YCF4 is part of a high molecular weight complex in C. reinhardtii that does not co-sediment with PSI subunits, whereas YCF3 is monomeric (103). Consistently we found that YCF4 accumulated in a
500 kDa complex (band 10, above PSI monomer) as well as at
220 kDa (band 16) and in low molecular weight fractions (regions below LHCIIs, bands 21–25). In contrast, maize YCF3 was found only in low molecular weight fractions, similar to that observed in Chlamydomonas (103). The contrasting BS/M accumulation patterns and oligomeric state of YCF3 and YCF4 suggest that they have different roles in PSI assembly.
We also quantified the expression of a maize homologue of Arabidopsis pale yellow green7 protein (PYG7; TC307174) with BS/M ratio of 1.53, which is similar to ratios observed for PSI subunits. PYG7 in A. thaliana is required for a specific accumulation of PSI subunits without affecting the accumulation of other photosynthetic complexes and was qualified as a PSI assembly factor (105). The protein contains a TPR motif, and it was only identified in low molecular weight complexes on 1D BN. Thylakoid-bound Rubredoxin A (RubA) was proposed to be specifically required for the assembly of the F(X) iron-sulfur cluster in PSI but was not needed for iron-sulfur clusters in other thylakoid proteins (e.g. Rieske protein). Homologues of RubA were detected in thylakoids of spinach and C. reinhardtii (106, 107). We identified a maize homologue of RubA (TC296284) expressed in both compartments (BS/M = 0.8).
A maize homologue of the Arabidopsis protein APE1 (acclimation of photosynthesis to environment 1) (TC299113) showed an equal distribution between BS and M membranes. A mutant in APE1 was isolated in a chlorophyll fluorescence screen for light acclimation mutants (108), and the protein was identified by proteomics in the thylakoid membrane of Arabidopsis (19).
Protein Degradation—
Protein degradation has a fundamental role in chloroplast development and maintenance as evidenced by strong phenotypes of chloroplast protease mutants in Arabidopsis (109, 110), and it is thus likely that proteases also play a role in BS/M differentiation. Chloroplast proteases in Arabidopsis include soluble ATP-dependent ClpP protease, thylakoid ATP-dependent FtsH metalloproteases, thylakoid SppA and EGY1, and members of the DegP family (109, 110).
Three maize accessions (TC287566, BE453757, and TC303156) with homology to thylakoid-localized SppA protease (111) accumulated preferentially in the M membranes (average BS/M ratio = 0.33). It has been suggested that this high light-inducible protease might target PSII and LHCII proteins in A. thaliana (111). We identified maize SppA mainly in a
500-kDa high molecular mass complex that could correspond to a dimeric form of a previously observed Arabidopsis SppA complex of 270 kDa (111). EGY1 was identified as a metalloprotease required in early chloroplast development (112), and we observed two EGY accessions (matching to EGY1 and -2) with high BS/M ratios (TC314564 and TC287793) of 3.23 and 4.9, respectively, suggesting substrates specifically enriched in BS thylakoids.
Various thylakoid-associated DegP protease are likely involved in degradation of the D1 protein (113–115). We identified lumenal DegP1 (TC294056) in low molecular weight fractions with a BS/M ratio of 0.56. Six accessions for FtsH proteases quantified by label-free analysis, iTRAQ, and 2D BN showed that FtsH protease is present in both BS and M chloroplasts with equal distribution (TC292185, TC292243, TC234343, TC286826, TC286827, and TC219258). FtsH is mainly present in a high molecular mass complex of
750 kDa and in
650- and
500-kDa complexes.
We did not identify any components of the ClpP protease complex, such as ClpP,R or ClpS (110). However, we did observe ClpC (TC281093), one of three closely related Clp chaperones, in a distinct major peak at
650 kDa (band 5) that could correspond to the predicted hexameric assembly. A putative oligopeptidase (TC294060) was identified with a high Mowse score (2200 in M and 103 in BS) and a BS/M ratio of 0.05. This is one of the most extreme ratios observed in our data set, suggesting a very specific function in M chloroplasts.
Redox and ROS Response Proteins—
Chloroplasts have an elaborate enzymatic system to detoxify ROS (116–119). M chloroplasts have high rates of linear photosynthetic electron transport, whereas BS chloroplasts mainly carry out cyclic electron flow. It is thus quite likely that there are substantial differences in the quantity and "quality" of ROS produced in these two chloroplast types. Indeed the majority of ROS defense proteins quantified in our study showed higher accumulation in M membranes than in BS membranes, similar to those observed for stromal M and BS maize proteins (8). For example, the thylakoid-bound APX (TC306001), the lumenal Peroxiredoxin Q (Prx-Q; TC306534), and two accessions for peripheral thylakoid protein 2-Cys Prx-B (TC304817 and TC304818) have BS/M ratios of 0.67, 0.14, 0.8, and 0.45, respectively. The stromal glutathione peroxidase 2 (GPX2; TC282772) had a similar BS/M ratio of 0.53. Interestingly PrxII-E (TC306147) seems to have a specialized BS function as indicated by its very high BS/M of 6.7. Neither the precise function nor the electron donor of PrxII-E is known in Arabidopsis (119).
The expression of ferritins in maize is regulated by abscisic acid in the case of Ferritin 2 (ZmFer2) and by oxidative stress in the case of ZmFer1 (for review, see Ref. 120). The preferential accumulation of Ferritin 1 (TC293195) in the M membrane is consistent with these observations and supports a role in iron chelation to reduce formation of reactive hydroxyl radicals by interaction of O2 with ferrous ions as suggested previously (120).
Calvin Cycle and C4 Markers—
Because the BS and M membrane fractions analyzed here were deliberately not washed with salts, abundant stromal enzymes involved in carbon fixation and C4-specific pathways could be identified and quantified in this study. As we mentioned earlier (Fig. 3B), these proteins provided an excellent test set to determine the consistency between the current study and our previous study regarding the BS and M maize stromal proteomes (8).
The specific Calvin cycle enzymes RBCL (NP_043033), RBCS (TC286735, TC286728, and TC286733), RCA (TC300568), and fructose-bisphosphate aldolase-2 (SFBA; TC285785) all showed strong BS accumulation. In agreement with our previous BS/M stromal analysis, we observed that the Calvin cycle enzymes involved in triose phosphate reduction, (phosphoglycerate kinase, NP539335, TC286022, and TC219625; glyceraldehyde-3-phosphate dehydrogenase A,B, TC292245 and TC292248; and triose-phosphate isomerase, TC279906) preferentially localized to M or were equally distributed over BS and M chloroplasts (Table I). Within the malate-pyruvate C4 shuttle, we observed strong M expression for PPDK (TC233444 and TC286559) and MDH (TC300313). Carbonic anhydrase (TC280282) was only identified in M membranes; this is consistent with activity assays (121). As expected the NADP-malic enzyme (TC280868) was only identified in the BS.
The activity of PPDK is regulated by a bifunctional kinase/phosphatase PPDK-regulatory protein (PPDK-RP) (122). This unusual enzyme uses ADP as source of Pi (123). Consistently we observed that adenylate kinase (ADK; DR786994, TC301869 and TC293517), which phosphorylates the AMP produced in this PPDK phosphorylation reaction, preferentially accumulated in the M. In agreement with our stromal analysis, we identified PPDK-RP (TC287496 in ZmGI v16.0) in M chloroplasts and not in BS chloroplasts consistent with the localization of its substrate PPDK. 1D BN profiles showed a co-migration of PPDK and PPDK-RP in high molecular mass complexes at >1000 kDa (band 1),
650 kDa (band 5), and
500 kDa (band 10) (not shown). The relative stoichiometry of PPDK and PPDK-RP appears to be conserved, strongly suggesting stable interaction. An additional
180-kDa (band 18) complex was observed for PPDK only, possibly representing homodimeric PPDK.
Translocation of Triose Phosphates and Starch—
M chloroplasts have to maintain high rates of export of phosphoenolpyruvate (PEP) (generated by PPDK) into the cytosol for oxaloacetate production by PEP carboxylase. Consistently we observed that the PEP/Pi translocator (PPT), located in the inner chloroplast envelope membrane, accumulated with a low BS/M ratio of 0.21 (TC286421). PPT homologues in Arabidopsis were shown to be expressed specifically in the mesophyll cells (PPT1) or veins (PPT2) and are required for formation of mesophyll tissue (124).
Interestingly we also identified two close maize homologues of inner envelope 2-oxoglutarate/malate translocator (Dit). TC282039 was closer to A. thaliana Dit1 (AT5G64290), TC282161 was closer to Dit2 (AT5G64280), and they preferentially accumulated in the M envelopes (BS/M = 0.2) and BS envelopes (BS/M = 1.7), respectively. This observation is consistent with previous observations of cell-specific transcript accumulation of Dit1,2 in sorghum (125) and maize (41). Consequently a two-translocator model for Dit1,2 was proposed (126) in which M-enriched Dit1 imports 2-oxoglutarate into chloroplast in exchange for malate, whereas BS enriched Dit2 imports malate and exports glutamate (for a review, see Ref. 127). In C3 plants, Dit1,2 serve to facilitate the photorespiratory pathway, and indeed the Dit2-null mutant is a classic photorespiratory mutant in Arabidopsis (125). Because photorespiration is strongly suppressed in maize (and other C4 species), the Dit translocators in maize essentially serve to drive malate from the M cell to the BS chloroplast.
The maize inner envelope maltose transporter Mex1 (TC302035) accumulated strongly in the BS membranes (BS/M ratio of 2.7). Arabidopsis Mex1 was shown to function as the main exporter of maltose to the cytosol (128) and is essential for carbohydrate export in Arabidopsis leaves at night. Its BS localization is consistent with previous localization of starch synthesis and degradation mechanisms in maize BS chloroplasts (see Ref. 8 for a discussion and references). Clearly, differentiation of envelope proteomes of BS and M chloroplasts participates in the adaptation to C4 photosynthesis.
Isoprenoid and Tetrapyrrole Pathways (Carotenoid, Plastoquinone, and Chlorophyll Biosynthesis)—
Within the initial steps of the carotenoid biosynthetic pathway, we found two isoforms of phytoene desaturase (TC306129 (BS/M = 1.41) and BG320144 (BS/M = 0.64)), a
-carotene desaturase (TC289473; BS/M = 0.42), and, further downstream, zeaxanthin epoxidase (TC292044; BS/M = 0.23). A β-hydroxylase (DT944642), a homologue of A. thaliana LUT5 (129), was identified in the M thylakoid. A maize homologue of Arabidopsis inner envelope membrane protein APG1 (TC304918) accumulated at higher levels in the M membranes. APG1 likely functions as an S-methionine methyltransferase in the methylation step of plastoquinone (PQ) biosynthesis (130).
We identified seven enzymes participating in the chlorophyl biosynthetic pathway. The initial steps in the pathway are carried out by soluble stromal enzymes and were not identified in the membrane samples in our current study. The most upstream enzyme in our study was envelope-associated Mg2+-protoporphyrin IX chelatase subunit H (CF244872 and TC287014) (CHLH, also named Gun5 in Arabidopsis), only identified in the M fraction with low scores. Mg-protoporphyrin IX monomethyl ester cyclase (TC299016), two steps further downstream, showed preferential BS accumulation (BS/M = 3.0). Homologues of this enzyme in Arabidopsis and C. reinhardtii are thylakoid- and envelope-associated CHL27 (131) and CRD1 (132), respectively. Two more steps further down in the pathway, we observed NADPH-protochlorophyllide oxidoreductase A (PORA; TC292871) with a BS/M ratio of 0.55. Geranylgeranyl reductase (TC311340) and chlorophyll synthase (TC300395) showed a BS/M ratio of 1.2 and 0.8, respectively. TC305891, a close homologue of the Arabidopsis (AT3G14110) FLU protein and responsible for negative feedback regulation within the tetrapyrrole pathway (133), accumulated with a BS/M ratio of 1. Finally maize lethal leaf spot protein (134), a pheophorbide a oxygenase (PaO; TC282596) involved in chlorophyll degradation, had a BS/M ratio of 1.26. Thus with exception of the CRD1/CHL27 homologue (see "Discussion"), enzymes in the chlorophyll metabolism were quite evenly distributed.
Jasmonic Acid Synthesis—
Jasmonic acid (JA) is an important hormone in plant development and defense, and we quantified the first two committed steps of JA synthesis. Phospholipase (Defective in Anther Dehiscence 1 protein or DAD1; TC293375), which performs the initial step in JA synthesis, showed a BS/M ratio of 0.38, whereas Lipoxygenase 2 (LOX2; TC298873), which carries out the second step in the JA biosynthetic pathway, also preferentially accumulated in M thylakoids (BS/M = 0.58).
The Plastoglobular Proteome—
Chloroplast plastoglobules (PGs) are thylakoid-associated lipoprotein particles; their size and number vary greatly during chloroplast differentiation, senescence, and abiotic stress (for a review, see Ref. 135). Analyses of PGs isolated from Arabidopsis thylakoids show that they have a unique proteome and are a transient storage site for e.g. prenyl quinones and tocopherol (vitamin E) (135). We quantified 11 maize accessions with homology to these Arabidopsis PG proteins. Two structural lipid coat proteins of the plastoglobulin family were expressed in the M (TC299786 and TC284028), whereas two others were preferentially expressed in the BS (TC281878 and TC307124). This family of proteins is believed to control PG size and stability, and in Arabidopsis the family members respond differentially to various environmental conditions and developmental stage (136). Additionally we identified four ABC1 kinases with ABC1K5 (TC281706) only found in M samples and ABC1K2 (TC296339) strongly expressed in the BS (BS/M = 4.8). The function of the kinases is unclear, but homologues in yeast mitochondria and in E. coli are believed to regulate quinone synthesis (for a brief discussion, see Ref. 23). The BS/M ratio of some other plastoglobuli components could be quantified and included a flavin reductase-like protein (TC305628) and an aldo/keto reductase-like protein (TC282614), both with preferential expression in M membranes (Table I).
Additional Quantified Proteins with Unknown Roles—
Numerous additional proteins without a known function or predicted functional domain were quantified, showing a wide expression between BS and M (Table I). Arabidopsis homologues for many of these proteins have been localized previously in the thylakoid or chloroplast envelope by proteomics (see the PPDB). The proteins with strong BS/M cell-specific accumulation patterns likely play a role in functional BS/M differentiation and represent new entry points to study C4 differentiation.
| DISCUSSION |
|---|
|
|
|---|
To integrate the quantitative and qualitative data obtained in this study, we summarized several of our findings in Fig. 6. This provides a schematic representation of the differential expression of photosynthetic and auxiliary thylakoid functions between the M and the BS as well as an integrated organization of the C4 carbon fixation and main chloroplast metabolic pathways. In the following section we discuss and conceptualize our findings particularly in the context of follow-up experiments to address mechanisms of BS/M differentiation and with relevance to molecular engineering of C4 photosynthesis (138, 139).
|
Many proteins are members of (small) gene families, and similar to our previous analysis of the BS- and M-specific stromal proteomes (8), we observed that different paralogues (e.g. for phosphoglycerate kinase) often show cell-specific accumulation. For instance, two LHCII members accumulated with higher BS/M ratios than the other 16 LHCII accessions. These BS enriched LHCII proteins could represent LHCII proteins preferentially interacting with PSI. In addition, several subunits of PSI (e.g. PsaE and PsaH) were represented by more than one accession, some of which showed differences in BS/M ratios. Based on observations of neofunctionalization between paralogues of PsaE and PsaH in Arabidopsis (140), this may suggest subtle differences in PSI composition between BS and M thylakoids, for instance to adapt to higher rates of cyclic electron flow in BS membranes. A high quality annotated maize genome sequence and additional experiments will be required to sort out these subtle but interesting differences.
Characterization of NDH Complexes and Associated Proteins: a Role for NDH in CO2 Concentration in Addition to Cyclic Electron Flow?—
Thylakoid-localized Ndh complex is proposed to play a role in chlororespiration and in cyclic electron flow around PSI in both C3 and C4 plants (Refs. 72, 141, and 142; for discussions, see Refs. 16 and 143). During chlororespiration in the dark (or light), the NDH complex oxidizes stromal metabolites and transfers these electrons to PSI (or IMMUTANS) in a process coupled to proton translocation to help generate ATP. The role of NDH in the light-driven cyclic electron flow is to accept electrons (likely via Fd rather than FNR (see Ref. 69)) from PSI and then cycle the electrons back to PSI via intersystem components in a process coupled to proton translocation to generate ATP (144). An alternative cyclic electron flow route does not involve NDH but requires PGR5, PGRL1, Fd, and PQ as well as the cytb6f complex (16, 82, 145, 146); all components were detected in maize thylakoids.
Our analysis identified eight of 11 chloroplast-encoded NDH subunits as well as homologues of four (putative) nucleus-encoded NDH subunits (NdhM, -N, and -O and PPL2) recently discovered in Arabidopsis (62, 69) and two putative maize homologues of Thermosynechococcus elongatus NdhL and NdhU. These 14 subunits accumulated with an average BS/M ratio of 3.0, which is in good agreement with Western blot data for NDH-H (16) and BS/M ratios observed for several other subunits (17, 18). These subunits likely form the L-shaped complex observed in cyanobacteria and complex I in mitochondria (Ref. 147 and references therein). In addition, we consistently found associated lumen-localized isomerases TLP20 and FKPB; together they likely function in protein folding, and the linkage to NDH possibly facilitates oxidation or reduction of lumenal proteins (148). The observed native masses of the NDH complex (550 kDa for the monomer and 1000 kDa for the dimer) and the various subcomplexes are consistent with published analyses (15, 17, 18). Three additional Arabidopsis proteins, CRR1, -4, and -6 (67, 149, 150), were shown to be involved in NDH biogenesis but were not considered bona fide NDH subunits, and we did not identity their homologues in the maize samples.
In cyanobacteria, specialized low CO2-induced thylakoid-localized NDH complexes (D3/F3 type) participate in carbon-concentrating mechanisms (CCMs) (151, 152). Recent analysis of their structural organization by single particle electron microscopy in T. elongatus (147) showed (i) an "L"-shaped complex composed of the connecting hydrophilic arm and the membrane-embedded subunits, ii) a "U"-shaped configuration with an additional hydrophilic arm, and (iii) a smaller complex on the lumenal side. The presence of the U-shaped NDH form is unlikely to be an artifact because milder purification conditions dramatically increased the population of U-shaped particles compared with L-shaped particles (147). In fact, the authors suggested that the U-shaped particle is the functional form of the NDH-1 in T. elongatus. The composition of the extra arm in the U-shaped NDH-1 complex is unknown, but genetics showed that ChpY (also named CupA) is an essential component in thylakoid-associated CCM (152). Surprisingly we identified an additional subcomplex in BS thylakoids with four identified proteins that was consistently found in higher molecular mass assemblies containing NDH subunits. Based on the presence of the U-shaped complex discussed above, we speculate that the 1000-kDa thylakoid NDH complex represents a functional homologue of the U-shaped NDH-1 complex in T. elongatus. Thus the additional abundant (peripheral) complex in the BS that we clearly observed on the 2D BN-PAGE gels as well as in the 1D BN analysis may be the functional equivalent of the extra "arm" involved in CO2 hydration in cyanobacteria. We speculate that this novel complex may be involved in hydration of CO2 into HCO3–, thus strongly reducing diffusion rates of inorganic carbon through chloroplast and cellular membranes and providing a reserve/buffer of inorganic carbon within the BS chloroplast (Fig. 6). This would prevent wasteful leakage of CO2 from the BS chloroplast, in particular when ME decarboxylation rates are higher than Rubisco carboxylation rates. A likely source of NADPH to perform this function would be the decarboxylation of malate by ME furnishing CO2 and NADPH in a 1:1 ratio. The benefits for C4 photosynthesis are clear, and this hypothesis is consistent with the low BS/M ratio for carbonic anhydrase (TC280282; only identified in the M with good scores). The preferential location of carbonic anhydrase in M chloroplasts may further be enhanced by the preferential location of lipid biosynthesis in M chloroplasts with a much lower internal CO2 concentration because HCO3– is needed as a co-substrate for acetyl-CoA carboxylase (153). Thus we suggest that maize may have retained a thylakoid-associated CMM that is best understood in cyanobacteria (152, 154), the progenitor of plant plastids.
Control and Regulation of Excitation Energy and ROS Defense—
The redox state, partial pressures of molecular oxygen, and the quality and quantity of the light absorption cross-sections are different between M and BS chloroplasts. Hence there must be BS- and M-specific differentiation to avoid light stress and optimize light harvesting capacity. This includes balancing excitation of PSI and PSII reaction centers through (de)phosphorylation-driven state transitions, non-photochemical quenching of excess excitation energy in the antennae, detoxification of triplet chlorophyll, and detoxification of the various ROS using enzymatic and non-enzymatic systems (116–119). We were able to quantify several membrane (-associated) proteins involved in these processes, and these data were complementary to findings of the BS and M stromal proteomes (8). In particular, the very low BS/M ratio for lumenal Prx-Q suggests a specific role in ROS defense in the M thylakoid lumen, whereas the very high BS/M ratio for PrxII-E on the stromal side of the thylakoid membrane suggests involvement with metabolic redox regulation in BS chloroplasts. The strong M enrichment of maize homologues of light stress proteins Ohp1,2, Lil4/Sep1, and PsbS suggests that their functions are related to PSII and/or that they are generally up-regulated under conditions of high linear electron transport. The preferential M accumulation of thylakoid kinases STN7 and STN8 and phosphoprotein TSP9 as well as lumenal isomerase TLP40 strongly suggests that M thylakoids retained phosphorylation-driven regulation of LHCII state transitions and PSII repair, whereas BS thylakoids have down-regulated these functions in the absence of significant levels of functional PSII.
Control of Protein Homeostasis in BS and M Chloroplast Differentiation—
To accommodate protein homeostasis in the two cell types, it is expected that the BS and M chloroplasts also display differences in the protein import and sorting, folding, assembly, and degradation machineries. Inspection of BS/M accumulation ratios for proteins involved in these different processes suggested indeed a number of divergences. High levels of import machinery in the inner envelope (Tic110,21,40) of M chloroplasts suggested an increased protein flux. Consistently, cpSRP45 involved in targeting of LCHC and chloroplast-encoded proteins (92) was enriched in M thylakoids. The low BS/M ratios for ribosomal and ribosome-associated proteins strongly suggest higher rates of translation (of thylakoid proteins) in M chloroplasts; this would be consistent with higher accumulation levels of PSII and high rates of damage and resynthesis of the D1 protein. This hypothesis is supported by protein 35S pulse labeling in the C4 plant Digitaria sanguinalis showing high rates of synthesis of D1 and D2 PSII subunits in the M, whereas only labeled RBCL is visible in BS chloroplasts (155).
Finally several thylakoid proteins were quantified that are believed to be involved in maturation or assembly of thylakoid complexes, such as lumenal isomerases. Of particular interest were maize homologues for Arabidopsis LPA1 and HCF136. Although we observed enrichment in the M membrane of LPA1 consistent with its role in synthesis of the D1 protein (100), HCF136 clearly involved in assembly of PSII reaction centers in Arabidopsis (101) and maize (4) showed a BS/M ratio close to 1, suggesting additional roles. Also interesting to note is that although PSII assembly factors Pyg7 and RubA have BS/M ratios similar to that of PSI itself, YCF3 and YCF4, have very different ratios. Thus the comparative analysis of BS and M thylakoids provides a good system to study the functions and specificities of various assembly factors.
Numerous studies in Arabidopsis indicated that thylakoid FtsH and Deg proteases are involved in the degradation of the D1 protein (113, 156, 157). However, we did not observe that the FtsH family was particularly enriched in either cell type, whereas lumenal DegP1 was nearly 2-fold higher in M thylakoids. Two maize homologues of Arabidopsis thylakoid protease EGY1 (112) and uncharacterized EGY2 showed more than 3–5-fold higher accumulation in BS membranes, whereas a putative oligopeptidase (TC294060) showed one of the lowest BS/M ratios observed in our data set, suggesting very specific functions in BS and M chloroplasts, respectively. Studies that follow cell type-specific accumulation of proteases along the developmental gradient could further explore the role of proteolysis in BS/M differentiation.
Cell-specific Differentiation of Chlorophyll, Carotenoid, and Quinone Biosynthetic Pathways—
The thylakoid membrane and its proteome contain a high number of various carotenoids and quinones (plastoquinone and phylloquinone) as well as tocopherol. The final steps in their biosynthesis occur at the thylakoids and/or PGs (135), and therefore one would expect that the abundance and/or activity of the enzymes in the BS or M membrane is proportional to demand. Demand is determined by the composition of the electron transport chain which suggests higher amounts of plastoquinone, carotenoids, and tocopherol in the M thylakoid. Indeed relevant biosynthetic enzymes were either only found in M membranes or found with low BS/M ratios. Although we quantified a number of PG proteins, it was particularly interesting to observe very divergent BS/M ratios in the family of ABC1 kinases. In particular ABC1K2 showed a very high BS/M ratio, whereas ABC1K5 was only found in M membranes with good scores. The function of the ABC1 kinase family is unknown, but homologues in E. coli and yeast are involved in regulation of (ubi)quinone synthesis (see Ref. 23). The observed BS/M differentiation provides an excellent opportunity to further analyze ABC1K and PG function.
In contrast, enzymes in chlorophyll synthesis and degradation were equally distributed over the two cell type with the exception of BS enriched CHL27, which requires molecular oxygen and NADPH in the cyclase reaction. Because the O2 concentration is lower in BS than in M, this higher BS accumulation may serve to compensate for slower reaction kinetics and is consistent with enhanced expression of homologue CRD1 in C. reinhardtii when grown under low O2 (158).
Regulation of BS/M Differentiation—
The molecular basis for C4 differentiation is still poorly understood, but it has been established that it includes nuclear transcriptional regulation through DNA regulatory elements, transcription factors, histone modifications, and likely also metabolic signals (5, 159). Nuclear gene expression is also influenced by signals from the plastid, but the actual signals have not been identified conclusively (160, 161). However, it is clear that multiple components within the plastid contribute to plastid-nucleus signaling (in both C3 and C4 plants), possibly integrated in a single pathway (via pentatricopeptide repeat protein GUN1) when the signal is "leaving" the plastid (162). These plastid components/contributions have been suggested to include (i) the redox state of the PQ pool and possible other redox-active components in the thylakoid and stroma, (ii) translational activity within the plastid, and (iii) intermediates of the tetrapyrrole pathway. It is conceivable that intraplastid signaling also involves phosphorylation cascades, various ROS species, and other small molecules.
We explored our data for possible evidence of differences in BS and M plastid-nucleus signaling, and indeed we did observe several components and processes that have been implicated in signaling. Within the tetrapyrrole/chlorophyll synthesis in particular CHL27 stands out for its high BS/M ratio. Clearly the BS and M differences in linear electron flow rates, tightly intertwined with differences in the ROS defense system, possibly provide a good starting point for differential plastid-nucleus signaling. Also redox equivalents imported via the malate shuttle and possibly redistributed via chlororespiration (involving the NDH complex) are likely to contribute to plastid-nucleus signaling. Although these components may help to complete the BS/M differentiation process and to enforce the differentiated state (the status quo), they are not the primary event in the differentiation process. The search for components involved in the primary steps in BS/M differentiation should involve leaf zones closer to the base of the leaves. It is interesting to note that when the plastid PsbA gene encoding for the D1 protein was deleted in tobacco the NDH complex was up-regulated 4-fold (163). This suggests that parts of regulatory networks or signal transduction chains involved in C4 chloroplast differentiation maybe be similar to those in C3 plants, also underlining the notion that genes necessary for C4 photosynthesis are present in C3 species (137). Finally two isoforms of haloacid dehalogenase were identified, one with no differential expression and the other with strong expression in the BS. In recent years, metabolism of the disaccharide trehalose in plants has received a lot of attention in particular because trehalose 6-phosphate is involved in the regulation of developmental and metabolic processes (164–167). A possible function of the trehalose pathway in BS and M differentiation has yet to be explored.
Conclusions and Challenges—
The current study complements our previous stromal analysis, completes a first quantitative overview of the differentiated state of M and BS chloroplasts, and shows a fascinating number of differences to meet the demand for carbon and energy as well as (re)distribution of metabolites between C4 photosynthetic compartments. The identification and quantification of hundreds of chloroplast proteins integrates observations from many studies over the last few decades. Moreover our study provides many new entry points to study the mechanisms and consequences of C4 differentiation in maize and proposes a new function for the NDH complex in C4 plants. A logical next step will be to study the kinetics of BS/M differentiation along the leaf gradient from base to tip; this is needed to establish the sequence of events. Our study also demonstrated the strength of label-free quantification using a highly sensitive, fast MS instrument with high mass accuracy, the LTQ-Orbitrap. Finally a completely sequenced maize genome, high quality assembly, and gene discovery will be needed to refine this analysis. Our MS/MS data set can help to annotate gene models as we demonstrated earlier on a small scale for Arabidopsis (20) and as demonstrated in large scale studies in other organisms such as Drosophila melanogaster (168).
| FOOTNOTES |
|---|
Published, MCP Papers in Press, May 2, 2008, DOI 10.1074/mcp.M800016-MCP200
1 The abbreviations used are: BS, bundle sheath; M, mesophyll; ME, NADP-malic enzyme; PS, Photosystem; cyt, cytochrome; LHC, light-harvesting complex; NDH, NAD(P)H dehydrogenase; FNR, ferredoxin-NAD(P)+ reductase; Fd, ferredoxin; ROS, reactive oxygen species; LTQ, linear ion trap triple quadrupole; PPDB, Plant Proteomics Database; BN, blue native; β-DM, β-n-dodecyl maltoside; FA, formic acid; cv, coefficient of variation; PPDK, pyruvate, phosphate dikinase; MDH, malate dehydrogenase; GAPB, glyceraldehyde dehydrogenase subunit B; RBCL, Rubisco large subunit; RBCS, Rubisco small subunit; RCA, Rubisco activase; PG, plastoglobule; PEP, phosphoenolpyruvate; PQ, plastoquinone; TPR, tetratricopeptide repeat; TAC, transcriptionally active chromosome; PSRP, plastid ribosome-associated protein; Rubisco, ribulose-bisphosphate carboxylase/oxygenase; iTRAQ, isobaric tags for relative and absolute quantification; 2D, two-dimensional; 1D, one-dimensional; EST, expressed sequence tag; TIGR, The Institute for Genomic Research; Mowse, molecular weight search; ZmGI, Z. mays Gene Index; FPR, false positive rate; BLAST, Basic Local Alignment Search Tool; OEC, oxygen-evolving complex; FKBP, FK506-binding protein; CRR, chlororespiratory reduction; PIFI, postillumination fluorescence increase; IM, IMMUTANS; MFP1, matrix attachment filament protein 1; LPA1, low PSII accumulation 1 protein; PYG7, pale yellow green7 protein; RubA, Rubredoxin A; APX, ascorbate peroxidase; Prx, Peroxiredoxin; PPDK-RP, PPDK-regulatory protein; ADK, adenylate kinase; PPT, PEP/Pi translocator; JA, jasmonic acid; CCM, carbon-concentrating mechanism. ![]()
2 B. Zybailov, unpublished data. ![]()
3 Q. Sun, B. Zybailov, and K. J. van Wijk, unpublished data. ![]()
* The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. ![]()
S The on-line version of this article (available at http://www.mcponline.org) contains supplemental material. ![]()
Present address: Dept. of Chemistry and Biochemistry, University of California, Los Angeles, CA 90095-1569. ![]()
** To whom correspondence should be addressed: Dept. of Plant Biology, Emerson Hall 332, Cornell University, Ithaca, NY 14853. Tel.: 607-255-3664; Fax: 607-255-3664; E-mail: kv35{at}cornell.edu
| REFERENCES |
|---|
|
|
|---|
pH-dependent pathways.
Biochim. Biophys. Acta 1541, 80
–90[Medline]This article has been cited by other articles:
![]() |
J. Joyard, M. Ferro, C. Masselon, D. Seigneurin-Berny, D. Salvi, J. Garin, and N. Rolland Chloroplast Proteomics and the Compartmentation of Plastidial Isoprenoid Biosynthetic Pathways Mol Plant, November 12, 2009; (2009) ssp088v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Suorsa, S. Sirpio, and E.-M. Aro Towards Characterization of the Chloroplast NAD(P)H Dehydrogenase Complex Mol Plant, November 1, 2009; 2(6): 1127 - 1140. [Abstract] [Full Text] [PDF] |
||||
![]() |
Z. Shen, P. Li, R.-J. Ni, M. Ritchie, C.-P. Yang, G.-F. Liu, W. Ma, G.-J. Liu, L. Ma, S.-J. Li, et al. Label-free Quantitative Proteomics Analysis of Etiolated Maize Seedling Leaves during Greening Mol. Cell. Proteomics, November 1, 2009; 8(11): 2443 - 2460. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Peng, Y. Fukao, M. Fujiwara, T. Takami, and T. Shikanai Efficient Operation of NAD(P)H Dehydrogenase Requires Supercomplex Formation with Photosystem I via Minor LHCI in Arabidopsis PLANT CELL, November 1, 2009; 21(11): 3623 - 3640. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Ishida, A. Takabayashi, N. Ishikawa, Y. Hano, T. Endo, and F. Sato A Novel Nuclear-Encoded Protein, NDH-Dependent Cyclic Electron Flow 5, is Essential for the Accumulation of Chloroplast NAD(P)H Dehydrogenase Complexes Plant Cell Physiol., February 1, 2009; 50(2): 383 - 393. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Sirpio, Y. Allahverdiyeva, M. Holmstrom, A. Khrouchtchova, A. Haldrup, N. Battchikova, and E.-M. Aro Novel Nuclear-encoded Subunits of the Chloroplast NAD(P)H Dehydrogenase Complex J. Biol. Chem., January 9, 2009; 284(2): 905 - 912. [Abstract] [Full Text] [PDF] |
||||
![]() |
P.-A. Christin, N. Salamin, E. A. Kellogg, A. Vicentini, and G. Besnard Integrating Phylogeny into Studies of C4 Variation in the Grasses Plant Physiology, January 1, 2009; 149(1): 82 - 87. [Full Text] [PDF] |
||||
![]() |
Q. Sun, B. Zybailov, W. Majeran, G. Friso, P. D. B. Olinares, and K. J. van Wijk PPDB, the Plant Proteomics Database at Cornell Nucleic Acids Res., January 1, 2009; 37(suppl_1): D969 - D974. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| All ASBMB Journals | Journal of Biological Chemistry |
| Journal of Lipid Research | ASBMB Today |