Proteomics, phylogenetics, and coexpression analyses indicate novel interactions in the plastid CLP chaperone-protease system.

The chloroplast chaperone CLPC1 unfolds and delivers substrates to the stromal CLPPRT protease complex for degradation. We previously used an in vivo trapping approach to identify interactors with CLPC1 in Arabidopsis thaliana by expressing a STREPII-tagged copy of CLPC1 mutated in its Walker B domains (CLPC1-TRAP) followed by affinity purification and mass spectrometry. To create a larger pool of candidate substrates, adaptors, or regulators, we carried out a far more sensitive and comprehensive in vivo protein trapping analysis. We identified 59 highly enriched CLPC1 protein interactors, in particular proteins belonging to families of unknown functions (DUF760, DUF179, DUF3143, UVR-DUF151, HugZ/DUF2470), as well as the UVR domain proteins EXE1 and EXE2 implicated in singlet oxygen damage and signaling. Phylogenetic and functional domain analyses identified other members of these families that appear to localize (nearly) exclusively to plastids. In addition, several of these DUF proteins are of very low abundance as determined through the Arabidopsis PeptideAtlas http://www.peptideatlas.org/builds/arabidopsis/ showing that enrichment in the CLPC1-TRAP was extremely selective. Evolutionary rate covariation indicated that the HugZ/DUF2470 family coevolved with the plastid CLP machinery suggesting functional and/or physical interactions. Finally, mRNA-based coexpression networks showed that all 12 CLP protease subunits tightly coexpressed as a single cluster with deep connections to DUF760-3. Coexpression modules for other trapped proteins suggested specific functions in biological processes, e.g., UVR2 and UVR3 were associated with extraplastidic degradation, whereas DUF760-6 is likely involved in senescence. This study provides a strong foundation for discovery of substrate selection by the chloroplast CLP protease system.

The CLP and PREP protease systems coordinate maturation and degradation of the chloroplast proteome in Arabidopsis thaliana.

A network of peptidases governs proteostasis in plant chloroplasts and mitochondria. This study reveals strong genetic and functional interactions in Arabidopsis between the chloroplast stromal CLP chaperone-protease system and the PREP1,2 peptidases, which are dually localized to chloroplast stroma and the mitochondrial matrix. Higher order mutants defective in CLP or PREP proteins were generated and analyzed by quantitative proteomics and N-terminal proteomics (terminal amine isotopic labeling of substrates (TAILS)). Strong synergistic interactions were observed between the CLP protease system (clpr1-2, clpr2-1, clpc1-1, clpt1, clpt2) and both PREP homologs (prep1, prep2) resulting in embryo lethality or growth and developmental phenotypes. Synergistic interactions were observed even when only one of the PREP proteins was lacking, suggesting that PREP1 and PREP2 have divergent substrates. Proteome phenotypes were driven by the loss of CLP protease capacity, with little impact from the PREP peptidases. Chloroplast N-terminal proteomes showed that many nuclear encoded chloroplast proteins have alternatively processed N-termini in prep1prep2, clpt1clpt2 and prep1prep2clpt1clpt2. Loss of chloroplast protease capacity interferes with stromal processing peptidase (SPP) activity due to folding stress and low levels of accumulated cleaved cTP fragments. PREP1,2 proteolysis of cleaved cTPs is complemented by unknown proteases. A model for CLP and PREP activity within a hierarchical chloroplast proteolysis network is proposed.

Autocatalytic Processing and Substrate Specificity of Arabidopsis Chloroplast Glutamyl Peptidase.

