Workflow for Large Scale Detection and Validation of Peptide Modifications by RPLC-LTQ-Orbitrap: Application to the Arabidopsis thaliana Leaf Proteome and an Online Modified Peptide Library

Detection of the Arabidopsis Proteome and Its Post-translational Modifications and the Nature of the Unobserved (Dark) Proteome in PeptideAtlas

This study describes a new release of the Arabidopsis thaliana PeptideAtlas proteomics resource (build 2023-10) providing protein sequence coverage, matched mass spectrometry (MS) spectra, selected post-translational modifications (PTMs), and metadata. 70 million MS/MS spectra were matched to the Araport11 annotation, identifying ∼0.6 million unique peptides and 18,267 proteins at the highest confidence level and 3396 lower confidence proteins, together representing 78.6% of the predicted proteome. Additional identified proteins not predicted in Araport11 should be considered for the next Arabidopsis genome annotation. This release identified 5198 phosphorylated proteins, 668 ubiquitinated proteins, 3050 N-terminally acetylated proteins, and 864 lysine-acetylated proteins and mapped their PTM sites. MS support was lacking for 21.4% (5896 proteins) of the predicted Araport11 proteome: the "dark" proteome. This dark proteome is highly enriched for E3 ligases, transcription factors, and for certain (e.g., CLE, IDA, PSY) but not other (e.g., THIONIN, CAP) signaling peptides families. A machine learning model trained on RNA expression data and protein properties predicts the probability that proteins will be detected. The model aids in discovery of proteins with short half-life (e.g., SIG1,3 and ERF-VII TFs) and for developing strategies to identify the missing proteins. PeptideAtlas is linked to TAIR, tracks in JBrowse, and several other community proteomics resources.

Compendium of Chromatographic Behavior of Post-translationally and Chemically Modified Peptides in Bottom-Up Proteomic Experiments

We have systematically evaluated the chromatographic behavior of post-translationally/chemically modified peptides using data spanning over 70 of the most relevant modifications. These retention properties were measured for standard bottom-up proteomic settings (fully porous C18 separation media, 0.1% formic acid as ion-pairing modifier) using collections of modified/nonmodified peptide pairs. These pairs were generated by spontaneous degradation, chemical or enzymatic treatment, analysis of synthetic peptides, or the cotranslational incorporation of noncanonical proline analogues. In addition, these measurements were validated using external data acquired for synthetic peptides and enzymatically induced citrullination. Working in units of hydrophobicity index (HI, % ACN) and evaluating the average retention shifts (ΔHI) represent the simplest approach to describe the effect of modifications from a didactic point of view. Plotting HI values for modified (y-axis) vs nonmodified (x-axis) counterparts generates unique slope and intercept values for each modification defined by the chemistry of the modifying moiety: its hydrophobicity, size, pKa of ionizable groups, and position of the altered residue. These composition-dependent correlations can be used for coarse incorporation of PTMs into models for prediction of peptide retention. More accurate predictions would require the development of specific sequence-dependent algorithms to predict ΔHI values.

Mapping the Arabidopsis thaliana proteome in PeptideAtlas and the nature of the unobserved (dark) proteome; strategies towards a complete proteome

This study describes a new release of the Arabidopsis thaliana PeptideAtlas proteomics resource providing protein sequence coverage, matched mass spectrometry (MS) spectra, selected PTMs, and metadata. 70 million MS/MS spectra were matched to the Araport11 annotation, identifying ∼0.6 million unique peptides and 18267 proteins at the highest confidence level and 3396 lower confidence proteins, together representing 78.6% of the predicted proteome. Additional identified proteins not predicted in Araport11 should be considered for building the next Arabidopsis genome annotation. This release identified 5198 phosphorylated proteins, 668 ubiquitinated proteins, 3050 N-terminally acetylated proteins and 864 lysine-acetylated proteins and mapped their PTM sites. MS support was lacking for 21.4% (5896 proteins) of the predicted Araport11 proteome - the 'dark' proteome. This dark proteome is highly enriched for certain ( e.g. CLE, CEP, IDA, PSY) but not other ( e.g. THIONIN, CAP,) signaling peptides families, E3 ligases, TFs, and other proteins with unfavorable physicochemical properties. A machine learning model trained on RNA expression data and protein properties predicts the probability for proteins to be detected. The model aids in discovery of proteins with short-half life ( e.g. SIG1,3 and ERF-VII TFs) and completing the proteome. PeptideAtlas is linked to TAIR, JBrowse, PPDB, SUBA, UniProtKB and Plant PTM Viewer.

