Literature DB >> 20110265

Proteomic analysis of pollination-induced corolla senescence in petunia.

Shuangyi Bai1, Belinda Willard, Laura J Chapin, Michael T Kinter, David M Francis, Anthony D Stead, Michelle L Jones.   

Abstract

Senescence represents the last phase of petal development during which macromolecules and organelles are degraded and nutrients are recycled to developing tissues. To understand better the post-transcriptional changes regulating petal senescence, a proteomic approach was used to profile protein changes during the senescence of Petuniaxhybrida 'Mitchell Diploid' corollas. Total soluble proteins were extracted from unpollinated petunia corollas at 0, 24, 48, and 72 h after flower opening and at 24, 48, and 72 h after pollination. Two-dimensional gel electrophoresis (2-DE) was used to identify proteins that were differentially expressed in non-senescing (unpollinated) and senescing (pollinated) corollas, and image analysis was used to determine which proteins were up- or down-regulated by the experimentally determined cut-off of 2.1-fold for P <0.05. One hundred and thirty-three differentially expressed protein spots were selected for sequencing. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) was used to determine the identity of these proteins. Searching translated EST databases and the NCBI non-redundant protein database, it was possible to assign a putative identification to greater than 90% of these proteins. Many of the senescence up-regulated proteins were putatively involved in defence and stress responses or macromolecule catabolism. Some proteins, not previously characterized during flower senescence, were identified, including an orthologue of the tomato abscisic acid stress ripening protein 4 (ASR4). Gene expression patterns did not always correlate with protein expression, confirming that both proteomic and genomic approaches will be required to obtain a detailed understanding of the regulation of petal senescence.

Entities:  

Mesh：

Year:  2010        PMID： 20110265      PMCID： PMC2826652          DOI： 10.1093/jxb/erp373

Source DB:  PubMed          Journal:  J Exp Bot        ISSN： 0022-0957            Impact factor:   6.992


Introduction

Many flowers have large, colourful petals (collectively called the corolla) that serve to attract pollinators. The longevity of flower petals varies considerably by species and is, in part, determined by the frequency of pollinator visitation (Ashman and Schoen, 1994). The maintenance of flower petals is costly in terms of metabolic energy and water loss, therefore it is advantageous for the plant to induce corolla senescence as soon as the flower has been successfully pollinated (Stead, 1992). The lifespan of most pollinated flowers is terminated by petal wilting followed by corolla abscission. During flower senescence, macromolecules are degraded and organelles are systematically dismantled so that nutrients can be remobilized to developing tissues. This process of cellular disassembly allows the plant to recover carbon, nitrogen, and phosphorus from the petals before the corolla is shed (Verlinden, 2003; Jones, 2004; Stead ; Chapin and Jones, 2007, 2009). Studies of both leaf and flower senescence have focused almost entirely on the identification of senescence-associated genes (reviewed in Buchanan-Wollaston ; Jones, 2004; Stead ; van Doorn and Woltering, 2008). While these studies have indicated that the timely initiation and progression of the senescence programme requires new gene transcription, senescence is also controlled at the post-transcriptional level (Thomas and Stoddart, 1980). Advances in mass spectrometry (MS) and 2-dimensional gel electrophoresis (2-DE) techniques and the availability of continually expanding EST and genomic databases have provided new opportunities for studying plant proteomes. Comparative proteomic analysis has been used to identify protein changes during plant development and following abiotic and biotic stresses (reviewed in Jorrin ). Recently, plant proteomics has been used to study various developmental processes including leaf senescence, fruit ripening, and early flower development (Wilson ; Dafny-Yelin ; Carp and Gepstein, 2007; Faurobert ; Hebeler ; Desclos ). Drought, heat, water logging, and salt-induced stress responses have also been investigated using comparative proteomic approaches (Askari ; Horvath-Szanics ; Huang ; Ahsan ; Hajheidari ; Ito ; Jiang ; Veerasamy ; Yang ). These proteomic studies have provided a high throughput means of confirming the involvement of biochemical pathways that have been identified in genomic experiments and have provided unique insights into the role of post-translational modifications regulating and executing plant development and stress responses. The premature senescence of floral tissues can have a detrimental impact on the quality, yield and subsequent value of agricultural and horticultural crops. It is therefore important to understand the molecular and biochemical mechanisms of senescence initiation, regulation, and execution in plants. Petunias provide an excellent model system for studies of flower senescence because they flower profusely and have large floral organs that are amenable to molecular and biochemical analysis. While the longevity of petunia corollas is under tight developmental control, petal senescence can be accelerated and synchronized by treating flowers with exogenous ethylene or by pollination (Jones ). Changes in dry weight and petal nutrient content indicate that remobilization is occurring during petunia corolla senescence (Verlinden, 2003; Chapin and Jones, 2007, 2009; Jones, 2008). The aim of this study was to identify proteins that increased and decreased in abundance during petunia corolla senescence. 2-DE, mass spectrometry, and bioinformatics tools were used to identify and functionally classify proteins that were differentially expressed in senescing and non-senescing petals. These studies provide insight into the dynamic changes occurring in the corolla proteome during senescence, and further support the idea that a combination of genomic and proteomic approaches will be necessary to elucidate molecular and biochemical mechanisms of senescence.

Materials and methods

Plant materials

Petunia×hybrida ‘Mitchell Diploid’ (MD) plants were used in all experiments. Seeds were sown in cell-packs on top of soil-less mix (Promix BX, Premier Horticulture, Quebec, Canada). All plants were established in the greenhouse after germination and plants were transferred to 16 cm pots after 4 weeks. Plants were fertilized at each watering with 150 mg l−1 Scott's Excel 15N-5P-15K (The Scotts Co., Marysville, OH), and a one-time treatment of Soluble Trace Element Mix (STEM, The Scotts Co., Marysville, OH) was applied 4 weeks after transferring to 16 cm pots. Growing conditions were 24/16 °C (day/night) with a 13 h photoperiod supplemented by high pressure sodium and metal halide lights. Petunia petals (collectively called the corolla) were collected from unpollinated and pollinated flowers at 0, 24, 48, and 72 h after flower opening. Three replicates, each containing eight corollas from at least three different plants, were collected at each time point. The eight corollas in a replicate were pooled prior to protein or RNA extraction.

Protein extraction and 2-DE separation of proteins

Approximately 2.5 g corolla tissue (eight corollas) was ground in liquid nitrogen and mixed with 8000 μl homogenization buffer [100 mM HEPES-KOH (pH 7.5), 5% (v/v) glycerol, 15 mM EGTA, 5 mM EDTA, 0.5% (v/v) polyvinylpyrrolidone, 3 mM dithiothreitol, 60 μl proteinase inhibitor cocktail (Sigma-Aldrich, St Louis, Mo)]. Total soluble proteins were extracted as described by Coaker . The resulting pellet was resuspended in 800 μl of rehydration buffer [8 M urea, 2 M thiourea, 4% (w/v) CHAPS, 0.2% (v/v) immobilized pH gradient (IPG) buffer (pH 5–8; Bio-Rad Laboratories, Hercules, CA), 1% (v/v) Triton X-100, 0.5% (w/v) DTT] and sonicated for 2 min in a bath sonicator (Model PC5; L&R Manufacturing Co., Kearny, NJ). After shaking at 120 rpm for 1.5 h at 25 °C, the protein concentration was determined with the Amersham Biosciences 2-D quantification kit (Amersham Biosciences Corp., Piscataway, NJ). Two hundred micrograms of total protein for each sample was used for passive rehydration of 11 cm IPG strips (pH 5–8; Bio-Rad Laboratories, Hercules, CA). Programmed isoelectric focusing (IEF) was performed with the Protean IEF Cell (Bio-Rad Laboratories, Hercules, CA) using the following preset conditions: start voltage 0 V, end voltage 8000 V, total volt-hours 35 000 V h at 20 °C. After IEF, IPG strips were immediately equilibrated in buffer 1 [6 M urea, 0.375 M TRIS (pH 8.8), 2% (w/v) SDS, 20% (v/v) glycerol, 2% (w/v) DTT] and buffer 2 (6 M uUrea, 0.375 M TRIS (pH 8.8), 2% (w/v) SDS, 20% (v/v) glycerol, 2.5% (w/v) iodoacetamide], each for 10 min. The equilibrated strips were then subjected to SDS-PAGE (12.5%) for second dimension separation. After running for 2 h under a constant 200 V, gels were stained with GelCode Blue Stain Reagent (Pierce, Rockford, IL).

Data analysis of 2-D gels