Chloroplast proteostasis is governed by a network of peptidases. As a part of this network, we show that Arabidopsis (Arabidopsis thaliana) chloroplast glutamyl peptidase (CGEP) is a homo-oligomeric stromal Ser-type (S9D) peptidase with both exo- and endo-peptidase activity. Arabidopsis CGEP null mutant alleles (cgep) had no visible phenotype but showed strong genetic interactions with stromal CLP protease system mutants, resulting in reduced growth. Loss of CGEP upregulated the chloroplast protein chaperone machinery and 70S ribosomal proteins, but other parts of the proteostasis network were unaffected. Both comparative proteomics and mRNA-based coexpression analyses strongly suggested that the function of CGEP is at least partly involved in starch metabolism regulation. Recombinant CGEP degraded peptides and proteins smaller than ∼25 kD. CGEP specifically cleaved substrates on the C-terminal side of Glu irrespective of neighboring residues, as shown using peptide libraries incubated with recombinant CGEP and mass spectrometry. CGEP was shown to undergo autocatalytic C-terminal cleavage at E946, removing 15 residues, both in vitro and in vivo. A conserved motif (A[S/T]GGG[N/G]PE946) immediately upstream of E946 was identified in dicotyledons, but not monocotyledons. Structural modeling suggested that C-terminal processing increases the upper substrate size limit by improving catalytic cavity access. In vivo complementation with catalytically inactive CGEP-S781R or a CGEP variant with an unprocessed C-terminus in a cgep clpr2-1 background was used to demonstrate the physiological importance of both CGEP peptidase activity and its autocatalytic processing. CGEP homologs of photosynthetic and nonphotosynthetic bacteria lack the C-terminal prosequence, suggesting it is a recent functional adaptation in plants.

Tissue-type specific accumulation of the plastoglobular proteome, transcriptional networks, and plastoglobular functions.

Plastoglobules are dynamic protein-lipid microcompartments in plastids enriched for isoprenoid-derived metabolites. Chloroplast plastoglobules support formation, remodeling, and controlled dismantling of thylakoids during developmental transitions and environmental responses. However, the specific molecular functions of most plastoglobule proteins are still poorly understood. This review harnesses recent co-mRNA expression data from combined microarray and RNA-seq information in ATTED-II on an updated inventory of 34 PG proteins, as well as proteomics data across 30 Arabidopsis tissue types from ATHENA. Hierarchical clustering based on relative abundance for the plastoglobule proteins across non-photosynthetic and photosynthetic tissue types showed their coordinated protein accumulation across Arabidopsis parts, tissue types, development, and senescence. Evaluation of mRNA-based forced networks at different coefficient thresholds identified a central hub with seven plastoglobule proteins and four peripheral modules. Enrichment of specific nuclear transcription factors (e.g. Golden2-like) and support for crosstalk between plastoglobules and the plastid gene expression was observed, and specific ABC1 kinases appear part of a light signaling network. Examples of other specific findings are that FBN7b is involved with upstream steps of tetrapyrrole biosynthesis and that ABC1K9 is involved in starch metabolism. This review provides new insights into the functions of plastoglobule proteins and an improved framework for experimental studies.

Targeted Profiling of Arabidopsis thaliana Subproteomes Illuminates Co- and Posttranslationally N-Terminal Myristoylated Proteins.

N-terminal myristoylation, a major eukaryotic protein lipid modification, is difficult to detect in vivo and challenging to predict in silico. We developed a proteomics strategy involving subfractionation of cellular membranes, combined with separation of hydrophobic peptides by mass spectrometry-coupled liquid chromatography to identify the Arabidopsis thaliana myristoylated proteome. This approach identified a starting pool of 8837 proteins in all analyzed cellular fractions, comprising 32% of the Arabidopsis proteome. Of these, 906 proteins contain an N-terminal Gly at position 2, a prerequisite for myristoylation, and 214 belong to the predicted myristoylome (comprising 51% of the predicted myristoylome of 421 proteins). We further show direct evidence of myristoylation in 72 proteins; 18 of these myristoylated proteins were not previously predicted. We found one myristoylation site downstream of a predicted initiation codon, indicating that posttranslational myristoylation occurs in plants. Over half of the identified proteins could be quantified and assigned to a subcellular compartment. Hierarchical clustering of protein accumulation combined with myristoylation and S-acylation data revealed that N-terminal double acylation influences redirection to the plasma membrane. In a few cases, MYR function extended beyond simple membrane association. This study identified hundreds of N-acylated proteins for which lipid modifications could control protein localization and expand protein function.

MET1 is a thylakoid-associated TPR protein involved in photosystem II supercomplex formation and repair in Arabidopsis.