Peptide retention time prediction for peptides with post-translational modifications: N-terminal (α-amine) and lysine (ε-amine) acetylation

Development of a peptide retention prediction model in reversed-phase chromatography is reported for acetylated peptides - both N-terminal (α-) and side chain of Lys (ε-amine) residues. Large-scale proteomic 2D LC-MS analyses of acetylated/non-acetylated tryptic digest of whole human cell lysate have been used to assemble representative retention data sets of 25,000+ modified/non-modified pairs. This allowed elucidating chromatographic behaviour of modified peptides in three different separation modes: high pH reversed-phase, HILIC separation on amide phase (first dimension of 2D) and reversed-phase separation with formic acid as ion-pairing modifier in the second dimension. On average, N-terminal acetylation increases peptide RP retention at acidic pH by 5 Hydrophobicity Index units (% acetonitrile). Acetylation of first lysine adds another 4.1%. The magnitude of the retention shift varies greatly depending on the number of modified amines, peptide length, and N-terminal peptide sequence. Large retention shifts have been observed for peptides with hydrophobic N-termini and specifically peptides carrying sequences characteristic for amphipathic helical structures - all in complete agreement with major sequence-specific features of RP retention mechanism. The utility of the modified Sequence Specific Retention Calculator model has been verified for the in-vivo N-terminally acetylated peptides detected by 2D LC-MS/MS analysis of a yeast tryptic digest. The effect of N-terminal acetylation was also evaluated for six different HILIC columns, strong cation- and strong anion exchange separations using previously acquired 2D LC-MS/MS data.

The Arabidopsis PeptideAtlas: Harnessing worldwide proteomics data to create a comprehensive community proteomics resource

We developed a resource, the Arabidopsis PeptideAtlas (www.peptideatlas.org/builds/arabidopsis/), to solve central questions about the Arabidopsis thaliana proteome, such as the significance of protein splice forms and post-translational modifications (PTMs), or simply to obtain reliable information about specific proteins. PeptideAtlas is based on published mass spectrometry (MS) data collected through ProteomeXchange and reanalyzed through a uniform processing and metadata annotation pipeline. All matched MS-derived peptide data are linked to spectral, technical, and biological metadata. Nearly 40 million out of ∼143 million MS/MS (tandem MS) spectra were matched to the reference genome Araport11, identifying ∼0.5 million unique peptides and 17,858 uniquely identified proteins (only isoform per gene) at the highest confidence level (false discovery rate 0.0004; 2 non-nested peptides ≥9 amino acid each), assigned canonical proteins, and 3,543 lower-confidence proteins. Physicochemical protein properties were evaluated for targeted identification of unobserved proteins. Additional proteins and isoforms currently not in Araport11 were identified that were generated from pseudogenes, alternative start, stops, and/or splice variants, and small Open Reading Frames; these features should be considered when updating the Arabidopsis genome. Phosphorylation can be inspected through a sophisticated PTM viewer. PeptideAtlas is integrated with community resources including TAIR, tracks in JBrowse, PPDB, and UniProtKB. Subsequent PeptideAtlas builds will incorporate millions more MS/MS data.

Confident Identification of Citrullination and Carbamylation Assisted by Peptide Retention Time Prediction

