Proteomic Analysis of a Fraction with Intact Eyespots of Chlamydomonas reinhardtii and Assignment of Protein Methylation.

Flagellate green algae possess a visual system, the eyespot. In Chlamydomonas reinhardtii it is situated at the edge of the chloroplast and consists of two carotenoid rich lipid globule layers subtended by thylakoid membranes (TM) that are attached to both chloroplast envelope membranes and a specialized area of the plasma membrane (PM). A former analysis of an eyespot fraction identified 203 proteins. To increase the understanding of eyespot related processes, knowledge of the protein composition of the membranes in its close vicinity is desirable. Here, we present a purification procedure that allows isolation of intact eyespots. This gain in intactness goes, however, hand in hand with an increase of contaminants from other organelles. Proteomic analysis identified 742 proteins. Novel candidates include proteins for eyespot development, retina-related proteins, ion pumps, and membrane-associated proteins, calcium sensing proteins as well as kinases, phosphatases and 14-3-3 proteins. Methylation of proteins at Arg or Lys is known as an important posttranslational modification involved in, e.g., signal transduction. Here, we identify several proteins from eyespot fractions that are methylated at Arg and/or Lys. Among them is the eyespot specific SOUL3 protein that influences the size and position of the eyespot and EYE2, a protein important for its development.

Introduction

Many motile algae exhibit a peculiar photo-behavior: movement toward or away from the light source depending on the light intensity and quality. This behavior is known as positive or negative phototaxis. To allow such precise movement responses, many flagellate algae of all major phylogenetic lineages have developed specialized optical devices (eyespots), which are antennae that determine the direction of incident light (Foster and Smyth, 1980; Kreimer, 1994, 2009). In some warnowiid dinoflagellates, this structure is extremely complex and called ocelloid. Its ultrastructure bears apparent resemblance to camera-type eyes of some animals and has recently been shown to be composed of different specialized cell organelles such as plastids, mitochondria, and vesicles. The ocelloid is probably homologous to simpler dinoflagellate eyespot types, most of which involve parts of the chloroplasts (Dodge, 1984; Kawai and Kreimer, 2000; Gavelis et al., 2015; Hayakawa et al., 2015). In green algae, the functional eyespot is also a composed “organelle”, involving local specializations from different subcellular compartments and highly ordered carotenoid-rich globules inside the chloroplast. It is peripherally located in the cell and readily observable by bright-field microscopy as an orange- to red-colored spot (Melkonian and Robenek, 1984; Kreimer, 2009). In Chlamydomonas reinhardtii, the eyespot typically consists of two highly ordered layers of such globules, each subtended by a thylakoid (see Figure for a schematic drawing). The outermost globule layer is attached to specialized areas of both the chloroplast envelope and the adjacent PM. The globule layers modulate the light intensity reaching the photoreceptors in the PM patch as the cell rotates around its longitudinal axis during forward swimming. They function as a combined absorption screen and quarter-wave interference reflector. Thereby, the contrast at the photoreceptors is increased up to eightfold, making the whole optical system highly directional (Foster and Smyth, 1980; Kreimer and Melkonian, 1990; Harz et al., 1992; Kreimer et al., 1992). Beside its function as a sensor for light direction and quality, the eyespot might have potential additional roles mainly for chloroplast function. There is, e.g., increasing evidence from both, ultrastructural and proteomic data, for a link between eyespot globules (EG) and plastoglobules (PG) and thereby for a more general role in chloroplast metabolism such as the biosynthesis of prenylquinones and carotenoids (Kreimer, 2009; Nacir and Bréhélin, 2013; Davidi et al., 2015). Whereas PG are directly connected with the thylakoids via the outer lipid leaflet of the TM (Austin et al., 2006), it is currently not known whether this is also true for the globules of the eyespot. Methods for the isolation of PG from green algae as well as EG and fragments are established (Kreimer et al., 1991; Renninger et al., 2001; Schmidt et al., 2006; Davidi et al., 2015).

Schematic drawing of the eyespot in Schematic longitudinal section through the eyespot in the region of the four stranded microtubular root (D4), which is important for eyespot positioning. The functional eyespot consists of local specializations of different compartments. The part of the plasma membrane (PM) overlying the eyespot globules (EG) is thickened and in close association with the inner and outer chloroplast envelope (IM/OM). The EG are in contact with the thylakoid membrane (TM) and arranged in layers, the outermost being also in contact with the chloroplast envelope. The possible link between plastoglobules (PG) and EG is indicated. (B) Separation of the cell homogenate in the first gradient using the established eyespot isolation method by Schmidt et al. (2006). (C) Separation of the cell homogenate in the first gradient using the here described (see Materials and Methods) modified cell rupture method for an extended (“e”) eyespot fraction; separation of fraction 2 (F2e) from these gradients in fraction 2A (F2Ae) and 2B (F2Be) in the final gradient. Note the increased amount of F2e and the decreased amount of F1e and thylakoid debris in comparison to the original cell rupture method.

Until 2005, only six proteins relevant to the structure and function of the eyespot of C. reinhardtii were identified, mainly based on genetic approaches. These included EYE2 and MIN1, two proteins important for eyespot assembly (Roberts et al., 2001; Dieckmann, 2003), two splicing variants of the abundant retinal binding protein COP, and the two unique seven-transmembrane domain photoreceptors, Channelrhodopsin1 (ChR1), and ChR2 serving as basis for the development of optogenetics (Deininger et al., 1995; Fuhrmann et al., 2001; Nagel et al., 2002; Sineshchekov et al., 2002; Suzuki et al., 2003; Hegemann and Nagel, 2013). As an in depth knowledge of the protein composition of this complex organelle is one prerequisite to understand its function at a molecular level, also proteomic approaches to a fraction strongly enriched in eyespot fragments were applied. Thereby, 203 different proteins covered with at least two different peptides were identified and some of them were shown to be phosphorylated (Schmidt et al., 2006; Wagner et al., 2008). Based on these proteomic approaches, recent functions in eyespot development have been demonstrated for SOUL3, one of the five SOUL heme-binding proteins in C. reinhardtii (Merchant et al., 2007; Schulze et al., 2013), Casein kinase 1 and the blue-light receptor Phototropin (Schmidt et al., 2006; Trippens et al., 2012). Additionally, evidence for a specialized localization of the alpha- and beta-subunit of the chloroplast ATP synthase and thereby probably function within the eyespot was found (Schmidt et al., 2007). The eyespots used in the proteomic approach by Schmidt et al. (2006) retained, however, only small parts of the eyespot membranes. Thus, it is well feasible that core membrane-associated eyespot proteins were not identified in this study. To further increase the understanding of eyespot related signaling processes in depth, knowledge of the protein composition of the membranes in its close vicinity is desirable. We thus developed a method that allows isolation of intact eyespots containing the two layers of lipid rich globules still associated with large parts of the eyespot membranes and the adjacent PM as well as chloroplast envelope areas and analyzed this fraction by a proteomic approach.