Photosystem II (PSII) requires constant disassembly and reassembly to accommodate replacement of the D1 protein. Here, we characterize Arabidopsis thaliana MET1, a PSII assembly factor with PDZ and TPR domains. The maize (Zea mays) MET1 homolog is enriched in mesophyll chloroplasts compared with bundle sheath chloroplasts, and MET1 mRNA and protein levels increase during leaf development concomitant with the thylakoid machinery. MET1 is conserved in C3 and C4 plants and green algae but is not found in prokaryotes. Arabidopsis MET1 is a peripheral thylakoid protein enriched in stroma lamellae and is also present in grana. Split-ubiquitin assays and coimmunoprecipitations showed interaction of MET1 with stromal loops of PSII core components CP43 and CP47. From native gels, we inferred that MET1 associates with PSII subcomplexes formed during the PSII repair cycle. When grown under fluctuating light intensities, the Arabidopsis MET1 null mutant (met1) showed conditional reduced growth, near complete blockage in PSII supercomplex formation, and concomitant increase of unassembled CP43. Growth of met1 in high light resulted in loss of PSII supercomplexes and accelerated D1 degradation. We propose that MET1 functions as a CP43/CP47 chaperone on the stromal side of the membrane during PSII assembly and repair. This function is consistent with the observed differential MET1 accumulation across dimorphic maize chloroplasts.

Discovery of a Unique Clp Component, ClpF, in Chloroplasts: A Proposed Binary ClpF-ClpS1 Adaptor Complex Functions in Substrate Recognition and Delivery.

Clp proteases are found in prokaryotes, mitochondria, and plastids where they play crucial roles in maintaining protein homeostasis (proteostasis). The plant plastid Clp machinery comprises a hetero-oligomeric ClpPRT proteolytic core, ATP-dependent chaperones ClpC and ClpD, and an adaptor protein, ClpS1. ClpS1 selects substrates to the ClpPR protease-ClpC chaperone complex for degradation, but the underlying substrate recognition and delivery mechanisms are currently unclear. Here, we characterize a ClpS1-interacting protein in Arabidopsis thaliana, ClpF, which can interact with the Clp substrate glutamyl-tRNA reductase. ClpF and ClpS1 mutually stimulate their association with ClpC. ClpF, which is only found in photosynthetic eukaryotes, contains bacterial uvrB/C and YccV protein domains and a unique N-terminal domain. We propose a testable model in which ClpS1 and ClpF form a binary adaptor for selective substrate recognition and delivery to ClpC, reflecting an evolutionary adaptation of the Clp system to the plastid proteome.

Structures, Functions, and Interactions of ClpT1 and ClpT2 in the Clp Protease System of Arabidopsis Chloroplasts.

Plastid ClpT1 and ClpT2 are plant-specific proteins that associate with the ClpPR protease. However, their physiological significance and structures are not understood. Arabidopsis thaliana loss-of-function single clpt1 and clpt2 mutants showed no visible phenotypes, whereas clpt1 clpt2 double mutants showed delayed development, reduced plant growth, and virescent, serrated leaves but were viable and produced seed. The clpt1 and clpt1 clpt2 mutants showed partial destabilization of the ClpPR complex, whereas clpt2 null mutants showed only marginal destabilization. Comparative proteomics of clpt1 clpt2 plants showed a proteostasis phenotype similar to viable mutants in ClpPR core subunits, indicating reduced Clp protease capacity. In vivo and in vitro assays showed that ClpT1 and ClpT2 can independently interact with the single ClpP ring and ClpPR core, but not with the single ClpR ring. We determined ClpT1 and ClpT2 structures (2.4- and 2.0-Å resolution) and detailed the similarities to the N-domains of bacterial ClpA/C chaperones. The ClpT structures suggested a conserved MYFF motif for interaction with the ClpPR core near the interface between the P- and R-rings. In vivo complementation showed that ClpT function and ClpPR core stabilization require the MYFF motif. Several models are presented that may explain how ClpT1,2 contribute to ClpPR protease activity.

ClpS1 is a conserved substrate selector for the chloroplast Clp protease system in Arabidopsis.