The chromatographic behavior of peptides carrying citrulline and homocitrulline residues in proteomic two-dimensional (2D) liquid chromatography-mass spectrometry (LC-MS) experiments has been investigated. The primary goal of this study was to determine the chromatographic conditions that allow differentiating between arginine citrullination and deamidation of asparagine based on retention data, improving the confidence of MS-based identifications. Carbamylation was used as a reference point due to a high degree of similarity between modification products and anticipated changes in chromatographic behavior. We applied 2D LC-MS/MS (a high-pH-low-pH reversed phase (RP), hydrophilic interaction liquid chromatography (HILIC)-low-pH RP, and strong cation exchange (SCX)-low-pH RP) to acquire retention data for modified-nonmodified peptide pairs in the four separation modes. Modifications of a standard protein mixture were induced enzymatically (PAD-2) or chemically (urea) for citrullination and carbamylation, respectively. Deamidation occurs spontaneously. Similar retention shifts were observed for all three modifications in a high-pH RP (decrease) and a low-pH RP (increase), thus limiting the applicability of this 2D LC combination. HILIC on bare silica and strong cation exchange separations have been probed to amplify the effect of charge loss upon citrullination, with SCX demonstrating the most differentiating power: the elimination of basic residues upon citrullination/carbamylation results in an ∼58 mM KCl retention decrease, while retention of deamidated products decreases slightly.

Potential Plasticity of the Mannoprotein Repertoire Associated to Mycobacterium tuberculosis Virulence Unveiled by Mass Spectrometry-Based Glycoproteomics

To date, Mycobacterium tuberculosis (Mtb) remains the world’s greatest infectious killer. The rise of multidrug-resistant strains stresses the need to identify new therapeutic targets to fight the epidemic. We previously demonstrated that bacterial protein-O-mannosylation is crucial for Mtb infectiousness, renewing the interest of the bacterial-secreted mannoproteins as potential drug-targetable virulence factors. The difficulty of inventorying the mannoprotein repertoire expressed by Mtb led us to design a stringent multi-step workflow for the reliable identification of glycosylated peptides by large-scale mass spectrometry-based proteomics. Applied to the differential analyses of glycoproteins secreted by the wild-type Mtb strain—and by its derived mutant invalidated for the protein-O-mannosylating enzyme PMTub—this approach led to the identification of not only most already known mannoproteins, but also of yet-unknown mannosylated proteins. In addition, analysis of the glycoproteome expressed by the isogenic recombinant Mtb strain overexpressing the PMTub gene revealed an unexpected mannosylation of proteins, with predicted or demonstrated functions in Mtb growth and interaction with the host cell. Since in parallel, a transient increased expression of the PMTub gene has been observed in the wild-type bacilli when infecting macrophages, our results strongly suggest that the Mtb mannoproteome may undergo adaptive regulation during infection of the host cells. Overall, our results provide deeper insights into the complexity of the repertoire of mannosylated proteins expressed by Mtb, and open the way to novel opportunities to search for still-unexploited potential therapeutic targets.

Protein methylome analysis in Arabidopsis reveals regulation in RNA-related processes

Protein methylation has been proposed as an important post-translational modification, which occurs predominantly on lysine and arginine residues. Recent discoveries have revealed that protein methylation is also present on non-histones besides histones, and plays critical roles in regulating protein stability and function. However, proteome-wide identification of methylated proteins in plants remains unexplored. Here, we present the first global survey of monomethyl arginine, symmetric and asymmetric dimethyl arginine, and monomethyl, dimethyl, trimethyl lysine modifications in the proteomes of 10-day-old Arabidopsis seedlings through a combination of immunoaffinity purification and mass spectrometry analysis. In total, we identified 617 methylation sites which mapped to 412 proteins, with 263 proteins harboring 381 lysine methylation sites and 149 proteins harboring 236 arginine methylation sites. Among them, 607 methylation sites on 408 proteins were novel findings. Motif analysis revealed that glycine preferentially flanked methylated arginine residues, whereas aspartate and glutamate enriched around mono- and dimethylated lysine sites. Methylated proteins were involved in a variety of metabolic processes, showing significant enrichment in RNA-related metabolic pathways including spliceosome, RNA transport, and ribosome. Our data provide a global view of methylated non-histone proteins in Arabidopsis, laying foundations for elucidating the biological function of protein methylation in plants. SIGNIFICANCE: Protein methylation has emerged as a common and important modification both in eukaryotes and prokaryotes. The identification of methylated sites/peptides is fundamental for further functional analysis of protein methylation. This study was the first proteome-scale identification of lysine and arginine methylation in plants. We found that methylation occurred widely on non-histone proteins in Arabidopsis and was involved in diverse biological functions. The results provide foundations for the investigation of the protein methylome in Arabidopsis and provide powerful resources for the functional analysis of protein methylation in plants.

