Proteome Analysis of the Inner Integument from Developing

Anatomical structure of J. curcas seed. ... (27) Six biological replicates were performed for each sample (Figure 2). Each biological replicate was co...
0 downloads 0 Views 1MB Size
Subscriber access provided by UZH Hauptbibliothek / Zentralbibliothek Zuerich

Article

Proteome analysis of the inner integument from developing seeds of Jatropha curcas L. Emanoella L. Soares, Mohibullah Shah, Arlete A. Soares, José Hélio Costa, Paulo C. Carvalho, Gilberto B. Domont, Fábio C.S. Nogueira, and Francisco A.P. Campos J. Proteome Res., Just Accepted Manuscript • DOI: 10.1021/pr5004505 • Publication Date (Web): 10 Jul 2014 Downloaded from http://pubs.acs.org on July 14, 2014

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Journal of Proteome Research is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

Proteome analysis of the inner integument from developing Jatropha curcas L. seeds

Emanoella L. Soares1#, Mohibullah Shah1#, Arlete A. Soares2, José H. Costa1, Paulo Carvalho4, Gilberto B. Domont3, Fábio C.S. Nogueira3* and Francisco A.P. Campos3* 1 Department of Biochemistry and Molecular Biology, Federal University of Ceara, Fortaleza, Brazil 2 Department of Biology, Federal University of Ceara, Fortaleza, Brazil 3 Proteomic Unit, Institute of Chemistry, Federal University of Rio de Janeiro, Rio de Janeiro, Brazil 4 Laboratory for Proteomics and Protein Engineering, Carlos Chagas Institute, Fiocruz – Parana, Brazil #Contributed equally *Correspondence: Prof. Francisco A.P. Campos, Department of Biochemistry and Molecular Biology, Federal University of Ceara, 60455-900 Fortaleza, Brazil Email: [email protected] Fax: +55-85-33669829 Prof. Fábio C. S. Nogueira, Proteomic Unit, Institute of Chemistry, Federal University of Rio de Janeiro, 21941-909 Rio de Janeiro, Brazil 1

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 33

Email: [email protected] Phone: +55-21-39388862 Fax: +55-21-39387266

2

ACS Paragon Plus Environment

Page 3 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

Keywords: label-free quantitative proteomics, programmed cell death, Jatropha curcas, oilseed, seed coat Abstract In this study, we performed a systematic proteomic analysis of the inner integument from developing seeds of Jatropha curcas and further explored the protein machinery responsible for generating the carbon and nitrogen sources to feed the growing embryo and endosperm. The inner integument of developing seeds was dissected into two sections called distal and proximal and proteins were extracted from these sections and from the whole integument and analyzed using an EASY-nanoLC system coupled to an ESI-LTQ-Orbitrap Velos mass spectrometer. We identified 1526, 1192, and 1062 proteins from the proximal, distal and whole inner integuments, respectively. The identifications include those of peptidases and other hydrolytic enzymes that play a key role in developmental programmed cell death (PCD), and proteins associated with the cell wall architecture and modification. As many of these proteins are differentially expressed within the integument cell layers, these findings suggest that the cells mobilize an array of hydrolases to produce carbon and nitrogen sources from proteins, carbohydrates, and lipids available within the cells. Not least, the identification of several classes of seed storage proteins in the inner integument provides additional evidence of the role of the seed coat as a transient source of reserves for the growing embryo and endosperm.

3

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 33

1. Introduction The seeds of J. curcas are considered to be a potential source of raw material for the production of biodiesel1 and two by-products of the oil extraction viz. the protein-rich press-cake and the seed coat, also offer opportunities for adding value to the crop. The former can potentially be used as animal feed and the latter, representing almost 40% of the seed weight, constitutes a potential source of biomass for the production of second-generation ethanol.2 The full exploitation of such potentials is hindered by a lack of knowledge of the chemical identity and biosynthesis of cell wall components of the seed coat and by the presence, in the seeds, of high concentrations of phorbol esters, tetracyclic diterpenoids known for their tumor-promoting activity, which represent serious risks for humans, animals and for the environment.3 In mature seeds, the highest concentrations of phorbol esters are found in the inner layer of the seed coat, called the tegmen, from where they diffuse to the endosperm.4,5 The seeds of Angiosperms are derived from the ovule, an organ that harbors the female gametophyte which in turn is surrounded by the nucellus and one or two integuments.6 Following fertilization of the egg cell and the central cell, the female gametophyte will differentiate into the diploid embryo and triploid endosperm respectively and the maternally derived diploid seed coat develops from the integuments of the ovule.7 The carbon and nitrogen sources necessary for the developing embryo and endosperm are provided by the long distance transport of metabolites, such as sucrose, amino acids and amides through the maternal vascular tissue. 8 While knowledge of the developmental biology of embryo and endosperm is expanding rapidly9,10 the cellular events underlying seed coat differentiation have received comparatively little attention.11,12 However, it is well established that during early to mid stages of development, the seed coat acts as a maternal conduit to convey nutrients for the developing embryo and endosperm13-15 and at maturity it provides protection, regulates germination and additionally promotes seed dispersal.16 4

ACS Paragon Plus Environment

Page 5 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

Details of the protein machinery responsible for producing the nitrogen and carbon sources for feeding the embryo and endosperm during seed development are elusive. Although differentiation of the integuments to seed coat is known to involve developmentally regulated programmed cell death (PCD) triggered by vacuole collapse,17,18 very little is known about the synthesis, subcellular location, and site of action of the peptidases, lipases, carbohydrases, and nucleases. These enzymes are needed to digest cell components, including the cell wall, so that the products may be used as energy sources and/or buildings blocks for biosynthetic reactions within the embryo and endosperm. Most of the studies related to the seed coat have been restricted to individual enzymes19-21 or to high throughput transcriptomic analysis.10 Recently Miernyk and Johnston11 published the very first in-depth proteomic analysis of the seed coat, providing a glimpse of the dynamic changes in the proteome of the testa of developing soybean seeds, but similar studies for the seeds of Euphorbiaceae are lacking. In this present work, we capitalized on a detailed description of the histological development of the seed coat of Jatropha curcas22 to guide our proteomic study within the layers of the inner integument during their midstage development phase. This allowed us to report with unprecedented details of the proteome of the enzymatic machinery responsible for executing the developmentally controlled PCD within this tissue.

2. Material and Methods 2.1. Plant materials and histological analysis Seeds (25 days after pollination) were collected from hand-pollinated three year-old plants that were field-grown and irrigated three days per week until the field capacity, according to the measurement of evapotranspiration by a Class A land tank. Histological analyses were performed as described by Rocha et al. 23 Briefly, the inner integument was isolated with a sharp 5

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 33