Whereas the plastid caseinolytic peptidase (Clp) P protease system is essential for plant development, substrates and substrate selection mechanisms are unknown. Bacterial ClpS is involved in N-degron substrate selection and delivery to the ClpAP protease. Through phylogenetic analysis, we show that all angiosperms contain ClpS1 and some species also contain ClpS1-like protein(s). In silico analysis suggests that ClpS1 is the functional homolog of bacterial ClpS. We show that Arabidopsis thaliana ClpS1 interacts with plastid ClpC1,2 chaperones. The Arabidopsis ClpS1 null mutant (clps1) lacks a visible phenotype, and no genetic interactions with ClpC/D chaperone or ClpPR core mutants were observed. However, clps1, but not clpc1-1, has increased sensitivity to the translational elongation inhibitor chloramphenicol suggesting a link between translational capacity and ClpS1. Moreover, ClpS1 was upregulated in clpc1-1, and quantitative proteomics of clps1, clpc1, and clps1 clpc1 showed specific molecular phenotypes attributed to loss of ClpC1 or ClpS1. In particular, clps1 showed alteration of the tetrapyrrole pathway. Affinity purification identified eight candidate ClpS1 substrates, including plastid DNA repair proteins and Glu tRNA reductase, which is a control point for tetrapyrrole synthesis. ClpS1 interaction with five substrates strictly depended on two conserved ClpS1 residues involved in N-degron recognition. ClpS1 function, substrates, and substrate recognition mechanisms are discussed.

Loss of plastoglobule kinases ABC1K1 and ABC1K3 causes conditional degreening, modified prenyl-lipids, and recruitment of the jasmonic acid pathway.

Plastoglobules (PGs) are plastid lipid-protein particles. This study examines the function of PG-localized kinases ABC1K1 and ABC1K3 in Arabidopsis thaliana. Several lines of evidence suggested that ABC1K1 and ABC1K3 form a protein complex. Null mutants for both genes (abc1k1 and abc1k3) and the double mutant (k1 k3) displayed rapid chlorosis upon high light stress. Also, k1 k3 showed a slower, but irreversible, senescence-like phenotype during moderate light stress that was phenocopied by drought and nitrogen limitation, but not cold stress. This senescence-like phenotype involved degradation of the photosystem II core and upregulation of chlorophyll degradation. The senescence-like phenotype was independent of the EXECUTER pathway that mediates genetically controlled cell death from the chloroplast and correlated with increased levels of the singlet oxygen-derived carotenoid ß-cyclocitral, a retrograde plastid signal. Total PG volume increased during light stress in wild type and k1 k3 plants, but with different size distributions. Isolated PGs from k1 k3 showed a modified prenyl-lipid composition, suggesting reduced activity of PG-localized tocopherol cyclase (VTE1), and was consistent with loss of carotenoid cleavage dioxygenase 4. Plastid jasmonate biosynthesis enzymes were recruited to the k1 k3 PGs but not wild-type PGs, while pheophytinase, which is involved in chlorophyll degradation, was induced in k1 k3 and not wild-type plants and was localized to PGs. Thus, the ABC1K1/3 complex contributes to PG function in prenyl-lipid metabolism, stress response, and thylakoid remodeling.

Correlation of mRNA and protein abundance in the developing maize leaf.

To help understand regulation of maize leaf blade development, including sink-source transitions and induction of C4 photosynthesis, we compared large-scale quantitative proteome and transcriptomes collected at specific stages along the developmental maize leaf blade gradient. Proteome data were based on label-free shotgun proteomics (spectral counting) and transcript data were based on RNA-seq using the same source materials, and had been published previously (Nat Genet, 42, 2010, 1060-1067; The Plant Cell, 22, 2010, 3509-3542). Transcript and protein abundance followed near normal distributions, in contrast with several studies with other organisms. Protein observability correlated with transcript abundance following a 'lazy step function' similar to that in bacteria and yeast. mRNA and protein abundance showed significant positive correlations (up to 0.8) for log-transformed length-weighted normalized spectral abundance factor (NSAF) and reads per kilobase of exon model per million mapped reads (RPKM) and non-weighted abundances (NadjSPC and COV) in dependence of function and development. Correlations were much weaker in the leaf 'sink-source' transition zone, i.e. the zone with massive investments in leaf chloroplast biogenesis and build-up of photosynthetic capacity. Clustering analyses of gene-specific protein-mRNA ratios revealed co-ordinated shifts in control points in gene expression along the leaf blade developmental gradient. The highest protein-mRNA ratio for each gene generally corresponded to leaf developmental stages in which the protein function was most important, with the exception of the 80S ribosome. Specific examples are discussed in the context of C4 photosynthesis, leaf development and sink-source transitions. This large-scale mRNA-protein correlation analysis in plants (maize) using label-free spectral counting for protein quantification and RNA-seq for mRNA abundance will provide a template for future mRNA-protein correlation studies.

Ontogeny of the maize shoot apical meristem.