Selective Sensing of Phosphorylated Peptides and Monitoring Kinase and Phosphatase Activity with a Supramolecular Tandem Assay

Simple tuning of a host:guest pair allows selective sensing of different peptide modifications, exploiting orthogonal recognition mechanisms. Excellent selectivity for either lysine trimethylations or alcohol phosphorylations is possible by simply varying the fluorophore guest. The phosphorylation sensor can be modulated by the presence of small (μM) concentrations of metal ions, allowing array-based sensing. Phosphorylation at serine, threonine, and tyrosine can be selectively sensed via discriminant analysis. The phosphopeptide sensing is effective in the presence of small-molecule phosphates such as ATP, which in turn enables the sensor to be employed in continuous optical assays of both serine kinase and tyrosine phosphatase activity. The activity of multiple different kinases can be monitored, and the sensor is capable of detecting the phosphorylation of peptides containing multiple different modifications, including lysine methylations and acetylation. A single deep cavitand can be used as a "one size fits all" sensor that can selectively detect multiple different modifications to oligopeptides, as well as monitoring the function of their post-translational modification writer and eraser enzymes in complex systems.

Cold and Heat Stress Diversely Alter Both Cauliflower Respiration and Distinct Mitochondrial Proteins Including OXPHOS Components and Matrix Enzymes

Complex proteomic and physiological approaches for studying cold and heat stress responses in plant mitochondria are still limited. Variations in the mitochondrial proteome of cauliflower (Brassica oleracea var. botrytis) curds after cold and heat and after stress recovery were assayed by two-dimensional polyacrylamide gel electrophoresis (2D PAGE) in relation to mRNA abundance and respiratory parameters. Quantitative analysis of the mitochondrial proteome revealed numerous stress-affected protein spots. In cold, major downregulations in the level of photorespiratory enzymes, porine isoforms, oxidative phosphorylation (OXPHOS) and some low-abundant proteins were observed. In contrast, carbohydrate metabolism enzymes, heat-shock proteins, translation, protein import, and OXPHOS components were involved in heat response and recovery. Several transcriptomic and metabolic regulation mechanisms are also suggested. Cauliflower plants appeared less susceptible to heat; closed stomata in heat stress resulted in moderate photosynthetic, but only minor respiratory impairments, however, photosystem II performance was unaffected. Decreased photorespiration corresponded with proteomic alterations in cold. Our results show that cold and heat stress not only operate in diverse modes (exemplified by cold-specific accumulation of some heat shock proteins), but exert some associations at molecular and physiological levels. This implies a more complex model of action of investigated stresses on plant mitochondria.

Separation of Methylated Histone Peptides via Host-Assisted Capillary Electrophoresis

Lysine methylation in protein is one important epigenetic mechanism that regulates diverse biological processes but is challenging to study due to the large variability in methylation levels and sites. Here, we show that supramolecular hosts such as calixarenes and cucurbiturils can be applied in the background electrolyte (BGE) of capillary electrophoresis (CE) for highly effective separation of post-translationally methylated histone peptides. The molecular recognition event causes a shift in the electrophoretic mobility of the peptide, allowing affinity measurement for binding between the synthetic receptor and various methylated lysine species. Successful separation of the H3 peptides carrying different methylation levels at the K9 position can be achieved using CX4 and CX6 as the BGE additives in CE, enabling monitoring of the activity of the histone lysine demethylase JMJD2E. This reveals the power of combining high resolution CE with synthetic hosts for study of protein methylation, and the method should be capable of analyzing complex biological samples for better understanding of the functions of histone methylation.

Proteins in aggregates functionally impact multiple neurodegenerative disease models by forming proteasome-blocking complexes