Beside phosphorylation, protein methylation is an important posttranslational modification involved in the regulation of protein stability, localization, activity, and protein–protein interactions (Reinders and Sickmann, 2007; Biggar and Li, 2015). The importance of protein methylation at Arg and Lys residues in regulating protein activity is also becoming apparent in vascular plants and algae (e.g., Deng et al., 2010; Blifernez et al., 2011; Werner-Peterson and Sloboda, 2013). Recently, Arg and Lys methylation was reported for several chloroplast proteins (Alban et al., 2014). We therefore assigned protein methylation sites in the previous and currently analyzed eyespot fractions. In this study, we identify 742 proteins of the extended eyespot fraction, including numerous membrane-associated proteins as intended beside other candidates for eyespot development and signaling. We also assign protein methylation sites on two proteins with already demonstrated important functions for the eyespot (SOUL3 and EYE2; Roberts et al., 2001; Boyd et al., 2011; Schulze et al., 2013) and on 23 other proteins from the eyespot fractions, among them the alpha-, beta- and I-subunits of the chloroplast ATP synthase and a 14-3-3 protein.

Materials and Methods

Isolation of the Different Eyespot Fractions

Growth of 20 L C. reinhardtii strain cw15 to late log-phase and isolation of fraction 1 (Figure ) enriched in eyespot fragments was done as previously described (Schmidt et al., 2006). The isolation of a fraction enriched in largely intact eyespots (F2Ae in Figure ) was achieved by reducing the power and duration of the ultrasonic treatment (Bandelin Sonopuls HD2070, Microtip HS 73) to 16% output and seven cycles (15 s each interrupted by 15 s of cooling) during cell rupture. In addition, the submersion depth of the sonotrode tip was reduced to 1.5 cm. Eyespot fragments were separated by discontinuous sucrose gradients buffered with gradient stock solution (GSS) consisting of 12.5 mL sample brought to 42% (w/v) sucrose, 10 mL 31.8% (w/w), and 10 mL 20.5% (w/w) sucrose, overlaid with 5.5 mL GSS. After centrifugation (100,000 g, 105 min, 4°C), the orange-green bands (F2e in Figure ) at the interface of 20.5 and 31.8% sucrose were collected and brought to 35% (w/v) sucrose. Fraction F2e was further purified by flotation centrifugation (100,000 g, 60 min, 4°C) on discontinuous sucrose gradients (15 mL sample, 3 mL 31.8% [w/w] sucrose, 18 mL 25% [w/w] sucrose, and 2 mL 20.5% [w/w] sucrose). The yellow–orange 20.5% sucrose fractions (F2Ae) were collected and brought to 32% (w/v) sucrose. For concentration, 10 ml of this fraction were layered on 26 mL 42% (w/w) sucrose overlaid by 2 mL GSS and centrifuged again (100,000 g, 45 min, 4°C). The concentrated fraction F2Ae was collected from the top of the gradient, directly extracted with chloroform:methanol:water (4:8:3) and the precipitated proteins were washed several times with methanol:chloroform (2:1, v/v) to remove lipids and pigments before they were dissolved in 2x SDS sample buffer. Cell harvesting, buffers and all other treatments were otherwise done as described by Schmidt et al. (2006).

Crude Extract Preparation and Electrophoretic Methods

Log phase cultures were supplemented directly prior to cell harvesting with 1 mM phenylmethylsulfonyl fluoride (PMSF). Cells were harvested by centrifugation (2000 g, 10 min, 4°C), and pellets were suspended in TNED buffer (20 mM Tris, 80 mM NaCl, 1 mM EDTA, and 1 mM DTT, pH 7.5) supplemented with 1 mM PMSF and complete protease inhibitor cocktail (Roche), according to the instructions of the supplier. In the case of the cell-wall less strain cw15, aliquots (300 μl) of the suspended cells were directly mixed with methanol:chloroform (2:1, v:v; 1200 μL) for parallel protein precipitation and lipid/pigment removal. Cell wall-possessing strains were homogenized by sonication. Protein solubilization, SDS-PAGE and protein staining with silver or colloidal Coomassie were conducted as described by Schmidt et al. (2006). Gel loading based on equal protein content and immunoblot analyses followed standard techniques. Incubations with the monoclonal anti-modified (monomethyl, asymmetric dimethyl, and symmetric dimethyl) methyl-arginine antibody 7E6 (anti-Rm; company Covalab SAS, France, 1:2,500) were done overnight at 4°C.

Electron Microscopy

For fixation, aliquots of F2Ae were diluted 1:1 with ice-cold 40% (w/v) BSA in GSS and mixed with cold glutaraldehyde (final concentration 3.1 or 6.3 %). The BSA/F2Ae gel pieces were sliced in small pieces and fixed overnight in 25 mM Hepes/NaOH (pH 7.8) and 3.5% glutaraldehyde at 4°C. Samples were washed three times for 10 min with 25 mM Hepes/NaOH (pH 7.8) prior to a 2 h incubation at 4°C with osmium tetroxide (1% in 25 mM Hepes/NaOH, pH 7.8), followed by four washing steps with Millipore water. Dehydration, embedding in Epon and staining of ultrathin sections with uranyl acetate/lead citrate was done as described (Reynolds, 1963; Renninger et al., 2001). EM negatives were scanned and processed with Photoshop (Adope Systems).

Mass Spectrometry (MS) Analysis

In-Gel Digestion and Nano HPLC Electrospray Ionization Tandem MS (LC-ESI-MS/MS)

The gel was dissected into 50 bands and each gel slice was subjected to a washing procedure followed by an in gel digestion with trypsin as previously described (Schmidt et al., 2006). The pellets containing the tryptic peptides were resuspended in 5 μL 5% (v/v) acetonitrile/0.1% (v/v) formic acid and subjected to nano LC-ESI-MS using an UltiMateTM 3000 nano HPLC (Dionex Corporation) with a flow rate of 300 nL/min coupled on-line with a linear ion trap electrospray-ionization (ESI) mass spectrometer (FinnigianTM LTQTM, Thermo Electron Corp.) as described previously (Schmidt et al., 2006). The mass spectrometer was cycling between one full MS and MS/MS scans of the four most abundant ions. After each cycle, these peptide masses were excluded from analysis for 3 min.

Data Analysis

Data analysis was done using the Proteome Discoverer software (Version 1.4) from Thermo Electron Corp. including the SEQUEST algorithm (Link et al., 1999). Searches were done for tryptic peptides allowing two missed cleavages. The software parameters were set to detect lysine modifications (+14.016 for monomethyl, +28.031 for dimethyl, and +42.047 for trimethyl) and arginine modifications (+14.016 for monomethyl and +28.031 for dimethyl) with a maximum of four modifications of one type per peptide. Further, detection of variable methionine oxidation (+15.995) was enabled. Peptide mass tolerance was set to 1.5 Da in MS mode. In MS2 mode, fragment ion tolerance was set up to 1 Da. The parameters for all database searches were set to achieve a false discovery rate (FDR) of not more than 1% for each individual analysis (Veith et al., 2009). All spectra used for the assignment of protein methylation were in addition manually validated and checked for the presence of the b- and/or y-type ions representing the methylation sites. In case that a methylation site is present on the first amino acid at the N-terminus of the peptide and therefore not present in the b- and y-type ions, these methylation sites were considered only when all other possible methylation sites in the peptide were validated as methylated or non-methylated, respectively. In case of trimethylated Lys residues, spectra were additionally screened for the presence of neutral loss events of trimethylamine (59 Da), which is a signature for trimethylation (Zhang et al., 2004; Erce et al., 2012). Data were searched against the C. reinhardtii database (Vs. 5.3.1), hosted by Phytozome (Vs. 9.1)[1], the mitochondrial database available from NCBI (NC001638, [gi:11467088]) and the chloroplast database[2]. Data from all runs were combined and further evaluated using an in house developed program. The peptide sequences of the gene models were compared to the NCBI protein database using BLAST (Altschul et al., 1997). For positive identification of both, protein and functional domain prediction, an internal cut-off E-value of 1e-05 was used. Transmembrane domain information was based on predictions by the programs TMHMM (Krogh et al., 2001), TMpred (Hofmann and Stoffel, 1993), and TopPred (von Heinjne, 1992). The GRAVY index was determined with ProtParam (Gasteiger et al., 2005). The mass spectrometry (MS) proteomics data have been deposited to the ProteomeXchange Consortium[3] via the PRIDE partner repository (Vizcaino et al., 2013) under the project description “Extended eyespot proteome and methylated proteins in the eyespot of C. reinhardtii strain cw15 (ID: PXD003254)”.