The maize (Zea mays) shoot apical meristem (SAM) arises early in embryogenesis and functions during stem cell maintenance and organogenesis to generate all the aboveground organs of the plant. Despite its integral role in maize shoot development, little is known about the molecular mechanisms of SAM initiation. Laser microdissection of apical domains from developing maize embryos and seedlings was combined with RNA sequencing for transcriptomic analyses of SAM ontogeny. Molecular markers of key events during maize embryogenesis are described, and comprehensive transcriptional data from six stages in maize shoot development are generated. Transcriptomic profiling before and after SAM initiation indicates that organogenesis precedes stem cell maintenance in maize; analyses of the first three lateral organs elaborated from maize embryos provides insight into their homology and to the identity of the single maize cotyledon. Compared with the newly initiated SAM, the mature SAM is enriched for transcripts that function in transcriptional regulation, hormonal signaling, and transport. Comparisons of shoot meristems initiating juvenile leaves, adult leaves, and husk leaves illustrate differences in phase-specific (juvenile versus adult) and meristem-specific (SAM versus lateral meristem) transcript accumulation during maize shoot development. This study provides insight into the molecular genetics of SAM initiation and function in maize.

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.

Tissue- and cell-type specific transcriptome profiling of expanding tomato fruit provides insights into metabolic and regulatory specialization and cuticle formation.

Tomato (Solanum lycopersicum) is the primary model for the study of fleshy fruits, and research in this species has elucidated many aspects of fruit physiology, development, and metabolism. However, most of these studies have involved homogenization of the fruit pericarp, with its many constituent cell types. Here, we describe the coupling of pyrosequencing technology with laser capture microdissection to characterize the transcriptomes of the five principal tissues of the pericarp from tomato fruits (outer and inner epidermal layers, collenchyma, parenchyma, and vascular tissues) at their maximal growth phase. A total of 20,976 high-quality expressed unigenes were identified, of which more than half were ubiquitous in their expression, while others were cell type specific or showed distinct expression patterns in specific tissues. The data provide new insights into the spatial distribution of many classes of regulatory and structural genes, including those involved in energy metabolism, source-sink relationships, secondary metabolite production, cell wall biology, and cuticle biogenesis. Finally, patterns of similar gene expression between tissues led to the characterization of a cuticle on the inner surface of the pericarp, demonstrating the utility of this approach as a platform for biological discovery.

Nucleoid-enriched proteomes in developing plastids and chloroplasts from maize leaves: a new conceptual framework for nucleoid functions.

Plastids contain multiple copies of the plastid chromosome, folded together with proteins and RNA into nucleoids. The degree to which components of the plastid gene expression and protein biogenesis machineries are nucleoid associated, and the factors involved in plastid DNA organization, repair, and replication, are poorly understood. To provide a conceptual framework for nucleoid function, we characterized the proteomes of highly enriched nucleoid fractions of proplastids and mature chloroplasts isolated from the maize (Zea mays) leaf base and tip, respectively, using mass spectrometry. Quantitative comparisons with proteomes of unfractionated proplastids and chloroplasts facilitated the determination of nucleoid-enriched proteins. This nucleoid-enriched proteome included proteins involved in DNA replication, organization, and repair as well as transcription, mRNA processing, splicing, and editing. Many proteins of unknown function, including pentatricopeptide repeat (PPR), tetratricopeptide repeat (TPR), DnaJ, and mitochondrial transcription factor (mTERF) domain proteins, were identified. Strikingly, 70S ribosome and ribosome assembly factors were strongly overrepresented in nucleoid fractions, but protein chaperones were not. Our analysis strongly suggests that mRNA processing, splicing, and editing, as well as ribosome assembly, take place in association with the nucleoid, suggesting that these processes occur cotranscriptionally. The plastid developmental state did not dramatically change the nucleoid-enriched proteome but did quantitatively shift the predominating function from RNA metabolism in undeveloped plastids to translation and homeostasis in chloroplasts. This study extends the known maize plastid proteome by hundreds of proteins, including more than 40 PPR and mTERF domain proteins, and provides a resource for targeted studies on plastid gene expression. Details of protein identification and annotation are provided in the Plant Proteome Database.

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.

Large-scale label-free quantitative proteomics of the pea aphid-Buchnera symbiosis.