Age-dependent neurodegenerative diseases progressively form aggregates containing both shared components (e.g., TDP-43, phosphorylated tau) and proteins specific to each disease. We investigated whether diverse neuropathies might have additional aggregation-prone proteins in common, discoverable by proteomics. Caenorhabditis elegans expressing unc-54p/Q40::YFP, a model of polyglutamine array diseases such as Huntington's, accrues aggregates in muscle 2–6 days posthatch. These foci, isolated on antibody-coupled magnetic beads, were characterized by high-resolution mass spectrometry. Three Q40::YFP-associated proteins were inferred to promote aggregation and cytotoxicity, traits reduced or delayed by their RNA interference knockdown. These RNAi treatments also retarded aggregation/cytotoxicity in Alzheimer's disease models, nematodes with muscle or pan-neuronal Aβ_1–42 expression and behavioral phenotypes. The most abundant aggregated proteins are glutamine/asparagine-rich, favoring hydrophobic interactions with other random-coil domains. A particularly potent modulator of aggregation, CRAM-1/HYPK, contributed < 1% of protein aggregate peptides, yet its knockdown reduced Q40::YFP aggregates 72–86% (P < 10⁻⁶). In worms expressing Aβ_1–42, knockdown of cram-1 reduced β-amyloid 60% (P < 0.002) and slowed age-dependent paralysis > 30% (P < 10⁻⁶). In wild-type worms, cram-1 knockdown reduced aggregation and extended lifespan, but impaired early reproduction. Protection against seeded aggregates requires proteasome function, implying that normal CRAM-1 levels promote aggregation by interfering with proteasomal degradation of misfolded proteins. Molecular dynamic modeling predicts spontaneous and stable interactions of CRAM-1 (or human orthologs) with ubiquitin, and we verified that CRAM-1 reduces degradation of a tagged-ubiquitin reporter. We propose that CRAM-1 exemplifies a class of primitive chaperones that are initially protective and highly beneficial for early reproduction, but ultimately impair aggregate clearance and limit longevity.

Global analysis of lysine acetylation suggests the involvement of protein acetylation in diverse biological processes in rice (Oryza sativa)

Lysine acetylation is a reversible, dynamic protein modification regulated by lysine acetyltransferases and deacetylases. Recent advances in high-throughput proteomics have greatly contributed to the success of global analysis of lysine acetylation. A large number of proteins of diverse biological functions have been shown to be acetylated in several reports in human cells, E.coli, and dicot plants. However, the extent of lysine acetylation in non-histone proteins remains largely unknown in monocots, particularly in the cereal crops. Here we report the mass spectrometric examination of lysine acetylation in rice (Oryza sativa). We identified 60 lysine acetylated sites on 44 proteins of diverse biological functions. Immunoblot studies further validated the presence of a large number of acetylated non-histone proteins. Examination of the amino acid composition revealed substantial amino acid bias around the acetylation sites and the amino acid preference is conserved among different organisms. Gene ontology analysis demonstrates that lysine acetylation occurs in diverse cytoplasmic, chloroplast and mitochondrial proteins in addition to the histone modifications. Our results suggest that lysine acetylation might constitute a regulatory mechanism for many proteins, including both histones and non-histone proteins of diverse biological functions.

Modified Clp protease complex in the ClpP3 null mutant and consequences for chloroplast development and function in Arabidopsis

