The Metabolic Status Drives Acclimation of Iron Deficiency Responses in Chlamydomonas reinhardtii as Revealed by Proteomics Based Hierarchical Clustering and Reverse Genetics*

Iron is a crucial cofactor in numerous redox-active proteins operating in bioenergetic pathways including respiration and photosynthesis. Cellular iron management is essential to sustain sufficient energy production and minimize oxidative stress. To produce energy for cell growth, the green alga Chlamydomonas reinhardtii possesses the metabolic flexibility to use light and/or carbon sources such as acetate. To investigate the interplay between the iron-deficiency response and growth requirements under distinct trophic conditions, we took a quantitative proteomics approach coupled to innovative hierarchical clustering using different “distance-linkage combinations” and random noise injection. Protein co-expression analyses of the combined data sets revealed insights into cellular responses governing acclimation to iron deprivation and regulation associated with photosynthesis dependent growth. Photoautotrophic growth requirements as well as the iron deficiency induced specific metabolic enzymes and stress related proteins, and yet differences in the set of induced enzymes, proteases, and redox-related polypeptides were evident, implying the establishment of distinct response networks under the different conditions. Moreover, our data clearly support the notion that the iron deficiency response includes a hierarchy for iron allocation within organelles in C. reinhardtii. Importantly, deletion of a bifunctional alcohol and acetaldehyde dehydrogenase (ADH1), which is induced under low iron based on the proteomic data, attenuates the remodeling of the photosynthetic machinery in response to iron deficiency, and at the same time stimulates expression of stress-related proteins such as NDA2, LHCSR3, and PGRL1. This finding provides evidence that the coordinated regulation of bioenergetics pathways and iron deficiency response is sensitive to the cellular and chloroplast metabolic and/or redox status, consistent with systems approach data.