spatula, dissected into two regions having the vascular bundle as reference – one of these internal to the vascular bundle faced the central cavity (proximal region) and the other external to the vascular bundle, facing the exotegmen (distal region) (Figure 1). The separation of these regions was facilitated by the large size of the inner integument and the easy visualization of the vascular bundle with the aid of a binocular microscope. The middle section between the proximal and distal regions, comprising the vascular bundle, was discarded in order to avoid crosscontamination. In addition to the proximal and distal regions we also collected the whole intact inner integument (here called total integument). For histological analysis tissues were fixed in Karnovsky solution,

24

dehydrated in an ethanol series of 10, 20, 30, 40, 50, 60, 70, 80, 90 and

100% (one hour per step) and slowly embedded in historesin LEICA for 20 days. Sections of 5 µm were prepared in a rotatory microtome LEICA 2065, stained with toluidine blue 0.5% in borax 0.12% followed by basic fuchsin 25 (0.5%) and mounted with Tissue Mount.

2.2. Sample preparation and protein determination After separation, the three tissues (whole integument, proximal and distal sections) were immediately frozen in liquid nitrogen and freeze-dried. Subsequently, samples were ground down to powder and 100 mg of each fraction were separately subjected to protein extraction according to Vasconcelos, et al.26 Briefly, powdered tissues were homogenized in pyridine buffer (50 mM pyridine, 10 mM thiourea and 1% SDS, pH 5.0) with polyvinyl−polypyrrolidone (PVPP) in a proportion of 1:40:2 (w/v/w), respectively. The homogenate was stirred for 4 hours at 4º C and centrifuged for 30 minutes at 10000 g. Proteins in the supernatant were precipitated overnight with cold 10% trichloroacetic acid (TCA) in acetone and centrifuged for 15 minutes at 10,000 g. Pellet was washed three times with cold acetone, dried under vacuum and dissolved in 7 M urea/2 M thiourea. Protein concentrations were determined by the Bradford assay.27 Six

6

ACS Paragon Plus Environment

Page 7 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

biological replicates were performed for each sample (Figure 2). Each biological replicate was constituted from tissues of approximately 60 seeds collected from approximately 20 fruits.

2.3. LC-MS/MS and data analysis Trypsin digestions were performed using 100 µg of the proteins from each of the three samples as described by Pinheiro et al.

28

Peptides were dissolved in 80 µL of 0.5% formic acid

and diluted five times with the same solution. 4 µL of the diluted samples were applied to an EASY II nano LC system (Proxeon Biosystem) coupled online to an ESI-LTQ Orbitrap Velos mass spectrometer (Thermo Fisher Scientific). Peptides were loaded in a trap column (150 µm × 2 cm) packed in-house with C-18 ReproSil 3 µm resin (Dr. Maisch) and eluted in an analytical column (100 µm x 15 cm) packed with the same material. Peptide separations were performed using a gradient from 100% of A (0.1% formic acid) to 35% of B (0.1% formic acid, 95% acetonitrile) for 150 min, followed by 35% to 90% of B for 15 min, and 90% for 5 min. MS1 spectra were acquired in a positive mode using the data-dependent automatic (DDA) survey MS scan. Each DDA consisted of a survey scan in the m/z range 300−2000 and resolution 60,000 with a target value of 1 × 10-6 ions. The ten most intense ions were subjected to MS2 acquisition in the LTQ using a normalized collision-induced dissociation (CID) of previously fragmented ions were dynamically excluded for 60 s. Raw data were viewed in Xcalibur v.2.1 (Thermo Scientific). Three raw files were generated for each biological replicate of the distal and proximal regions while two raw files were generated for each biological replicate of the total integument, representing three and two technical replicates, respectively (Figure 2). Raw files were converted into MS2 files using PatternLab’s RawReader. ProLuCID v1.329 was used to perform peptide spectrum matching against the J. curcas proteins database30 downloaded from http://www.kazusa.or.jp/jatropha/ 7

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 33