The plastid ClpPRT protease consists of two heptameric rings of ClpP1/ClpR1/ClpR2/ClpR3/ClpR4 (the R-ring) and ClpP3/ClpP4/ClpP5/ClpP6 (the P-ring) and peripherally associated ClpT1/ClpT2 subunits. Here, we address the contributions of ClpP3 and ClpP4 to ClpPRT core organization and function in Arabidopsis (Arabidopsis thaliana). ClpP4 is strictly required for embryogenesis, similar to ClpP5. In contrast, loss of ClpP3 (clpp3-1) leads to arrest at the hypocotyl stage; this developmental arrest can be removed by supplementation with sucrose or glucose. Heterotrophically grown clpp3-1 can be transferred to soil and generate viable seed, which is surprising, since we previously showed that CLPR2 and CLPR4 null alleles are always sterile and die on soil. Based on native gels and mass spectrometry-based quantification, we show that despite the loss of ClpP3, modified ClpPR core(s) could be formed, albeit at strongly reduced levels. A large portion of ClpPR subunits accumulated in heptameric rings, with overaccumulation of ClpP1/ClpP5/ClpP6 and ClpR3. Remarkably, the association of ClpT1 to the modified Clp core was unchanged. Large-scale quantitative proteomics assays of clpp3-1 showed a 50% loss of photosynthetic capacity and the up-regulation of plastoglobules and all chloroplast stromal chaperone systems. Specific chloroplast proteases were significantly up-regulated, whereas the major thylakoid protease (FtsH1/FtsH2/FtsH5/FtsH8) was clearly unchanged, indicating a controlled protease network response. clpp3-1 showed a systematic decrease of chloroplast-encoded proteins that are part of the photosynthetic apparatus but not of chloroplast-encoded proteins with other functions. Candidate substrates and an explanation for the differential phenotypes between the CLPP3, CLPP4, and CLPP5 null mutants are discussed.

Defining the core proteome of the chloroplast envelope membranes

High-throughput protein localization studies require multiple strategies. Mass spectrometric analysis of defined cellular fractions is one of the complementary approaches to a diverse array of cell biological methods. In recent years, the protein content of different cellular (sub-)compartments was approached. Despite of all the efforts made, the analysis of membrane fractions remains difficult, in that the dissection of the proteomes of the envelope membranes of chloroplasts or mitochondria is often not reliable because sample purity is not always warranted. Moreover, proteomic studies are often restricted to single (model) species, and therefore limited in respect to differential individual evolution. In this study we analyzed the chloroplast envelope proteomes of different plant species, namely, the individual proteomes of inner and outer envelope (OE) membrane of Pisum sativum and the mixed envelope proteomes of Arabidopsis thaliana and Medicago sativa. The analysis of all three species yielded 341 identified proteins in total, 247 of them being unique. 39 proteins were genuine envelope proteins found in at least two species. Based on this and previous envelope studies we defined the core envelope proteome of chloroplasts. Comparing the general overlap of the available six independent studies (including ours) revealed only a number of 27 envelope proteins. Depending on the stringency of applied selection criteria we found 231 envelope proteins, while less stringent criteria increases this number to 649 putative envelope proteins. Based on the latter we provide a map of the outer and inner envelope core proteome, which includes many yet uncharacterized proteins predicted to be involved in transport, signaling, and response. Furthermore, a foundation for the functional characterization of yet unidentified functions of the inner and OE for further analyses is provided.

Defining the boundaries and characterizing the landscape of functional genome expression in vascular tissues of Populus using shotgun proteomics

Current state-of-the-art experimental and computational proteomic approaches were integrated to obtain a comprehensive protein profile of Populus vascular tissue. This featured: (1) a large sample set consisting of two genotypes grown under normal and tension stress conditions, (2) bioinformatics clustering to effectively handle gene duplication, and (3) an informatics approach to track and identify single amino acid polymorphisms (SAAPs). By applying a clustering algorithm to the Populus database, the number of protein entries decreased from 64,689 proteins to a total of 43,069 protein groups, thereby reducing 7505 identified proteins to a total of 4226 protein groups, in which 2016 were singletons. This reduction implies that ∼50% of the measured proteins shared extensive sequence homology. Using conservative search criteria, we were able to identify 1354 peptides containing a SAAP and 201 peptides that become tryptic due to a K or R substitution. These newly identified peptides correspond to 502 proteins, including 97 previously unidentified proteins. In total, the integration of deep proteome measurements on an extensive sample set with protein clustering and peptide sequence variants provided an exceptional level of proteome characterization for Populus, allowing us to spatially resolve the vascular tissue proteome.

Subunit stoichiometry, evolution, and functional implications of an asymmetric plant plastid ClpP/R protease complex in Arabidopsis