Miscellaneous

Protein content was determined according to Neuhoff et al. (1979) with BSA as standard. Chlorophyll and carotenoids were determined as described by Lichtenthaler (1987). Spectra were recorded in 90% methanol with an Ultrospec 2100 pro (GE Healtcare Life Science). Images of gradients, gels and blots were captured with a Coolpix 990 (Nikon) and processed with Photoshop (Adobe Systems).

Results

Isolation and Characterization of a Fraction Enriched in Intact Eyespots

Based on a previously established method for the isolation of eyespot fragments from C. reinhardtii (Schmidt et al., 2006), we developed a procedure that allows the isolation of largely intact eyespot apparatuses. They contain the two layers of lipid rich globules with large parts of the eyespot associated membranes and adjacent PM areas. As a visual marker for enrichment of eyespots within the sucrose gradients, we first used the conspicuous yellow/orange to red color of the EG (Figures ). We noticed a deep orange-red fraction enriched in highly purified eyespot fragments (named F1; Schmidt et al., 2006), as well as an additional weak yellow–orange fraction (F2) apart from the bulk of the cell and chloroplast debris. F2 accumulates at the interface between 20.5 and 31.8% sucrose. At these sucrose concentrations, enrichment of largely intact eyespots with a high proportion of eyespot-associated membranes was described for the green alga Spermatozopsis similis (Kreimer et al., 1991). Thus, it seemed possible that the Chlamydomonas F2 fraction that was always present only in minute amounts may also bear intact eyespots. We therefore tried to optimize the yield of fraction F2 by variation of the ultrasonic treatment used for cell rupture. Best results were obtained upon a 50% reduction of the output power and cycle number (32 to 16% and 14 to 7 cycles) in combination with a reduction of the submersion depth of the sonotrode tip to 1.5 cm. For a clear differentiation between the different eyespot preparations, all fractions obtained by this modified cell rupture/purification strategy are characterized by the subscripted “e”, which stands for extended eyespot. Although these reductions result in a clear loss in general cell rupture and thereby total yield of F1e, the relative amount of F2e, as judged by its color, was increased compared to the standard procedure (Figures ). F2e was further purified by an additional gradient centrifugation step to reduce contamination by thylakoids and other cell organelles. In this subsequent gradient, F2e splits up into two fractions: an orange–yellow (F2Ae) and a greenish-yellow (F2Be) fraction (Figure , final gradient). F2Ae was finally concentrated by a floating centrifugation step. The total protein yields of F2Ae varied between 28 and 46 μg when starting with 20 L of culture.

To verify enrichment of eyespots and to judge their intactness and purity, F2Ae was analyzed by transmission electron microscopy (Figure ). Indeed, a significant number of well-defined eyespots with intact globule layers are enriched in this fraction (Figures ). The close packing of the globules is largely preserved and the layers are still subtended by a thylakoid. The average size of the globules in these layers is well in the in situ range of 80–130 nm (Melkonian and Robenek, 1984). Analyses of 100 globules from different eyespots and sections yielded values of 114 ± 28 nm for the width and 121 ± 30 nm for the height, respectively. Only a few EG were enlarged, probably by fusion. In addition, large parts of the chloroplast envelope, the PM and a fuzzy appearing material observed between the chloroplast envelope and the PM in the eyespot region in situ were frequently present (Figures ; Melkonian and Robenek, 1984; Kreimer, 2009). Additionally, less well-preserved eyespots predominantly with enlarged, possibly fused, globules were also observed (Figure ). Along with the intact eyespots, still possessing all eyespot membranes, varying amounts of the adjacent PM and chloroplast envelope were co-isolated (Figures ). Depending on their length, these membranes partially enwrapped the isolated eyespots as well as cell components not belonging to the functional eyespot, but located in situ in its close proximity like, e.g., grana thylakoids and ribosomes (Figure ). Major free contaminants were individualized fused globules associated with membranes, membrane vesicles of unknown origin as well as stroma and grana thylakoids (Figures ). Occasionally, also mitochondrial fragments were detected (Figure ). In summary, transmission electron microscopy confirmed enrichment of intact eyespots associated with all eyespot membranes. The gain in intactness of eyespots seems, however, accompanied by an inevitable increase of non-eyespot contaminants.

