Proteomic and Metabolic Pro
- 格式:pdf
- 大小:517.85 KB
- 文档页数:12
Proteomic and Metabolic Profiling of Rice Suspension Culture Cells as a Model to Study Abscisic Acid Signaling Response Pathways inPlantsSushma R.Rao,†Kristina L.Ford,Andrew M.Cassin,Ute Roessner,John H.Patterson,‡andAntony Bacic*Australian Centre for Plant Functional Genomics,School of Botany,University of Melbourne,VIC 3010,AustraliaReceived July 28,2010Rice (Oryza sativa cv Taipei 309)suspension culture cells (SCCs)were used as a simple,single cell model system to gain insights into the complex abscisic acid (ABA)signaling response pathways in plants.Following system establishment involving morphological observations and transcript profiling of genes known to be ABA responsive in planta,a comprehensive proteomic and metabolomic study was performed.A total of 759buffer-soluble proteins that included 3284peptides categorized into 656protein families are ing iTRAQ,only 36of these proteins showed statistically significant changes in abundance in response to ABA.In addition,a GC -MS based metabolite profiling study allowed the identification of 148metabolites that included 25amino acids (AAs),45organic acids (OAs),35sugars,19fatty acids,2polyamines,4sterols,5sugar acids,4sugar alcohols,and 9miscellaneous compounds.Of these,only 11(8.8%)changed in a statistically significant manner in response to ABA treatment.These studies provide important insights into plant responses to ABA at the protein and metabolite level.Keywords:Oryza sativa (rice)suspension culture cells •ABA •proteomics •metabolomicsIntroductionThe plant hormone abscisic acid (ABA)is a mediator of cellular responses during growth and development as well as to environmental stresses,including drought.ABA is produced in the roots in response to a water deficit and this signal moves throughout the plant where it is perceived by ABA receptors.This interaction activates a complex web of signaling pathways that include protein phosphorylation,activation of G-proteins and modulation of RNA metabolism.1,2The pH and redox status of the cell are crucial factors in mediating ABA signal transduction.3These signals lead to appropriate physiological,metabolic,transcriptional and translational modifications.Studying ABA signal transduction at the molecular level is helping to identify new targets and pathways involved in these responses.Part of the rationale for the present study was to evaluate the utility of SCCs,given they are a simple single cell type that can be generated in high yield,as a model system to study plant ABA responses using high throughput ’omics technologies.Advances in proteomic and metabolomic techniques have allowed the large-scale profiling of samples from many species.The ability to analyze and compare dynamic changes quanti-tatively in a plant proteome in response to a stimulus,using the iTRAQ approach is a significant added advantage of modern proteomics.4Proteome profiling has been carried out in various plant species such as Arabidopsis mitochondrion,5Pinus ra-diata ,6Holm oak leaf,7wheat lemma 8and leaf.9Comprehensive comparative studies of protein profiling during either develop-ment or stress responses,and functional analysis and charac-terization of regulatory processes are now needed to under-stand plant complex signaling and response mechanisms.Likewise,with the detection and relative quantification of many of the major metabolites,it should be possible to link them to biochemical pathways and understand their role.It is predicted that there are around 200000metabolites in the plant king-dom.10Metabolite profiling has been used to study tempera-ture,water,salinity,sulfur,phosphorus,oxidative and heavy metal stress responses in a number of plant species 11-13and have resulted in the identification of previously unknown metabolic responses and networks.In this study,GC -MS based metabolite profiling was used,as its utility has been demon-strated in the analysis of metabolites from Solanum tubero-sum ,14Lycopersicon esculentum ,15Hordeum vulgare ,12Lotus japonicus ,16Oryza sativa 17,18and Arabidopsis thaliana.19The integrated metabolomic and quantitative proteomic approach adopted in this study reveals some of the known components of the ABA signaling process and also suggests the possible involvement of other factors.The time-course design enabled*To whom correspondence should be addressed.Antony Bacic.E-mail:abacic@.au.Phone.:+613-83445041.Fax:+913-93471071.Web-site:.au.†Current address:Sushma Rao,Children’s Medical Research Institute,214Hawkesbury Road,Westmead,NSW 2145.‡Current address:John Patterson,Viewbank College,Warren Rd,View-bank VIC 3084.10.1021/pr100788m2010American Chemical SocietyJournal of Proteome Research 2010,9,6623–66346623Published on Web10/20/2010the temporal tracking of changes together with identification of relevant proteins and metabolites in response to ABA. Materials and MethodsPlant Material and ABA Treatment.All chemicals were from Sigma Aldrich unless otherwise specified.Rice(Oryza sativa cv Taipei309)suspension cell cultures(SCCs),a gift from the late Prof.Hans Kende,Michigan State University,East Lansing, MI,were grown in Murashige and Skoog medium lacking IAA, kinetin or sucrose(MP Biomedicals).This was supplemented with6%(w/v)sucrose,60mM2,4-dichlorophenoxyacetic acid (2,4-D)and 2.5mM2-(N-morpholino)ethanesulfonic acid (MES).The pH of the medium was adjusted to5.6with1M KOH.The cultures were maintained in250mLflasks on shakers running at100rpm at a temperature of23°C in the dark.The cells were subcultured every7days by transferring50mL of cells and media into fresh MS medium.Rice SCCs(1L)were grown to their midlog phase(7d after subculture).Cell growth curves were obtained by monitoring the fresh weight of cells for2weeks in the presence or absence of100µM ABA(added from a100mM stock,dissolved in ethanol).At48h intervals, the fresh weight of cells was recorded for control cells and cells from cultures supplemented with either100µL ethanol or100µL MS media(3replicates in separateflasks used for each treatment).In each experiment,the cells were harvested and weighed in triplicate(Supporting Information I).Cells in the exponential phase of their growth(7days)were used for all further experimentation.For proteomic and metabolomic experiments,the parental culture(1L)was divided into4separate cultureflasks(250mL each).Oneflask from each set was used as the untreated control.The3remainingflasks were treated with100µM ABA for0.5,2,and6h,respectively.Three individual biological replicate parentalflasks were grown and divided as described for each time point.The triplicates were harvested byfiltering the SCCs through a vacuumfilter setup with a cellulosefilter disc(Whatman)to remove the media.The cells were then snap frozen in liquid nitrogen.Biological replicates,three for pro-teomics andfive for metabolomics were used.Real Time-Q-PCR.Cells treated with100µM ABA for0 (without ABA),2,4,6,and8h were used for the Q-PCR experiments.Separateflasks were used for each time sampling point.The primer pairs indicated in Supplementary Table4 (Supporting Information)were used for RT-QPCR experiments carried out at the University of Adelaide with the assistance of Dr.Neil Shirley(ACPFG,.au/).The size and the identity of the products were confirmed by sequencing. Q-PCR analysis was carried out as per the method described by Burton et al.20Three replicate PCRs for each of the cDNAs were included in every run.The data obtained were normalized with respect to the control genes,actin(OsActin),ribosomal (10S)RNA(Os10SRNA),tubulin(OsTubulin)and glyceralde-hyde-3-phosphate dehydrogenase(OsGAPDH)and analyzed.Protein Extraction.The frozen samples were ground in liquid nitrogen to afine powder and suspended in an extraction buffer containing Tris.HCl(50mM pH8.5),MgCl2(15mM), NaCl(75mM), -glycerophosphate(15mM),para-nitrophe-nylphosphate(pNPP)(15mM),NaF(1mM),EDTA(0.25M), PMSF(1mM),Na3VO4(0.5mM),DTT(1mM),aprotinin(40 g/L),leupeptin(20g/L)and pepstatin A(20g/L)at4°C. Samples were vortexed and centrifuged at3000×g for30min. The supernatant was collected and further centrifuged at 100000×g for30min.Proteins in the membrane-depleted supernatant were collected and concentrated by mixing with six volumes of-20°C equilibrated acetone,left at-20°C for 2.5h and centrifuged at3000×g for10min.The protein pellet was resuspended in triethylammonium bicarbonate(0.5M,pH 8.5)containing0.1%SDS.Soluble protein concentration was determined by the Bradford assay using bovine serum albumin (BSA)as a standard.21The acetone precipitated proteins(40µg)were reduced with tris-(2-carboxylethyl)phosphine(5mM), at60°C for1h.Cysteine residues were blocked by incubating proteins in methyl methane thiosulphate(MMTS)(90mM)in isopropanol for10min at RT.Proteins were then digested with trypsin(4µg;1:10w/w,sequencing grade,Promega)for16h at37°C.iTRAQ Labeling.iTRAQ tags(Applied Biosystems)were resuspended in ethanol(70µL)and added to the peptide mixtures.One set of iTRAQ tags were used per biological replicate(replicate1:control114;0.5h115;2h116;6h117, replicate2:control115;0.5h116;2h117;6h114,replicate3: control116;0.5h117;2h114;6h115).iTRAQ tagging was carried out by incubation at RT for1h and resulting labeled peptide mixtures were pooled and lyophilized for2D-LC-MS/ MS.Offline Strong Cation Exchange of the iTRAQ-Labeled Peptides.The iTRAQ tagged peptide sample was resuspended in0.1%formic acid and passed over a C18SepPak column (Waters,Milford,MA)and concentrated under vacuum to a final volume of100µL.The concentrated tryptic peptides were separated on a PolySULFOETHYL Aspartamide SCX column (4.6mm×200mm,5µm,300Å,PolyLC Inc.,Columbia,MD) attached to an Agilent1100series HPLC system(Agilent Technologies,Palo Alto,CA)with the following separation gradient:buffer A(25%(v/v)acetonitrile in5mM phosphate buffer,pH3)for10min,then up to100%buffer B(300mM potassium chloride,25%(v/v)acetonitrile in5mM phosphate buffer,pH5)over30min at aflow rate of0.7mL/min with0.5 min fractions being collected in a96-well plate.Inline Reversed-Phase LC-ESI-MS/MS.Fractions obtained from SCX-HPLC were reduced under vacuum and resuspended in0.1%formic acid(60µL),filtered through a minisart membrane(0.2µm;Sartorius Stadim Biotech,Aubagne,France) and one-quarter of each fraction was loaded onto a reversed-phase precolumn(300µm×5mm Zorbax300SB-C18;Agilent Technologies,Palo Alto,CA)attached to a Shimadzu Promi-nence nano LC system(Shimadzu Corporation,Kyoto,Japan). The precolumn was washed with0.1%formic acid in5% acetonitrile for15min before placing in-line with a75µm i.d.×150mm Zorbax300SB-C18(Agilent Technologies,Palo Alto, CA)reversed-phase column.Peptides were eluted using a gradient of5-65%(v/v)acetonitrile in0.1%formic acid over 60min,at aflow rate of0.25µL min-1.Peptides were analyzed via electrospray ionization(ESI)on a QSTAR Elite hybrid quadrupole time-of-flight mass spectrometer(Applied Biosys-tems/MDS Sciex,Foster City,CA).Each SCX-HPLC fraction was chromatographed and analyzed three times.The MS was operated in the positive ion mode,ion source voltage of2200V,using10µm uncoated SilicaTips(New Objective,Woburn,MA).Analyst QS2.0software(Applied Biosystems/MDS Sciex,Foster City,CA)was used to collect data in a data-dependent acquisition mode for the three most intense ions fulfilling the following criteria:m/z between450 and2000;ion intensity40counts;and charge state between +2and+5.After MS/MS analysis,these ions were dynamically excluded for18s,using a mass tolerance of50mDa.MS scansresearch articles Rao et al. 6624Journal of Proteome Research•Vol.9,No.12,2010were accumulated for0.5s,and MS/MS scans were collected in automatic accumulation mode for a maximum of2s.Mass and charge state-dependent rolling collision energy was used and the MS instrument was calibrated daily with[Glu]-fibrinopeptide B(Sigma-Aldrich,St.Louis,MO).Proteomic Data Analysis.Peak lists from the MS/MS spectra were made using ProteinPilot software version2.0.1(Applied Biosystems/MDS Sciex Foster City,CA).The peak list was searched against Oryza sativa proteins downloaded from NCBI in July2009using MASCOT22and the Paragon Algorithm (Applied Biosystems/MDS Sciex).23The rationale for using two separate algorithms was to reduce the false positive rates of the peptides identified.The false positive rate for this study was calculated using randomized version of the rice protein database from NCBI and was found to be0.2%.The MASCOT parameters were:enzyme:trypsin,fixed modifications:iTRAQ4plex(N-term);iTRAQ4plex(K);Meth-ylthiol(C),variable modifications:iTRAQ4plex(Y),MS peptide tolerance:0.25Da,MS/MS tolerance:0.15Da,number of missed cleavages:up to1.The Paragon Algorithm parameters were:sample type:iTRAQ4plex(peptide labeled);Cys Alkyla-tion:MMTS;Digestion:Trypsin;Search effort:Thorough ID. The outputs from both search algorithms were combined and only proteins with2or more peptides with a P<0.05in both search algorithms were reported.The reporter ion peak areas generated in ProteinPilot were used for quantification.Any peptide with a reporter ion peak area of less than20was removed from quantification.When a peptide was detected more than once,the peak area for each reporter ion was summed,each peptide was then normalized by the sum of its channel intensities(114,115,116,and117). Peptides were ignored when the normalized peptide value was more than2standard deviations from the calculated mean of the protein the peptide matched to.The mean was then calculated for proteins with3or more peptides that fulfilled the above criteria.The time points were normalized to the control and again by the average protein ratio for each time point.The results for each replicate were combined and only those proteins seen in all3replicates are reported.This analysis therefore focused primarily on proteins that showed changes in abundance in response to ABA,rather than proteins that may have become differentially post-translationally modified in response to ABA.Proteins were classified into their respective functional categories as specified either in NCBI or Uniprot (Universal protein resource:/).Extraction,Derivatization,Profiling and Statistical Analysis of Metabolites.The harvested rice SCCs were frozen in liquid nitrogen,homogenized to afine powder which was used for metabolite extraction using a modified method of Jacobs et al.18 To account for biological and technical variability,5biological replicates(SCCs grown in separateflasks)were analyzed for each particular time point.15Approximately70-100mg of the sample(accurate weights were recorded)was extracted in100% methanol(500µL)and a polar internal standard(20µL of0.2 mg/mL ribitol/norleucine in water)was added.The mixture was extracted for15min at70°C and water(500µL)was added and vortexed.The mixture was centrifuged at14000rpm for 10min.The dried polar residue was derivatized using both the trimethylsilyl(TMS)and tributyldimethylsilyl(TBS)methods.14,15 as described by Roessner et al.12and the samples(1µL)were injected onto a GC column in the splitless mode and run on a GC-MS system comprising of an AS3000auto sampler,an ultra trace GC and a DSQ quadrupole MS(Thermo Electron Cor-poration,Madison,WI)using the conditions described in Roessner et al.12The analysis of TBS samples was performed as described in Jacobs et al.18and the mass spectra were recorded at two scans per second over a mass range of mass-to-charge(m/z)ratio of70to600atomic mass units(amu). Both total ion chromatograms(TIC)and mass spectra were evaluated using the Xcalibur program(ThermoFinnigan)and the resulting data prepared and normalized as described in Roessner et al.14Mass spectra of eluting compounds were identified using a combination of an in-house constructed mass spectra library of authentic standards,the public domain mass spectra library of the Max-Planck-Institute for Molecular Plant Physiology,and the commercial mass spectra library of the National Institute of Standards and Technology.14,18,24,25Where possible,all matching mass spectra were additionally confirmed by determination of the retention time and mass spectra by analysis of authentic standards.Resultant relative response ratios were normalized per gram extracted fresh weight as previously described.12Metabolite profile data is presented as fold changes with respect to the control which is set at1. Differences between the treated samples were considered as significant when the p-value using student’s t test was<0.05. The Bonferroni correction was applied to reduce the number of false positives and only metabolites significant after this correction were considered as most important.26,27 ResultsTranscriptomic Profiling.To establish the validity of the SCC system as a model to examine ABA responses,the effect of exposure to100µM ABA on transcript abundance(see Figure 1)growth rate and morphology(Supporting Information I)was initially examined.The expression levels of four different genes known to be ABA responsive in planta,myosin heavy chain-like protein(OsMyo1),28chlorophyll a/b-binding protein pre-cursor(OsCab26),28low temperature induced protein19 (OsLIP19)28and putative sphingosine kinase(OsSPHK)28were analyzed by QPCR.Expression was normalized to the control genes,OsActin,Os10SRNA,OsTubulin and OsGAPDH(Figure 1).Three of the four genes showed an ABA-dependent increase in transcript abundance.The increase was transient and either at the2h(OsCab26;Figure1B)or the4h(OsMyo1;Figure1A, OsLIP19;Figure1C)time point after ABA exposure.OsSPHK did not respond to ABA treatment(Figure1D)and surprisingly decreased relative to the no-ABA control.Measurement of the growth of SCCs following ABA exposure also suggested an effect that presented as an initial lag in growth caused by a loss in turgor but following osmotic adjustment the cells entered exponential phase and reached a similar biomass to the controls after18days culture(see Supplementary Figure I, Supporting Information).On the basis of these data,proteomic and metabolomic profiling was conducted at early time points (0.5,2.0,and6.0h)and is described below.Proteomic Profiling.A total of3284peptides(Supplemen-tary Table1,Supporting Information)were detected that corresponded to759proteins(minimal set)identified with two or more peptides and an additional583proteins if relying on only one peptide for identification.Only those proteins with2 or more peptide matches and with a p<0.05(in both Mascot and Paragon)were reported for greater confidence.The list of peptides identified and the set of proteins are presented in Supplementary Table1(Supporting Information).Figure2 illustrates the number of proteins identified based on the number of peptides that matched to that particular protein.AOryza sativa Cell Suspension Culture Proteomics research articlesJournal of Proteome Research•Vol.9,No.12,20106625pie-chart of the functional classifications of these proteins is presented in Figure 3.One of the interesting features of this data set was the identification of proteins that are predicted to have a relatively low abundance (including,for example,Skp1,pathogenesis related protein,14-3-3proteins and the cell division cycle protein 48),alongside proteins with a predicted,relatively high abundance (including,for example,enzymes of primary metabolism,for example,glycolytic and citric acid cycle pathways).The largest functional classification group,representing 38%of the proteins identified,was proteins involved in primary metabolic processes,including glycolysis and the citric acid cycle.The second largest group was proteins with an unknown function (18%)that were unable to be matched to any other proteins of known function,after BLAST searching against the NCBInr database.Proteins involved in stress responses con-stituted 7%of the data set.Translation factors,elongation factors and translation initiation factors were categorized into a single group “protein biosynthesis”,which included 6%of the proteins identified.Five percent of the proteins were involved in amino acid metabolism and another 4%included proteins that were involved in DNA or RNA binding categorized as “nucleic acid binding”proteins.Proteins involved in protein folding and transport made up to 3%each of identified proteins.Proteins responsible for maintaining the cell homeo-stasis by regulating the redox status constituted 3%and those involved in nucleotide biosynthesis included another 2%of the proteins.Other groups included proteins involved in the production of ATP -energy proteins (1%),proteins involved in nucleotide binding (1%),protein modification (1%),transcrip-tion (1%),cell cycle (<1%),cellular organization (1%),nucleo-some assembly (1%),protein degradation (<1%)and storage proteins contributed toward 1%of the data set.Relative Protein Abundance.Of the 759proteins identified,92proteins were able to be quantified after collating the data from the three biological replicates over a time period of ABAFigure 1.Normalized expression levels of OsMyo1,OsCab26,OsLIP19and OsSPHK gene transcripts.Q-PCR was used to examine the transcript abundance of the above genes following exposure to ABA over 0,2,4,6,and 8h.Light gray bars indicate expression levels of genes after 100µM ABA treatment and the dark gray bars are the levels of expression of control samples (treated with 100µL ethanol).Error bars represent (SD,n )3replicates.research articlesRao et al.6626Journal of Proteome Research •Vol.9,No.12,2010treatment at 0.5,2,and 6h relative to an untreated control,using the stringent criteria indicated in the methods section.The 92quantified proteins are listed in Supplementary Table 2(Supporting Information).The distribution of protein abundances is depicted in Figure 4.Eleven proteins (11.9%),which included phosphoenolpyru-vate carboxylase,sucrose synthase 2and OsGDI1(the three with the greatest changes)and aminopeptidase C and aconitate hydratase,show increased abundances of >1.25fold at at least one time-point and just 1protein (OSIGBa0148P16.5;1%)shows a decrease of <0.8fold at the 0.5h time point.A BLAST search of this protein revealed it to be similar to an unchar-Figure 2.Proteins assigned based upon the number of peptides identified.The graph shows the number of proteins assigned based on the number of peptides identified.The x -axis represents the number of peptides assigned to a single protein.The y -axis represents the number of proteins identified.Figure 3.Functional classification of proteins identified.The pie chart includes the 759proteins identified (see Supplementary Tables 1and II,Supporting Information)on the basis of two or more peptides and classified based on known and predicted functions as specified either in NCBI or Uniprot.Figure 4.Relative protein abundances.Box plots showing distribution of iTRAQ protein ratios on a log scale across the time-course of ABA treatment.All 92proteins are represented in this figure.Boxes define the 25th and 75th percentiles of the popula-tion.Asterisks indicate extreme outlying values (*)phosphe-nolpyruvate carboxylase,**)sucrose synthase 2,***)OsG-DI1).Oryza sativa Cell Suspension Culture Proteomicsresearch articlesJournal of Proteome Research •Vol.9,No.12,20106627acterized protein,Os04g0386600.It is however important to note that the above-mentioned proteins showed consistent increases considering the average values of the three biological replicates at the three different time points.However,25other proteins after a student t test showed a p value of <0.02with consistent changes across all three biological replicates were also of interest.Most of these proteins included primary metabolic enzymes and proteins involved in oxidative stress responses described in the Discussion (see ABA-responsive metabolites,proteins and pathways).Metabolite Profiling.The metabolite responses of rice SCCs to ABA (100µM)was examined over a time period of 0.5h,2and 6h relative to an untreated control.Metabolite levels were determined using GC-MS and a total of 148metabolites were identified which included 25amino acids (AAs),45organic acids (OAs),35sugars,19fatty acids,2polyamines,4sterols,5sugar acids,4sugar alcohols and 9miscellaneous compounds (Supplementary Table 3,Supporting Information).This is comparable to the 129compounds identified in the analysis of rice leaves.18In part,this may reflect the metabolic com-plexity of the SCC system under examination,but also,the increased maturity of GC-MS based system and compound libraries that were used for searching the data to identify metabolites.The identified metabolites and their fold-changes relative to the control samples are presented in Supplementary Table 3(Supporting Information)with a standard error to the mean ((SE).Although many metabolites showed a p-value <0.05after a Student’s t test,the Bonferroni correction was used as the threshold to reduce false discovery rate.The values that were found to be statistically significant after Bonferroni correction 26,27are presented in bold font in Supplementary Table 3.A surprisingly small number of metabolites (13;8.78%)showed statistically significant changes in response to ABA at one or more time points,indicating that the response was quite specific.Of the 13significantly altered metabolites,8displayed significant decreases at 0.5h,2at 2h,while 10were signifi-cantly changed after 6h (includes metabolites that showed significant changes at more than one time point).The distribution of metabolite abundances are depicted in Figure 5.Unlike the protein abundance distribution,no single metabolites displayed consistently altered abundance at more than one time-point.Some of the outlier metabolites (most significantly changed)include glucuronate,erythritol,2-oxo-glutarate,galactose,aspartate and fructose-6-phosphate (Figure 5).The abundance of fumarate decreased significantly at both the 2h (0.7fold)and 6h (0.2fold)time points.This compound is one of the core metabolites of the citric acid cycle.Notably,three of the four other metabolites identified that are part of the citric acid cycle,citrate,isocitrate,malate and succinate also displayed decreased abundance following exposure to ABA,although none of these decreases were statistically significant.Tyrosine (Tyr)showed a decrease in abundance at the 0.5h (0.7fold)and 6h (0.6fold)time points.This metabolite can be decarboxylated to form tyramine through the action of a tyrosine decarboxylase (EC 4.1.1.28).29Three other organic acids displayed significant decreases in abundance during ABA treatment,4-hydroxybenzoate,itaconate and propanoic acid.4-Hydroxybenzoate,involved in the degradation of benzoate,displays a similar pattern of decreased abundance as Tyr.Benzoate also displayed a decrease in abundance following exposure to ABA,although it was not statistically significant.The abundance of itaconate significantly decreased at the 0.5h (0.5fold)and 6h (0.3fold)time points.cis -Aconitate is converted to itaconate releasing CO 2by the action of the enzyme cis -aconitate carboxylyase.29cis -Aconitate also dis-played a large decrease at the 6h time points,however,it increased at the 0.5h time point.Mannitol,a sugar alcohol,involved in fructose (Fru)and mannose metabolism 29was slightly reduced in abundance at the 0.5and 2h time point but displayed a significant (<0.6fold)reduction in abundance at the 6h time point.Fru-6-phosphate (Fru-6-P)increases rapidly in the first 0.5h (6-fold)and then gradually decreases over the 2and 6h time points.DiscussionRice SCCs were established as a model system to study ABA responses by studying known ABA responsive genes.Three of the four in planta ABA responsive genes examined,OsMyo1,OsCab26and OsLIP19,responded to ABA treatment within 4h of application.As reported by Lin et al,28our results are consistent in showing that the OsMyo1transcript responded to ABA in a manner similar to whole plants by increasing within 4h.28Cab26is both negatively regulated by ABA in soybean cotyledons during embryogenesis 30and induced by ABA in rice shoots.28In rice SCCS the OsCab26transcript levels increased at 2h of ABA treatment.These reports suggest that the Cab genes are differentially expressed during different stages of plant development and are under different mechanisms of regulation.We also observed the induction of OsLIP19at 2h further illustrating that SCCS can mimic whole plants.In cDNA macroarray and RT-PCR experiments OsSPHK transcript did increase after 2h of ABA treatment.28From the Q-PCR data obtained (Figure 1D),there was no increase in the transcript abundance of OsSPHK over the time-course of ABA treatment.This could be because the transcript was activated within 2h of treatment (first time-point analyzed in the Q-PCR)with ABA similar to that previously observed for rice shoots.28These results indicate that the SCCs respond to ABA treatment in a similar,but not always identical manner to that reported in planta .We therefore concluded it was a useful model to apply an ’omics technology approach to studying the proteome and metobolome responses to ABA.Figure 5.Relative metabolite abundances.Box plots showing distribution of relative metabolite fold change values on a log scale across the time-course of ABA treatment.All 148metabo-lites are represented in this figure.Boxes define the 25th and 75th percentiles of the population.research articlesRao et al.6628Journal of Proteome Research •Vol.9,No.12,2010。