The caseinolytic protease (Clp) protease system has been expanded in plant plastids compared with its prokaryotic progenitors. The plastid Clp core protease consists of five different proteolytic ClpP proteins and four different noncatalytic ClpR proteins, with each present in one or more copies and organized in two heptameric rings. We determined the exact subunit composition and stoichiometry for the intact core and each ring. The chloroplast ClpP/R protease was affinity purified from clpr4 and clpp3 Arabidopsis thaliana null mutants complemented with C-terminal StrepII-tagged versions of CLPR4 and CLPP3, respectively. The subunit stoichiometry was determined by mass spectrometry-based absolute quantification using stable isotope-labeled proteotypic peptides generated from a synthetic gene. One heptameric ring contained ClpP3,4,5,6 in a 1:2:3:1 ratio. The other ring contained ClpP1 and ClpR1,2,3,4 in a 3:1:1:1:1 ratio, resulting in only three catalytic sites. These ClpP1/R1-4 proteins are most closely related to the two subunits of the cyanobacterial P3/R complex and the identical P:R ratio suggests conserved adaptation. Furthermore, the plant-specific C-terminal extensions of the ClpP/R subunits were not proteolytically removed upon assembly, suggesting a regulatory role in Clp chaperone interaction. These results will now allow testing ClpP/R structure-function relationships using rationale design. The quantification workflow we have designed is applicable to other protein complexes.

The workflow for quantitative proteome analysis of chloroplast development and differentiation, chloroplast mutants, and protein interactions by spectral counting

This chapter outlines a quantitative proteomics workflow using a label-free spectral counting technique. The workflow has been tested on different aspects of chloroplast biology in maize and Arabidopsis, including chloroplast mutant analysis, cell-type specific chloroplast differentiation, and the proplastid-to-chloroplast transition. The workflow involves one-dimensional SDS-PAGE of the proteomes of leaves or chloroplast subfractions, tryptic digestions, online LC-MS/MS using a mass spectrometer with high mass accuracy and duty cycle, followed by semiautomatic data processing. The bioinformatics analysis can effectively select best gene models and deals with quantification of closely related proteins; the workflow avoids overidentification of proteins and results in more accurate protein quantification. The final output includes pairwise comparative quantitative analysis, as well as hierarchical clustering for discovery of temporal and spatial patterns of protein accumulation. A brief discussion about potential pitfalls, as well as the advantages and disadvantages of spectral counting, is provided.

Structural and metabolic transitions of C4 leaf development and differentiation defined by microscopy and quantitative proteomics in maize

C(4) grasses, such as maize (Zea mays), have high photosynthetic efficiency through combined biochemical and structural adaptations. C(4) photosynthesis is established along the developmental axis of the leaf blade, leading from an undifferentiated leaf base just above the ligule into highly specialized mesophyll cells (MCs) and bundle sheath cells (BSCs) at the tip. To resolve the kinetics of maize leaf development and C(4) differentiation and to obtain a systems-level understanding of maize leaf formation, the accumulation profiles of proteomes of the leaf and the isolated BSCs with their vascular bundle along the developmental gradient were determined using large-scale mass spectrometry. This was complemented by extensive qualitative and quantitative microscopy analysis of structural features (e.g., Kranz anatomy, plasmodesmata, cell wall, and organelles). More than 4300 proteins were identified and functionally annotated. Developmental protein accumulation profiles and hierarchical cluster analysis then determined the kinetics of organelle biogenesis, formation of cellular structures, metabolism, and coexpression patterns. Two main expression clusters were observed, each divided in subclusters, suggesting that a limited number of developmental regulatory networks organize concerted protein accumulation along the leaf gradient. The coexpression with BSC and MC markers provided strong candidates for further analysis of C(4) specialization, in particular transporters and biogenesis factors. Based on the integrated information, we describe five developmental transitions that provide a conceptual and practical template for further analysis. An online protein expression viewer is provided through the Plant Proteome Database.

Arabidopsis stromal 70-kDa heat shock proteins are essential for chloroplast development