Iron is a crucial cofactor in numerous redox-active proteins operating in bioenergetic pathways including respiration and photosynthesis. Cellular iron management is essential to sustain sufficient energy production and minimize oxidative stress. To produce energy for cell growth, the green alga Chlamydomonas reinhardtii possesses the metabolic flexibility to use light and/or carbon sources such as acetate. To investigate the interplay between the iron-deficiency response and growth requirements under distinct trophic conditions, we took a quantitative proteomics approach coupled to innovative hierarchical clustering using different "distance-linkage combinations" and random noise injection. Protein co-expression analyses of the combined data sets revealed insights into cellular responses governing acclimation to iron deprivation and regulation associated with photosynthesis dependent growth. Photoautotrophic growth requirements as well as the iron deficiency induced specific metabolic enzymes and stress related proteins, and yet differences in the set of induced enzymes, proteases, and redoxrelated polypeptides were evident, implying the establishment of distinct response networks under the different conditions. Moreover, our data clearly support the notion that the iron deficiency response includes a hierarchy for iron allocation within organelles in C. reinhardtii. Importantly, deletion of a bifunctional alcohol and acetaldehyde dehydrogenase (ADH1), which is induced under low iron based on the proteomic data, attenuates the remodeling of the photosynthetic machinery in response to iron deficiency, and at the same time stimulates ex-pression of stress-related proteins such as NDA2, LHCSR3, and PGRL1. This finding provides evidence that the coordinated regulation of bioenergetics pathways and iron deficiency response is sensitive to the cellular and chloroplast metabolic and/or redox status, consistent with systems approach data. The green alga Chlamydomonas reinhardtii has an enormous metabolic versatility (1) and possesses the flexibility to grow in the presence of different carbon sources. It may use carbon dioxide (CO 2 ) for photoautotrophic, acetate for heterotrophic, and both carbon sources for mixotrophic growth. In this alga CO 2 is fixed via the Calvin Benson Bassham cycle (2), while acetate can be taken up, converted to acetyl-CoA, and enter the glyoxylate cycle where it may be incorporated into C4 acids (3). In addition to the use of acetate as a source of energy and carbon backbone for biosynthetic processes, acetate can control respiration and photosynthesis in conjunction with the light intensity and CO 2 availability (4 -6). Moreover, acclimation responses to iron-and copper-deficiencies significantly vary in photoautotrophic versus heterotrophic conditions (7)(8)(9)(10), indicating that the metabolic status of the cells influence overall cellular acclimation responses.
Transition metals like copper, manganese, and iron possess the ability to donate and accept electrons, making these metals suitable cofactors in enzymes that catalyze redox reactions. In particular, iron is used as a cofactor in numerous biochemical pathways and is therefore an essential nutrient. Cells require relatively high levels of iron because it is present in heme-, iron-sulfur and other proteins that function in respiratory and photosynthetic energy transducing. Correspondingly, in eukaryotic cells, the mitochondrion is a major ironutilizing compartment. It is well established that iron is transported into mitochondria for heme synthesis and ironsulfur cluster assembly. This is required for the formation of a functional respiratory electron transport machinery (11). Therefore, mitochondrial metabolism in mammals, fungi and plants is significantly affected under iron deficiency, as demonstrated by a number of studies (12)(13)(14). In plants, the chloroplasts are a primary target of iron deficiency. Changes in chloroplast structure, photosynthetic capacity and the composition of thylakoid membranes have been described for plants deprived of iron (15)(16)(17)(18)(19)(20)(21).
Plants have devised various strategies for acquiring iron (22). Generally, iron deficiency leads to the activation of the iron uptake systems in photosynthetic organisms. For example, the accumulation of the ferroxidase, a component of the high affinity iron uptake system in C. reinhardtii, is very rapidly enhanced when iron becomes limiting (23). Inactivation of IRT1, the most prevalent Fe 2ϩ transporter in Arabidopsis thaliana leads to a dramatic iron deficiency that is reflected by chlorosis (24 -26). Despite the evolution of elaborate ironuptake mechanisms in plants, iron deficiency-induced chlorosis remains a major agricultural problem (27,28).
The global impact of iron deficiency on photosynthetic productivity has been also shown in vast ocean regions, which are severely limited for iron (29,30). Generally, one can conclude that photosynthesis in the oceans and on land can occur in environment where iron availability is restricted.
Photosystem I (PSI) is a prime target of iron deficiency as it contains 12 atoms of iron per core complex. In algae, the degradation of PSI is also linked to remodeling of PSI-associated light-harvesting antenna (LHCI) (31)(32)(33). Cyanobacteria respond to iron deficiency by degradation of light harvesting phycobilisomes (34) and induction of the "iron-stress-induced" gene isiA. The ISIA protein, which has significant sequence similarity with CP43, a chlorophyll a-binding protein of photosystem II (PSII; (35,36), forms a ring of 18 molecules around a PSI trimeric reaction center, as shown by electron microscopy (37,38). The overall reorganization of the PSI complex from 900 kDa into 1.7 MDa complex highlights the large adaptive nature of the cellular response to iron deficiency, which helps to optimize the architecture of the photosynthetic apparatus to conditions in which iron is a limiting factor.
The marine diatom Thalassiosira oceanica shows a remarkable retrenchment of cellular metabolism and remodeling of bioenergetic pathways in response to iron availability (39). Low iron triggers a reduction in the level of iron-rich photosynthetic proteins while iron-rich mitochondrial proteins are preserved. Furthermore, iron deprivation causes a remodeling of the photosynthetic machinery resulting in the adjustment of light energy use to an overall decline in the level of photosynthetic electron transport complexes (39). These responses, reported for green algae such as C. reinhardtii (31,40,41), are important for minimizing photo-oxidative stress and optimizing photosynthetic function. As observed for T. oceanica, under conditions of low iron availability (in the presence of organic carbon) a hierarchy of iron allocation responses in C. reinhardtii result in the down-regulation of iron-rich photosyn-thetic complexes while iron-rich mitochondrial complexes remain stable (41). Notably, under photoautotrophic and mixotrophic conditions C. reinhardtii displays distinct iron deprivation responses, suggesting that the cell's response to iron deficiency is also dependent on trophic conditions (7)(8)(9). Thus bioenergetics pathways are remodeled in response to iron availability as well as to the type of carbon source available. Moreover, recent data has indicated that the regulation of iron-induced remodeling of the photosynthetic apparatus is linked to energy metabolism. Depletion of Proton Gradient Regulation Like1 protein (PGRL1) in C. reinhardtii has revealed a decreased efficiency of cyclic electron transfer under low iron conditions resulting in higher vulnerability toward iron deprivation (42).
It was our aim to generate a more comprehensive picture of how the proteome of C. reinhardtii varies in response to low iron under distinct trophic conditions and how these changes compare with differences observed for cells grown under photoautotrophic and photoheterotrophic iron replete conditions. Quantitative proteomics in conjunction with a novel hierarchical clustering approach revealed information about the responses of C. reinhardtii to low iron conditions and the iron requirements of photoautotrophic growth. These analyses provide novel insights into the relationships between protein networks required for photosynthesis and iron deprivation-elicited stress responses; these studies are providing the knowledge required for modulating the level of available iron to improve the photosynthetic performance of plants (43,44).

EXPERIMENTAL PROCEDURES
Strains and Culture Conditions-For all experiments the cell wallless wild type strain CW15 was used. Cells were grown photoheterotrophically (PH) in Tris-Acetate-Phosphate (TAP) medium (45) and photoautotrophically (PA) in minimal medium (HSM) (45) to a cell density of 1-2 ϫ 10 6 cells/ml at a light intensity of 60 E/m 2 /s. PHand PA-grown cells were isotopically labeled with 15 NH 4 Cl (Cambridge Isotope Laboratories Inc.) in iron-replete conditions (18 M iron) and nonlabeled ( 14 NH 4 Cl) in iron-deficient conditions (0.1 M iron). The cells were treated as a "steady state" culture and harvested from a liquid culture after 5 days (PH) or 10 days (PA) of iron deficiency. For the comparison of PA to PH (both in 18 M iron), for swapping experiments PH-or PA-grown cells were isotopically labeled.
For data shown in Fig. 4, the adh1 mutant (46) and its wild type parental strain CC125 (cw ϩ , mt ϩ , nit Ϫ ) were used. The two strains were grown in TAP medium at a light intensity of 60 E m Ϫ2 s Ϫ1 until they reached a cell density of about 2-3 ϫ 10 6 cells/ml (0 time point). Cells were then pelleted at 4000g for 2 min, washed twice and re-suspended in iron-free TAP medium to a cell density 1 ϫ10 6 cells/ml. Samples for protein, cell counting, chlorophyll, and fluorescence emission spectroscopy at 77K were collected at day 1 and 2 after the transfer.
Whole Cells and Isolation of Chloroplasts and Mitochondria-Cells were pelleted at 2500 ϫ g (Beckmann Coulter J 20 XP) and chloroplasts and mitochondria were isolated as described in Terashima et al. (47).
Determination of the Protein Amount and Immunodetection-The determination of the protein amount was carried out using the Pierce® Protein assay kit (Thermo Scientific, Perbio Science GmbH, Bonn, Germany). The absorbance of the BCA complex was detected at 562 nm after incubation of the protein with BCA at 55°C for 30 min. Protein concentrations were calculated based on a BSA calibration curve. Immunodetection was carried out as described in Terashima et al. (47). Antibodies against ADH1 1 were purchased from Agrisera (ref #AS10 748); anti-NDA2 antibodies were a kind gift of Dr. Pierre Cardol, University of Liè ge (Belgium).
Protein Analysis and Liquid Chromatography-Tandem MS (LC-MS/ MS) Analysis of Proteins-Isolated chloroplasts and whole cell samples were mixed at an equal amount of labeled/unlabeled 1:1 protein.
For isolated mitochondria, concentrations of mitochondrial proteins in samples stemming from iron sufficient or iron deficient conditions were enumerated by Western blot analyses (not shown) using the COX2b antibody mentioned above (section above). To account for differences in enrichment of COX2b in the distinct preparations, isolated mitochondria were mixed in the ratio labeled/unlabeled 2.5:1 protein. The mixed samples were fractionated by SDS-PAGE. A different number of bands was excised and digested tryptically as described in Naumann et al. (40). A schematic overview of all mixed preparations and performed measurements is shown in Fig. 1. The LC-MS/MS measurements were carried out as described in (47) (49) with 32 bit precision. Standard proteomics data identifications were performed using Proteomatic (50) as a data processing pipeline tool providing a graphical user interface to the later mentioned tools and pymzML (51) for direct access to mzML files. mzML files associated with this manuscript may be downloaded from peptideatlas.org with the database tag 'Hoehner_2013Ј. Files in mgf format were used as input for database search algorithms. Mgf files were generated using pymzML (51).
Enhancing Proteomics Data-The detailed procedure of enhancing proteomics data by means of LC-MS/MS alignment and quantifications based on identifications performed in different LC-MS/MS runs within a well-defined elution window will be described in more detail elsewhere (Barth & Fufezan, manuscript in prep.). Briefly, all LC-MS/MS runs were aligned using a combination of the qualities of peptide identifications and peptide quantifications. Peptide identifications were performed by OMSSA (version 2.1.9) (52) and X! Tandem (2010.12.01) (53). Identification parameters were used as described in (47) with the following changes. Precursor mass tolerance was set to 5 ppm and oxidations of methionine and tryptophan (single and double) were set as variable modifications. For OMSSA a linear precursor charge dependence and enzyme specificity was used. For X!Tandem noise suppression was enabled. Protein databases JGI4.3 AUGUSTUS10.2, NCBI BK000554.2, and NCBI NC_001638.1 were used. The resulting database contained 17195 nonredundant protein sequences, with additional 1198 common contaminants. Decoy protein sequences were generated by randomly shuffling tryptic peptides while retaining the redundancy of nonproteotypic peptides. Qvality (version 2.02, (54)) was used for statistical validation of peptide identifications using a posterior error probability (PEP) threshold of less than 0.05. Quantifications were performed using pyQMS (Niehues & Fufezan, manuscript in prep). Preliminary data were quantified using qTrace (47). Ratios obtained by both quantification algorithms were very similar.
Isotopologs, their elution profiles and their elution windows could be assigned to a given peptide. Confident quantifications of the peptides were possible by retaining the strict rule that a peptide was only considered if it could be identified within the aligned retention time window allowing a deviation of Ϯ 30 s. All data analysis were performed using piqDB (Fufezan and coworkers, unpublished), which is very similar to SuperHirn (55) incorporating quality of the identifications and quantifications of a given peptide into the definition of its accurate mass time tag, that is in our case into the definition of the retention time windows.
with ratios based on few MS1 scans are weighted less. Protein ratios and standard deviations of all proteins in all experiments (log 2 ) are shown in supplemental Table S1.
Hierarchical Clustering of Proteomics Data-Enhanced proteomics data was clustered using the Python (56) and Python modules numpy (57) and scipy (58) combined in the Python module pyGCluster (Jaeger et al., manuscript in prep). Briefly, agglomerative hierarchical clustering for different distance-linkage combinations was performed using a pool of (1) six methods to define the distance among clusters (linkage): "complete," "average," "weighted," "centroid," "median," "ward," and (2) two functions to calculate the distance between clusters (metric): "Euclidean" and "correlation." Because the ward, centroid and median linkage methods cannot be used with the correlation function, a total of nine distance-linkage combinations were evaluated. Data were clustered independently 250,000 times. During each iteration step, each ratio was altered by randomly sampling a ratio from a Gaussian distribution described by the original ratio mean and its standard deviation. This step injects a high degree of noise into the data set, ensures stringent hierarchical clustering and allows discriminating between core elements of a well-defined high frequency cluster in, e.g. Euclidian space and the adjacent less frequent elements. A cluster is a group of genes or proteins that are similar given the distance definition. For example using Euclidian or correlation distance will group proteins by similarity in ratios over all conditions or by similar trends in ratios comparing the conditions, respectively. Because ratios, i.e. (0, . . . ,ϱ) are not well suited to estimate distances, the protein ratios were log 2 transformed before the clustering.
Obviously, creating a "consensus tree" is neither feasible nor informative using such an approach. The overall frequency in which a certain cluster, i.e. a certain set of proteins, is found is more informative because its position, i.e. clade, can occur on different levels of the tree. In this approach, only clusters with a frequency higher than 0.5% in at least one of the nine hierarchical clustering assays were retained for further analysis. The 0.5% threshold was chosen for the sake of simplicity, i.e. only the top 490 clusters were analyzed further and grouped into communities.
Clusters, their frequency, their relation to each other and the communities that they form are illustrated as a node map using NetworkX (59) and gephi (60). Clusters are shown as nodes and their relations are shown as edges connecting related nodes, i.e. those sharing same proteins. Therefore, the node map shows the high frequency clusters from a combination of multiple hierarchical clustering assays using high number repetitions with a standard deviation weighted degree of variation in each dimension (condition) of each object (protein) during each iteration, and the relation among those high frequency clusters. The communities that are found represent welldefined groups of proteins that show cotranslational regulation. Such node maps are referred to as cNode maps.
Changes in protein levels were plotted for each community as a heat map using the Python scripting language and pyGCluster (Jaeger et al., manuscript in prep). The log 2 transformed protein ratios are color-coded and the standard deviation is illustrated by the size of the plotted box in the heat map (see text). The cFrequency is the sum of all cluster frequencies in which the given protein is within the whole community. The cluster frequencies themselves were calculated as the median of the frequencies for all nine assays. The cFrequency illustrates therefore the major elements of each community and the relative frequency among communities. For more detailed information on frequencies for each clustering assay, see supplemental Table S2, and for more details concerning pyGCluster, see Jaeger et al. (2013, manuscript in prep).
Fluorescence Emission Spectroscopy at 77K-The samples were used at a chlorophyll concentration of 5 g/ml, as described in (42).
Fluorescence emission spectra were recorded as an average of four scans within the window of 600 -800 nm at a scan speed of 500 nm/min and a slit width of 3 nm. Excitation of the sample was at 435 nm.
Oxygen Exchanges-Oxygen uptake and evolution rates were measured as described in (42). The measurements were performed at 1000 E/m 2 /s. After the measurement, 500 l of the sample was removed from the cuvette for determining cell number to allow for the calculation of pM O 2 /cell.

RESULTS
Iron deficiency responses in C. reinhardtii are influenced by the trophic status of the cells (7)(8)(9)41). Yet, although many aspects of the cellular iron deficiency response are driven by low iron availability, other properties appear to be governed by the cell's metabolic status. To dissect these responses and obtain more detailed insights into these regulatory processes we performed an in-depth quantitative proteomics approach using metabolic labeling. To this end iron sufficient and deficient cells were grown under photoheterotrophic (PH) or -autotrophic (PA) growth conditions in 15 N or 14 N labeled media, respectively. The iron deficiency status was verified using low temperature fluorescence emission spectroscopy and immunoblotting using antibodies directed against marker proteins, diminished in response to iron deprivation (see supplemental Fig. S1). Moreover PH grown iron sufficient cells were labeled with 15 N and compared with 14 N labeled cells grown under PA conditions. Here also a swapping experiment was performed to detect and avoid isotope effects. Quantitative analysis showed very good correlation among the label swap experiments (see supplemental Fig. S2). 15 N labeled iron sufficient cells or isolated chloroplasts were mixed on equal protein basis with the respective complement from 14 N labeled iron deficient or iron sufficient material (Fig. 1). In addition mitochondria were isolated from 15 N labeled iron sufficient and 14 N labeled iron deficient cells grown under PH conditions and also mixed on equal protein basis (Fig. 1). Mixed samples were fractionated by SDS-PAGE and protein bands were excised from the gels. On average about 56 bands were cut per gel. SDS-PAGE bands were in-gel digested with trypsin and peptides were analyzed by LC-MS/MS. We completed and processed two independent biological replicates. For the whole cell and chloroplast sample sets we investigated 337 and 523 excised and tryptically digested protein bands, respectively. Moreover, for isolated mitochondria in total 124 excised and proteolytically digested protein bands were examined by LC-MS/MS. The mass spectrometric analyses resulted in the identification of peptides and proteins and subsequently allowed the quantitation of protein ratios. The results stemming from the two independent biological replicates are summarized in Fig. 2 and are plotted as number of quantified proteins and their respective proteins ratios starting from the lowest to the highest observed ratios. As shown in Fig. 2, from the six independent experimental set-ups at minimum 1316 and at maximum 1944 proteins could be quanti-fied. These high numbers can be attributed to the enhancement procedure using piqDB, (Fufezan and coworkers, unpublished), which is similar to SuperHirn (55). Consequently, such large amount of protein ratios provided a good basis for further exploration. Thus coregulation analyses of the proteomics data was performed using a novel and advanced hierarchical clustering approach, that injects noise into the data set and uses a large number of iterations, using pyGCluster (Jaeger et al., manuscript in prep). Briefly, pyGCluster uses the uncertainty of the data at hand by including the standard deviation (std) into the agglomerative hierarchical clustering approach. For each iteration, a new data set is generated in which each data point is randomly sampled out of the Gauss distribution described by its mean and std. The generated data set is described by the original data yet it will cluster differently. pyGCluster then applies nine different distance-linkage combinations for each of those perturbed data sets to ultimately, after a high number of iterations, isolate the high frequency clusters and merge those into coregulated communities. The advantages lie in less user bias and most importantly in reducing the data set at hand to very few proteins that show coregulations although heavily perturbed based on the uncertainty of their data.
Quantitative mass spectrometric data from whole cell extracts (WC) and isolated chloroplasts (C) have been clustered using pyGCluster (see Material and Methods). It is expected that chloroplast localized proteins harbor similar trends in the six settings compared, which is indeed observed for all chloroplast proteins found in the clusters. Notably among the 109 proteins clustered, 53 were reported to be part of a chloroplast proteome (61) (supplemental Table S3). For nonchloroplast proteins, tendencies seen in whole cells should be representative, whereas changes observed in enriched chloroplast fractions might be biased by the organelle enrichment and purity, yet, jointed clustering might be reached because of these properties. The results are illustrated as cNode maps (Fig. 3) that highlight well-defined high frequency clusters of proteins, their relation to each other and how they are grouped into communities. These communities are colored coded and labeled with roman numerals in the cNode map. A total of 23 communities have been identified. The proteins within each community are illustrated as heat maps, a selection of heat maps is shown in

FIG. 1. Schematic overview of the performed measurements for protein identification and quantification via stable 15 N isotopic labeling.
Whole cells grown in iron deficiency ( 14 N (ϪFe)) and iron sufficiency ( 15 N (ϩFe)) were harvested from photoheterotrophic conditions (PH) and photoautotrophic conditions (PA). Mixed samples of ( 15 N (ϩFe)) and ( 14 N (ϪFe)), which were separated via one-dimensional SDS-PAGE, were cut into 56 bands for PH and 54 and 59 bands for PA for two biological replicates, respectively. Isolated chloroplasts of a mix of ( 15 N (ϩFe)) and ( 14 N (ϪFe)) samples resulted in 54 and 56 bands for PH, 52 bands for PA and for isolated mitochondria in 62 bands for PH for two biological samples, respectively. The mix of samples from 14 N PA and 15 N PH in iron-sufficient conditions resulted in 56 bands for whole cells and 55 bands for isolated chloroplasts for two biological samples, respectively. For whole cells from iron-sufficient conditions also a label swap experiment was conducted ( 15 N PA and 14 N PH). Mitochondria could not be isolated from PA of ( 15 N (ϩFe)) and ( 14 N (ϪFe)). Therefore no measurements were carried out for isolated mitochondria from PA under these conditions as well as for the comparison of 14 N PA to 15 Table T2. The cFrequency of a protein is equal to the sum of all cluster frequencies the protein is part of, within its community. Heat maps are subdivided to illustrate the separation of subcommunities visible as branches in the cNode map. Fig. 3. All heat maps can be found in the supplemental Fig. S3 (communities 0 -XXII). The heat maps also show the cFrequencies of each protein within a community. The cFrequencies are the sum of all cluster frequencies the protein is part of within a community. Because each cluster is found with different frequencies by the 9 hierarchical clustering assays performed here, the median frequency of those was used to represent each cluster. The colors in the heat maps between purple and yellow represent reduced and increased protein levels, respectively, comparing iron deficient (ϪFe) and iron sufficient (ϩFe) conditions or PA and PH settings. Ratios (i.e. -Fe/ϩFe or PA/PH) are log 2 transformed and the uncertainty, i.e. the standard deviation of a ratio is reflected in its box size. The higher the standard deviation and therefore the higher the uncertainty of a given ratio, the smaller is the box. This is illustrated in the legend.
Subunits of the cyt b 6 f Complex and Photosystem I As Well As Proteins Involved in Chlorophyll Biosynthesis Exhibit a Uniform Down-regulation Under Iron Deprivation-Community III (orange, Fig. 3 top right) contains 17 proteins, which are all down-regulated under low iron. This cluster is dominated by four subunits of the cytochrome b 6 f complex (cyt b 6 f complex, cytf, Rieske, cytb 6 and subunit IV) and five enzymes involved in chlorophyll biosynthesis (the geranyl-geranylreductase (GGR), uroporphyrinogen-III decarboxylase (UROD1), light-dependent protochlorophyllide reductase (POR), copper target homolog 1 (CTH1), and porphobilinogen deaminase (PBGD1)). This cluster supports earlier notions that cyt b 6 f complex and the chlorophyll biosynthesis pathway are prime targets of iron-deficiency (31,(62)(63)(64)(65). A light-harvesting-like protein LIL (Lhl3) shows the same tendency as the chlorophyll biosynthesis proteins. Another photosynthetic protein complex that is preferentially down-regulated under iron-deprivation is photosystem I. In community IV (yellow, Fig. 3 left bottom) seven out of eight proteins belong to the photosystem I complex (PsaA, B, D, E, F, L, H (Cre07.g330250)). The only nonphotosystem I protein in this cluster is a protein of unknown function (Cre05.g244500). All photosystem I complex subunits remain unchanged or slightly increase when iron replete PA/PH conditions are compared.
Core subunits of photosystem II (PSII) and the ATPase and many light-harvesting proteins are not among the clustered proteins. The polypeptides belonging to PSII and ATPase multiprotein complexes as well as to the light-harvesting complement show only moderate down-regulation under low iron availability (supplemental Fig. S4A). Moreover it appears that the impact on these complexes under low iron is generally more pronounced under PA rather than under PH conditions. In contrast, light-harvesting proteins that are associated to PSII, with the exception of Lhcb4 and Lhcb5, and listed in supplemental Fig. S4A are not reduced in amounts by iron deficiency. These data correlate very well with published results (41) and underpin the differential regulation of photosynthetic complex proteins in response to iron deprivation.
Clustering Reveals Response Networks of Acclimation to Iron Deficiency-A complete different picture can be seen in community VI (green, Fig. 3 top right). Here proteins are increased in expression under iron deprivation and mainly decreased under PA/PH iron replete conditions. The protein with the highest cFrequency belongs to the reticulon protein family (Cre06.g308950.t1.1), such proteins are found to be associated with the endoplasmic reticulum (ER) (66). The Chlamydomonas protein possesses a characteristic dilysine motif at the C-terminal and may be associated with ER channel-like membrane complexes as found for the reticulon NSP-A that is connected with an ER Ca 2ϩ -ATPase (66). Other members of this community are a 6-phosphogluconate dehydrogenase (GND1a, Cre12.g526800.t1.1), a putative glutathione S-transferase (Cre12.g538100.t1.1) and a dynamin related GTPase (DRP1), followed by two proteases, ASP, an aspartic-type endopeptidase and CEP (a candidate serine carboxypeptidase) and a putative mitochondrial phosphate carrier 1 (MPC1). Lower cFrequency members are a septin (SEP1), a granule bound starch synthase I (STA2), GSH1, a gamma-glutamylcysteine synthetase, a flagellar associated protein (AGG3/AGG4; star behind the name tag in Fig. 3 indicates a merge tag with models Cre10.g456100.t1.1/ Cre10.g456050.t1.1 because of non proteotypic peptides matching on both protein sequences), Cre12.g547300.t1.2 (a putative ABC transporter), as well as two proteins of unknown function (Cre08.g382950.t1.1 and Cre01.g010400.t1.1).
Overall, members of community VI are involved in bioenergetic remodeling and stress responses as many metabolic enzymes, proteases, and redox related enzymes dominate in this assembly. Importantly the clustering reveals that certain responses are distinct among the different conditions analyzed, a scenario that is also observed in community XI.
Three proteins are represented in an adjacent subcommunity. These proteins exhibit a distinctive up-regulation under PA/PH (ϩFe) and a diminishment under iron-deficiency. The proteins present in this cluster are a stomatin-like protein (Cre33.g782550.t1.1), a NADP malic enzyme (MME5) as well as a protein that contains a thioredoxin domain and four putative EF-hand motifs (Cre03.g202950.t1.2). According to their regulation, these proteins might be important for phototrophic growth.
Proteins of another subcommunity show a uniform downregulation under PA growth conditions (ϪFe/ϩFe) and exhibit a slight up-regulation in the other two growth settings. Out of seven proteins present in this community, three proteins belong to the glycine cleavage multienzyme complex (GCSH, GCST and GCSP). This complex has glycine decarboxylase activity and is required for efficient photorespiration as demonstrated by analysis of an Arabidopsis thaliana mutant devoid of glycine decarboxylase activity (70). Another photorespiratory enzyme in this cluster is the serine glyoxylate aminotransferase (Cre01.g005150.t1.1) (71). Other members of this assembly are a N-acetyl-gamma-glutamyl-phosphate reductase (ARGC1), the large subunit of carbonyl phosphate synthase (Cre28.g776100.t1.1), and a putative 10-formyltetrahydrofolate synthetase (FTHFS) (Cre13.g566000.t1.1). The latter enzyme is involved in formate dependent synthesis of 10-formyltetrahydrofolate, a precursor for 5,10-methylene tetrahydrofolate that is required for synthesis of serine and thus also linked to the photorespiratory process (72). 5,10 methylene tetrahydrofolate is also produced by glycine decarboxylase activity which is then used as the main single carbon donor for synthesis of serine (73). Yet, the synthesis via the 10-formyltetrahydrofolate synthetase provides independent 5,10 methylene tetrahydrofolate for the synthesis of serine (74). Interestingly using 13 C nuclear magnetic resonance, interactions among serine hydroxymethyltransferase, tetrahydrofolate synthase, and glycine decarboxylase complex activities were shown (74). Accordingly our finding of coclustering of the glycine cleavage multienzyme complex, theserineglyoxylateaminotransferaseand10-formyltetrahydrofolate synthetase supports the assumption of coordinated action in serine biosynthesis and reveal that the process of photorespiration is active under PH (ϪFe/ϩFe) and PA/PH (ϩFe), whereas it appears to be not active under PA (ϪFe/ϩFe).
Respiratory Complexes are More Severely Down-regulated Under Autotrophic Low Iron Conditions-The community XIII (blue, Fig. 3 bottom middle) consists of 12 proteins. Among those, seven proteins are subunits of mitochondrial multiprotein respiratory protein complexes (three subunits of complex I, three subunits of complex III, and one subunit of complex IV). Considering the whole cell data, subunits of the respiratory complexes are slightly induced under PH (ϪFe/ϩFe) and PA/PH (ϩFe), but down-regulated under PA (ϪFe/ϩFe). Looking at the entire list of mitochondrial respiratory proteins organized in complexes I to IV as well as in ATPase in whole cells or isolated mitochondria (supplemental Fig. S4B), it appears that, except for complex II proteins and subunits of the mitochondrial ATPase, protein expression has a stronger impact under PA (ϪFe/ϩFe) or PA/PH (ϩFe) than under iron deficient PH growth condition. The rather minor down-regulation of respiratory complexes under PH low iron availability has been already described (41) and is confirmed by the current study. However a new aspect is that autotrophic growth under low iron induces a remodeling of the respiratory machinery.
Deletion of ADH1 Attenuates the Remodeling of the Photosynthetic Machinery in Response to Iron Deficiency and Stimulates Expression of Stress-related Proteins-To substantiate our quantitative proteomic data, we investigated the up-regulation of ADH1 by immunochemical methods and analyzed the iron-deficiency response in an ADH1 deletion mutant that has been recently reported (46). To analyze the induction of ADH1 under low iron conditions and possible acclimation phenotype of the ADH1 deletion mutant in response to the iron-deficiency response we shifted wild type and mutants grown from iron-sufficient acetate-containing medium for 2 days into iron-free culture solution. Such a shift causes a pronounced remodeling of the photosynthetic complexes leading to functional uncoupling of PSI associated LHCI with PSI and diminishment of PSI (31). After shifting to iron-free medium, growth performance of wild type and ADH1-deficient cells was comparable (4.53 Ϯ 0.46 ϫ 10 6 cells/ml for WT versus 3.93 Ϯ 0.25 ϫ 10 6 cells/ml for adh1 after 2 days from the shift; Fig. 4A, upper panel). Interestingly, the decrease in chlorophyll per cell content associated with the iron-starvation response appeared to be less pronounced in the mutant devoid of ADH1 (Fig. 4A, lower panel). Low temperature fluorescence spectroscopy of wild type cells after 2 days of iron-deficient growth showed a maximum emission at 705 nm whereas iron sufficient cells peaked at 714 nm (Fig. 4B), revealing that PSI and LHCI are functionally uncoupled under low iron as described (31). Interestingly, this shift in maximal fluorescence emission is not observed in adh1 mutant cells grown in iron-free medium for 2 days. To obtain further insights into the acclimation to iron-deficiency, cells were harvested from the outlined conditions and fractionated by SDS-PAGE, immunoblotted, and analyzed by a set of different antibodies (Figs. 4C, 4D). Analysis of ATPb protein levels by anti-ATPb antibodies was used as a loading control. As expected the immunoblot data revealed for wild type a downregulation of PSI and PSII as visualized by decreases of PSAD and PSBA in day 1 and 2 of iron deficiency, respectively. On contrary this down-regulation is less pronounced in the ADH1 deletion mutant, particularly, PsbA appears to be two-times more abundant after 2 days of iron-deficiency as compared with the wild type (Fig. 4D). Additionally down-regulation of LHCA3 is less pronounced in adh1 (Fig. 4C), which likely explains the minor blue shift of the low temperature fluorescence emission spectrum in comparison to wild type, as degradation of LHCA3 correlates with changes in low temperature fluorescence emission of PSI-LHCI (40). Interestingly, degradation of FDX1 on the other hand is paralleled in wild type and adh1 mutant. In contrast, the iron-deficiency induced expression of NDA2, LHCSR3 and PGRL1 is significantly stimulated in adh1 (Fig. 4D). Notably, iron-deficiency stimulated expression of the three proteins is also observed in the quantitative proteomics data (supplemental Fig. S4C), thereby confirming the mass spectrometric data. DISCUSSION In this work we compared iron deficiency responses under different metabolic conditions with data obtained from iron sufficient cells grown under PA and PH settings. Taking advantage of quantitative mass spectrometric data stemming from these distinct conditions, protein co-expression and the resulting network analyses provided new insights into coordinated remodeling of bioenergetics pathways and iron deficiency responses in C. reinhardtii and revealed interdependencies with the cellular metabolic status. Analyses of an adh1 knockout mutant suggest that the low iron regulatory acclimation networks are particularly sensitive to the chloroplast metabolic and/or redox status.
Iron Deficiency Uniformly Down-regulates Assemblies of Proteins Connected to PSI, the cyt b 6 f Complex and Chlorophyll Biosynthesis-A hallmark of the iron deficiency response in cyanobacteria, algae and higher plants is the down-regulation of PSI (see above). This is also reflected in community IV shown in Fig. 3, where most of its members are PSI subunits. Moreover, the down-regulation of four subunits belonging to the cyt b 6 f complex can be observed in community III (Fig. 3). Down-regulation of the cyt b 6 f complex under iron deprivation is expected (31). Nevertheless, coclustering of subunits belonging to a common multiprotein complex as already seen for PSI subunits underscores the robustness of the quantitative data and the strength of the clustering algo- rithm, pyGCluster. A protein that clusters with the subunits of the cyt b 6 f complex is CTH1. CTH1 is a diiron protein that functions as the aerobic cyclase in chlorophyll biosynthesis and is another key chloroplast target under iron deficiency (31,63,64,75,76). The diminishment of CTH1 under iron deprivation is likely responsible for the accumulation of Mgprotoporphyrin IX and Mg-protoporphyrin IX monomethyl ester under such conditions (62,63). The regulation of CTH1 goes in line with the down-regulation of GGR, POR, PBGD1, and UROD1 that are also related to the tetrapyrole biosynthesis process. The impact in chlorophyll biosynthesis, however, has a very distinct effect on chlorophyll binding proteins and protein complexes as discussed above. Remarkably, LHCSR3 is upregulated under PH iron deficiency ( Fig. 4D; supplemental Fig.  S4C) as well as other LHCII polypeptides (supplemental Fig.  S4A) despite a decrease in chlorophyll biosynthesis. Thus the alteration of the abundance of chlorophyll binding proteins is more likely to be a regulatory mechanism than an answer to chlorophyll deficiency under iron deprivation in agreement with previous conclusions (31,41).

Remodeling of Bioenergetics Pathways Under Iron-deficiency is Modulated by the Cellular Metabolic Status-As
found in an earlier study (41) mitochondrial protein complexes are only marginally affected by iron deficiency in C. reinhardtii, although this impact is stronger under PA than under PH growth conditions, which is also reflected in respiratory activity (Table I, supplemental Fig. S4B). Respiratory complex I itself contains at least 26 atoms of iron (77), which is more than the entire photosynthetic machinery possesses, yet this complex as well as the other iron rich respiratory electron transfer complexes remain unaffected by iron deficiency (supplemental Fig. S4B). Down-regulation of the photosynthetic electron transfer complexes PSI as well as cyt b 6 f is therefore a result of a coordinated regulation. Notably, a bifunctional alcohol and acetaldehyde dehydrogenase (ADH1) is induced under iron deprivation (Fig. 3, Fig. 4C). Interestingly expression of ADH1 follows a day/night cycle and was mainly insensitive to oxygen availability ((78); Barth and Fufezan, unpublished results). A proteomic study found the chloroplast localized ADH1 to be 1.5 fold up-regulated under anaerobic conditions (47). Moreover, ADH1 is also induced under copper and zinc deficiency (79). As discussed for anoxic growth in C. reinhardtii (46), ADH1 could oxidize two molecules of NAD(P)H per acetyl-CoA under low iron availability and thereby modulate the redox poise of the chloroplast by re-plenishing it with NAD(P) ϩ . The induction of cyclic photosynthetic electron flow (CEF) (42) would likewise decrease the flow of electrons into the NADPH pool, acidify the thylakoid lumen and activate energy dependent nonphotochemical quenching, going along with the described induction of LHCSR3 (80), particularly under PH settings of low iron availability (8,41) ; Fig. 4D). Most interestingly deletion of ADH1 leads to an increased expression of LHCSR3, PGRL1 and NDA2 after the onset of iron deprivation as compared with wild type (Fig. 4D), indicating that the redox poise of the chloroplast in adh1 knockout strain is more reducing than in the wild type. Beside increase of PGRL1 also induction of NDA2 can be seen as a response to diminish the chloroplast redox poise, as NDA2 can lower the stromal redox poise by oxidizing NADPH and reducing plastoquinone (81). The attenuated degradation of PSI and PSII and the absence of the low iron induced low temperature fluorescence emission shift of PSI-LHCI in the ADH1 deletion mutant could be explained by the fact that stress protection in form of LHCSR3, PGRL1 and NDA2 is faster induced and therefore more protective. Alternatively, iron deficiency induced proteases could be less active and/or less induced in the adh1 mutant. On the other hand, iron deficiency induced degradation of FDX1 is unaltered when wild type and adh1 are compared (Fig. 4D), indicating that the iron deficiency response is active in the mutant. Summarizing these considerations, we can at least conclude that the deletion of ADH1 has a profound impact on the remodeling of the photosynthetic apparatus as well as in the iron dependent regulation of LHCSR3, PGRL1 and NDA2 expression. Similarly, the activation of a metabolic by-pass resulting in glycerol synthesis had been shown for the adh1 mutant placed under dark anoxic conditions. This new metabolic pathway has been suggested to compensate for the inability of adh1 to synthesize ethanol and concomitantly reoxidize NAD(P)H (46). ADH1 clustered with other metabolic enzymes namely 6-phosphogluconate dehydrogenase (GND1a) and glyceraldehyde-3-phosphate dehydrogenase (GAP1, plastidic) and shared up-regulation under low iron and down-regulation by PA/PH under iron sufficiency (Fig. 3). GND1a catalyzes the second step of the oxidative Pentose-Phosphate-Pathway (OPPP). The glycolytic enzyme GAP1 catalyzes the first step of the second phase of glycolysis and was found in vascular plants to be induced in expression under anaerobiosis and heat stress (82) as well as under iron-deficiency (83). Impaired iron availability clearly decreases photosynthetic activity, so that plants as well as algae probably enter into an anaerobic growth mode. The increased expression of GAP1 probably mirrors increased request to synthesize ATP under conditions where ATP synthesis is diminished because of decreased photosynthetic and respiratory activities. In line, it was found that specific activities of GAP and GND1 were greater in leaf tissues unable to generate reducing equivalents and ATP by photosynthesis (84). Thus induction of GAP1 and GND1 in C. reinhardtii under low iron likely reflects the necessity to generate NADH, NADPH, and ATP from the gluconeogenesis and successive oxidation of glucose. Under such circumstances ADH1 is required to oxidize NAD(P)H to allow advancement of glucose oxidizing by replenishing NAD(P) ϩ . Its deletion, under conditions when GAP1 and GND1 are induced, on the other hand will therefore certainly increase the cellular redox poise, supporting the conclusions made above. As ADH1, GAP1, and GND1a are chloroplast localized (47), it is likely that the redox poise of the chloroplast is involved in communication with the nucleus to modulate gene expression. Interestingly, transcripts of genes encoding ADH1, GAP1, and GND1a are down-regulated by iron-deficiency (9), whereas the protein amounts are up-regulated (herein; (9)), thus strengthening our interpretation that these proteins are in the same regulatory circuit. Moreover, the divergence of transcript and protein expression may implicate post-transcriptional regulation for expression control. It is tempting to speculate that this may occur via the cytosolic translational machinery and the redox poise of the chloroplast. Besides remodeling of bioenergetics pathways, the low iron response also induced proteases, redox-and oxidative stress related enzymes. The Iron Deficiency Response Induces Proteases, Redoxand Oxidative Stress Related Enzymes -Nonetheless Differences Between PH and PA Conditions are Evident-In the iron deficiency response of C. reinhardtii several proteases are up-regulated. This up-regulation is especially pronounced under PH conditions (Fig. 3), as seen for a putative metallopeptidase (Cre01.g023350.t1.2) and a serine-type protease (Cre10.g423300. t1.1) as well as a peptidase M (Cre20.g758550.t1.2, supplemental Fig. S3 community XVIII). Two other proteases, a cysteine protease, CEP, and an aspartic acid protease, ASP, clustered together and are also induced under PA iron-deficiency but down-regulated under PA/PH iron sufficient conditions. Notably, CEP and ASP proteases are chloroplast-localized (47). Generally plant proteases are involved in many cellular functions, including chloroplast related properties (85). However, the precise functions of the four proteases found in our protein assemblies are unknown. For strain C. reinhardtii W80, a cysteine protease was isolated that was also induced under photo-oxidative stress conditions (86). Photo-oxidative stress in the chloroplast of iron-deficient photoheterotrophic cells can be conceived, as PSI and the cytochrome b 6 f complex are down-regulated, whereas the PSII related light-harvesting protein complement remains stable or either increases or decreases in amount as in the case of LHCBM1 and LHCBM6/3/7 or LHCBM4, respectively (supplemental Fig. S3 community XVII). It was reported that under PH growth and iron deprivation the light-harvesting antenna was still functionally coupled to PSII, which caused significant photoinhibition (41). Thus induction of the proteases could be related to the light-induced photo-oxidative stress that is particular for iron deficient PH grown cells and potentially involved in the degradation of Fe-containing proteins (e.g. FDX1, FAB2, cytf etc.) At the same time respiratory activity of these cells, as measured by oxygen uptake activity in the dark, is only 1.5 times lower as compared with iron sufficient PH grown cells (Table I). In contrast, respiratory activity of iron deficient PA grown cells is four times lower in regard to iron sufficient cells. Importantly the oxygen evolution activity of iron-sufficient PA grown cells is increased almost twofold compared with the activity of PH grown cells, whereas iron deficient PA grown cells produce less oxygen as compared with iron deficient PH cells (Table I). Thus iron deficient PH grown cells possess a higher photosynthetic activity as well as a higher respiratory rate than iron deficient PA grown cells. Apparently, PH grown cells suffer a higher degree of photo-oxidative stress for their pronounced bioenergetic capacity. Notably, we found proteins with a flavodoxin domain (AGG3 (Cre10.g456000.t1.1) and particularly AGG4 (Cre10.g456050.t1.1)) that are highly up-regulated under PH low iron conditions (Fig. 3, supplemental Fig. S3 community XVIII). AGG3 and AGG4 were found to be enriched in the flagella proteome and the membranematrix fraction of C. reinhardtii (87). AGG3 was shown to be involved in mediating the orientation of the alga toward a directional light source (88). It was concluded that downregulation of AGG3 blocks positive phototaxis (88). It is unclear what impact overexpression of AGG3 would have, also the concrete function of AGG4, a close paralog of AGG3, is currently unknown but it is tempting to speculate that increase of AGG3 and AGG4 under iron deficiency impacts the phototaxis of the alga, thereby contributing to the overall strategy to minimize photooxidative stress and optimize photosynthetic function. Interestingly, under iron deficiency, flavodoxin substitutes for ferredoxin as electron acceptor of PSI in diatoms and cyanobacteria (89 -92). However, because of their localization and low sequence similarity to photosynthetic flavodoxins, AGG3 and AGG4 are rather involved in regulation of phototaxis than in photosynthesis.
Having acetate as a fueling source, the iron deprived PH grown cells produce energy via respiration and apparently invest in a costly stress response to maintain photosynthetic activity. In line with this argumentation, a Peptidyl-prolyl cistrans isomerases (CYN37), a translation factor for expression of the chloroplast-encoded psbA gene (TBA1) (67), a cytochrome c peroxidase (CCPR1) and a thioredoxin o (TRXo) were induced solely under PH iron-deficiency and clustered together (Fig. 3, XI). Notably these proteins are also induced under PA iron plus conditions. CYN37 shares highest similar-ity with Arabidopsis TLP38, a luminal chloroplast protein (93). PPIases catalyze the cis-trans isomerization of the peptidylprolyl bond during protein folding, a rate-limiting step in this process (94). Interestingly, Vener and colleagues (95) found redox-dependent activity of chloroplast PPIase implying a link to the redox poise of the stroma and suggesting that protein folding catalysis under conditions of active photosynthesis is important. The finding of a constant up-regulation of CYN37 under iron deficiency stress or PA conditions underpins this notion. TBA1 encodes a putative oxidoreductase required for translation of the chloroplast psbA mRNA (67). Up-regulation of chloroplast localized TBA1 may reflect high turnover of the D1 protein, thereby requiring increased translation of psbA RNA. The oxidoreductase domain may connect TBA1 activity to the redox poise of the stroma and implicate a regulation similar to CYN37. The mitochondrial cytochrome c peroxidase (CCPR1) is located at the intermembrane space, reduces hydrogen peroxide to water using the electrons provided by cytochrome c and is thereby involved in scavenging of reactive oxygen species in the mitochondria (96). The mitochondrial thioredoxin o (TRXo) that is found in the same cluster is another player in this scavenging system. It acts as an electron donor to mitochondrial type II peroxiredoxin F, which has been shown to be essential for redox homeostasis and root growth of A. thaliana under stress (97). Correspondingly, the mitochondrial manganese containing superoxide dismutases (MSD1 and MSD2) are strongly induced under PH low iron, likely combating superoxide formed by the respiratory electron transport (supplemental Fig. S3 community XVIII, supplemental Fig. S4D). Beside these responses, Fig. 3 reveals other redox and oxidative stress related proteins that are up-regulated under PH and PA iron-deficiency, though up-regulation of these proteins is stronger under PH low iron. These proteins are represented by a putative glutathione Stransferase (Cre12.g526800.t1.1) and GSH1, the gamma-glutamylcysteine synthetase. This enzyme catalyzes the second step of glutathione (GSH) biosynthesis by adding glycine to the C-terminal site of ␥-glutamylcysteine to form GSH (98). GSH operates as a redox buffer but contributes also to other cellular functions; in particular it participates in heavy metal detoxification and in ROS scavenging (99). GSH is involved in the generation of reduced ascorbate, another cellular redox buffer (100) and antioxidant (101), which is required e.g. for the detoxification of ROS produced by PSI (102,103). Likewise monodehydroascorbate reductase (MDAR), which also participates in regeneration of ascorbate (104) by reducing monodehydroascorbate using NADPH, is up-regulated under PH low iron growth condition (Fig. 3). The putative glutathione S-transferase (Cre12.g526800.t1.1) might be functionally related to the glutathione S-transferase, GST1. GST1 expression is induced as an acclimation response to singlet oxygen treatment and its overexpression was sufficient to enforce resistance against singlet oxygen stress (105). Ledford and colleagues (105) discuss that GST1 could function as lipid peroxidase, thereby increasing lipid peroxidase activity of the cells, which in turn would protect against reactive oxygen species (ROS) and/or participate in signaling activity.
Interestingly another glutathione S-transferase (Cre01. g064400.t1.1) and a DegP protease (DEG5A) are induced under PA iron-deprivation (Fig. 3), indicating functional relation between these two but on the contrary distinct tasks for the one induced under PH iron deficiency.
Carbon-concentrating Mechanism, Photorespiration, and Arginine Biosynthesis are Linked Under Photoautotrophic Growth Conditions-A common motif of regulation for growth under photoautotrophic condition is the induction of carbonconcentrating mechanism (CCM) related proteins. In C. reinhardtii and other aquatic unicellular organism, CCM is induced under low CO 2 availability to improve photosynthesis by increasing the concentration of intracellular inorganic carbon to elevate the CO 2 concentration at the site of ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco; for reviews see (106,107)). Rubisco catalyzes the first step of photosynthetic CO 2 fixation and beside the carboxylation of its substrate ribulose-1,5-bisphosphate (RuBP), it may also oxygenate it. Oxygenation produces the toxic component 2-phosphoglycolate (2-PG), which is recycled to 3-phosphoglycerate (3-PGA) during the process of photorespiration. Consequently an increase in the CO 2 to O 2 ratio will decrease photorespiration. Among the proteins that are found to be up-regulated under PA conditions compared with PH (communities VI and XI, Fig. 3) and found in a shared cluster are LCIC and LCIB. Studies of mRNA expression in response to low CO 2 revealed that transcripts of lcib and lcic are among the most induced mRNA molecules recorded (108,109). Genetic analyses indicated that LCIB is involved in the accumulation of inorganic carbon in the chloroplast (108) and that correspondingly to their co-expression, LCIB and LCIC form a 350 kDa hexameric complex in the chloroplast (110). Induction of CCM proteins is also strongly supported by the fact that CAH3 (for review see (106,107) is significantly increased under iron sufficient PA settings. In accordance, LCIB/LCIC and CAH3 were found in CCM clusters when transcriptomewide changes in relation to CO 2 and the CO 2 -CCM regulators CIA5/CCM1 were analyzed (111). Despite the fact that C. reinhardtii possesses a carbon-concentrating mechanism, photorespiratory products have been detected after shifting cells to low CO 2 by using gas chromatography-coupled time of flight mass spectrometry (112). Accordingly, transcript levels for the photorespiratory enzymes, including HRP1, were increased in low CO 2 grown cells (109,111). HPR1 is involved in photorespiration as explained above, and is found in a coregulated community with LCIB and LCIC (Fig. 3, community XI). Thus a close connection between photorespiration and the CCM can be interfered as already suggested (109,111). Coclustering of the glycine cleavage multienzyme complex, the serine glyoxylate aminotransferase and 10-formyltetrahydrofolate synthetase, as suggested above, indicates coordinated action in serine biosynthesis and supports that photorespiration is more active under PH low iron and PA/PH as compared with PH low iron. Other members of this assembly are a N-acetyl-gamma-glutamyl-phosphate reductase (ARGC1) and the large subunit of carbonyl phosphate synthase (Cre28.g776100.t1.1) (CPS), two enzymes that are involved in arginine biosynthesis. It is of note that these enzyme are all found in the chloroplast proteome of C. reinhardtii (47) and Arabidopsis thaliana (113,114). Knock-out of the small or large subunit of CPS resulted in a reticulate leaf phenotype that was correlated with a defect in mesophyll development (115), implicating that CPS function and arginine may play a role in these developmental processes. CPS is involved in the conversion of glutamine and bicarbonate into carbamoyl phosphate (CP) and glutamate (116). CP formation is the first committed step in the biosynthesis of arginine and in de novo biosynthesis of pyrimidines (117). Because CPS, some CCM proteins, HPR1, three subunits of the glycine cleavage multienzyme complex, the serine glyoxylate aminotransferase and 10-formyltetrahydrofolate synthetase are coregulated, the question arises about the connection between these biological processes and the arginine and pyrimidine biosynthesis. During photorespiration, the formation of glycine from glyoxylate via the glutamate:glyoxylate aminotransferase (GGT) requires glutamate as a cosubstrate (71). Thus the creation of CP as a committed step in the arginine biosynthetic pathway could participate in glutamate production as substrate for the photorespiratory process, linking thereby arginine biosynthesis and photorespiration. This goes hand in hand with the observed increase in 2-oxogluterate, which is produced via GGT action induced after shifting C. reinhardtii cells into low CO 2 conditions (112), whereas glutamate is decreased because it is steadily consumed. Notably arginine accumulates under stress in A. thaliana (118,119) and is considered an important intermediate for nitrogen storage (120). In line our data suggest that arginine biosynthesis, as well as photorespiration, is characteristic for cells actively performing photosynthesis.
Photoautotrophic Growth Regulation May Require Calcium and Redox Signaling-Beside proteins involved in CCM, photorespiration or arginine biosynthesis as presented in community XI, other proteins were particularly up-regulated under PA iron sufficient conditions (Fig. 3). These proteins are a stomatin-like protein (Cre33.g782550.t1.1), a protein that contains a thioredoxin domain and four putative EF-hand motifs (Cre03.g202950.t1.2) as well as a NADP malic enzyme (MME5). MME5 is a chloroplast localized NADP-malic enzyme (61) that either operates in decarboxylating malate to pyruvate or in catalyzing the carboxylation of pyruvate to malate thereby generating NADPH or oxidizing it to NADP ϩ . As the direction is dependent on the chloroplast redox poise, the enzyme could thereby be important for the fine-tuning of the chloroplast redox status.
In mammalians, stomatin and its paralogs are described, as homo-oligomeric, lipid raft-associated, integral membrane proteins that may interact with and modulate various ion channels and transporters (121). More particularly, stomatin was found to interact with a Ca 2ϩ ATPase (122). However, whether a link between C. reinhardtii stomatin and calcium exists remains to be shown.
Four EF hands domains were found in Cre03.g202950.t1.2. EF hands present a common motif for proteins involved in the binding of Ca 2ϩ . The EF is a helix-loop-helix motif, which mostly occurs in pairs (not identical; up to 6 pairs in one protein), allowing high affinity binding to Ca 2ϩ . The pairing is important for the Ca 2ϩ -binding mechanism and Ca 2ϩ -induced conformational change, translating the Ca 2ϩ concentration into an output response (123). EF-hand containing proteins operate in all aspects of cell function and more than 3000 EF-hand related entries can be found in the NCBI RefSeq database. Likewise thioredoxin domains are contained in many proteins andenzymaticallythioredoxinactsasaproteindisulphideoxidoreductase that catalyzes the reversible oxidation of two cysteine thiol groups to a disulfide, thereby transferring two electrons and two protons (124). Thus combining EF-hand motifs with a thioredoxin domain would link calcium and redox signaling. From network analyses it distinctly appears that thioredoxins are central to chloroplast redox signaling (125). Recent data revealed that chloroplast calcium signaling may also contribute to retrograde signaling between chloroplast and the nucleus. In C. reinhardtii, the expression of the nuclear encoded LHCSR3 is regulated via the thylakoid membrane protein CAS (Ca 2ϩ -sensing receptor) and Ca 2ϩ (126). In A. thaliana, CAS is involved in chloroplast mediated activation of plant immune signaling (127) and in stomatal regulation in response to elevations of external Ca 2ϩ through the modulation of cytoplasmic Ca 2ϩ dynamics (128,129), thereby regulating the availability of carbon dioxide and thus indirectly also the redox poise of the chloroplast stroma. Importantly, CAS has been also implicated in regulation of cyclic electron transfer (130). Here consistently a link between redox poise and calcium can be anticipated.
Notably the induction of a Peptidyl-prolyl cis-trans isomerase (CYN37) and a redox-related translation factor for expression of the chloroplast-encoded psbA gene (TBA1) (67), under PA/PH conditions, as already discussed for PH iron deficiency, support the conclusion that photoautotrophic growth regulation is related to the chloroplast redox poise. Authors declare no competing commercial interests.