Characterization of the extended eyespot fraction 2A (F2A Overview (A) and details (B–F). (A) Black arrows indicate well-preserved eyespots, whereas the white arrow indicates structurally less well-preserved eyespots with fused EG; scale bar: 600 nm. (B,C) Cross-sections through isolated, double-layered eyespots associated with eyespot and other membranes. The globule layers are subtended by a thylakoid (black asterisk). Small arrowheads: plasma membrane area overlying the EG. Note the close association with the chloroplast envelope membranes and the thickened, electron-dense appearance of the plasma membrane in this region. In addition, the fuzzy fibrillar material typically observed in situ between the plasma membrane and chloroplast envelope in the region of the eyespot apparatus (Kreimer, 2009) is preserved. Small arrows point to normal plasma membrane regions tending to enwrap the globule layers and not eyespot-associated structures like a thylakoid stack (white asterisk) or ribosomes (large arrowheads); scale bars: 150 nm. (D) Tangential section through a globule layer. Arrowheads indicate membranes enclosing the eyespot; scale bar: 250 nm. (E,F) Co-isolated free mitochondrial fragment (E) and free thylakoid stack (F); scale bars: 400 nm.

Spectral analysis of the pigments in fraction F2Ae was performed for further characterization. The visual appearance of F2Ae already indicates that carotenoids are strongly enriched. In vegetative cells of C. reinhardtii, carotenoids are found in the thylakoids, the chloroplast envelope and EG, whereas chlorophyll is solely present within thylakoids. A comparison of the amount of chlorophyll present in the crude extract (CE) with the amount of chlorophyll in fraction F2Ae revealed that 0.015% of the total chlorophyll applied to the gradients remained there. The average carotenoid:chlorophyll ratio measured as the absorbance ratio 478 to 680 nm was 40 ± 12.6. These values indicate effective removal of the majority of free thylakoids not associated with the EG layers. Compared to fraction 1, however, these values are significantly higher (<0.0005% and 60–70; Schmidt et al., 2006). Spectra of methanol extracts of fraction F2Ae were identical to those of fraction F1 in the range between 350 and 525 nm and revealed a typical carotenoid spectrum (Figure ). Major absorption peaks where at 446 and 470 nm and a shoulder at 423 nm.

Spectral analysis and SDS-PAGE of the different eyespot fractions. (A) Normalized absorption spectra of fractions F1 and F2Ae in 90% methanol in comparison to that of a crude cell extract (CE). (B) Comparison of the protein patterns of fractions F1, F2Ae, F2Be, and F2e. Total proteins (10 μg) of the different gradient fractions were separated by SDS-PAGE (11%) and stained with silver. The positions of molecular mass markers are indicated on the left in kilodaltons.

The above analyses clearly shows that eyespot purity of the extended eyespot fraction F2Ae was reduced compared to that of fraction F1 (Schmidt et al., 2006). Nonetheless, the significantly increased amount of eyespot associated membranes as well as of membrane areas adjacent to them in fraction F2Ae offers the possibility to identify novel putative eyespot proteins and those present in close vicinity to this visual system via a proteomic approach. To get a first impression of possible differences in the protein patterns of the eyespot fractions, a comparative SDS-PAGE analysis was conducted (Figure ). As expected, the general protein pattern of the fractions enriched in eyespot fragments (F1) and intact eyespots (F2Ae) were similar to some extent indicating an overlap of a majority of proteins. In the complex protein patterns, however, also subtle but clear differences were resolved (see, e.g., the molecular mass ranges between 35 and 45 kDa and below 25 kDa). Their protein patterns clearly differed also from those fractions with higher chlorophyll content (F2e and F2Be). The protein pattern of fraction F2Ae was highly reproducible in independent purifications [Supplemental Figure S1A (Data Sheet 1)].

The Proteome of the Extended Eyespot

To identify proteins of the enriched F2Ae eyespot fraction by MS/MS, proteins from two independent eyespot isolations were combined, separated by SDS-PAGE and the gel was sliced into 50 pieces [Supplemental Figure S1B (Data Sheet 1)]. Following in-gel digestion with trypsin, the generated peptide fragments were subjected to LC-ESI-MS/MS analyses using a linear ion trap mass spectrometer (see Materials and Methods). In total, 742 proteins were identified with at least two different peptides and a FDR of ≤1. The majority of the 203 proteins identified in the former fraction F1 with at least two different peptides (Table 1 in Schmidt et al., 2006) were positively confirmed to be present in the extended eyespot of F2Ae [Supplemental Table S1 (Data Sheet 3) along with the currently identified peptides and values]. Supplemental Table S2 (Data Sheet 4) lists the 27 missing candidates, 10 of them (category ribosomes, translation, and DNA related) being obvious contaminants. Reasons for the differences are discussed later.

A selection of promising novel candidates with regard to eyespot development, signal transduction as well as membrane association and transporters is presented in Table . This selection is based on current knowledge about the eyespot apparatus, including its different complex associated subcellular structures as well as signaling pathways important for photoorientation, its development including mutant analysis as well as information on PG. Detailed information for all identified peptides of these selected candidates is given in Supplemental Table S3 (Data Sheet 5). We think that these categories enclose the primarily most interesting candidates to be studied functionally in the future. Supplemental Table S4 (Data Sheet 6) lists all further identified novel proteins along with their peptides from F2Ae including also the putative contaminants. Among the novel candidates (Table ), a protein involved in eyespot assembly, the Ser/Thr kinase EYE3 (Boyd et al., 2011), a retinal pigment epithelial membrane receptor as well as a protein with homologies to the eye pigment precursor transporter protein family (Vopalensky et al., 2012) were found. Moreover, two SOUL heme-binding proteins (SOUL2 and SOUL5) were identified. The so far known SOUL3 eyespot protein influences size and position of the eyespot (Schulze et al., 2013). Excitation of the ChR photoreceptors initiates a Ca2+-based signaling cascade toward the flagella in the eyespot (reviewed by Hegemann and Berthold, 2009). Further, rapid reversible protein phosphorylation in isolated eyespots is strongly affected by free Ca2+ concentrations between 10-8 and 10-7 M (Linden and Kreimer, 1995) and both ChRs as well as SOUL3 are targets of kinases (Wagner et al., 2008). Therefore, the identification of six novel calcium-sensing and binding proteins as well as 11 different kinases and four phosphatases could be of high interest for understanding the signaling cascade network in the eyespot. Notably, a Ca2+-dependent protein kinase 1 was identified with 11 different peptides. Moreover, two 14-3-3 proteins that interact with phosphorylated proteins and are key regulators in many vital cellular processes including signal transduction (Denison et al., 2011) were identified with eight and ten different peptides, respectively. Of note is also the fact that some of the kinases belong to the ABC1 kinase family such as EYE3 (Table ). Homologs of two members of this family identified in the proteome of eyespot fraction F1 are present in PG of Arabidopsis thaliana and the beta-carotene PG of Dunaliella bardawil (Ytterberg et al., 2006; Lundquist et al., 2012; Davidi et al., 2015). Two of the newly identified ABC1 kinases with an AarF domain from the F2Ae fraction [Cre09.g407800.t1.3 (AKC1), Cre13.g570350.t1.3 (AKC4)] are also present in Arabidopsis PG (Lundquist et al., 2012). Additionally, AKC1 has a homolog in the plastoglobule proteome of D. bardawil (Davidi et al., 2015).

Among the novel candidates in the group of membrane associated/structural proteins are three PAP-fibrillin domain-containing proteins. Eight proteins of this group important for stabilization of PG and EG (Renninger et al., 2001; Ytterberg et al., 2006) were already identified in our previous analysis (Schmidt et al., 2006). In addition, Vesicle inducing protein in plastids (VIPP1) that is involved in the biogenesis of thylakoids (Rütgers and Schroda, 2013) has been also found with 11 different peptides (Table ). VIPP1 was detected recently in an algal plastoglobule proteome and it was postulated that PG evolved from eyespot lipid droplets (Davidi et al., 2015). In addition, two fasciclin-like proteins with three and four FAS1 domains, respectively, might be of interest. FAS1 domains are present in many secreted, membrane-anchored proteins. Five transporters were found as well, underlining the enrichment of eyespot associated membranes (Table ). In this group, identification of three P-type PM ATPases and a PM type Calmodulin-binding Ca2+-transporting ATPase is of special interest. These types of proteins may be important for the resting membrane potential at the PM and thereby eventually also affect the excitation of the cell through the ChRs. As the currents of these photoreceptors are carried by Ca2+, its extrusion in the eyespot region following ChR excitation will be important for signaling and adaptation. Due to the extreme close vicinity of the PM and the chloroplast envelope in the eyespot, also the chloroplast K+/H+-efflux antiporter 2 (KEA2) might affect the resting membrane potential in the eyespot region. KEA2, as one of the PM-type H+-ATPases, was found with 11 different peptides. These different types of transporters thus might probably have important indirect functions in the context of eyespot related signaling.

Assignment of Protein Methylation

In recent years, increasing evidence points to the importance of diverse posttranslational protein modifications beside reversible phosphorylation in the regulation of protein stability, localization, activity and protein–protein interactions in cell organelles (e.g., Lehtimäki et al., 2015). One of them is methylation. In a recent study, 23 chloroplast proteins have been shown to be methylated at Lys and/or Arg residues (Alban et al., 2014). Methylation of these amino acids increases their basicity and hydrophobicity without altering their charge (Rice and Allis, 2001). As a part of the functional eyespot involves local specializations of the chloroplast (Figure ), it was of interest to find out whether methylated proteins are present in the eyespot fractions. Prior to an analysis via MS, we checked with a commercially available anti-methyl-Arg specific antibody recognizing monomethyl, asymmetric dimethyl, and symmetric dimethyl-Arg whether methylated proteins are still detectable after the isolation of eyespots (Figure ). Immunoblot analysis of the proteins from fraction F1 with this antibody revealed four clearly and several weakly labeled protein bands. This result demonstrates that methylated proteins are still present in the eyespot fraction after the isolation procedure. Position, number, and intensity of labeled bands in a CE clearly differed, indicating that a specific subset of methylated proteins may occur in the eyespot fractions.

Western blot analysis of proteins from the eyespot fraction F1 and a CE with the anti-methyl-arginine specific antibody 7E6. Proteins (8 μg) were separated by 11% SDS-PAGE, transferred to a PVDF membrane, and either analyzed with a monoclonal anti-methyl-arginine antibody (Anti-Rm; 1:2,500) or stained with Coomassie Brilliant Blue R250 (Coomassie).

Mass spectrometry analysis of the peptides from both eyespot fractions resulted in the identification of 25 methylated proteins (Tables and ) where the methylation sites could be specified after additional manual evaluation of the spectra [Supplemental Figure S2 (Data Sheet 2), Supplemental Table S5 (Data Sheet 7)]. These proteins were characterized by 36 different methylated peptides and nine additional overlapping peptides having variations in the methylation status and/or sites and/or in oxidized versus non-oxidized Met. Six proteins are methylated at Lys and Arg residues, 17 proteins only at Lys residues, and two only at Arg residues. In total, we detected 10 Arg sites (nine mono and one dimethyl) and 42 Lys methylation sites (23 mono-, 14 di-, and 8 trimethylations; two of the residues were found mono- as well as dimethylated and one further residue di- as well as trimethylated; Tables and ). It should be mentioned that the identification of trimethylated Lys sides by LC-ESI-MS/MS is hampered by the fact that modification by acetylation is very close in mass (42.04695 vs. 42.01056 Da; Alban et al., 2014). Identification of neutral loss of trimethylamine (59 Da) being specific for trimethylation allows to discriminate (Zhang et al., 2004; Erce et al., 2012; Alban et al., 2014). All seven peptides (Tables and ) with potential trimethylated sites were therefore screened for neutral loss events by analyzing the spectra manually. In all corresponding spectra, b- and y-type ions of neutral loss fragments were found [Supplemental Figure S2 (Data Sheet 2)], indicating that there is trimethylation and not acetylation. Moreover, 14 methylated proteins were identified that are nucleic acid-related or considered as obvious contaminants. They are listed in Supplemental Table S6 (Data Sheet 8).

Interestingly, two proteins connected to eyespot development and positioning, EYE2 and SOUL3, as well as one of the 14-3-3 proteins were detected among the methylated proteins. In the 14-3-3 protein, a methylated Arg and a dimethylated Lys residue were detected. As in the chloroplast of Arabidopsis (Alban et al., 2014), the Chlamydomonas chloroplast ATP synthase subunits alpha and beta both contain methylation sites. In addition, subunit I of the CF0 membrane part of ATP synthase bears a trimethylation site. Also, other TM associated components bear methylation sites such as photosystem I and II components or the light harvesting complex protein Lhca7. N-terminal methylation of core light-harvesting complex proteins was shown before in purple photosynthetic bacteria (Wang et al., 2002). Furthermore, the photosystem I assembly protein Ycf4 was methylated at one Lys residue. In addition, methylated peptides of Ycf10, a transporter of the inner chloroplast envelope possibly involved in H+ extrusion into the cytosol in algae (Spalding, 2008), and the Ferredoxin-NADP reductase were found. This enzyme of C. reinhardtii was already characterized earlier as a methylated protein (Decottignies et al., 1995). It possesses three methylated Lys residues (Lys109, 115, and 161; as the first 26 amino acids were missing in the study from 1995, position of the Lys residues were changed according to the actual sequence of the protein). These positions are remarkably similar to those we obtained (Lys113, Lys162) and coincide in one case (Lys109). Additionally, several proteins of yet unknown functions bear methylation sites (Table ). One protein, Cre01.g000900.t1.2 with two DUF1350 domains, was methylated only at Arg residues.

Discussion

The functional eyespot apparatus involves parts of numerous subcellular compartments including different membranes, proteins for its development and positioning, as well as proteins known from PG that seem to be involved in its structural organization and preservation. In addition, proteins involved in signaling, adaptational responses, and biochemical pathways, like, e.g., retinal biosynthesis, are part of this light sensitive organelle. Due to its high ultrastructural complexity and its diverse associated processes, selection criteria for proteins associated specifically with the eyespot cannot be derived on the basis of a simple routine workflow procedure. Moreover, conserved targeting sequences for the eyespot are not evident so far and some proteins in the eyespot also occur in additional compartments such as Casein kinase 1 or Phototropin (Schmidt et al., 2006; Trippens et al., 2012). Comparable green algal eyespot proteomes (beside the one in Schmidt et al., 2006) are still missing, but algal lipid droplets and higher plant plastoglobule proteomes are known and can be used for comparison (e.g., Lundquist et al., 2012; Davidi et al., 2015). Moreover, mutant screens for the eyespot are available (e.g., Boyd et al., 2011). All these literature based information has been used for the selection of highly favorable candidates for the extended eyespot proteome listed in Table . Thereby, we cannot exclude that we have missed certain components or that some false positives were included. For example, light signaling includes reversible phosphorylation, but it is hard to judge if a predicted kinase or phosphatase (without any further information) may be due to a contamination. Thus, additional functional studies as they were done for Casein kinase 1 (Schmidt et al., 2006) and Phototropin (Trippens et al., 2012) that is also present in the cytosol and flagella of C. reinhardtii, or the eyespot-specific SOUL3 (Schulze et al., 2013) and the newly identified EYE3 (Table ; Boyd et al., 2011) are crucial. The present listed candidates in Table provide, however, a good basis for starting such functional tests.

Eyespot components involved in signal transduction are mainly in connection with the light signaling pathway. In the phototactic response of C. reinhardtii extracellular Ca2+ and Ca2+ fluxes are intricately involved. Both ChRs are directly light gated ion channels, which conduct Ca2+ under physiological conditions. Their excitation initiates fast inward-directed complex currents in the eyespot region, which finally produce a Ca2+ dependent alteration of the flagella beating and thereby leads to the steering of the cell toward or away from the light source (for reviews, see Witman, 1993; Hegemann and Berthold, 2009; Kreimer, 2009). Both ChRs as well as SOUL3, a protein important for the size and position of the eyespot, are targets of kinases (Wagner et al., 2008). Furthermore, in isolated green algal eyespots an increase in the free Ca2+ concentration from 10-8 to 10-7M is known to affect protein phosphorylation and two protein phosphatases with PP2Cc domains are among the dominant proteins in the eyespot proteome (Linden and Kreimer, 1995; Schmidt et al., 2006). Thereby, reversible phosphorylation as well as Ca2+ signaling in the eyespot region are primary events after excitation of the ChRs. In this context, the increase in kinases and phosphatases as well as in calcium sensing components in the extended eyespot proteome is of high interest for future studies. Especially the identification of the Calcium dependent protein kinase 1 (CDPK 1) is noteworthy, as no CDPKs have been detected in the core eyespot proteome (Schmidt et al., 2006). The Chlamydomonas genome encodes 14 CDPKs; three of them (1, 3, and 11) are found in the flagellar proteome (Pazour et al., 2005; Liang and Pan, 2013). Although CDPK1 has important functions in the flagella, it is also strongly present in the cell body (Liang et al., 2014; Motiwalla et al., 2014). A similar situation is found for Casein kinase 1. This kinase is specifically enriched in both, flagella and the eyespot, and has been shown to be involved in circadian phototaxis as well as in different flagella functions (Schmidt et al., 2006; Wirschell et al., 2011; Boesger et al., 2012). Similarly, the identification of additional members of the ABC kinase family in the extended eyespot proteome could be of interest for eyespot function and development. One of them, EYE3, is already known to be important for eyespot development. Mutants in EYE3 lack a visible eyespot (Lamb et al., 1999; Boyd et al., 2011). Four members of this group found in the eyespot proteomes (Cre13.g581850.t1.2; Cre03.g158500.t1.1; Cre09.g407800.t1.3; Cre13.g570350.t1.3) have homologs in PG of A. thaliana. Two of them, which were newly identified in the extended eyespot proteome, are known as AKC1 and AKC4. In higher plants, six of in total eight members of the ABC1-like kinases found in chloroplasts are associated with PG. They play a role in regulation of phylloquinone biosynthesis, redox recycling under high light and probably have other not yet known regulatory functions (Ytterberg et al., 2006; Lundquist et al., 2012; Martinis et al., 2013; Spicher and Kessler, 2015). With two 14-3-3 proteins, which interact with phosphoproteins, members of another protein group often involved in diverse signaling networks including those related to light perception were identified. For example, in higher plants 14-3-3 proteins interact with the blue light receptor Phototropin 1 (Sullivan et al., 2009). In this context, it is interesting that Phototropin in C. reinhardtii has recently been shown to affect eyespot size, the content of ChR1 as well as the sign of phototaxis (Trippens et al., 2012).

The procedure for purifying intact eyespots resulted also in an increase in transporters and membrane-associated proteins. Especially the presence of three P-type PM ATPases and a PM type Calmodulin-binding Ca2+-transporting ATPase is of special interest, as, e.g., activation of PM ATPases could affect the speed of recovery of the resting membrane potential following excitation of the ChRs. This might be important for, e.g., desensitization and dark recovery of the cell (Govorunova et al., 1997). Also lowering of the cytosolic free Ca2+ concentration in the eyespot region following activation of the ChRs is important for signaling and adaptation. The function of the Ca2+-ATPase might be complemented by light-induced Ca2+ uptake into the chloroplast, which has been demonstrated for isolated chloroplasts and appears to be important for the regulation of the sign of phototaxis in C. reinhardtii (Kreimer et al., 1985; Takahashi and Watanabe, 1993). Additionally, the discovery of two fascilin-like proteins with four (Cre12.g492650.t1.2) and three (Cre12.g492600.t1.2) FAS1 domains might be functionally relevant. These are thought to represent ancient cell adhesion domains. As the MORN-repeat motif, which has a critical role in several proteins with functions in the organization of membranous and cytoskeletal structures, these domains may be important for the close contact of the different membrane types in the eyespot region. A MORN-repeat protein is present in the proteomes of both eyespot fractions (Schmidt et al., 2006). It should, however, be noted that for all proteins discussed so far, a role in eyespot development, ultrastructure, positioning or eyespot related signaling needs to be proven experimentally in the future. While the extended eyespot proteome revealed numerous novel eyespot proteins, some were also “lost” compared to the already published one (Schmidt et al., 2006). Thereby, one should consider that the former bioinformatics analysis of the proteome data was based on a so-called probability score. In the current analysis, we are applying a FDR of ≤1%, which is more restrictive. On the other hand, several of the missing candidates (10 out of 27) are putative contaminants. It may also be that the different cell rupture and eyespot isolation methods yielding either eyespot fragments (Schmidt et al., 2006) or the intact structure (present study), have different contaminants associated.

In the past years, it turned out that cellular signal transduction is not only frequently mediated by phosphorylation at Tyr, Ser and Thr but also by methylation at Arg and Lys (Biggar and Li, 2015). Protein methylation is meanwhile considered as an integral part of cellular biology. Protein Lys methyltransferases (PKMTs, often containing a so-called SET domain) and protein Arg methyltransferases (PRMTs) have been well characterized, while Lys and Arg demethylases are less well studied so far or still under debate in case of Arg demethylases (Alban et al., 2014). In A. thaliana, six of the SET-domain containing PKMT’s are known or predicted to be targeted to the chloroplast (Table 1 in Alban et al., 2014). A search for homologous PKMTs in the C. reinhardtii genome, which was based on the Arabidopsis candidates, revealed five putative protein Lys methyltransferases. They have either a Rubisco-lysine N-methyltransferase domain [Cre16.g661350.t1.2 (RMT1), Cre12.g524500.t.1.2 (RMT2), and Cre12.g503800.t1.1 (RMT5)], a SET domain (Cre16.g649700.t1.1) or both (Cre12.g541777.t1.1). In case of the ATXR5 protein in Arabidopsis, which is supposed to be a histone-lysine methyltransferase, only a limited number of proteins with a very weak homology was found in the C. reinhardtii genome. None of them has a methyltransferase domain or a SET domain. Another distantly related group of PKMT’s predicted for Arabidopsis chloroplasts (Table 1 in Alban et al., 2014) belongs to the seven-beta strand methyltransferases. One is a ribosomal protein L11 methyltransferase-like protein (PrmA-like) and is also found in Chlamydomonas (Cre09.g396735.t2.1). The other is a S-adenosyl-L-methionine-dependent methyltransferase (CaMKMT-like) and has only very limited homology (E-value: 5.6E-9) to a putative N2-dimethylguanosine tRNA methyltransferase (Cre07.g347800.t1.1) in Chlamydomonas. In the case of the PRMT’s, five enzymes are known in C. reinhardtii (summarized in Werner-Peterson and Sloboda, 2013). PRMTs 1, 2, 3, and 6 are Type I enzymes, which produce asymmetric dimethyl Arg residues, whereas PRMT5 is a Type II enzyme and produce symmetric dimethyl Arg’s. While PRMT1 is known to be localized in flagella (Werner-Peterson and Sloboda, 2013), nothing is known about the cellular localization of the other proteins. In Arabidopsis (Alban et al., 2014), only PRMT7 is predicted to be localized in the chloroplast. In C. reinhardtii, the closest two homologs to this protein are an uncharacterized PRMT (Cre02.g14626.t1.1) and a Protein/Histone-arginine N-methyltransferase [Cre.g558100.t1.2 (PRMT2)]. It seems likely that the methylation events of eyespot proteins located within the chloroplast part of this structure are catalyzed by chloroplast localized methyltransferases. This may even hold true for the methylated 14-3-3 protein, since members of this protein group are also reported to be present in the chloroplast (Sehnke et al., 2000, 2002).

In this context, methylation of known eyespot proteins and proteins forming complexes with eyespot proteins is of interest. Remarkably all of these proteins are localized in chloroplast parts of the eyespot. EYE2 is localized in the chloroplast envelope and important for the formation of the EG layer (Boyd et al., 2011). The photosystem I assembly protein Ycf4, which was found to be methylated at one Lys residue, is known to be in a complex with the eyespot protein COP2 (Ozawa et al., 2009) as well as certain PSI subunits including PsaB, which was also identified as a methylated protein in our analysis. Also the ATP synthase subunits alpha and beta, which are part of the soluble CF1 subunit, are of interest. Both subunits were found as methylated proteins in the chloroplast of Arabidopsis (Alban et al., 2014) and they are also present as methylated proteins in the Chlamydomonas eyespot. Former studies (Schmidt et al., 2007) showed that the alpha and beta subunits present in the eyespot fraction of Chlamydomonas are resistant against thermolysin treatment, which is not the case for the thylakoid-localized alpha and beta subunits of the ATP synthase and other membrane-associated proteins in the eyespot fraction. Therefore, it was assumed that a significant proportion of these ATP-synthase subunits have a specialized localization and function within the eyespot, possibly between the EG. Blue native PAGE of thermolysin-treated eyespots followed by SDS-PAGE even revealed that the alpha and beta subunits are present in conjunction with the gamma-subunit in a thermolysin resistant complex (Schmidt et al., 2007). Thus, methylation of these subunits is of special interest with regard to the conservation of methylation sites compared to Arabidopsis chloroplasts but also with regard to potential special functional features of methylation within the eyespot. In Arabidopsis, the alpha subunit reveals a methylation site on Arg141. In Chlamydomonas eyespots, we detected two methylation sites in this highly conserved protein at Arg161 and Lys470. Although an Arg is present at position 161 in Arabidopsis, it is not methylated there. On the other hand, at position 141 in C. reinhardtii also a non-methylated Arg is present. The methylated Lys at position 470 is not conserved in Arabidopsis. For the beta subunit, we identified six methylated residues in C. reinhardtii (Lys96, Lys104, Lys191, Arg418, Lys447, and Lys460) compared to two sites in the A. thaliana chloroplast (Alban et al., 2014). For three positions (Lys96, Lys104, Lys460) the corresponding Lys residues are not conserved in A. thaliana. Lys191 and Arg418 are conserved, but not methylated in Arabidopsis. One of the methylation sites (Lys447), however, fits exactly to one of the two methylation sites that were identified in the ATP synthase subunit beta in the chloroplast of A. thaliana (Lys447, Arg52; Alban et al., 2014). It remains to be studied in the future whether the higher methylation status and rather distinct methylation pattern (with one exception) of these two subunits in the eyespot compared to the thylakoids is possibly related to their specialized location in the eyespot. The increased hydrophobicity without altering the charge of the methylated amino acid residues (Rice and Allis, 2001) might be beneficial for their localization between the highly hydrophobic EG.

Among the methylated proteins found in the eyespot fractions is also the SOUL heme binding protein SOUL3 that is localized at the EG and important for a correct positioning of the eyespot and its size (Schulze et al., 2013). It was shown before that SOUL3 is phosphorylated at two sites (Wagner et al., 2008). A crosstalk between methylation–phosphorylation events has been recently found (reviewed in Biggar and Li, 2015). If the methylation and phosphorylation sites are neighbored, this is called “methylation–phosphorylation switch”. Moreover, a methylation–methylation crosstalk has been described. In case of SOUL3, a “methylation–phosphorylation switch” seems rather unlikely as the relevant sites are not directly neighbored. A methylation-methylation crosstalk, however, may exist. SOUL3 methylation sites are at Arg95 and Lys102 while the phosphorylation sites are at Thr42 and Ser44. Notably, both posttranslational modifications are situated at the N-terminal part of SOUL3 while the SOUL/HBP domain is situated at the C-terminus. Methylation of SOUL3 may be involved in regulation of protein–protein interactions. For example, dimethylation of p53 promotes its association with the co-activator protein p53-binding protein 1. Lys methylation can, however, also block protein–protein interactions as demonstrated for the MAPK kinase kinase 2 and its interaction with the serine/threonine protein phosphatase 2A complex (reviewed in Hamamoto et al., 2015). SOUL3 (approximately 45 kDa) has been found to be present in complexes during the day (approximately 100 kDa) and the night phase (approximately 270 kDa) over a circadian cycle (Schulze et al., 2013). It has to be investigated whether methylation affects complex formation of SOUL3 and whether its methylation status changes over a circadian cycle.

The knowledge about the protein composition of the eyespot proteomes as well as about posttranslational modifications such as methylation and phosphorylation provides now an efficient basis for further functional studies. For example, potential effects of protein methylation on protein–protein interactions, protein stability, or subcellular localization can now be put under focus. Signaling components such as EYE2, a protein of the thioredoxin superfamily, or the 14-3-3 protein might recruit interaction partners or induce/inhibit signaling pathways upon methylation. As mentioned before for the ATPase, it is also well visible that eyespot proteins may rely on methylation to increase their hydrophobicity profile especially when situated in between the EG. These emerging functional possibilities will need further attention in future approaches in order to understand this highly complex primordial visual system of a unicellular organism in depth with regard to its structural properties and signaling pathways.

Author Contributions

GK and MM designed the research. DB, NE, VW, and WW performed the experiments. NE, SG, and GK did the EM analyses; GK, MM, and VW wrote the article with contributions and edits from all other authors.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Table 1

Functional categorization of newly identified proteins from the extended eyespot fraction F2Ae.

Transcript name (Phytozome)	No. of different peptides	Function and/or homologies	TMDsa
Proteins important for eyespot development/retina related proteins
Cre02.g105600.t1.3	4	Eyespot assembly protein EYE3, ABC1 kinase family	(+)
Cre12.g547300.t1.3	3	Eye Pigment Precursor Transporter (EPP) family protein	+
Cre12.g488350.t1.3	2	Retinal pigment epithelial membrane receptor	(+)
SOUL heme-binding proteins
Cre13.g566850.t1.2	2	Similar to SOUL2	(+)
Cre06.g292400.t1.3	2	Similar to SOUL5	(+)
14-3-3 proteins
Cre12.g559250.t1.2	10	14-3-3 protein	-
Cre06.g257500.t1.2	8	14-3-3 protein	-
Kinasesb
Cre17.g705000.t1.2	11	Calcium-dependent protein kinase 1	(+)
Cre09.g407800.t1.3	10	ABC1/COQ8 ser/thr kinase with an AarF (predicted unusual kinase) domain	(+)
Cre16.g663200.t1.3	10	Cyclic nucleotide dependent protein kinase	(+)
g8097.t1	9	Predicted protein with AarF (predicted unusual protein kinase) domain	(+)
g8097.t2c	5	Predicted protein with AarF (predicted unusual protein kinase) domain	(+)
Cre13.g570350.t1.3; Cre13.g570350.t2.1	6	ABC-1-like kinase with an AarF (predicted unusual protein kinase) domain	+
g13907.t1	4	Predicted protein with AarF (predicted unusual protein kinase) domain	(+)
Cre06.g307100.t1.3	3	ABC1/COQ8 ser/thr kinase with an AarF domain	(+)
Cre16.g657350.t1.2	3	Protein with catalytic domain of serine/threonine protein kinases	(+)
Cre03.g168150.t1.2	2	Protein with catalytic domain of tyrosine kinase	(+)
Cre13.g571700.t1.3	2	Protein with catalytic domain of serine/threonine protein kinases	(+)
g17359.t1	2	Protein with catalytic domain of serine/threonine protein kinases	(+)
Phosphatasesb
Cre06.g292550.t1.2	3	Protein phosphatase 1	-
Cre09.g388750.t1.2	3	Phosphoinositide phosphatase	+
Cre06.g257850.t1.2	2	Serine/threonine protein phosphatase	-
Cre01.g030200.t1.2	2	Protein phosphatase 2C-like protein	-
Calcium-sensing and binding proteins
Cre11.g468450.t1.2	6	Similar to centrin	-
Cre14.g615750.t1.1	4	Protein with EF-hand, calcium binding motif	(+)
Cre12.g559450.t1.3	4	Protein with a C2 domain (found in kinases and membrane trafficking proteins)	+
Cre03.g150300.t1.2	3	Protein with EF-hand, calcium binding motif	-
Cre03.g178150.t1.1	2	Similar to calmodulin	-
Cre15.g641250.t1.2	2	Protein with EF-hand, calcium binding motif	-
Membrane-associated/structural proteins, proteins with PAP-fibrillin domain
Cre13.g583550.t1.2	11	VIPP1, Vesicle inducing protein in plastids 1	-
Cre03.g197650.t1.2	8	Protein with PAP-fibrillin domain	(+)
Cre12.g502250.t1.2	4	Protein with PAP-fibrillin domain	+
Cre12.g492600.t1.2	3	Fasciclin-like protein	+
Cre12.g492650.t1.2	2	Fasciclin-like protein	+
Cre11.g478850.t1.2	2	Protein with PAP-fibrillin domain	(+)
Transporter
Cre10.g459200.t1.2	11	Plasma membrane-type proton ATPase	+
Cre04.g220200.t2.1; Cre04.g220200.t1.2; Cre04.g220200.t3.1	11	K+/H+-efflux antiporter 2 (KEA2, chloroplast inner envelope)	+
Cre09.g388850.t1.1	4	Calmodulin binding calcium-transporting ATPase (P-type/plasma membrane)	+
Cre03.g164600.t1.2;	4	Plasma membrane hydrogen ATPase	+
Cre03.g165050.t1.2	2	Plasma-membrane proton-efflux P-type ATPase	+

Table 2

Functional categorization of identified methylated proteins from the eyespot fractions F2Ae and F1 (Schmidt et al., 2006).

Transcript name (Phytozome)	Function and/or homologies	Methylated peptide	z	Xcorr	x-times found
Proteins important for eyespot development
Cre16.g666550.t1.2	SOUL3	QRQAFIMNDTCRmFLATDLKm2a	3	3.80	2
		–>QRQAFIMoNDTCRmFLATDLKm2b	3	3.51	1
Cre12.g509250.t1.1	EYE2, no eyespot	LTDDELIALVNSDPDLDKmb	2	4.40	1
14-3-3 proteins
Cre12.g559250.t1.2	14-3-3 protein	DNLTLWTSDMoQDPAAGDDRmEGADMo Km2VEDAEPa	3	4.37	2
Photosynthesis/electron transport/light harvesting
Cp genome	AtpB, ATP synthase subunit beta	FLSQPFFVAEVFTGSPGKm2YVSLAETIEG FGKa,b	3	5.97	6
		–>FLSQPFFVAEVFTGSPGKmYVSLAETIE GFGKma	3	3.81	1
		ELQDIIAILGLDELSEEDRmb	2	3.51	2
		GMEVVDTGKm2PLSVPVGKa	2	2.91	1
		–>GMoEVVDTGKPLSVPVGKmb	2	3.09	1
		–>GMoEVVDTGKmPLSVPVGKb	2	2.81	1
		TVLIMoELINNIAKmb	2	2.73	1
Cp genome	AtpA, ATP synthase subunit alpha	SYLANSYPKm2YGEILRa	2	2.93	1
		SVYEPLATGLVAVDAMoIPVGRmb	2	3.38	1
Cp genome	Atpl, CF0 ATP synthase subunit I	YVEPAAFLLPINVLEDFTKm3PLSLSFRa	3	3.11	1
Cre06.g261000.t1.2	10 kDa PS II polypeptide	GKm2GYGVYRb	1	2.26	2
		YEDKmYGANVDGYSPIYTPDLWTESGDSYTLGTKa	3	7.16	1
Cp genome	PetA, Cytochrome f	KmYSEMVVPILSPDPAKmb	2	2.60	1
		–>YSEMoVVPILSPDPAKmb	2	2.60	1
Cre16.g687900.t1.2	Lhca7, light-harvesting protein of PS I	NPGSQADGSFLGFTEEFKma	2	3.48	5
Cp genome	PsbD, PS II D2 protein	AAEDPEFETFYTKma	2	3.40	1
Cp genome	PsaB, PS I P700 chlorophyll a apoprotein A2	GYWQELIETLVWAHEKmTPLANLVYWKma	3	3.40	1
Cp genome	Ycf4, PS I assembly protein	EIEKmQASELANFLQVSLEAb	2	2.77	1
Transporter
Cp genome	CemA (Ycf10), inner envelope protein	FLKm3QLFSDVDNLVIQEYRa	3	3.10	1
		GSLDSIKmNKm3DISKa	3	3.00	1
g11711.t1	Similar to ATPase components of ABC transporters	LQTTKIGMLSEGQKm2SRa	2	3.28	1
Ferredoxin and thioredoxin-related proteins
Cre11.g476750.t1.2	Ferredoxin-NADP reductase	KmGLCSNFLCDATPGTEISMoTGPTGKa	2	3.61	1
		–>KmGLCSNFLCDATPGTEISMTGPTGKa	2	3.42	1
		IPFWEGQSYGVIPPGTKmINSKm2a	3	3.82	1
		–>IPFWEGQSYGVIPPGTKINSKm3a	3	3.84	1
(Lipid) Metabolism
g11946.t1, g11946.t2	Similar to Cytochrome b5 reductase	APDYSQGEVSGLLKm2a	2	3.32	3
Cre07.g349700.t1.2	Similar to 3-beta hydroxysteroid dehydrogenase/isomerase	ALVRDVSKmATSGSGLLAGVGSTTEVVRb	3	3.93	2
		–>ALVRmDVSKATSGSGLLAGVGSTTEVVRb	3	3.70	1
Cre01.g017100.t1.3	Similar to proteins with a acylglycerol/acyl-transferase domain	WFESFGAVKASPMAAFRm2LLRa	3	3.87	1

Table 3

Methylated proteins of unknown function from the eyespot fractions F2Ae and F1 (Schmidt et al., 2006).

Transcript name (Phytozome)	Function and/or homologies	Methylated peptide	z	Xcorr	x-times found
Cre06.g263250.t1.3	No significant hit in NCBI BLASTp	AAVADATGAASSAAADAKm2a	2	5.71	7
		AAVADATGAASSAATDAKm2a	2	4.20	3
Cp genome	ORF1995 unknown protein	MALEDLSKm3WKm3a,b	2	2.71	5
		SFDITSMTTTLPFYAGWDESLKm2a	2	4.61	3
		–>SFDITSMoTTTLPFYAGWDESLKm2a	2	3.59	1
Cre01.g000900.t1.2	Similar to conserved plant/cyanobacterial proteins of unknown functions, contains two DUF1350 domains	LATVAGQLGVSAATAPLEELSRma,b	2	4.03	4
		FKDDSLDDTNNLVQLLQGSSSVGEVLDLTVRmb	3	4.35	1
Cp genome	ORF2971 unknown protein	VAMoLAELSLSNLSAKm3LDMITDLLVIIDSVRma	3	3.31	1
		MoGQRmKmSQITLLEKma	2	2.52	1
g2947.t1	No significant hit in NCBI BLASTp	ADGAAATATTAATGVLGAGFAKmADEAAASATTAATGVLGAGFAKm2a	3	3.67	1
g14174.t1	No significant hit in NCBI BLASTp	GLGDVVGMKm3GPAAEINNGRa	2	3.33	1
Cre10.g438450.t1.3	No significant hit in NCBI BLASTp	GWGKm2LPDSGAALPAFLYKmHVLKma	2	3.13	1