70 kDa heat shock proteins (Hsp70s) act as molecular chaperones involved in essential cellular processes such as protein folding and protein transport across membranes. They also play a role in the cell's response to a wide range of stress conditions. The Arabidopsis family of Hsp70s homologues includes two highly conserved proteins, cpHsc70-1 and cpHsc70-2 which are both imported into chloroplasts (Su and Li in Plant Physiol 146:1231-1241, 2008). Here, we demonstrate that YFP-fusion proteins of both cpHsc70-1 and cpHsc70-2 are predominantly stromal, though low levels were detected in the thylakoid membrane. Both genes are ubiquitously expressed at high levels in both seedlings and adult plants. We further show that both cpHsc70-1 and cpHsc70-2 harbour ATPase activity which is essential for Hsp70 chaperone activity. A previously described T-DNA insertion line for cpHsc70-1 (DeltacpHsc70-1) has variegated cotyledons, malformed leaves, growth retardation, impaired root growth and sensitivity to heat shock treatment. In addition, under stress conditions, this mutant also exhibits unusual sepals, and malformed flowers and sucrose concentrations as low as 1% significantly impair growth. cpHsc70-1/cpHsc70-2 double-mutants are lethal. However, we demonstrate through co-suppression and artificial microRNA (amiRNA) approaches that transgenic plants with severely reduced levels of both genes have a white and stunted phenotype. Interestingly, chloroplasts in these plants have an unusual morphology and contain few or no thylakoid membranes. Our data show that cpHsc70-1 and cpHsc70-2 are essential ATPases, have overlapping roles and are required for normal plastid structure.

Megadalton complexes in the chloroplast stroma of Arabidopsis thaliana characterized by size exclusion chromatography, mass spectrometry, and hierarchical clustering

To characterize MDa-sized macromolecular chloroplast stroma protein assemblies and to extend coverage of the chloroplast stroma proteome, we fractionated soluble chloroplast stroma in the non-denatured state by size exclusion chromatography with a size separation range up to approximately 5 MDa. To maximize protein complex stability and resolution of megadalton complexes, ionic strength and composition were optimized. Subsequent high accuracy tandem mass spectrometry analysis (LTQ-Orbitrap) identified 1081 proteins across the complete native mass range. Protein complexes and assembly states above 0.8 MDa were resolved using hierarchical clustering, and protein heat maps were generated from normalized protein spectral counts for each of the size exclusion chromatography fractions; this complemented previous analysis of stromal complexes up to 0.8 MDa (Peltier, J. B., Cai, Y., Sun, Q., Zabrouskov, V., Giacomelli, L., Rudella, A., Ytterberg, A. J., Rutschow, H., and van Wijk, K. J. (2006) The oligomeric stromal proteome of Arabidopsis thaliana chloroplasts. Mol. Cell. Proteomics 5, 114-133). This combined experimental and bioinformatics analyses resolved chloroplast ribosomes in different assembly and functional states (e.g. 30, 50, and 70 S), which enabled the identification of plastid homologues of prokaryotic ribosome assembly factors as well as proteins involved in co-translational modifications, targeting, and folding. The roles of these ribosome-associating proteins will be discussed. Known RNA splice factors (e.g. CAF1/WTF1/RNC1) as well as uncharacterized proteins with RNA-binding domains (pentatricopeptide repeat, RNA recognition motif, and chloroplast ribosome maturation), RNases, and DEAD box helicases were found in various sized complexes. Chloroplast DNA (>3 MDa) was found in association with the complete heteromeric plastid-encoded DNA polymerase complex, and a dozen other DNA-binding proteins, e.g. DNA gyrase, topoisomerase, and various DNA repair enzymes. The heteromeric >or=5-MDa pyruvate dehydrogenase complex and the 0.8-1-MDa acetyl-CoA carboxylase complex associated with uncharacterized biotin carboxyl carrier domain proteins constitute the entry point to fatty acid metabolism in leaves; we suggest that their large size relates to the need for metabolic channeling. Protein annotations and identification data are available through the Plant Proteomics Database, and mass spectrometry data are available through Proteomics Identifications database.