September 2012, combined with J. curcas chloroplast genome encoded proteins, downloaded from NCBI September 2012. Search parameters were: semi tryptic hydrolysis, two missed cleavages, oxidation of methionine as variable and carbamidomethylation as fixed modifications, and a peptide tolerance of 50 ppm. Search results were subsequently filtered/processed through Search Engine Processor31 to achieve a list of identified proteins with 1% false discovery rate (FDR) at protein level. The validity of the Peptide Sequence Matches (PSM) was assessed using the Search Engine Processor v. 2.2.0.2 (SEPro). Identifications were grouped by the charge state (+2 and > +3) and then by tryptic status (tryptic and semi-tryptic), resulting in four distinct subgroups. For each result, the ProLuCID XCorr, DeltaCN and ZScore values were used to generate a Bayesian discriminating function. A cutoff score was established to accept a false-discovery rate (FDR) of 1% based on the number of decoys. This procedure was independently performed on each data subset, resulting in a false-positive rate that was independent of tryptic status or charge state. A minimum sequence length of 6 amino acids was required. Results were post-processed to accept PSM with less than 8 ppm. Proteins were grouped according to the maximum parsimony by the bipartite graph approach.32 Proteins appearing in at least two biological replicates, considered as representatives of that biological sample, were blasted against NCBI non redundant database using Blast2GO annotation tool33 and used for downstream analysis. Spectral counts of the identified proteins were obtained using the PatternLab computational environment.34 Gene ontology (GO) annotation was done using AgBase tools and database (http://agbase.msstate.edu/index.html). Goanna35 was used to retrieve the GO annotations assigned on the basis of sequence similarities and Plant GOSlim was introduced to summarize the sub-categories of the identified proteins.

8

ACS Paragon Plus Environment

Page 9 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

Differentially expressed proteins between the distal and proximal regions were pinpointed by using PatternLab’s revised TFold module 36, that uses the spectral count method. Spectral counting is a label-free semi-quantitative measure of the protein abundance and defined as, the total number of spectra identified for a protein.

37

Briefly for a given FDR bound, the

revised TFold uses a theoretical FDR estimator (in this case, the Benjamini-Hochberg) to maximize the number of identifications that satisfy both a fold-change cutoff that varies with the t-test P-value as a power law and a stringency criterion that aims to detect lowly abundant proteins. Fold changes were calculated by dividing the average spectral counts obtained in the distal region by those obtained in the proximal region.

2.4 RNA isolation, cDNA synthesis and reverse transcription quantitative real-time PCR (RT-qPCR) Total RNA from proximal and distal regions was isolated from 300mg of tissue, using the Plant RNA Purification Reagent (Invitrogen®) according to the manufacturer’s instructions. RNA was digested with DNase I (Invitrogen) for removal of genomic DNA. The first-strand cDNA synthesis was performed with 0.5 µg of DNase I treated-RNA using the ImProm-II™ Reverse Transcription System (Promega) according to the manufacturer’s instructions and was stored at -20 °C until use. RT-qPCR was carried out using the Fast SYBR Master Mix (Applied Biosystems) with a Mastercycler® ep realplex (Eppendorf AG, Hamburg) under the following conditions: 10 min at 95◦C for activation of the enzyme, followed by 40 cycles of denaturation at 95◦C for 15 s, annealing at 60◦C for 20 s, and extension at 60◦C for 20 s. A melting curve additional step was carried out to check the specificity of the primers. Each reaction was performed with 20 ng of cDNA and 400 nM of each primer in a final volume of 20 µL. The sequences coding for the target genes (arabinofuranosidase - Jcr4S00428.20; aspartic peptidase - Jcr4S06063.10; cysteine 9

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 33

peptidase - Jcr4S01104.40; cysteine peptidase inhibitor - Jcr4S00024.130; serine peptidase Jcr4S01605.60; serpin - Jcr4S00079.140 and subtilisin - Jcr4S01651.50, Jcr4S01323.10) were retrieved from the genome of J. curcas (www.kazusa.or.jp/jatropha/) to design specific primer sets (Supplementary Table I) using the Perl Primer v1.1.19 software.38 Three reference genes (glyceraldehyde-3-phosphatedehydrogenase, polyubiquitin-3 and ubiquitin conjugating enzyme) were used for data normalization.23 Three biological replicates were performed for each sample, each one with two technical replicates in quadruplicate (n=8). The analysis of relative expression was carried out by the 2–∆∆CT method.39

3. Results and Discussion 3.1 Histology of seed development The mature ovule of J. curcas is anatropous, crassinucellar, bitegmic, and the vascular bundles present in the outer integument extend by postchalazal branching from the chalaza through the inner integument.22,40 Both integuments are composed of three layers: endotegmen, mesotegmen and exotegmen in the inner integument and endotesta, mesotesta and exotesta in the outer integument. During development the integuments will differentiate to form the seed coat. At the mid stage of seed development, 25 days after pollination (DAP), different parts of the inner integument are experiencing different changes in its cell features indicative that PCD is occurring in a gradual way (Figure 1). These observations led us to select this stage for proteomic analysis to see the differences in the protein profiles in different regions of the inner integument. At 25 DAP the embryo is at the heart stage, the cellularization of the endosperm is well underway (Figure 1A) and the cell layers closest to the exotegmen are small and become progressively larger and vacuolated as they are positioned closer to the endosperm and those in the immediate vicinity of the endosperm collapse, forming a layer of cell debris (Figure 1B). In 10

ACS Paragon Plus Environment

Page 11 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

later stages, when cellularization of the nuclear endosperm is complete and the embryo is fully developed, the mesotegmen cells, which are internal to the vascular bundle, are completely consumed, the cells external to the vascular bundle are crushed and those from the exotegmen become sclerified forming the main mechanical barrier of the mature seed coat. At this stage the three layers of the outer integument are differentiated. Figure 1B indicates that cells closer to the exotegmen and those closer to the endosperm are at different metabolic states and are experiencing PCD at different degrees, suggesting the possible existence within the inner integument of a concentration gradient of the proteins involved with the execution of PCD (e.g. peptidases, lipases and carbohydrases). To further investigate this, we took advantage of the large size of the developing seeds to isolate the whole inner integument and to obtain from it proximal and distal sections to the developing endosperm (Figure 1A) and to perform an indepth proteomic analysis of each of the three parts.

3.2 Protein identification and functional classification Our results permitted the identification of 1752, 1405, and 1251 proteins from distal, proximal and total integument, representing 1169, 912, and 823 protein groups, respectively, at a 1% FDR. Supplementary Table II summarizes details about the number of MS, MS2, PSMs, peptides and unique peptides for all tissues and biological/technical replicates. The proteins identified in at least two biological replicates were considered as highly confident; these identifications comprise 1526 (distal), 1192 (proximal), and 1062 (total integument) proteins summing up a total of 1770 (Supplementary Table III). The identification of sucrose synthases and cell wall and apoplastic invertases, which are key players in the control of nutrient supply to the embryo and endosperm,41 indicated that the inner integument of developing J. curcas seeds is metabolically very active. This is further supported by the identification of proteins belonging to major biochemical pathways, such as those related to amino acid and protein biosynthesis, and a 11

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 33

number of transporters (Supplementary Table III). Additionally, the role of the inner integument as a transient storage source of nutrients to the developing embryo and endosperm is supported by the identification of several classes of seed storage proteins such as nutrient reservoir, globulins 11S, legumins, and glutelins (Table 1). Although a number of proteins involved in secondary metabolism were identified, we were unable to identify any casbene synthase or other protein that may be involved in the biosynthesis of phorbol esters, the major toxic components of J. curcas seeds. This result is in agreement with recent findings by Nakano et al.42 which showed that in J. curcas, casbene synthase is expressed in seedlings, mature leaves and in the flesh of developing fruits, but not in developing seeds. Another recent proteome analysis of plastids from developing J. curcas seeds28 reported no evidence for the presence of casbene synthase among more than 1100 proteins. These results indicate that the phorbol esters that are found in the seeds are probably translocated either from the leaves or from the roots, but experimental support for this hypothesis is lacking. Despite the heterogeneity in the total number of the proteins identified in each of the three samples analyzed (Supplementary Table III), we did not observe significant differences in the distribution of these proteins into different Gene Ontology categories (Supplementary Table IV). However, Gene Ontology annotation of the proteins unique to distal and proximal region of the integument revealed a clear difference in their distribution to some of the subcategories of GO Molecular Function and Cellular Component (Supplementary Table V and Figure 3). Proteins related to hydrolase activity of the GO Molecular Function (Figure 3A) were found to comprise the largest functional class in the proximal region as compared to the distal region where proteins related to protein binding (Figure 3A) make up the largest functional class. This indicates that most of the metabolic activity located in the proximal region is geared towards the reallocation to the embryo and endosperm of the catabolic products produced by the digestion of proteins, carbohydrates, lipids, and nucleic acids within the cells undergoing PCD. Moreover, 12

ACS Paragon Plus Environment

Page 13 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

the percentage of the cell wall sub-category from the GO Cellular Component (Figure 3B) was greater in the proximal region than in the distal region. The proteins mapping within this category are mostly related to the processes of cell wall degradation which further supports a metabolism geared towards the digestion of carbohydrate polymers and proteins that comprise the cell wall.

3.3 Differentially expressed proteins A PCA plot was performed using PatternLab´s Buzios module43; the plot confirms the proximal and distal proteomic profiles to be very distinct (Supplementary Figure I). We used PatternLab’s Approximately Area Proportional Venn Diagram module to pinpoint proteins uniquely identified in distal and proximal regions of the inner integument. The analysis only considered proteins found in two or more biological replicates from a particular sample. Differentially expressed proteins found in two or more biological replicates identified in the distal and proximal regions of the inner integument were discriminated by PatternLab’s TFold module using a q-value of 0.05 (Figure 4). Only proteins that satisfy all statistical tests (blue dots in Figure 4) were considered, and those dots found in the upper section of the plot are upregulated in the distal sample; likewise, those found in the lower section are up-regulated in the proximal sample. We recall that the TFold module uses a theoretical FDR estimator to maximize identifications satisfying both a fold-change cutoff that varies with the t-test p-value as a power law and a stringency criterion that aims to fish out proteins of low abundance that are likely to have had their quantitation compromised. Differentially expressed proteins between distal and proximal regions are present in Supplementary Table VI.

13

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 33

3.4 Functional classes related to programmed cell death 3.4.1 Peptidases and peptidase inhibitors Of all of the identified proteins, almost 10% belong to one of the four major mechanistic classes of peptidases and their inhibitors or are catalytic subunits of the proteasome or ubiquitination proteins (Supplementary Table VII). The identification of proteasome components and ubiquitination proteins indicates a proper functioning of the proteasome, a conditio sine qua non for the occurrence of PCD in plant cells.44 Many of the identified peptidases are known to have a role in PCD. For example, ϒ-vacuolar processing enzyme (ϒ-VPE), only identified in the distal region, is a homologue to a VPE also identified in castor seeds during the proteomic analysis of nucellus tissue undergoing PCD.45 The VPEs are a sub-class of cysteine peptidases (CP) which possess caspase-1-like activity and are used by the vacuolar-collapse system to promote developmentally controlled PCD in plant tissues.17 The KDEL-tailed CP is another type of CP that has been found to be a hallmark of PCD in plant tissues undergoing PCD.23,45-47 Here we identified a KDEL-tailed CP (Jcr4S01104.40) having a statistically higher expression in the proximal region (Supplementary Table VI); this evidence supports the hypothesis that these proteinases are released into the cytoplasm following the vacuole collapse triggered by the action of the VPEs.21 The spatial and temporal differences in expression of the VPEs and KDEL-tailed CP in the inner integument of J. curcas seeds underline the different roles of these CP in PCD and are in line with recent findings by Rocha et al.23 which indicated that during seed development, transcripts for ϒ-VPE are detected earlier than those for the KDEL-tailed CP. Our results disclose identifications of several serine peptidases (SP) and serine peptidase inhibitors (SPI) (Supplementary Table VII) of which 10 are differentially expressed (Supplementary Table VI). Notably, three subtilisin-like serine peptidases presented spectralcounting fold changes greater than 10 fold higher in the proximal region. These subtilisin-like 14

ACS Paragon Plus Environment

Page 15 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

serine peptidases are characterized by an aspartate, histidine and serine catalytic triad, have been shown to display caspase like activity, and were seen to be associated with cell death in several plant species.48 In order to further confirm these results, we used qPCR to measure the relative expression of transcripts coding for selected peptidase and peptidase inhibitors in distal and proximal sections of the inner integument. As shown in Table 2, the relative expression of the transcripts analyzed, matches closely the relative abundance of the proteins coding for these transcripts as determined in our proteome analysis. Identification of this remarkable number and diversity of peptidases belonging to the four mechanistic classes of peptidases in tissues undergoing PCD underlines the importance of proteolysis for providing amino acids for protein synthesis in the filial tissues as suggested by Gallardo et al.49 Moreover, such diversity highlights the role of seed integuments in the developmental biology of J. curcas seeds. Miernyk and Johnston11 and Gallardo et al.49 have identified an abundance of peptidase in several stages of the development of the seed coat of soybean and Medicago truncatula, but the possibility that the abundance and variety of these peptidases in the developing seed coat may be related to PCD was not discussed.

3.4.2 Proteins related to cell wall metabolism Besides the peptidases referred to above, we identified proteins involved in the architecture of the cell wall such as several expansins and extensins and various other proteins associated with cell wall modification and degradation processes, such as arabinofuranosidases,

lysosomal

alpha-mannosidases,

pectinesterase,

alpha-

alpha-Land

beta-

glucosidases, beta-xylosidases, alpha-galactosidases, glucanases, polygalacturonases, among others (Supplementary Table III). Both beta-glucosidase and xylosidase are known to be involved in the breakdown or modification of cell wall hemicelluloses and pectins.50,51 Alpha-dgalactosidase is associated with the mobilization of galactomannan, while beta-galactosidase is 15

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 33

involved in the degradation of the galactan side chain.51 Pectinesterases also called pectin methylesterases (PMEs) catalyse the demethyl-esterification of the pectin domains and make it feasible to be degraded by the action of polygalactorunases.52 Additionally, we identified pectinesterase inhibitors (Supplementary Table III) which modulate pectinesterase activity.50 Similarly the enzyme alpha-arabinofuranosidase identified here, in combined action with betagalactosidase and beta-glucuronidase (not identified here) is known to degrade the carbohydrate moieties of arabinogalactan proteins (AGPs).51 Our results show that the expression of many of these proteins is differentially regulated, as is the case of alpha-L-arabinofuranosidase, pectinesterase inhibitor and polygalacturonase that are more expressed in the proximal section than in the distal section of the inner integument (Supplementary Table VI). As measured by qPCR, the relative levels of transcripts coding for the alpha-L-arabinofuranosidase in proximal and distal sections of the inner integument (Table 2), matches closely the relative abundance of the protein coded by this transcript as determined in our proteome analysis (Supplementary Table VI). Additionally, several other proteins such as glucosidases, xylanases, mannosidases and pectinesterase, among others, were identified only in the proximal region of the inner integument (Supplementary Table III). These results indicate that along other possible roles, these polysaccharide acting enzymes together with different kind of peptidases, are acting to liberate all possible carbon and nitrogen sources available in the cell, including those that belong to the cell wall.

3.5.3. Proteins related to lipid catabolism It is noteworthy to point out that our results disclose several proteins involved in the degradation of lipids, viz.: phospholipase D, lipoxygenases, acyl-CoA oxidases, 3-ketoacyl thiolase, among others (Supplementary Table III). Phospholipase is an enzyme that catalyses the 16

ACS Paragon Plus Environment

Page 17 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

hydrolysis of structural phospholipids into hydroperoxides leading to the loss of membrane integrity, an event associated with PCD.53 This enzyme is known to be involved in the largescale degradation of lipids associated with cell death54 and sometimes contributes to caspasedependent cell death signaling.55 Fatty acids resulting from the degradation of the cell membrane are oxidized by the action of lipoxygenases56 and acyl hydrolases. The appearance of these enzymes along with other hydrolases, reinforce the occurrence of PCD in this maternal tissue.

3.5.3 Comparison with other proteomic analyses of J. curcas So far only a few proteomic studies have been published for J. curcas, concentrating mainly on endosperm, embryo, oilbodies and seed kernel, that collectively produced limited information57-63. We have recently published a first in depth proteomic analysis of J. curcas by analyzing the proteome of plastids isolated from the endosperm of developing J. curcas seeds, which resulted into the identification of 923 proteins28. Comparing the integument and plastids proteomes, we have identified 2153 proteins among which 1230 are unique to integument while 383 are unique to plastids. Proteins unique to the integument proteome consist mainly of different kinds of hydrolases like peptidases and glycosidases, which are involved in the degradation of complex molecules to liberate nutrients to the developing endosperm and embryo. On the other hand, proteins unique to the plastids proteome consist of proteins involved in the lipids metabolism. The other proteomic analyses of J. curcas present in literature57-63 collectively identified 209 non redundant proteins. Comparison of our plastids and integument proteomes to these studies showed, that out 209, only 42 proteins were not identified in our results, which is probably due to the use of different tissues and experimental conditions in each of these different studies.

17

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 33

4. Conclusion The data presented here demonstrate the role of the inner integument of developing seeds of J. curcas in providing carbon and nitrogen sources for the growing embryo and endosperm. The identification of several classes of peptidases, particularly of ϒ-VPE and KDEL-tailed CP highlights the role of developmental PCD in the developmental biology of seeds of J. curcas. Additionally, the demonstration of the differential expression of these peptidases and other hydrolases within cell layers of the inner integument after the triggering of PCD, suggests that the cells mobilize a whole array of hydrolases to produce carbon and nitrogen sources from proteins, carbohydrates and lipids available within the cells. Not least, the identification of several classes of seed storage proteins in the inner integument is an additional evidence of the role of the seed coat as a transient source of reserves for the growing embryo and endosperm.

Acknowledgements

This work was supported by grants from PETROBRAS, Banco do Nordeste, FUNCAP, CNPq, CAPES and TWAS. Data availability Mass spectrometry raw data and Search Engine Processor (SEPro) file are available at: Chorus repository website (chorusproject.org; Project Name: Proteome analysis of the inner integument from developing Jatropha curcas L seeds) and http://yoda.iq.ufrj.br/mohib2014/.

18

ACS Paragon Plus Environment

Page 19 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

5. References (1) Contran, N.; Chessa, L.; Lubino, M.; Bellavite, D.; Roggero, P. P.; Enne, G. State-of-the-art of the Jatropha curcas productive chain: from sowing to biodiesel and by-products. Ind. Crop Prod. 2013, 42, 202-215. (2) Wever, D.-A. Z.; Heeres, H. J.; Broekhuis, A. A. Characterization of Physic nut (Jatropha curcas L.) shells. Biomass Bioenergy 2012, 37, 177-187. (3) Gressel, J. Transgenics are imperative for biofuel crops. Plant Sci. 2008, 174 (3), 246-263. (4) King, A. J.; Montes, L. R.; Clarke, J. G.; Affleck, J.; Li, Y.; Witsenboer, H.; van der Vossen, E.; van der Linde, P.; Tripathi, Y.; Tavares, E.; Shukla, P.; Rajasekaran, T.; van Loo, E. N.; Graham, I. A. Linkage mapping in the oilseed crop Jatropha curcas L. reveals a locus controlling the biosynthesis of phorbol esters which cause seed toxicity. Plant Biotechnol. J. 2013, 11 (8), 986-96. (5) He, W.; King, A. J.; Khan, M. A.; Cuevas, J. A.; Ramiaramanana, D.; Graham, I. A. Analysis of seed phorbol-ester and curcin content together with genetic diversity in multiple provenances of Jatropha curcas L. from Madagascar and Mexico. Plant Physiol. Biochem. 2011, 49 (10), 1183-1190. (6) Boesewinkel, F. D.; Bouman, F. The seed: structure. In Embryology of Angiosperms, Johri, B., Ed. Springer Berlin Heidelberg: 1984; pp 567-610. (7) Ingram, G. C. Family life at close quarters: communication and constraint in angiosperm seed development. Protoplasma 2010, 247 (3-4), 195-214. (8) Patrick, J. W.; Offler, C. E. Compartmentation of transport and transfer events in developing seeds. Journal of Experimental Botany 2001, 52, 551-564. (9) Harada, J. J.; Pelletier, J. Genome-wide analyses of gene activity during seed development. Seed Sci. Res. 2012, 22, (Supplement S1), S15-S22. (10) Belmonte, M. F.; Kirkbride, R. C.; Stone, S. L.; Pelletier, J. M.; Bui, A. Q.; Yeung, E. C.; Hashimoto, M.; Fei, J.; Harada, C. M.; Munoz, M. D.; Le, B. H.; Drews, G. N.; Brady, S. M.; Goldberg, R. B.; Harada, J. J. Comprehensive developmental profiles of gene activity in regions and subregions of the Arabidopsis seed. Proc. Natl. Acad. Sci. U S A 2013, 110 (5), 435-44. (11) Miernyk, J. A.; Johnston, M. L. Proteomic analysis of the testa from developing soybean seeds. J. Proteomics 2013, 89, 265-72. (12) Verdier, J.; Dessaint, F.; Schneider, C.; Abirached-Darmency, M. A combined histology and transcriptome analysis unravels novel questions on Medicago truncatula seed coat. J. Exp. Bot. 2013, 64 (2), 459-70. (13) Heim, U.; Weber, H.; Baumlein, H.; Wobus, U. A sucrose-synthase gene of Vicia faba L.: expression pattern in developing seeds in relation to starch synthesis and metabolic regulation. Planta 1993, 191 (3), 394-401. 19

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 33

(14) Weber, H.; Borisjuk, L.; Heim, U.; Buchner, P.; Wobus, U. Seed coat-associated invertases of fava bean control both unloading and storage functions: cloning of cDNAs and cell typespecific expression. Plant Cell 1995, 7 (11), 1835-1846. (15) Rochat, C.; Boutin, J.-P. Metabolism of phloem-borne amino acids in maternal tissues 2 Pisum sativum L. J. Exp. Bot. 1991, 42 (2), 207-214. (16) Dean, G.; Cao, Y.; Xiang, D.; Provart, N. J.; Ramsay, L.; Ahad, A.; White, R.; Selvaraj, G.; Datla, R.; Haughn, G. Analysis of gene expression patterns during seed coat development in Arabidopsis. Mol. Plant 2011, 4 (6), 1074-91. (17) Hara-Nishimura, I.; Hatsugai, N. The role of vacuole in plant cell death. Cell Death Differ. 2011, 18 (8), 1298-304. (18) Haughn, G.; Chaudhury, A. Genetic analysis of seed coat development in Arabidopsis. Trends Plant Sci. 2005, 10 (10), 472-7. (19) Nakaune, S.; Yamada, K.; Kondo, M.; Kato, T.; Tabata, S.; Nishimura, M.; HaraNishimura, I. A Vacuolar Processing Enzyme, δVPE, is involved in seed coat formation at the early stage of seed development. Plant Cell 2005, 17 (3), 876-887. (20) Wan, L.; Xia, Q.; Qiu, X.; Selvaraj, G. Early stages of seed development in Brassica napus: a seed coat-specific cysteine proteinase associated with programmed cell death of the inner integument. Plant J. 2002, 30 (1), 1-10. (21) Trobacher, C. P.; Senatore, A.; Holley, C.; Greenwood, J. S. Induction of a ricinosomalprotease and programmed cell death in tomato endosperm by gibberellic acid. Planta 2013, 237 (3), 665-79. (22) Singh, R. P. Structure and development of seeds in Euphorbiaceae: Jatropha species. Beitr Biol. Pflanzen 1970, 47, 79-90. (23) Rocha, A. J.; Soares, E. L.; Costa, J. H.; Costa, W. L.; Soares, A. A.; Nogueira, F. C.; Domont, G. B.; Campos, F. A. Differential expression of cysteine peptidase genes in the inner integument and endosperm of developing seeds of Jatropha curcas L. (Euphorbiaceae). Plant Sci. 2013, 213, 30-7. (24) Karnovsky, M.J. A formaldehyde-glutaraldehyde fixative of high osmolality for use in electron microscopy. J. Cell Biol. 1965, 27, 137–138. (25) Junqueira, C.U. O uso de cortes finos de tecidos na Medicina e Biologia. Meios e Métodos. 1990, 167–171. (26) Vasconcelos, É. A. R.; Nogueira, F. C. S.; Abreu, E. F. M.; Gonçalves, E. F.; Souza, P. A. S.; Campos, F. A. P. Protein extraction from cowpea tissues for 2-D gel electrophoresis and MS analysis. Chromatographia 2005, 62 (7-8), 447-450. (27) Bradford, M. M. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem. 1976, 72, 248-54. 20

ACS Paragon Plus Environment

Page 21 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

(28) Pinheiro, C. B.; Shah, M.; Soares, E. L.; Nogueira, F. C.; Carvalho, P. C.; Junqueira, M.; Araujo, G. D.; Soares, A. A.; Domont, G. B.; Campos, F. A. Proteome analysis of plastids from developing seeds of Jatropha curcas L. J. Proteome Res. 2013, 12 (11), 5137-45. (29) Xu, T.; Venable, J. D.; Park, S. K.; Cociorva, D.; Lu, B.; Liao, L.; Wohlschlegel, J.; Hewel, J.; Yates, J. R. I. ProLuCID, a fast and sensitive tandem mass spectra-based protein identification program. Mol Cell Proteomics 2006, 5, S174. (30) Hirakawa, H.; Tsuchimoto, S.; Sakai, H.; Nakayama, S.; Fujishiro, T.; Kishida, Y.; Kohara, M.; Watanabe, A.; Yamada, M.; Aizu, T.; Toyoda, A.; Fujiyama, A.; Tabata, S.; Fukui, K.; Sato, S. Upgraded genomic information of Jatropha curcas L. Plant Biotechnol. J. 2012, 29 (2), 123130. (31) Carvalho, P. C.; Fischer, J. S.; Xu, T.; Cociorva, D.; Balbuena, T. S.; Valente, R. H.; Perales, J.; Yates, J. R., 3rd; Barbosa, V. C. Search engine processor: filtering and organizing peptide spectrum matches. Proteomics 2012, 12 (7), 944-9. (32) Zhang, B.; Chambers, M. C.; Tabb, D. L., Proteomic Parsimony through Bipartite Graph Analysis Improves Accuracy and Transparency. Journal of Proteome Research 2007, 6 (9), 3549-3557. (33) Conesa, A.; Gotz, S. Blast2GO: A comprehensive suite for functional analysis in plant genomics. Int J. Plant Genomics 2008, 2008, 619832. (34) Carvalho, P. C.; Fischer, J. S. G.; Xu, T.; Yates, J. R.; Barbosa, V. C., PatternLab: From Mass Spectra to Label-Free Differential Shotgun Proteomics. In Current Protocols in Bioinformatics, John Wiley & Sons, Inc.: 2012. (35) McCarthy, F. M.; Wang, N.; Magee, G. B.; Nanduri, B.; Lawrence, M. L.; Camon, E. B.; Barrell, D. G.; Hill, D. P.; Dolan, M. E.; Williams, W. P.; Luthe, D. S.; Bridges, S. M.; Burgess, S. C. AgBase: a functional genomics resource for agriculture. BMC Genomics 2006, 7, 229. (36) Carvalho, P. C.; Yates, J. R., 3rd; Barbosa, V. C. Improving the TFold test for differential shotgun proteomics. Bioinformatics 2012, 28 (12), 1652-4. (37) Lundgren, D. H.; Hwang, S. I.; Wu, L.; Han, D. K. Role of spectral counting in quantitative proteomics. Expert Rev Proteomics 2010, 7 (1), 39-53. (38) Marshall, O.J. Perl Primer: cross-platform, graphical primer design for standard, bisulphite and real-time PCR. Bioinformatics. 2004, 20, 2471–2472. (39) Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2-∆∆CT method. Methods. 2001, 25, 402–408. (40) Tokuoka, T.; Tobe, H. Ovules and seeds in Crotonoideae (Euphorbiaceae): structure and systematic implications. Bot. Jahrb Syst. 1998, 120 (2), 165-86. (41) Weber, H.; Borisjuk, L.; Wobus, U. Molecular physiology of legume seed development. Annu. Rev. Plant Biol. 2005, 56, 253-79. 21

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 33

(42) Nakano, Y.; Ohtani, M.; Polsri, W.; Usami, T.; Sambongi, K.; Demura, T. Characterization of the casbene synthase homolog from Jatropha (Jatropha curcas L.). Plant Biotech. J. 2012, 29 (2), 185-189. (43) Aquino, P. F.; Fischer, J. S. G.; Neves-Ferreira, A. G. C.; Perales, J.; Domont, G. B.; Araujo, G. D. T.; Barbosa, V. C.; Viana, J.; Chalub, S. R. S.; Lima de Souza, A. Q.; Carvalho, M. G.; Leão de Souza, A. D.; Carvalho, P. C. Are gastric cancer resection margin proteomic profiles more similar to those from controls or tumors? J. Proteome Res. 2012, 11, 5836–5842. (44) Hatsugai, N.; Iwasaki, S.; Tamura, K.; Kondo, M.; Fuji, K.; Ogasawara, K.; Nishimura, M.; Hara-Nishimura, I. A novel membrane fusion-mediated plant immunity against bacterial pathogens. Genes Dev. 2009, 23 (21), 2496-506. (45) Nogueira, F. C.; Palmisano, G.; Soares, E. L.; Shah, M.; Soares, A. A.; Roepstorff, P.; Campos, F. A.; Domont, G. B. Proteomic profile of the nucellus of castor bean (Ricinus communis L.) seeds during development. J. Proteomics 2012, 75 (6), 1933-9. (46) Helm, M.; Schmid, M.; Hierl, G.; Terneus, K.; Tan, L.; Lottspeich, F.; Kieliszewski, M. J.; Gietl, C. KDEL-tailed cysteine endopeptidases involved in programmed cell death, intercalation of new cells, and dismantling of extensin scaffolds. Am. J. Bot. 2008, 95 (9), 1049-62. (47) Greenwood, J. S.; Helm, M.; Gietl, C. Ricinosomes and endosperm transfer cell structure in programmed cell death of the nucellus during Ricinus seed development. Proc. Natl. Acad. Sci. U. S. A. 2005, 102 (6), 2238-43. (48) Vartapetian, A. B.; Tuzhikov, A. I.; Chichkova, N. V.; Taliansky, M.; Wolpert, T. J. A plant alternative to animal caspases: subtilisin-like proteases. Cell Death Differ. 2011, 18 (8), 1289-97. (49) Gallardo, K.; Firnhaber, C.; Zuber, H.; Hericher, D.; Belghazi, M.; Henry, C.; Kuster, H.; Thompson, R. A combined proteome and transcriptome analysis of developing Medicago truncatula seeds: evidence for metabolic specialization of maternal and filial tissues. Mol Cell Proteomics 2007, 6 (12), 2165-79. (50) Day, A.; Fenart, S.; Neutelings, G.; Hawkins, S.; Rolando, C.; Tokarski, C. Identification of cell wall proteins in the flax (Linum usitatissimum) stem. Proteomics 2013, 13 (5), 812-25. (51) Minic, Z. Physiological roles of plant glycoside hydrolases. Planta 2008, 227 (4), 723-40. (52) Pelloux, J.; Rustérucci, C.; Mellerowicz, E. J. New insights into pectin methylesterase structure and function. Trends Plant Sci. 2007, 12 (6), 267-277. (53) Rubinstein, B. Regulation of cell death in flower petals. Plant Mol. Biol. 2000, 44 (3), 30318. (54) Laxalt, A. M.; ter Riet, B.; Verdonk, J. C.; Parigi, L.; Tameling, W. I.; Vossen, J.; Haring, M.; Musgrave, A.; Munnik, T. Characterization of five tomato phospholipase D cDNAs: rapid and specific expression of LePLDbeta1 on elicitation with xylanase. Plant J. 2001, 26 (3), 23747. 22

ACS Paragon Plus Environment

Page 23 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

(55) Iakimova, E. T.; Michaeli, R.; Woltering, E. J. Involvement of phospholipase D-related signal transduction in chemical-induced programmed cell death in tomato cell cultures. Protoplasma 2013, 250 (5), 1169-83.

(56) Siedow, J. N. Plant Lipoxygenase: Structure and Function. Annu. Rev. Plant Physiol. Plant Mol. Biol. 1991, 42 (1), 145-188. (57) Liang, Y.; Chen, H.; Tang, M. J.; Yang, P. F.; Shen, S. H. Responses of Jatropha curcas seedlings to cold stress: photosynthesis-related proteins and chlorophyll fluorescence characteristics. Physiol Plant 2007, 131, 508-17. (58) Liu, H.; Liu, Y. J.; Yang, M. F.; Shen, S. H. A comparative analysis of embryo and endosperm proteome from seeds of Jatropha curcas. J Integr Plant Biology 2009, 51, 850-7. (59) Yang, M. F.; Liu, Y. J.; Liu, Y.; Chen, H.; Chen, F.; Shen, S. H. Proteomic analysis of oil mobilization in seed germination and postgermination development of Jatropha curcas. Journal of Proteome Research 2009, 8, 1441-51. (60) Popluechai, S.; Froissard, M.; Jolivet, P.; Breviario, D.; Gatehouse, A. M.; O'donnell, A. G.; Chardot, T.; Kohli, A. Jatropha curcas oil body proteome and oleosins: Lform JcOle3 as a potential phylogenetic marker. Plant Physiol Biochem. 2011, 49, 352-6. (61) Liu, H.; Yang, Z.; Yang, M.; Shen, S. The differential proteome of endosperm and embryo from mature seed of Jatropha curcas. Plant Science 2011, 181, 660-6. (62) Liu, H.; Wang, C.; Komatsu, S.; He, M.; Liu, G.; Shen, S. Proteomic analysis of the seed development in Jatropha curcas: from carbon flux to the lipid accumulation. J Proteomics 2013, 91, 23-40. (63) Booranasrisak, T.; Phaonakrop, N.; Jaresitthikunchai, J.; Virunanon, C.; Roytrakul, S.; Chulalaksananukul, W. Proteomic evaluation of free fatty acid biosynthesis in Jatropha curcas L. (physic nut) kernel development. African Journal of Biotechnology 2013, 12, 3132-42.

23

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 33

Tables Table 1: Classes of seed storage proteins identified in the inner integument of Jatropha curcas seeds. Table 2: Relative expression of selected genes, identified in our proteomic analysis of the inner integument, quantified through qPCR. Positive values represent higher expression in the distal region and negative values represent higher expression in the proximal region. Figures Figure 1: Anatomical structure of J. curcas seed. A. Seed structure 25 DAP. B. Magnified representation of the square area in A showing cell features of the distal and proximal region. Legends: D, distal region; Em, embryo; En, endosperm; II, inner integument; OI, outer integument; P, proximal region; arrowheads indicate cell remnants; arrows indicate vascular bundles. Figure 2: Work flow of the experimental design of the proteomic analysis of the inner integument from J. curcas seeds. Red bar, Proximal region; Yellow bar, Distal region; Green bar, Total inner integument. Figure 3: Gene ontology (GO) categories of the proteins unique to distal and proximal regions of the integument. Plant GO Slim was used to summarize the sub-categories of the identified proteins. Y axes represent percent of proteins. A: GO-Molecular function. B: GO-Cellular component. C: GO-Biological process. Figure 4: TFold analysis pinpointing differentially expressed proteins when comparing the proteomics profiles of the distal and proximal integuments. Each protein is mapped as a dot on the plot according to its –Log2 (p-value) (x-axis) and Log2 (Fold change) (y-axis). Red dots are proteins that satisfy neither the variable fold-change cutoff nor the FDR cutoff α= 0.05. Green dots are those that satisfy the fold-change cutoff but not α. Orange dots are those that satisfy both 24

ACS Paragon Plus Environment

Page 25 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

the fold-change cutoff and α, but are lowly abundant proteins and therefore most likely have their quantitation compromised. Finally, blue dots are those that satisfy all statistical filters. Dots found in the upper section of the plot are up-regulated in the distal sample; likewise, those found in the lower section are up-regulated in the proximal sample.

Supporting Information - this material is available free of charge via http://pubs.acs.org. Supplementary Material Supplementary Figure I: Principal component analysis generated using PatternLab´s Buzios Module. Each blue dot is representative of one analysis from the distal sample; likewise, each orange dot originates from an analysis obtained from a proximal sample. Supplementary Table I: Primer sequences of the target genes used for gene expression analysis. Supplementary Table II: Details about the number of MS, MS2, PSM, peptides and unique peptides. Supplementary Table III: List of the proteins appearing in at least two biological replicates identified in the distal and proximal regions and intact inner integument of developing J. curcas seeds. Supplementary Table VI: Gene ontology (GO) annotation of the proteins identified in the distal and proximal regions and intact inner integument of developing J. curcas seeds. Plant GO Slim was used to summarize sub-categories of the identified proteins.

25

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 33

Supplementary Table V: Gene ontology (GO) annotation of the proteins unique to the distal and proximal regions of the inner integument of developing J. curcas seeds. Plant GO Slim was used to summarize the sub-categories of the proteins. Supplementary Table VI: Differentially expressed proteins identified from distal and proximal regions of the inner integument of developing J. curcas seeds. In the column T, positive values show higher expression in distal region while negative values show higher expression in the proximal region of the inner integument. Supplementary Table VII: Peptidases and peptidase inhibitors identified in distal and proximal regions and intact inner integument of developing J. curcas seeds.

26

ACS Paragon Plus Environment

Page 27 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48

Journal of Proteome Research

Table 1: Classes of seed storage proteins identified in the inner integument of Jatropha curcas seeds.

Protein ID

Description

R. communis (UNIPROT)

Jcr4S01636.40 Jcr4S00279.60 Jcr4S01636.60 Jcr4S01617.40 Jcr4S08024.20 Jcr4S00423.10 Jcr4S03933.20 Jcr4S01793.30

11S globulin seed storage protein 2-like glutelin type-A precursor legumin B precursor nutrient reservoir nutrient reservoir nutrient reservoir nutrient reservoir late embryogenesis abundant protein

B9SW16 B9T5E7 Q9M4Q8 B9SYN7 B9SYN7 B9SYN7 B9SYN7 B9T526

A. thaliana (TAIR) AT1G03890.1 AT5G44120.3 AT5G44120.3 AT2G28680.1 AT1G07750.1 AT2G28680.1 AT1G07750.1 AT2G44060.1

Spectrum count Distal Proximal Total region region Integument 3 13 5 18 28 46 13 25 6 246 45 66 165 33 43 187 42 42 193 52 45 60 41 25

27

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 33

Table 2: Relative expression of selected genes, identified in our proteomic analysis of the inner integument, quantified through qPCR. Positive values represent higher expression in the distal region and negative values represent higher expression in the proximal region. Protein ID Jcr4S00428.20 Jcr4S01104.40 Jcr4S00024.130 Jcr4S01605.60 Jcr4S01651.50 Jcr4S01323.10 Jcr4S00079.140

Description Arabinofuranosidase Cysteine peptidase (KDEL) Cysteine peptidase inhibitor Serine peptidase Subtilisin Subtilisin Serpin

Transcript relative expression Mean Fold Change SD CV (%) - 4.14 1.67 40.46 - 16.64 4.61 27.69 - 3.51 1.21 34.36 - 6,55 3.86 58.96 - 48.56 15.78 32.49 - 20.05 4.44 22.16 1.23 0.18 14.91

28

ACS Paragon Plus Environment

Page 29 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

Anatomical structure of J. curcas seed. A. Seed structure 25 DAP. B. Magnified representation of the square area in A showing cell features of the distal and proximal region. Legends: D, distal region; Em, embryo; En, endosperm; II, inner integument; OI, outer integument; P, proximal region; arrowheads indicate cell remnants; arrows indicate vascular bundles. 187x71mm (150 x 150 DPI)

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Work flow of the experimental design of the proteomic analysis of the inner integument from J. curcas seeds. Red bar, Proximal region; Yellow bar, Distal region; Green bar, Total inner integument. 151x190mm (150 x 150 DPI)

ACS Paragon Plus Environment

Page 30 of 33

Page 31 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

Gene ontology (GO) categories of the proteins unique to distal and proximal regions of the integument. Plant GO Slim was used to summarize the sub-categories of the identified proteins. Y axes represent percent of proteins. A: GO-Molecular function. B: GO-Cellular component. C: GO-Biological process. 255x162mm (150 x 150 DPI)

ACS Paragon Plus Environment

Journal of Proteome Research

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

TFold analysis pinpointing differentially expressed proteins when comparing the proteomics profiles of the distal and proximal integuments. Each protein is mapped as a dot on the plot according to its –Log2 (pvalue) (x-axis) and Log2 (Fold change) (y-axis). Red dots are proteins that satisfy neither the variable foldchange cutoff nor the FDR cutoff α= 0.05. Green dots are those that satisfy the fold-change cutoff but not α. Orange dots are those that satisfy both the fold-change cutoff and α, but are lowly abundant proteins and therefore most likely have their quantitation compromised. Finally, blue dots are those that satisfy all statistical filters. Dots found in the upper section of the plot are up-regulated in the distal sample; likewise, those found in the lower section are up-regulated in the proximal sample. 254x151mm (101 x 101 DPI)

ACS Paragon Plus Environment

Page 32 of 33

Page 33 of 33

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

162x86mm (150 x 150 DPI)

ACS Paragon Plus Environment