Protein profiles of pollinated corollas were compared to unpollinated corollas at the same developmental age to identify differentially expressed proteins for sequence analysis. Proteins up-regulated during senescence were defined as those where the pixel intensities associated with a spot increased in pollinated corollas compared to unpollinated corollas, or those proteins that were detected in pollinated corollas but not in unpollinated corollas. Down-regulated proteins were defined as those where the pixel intensities decreased in unpollinated corollas compared to the pollinated corollas, or proteins that were detected only in unpollinated corollas. Three 2-D gels from each time point, representing independent replications of eight corollas each, were analysed using PDQuest v. 7.40 (Bio-Rad Laboratories, Hercules, CA, USA). For each time point the three replicate gel images for pollinated and unpollinated corollas were imported as a dataset, and a synthetic Gaussian image from one of three replicate gel images was created as a reference for all subsequent quantification and protein analyses. The quantities of protein spots on each image were normalized based on the total density in the gel image. A scaling factor was used to give a more meaningful value and the normalized quantity was multiplied by parts per million (ppm). Protein spots were visually inspected to validate all automated analyses. Average intensities and correlation coefficients between gels were obtained. The distribution of protein intensities and residuals was inspected using the Univariate procedure of SAS (Version 9.0, SAS Institute Inc., Cary, NC, USA). Transformation of the data was necessary in order to approximate a normal distribution. Both loge [LN(value+1)] and log10 [log(value+1)] transformations of the data approximated normal distributions, and the loge (LN) transformation was subsequently applied across all quantitative data prior to statistical analysis. Analysis of variance (ANOVA) was conducted using the general linear model procedure in the SAS software. The statistical model included main effects for protein, treatment (pollinated or unpollinated), time, and replicate. All factors were considered fixed. Sufficient degrees of freedom were available to test all two-way interactions which were included in the statistical model for ANOVA. The Mean Square Error from the ANOVA was used to estimate an experiment-wide cut-off for significance, calculated as fold-change using the relationship e[t×(Sqrt(MSE/n)]/2 where t was 1.96 for P=0.05, 2.576 for P=0.01, and 3.291 for P=0.001; where n equals the number of replicates (Kerr ).

Identification of proteins by mass spectrometry

Mass spectrometry was conducted at the Cleveland Clinic Proteomics Laboratory (Cleveland Clinic Foundation, Cleveland, OH), and trypsin digests of the proteins of interest were prepared as described by Kinter and Sherman (2000). Protein spots were excised from the SDS-PAGE gels, divided into a number of smaller pieces and washed/destained in 50% ethanol/5% acetic acid. The gel pieces were dehydrated in acetonitrile and dried in a Speed-Vac. Digestion was carried out overnight at room temperature with 5 μl of 20 ng μl−1 trypsin in 50 mM ammonium bicarbonate. The resulting peptides were extracted twice with 50% acetonitrile/5% formic acid, and the extracts were combined and evaporated to <30 μl for LC-MS analysis. A ThermoScientific LCQ ion trap mass spectrometer (San Jose, CA) interfaced with a self-packed Phenomenex Jupiter C18 reversed-phase capillary column was used for the analysis. Two μl volumes of the trypsin digested extract were injected and peptides were eluted from the column by a 30 min binary gradient with 0.05 M acetic acid and acetonitrile as the mobile phases at a flow rate of 250 nl min−1. The microelectrospray ion source was operated at 2.5 kV. Tryptic digests were analysed using the multi-task capability of the instrument to acquire full scan mass spectra to determine peptide molecular weights and product ion spectra to determine amino acid sequences in successive instrument scans. All collisionally induced dissociation (CID) spectra were searched against the National Center for Biotechnology Information (NCBI) green plant protein database using the Mascot search program (Matrix Science, Boston, MA). An ion score greater than 25 was used as a cut-off for consideration, and all identifications were based on at least two peptides that could then be verified by manual inspection of the spectra. If a protein could not be identified using this method, the CID spectra were manually interpreted to determine the specific amino acid sequence of the peptides. The peptide sequence data were then subjected to either FASTA searches performed on the University of Virginia website (http://fasta.bioch.virginia.edu/fasta_www2/fasta_list2.shtml) to identify the homologous NCBI sequences or to FASTA searches of translated EST databases using the search program Sequest (ThermoScientific, San Jose, CA). The translated EST databases were created by The Ohio State University Molecular and Cellular Imaging Center (Wooster, OH) and consisted of Solanaceae ESTs translated in six reading frames. Petunia ESTs were obtained from the SOL Genomics Network (SGN) (http://sgn.cornell.edu/; Mueller ) and Nicotiana benthamiana and Solanum tuberosum ESTs were provided by Dr Sophien Kamoun (Sainsbury Laboratory, John Innes Centre, UK). Xcorr values for each identified peptide were filtered according to the peptide charge states Xcorr ≥1.5 for +1, Xcorr ≥2.0 for +2, and Xcorr ≥2.5 for +3. To validate all identifications, the interpreted peptide sequences were manually aligned with the homologous NCBI protein sequences (see Supplementary Fig. S1 at JXB online).

Protein expression clustering and functional classification

Cluster (version 2.11) hierarchical clustering software (http://rana.lbl.gov/EisenSoftware.htm) was used to identify the main classes of protein expression (Eisen ). Protein expression patterns were clustered as described in the software user manual. Briefly, input data were normalized by dividing the loge values by the intensity at 0 h. Weighted Pearson (centred) correlation was chosen to define the similarity matrix. A dendrogram was generated by running the average-linkage clustering algorithm. The clustering image was visualized using the TreeView software (http://rana.lbl.gov/EisenSoftware.htm). Protein sequences were searched against the NCBI databases using the tBLASTn program to identify the most similar Arabidopsis genes. The locus identifier for this gene was then used to search the TAIR gene ontology (GO) annotation database (http://www.arabidopsis.org/tools/bulk/go/index.jsp) to identify the putative biological function of the protein of interest (Berardini ).

Real-time PCR analysis of gene expression

Total RNA was extracted using TRIzol reagent (Invitrogen, Carlsbad, CA) and treated with RQ1 RNase-free DNase (Promega, Madison, WI). cDNA was synthesized from 2 μg total RNA using the Omniscript RT Kit (Qiagen, Valencia, CA). The gene sequences for selected proteins were identified from the SGN petunia EST database (Mueller ) and gene-specific primers were designed using Primer Quest software available through IDT SciTools (http://www.idtdna.com/Scitools/Applications/Primerquest/) (Table S1). If petunia sequences corresponding to the peptides could not be identified from the EST database, degenerate primers were constructed to the most distal peptide sequences to allow for RT-PCR amplification of the largest amount of sequence from the corresponding gene. The PCR amplicon was then sequenced and Primer Quest software was used to design primers for real-time PCR. Quantitative PCR was performed in a 20 μl reaction volume using B-R SYBR Green SuperMix for iQ (Quanta BioSciences, Inc., Gaithersburg, MD). One microlitre cDNA was used as template, and all reactions were performed in triplicate. After a 3 min activation step at 95 °C, qPCR was conducted for 40 cycles of 94 °C for 10 s, 60 °C for 30 s, and 72 °C for 30 s using the iQ5 Thermocycler (Bio-Rad, Hercules, CA). The annealing temperature for some primer pairs varied and the specific temperatures are listed in Supplementary Table S1 at JXB online. Melt curves were generated to check amplification specificity. The absolute starting quantity of the target genes was normalized to that of PhACTIN1 in each cDNA sample to calculate relative target gene expression. PhACTIN1 had very consistent transcript abundance in corollas from 0 h to 72 h after pollination (ML Jones, unpublished results), and is not the same petunia actin that showed senescence-related increases in protein abundance in this study.

Results

Protein profiling during petunia corolla senescence

Pollination accelerated the senescence of petunia corollas, and flowers were visibly wilted by 48 h after pollination (48 P) (Fig. 1A). Corollas at 72 h after pollination (72 P) were severely wilted, but were not yet drying at the petal margins. Unpollinated petunia corollas from 24–72 h after flower opening (24 U, 48 U, and 72 U) were fully turgid, showed no senescence symptoms and were visually indistinguishable from flowers at anthesis (0 U). The 2-D separation of total soluble proteins from petunia corollas resulted in the detection of approximately 600 distinct protein spots in each gel (Fig. 1B). The molecular masses (Mr) of the proteins ranged from 15 kDa to 110 kDa (Fig. 2).
Fig. 1.

Protein changes during petunia corolla development and senescence. (A) Pollination accelerates petunia flower senescence. Petunia×hybrida ‘Mitchell Diploid’ (MD) unpollinated (U) flowers at 0, 24, 48, and 72 h and pollinated (P) flowers at 24, 48, and 72 h. (B) 2-D gel profiling of protein changes corresponding to corolla senescence. Representative 2-D gels of petunia corolla proteome variation during pollination-induced senescence. The first dimension was performed using 200 μg total soluble proteins on linear gradient IPG strips with pH 5–8. Second dimension separation was conducted on 12% SDS-PAGE gels and proteins were visualized using GelCode blue staining. Approximately 600 spots were detected on each gel. Three biological replicates at each time point were conducted.

Fig. 2.

Representative 2-DE comparisons of unpollinated and pollinated gels at 72 h. (A) All of the down-regulated proteins are indictated with arrows on the 72 U gels. (B) All of the up-regulated proteins are indicated with arrows on 72 P gels. The spot numbers are the same as the sample numbers (without the dash in the middle) in Tables 2 and 3. The assigned identities of up- and down-regulated proteins and their expression data are also shown in Tables 2 and 3.

Protein changes during petunia corolla development and senescence. (A) Pollination accelerates petunia flower senescence. Petunia×hybrida ‘Mitchell Diploid’ (MD) unpollinated (U) flowers at 0, 24, 48, and 72 h and pollinated (P) flowers at 24, 48, and 72 h. (B) 2-D gel profiling of protein changes corresponding to corolla senescence. Representative 2-D gels of petunia corolla proteome variation during pollination-induced senescence. The first dimension was performed using 200 μg total soluble proteins on linear gradient IPG strips with pH 5–8. Second dimension separation was conducted on 12% SDS-PAGE gels and proteins were visualized using GelCode blue staining. Approximately 600 spots were detected on each gel. Three biological replicates at each time point were conducted. Representative 2-DE comparisons of unpollinated and pollinated gels at 72 h. (A) All of the down-regulated proteins are indictated with arrows on the 72 U gels. (B) All of the up-regulated proteins are indicated with arrows on 72 P gels. The spot numbers are the same as the sample numbers (without the dash in the middle) in Tables 2 and 3. The assigned identities of up- and down-regulated proteins and their expression data are also shown in Tables 2 and 3.
Table 2.

Functional classification of petunia proteins that were up-regulated during pollination-induced corolla senescence

Sample no.aPutative protein assignmentbMatched peptidescMatched speciesAccession/EST no.dObser. Mr/pIeTheor. Mr/pIfPeptide coveragegMean of intensitiesh
48 Ui48 Pj72 U72 P
Defence response and response to stress
50-9Plasma membrane polypeptide16 (4,12)Nicotiana tabacum276499224/5.122.9/5.062%ud k1068ud470
14-23Manganese superoxide dismutase8 (4,4)Gossypium hirsutum321935326/6.922.1/8.553%271701111984
14-52Manganese superoxide dismutase7 (4,3)Gossypium hirsutum321935328/6.622.1/8.545%3635631211078
36-16Iron superoxide dismutase8 (6,2)Solanum lycopersicum3341330326/6.027.9/6.634%3717661861106
14-51Superoxide dismutase; Nectarin 1 precursor6 (0,6)Nicotiana plumbaginifolia609082930/5.824.8/7.729%5574558732463
36-12Annexin17 (11,6)Capsicum annuum107166037/6.335.8/5.749%137513455791477
36-23Similar to NtPRp27-like protein11 (11,0)Petunia×hybridaSGN-E52081224/7.729.5/7.9na156108235898
50-4Auxin-responsive-like protein; DFL19 (8,1)Arabidopsis thaliana1523965360/6.868.9/5.516%ud364ud240
36-10Mitochondrial formate dehydrogenase precursor20 (17,3)Solanum tuberosum1199152744/7.242.0/6.656%2291246811732717
50-7Similar to abscisic stress ripening protein2 (2,0)Petunia×hybridaSGN-E52444428/5.926.7/4.8naud202ud101
14-9UTP-glucose-1-phosphate uridyltransferase21 (13,8)Solanum tuberosum2159951/6.951.8/5.446%ud253ud364
Electron transport or energy pathways
36-5Succinate dehydrogenase flavoprotein alpha subunit15 (10,5)Arabidopsis thaliana884373449/6.369.7/5.931%894766250502
14-221,4-benzoquinone reductase-like protein8 (4,4)Arabidopsis thaliana2153948127/6.521.8/6.045%3885162101014
36-171,4-benzoquinone reductase-like protein16 (5,11)Arabidopsis thaliana2153948126/6.521.7/6.072%8988786541322
36-30Similar to plastocyanin-like domain-containing protein; early nodulin-like protein 2 precursor3 (1,2)Petunia×hybridaSGN-E52855016/5.629.9/9.4na171582ud1067
14-25Photosystem II 23 kDa protein5 (0,5)Solanum lycopersicum1931722/5.727.8/8.335%ud753ud576
Proteolysis
14-50Cysteine proteinase P218 (7,1)Petunia×hybrida94508130/5.639.3/5.823%1670215012974919
36-15Proteasome delta-subunit10 (8,2)Nicotiana tabacum174335627/5.925.1/5.247%291509319640
14-27Serine carboxylase II-2 (N-term truncated)8 (1,7)Oryza sativa5577386120/7.153.4/5.518%444431361062
36-13Glucose acyltransferase; serine carboxypeptidase-like protein7 (3,4)Solanum berthaultii410170335/5.853.4/5.118%221605144719
36-24Ubiquitin-conjugating enzyme E2; ubiquitin-protein ligase3 (2,1)Arabidopsis thaliana3069365621/6.716.6/6.224%ududud927
Carbohydrate metabolism
44-7Beta-xylosidase 2; LeXYL213 (8,5)Solanum lycopersicum3735970865/6.268.9/8.020%3513456524356521
14-13Beta-xylosidase 2; LeXYL1(C-terminally truncated)16 (8,8)Solanum lycopersicum3735970843/6.468.9/8.035%ud925ud2601
36-3Beta-xylosidase 1; AtBXL17 (2,5)Arabidopsis thaliana1523986773/6.283.5/8.811%ud404ud575
14-2Alpha-mannosidase8 (3,5)Arabidopsis thaliana1017713071/6.7118/8.37%1593942021187
14-29Vacuolar invertase (N or C-terminally truncated)14 (8,6)Nicotiana tabacum2989306448/5.571.2/5.419%89610317892828
36-22Vacuolar invertase (N-terminally truncated)5 (3, 2)Nicotiana tabacum2989306422/5.671.2/5.47%ududud455
36-21Vacuolar invertase precursor (N-terminally truncated)2 (0,2)Nicotiana tabacum2989306424/6.971.2/5.44%ud361ud765
14-11Putative beta-galactosidase (C-terminally truncated)13 (8,5)Solanum lycopersicum793962346/6.793.2/6.814%166825307773531
50-3Cytosolic phosphoglucomutase21 (15,6)Solanum tuberosum825062459/6.763.5/6.044%86514865451315
Lipid metabolism
14-6Lipoxygenase (N-terminally truncated)21 (9,12)Nicotiana tabacum89934462/5.997.6/5.524%1695901262226
50-2Lipoxygenase8 (5,3)Nicotiana tabacum89934465/5.797.6/5.512%ud560ud368
14-20Putative GDSL- motif lipase4 (1,3)Vitis vinifera3778982528/6.119.1/5.525%2101843951078
36-14Proline-rich protein; SGNH plant-lipase-like8 (1,7)Oryza sativa5091054739/5.540.1/5.733%39211925391478
Cytoskeleton organization and biogenesis
14-30Actin20 (17,3)Arabidopsis thaliana2046586547/5.641.6/5.360%100126100674
50-10Actin-depolymerizing factor 13 (3,0)Petunia×hybrida1490621922/6.316.0/5.819%ud656ud1362
51-14Actin-depolymerizing factor 27 (7,0)Petunia×hybrida1490621021/6.916.0/5.852%353810182551
Nucleotide/nucleoside/nucleobase metabolism
14-12Endonuclease8 (4,4)Solanum tuberosum5065759643/5.734.4/5.632%ud421ud1549
Amino acid metabolism
36-32Methylene-tetra hydrofolate reductase; MTHFR1(C-terminally truncated)4 (1,3)Arabidopsis thaliana3069509715/5.645.6/6.210%ud271ud659
Nitrogen metabolism
36-33Beta-Ureidopropionase; PYD33 (2,1)Arabidopsis thaliana3069800948/6.445.5/5.913%ud164ud432
Unclassified
36-27Putative acid phosphatase (N-terminally truncated)4 (0,4)Oryza sativa5529647719/6.032.3/5.49%ud154ud511
36-19Putative Kunitz-type proteinase inhibitor3 (1,2)Petunia×hybridaSGN-E53618224/5.723.8/7.0na2166593711010
36-20Putative Kunitz-type proteinase inhibitor5 (2,3)Petunia×hybridaSGN-E53618223/6.123.8/7.0na251637981720
14-24Putative aluminium-induced protein6 (1,5)Arabidopsis thaliana2153724624/6.527.5/6.425%186569255485
36-28Hypothetical protein6 (5,1)Petunia×hybridaSGN-E52163918/5.624/9.4na368751289960
50-11S-adenosyl-L-methionine:benzoic acid/salicylic acid carboxyl methyltransferase; PhBSMT221 (21,0)Petunia×hybrida2862949740/6.540.7/5.748%353810182551

Sample number indicates the spot sample for in-gel digestion in the sequencing reports. These sample numbers (without dash) were also used to indicate spots on the 2-D gels in Fig. 2.

Proteins were classified based on the biological processes according to homology to genes in the TAIR Gene Ontology database.

Matched peptides indicate total number of peptides that matched to other proteins. The first and the second number in the parenthesis indicate the peptides with exact and homologous matches, respectively.

Accession/EST number indicates the sequence in the searched database identified from either NCBI non-redundant protein database or SGN petunia EST database (http://www.sgn.cornell.edu/index.pl).

Obser. Mr/pI, observed Mr/pI for each protein was calculated from the 2-D gels with Image PDQuest 7.4.0 software according to standard marker proteins.

Theor. Mr/pI, theoretical Mr/pI of the matched proteins.

Peptide coverage refers to the percentage of sequence coverage of the matched protein; no available (na) data when EST database was used.

Mean of intensities is the average of three replicated spot normalized intensities at 48 U, 48 P, 72 U, and 72 P.

U is hours after flowers are fully opened without pollination (unpollinated).

P is hours after flowers are fully opened with pollination.

ud, spot intensity was undetectable on the gel.

Table 3.

Functional classification of petunia proteins that were down-regulated during pollination-induced corolla senescence

Sample no.aPutative protein assignmentbMatched peptidescMatched speciesAccession/EST No.dObser. Mr/pIeTheor. Mr/pIfPeptide coveragegMean of intensitiesh
48 Ui48 Pj72 U72 P
Defence response and response to stress
44-123-Phosphoshikimate 1-carboxy vinyltransferase19 (19,0)Petunia×hybrida11417649/6.555.5/8.039%24425811665196
51-11Putative pyridoxine biosynthesis isoform B11 (11,0)Nicotiana tabacum4639927136/6.633.1/5.932%1100465671411
49-24Copper-zinc superoxide dismutase5 (3,2)Solanum lycopersicum1344591818/6.014.8/5.345%3037151026292511
Amino acid metabolism
14-31Methionine synthase22 (16,6)Solanum tuberosum843954580/6.684.7/5.932%349314771322230
44-1Methionine synthase2 (2, 0)Solanum tuberosum843954573/5.784.7/5.95%576ud k322ud
44-2Methionine synthase43 (24, 19)Solanum tuberosum843954581/6.784.7/5.963%7428348542301870
44-3Methionine synthase24 (17,7)Solanum tuberosum843954580/6.684.7/5.930%1556739363325
14-37S-adenosylmethionine synthetase19 (19,0)Petunia×hybrida55950645/6.243.8/5.573%276417862883889
14-36S-adenosylmethionine synthetase21 (14,7)Petunia×hybrida572659445/6.043.7/5.471%17485382449220
44-13S-adenosylmethionine synthetase16 (8,8)Nicotiana tabacum3334051747/6.543.1/5.855%16316811520888
44-14S-adenosylmethionine synthetase12 (10,2)Nicotiana tabacum3334051748/6.643.1/5.837%18596761464486
49-18Hydroxymethyltransferase28 (12,16)Arabidopsis thaliana224474955/7.951.7/6.865%1968593476299
Cytoskeleton organization and biogenesis
14-34Alpha tubulin17 (17,0)Physcomitrella patens2539654548/5.650.1/5.040%8967571503431
49-21Actin-depolymerizing factor 17 (7,0)Petunia×hybrida1490621923/5.816.0/5.846%392ud308ud
49-22Actin-depolymerizing factor 210 (10,0)Petunia×hybrida1490621022/6.416.5/5.876%604ud386ud
Carbohydrate metabolism
51-13Xyloglucan endotransglucosylasehydrolase; XTH712 (3,9)Solanum lycopersicum4279546637/6.833.4/7.647%291019511650409
44-15Ribulose bisphosphate carboxylase large chain precursor16 (15,1)Petunia×hybrida13201154/6.953.0/6.634%5642421429041441
44-26Ribulose 1,5-bisphosphate carboxylase small chain precursor8 (8,0)Petunia×hybrida13208316/6.620.4/8.352%282212122384799
44-21Caffeoyl CoA 3-O-methyltransferase15 (10,5)Betula platyphylla5763962933/5.627.8/5.366%2806192038471187
44-22Caffeoyl CoA 3-O-methyltransferase13 (8,5)Solanum tuberosum2474596932/5.827.2/5.366%464ud499ud
44-23Caffeoyl CoA 3-O-methyltransferase15 (8,7)Solanum tuberosum2474596931/5.827.2/5.377%297996321561110
Nitrogen metabolism
44-18Glutamine synthetase shoot isozyme17(9,8)Oryza sativa5091251142/5.939.2/5.558%4066200736271869
Electron transport or energy pathways
44-4Transketolase 110 (9,1)Capsicum annuum355981477/6.480/6.213%1154553935457
49-20Photosystem II 23 kDa protein14 (6,8)Solanum lycopersicum1931727/5.627.8/8.362%819ud697ud
Unclassified
51-12Putative aluminium-induced protein9 (3,6)Arabidopsis thaliana1523999333/6.827.5/6.434%510ud382ud
44-5Similar to Acyl activating enzyme 112 (2,0)Petunia×hybridaSGN-E52340066/7.313.2/8.5na1161481842493
14-40S-adenosyl-L-methionine:benzoic acid/ salicylic acid carboxyl methyltransferase; PhBSMT217 (17,0)Petunia×hybrida2862949740/6.240.7/5.737%5071201824001745

Sample number indicates the spot sample for in-gel digestion in the sequencing reports. These sample numbers (without dash) were also used to indicate spots on the 2-D gels in Fig. 2.

Proteins were classified based on the biological processes according to homology to genes in the TAIR Gene Ontology database.

Matched Peptides indicate total number of peptides that matched to other proteins. The first and the second number in the parenthesis indicate the peptides with exact and homologous matches respectively.

Accession/EST number indicates the sequence in the searched database identified from either NCBI non-redundant protein database or SGN petunia EST database (http://www.sgn.cornell.edu/index.pl).

Obser. Mr/pI, observed Mr/pI for each protein was calculated from the 2-D gels with Image PDQuest 7.4.0 software according to standard marker proteins.

Theor. Mr/pI, theoretical Mr/pI of the matched proteins.

Peptide coverage refers to the percentage of sequence coverage of the matched protein; no available (na) data when EST database was used.

Mean of intensities is the average of three replicated spot normalized intensities at 48 U, 48 P, 72 U, and 72 P.

U is hours after flowers are fully opened without pollination (unpollinated).

P is hours after flowers are fully opened with pollination.

ud, spot intensity was undetectable on the gel.

Pixel intensity data of protein spots were not normally distributed, though both natural log and log base 10 transformations were effective at creating a normal distribution of pixel intensity data. Analysis of variance (ANOVA) indicated that the main effects for differences between proteins (P <0.0001), pollination treatment (P <0.0001), and time points (P <0.0001) were all highly significant. The replication effect and all two-way interactions with replicate were non-significant (0.760

1.00). The treatment by protein and time point by protein interactions were both highly significant (P <0.0001), indicating differential expression of individual proteins due to pollination and time. Based on the mean squared error associated with the ANOVA, the experimental-wide cut-off for differentially expressed spots was 2.1-fold at P <0.05, 3.4-fold at P <0.01, and 3.6-fold at P <0.001 (Kerr ). This analysis supplemented the identification of quantitative differences between treatments (pollinated versus unpollinated) by a t test (P ≤0.05) using the PDQuest software. In addition, qualitative differences (presence/absence) were also detected. No significant protein changes were identified at 24 h when comparing protein profiles in pollinated (24 P) and unpollinated (24 U) corollas (Table 1). By 48 h after flower opening, 72 proteins were up-regulated and 41 proteins were down-regulated in pollinated corollas by greater than 2.1-fold (P <0.05), suggesting that about 20% of the abundant proteins were differentially regulated. One-hundred-and-twelve proteins were up-regulated and 57 proteins were down-regulated at 72 h (P <0.05). This represented approximately 30% of the abundant proteins detected on the 2-D gels.

Table 1.

Profiling of differentially expressed proteins during petunia corolla senescence

Time points and treatments
Differentially expressed proteinsa24 P versus 24 U
48 P versus 48 U
72 P versus 72 U
No.Total (%)No.Total (%)No.Total (%)
Up-regulated proteinsb00731311220
Down-regulated proteinsc004175710
Total001142016930

Significant differences among treatments within a time point were detected by the t-test (P ≤0.05).

Proteins were defined as up-regulated if they were not detected in unpollinated petals and were detected in pollinated petals or if they increased in abundance in pollinated petals.

Proteins were defined as down-regulated if they were detected in unpollinated but not in pollinated petals or if they decreased in abundance in pollinated petals.

Profiling of differentially expressed proteins during petunia corolla senescence Significant differences among treatments within a time point were detected by the t-test (P ≤0.05). Proteins were defined as up-regulated if they were not detected in unpollinated petals and were detected in pollinated petals or if they increased in abundance in pollinated petals. Proteins were defined as down-regulated if they were detected in unpollinated but not in pollinated petals or if they decreased in abundance in pollinated petals. From the differentially expressed proteins, 133 unique spots were selected for sequencing by tandem mass spectrometry (Fig. 2A, B). Forty-six up- and 27 down-regulated spots were identified as single proteins (Tables 2, 3). Sixty-nine per cent of these proteins were identified based on peptide matches from Solanaceae sequences, with 26% of those being from petunia. A total of 47 spots included peptides that corresponded to multiple proteins (data not shown). Since the expression pattern could not unambiguously be determined for each individual protein, the data for spots with multiple protein components are not reported in this paper. Thirteen (11 up- and two down-regulated) spots could not be identified from any database (data not shown). The experimental Mr of many of the differentially regulated proteins was similar to the theoretical Mr of the proteins from the matched species (Tables 2, 3). Some protein spots had a slightly higher Mr than the theoretical Mr, and these mass differences may be due to various post-translational modifications during senescence. More striking was the number of proteins that had an experimental Mr that was much less than the theoretical Mr. A number of these proteins were truncated at the N- or C-termini, as the peptides identified by MS only aligned with a portion of the protein. These may represent proteins that have been degraded during senescence or proenzymes that have been cleaved and activated during the progression of senescence. While protease inhibitors were included in the protein extraction buffers, truncated proteins may also have resulted from in vitro protein degradation. Proteins that were identified as N- or C-terminally truncated due to MS peptide alignments and low observed Mr included the up-regulated proteins serine carboxylase II-2 (14-27), beta-xylosidase 2 (14-13), vacuolar invertase (14-29, 36-22), beta-galactosidase (14-11), lipoxygenase (14-6), methylenetetrahydrofolate reductase (36-32), cysteine protease (14-50), and acid phosphatase (36-27). Many of these truncated proteins were detected only in senescing corollas (i.e. 48 P and 72 P). None of the down-regulated proteins were identified by the peptide coverage as being N- or C-terminally truncated.

Functional classification of identified proteins based on biological processes

The 46 up-regulated proteins were grouped into 10 functional classifications and the 27 down-regulated proteins were grouped into seven functional classifications (Tables 2, 3). The up-regulated proteins were classified as defence and stress response (11), carbohydrate metabolism (9), electron transport or energy pathways (5), proteolysis (5), lipid metabolism (4), cytoskeleton organization and biogenesis (3), nucleotide/nucleoside/nucleobase metabolism (1), amino acid metabolism (1), nitrogen metabolism (1) or unclassified proteins (6) (Fig. 3A). The down-regulated proteins were classified as amino acid metabolism (9), carbohydrate metabolism (6), defence and stress response (3), cytoskeleton organization and biogenesis (3), electron transport or energy pathways (2), nitrogen metabolism (1), or unclassified proteins (3) (Fig. 3B).
Fig. 3.

Biological processes classification of up- (A) and down-regulated (B) proteins. Forty-six up- and 26 down-regulated proteins were assigned to putative biological process categories according to the TAIR GO annotation database (http://www.arabidopsis.org/tools/bulk/go/index.jsp).

Biological processes classification of up- (A) and down-regulated (B) proteins. Forty-six up- and 26 down-regulated proteins were assigned to putative biological process categories according to the TAIR GO annotation database (http://www.arabidopsis.org/tools/bulk/go/index.jsp).

Clustering analysis of protein expression patterns

During petal senescence, certain groups of proteins will probably be co-regulated. Hierarchical clustering revealed three groups, which included 24 U and 24 P, 48 U and 72 U, and 48 P and 72 P; based on expression similarities (Fig. 4) and corresponding to developmental stages from flower opening to senescence. The two main clusters included up- and down-regulated proteins. There were eight and six subclusters of up-regulated and down-regulated proteins, respectively. The relative abundance of a representative protein from each subcluster is shown in Fig. 5. The largest up-regulated subcluster included proteins (36-13 to 44-7) with the highest abundance at 48 and 72 h after pollination. Many of these represented senescence-specific proteins that were detected only in 48 P and 72 P corollas. There were also a few proteins that were detected only at 72 h after pollination.
Fig. 4.

Hierarchical clustering analysis of the 73 differentially expressed proteins. Protein expression levels were presented as the Log ratio relative to the 0 h unpollinated reference. Colours ranging from green to red represent protein expression from the highest level of down-regulation to the highest level of up-regulation, respectively. Black colour indicates no change compared to 0 h. Two main clusters were formed, representing up- and down-regulated proteins during senescence. The similarities of protein expression patterns represented by the distance of tree branches are shown on the left side. The spot sample numbers with putative protein identification are indicated on the right side of the heat map.

Fig. 5.

Changes in relative protein abundance during corolla development and senescence. Protein expression patterns of each of the 14 subclusters from the hierarchical tree in Fig. 4 are shown by one representative protein. Spots that were up- or down-regulated are indicated with arrows on the 2-D gels. The solid and dotted lines indicate the protein expression changes in pollinated and unpollinated corollas from 0 h to 72 h, respectively.

Hierarchical clustering analysis of the 73 differentially expressed proteins. Protein expression levels were presented as the Log ratio relative to the 0 h unpollinated reference. Colours ranging from green to red represent protein expression from the highest level of down-regulation to the highest level of up-regulation, respectively. Black colour indicates no change compared to 0 h. Two main clusters were formed, representing up- and down-regulated proteins during senescence. The similarities of protein expression patterns represented by the distance of tree branches are shown on the left side. The spot sample numbers with putative protein identification are indicated on the right side of the heat map. Changes in relative protein abundance during corolla development and senescence. Protein expression patterns of each of the 14 subclusters from the hierarchical tree in Fig. 4 are shown by one representative protein. Spots that were up- or down-regulated are indicated with arrows on the 2-D gels. The solid and dotted lines indicate the protein expression changes in pollinated and unpollinated corollas from 0 h to 72 h, respectively.

Protein isoforms were identified from multiple spots

Forty-four per cent of the proteins in Tables 2 and 3 were identified in multiple spots with different pI and/ or Mr. Thus almost half of the differentially expressed proteins corresponded to various protein isoforms. Protein isoforms may be produced from different members of multigenic protein families or they may result from alternative splicing or post-translational modifications. In most instances, the multiple spots representing a protein function had the same pattern of either up- or down-regulation during the progression of corolla senescence. The proteins manganese superoxide dismutase (14-23 and 14-52), 1,4-benzoquinone reductase-like protein (14-22 and 36-17), lipoxygenase (14-6 and 50-2), Kunitz-type protease inhibitor family protein (36-19 and 36-20), vacuolar invertase (14-29, 36-21, and 36-22) and beta-xylosidase 2 (PhXYL2; 44-7 and 14-13) were all up-regulated during senescence (Table 2). The spots identified as methionine synthase (14-31, 44-1, 44-2, and 44-3), caffeoyl CoA 3-O-methyltransferase (44-21, 44-22, and 44-23), and S-adenosylmethionine synthetase (SAMS) (14-36, 14-37, 44-13, and 44-14) were all down-regulated during corolla senescence. By contrast, S-adenosyl-L-methionine:benzoic acid/salicylic acid carboxyl methyltransferase (PhBSMT2; 14-40 and 50-11), actin depolymerizing factor 1 (PhADF1; 49-21 and 50-10), actin depolymerizing factor 2 (PhADF2; 49-22 and 51-14) and photosystem II protein (14-25 and 49-20) had isoforms with opposite expression patterns (Tables 2, 3; Fig. 6).
Fig. 6.

Relative protein abundance of selected proteins and the corresponding gene expression in petunia corollas after pollination. Protein abundance of all isoforms corresponding to the protein of interest are presented as the spot volume from the replicated 2-DE gels (n=3) and gene expression was quantified using real-time PCR (n=3) as described in the methods. Error bars represent ±SD.

Functional classification of petunia proteins that were up-regulated during pollination-induced corolla senescence Sample number indicates the spot sample for in-gel digestion in the sequencing reports. These sample numbers (without dash) were also used to indicate spots on the 2-D gels in Fig. 2. Proteins were classified based on the biological processes according to homology to genes in the TAIR Gene Ontology database. Matched peptides indicate total number of peptides that matched to other proteins. The first and the second number in the parenthesis indicate the peptides with exact and homologous matches, respectively. Accession/EST number indicates the sequence in the searched database identified from either NCBI non-redundant protein database or SGN petunia EST database (http://www.sgn.cornell.edu/index.pl). Obser. Mr/pI, observed Mr/pI for each protein was calculated from the 2-D gels with Image PDQuest 7.4.0 software according to standard marker proteins. Theor. Mr/pI, theoretical Mr/pI of the matched proteins. Peptide coverage refers to the percentage of sequence coverage of the matched protein; no available (na) data when EST database was used. Mean of intensities is the average of three replicated spot normalized intensities at 48 U, 48 P, 72 U, and 72 P. U is hours after flowers are fully opened without pollination (unpollinated). P is hours after flowers are fully opened with pollination. ud, spot intensity was undetectable on the gel. Functional classification of petunia proteins that were down-regulated during pollination-induced corolla senescence Sample number indicates the spot sample for in-gel digestion in the sequencing reports. These sample numbers (without dash) were also used to indicate spots on the 2-D gels in Fig. 2. Proteins were classified based on the biological processes according to homology to genes in the TAIR Gene Ontology database. Matched Peptides indicate total number of peptides that matched to other proteins. The first and the second number in the parenthesis indicate the peptides with exact and homologous matches respectively. Accession/EST number indicates the sequence in the searched database identified from either NCBI non-redundant protein database or SGN petunia EST database (http://www.sgn.cornell.edu/index.pl). Obser. Mr/pI, observed Mr/pI for each protein was calculated from the 2-D gels with Image PDQuest 7.4.0 software according to standard marker proteins. Theor. Mr/pI, theoretical Mr/pI of the matched proteins. Peptide coverage refers to the percentage of sequence coverage of the matched protein; no available (na) data when EST database was used. Mean of intensities is the average of three replicated spot normalized intensities at 48 U, 48 P, 72 U, and 72 P. U is hours after flowers are fully opened without pollination (unpollinated). P is hours after flowers are fully opened with pollination. ud, spot intensity was undetectable on the gel. Relative protein abundance of selected proteins and the corresponding gene expression in petunia corollas after pollination. Protein abundance of all isoforms corresponding to the protein of interest are presented as the spot volume from the replicated 2-DE gels (n=3) and gene expression was quantified using real-time PCR (n=3) as described in the methods. Error bars represent ±SD.

Transcript abundance does not always correspond with protein expression

Genes encoding 10 selected proteins were cloned by RT-PCR. These included the senescence-specific proteins, auxin-response-like protein (PhAUX; 50-4), UTP-glucose-1-phosphate uridyltransferase (PhUGP; 14-9); abscisic stress ripening protein (PhASR4; 50-7), and beta-ureidopropionase (PhBUP; 36-33); and the protein isoforms of beta-xylosidase 1 (PhXYL1; 36-3), beta-xylosidase 2 (PhXYL2; 44-7, 14-13), lipoxygenase (PhLOX; 50-2, 14-6), vacuolar invertase (PhINV; 14-29, 36-22, 36-21), actin depolymerizing factor 1 (PhADF1; 49-21, 50-10), and actin depolymerizing factor 2 (PhADF2; 49-22, 51-14). Of those proteins that were up-regulated during corolla senescence, only 60% of their corresponding genes had similar senescence-related increases in transcript abundance (Fig. 6). Gene expression patterns for PhXYL2 (encoding proteins 44-7 and 14-13) were opposite to that observed for protein abundance, with higher transcript levels in unpollinated corollas. A large decrease in mRNA levels was detected by 24 h after pollination with little subsequent change from 48 h to 72 h. Three proteins were identified as vacuolar invertase (14-29, 36-21, and 36-22). The five peptides identified in spot 14-29 were identical to those identified in spot 36-22, suggesting that they may represent a single protein that is encoded by the same gene (PhINV). The accumulation of PhINV transcripts also had an opposite pattern to that of the protein spots. Transcript levels were highest at 0 h and decreased after pollination. Four proteins were identified as an actin depolymerizing factor. Two of these proteins corresponded to the petunia actin depolymerizing factor PhADF1 (Genbank no. AF183903) and two corresponded to petunia actin depolymerizing factor PhADF2 (Genbank no. AF183904) (Mun ). The PhADF1 spot 50-10 was not detectable in unpollinated corollas (0 h to 72 h) or at 24 h after pollination, but was newly detected in senescing corollas at 48 h and 72 h after pollination (Fig. 6). By contrast, the PhADF1 spot 49-21 was detected at basal levels in non-senescing corollas (0 to 72 U and 24 P) but was undetectable in 48 P and 72 P corollas. Similar patterns were observed for the PhADF2 spots (49-22 and 51-14). PhADF1 and PhADF2 transcripts were neither up- nor down-regulated by pollination, and showed very little change from 0 h to 72 h.

Discussion

Two-dimensional gel electrophoresis followed by mass spectroscopy was used successfully to identify proteins that were differentially regulated during petunia corolla senescence. While large-scale gene expression changes can be detected early in the senescence programme (24 h; ML Jones, unpublished results), protein differences were not detected until the mid-senescence stage (48 h) when corollas were starting to show visible symptoms of wilting. Proteomic techniques have been used in recent publications to investigate the post-transcriptional and post-translational regulation of senescence in leaves, but this is the first proteomic analysis of petal senescence. While a recent proteomic analysis of leaf senescence in Arabidopsis indicated that more than 50% of the soluble proteins that were visualized by 2-DE declined in abundance during the progression of senescence (Carp and Gepstein, 2007), our analysis of petals indicated that the largest percentage of differentially regulated proteins increased in abundance during senescence. At the mid-senescence stage (48 h after pollination) 64% of the differentially regulated proteins increased and at the late senescence stage (72 h after pollination) 66% of the differentially expressed proteins were up-regulated. Genomic analyses of petal senescence have also reported a greater number of senescence up-regulated compared to down-regulated genes (Breeze ; Hoeberichts ). Functional classification of the differentially expressed proteins revealed that the largest number of up-regulated proteins fell into the single biological processes category of stress and defence response genes. Genomic studies of leaves and petals have also reported that a large percentage of the senescence up-regulated genes encode proteins putatively involved in plant–pathogen interactions and response to stress (primarily reactive oxygen species metabolism) (Bhalerao ; van Doorn ; Breeze ; Buchanan-Wollaston ; Hoeberichts ; Price ). It has been hypothesized that the pathogenesis-related and stress-response proteins may protect these tissues from pathogen attack or the accumulation of damaging reactive oxygen species (ROS) during senescence to allow for maximum degradation of cellular constituents and remobilization of nutrients (Buchanan-Wollaston ). Some defence-related proteins were also identified as being down-regulated by pollination and may function during early petal development. Similarly, a large number of stress-related proteins were identified from rose petals in a proteomic analysis of early rose flower development from the tight bud to the fully open flower stages (Dafny-Yelin ). While it is well known that ethylene is the primary hormone regulating petunia flower senescence, a number of the senescence up-regulated proteins identified in this study were putative orthologues of stress response genes that are regulated by other hormones including abscisic acid (ABA) and auxin. Protein 50-4 was identified as a putative orthologue of Arabidopsis DFL1 (Dwarf in Light), a member of the auxin-responsive GH3 gene family. The IAA/AUX, SAUR, and GH3 gene family members are auxin primary response genes and DFL1 functions as an indole-3-acetic acid amido synthetase involved in auxin signal transduction (Nakazawa ). A GH3 gene (CcGH3) in pepper is highly expressed during the late stage of fruit ripening rather than the early stage when auxin levels are highest. Promoter analysis shows that CcGH3 has both auxin and ethylene responsive elements and that CcGH3 expression is mainly responsive to ethylene during fruit ripening (Liu ). These results suggest that CcGH3 may be regulated by ethylene or that signal crosstalk between auxin and ethylene during fruit development and ripening is involved in its regulation (Liu ). While auxin treatment can accelerate senescence in some flowers, this response is mediated by ethylene. It is therefore possible that the auxin response-like protein (50-4) identified in our study plays a role in mediating the ethylene-induced senescence of pollinated flowers. A senescence-specific petunia protein (50-7) was identified as having homology to tomato abscisic stress ripening protein 4 (LeASR4). The ASR gene family includes ASR1, ASR2, ASR3, and ASR4. ASR1, ASR3, and ASR4 are induced by dehydration, ABA, and cold stresses, whereas ASR2 is ABA-independent and is specifically induced by dehydration (Doczi ). Overexpression of tomato ASR1 or lily ASR1 in transgenic tobacco or Arabidopsis enhances tolerance to drought and salt stresses (Goldgur ). Expression analysis using promoter::GUS fusion lines show that DS2 not only responds to dehydration, but it is also expressed in pollen, flowers, and fruits, suggesting that it is regulated by both developmental cues and stress stimuli (Doczi , 2005). The grape ASR1 orthologue and Solanum DS2 have been hypothesized to function as transcription factors (Silhavy ; Cakir ). The PhASR4 protein that was identified by 2-DE was detected only in senescing corollas at 48 h and 72 h after pollination, while PhASR4 transcript abundance was greatly increased at 24 h after pollination and decreased at 48 h and 72 h. To our knowledge, ASR genes have not previously been identified in senescing petals. This early senescence induction of PhASR4 gene expression supports a potential role in senescence signalling. The senescence of petunia corollas is accompanied by a decrease in fresh weight and an increase in ABA content (Chang ). Changes in the relative abundance of the ASR/DS2-type protein identified in our experiment may be the result of increased petal ABA levels or petal dehydration. Whether this protein plays a role in regulating the senescence response or is a side-effect of senescence is unknown. Up-regulated proteins in the biological classes of carbohydrate, lipid, nucleotide, amino acid, and nitrogen metabolism primarily included enzymes involved in catabolic processes that putatively function in the large-scale degradation of macromolecules and organelles that occurs during senescence to allow for nutrient remobilization. Proteins within these categories could collectively be referred to as remobilization proteins and make up 46% of the up-regulated proteins identified in this study. Most of the changes in gene expression identified using microarray and differential display techniques are also related to genes involved in remobilization, in keeping with the primary function of the senescence programme (Jones, 2004; Stead ; van Doorn and Woltering, 2008). One major difference that can be seen by comparing the genes that were up-regulated in a recent study of wallflower petal senescence is the absence of transporter proteins in our proteomic study (Price ). This is probably due to our protein extraction procedure, which does not adequately isolate membrane bound proteins. The proteomic techniques employed in this work are most useful for identifying highly abundant soluble proteins involved in the execution phase of senescence, while complementary genomic approaches are required to identify regulatory proteins involved in the initiation phase of senescence. Proteins within the functional categories of amino acid, carbohydrate, and nitrogen metabolism were also observed for the down-regulated protein spots. Many of these down-regulated proteins may be anabolic enzymes that catalyse synthesis reactions that cease with the onset of senescence. It should also be noted that some of the proteins that were detected only at 0 h and 24 h after pollination may represent proteins that have changed location on the gel due to a post-translational modification that resulted in a change in the Mr and/or pI of the protein. The next largest single functional category of up-regulated proteins after defence and stress response was carbohydrate metabolism. Microarray analyses of petal senescence in Alstroemeria and Dianthus caryophyllus (carnation) also revealed that a larger number of up-regulated genes were putatively involved in carbohydrate metabolism than were involved in lipid, nucleic acid or protein metabolism (Breeze ; Hoeberichts ). Many of the petunia proteins identified within this category have a putative function in cell wall degradation, while down-regulated proteins were also identified that might be involved in the synthesis of cell wall components. Cell wall degradation is a component of senescence in some flower species. Sandersonia flowers show very little cell wall degradation, while considerable degradation of the cell walls has been reported during the senescence of Hemerocallis (daylily), Ipomoea nil (morning glory), Iris, carnation, and petunia petals (Winkenbach, 1970; de Vetten ; Panavas ; O'Donoghue , 2009; van Doorn ). A recent study of flower development in Petunia×hybrida ‘Mitchell Diploid’ reported that the metabolism of cell wall-associated polymeric galactose was a major feature of flower opening and senescence (O'Donoghue ). Our 2-D profiling identified a β-galactosidase protein (14-11) that was highly abundant in non-senescing corollas and which increased in abundance at 48 h and then again at 72 h after pollination. The metabolism of cell wall galactose catalysed by β-galactosidases plays a role in cell wall expansion during growth and in the loss of cell wall adhesion during senescence. β-galactosidase enzyme activity increases during flower bud opening and continues to increase during the progression of petal senescence (O'Donoghue ). β-galactosidases are encoded by a multigene family and two β-galactosidase genes were recently cloned from petunia. PhGal1 is expressed at high levels in flower buds before flower opening and PhGAL2 is expressed at a similar level in open flowers and during petal senescence (O'Donoghue ). In our study, the abundance of the putative β-galactosidase protein 14-11 was highly up-regulated during petal senescence. The peptide sequences from 14-11 do not match the predicted amino acid sequences of either PhGAL1 or PhGAL2, and 14-11 seems to represent a third member of the petunia β-galactosidase gene family with a role in cell wall galactose mobilization during petal senescence. Beta-xylosidases are also involved in cell wall polysaccharide disassembly or modification and may function in cell wall degradation (Itai ; Hayama ). Multiple differentially regulated spots were identified as beta-xylosidases (36-3, 44-7, and 14-13). Proteins 44-7 and 14-13 were orthologues of the tomato beta-D-xylosidase LeXYL2, while 14-13 was most homologous to the Arabidopsis beta-D-xylosidase AtBXL1(Goujon ; Itai ). Studies in tomato have shown that the different members of the beta-D-xylosidase family are differentially regulated during fruit development. LeXYL1 mRNA levels increase at the late fruit ripening stage, while LeXYL2 mRNA levels are higher in premature fruit and decrease at the late ripening stage (Itai ). The increased transcript abundance of PhXYL2 in unpollinated flowers at 0 h suggests that it may play a primary role in cell wall loosening and expansion during flowering opening and that the up-regulation of the C-terminally truncated form of the PhXYL2 protein (14-13) during senescence may result from targeted degradation of the protein after pollination. By contrast, PhXYL1 may catalyse cell wall degradation specifically during petunia corolla senescence. A decrease in the total protein, RNA, and DNA content of senescing petals has been reported in many flowers including petunia (Winkenbach, 1970; Xu and Hanson, 2000; Jones, 2004; Jones ; Langston ). While this change could be the result of decreased protein and nucleic acid synthesis, large increases in protease and nuclease activity suggest that large-scale degradation is occurring (Stephenson and Rubinstein, 1998; Xu and Hanson, 2000; Jones ; Langston ; Pak and van Doorn, 2005). While the integrity of cellular membranes must be maintained during senescence to allow for maximum cellular recycling, the later stages of petal senescence are accompanied by increased activity of enzymes involved in phospholipid and fatty acid catabolism (Borochov ). In support of this role in nutrient remobilization, a number of senescence up-regulated proteins were identified with a putative role in protein, nucleic acid or lipid degradation. These included putative proteases (14-50, 36-13, and 14-27), endonucleases (14-12), lipoxygenases (50-2 and 14-6), and lipases (14-20). These types of hydrolytic enzymes have previously been reported in petunia (Borochov ; Xu and Hanson, 2000; Jones, 2004; Jones ; Langston ) and their putative roles in senescence have been discussed in a number of recent reviews (Hopkins ; Reid and Chen, 2007; van Doorn and Woltering, 2008; Jones ). While senescence up-regulated proteases may be involved in large-scale protein degradation, proteins involved in the ubiquitin–proteasome pathway appear to be involved in selective protein degradation during senescence (Buchanan-Wollaston ; Stead ; van Doorn and Woltering, 2008). In plants, ubiquitin- and proteasome-dependent proteolysis is involved in the regulation of transcription factors and ER-associated protein degradation during development and in response to stress. It also directly impacts ethylene signalling during senescence by regulating the degradation of the EIN3/EIL transcription factors and the ethylene receptors (ETR) (Guo and Ecker, 2003; Kevany ). Increased expression of genes putatively involved in protein ubiquitination and those encoding protein subunits of the proteasome have been identified in senescing Alstroemeria, Narcissus pseudonarcissus, carnation, morning glory, Mirabilis jalapa (four o'clocks), and petunia petals (Hunter ; Breeze ; Hoeberichts ; Xu a, b; Yamada ). A recent proteomic analysis of leaves also identified a proteasome beta subunit A1 that was up-regulated during nitrate starvation-induced senescence (Desclos ). Our proteomic analysis identified an ubiquitin-conjugating enzyme family protein (E2) (36-24) and a delta-subunit of the 20S proteasome (36-15) that were up-regulated by pollination. MS analysis revealed that a number of the up-regulated spots contained truncated proteins that had a much lower than predicted Mr. The appearance of these truncated proteins in senescing petals may be explained by the selective degradation of proteins that are no longer needed during the execution phase of senescence. The mass spectra for the peptides in band 14-29 (vacuolar invertase) support the theory that proteins are ubiquitinated during senescence. A peptide with a molecular mass of 1744 Da was detected that contained the sequence 231YSGNPVMoVPPPGIGVK246 plus an additional 114 Da, and the same peptide in its unmodified form (1630 Da) was also detected in band 14-29 (data not shown). The additional 114 Da occurred on the C-terminal fragment of the GVK residue. This mass increase can be explained by ubiquitination of the K residue, as the addition of two glycines to the K would account for the 114 Da increase. Down-regulating the expression of a RING domain E3 protein delays visual symptoms of petal wilting in petunia by 2 d, suggesting that the ubiquitin–proteasome pathway does play a regulatory role in senescence (Xu a). Actin-depolymerizing factor (ADF) is one of the small actin-binding proteins that modulates actin cytoskeleton dynamics in eukaryotic cells (Mun ; Thomas ). In our 2-DE experiments, two proteins encoded by PhADF1 (50-10 and 49-21) and two proteins encoded by PhADF2 (50-14 and 49-22) were identified. Previous studies have shown that the petunia PhADF1 and PhADF2 genes are expressed in all tissues except pollen, with the highest transcript abundance in petals (Mun ). In our proteomic analysis the two PhADF1 isoforms and the two PhADF2 isoforms had opposite expression patterns. MS analysis indicated that the down-regulated PhADF1 and PhADF2 isoforms were acetylated at the second amino acid of the N-terminus. We do not know if the two up-regulated proteins were also acetylated, as that same peptide was not sequenced. The down-regulated PhADFs had a higher observed Mr than the PhADF proteins that were up-regulated, supporting the theory that they may represent different stages in the post-translational modification of the ADF proteins. Although there is no evidence to indicate that acetylation regulates ADF function, the activity of the ADF protein and its interaction with actin is known to be regulated by pH, phosphorylation, and phosphoinositide (Lopez ). The stabilization of actin depolymerization has been shown to induce PCD in plants (Thomas ), and the post-transcriptional regulation of PhADF1 and PhADF2 may trigger or modulate corolla senescence. While most studies of senescence have focused on gene expression changes, senescence is also controlled by post-transcriptional and post-translational regulation (Thomas ). Treating flowers with the protein synthesis inhibitor, cycloheximide, significantly delays corolla senescence and provides evidence for the role of newly translated proteins in the senescence program (Wulster ; Lay Yee ; van Doorn ). In this study, ten genes corresponding to senescence up-regulated proteins were cloned to determine patterns of gene expression during petal senescence. Only six of the ten genes (60%) had senescence up-regulated increases in transcript abundance. Genes encoding PhXYL2 and PhINV showed an opposite expression pattern, with transcript abundance decreasing during senescence. Even more interestingly, the genes encoding PhADF1 and PhADF2 showed little change in relative mRNA abundance during senescence, even though protein isoforms of ADF1 and ADF2 were identified that were up- and down-regulated during petal senescence. In Arabidopsis, where extensive microarray data are available to complement proteomic analyses, a poor correlation between protein and gene expression has been reported (Jiang ). Similar to what was observed in petals, Carp and Gepstein (2007) reported that only approximately 60% of the up-regulated proteins identified during leaf senescence had corresponding increases in mRNA abundance.

Conclusion

In this study, a proteomic analysis of petunia petals identified 46 proteins that were up-regulated and 27 proteins that were down-regulated during senescence. The proteins identified were implicated in a wide range of biological processes with the majority having a putative role in stress and defence responses or catabolism and remobilization. This is the first proteomic analysis of petal senescence and as such it provides new insights into the pathways that execute senescence and the post-transcriptional regulation of senescence. While our MS analysis did not allow for the precise identification of specific protein modifications, 37% of the senescence up-regulated proteins had higher than predicted Mr and changes in the pI that indicate potential post-translational modifications. This study provides a good starting point for further genetic analysis of the functional role of these proteins in petal senescence.

Supplementary data

Supplementary data are available at JXB online. lists the primers and Tm used for real-time RT-PCR analysis. is a flow chart outlining the methodology for mass spectroscopy and protein identification.
  54 in total

1.  Plant responses to ethylene gas are mediated by SCF(EBF1/EBF2)-dependent proteolysis of EIN3 transcription factor.

Authors:  Hongwei Guo; Joseph R Ecker
Journal:  Cell       Date:  2003-12-12       Impact factor: 41.582

2.  Increases in DNA fragmentation and induction of a senescence-specific nuclease are delayed during corolla senescence in ethylene-insensitive (etr1-1) transgenic petunias.

Authors:  Brennick J Langston; Shuangyi Bai; Michelle L Jones
Journal:  J Exp Bot       Date:  2004-10-08       Impact factor: 6.992

3.  Comparative proteomic analysis of NaCl stress-responsive proteins in Arabidopsis roots.

Authors:  Yuanqing Jiang; Bo Yang; Neil S Harris; Michael K Deyholos
Journal:  J Exp Bot       Date:  2007-10-04       Impact factor: 6.992

Review 4.  Plant proteome analysis: a 2006 update.

Authors:  Jesús V Jorrín; Ana M Maldonado; Ma Angeles Castillejo
Journal:  Proteomics       Date:  2007-08       Impact factor: 3.984

5.  Endoglycanase-Catalyzed Degradation of Hemicelluloses during Development of Carnation (Dianthus caryophyllus L.) Petals.

Authors:  N C de Vetten; D J Huber; K C Gross
Journal:  Plant Physiol       Date:  1991-03       Impact factor: 8.340

6.  Two closely related cDNAs encoding actin-depolymerizing factors of petunia are mainly expressed in vegetative tissues.

Authors:  J H Mun; H J Yu; H S Lee; Y M Kwon; J S Lee; I Lee; S G Kim
Journal:  Gene       Date:  2000-10-31       Impact factor: 3.688

7.  Programmed cell death during pollination-induced petal senescence in petunia.

Authors:  Y Xu; M R Hanson
Journal:  Plant Physiol       Date:  2000-04       Impact factor: 8.340

8.  Gene expression during anthesis and senescence in Iris flowers.

Authors:  W G van Doorn; P A Balk; A M van Houwelingen; F A Hoeberichts; R D Hall; O Vorst; C van der Schoot; M F van Wordragen
Journal:  Plant Mol Biol       Date:  2003-12       Impact factor: 4.076

9.  Proteomic analysis of resistance mediated by Rcm 2.0 and Rcm 5.1, two loci controlling resistance to bacterial canker of tomato.

Authors:  Gitta L Coaker; Belinda Willard; Michael Kinter; Eric J Stockinger; David M Francis
Journal:  Mol Plant Microbe Interact       Date:  2004-09       Impact factor: 4.171

10.  Cluster analysis and display of genome-wide expression patterns.

Authors:  M B Eisen; P T Spellman; P O Brown; D Botstein
Journal:  Proc Natl Acad Sci U S A       Date:  1998-12-08       Impact factor: 11.205
View more
  16 in total

Review 1.  Using proteomics to study sexual reproduction in angiosperms.

Authors:  Ján A Miernyk; Anna Preťová; Adela Olmedilla; Katarína Klubicová; Bohuš Obert; Martin Hajduch
Journal:  Sex Plant Reprod       Date:  2010-09-10

Review 2.  Organisation and regulation of the cytoskeleton in plant programmed cell death.

Authors:  A Smertenko; V E Franklin-Tong
Journal:  Cell Death Differ       Date:  2011-05-13       Impact factor: 15.828

3.  Labellum transcriptome reveals alkene biosynthetic genes involved in orchid sexual deception and pollination-induced senescence.

Authors:  Filipa Monteiro; Mónica Sebastiana; Andreia Figueiredo; Lisete Sousa; Helena C Cotrim; Maria Salomé Pais
Journal:  Funct Integr Genomics       Date:  2012-06-17       Impact factor: 3.410

4.  Involvement of three annexin genes in the ripening of strawberry fruit regulated by phytohormone and calcium signal transduction.

Authors:  Jingxin Chen; Linchun Mao; Hongbo Mi; Wenjing Lu; Tiejin Ying; Zisheng Luo
Journal:  Plant Cell Rep       Date:  2016-01-02       Impact factor: 4.570

5.  Proteomes and Ubiquitylomes Analysis Reveals the Involvement of Ubiquitination in Protein Degradation in Petunias.

Authors:  Jianhang Guo; Juanxu Liu; Qian Wei; Rongmin Wang; Weiyuan Yang; Yueyue Ma; Guoju Chen; Yixun Yu
Journal:  Plant Physiol       Date:  2016-11-03       Impact factor: 8.340

6.  Interaction of plant growth regulators and reactive oxygen species to regulate petal senescence in wallflowers (Erysimum linifolium).

Authors:  Faezah Mohd Salleh; Lorenzo Mariotti; Natasha D Spadafora; Anna M Price; Piero Picciarelli; Carol Wagstaff; Lara Lombardi; Hilary Rogers
Journal:  BMC Plant Biol       Date:  2016-04-02       Impact factor: 4.215

7.  The scutellum of germinated wheat grains undergoes programmed cell death: identification of an acidic nuclease involved in nucleus dismantling.

Authors:  Fernando Domínguez; Javier Moreno; Francisco Javier Cejudo
Journal:  J Exp Bot       Date:  2012-08-09       Impact factor: 6.992

8.  Mineral nutrient remobilization during corolla senescence in ethylene-sensitive and -insensitive flowers.

Authors:  Michelle L Jones
Journal:  AoB Plants       Date:  2013-05-09       Impact factor: 3.276

9.  Comparisons of protein profiles of beech bark disease resistant and susceptible American beech (Fagus grandifolia).

Authors:  Mary E Mason; Jennifer L Koch; Marek Krasowski; Judy Loo
Journal:  Proteome Sci       Date:  2013-01-14       Impact factor: 2.480

10.  Annexin5 plays a vital role in Arabidopsis pollen development via Ca2+-dependent membrane trafficking.

Authors:  Jingen Zhu; Xiaorong Wu; Shunjie Yuan; Dong Qian; Qiong Nan; Lizhe An; Yun Xiang
Journal:  PLoS One       Date:  2014-07-14       Impact factor: 3.240
View more

北京卡尤迪生物科技股份有限公司 © 2022-2023.