Many insects are nutritionally dependent on symbiotic microorganisms that have tiny genomes and are housed in specialized host cells called bacteriocytes. The obligate symbiosis between the pea aphid Acyrthosiphon pisum and the γ-proteobacterium Buchnera aphidicola (only 584 predicted proteins) is particularly amenable for molecular analysis because the genomes of both partners have been sequenced. To better define the symbiotic relationship between this aphid and Buchnera, we used large-scale, high accuracy tandem mass spectrometry (nanoLC-LTQ-Orbtrap) to identify aphid and Buchnera proteins in the whole aphid body, purified bacteriocytes, isolated Buchnera cells and the residual bacteriocyte fraction. More than 1900 aphid and 400 Buchnera proteins were identified. All enzymes in amino acid metabolism annotated in the Buchnera genome were detected, reflecting the high (68%) coverage of the proteome and supporting the core function of Buchnera in the aphid symbiosis. Transporters mediating the transport of predicted metabolites were present in the bacteriocyte. Label-free spectral counting combined with hierarchical clustering, allowed to define the quantitative distribution of a subset of these proteins across both symbiotic partners, yielding no evidence for the selective transfer of protein among the partners in either direction. This is the first quantitative proteome analysis of bacteriocyte symbiosis, providing a wealth of information about molecular function of both the host cell and bacterial symbiont.

Strategies for transcriptome analysis in nonmodel plants.

Even with recent reductions in sequencing costs, most plants lack the genomic resources required for successful short-read transcriptome analyses as performed routinely in model species. Several approaches for the analysis of short-read transcriptome data are reviewed for nonmodel species for which the genome of a close relative is used as the reference genome. Two approaches using a data set from Phytophthora-challenged Rubus idaeus (red raspberry) are compared. Over 70000000 86-nt Illumina reads derived from R. idaeus roots were aligned to the Fragaria vesca genome using publicly available informatics tools (Bowtie/TopHat and Cufflinks). Alignment identified 16956 putatively expressed genes. De novo assembly was performed with the same data set and a publicly available transcriptome assembler (Trinity). A BLAST search with a maximum e-value threshold of 1.0 × 10(-3) revealed that over 36000 transcripts had matches to plants and over 500 to Phytophthora. Gene expression estimates from alignment to F. vesca and de novo assembly were compared for raspberry (Pearson's correlation = 0.730). Together, alignment to the genome of a close relative and de novo assembly constitute a powerful method of transcriptome analysis in nonmodel organisms. Alignment to the genome of a close relative provides a framework for differential expression testing if alignments are made to the predefined gene-space of a close relative and de novo assembly provides a more robust method of identifying unique sequences and sequences from other organisms in a system. These methods are considered experimental in nonmodel systems, but can be used to generate resources and specific testable hypotheses.

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.

Sequence comparison of prefrontal cortical brain transcriptome from a tame and an aggressive silver fox (Vulpes vulpes).

BACKGROUND: Two strains of the silver fox (Vulpes vulpes), with markedly different behavioral phenotypes, have been developed by long-term selection for behavior. Foxes from the tame strain exhibit friendly behavior towards humans, paralleling the sociability of canine puppies, whereas foxes from the aggressive strain are defensive and exhibit aggression to humans. To understand the genetic differences underlying these behavioral phenotypes fox-specific genomic resources are needed. RESULTS: cDNA from mRNA from pre-frontal cortex of a tame and an aggressive fox was sequenced using the Roche 454 FLX Titanium platform (> 2.5 million reads & 0.9 Gbase of tame fox sequence; >3.3 million reads & 1.2 Gbase of aggressive fox sequence). Over 80% of the fox reads were assembled into contigs. Mapping fox reads against the fox transcriptome assembly and the dog genome identified over 30,000 high confidence fox-specific SNPs. Fox transcripts for approximately 14,000 genes were identified using SwissProt and the dog RefSeq databases. An at least 2-fold expression difference between the two samples (p < 0.05) was observed for 335 genes, fewer than 3% of the total number of genes identified in the fox transcriptome. CONCLUSIONS: Transcriptome sequencing significantly expanded genomic resources available for the fox, a species without a sequenced genome. In a very cost efficient manner this yielded a large number of fox-specific SNP markers for genetic studies and provided significant insights into the gene expression profile of the fox pre-frontal cortex; expression differences between the two fox samples; and a catalogue of potentially important gene-specific sequence variants. This result demonstrates the utility of this approach for developing genomic resources in species with limited genomic information.