elucidation of complex nature of peg induced drought ... · elucidation of complex nature of peg...
TRANSCRIPT
ORIGINAL RESEARCHpublished: 29 September 2016doi: 10.3389/fpls.2016.01466
Frontiers in Plant Science | www.frontiersin.org 1 September 2016 | Volume 7 | Article 1466
Edited by:
Paul Christiaan Struik,
Wageningen University and Research
Centre, Netherlands
Reviewed by:
Abu Hena Mostafa Kamal,
University of Texas at Arlington, USA
Sanjay K. Gupta,
University of Minnesota, USA
Supratim Basu,
New Mexico Consortium, USA
*Correspondence:
Chandra S. Nautiyal
†These authors have contributed
equally to this work.
Specialty section:
This article was submitted to
Crop Science and Horticulture,
a section of the journal
Frontiers in Plant Science
Received: 23 July 2016
Accepted: 14 September 2016
Published: 29 September 2016
Citation:
Agrawal L, Gupta S, Mishra SK,
Pandey G, Kumar S, Chauhan PS,
Chakrabarty D and Nautiyal CS (2016)
Elucidation of Complex Nature of PEG
Induced Drought-Stress Response in
Rice Root Using Comparative
Proteomics Approach.
Front. Plant Sci. 7:1466.
doi: 10.3389/fpls.2016.01466
Elucidation of Complex Nature ofPEG Induced Drought-StressResponse in Rice Root UsingComparative Proteomics ApproachLalit Agrawal, Swati Gupta †, Shashank K. Mishra †, Garima Pandey, Susheel Kumar,
Puneet S. Chauhan, Debasis Chakrabarty and Chandra S. Nautiyal *
Division of Plant Microbe Interactions, Council of Scientific and Industrial Research-National Botanical Research Institute,
Lucknow, India
Along with many adaptive strategies, dynamic changes in protein abundance seem
to be the common strategy to cope up with abiotic stresses which can be best
explored through proteomics. Understanding of drought response is the key to decipher
regulatory mechanism of better adaptation. Rice (Oryza sativa L.) proteome represents
a phenomenal source of proteins that govern traits of agronomic importance, such
as drought tolerance. In this study, a comparison of root cytoplasmic proteome was
done for a drought tolerant rice (Heena) cultivar in PEG induced drought conditions.
A total of 510 protein spots were observed by PDQuest analysis and 125 differentially
regulated spots were subjected for MALDI-TOFMS-MS analysis out of which 102 protein
spots identified which further led to identification of 78 proteins with a significant score.
These 78 differentially expressed proteins appeared to be involved in different biological
pathways. The largest percentage of identified proteins was involved in bioenergy and
metabolism (29%) and mainly consists of malate dehydrogenase, succinyl-CoA, putative
acetyl-CoA synthetase, and pyruvate dehydrogenase etc. This was followed by proteins
related to cell defense and rescue (22%) such as monodehydroascorbate reductase
and stress-induced protein sti1, then by protein biogenesis and storage class (21%)
e.g. putative thiamine biosynthesis protein, putative beta-alanine synthase, and cysteine
synthase. Further, cell signaling (9%) proteins like actin and prolyl endopeptidase, and
proteins with miscellaneous function (19%) like Sgt1 and some hypothetical proteins
were also represented a large contribution toward drought regulatory mechanism in rice.
We propose that protein biogenesis, cell defense, and superior homeostasis may render
better drought-adaptation. These findings might expedite the functional determination of
the drought-responsive proteins and their prioritization as potential molecular targets for
perfect adaptation.
Keywords: drought, proteomics, rice root, tolerant cultivar, 2-DE, MALDI-MS/MS, SOTA analysis
Agrawal et al. Drought Affected Rice Root Proteome
INTRODUCTION
Rice (Oryza sativa L.) is an important staple cereal crops rankedsecond after maize worldwide and considered as a primary sourceof food for more than half the world’s population includingAsia. India is the second largest producer of rice after Chinaand account for about 20% of all world rice production. Rice isexclusively grown for human consumption and therefore there isimmense pressure for high production to feed the largely growingworld population (Ahsan et al., 2008). Climate related diverseabiotic stresses (drought, flood, salt, cold etc.) are the principalsources of risk and uncertainties in agriculture and causes widefluctuations in agricultural output. While attempts have beenmade to reduce the adverse effects of weather on agriculturethrough scientific research and technology development, theperformance of agriculture, especially in developing countries,still depends largely on the weather. Drought is the majorabiotic stress limiting the productivity of rice as it hampers plantgrowth and development, and shrinks harvest size (Subba et al.,2013). In Asia, drought, a major abiotic stress, is responsible foraffecting about 20% of the total rice-growing area (Pandey andBhandari, 2008). Rice germplasm exhibit a tremendous geneticsource for controlling agronomical importance characters suchas drought tolerance. Moreover, relatively smaller genome sizeand well-organized database make it a choice of model plant formonocots to study the physiological, biochemical and molecularaspects during development, abiotic as well as biotic stressconditions (Atkinson and Urwin, 2012; Narsai et al., 2013). Thebetter understanding of drought tolerance mechanism and theefficiency to develop drought tolerant varieties can be correlatedwith identification of trait associated genes (Tuberosa and Salvi,2006; Lafitte et al., 2007; Sreenivasulu et al., 2007). However, itremains a challenging task due to its highly complex network ofseveral metabolic pathways (Price et al., 2002). Understandingthe mechanism of dehydration response is the key to decipherthe regulatory mechanism of better adaptation (Subba et al.,2013).
Plant roots are important organs to uptake soil water andnutrients, perceiving and transducing of soil water deficit signalsto shoot (Davies and Zhang, 1991; Moumeni et al., 2011) whichfurther triggers an array of physiological, morphological and
Abbreviations: SAT, serine acetyl transferase; CS, cystein synthesae; CGL,Cystathionine γ-lyase; CBS, Cystathionine β-synthase; MS, methionine synthase;SAM, S-adenosylmethionine; SAMS, S-adenosylmethionine synthesae; G-6-P,Glucose 6 phosphate; PFK, pyrophosphate-dependent phosphofructokinase;ALD, aldehyde dehydrogenase; TPI, triose phosphate isomerase; GAPDH,glyceraldehyde 3-phosphate dehydrogenase; PGK, phosphoglycerate kinase; PGM,phosphoglycerate mutase; ENO, enolase; PK, pyruvate kinase; GLO1, glyoxalase 1;PDH, Pyruvate dehydrogenase; CoAS, Co-A synthesae; AH, aconitate hydratase;IDH, isocitrate dehydrogenase; SucCoA, Succinyl Co-A synthese; SD, succinateDehydrogenase; FR, fumarase; MD, malate dehydrogenase; GR, glutathionereductase; DHAR, dehydroascorbate reductase; APx, ascorbate peroxidase;MDAR, monodehydroascorbate reductase; POX, peroxidases; SOD, superoxidedismutase; HSPs, heat shock proteins; ACT, actin; SGT1, suppressor of the G2allele of skp1; PR, pathogenesis related protein; OAC, oryzain alpha chain; NiBP, Nibinding protein; STI1, stress inducible protein; PEP, phosphoenolpyruvate; HYPO,hypothetica.
molecular responses in the whole plant. Understanding theroot responses through transcriptomics with parallel insightson protein level changes has always been an area of interestand mostly relies on comparative studies of diverse geneticbackground under drought (Sengupta et al., 2011). Roots aresupposed to be the primary site for sensing drought stress toinitiate signaling cascade at molecular level in responses todrought. Plants can successfully adopt a tolerance mechanismagainst environmental stress by regulating gene expressionthrough transcriptional changes in regulatory and functionalproteins (Périn et al., 2007). In a study, Kawasaki et al. (2001)observed that gene expression profiling of rice under salt stressmay result in up and downregulation of∼10% of the transcripts.However, the gene expression studies do not offer insights intothe quantity and quality of the proteins as the amount ofproteins is not always correlated to that of mRNA. The post-translational modifications like phosphorylation, glycosylationand removal of signal peptides enable proteins for activities andsubcellular localization (Kawasaki et al., 2001). As a consequence,understanding the stress mechanism at protein level may providesome more useful insights to study the actual biotic stressresponse. In this context, evolution of proteomics is playing acrucial role as necessary and complementary approach to addresssuch issues in the post-genomic era (Zivy and de Vienne, 2000;van Wijk, 2001). Serious effort has been made for functionalidentification of tissues, organs, and development specific ricegene under environmental changes like biotic and abioticstresses using systematic studies in proteomic analysis (Khanand Komatsu, 2004; Komatsu, 2005; Ali and Komatsu, 2006)and knowledge of these proteins would help in understandingthe stress tolerance related molecular mechanism in rice at thetranslation level.
We report here the comparative proteomic analysis ofrice root for PEG simulated drought responsiveness in timedependent manner. Proteins were separated by two-dimensionalgel electrophoresis (2-DE) followed by PD Quest analysis andidentified by MALDI-Mass Spectrometry (MS) using availableproteome databases.
MATERIALS AND METHODS
Plant Material CollectionAfter screening of several local rice varieties, a drought tolerantvariety Heena cultivated in northern Indian was selected for thepresent study. Three-weeks-old seedlings of rice were stressed for7 days using 20% polyethylene glycol PEG-6000 in the nutrientsolution. Exposure to PEG-6000 solutions is supposed to mimicdrought stress with limited metabolic interferences. Effectivescreening of large sets of germplasm for drought tolerance hasbeen carried out by PEG-based in vitro screening as a suitablemethod with good accuracy (Muscolo et al., 2014). Roots samplesfrom three independent biological replicates were harvested onfirst, third, and seventh day and immediately flash frozen inliquid nitrogen and stored independently at −80◦C for furtheranalyses.
Frontiers in Plant Science | www.frontiersin.org 2 September 2016 | Volume 7 | Article 1466
Agrawal et al. Drought Affected Rice Root Proteome
Extraction of Rice Root Proteins and2-Dimensional Gel Electrophoresis (2DE)AnalysisRice roots harvested in three replications were pooled tonormalize the effect of variations in the biological replicates andused for extraction of soluble proteins. Root tissue (5 g) waspulverized to a fine powder in liquid nitrogen and suspended in10 ml lysis buffer containing 0.5 M Tris HCl (pH 7.5), 0.1 M KCl,50 mM EDTA (pH 8.0), 0.7 M Sucrose, 10 mM thiourea, 2 mMPMSF, and 2% β-mercaptoethenol. Suspension was mixed with10 ml of Tris saturated phenol (pH 8.0) for 30 min with gentleshaking and centrifuge at 9000 × g for 10 min. Upper phasewas transferred in separated vial and 0.1M ammonium acetate(dissolved in methanol) added in upper phase and incubatedat −20◦C for overnight. Next day, the solution was centrifugedagain at 9000 × g and supernatant was discarded. Pellet waswashed with 100% chilled acetone twice and then air dried pelletwas dissolved in rehydration buffer [8 M urea, 2 M thiourea, 4%(w/v) CHAPS, 20 mMDTT, 0.5% (v/v) Pharmalyte (4–7)]. Totalroot protein concentration was determined by Bradford assay(Bio-Rad, USA). Isoelectric focusing was carried out with 250 µgof protein as given by Agrawal et al. (2008). Briefly, the 13-cmIEF strips (pH 4–7) were rehydrated with protein for 16 h andelectrofocused using an IPGphor system (Bio-Rad, USA) at 20◦Cup to 25,000 Vh. This focused strips were subjected to reductionwith 1% (w/v) DTT in 10 mL of equilibration buffer [6 M urea,50 mM Tris-HCl (pH 8.8), 30% (v/v) glycerol and 2% (w/v) SDS]and followed by alkylation using 2.5% (w/v) iodoacetamide insame equilibration buffer. The strips were then loaded on topof 12.5% polyacrylamide gels for SDS-PAGE. The gels were fixedand stained with a silver stain plus kit as per protocol (Bio-Rad,USA).
Image Acquisition and Data AnalysisThe silver stained gel obtained after 2DE were scanned for imageacquisition using the Bio-Rad FluorS system equipped with a12-bit camera. Images from three biological replicate 2-DE gelswere taken for the data analysis in PDQuest version 8.0.1 (Bio-Rad). The parameters like protein spot quality, molecular mass,and pI of individual protein was assessed as described earlierby Agrawal et al. (2013). The spots showing reproducibility inquality and quantity in at least two of the three replicate gelswere taken in consideration during analysis. A normalized spotvolume for protein quantification was obtained by PDQuestsoftware using the total spot volume normalization procedureto avoid experimental variations in gels due to protein load orstaining.
Protein IdentificationThe silver stained protein spots were digested forMALDIMS/MSwere digested with trypsin (Promega Corporation, MA, USA)as per manufacturer instruction. The stained protein spots wereexcised manually, washed thoroughly with water and crushedin small pieces. These gel pieces were washed with 15 mMpotassium ferricynide and 50 mM sodium thiosulfate (1:1) fordestaining. Then gel pieces were dehydrated with solution A
[100% Acetonitrile (ACN): 50 mM Ammonium bicarbonate(ABC) in 2:1], subsequently rehydrated with 25 mM ABC.These steps were repeated until gel pieces become white. Finalsupernatant was discarded and gel pieces were dried in speedvac.After this gel pieces were rehydrated with 0.2 µg/µL trypsin and50mM ammonium bicarbonate and incubated at 37◦C overnightfor digestion. Peptides were extracted from gel slices with 1%TFA in 50% ACN twice. The peptide solution was concentrateddown to 5µL by vacuum centrifuge. Peptide solution mixed with5 mg/ml CHCA (α-Cyano-4-hydroxycinnamic acid) matrix andspotted onto MALDI plate.
MALDI-TOF-TOF AnalysisA 4800 Proteomics Analyzer (Applied Biosystems, USA) withTOF/TOF optics was used for all MALDI-MS and MS/MSapplications. Samples for MALDI were prepared by adding 0.5µl matrix solution (5 mg/mL a-Cyano-4-hydroxycinnamic acidin 50% ACN containing 0.1% TFA) to 0.5 µl trypsin digestedprotein sample and left for air dry at room temperature onstainless steel 384 well-target plate after spotting. The platecontaining spot was inserted in the mass spectrometer andsubjected tomass analysis. Amixture of four proteins angiotensinI, Glu-fibrino-peptide B, ACTH (1e17), and ACTH (18e39) wasused to calibrated the mass spectrometer externally. Further,the instrument was externally calibrated with fragment ofGlufibrino-peptide B for MS/MS experiments.
The monoisotopic peptide masses obtained from MALDI-TOF/TOF were analyzed by the 4000 Series Explorer softwareversion 3.5. Spectra were collected in a data dependent mode.For each spot total 500 laser shots were accumulated to generateions for MS analysis and 30 most abundant ions were usedfor MS/MS analysis. All exported spectra as Mascot GenericFormat (MGF) files from MALDI platforms were acquired inlinear positive mode, smoothed by the Savitzky–Golay algorithmand the baseline subtracted for further analysis. Peak detectioncriterion for peak to be considered was set to minimum S/N =
10. Protein identification was performed using Mascot software(http://www.matrixscience.com) NCBInr databases 20160522(87959973 sequences; 32280238673 residues). The databasesearch criteria were as follows: taxonomy,O. sativa (rice) (172275sequences), peptide tolerance, ±100 ppm, MS/MS tolerance,±0.2 Da; peptide charge+1; maximum allowed missed cleavage,1; fixed modification, cysteine carbamidomethylation; variablemodification, methionine oxidation; instrument type, MALDI-TOF/TOF. Protein scores were measured from sum of the seriesof peptide scores as a non-probabilistic basis for ranking proteinhits. The only protein spots whose MOWSE score was abovethe significant threshold level (Table S3) determined by Mascotalong with the number of peptides matched and % coverageof matched protein were considered to represent a positiveidentification. In all the protein identifications, probability scoreswere greater than the score fixed by Mascot as significant with ap < 0.05.
The similar expression pattern of the identified differentialproteins was determined by SOTA (self-organizing treealgorithm) clustering on the log transformed fold inductionexpression values using Multi Experiment Viewer (MEV)
Frontiers in Plant Science | www.frontiersin.org 3 September 2016 | Volume 7 | Article 1466
Agrawal et al. Drought Affected Rice Root Proteome
pH4 pH7
Control Day-1 Day-3 Day-7
pH4 pH7 pH4 pH7 pH4 pH7
pH4 pH7
A
B
C
FIGURE 1 | 2-DE analysis of root proteome of rice. (A) Rice plant showing different stages of drought. (B) Proteins were extracted from rice root tissue and equal
amounts (250 µg) of proteins were separated by 2-DE as described in Materials and Methods Section. (C) Three replicate silver-stained gels for each stage were
computationally combined using PDQuest software and one representative master standard gel image was generated.
software (The Institute for Genomic Research, TIGR) of rice.The clustering was performed with Pearson correlation asdistance with 10 cycles and maximum cell diversity of 0.8(Romijin et al., 2005). SOTA is a neural network that growsadopting a binary tree topology, where hierarchical cluster areobtained as result with the accuracy and robustness (Herreroet al., 2001). Each branch summarizes the patterns of allsimilarly expressed proteins with centroid expression profile of agroup.
Expression Analysis of Some SelectedCandidate Genes Using qRT PCRReal time-PCR was performed in 20 µl for a set of selected genesusing Fast SYBR Green PCR Master Mix (Agilent Technologies,USA). The list of selected genes and oligonucleotide primers(Sigma-Aldrich, USA) used for each gene are listed in Table S1.Oligonucleotide primers for vetiver ubiquitin gene were used asthe internal control for establishing equal amounts of cDNA inall reactions. The reactions were performed using the followingcycle conditions, an initial 94◦C for 2 min, followed by 30 cyclesof 94◦C for 30 s, 60◦C for 30 s, and 72◦C for 30 s, and the final5 min extension at 72◦C. After obtaining the ct-value for each
reaction, the relative expression was calculated by 2ˆ-delta Ctmethod.
RESULT AND DISCUSSION
Dehydration-Induced Changes and 2-DEAnalysisTo understand the drought tolerance mechanism in rice, wecarried out the comparative proteomic analysis of roots of thedrought tolerance Heena cultivar of rice at various time pointsafter drought induction in 10% PEG (Figure 1A). 2DE wasperformed at pH 4–7 IpG strip from 21 days old hydroponicallygrown rice root tissue in three replicates and images wereanalyzed by the PDQuest software version 8.0.1 as describedabove (Figure 1B). The mean value of the high-quality spots wasused as the spot quantity on the master gel (Figure 1C, Table S2).For comparative proteomics study minimum 2.5 fold changein protein expression either increase or decrease at any oneof the stages was considered for differential spot identification.A criterion of p < 0.01 was used to determine the significantdifference for analyzing the parallel spots between genotypewith analysis of one-way variance (ANOVA). A total of 125
Frontiers in Plant Science | www.frontiersin.org 4 September 2016 | Volume 7 | Article 1466
Agrawal et al. Drought Affected Rice Root Proteome
8717
7510
8703
8616
8601
8506
8412
8108
8003
7712
7614
7611
7609
7409
7314
7711
7111
7103
6607
65186414
6412
6201
6104
5717
5711
5504
5418
54095307
6206205205
4815 4808
4803
4510
4107
42024203 4207 4213
4302
37033618
3518
3516
3509
3413
3409
3318
43023413
3317
3314
3115 3111
2617
2610
2220
2201
261261
2107
1617
1007
12083218
2218220
7524
5107
2214
2313
3113
4315
6105
67106713
8408
3312
3417
641
3417
44144416
7617
22121215
12111212
1512
1001003
7619
2221
7210
2516
4104107
4112
8607
3416
4607
7
510
4107
22 5105
5718
5108
5011
831
6
3418
pH4 pH7
92
57
35
28
21
KDa
113
FIGURE 2 | A higher level master image was created in silico by
PDQuest from three replicate gels of each time point. The boxed areas
marked with dotted lines represent the zoomed-in gel sections in Figure 3.
The numbers correspond with the spot IDs listed in Table 1.
differentially expressed drought responsive protein spots fromall time points were subjected to MALDI-TOF/TOF analysis.Out of which 78 spots were identified significantly as indicatedby arrows on the higher level matchset image generated by PDQuest (Figure 2) and summarized in Table S3. Some typical gelregions representing protein spots with altered expression areenlarged and shown in Figure 3. These protein spots actuallyaccount for 45 distinct proteins (Table 1), suggesting 57% uniqueprotein identifications, while the remaining 43% of the identifieddifferentially expressed proteins either correspond to post-translationally modified forms or may be members of multigenefamilies. Of the 78 identified differentially expressed proteins,30 protein spots were clearly up-regulated, and 26 were down-regulated, while 22 of the protein spots showed a mixed patternof development stage dependent expression.
Functional Distribution of DehydrationResponsive ProteinsTo understand the function of the proteins associated withdrought induced changes in the rice roots, the differentiallyexpressed proteins were sorted into different functionalcategories. Seventy-eight identified differentially expressedproteins could be assigned to five functional classes based ontheir putative roles in the drought-response (Figure 4, Table 1).In case of Heena the largest percentage of the identified proteinswas involved in bioenergy and metabolism (29%), cell defenseand rescue (22%), protein biogenesis and storage (21%), cellsignaling 9%, and miscellaneous (19%; Figure 4). In a numberof cases, many proteins were represented by multiple isoelectric
4808
3312
3317
4808 4808 4808
3312 3312 3312
3317 3317 3317
87178717 8717 8717
7711
7614
7609
22201215
1211
1212 1212
1211
1215 2220
1212 1212
22201215
1211
22201215
1211
7609
7711
7614
7609 7609
7614
7711
7614
7711
6201
5718571857185718
6201 62016201
Control Day 1 Day 3 Day 7
G
H
E
D
C
B
A
F
FIGURE 3 | Some of the differentially expressed root proteins
represented in the enlarged gel sections (A–H) corresponds to the
marked boxed areas in Figure 2.
forms (Table 1), suggesting the possible post-translationalmodification of the candidate protein.
Bioenergy and Metabolism (BEM)The BEM class contained 24 proteins which appeared asdifferential spots. ATP synthase subunit d (OsC-3113), malatedehydrogenase, mitochondrial-like (OsC-5107), succinyl-CoAligase [ADP-forming] subunit beta (OsC-3218), putative acetyl-CoA synthetase (OsC-8316), V-type proton ATPase subunitB 1-like (OsC-3618), glutamate dehydrogenase 2 (OsC-8412),putative succinate dehydrogenase flavoprotein alpha subunit(OsC-8616) were upregulated, while pyruvate dehydrogenase E1component subunit beta-2 (OsC-3417), malate dehydrogenase,
Frontiers in Plant Science | www.frontiersin.org 5 September 2016 | Volume 7 | Article 1466
Agrawal et al. Drought Affected Rice Root Proteome
TABLE 1 | List of differentially expressed rice root proteins identified by MS/MS.
S. No. Spot IDa Protein name GI no. Protein No. of % coverage Expression patternc Experimental Theoritical
scoreb peptides kDa/PI kDa/PI
BIOENERGY AND METABOLISM (BEM)
1. OsC-3417 Pyruvate dehydrogenase E1
component subunit beta-2,
mitochondrial
GI:115480067 48 2 7 40.3/5.3 64.4/5.3
2. OsC-4803 Pyruvate dehydrogenase E1
component subunit beta-2,
mitochondrial-like [Setaria italica]
GI:115434904 100 9 12 93.9/5.1 97.6/5.5
3. OsC-2107 Lecithin:cholesterolacyltransferase
family protein, expressed [Oryza
sativa Japonica Group]
GI:218194989 41 1 3 33.6/9.4 29.1/5.2
4. OsC-3113 ATP synthase subunit d,
mitochondrial-like isoform X2 [Setaria
italica]
GI:115476908 58 3 20 19.7/6.2 9.0/5.4
5. OsC-5107 Malate dehydrogenase,
mitochondrial-like [Oryza
brachyantha]
GI:115456241 56 2 7 15.9/3.2 29.2/5.8
6. OsC-7524 Malate dehydrogenase,
mitochondrial-like [Oryza
brachyantha]
GI:115438875 394 6 20 35.6/8.7 76.5/6.3
7. OsC-7711 Sucrose synthase 1-like isoform X2
[Oryza brachyantha]
GI:115453437 170 6 7 93.3/5.9 90.6/6.4
8. OsC-3218 Succinyl-CoA ligase [ADP-forming]
subunit beta, mitochondrial-like
[Oryza brachyantha]
GI:115447367 363 11 26 45.4/5.9 47.7/5.5
9. OsC-8316 Putative acetyl-CoA synthetase
[Oryza sativa Japonica Group]
GI:49388286 95 4 6 78.5/5.6 57.6/6.4
10. OsC-1007 Succinate dehydrogenase
[ubiquinone] flavoprotein subunit,
mitochondrial; [Oryza sativa Japonica
Group]
GI:75135397 59 2 3 69.4/6.6 28.4/4.9
11. OsC-8616 Putative succinate dehydrogenase
flavoprotein alpha subunit [Oryza
sativa Japonica Group]
GI:34394418 49 2 3 69.4/6.6 86.3/6.6
12. OsC-2617 Probable nucleoredoxin 1-1; Short,
OsNrx1-1 [Oryza sativa Japonica
Group]
GI:115453457 56 1 2 64.0/4.9 80.4/5.1
13. OsC-3314 Full, Fructokinase-2; [Oryza sativa
Japonica Group]
GI:115474481 90 3 9 35.8/5.0 54.8/5.4
14. OsC-3409 Putative transaldolase [Oryza sativa
Japonica Group]
GI:115441963 72 3 7 46.5/6.0 69.5/5.4
15. OsC-3618 V-type proton ATPase subunit B
1-like [Oryza brachyantha]
GI:115468606 236 14 31 54.1/5.0 79.4/5.3
(Continued)
Frontiers in Plant Science | www.frontiersin.org 6 September 2016 | Volume 7 | Article 1466
Agrawal et al. Drought Affected Rice Root Proteome
TABLE 1 | Continued
S. No. Spot IDa Protein name GI no. Protein No. of % coverage Expression patternc Experimental Theoritical
scoreb peptides kDa/PI kDa/PI
16. OsC-4302 Spermidine synthase 1 [Oryza sativa
Japonica Group]
GI:6468656 76 2 9 27.9/5.1 50.3/5.5
17. OsC-5418 Photosystem II subunit D1, partial
(chloroplast) [Polyscias fruticosa]
GI:345105629 46 2 100 30.9/5.11 62.2/5.9
18. OsC-7103 Photosystem II subunit D1, partial
(chloroplast)
GI:345105629 45 2 100 17.2/5.1 28.9/6.1
19. OsC-5504 Enolase GI:90110845 737 11 42 48.2/5.4 70.3/5.7
20. OsC-7510 Pyrophosphate–fructose
6-phosphate 1-phosphotransferase
subunit beta-like [Oryza brachyantha]
GI:115467370 180 5 10 61.9/6.0 78.0/6.4
21. OsC-8108 Carbonic anhydrase [Oryza sativa] GI:3345477 44 2 10 29.4/8.4 35.4/6.8
22. OsC-8412 Glutamate dehydrogenase 2 [Oryza
sativa Japonica Group]
GI:81686712 265 4 12 44.9/6.3 62.0/6.8
CELL DEFENSE AND RESCUE (CDR)
23. OsC-3416 TPA: class III peroxidase 72 precursor GI:55701011 49 2 6 37.8/5.5 59.6/5.3
24. OsC-4416 TPA: class III peroxidase 72 precursor GI:55701011 42 2 6 37.8/5.5 66.5/5.5
25. OsC-6412 TPA: class III peroxidase 72 precursor
[Oryza sativa Japonica Group]
GI:55701011 60 3 11 37.8/5.5 67.2/5.9
26. OsC-8506 TPA: class III peroxidase 72 precursor GI:38426301 71 4 12 51.7/6.5 71.1/6.5
27. OsC-3703 Heat shock protein 90 GI:39104468 76 2 3 80.4/4.9 89.1/5.2
28. OsC-8607 Heat shock protein STI-like
[Brachypodium distachyon]
GI:115447567 140 5 10 65.1/6.0 84.3/6.6
29. OsC-4112 24.1 kDa heat shock protein,
mitochondrial-like [Oryza
brachyantha]
GI:115448791 61 3 10 24.0/6.8 29.5/5.7
30. OsC-4315 Glyoxalase I [Oryza sativa Japonica
Group]
GI:115475151 185 10 36 35.1/5.5 55.1/5.6
31. OsC-5409 Cytosolic
monodehydroascorbatereductase
[Oryza sativa Japonica Group]
GI:4666287 364 5 17 46.7/5.6 61.6/5.7
(Continued)
Frontiers in Plant Science | www.frontiersin.org 7 September 2016 | Volume 7 | Article 1466
Agrawal et al. Drought Affected Rice Root Proteome
TABLE 1 | Continued
S. No. Spot IDa Protein name GI no. Protein No. of % coverage Expression patternc Experimental Theoritical
scoreb peptides kDa/PI kDa/PI
32. OsC-7111 GSH-dependent
dehydroascorbatereductase 1 [Oryza
sativa Japonica Group]
GI:6939839 69 1 7 23.7/5.65 35.8/6.3
33. OsC-3418 Monodehydroascorbatereductase
[Oryza sativa Japonica Group]
GI:42407947 352 7 24 46.7/5.3 65.2/5.3
34. OsC-7314 Putative r40c1 protein—rice [Oryza
sativa Japonica Group]
GI:115452789 350 8 27 42.2/6.2 57.1/6.4
35. OsC-7611 Putative stress-induced protein sti1
[Oryza sativa Japonica Group]
GI:49388654 76 2 4 65.1/6.0 84.4/6.3
36. OsC-7617 Putative stress-induced protein sti1
[Oryza sativa Japonica Group]
GI:49388654 126 2 4 65.1/6.0 85.8/6.4
37. OsC-8601 Putative stress-induced protein sti1
[Oryza sativa Japonica Group]
GI:115447567 90 3 6 65.1/6.0 84.2/6.5
38. OsC-2313 Pathogenesis-related protein 1-like
[Oryza brachyantha]
GI:115474481 472 8 52 35.8/5.0 53.7/5.0
39. OsC-2212 Root specific pathogenesis-related
protein 10 [Oryza sativa Japonica
Group]
GI:38678114 179 8 41 17.0/4.8 41.9/4.9
PROTEIN BIOGENESIS AND STORAGE (PBS)
40. OsC-4607 Tricin synthase 1-like [Oryza
brachyantha]
GI:115477092 46 2 10 27.9/5.1 84.4/5.6
41. OsC-1208 Ribosome inactivating protein,
expressed [Oryza sativa Japonica
Group]
GI:29150390 36 1 2 55.6/9.7 46.2/4.9
42. OsC-1211 Elongation factor 1 beta’ [Oryza sativa
Japonica Group]
GI:218161 74 2 9 23.8/4.8 44.9/4.9
43. OsC-7712 Putative elongation factor 2 [Oryza
sativa Japonica Group]
GI:49387779 58 4 5 94.9/5.8 94.6/6.4
44. OsC-2220 Regulation of nuclear pre-mRNA
domain-containing protein 1B-like
isoform X2 [Oryza brachyantha]
GI:20160712 41 1 1 61.2/52.3 48.6/5.0
45. OsC-3111 Proteasome subunit beta type 3 GI:49388033 86 3 15 23.1/5.1 35.7/5.4
46. OsC-5205 Proteasome subunit alpha type-1-like
[Setaria italica]
GI:115444057 184 10 47 29.8/5.3 47.8/5.7
47. OsC-3516 Alpha-tubulin GI:1136120 79 1 3 50.4/4.8 69.5/5.4
48. OsC-4202 Putativetriosephosphateisomerase,
chloroplast precursor [Oryza sativa
Japonica Group]
GI:50725810 126 3 13 32.7/6.9 39.6/5.5
(Continued)
Frontiers in Plant Science | www.frontiersin.org 8 September 2016 | Volume 7 | Article 1466
Agrawal et al. Drought Affected Rice Root Proteome
TABLE 1 | Continued
S. No. Spot IDa Protein name GI no. Protein No. of % coverage Expression patternc Experimental Theoritical
scoreb peptides kDa/PI kDa/PI
49. OsC-4207 Cysteine synthase; GI:84028195 720 12 48 33.9/5.3 47.6/5.6
50. OsC-4815 Ubiquitin-activating enzyme E1 2,
putative, expressed [Oryza sativa
Japonica Group]
GI:108862075 45 2 1 109.3/5.2 97.5/5.5
51. OsC-7609 5-
methyltetrahydropteroyltriglutamate-
homocysteine methyltransferase,
putative, expressed
GI:108862990 95 4 8 79.2/7.1 87.9/6.3
52. OsC-7614 5-
methyltetrahydropteroyltriglutamate-
homocysteine methyltransferase,
putative, expressed
GI:108862990 299 7 15 84.6/6.0 88.8/6.4
53. OsC-8703 5-
methyltetrahydropteroyltriglutamate-
homocysteine methyltransferase,
putative, expressed
GI:108862990 91 4 9 48.2/7.1 88.9/6.5
54. OsC-3318 Putative thiamine biosynthesis protein GI:27261025 105 4 16 37.1/5.4 52.3/5.3
55. OsC-7409 Putative beta-alanine synthases
[Oryza sativa Japonica Group]
GI:22775671 54 2 5 46.0/6.0 62.5/6.4
SIGNALING (S)
56. OsC-2218 Actin GI:148886771 327 11 36 42.1/5.2 41.1/5.1
57. OsC-2221 Actin-1-like isoform X2 [Setaria italic] GI:115454971 638 9 33 41.8/5.2 46.5/5.0
58. OsC-5105 Actin-depolymerizing factor 4 GI:75243284 68 5 27 16.0/5.7 29.5/5.8
59. OsC-2214 Actin-depolymerizing factor 3 [Zea
mays]
GI:115489014 93 4 23 17.0/4.8 46.5/5.0
60. OsC-2201 14-3-3-like protein GF14-C; AltName:
G-box factor 14-3-3 homolog C
GI:115476520 83 2 11 28.9/4.7 44.2/5.0
61. OsC-4203 Putative membrane protein GI:9998903 120 4 13 31.5/5.0 47.0/5.5
62. OsC-5711 Prolylendopeptidase-like [Oryza
brachyantha]
GI:115470116 67 4 6 91.5/5.7 89.8/5.8
UNKNOWN (UK)
63. OsC-1215 Oryzain alpha chain; Flags: Precursor
[Oryza sativa Japonica Group]
GI:38345906 74 2 5 51.3/5.0 48.0/4.9
(Continued)
Frontiers in Plant Science | www.frontiersin.org 9 September 2016 | Volume 7 | Article 1466
Agrawal et al. Drought Affected Rice Root Proteome
TABLE 1 | Continued
S. No. Spot IDa Protein name GI no. Protein No. of % coverage Expression patternc Experimental Theoritical
scoreb peptides kDa/PI kDa/PI
64. OsC-2516 Sgt1 [Oryza sativa] GI:6581058 262 4 14 41.0/4.9 71.5/5.1
65. OsC-3115 Putative probable submergence
induced, nickel-binding protein 2A
GI:49388033 53 3 9 29.7/5.1 36.5/5.3
66. OsC-3312 Hypothetical protein OsI_37864 GI:125536157 45 2 6 28.9/5.0 56.1/5.4
67. OsC-3413 Hypothetical protein GI:54290369 45 1 9 17.2/9.2 59.8/5.4
68. OsC-3509 3’-N-debenzoyl-2’-deoxytaxol
N-benzoyltransferase-like [Oryza
brachyantha]
GI:115438576 129 3 8 46.1/5.1 69.9/5.4
69. OsC-3518 Hypothetical protein OsI_15081 GI:218194450 57 2 1 214.5/7.8 71.7/5.4
70. OsC-4107 Hypothetical protein OsJ_00565
[Oryza sativa Japonica Group]
GI:125569214 155 3 23 17.9/5.8 28.9/5.5
71. OsC-4213 Unknown protein [Oryza sativa
Japonica Group]
GI:30017570 60 3 10 35.7/9.6 41.9/5.7
72. OsC-4510 Hypothetical protein OsI_31140 GI:125563499 148 5 12 51.0/5.7 70.4/5.6
73. OsC-4808 Hypothetical protein OsI_02498 GI:218188499 53 2 2 96.9/5.5 102.3/5.6
74. OsC-6105 Hypothetical protein OsJ_22407
[Oryza sativa Japonica Group]
GI:222636106 45 1 3 30.5/7.4 36.2/6.0
75. OsC-6201 Hypothetical protein OsJ_04535 GI:125573095 46 1 8 29.6/10.3 43.2/5.8
76. OsC-6607 PREDICTED: uncharacterized protein
LOC100818188 isoform X1 [Glycine
max]
GI:571541740 45 7 15 43.6/5.4 82.9/5.9
77. OsC-8408 Hypothetical protein OsI_23019 GI:218198209 118 7 23 41.3/6.6 60.9/6.7
78. OsC-8717 Hypothetical protein
SORBIDRAFT_03g034200
GI:115456914 806 13 21 94.9/5.8 95.6/6.5
aSpot number as given on the two dimensional gel images. The first letters (Os) represent the source plant Oryza sativa followed by the fraction cytoplasm (c).bThe significance score (P < 0.05) of a protein, as produced by the Mascot algorithm.cProtein expression profile represents the average change in spot density among different stages of drought stress. The data were taken in terms of fold expression with respect to the
control value and were log transformed to the base 2 in order to level the scale of expression and to reduce the noise.
Frontiers in Plant Science | www.frontiersin.org 10 September 2016 | Volume 7 | Article 1466
Agrawal et al. Drought Affected Rice Root Proteome
Bioenergy &
Metabolism
29%
Cell Defence
& Rescue
22%
Protein
biogenesis &
storage
21%
Cell
signalling
9%
Miscellenious
19% B
P
C
C
P
P
C
M
C
FIGURE 4 | Distribution of 78 identified differentially expressed rice
root proteins in five functional classes based on their putative
functions assigned them using protein database.
mitochondrial-like (OsC-7524), succinate dehydrogenase[ubiquinone] flavoprotein subunit, mitochondrial (OsC-1007),probable nucleoredoxin (OsC-2617), fructokinase-2 (OsC-3314),spermidine synthase 1(OsC-4302), pyruvate dehydrogenaseE1 component (OsC-4803), lecithin:cholesterolacyl transferasefamily protein (OsC-2107), and enolase (OsC-5504) weredownregulated. Other proteins like sucrose synthase 1-likeisoform X2 (OsC-7711), putative transaldolase (OsC-3409),pyrophosphate-fructose 6-phosphate 1-phosphotransferase(OsC-7510) and carbonic anhydrase (OsC-8108) showed mixedpattern of expression.
Alteration in bioenergy metabolism under drought can beaddressed by increase in abundance of enzymes involved inacetyl CoA synthesis, TCA cycle, transport proteins, synthesisof osmoprotectants, and protein metabolism. Acetyl-CoAsynthetase (OsC-8316) showed a 43 fold increase in abundanceon first day and 66 fold on the third day of drought induction(Table 1); however the expression remains similar as control onthe 7th day. The enzyme catalyzes the formation of acetyl CoAfrom acetate utilizing inorganic phosphate (Lin and Oliver, 2008)and can be employed in the TCA cycle during aerobic respirationfor energy production and electron carriers. It is to be notedhere that drought stress resulted in downregulation of pyruvatedehydrogenase (OsC-3417), hence the upregulation of Acetyl-CoA synthetase possibly provided an alternate source of energyproduction (Schwer et al., 2006).
The TCA cycle enzymes like succinyl co-A ligase (OsC-3218),succinate dehydrogenase (OsC-1007) and malate dehydrogenase(OsC-7524) were upregulated. Succinyl-CoA ligase catalyzesthe reversible inter conversion of succinyl-CoA to succinate(Cavalcanti et al., 2014). There was a 20 fold increase in theenzyme abundance at first day (Table 1). The result corroborateswith earlier studies which reported increase in enzyme expressionin wheat roots in response to aluminum stress (Drummondet al., 2001). Succinate dehydrogenase catalyzes the oxidationof succinate to form fumarate (Singer et al., 1973; Jardim-Messeder et al., 2015). The enzyme is an important componentof TCA cycle as well as the electron transport chain (Popova
et al., 2010; Huang and Millar, 2013). Two isozymes of succinatedehydrogenase were found as differential spots, of which one(OsC-1007) was downregulated while the other (OsC-8616) wasupregulated showing 28 fold increase in activity on the thirdday (Table 1). The transcription of succinate dehydrogenaseflavoprotein subunit has been found to be up-regulated in Ilexparaguariensis leaves in response to water deficit and ABAapplication (Acevedo et al., 2013). Furthermore, its increasedactivity results in root elongation (Huang et al., 2012), andalso has a role in plant adaptation toward stress and especiallyby combating reactive oxygen species (Pastore et al., 2007).Of the two isozymes of malate dehydrogenase identified, onewas downregulated (OsC-5107), while the other isozyme (OsC-7524) showed 5 and 11 fold increase in activity on third andseventh day of drought stress, respectively (Table 1). ATPasesare integral transport proteins that help in the hydrolysisof ATP as well as movement of protons across membranesto generate electrochemical gradients (Palmgren and Harper,1999). The action of ATPase can influence stress mechanism bychanging the membrane potential and proton gradient (Elmoreand Coaker, 2011). Transcript and protein levels of ATPasehave been reported to increase under salt stress conditions(Batelli et al., 2007). Interestingly V-type proton ATPase subunitB 1-like protein (OsC-3618) and ATP synthase subunit dmitochondrial-like isoform X2 (OsC-3113) were upregulated.While, approximately 8 fold increases in activity of V-type protonATPase was noted on third and seventh day; the activity ofmitochondrial ATP synthase increased to approximately 33 foldin first day and 58 fold in third day (Table 1).
Mitochondrial pyruvate dehydrogenase (OsC-4803) is anintermediate enzyme that links glycolysis to the citric acid cycle.The enzyme was downregulated under drought stress and resultssubstantiate the earlier reports (Simova-Stoilova et al., 2015).Fructokinase-2 (OsC-3314) is mainly involved in carbohydratemetabolism, more specifically, sucrose and fructose metabolism(Odanaka et al., 2002). Previously rice fructokinase showedelevated levels of expression during infection with the fungus,Magnapor theoryzae (Ryu et al., 2006). Sucrose synthase (OsC-7711), a major protein of energy metabolism, plays an importantrole in controlling the mobilization of sucrose into variouspathways essential for various metabolic, structural, and otherstorage functions of the plant cell (Hesse and Willmitzer, 1996).The protein expression followed a miscellaneous pattern with 2.7fold abundance on third day (Table 1), indicating an alteration inenzyme function more toward osmotic adjustment and storagecompounds (Hasibeder et al., 2015).
Cell Defense and Rescue (CDR)Seventeen differentially expressed proteins related to celldefense and rescue were identified, out of which pathogenesis-related protein 1 (PR1, OsC-2313), one isoform of TPA:class III peroxidase 72 precursor (OsC-4416), glyoxalase I(OsC-4315), putative r40c1 protein (OsC-7314), stress-inducedprotein sti1 (OsC-7617, OsC-8601), and monodehydroascorbatereductase (OsC-3418) were upregulated whereas root specificpathogenesis-related protein 10 (OsC-2212), 24.1 kDa heatshock protein (OsC-4112), isozyme of monodehydroascorbate
Frontiers in Plant Science | www.frontiersin.org 11 September 2016 | Volume 7 | Article 1466
Agrawal et al. Drought Affected Rice Root Proteome
reductase (OsC-5409), TPA: class III peroxidase 72 precursor(OsC-6412 and OsC-3416), isozyme of stress-induced proteinsti1 (OsC-7611) and heat shock protein STI (OsC-8607) weredownregulated, while heat shock protein (HSP)-90 (OsC-3703)and GSH-dependent dehydroascorbatereductase 1 (OsC-7111)showed miscellaneous pattern of expression.
The ROS scavengers provide cells with an efficient machineryfor detoxifying O•−
2 and H2O2 and constitute the first lineof defense. A considerable upregulation of these enzymes wasa major tolerance mechanism adopted by the variety. Theplant class III peroxidases are known to basically involvedin metabolism of ROS along with some other function likeauxin metabolism, lignin and suberin formation and cross-linking of cell wall components (Almagro et al., 2009). Theclass III peroxidase (OsC-4416) expression gradually increasedup to 11 fold from first to seventh day (Table 1). Since theprotein is involved in lignin biosynthesis, a gradual increasein expression suggests enhanced lignin production. Increasedamount of lignin builds the mechanical strength of cell walland thereby protects roots against the dry soil. In addition,cell wall modification is also used to minimize water loss andcell dehydration, thus helping plants to resist and recover fromdrought upon availability of water (Yoshimura et al., 2008).An 11 fold increase in expression was observed in glyoxalase1 (OsC-4315), which is a component of the glyoxalase systemthat carries out the detoxification of methylglyoxal and the otherreactive aldehydes produced as a normal part of metabolism(Vander, 1989). Glyoxalase 1 has been identified as a saltinduced proteome in rice roots and plays an important role inmethylglyoxal detoxification in salt tolerant rice, thus ensuringredox homeostasis (El-Shabrawi et al., 2010). Putative r40c1protein showed increased abundance in response to drought.The exact function of the protein is unknown; however itsupregulation has been correlated to drought tolerance (Kumaret al., 2015). The expression of monodehydroascorbate reductase(OsC-5409) was increased by 9 fold on seventh day. It is acritical component of the ROS scavenging system in rice, andits expression can improve drought tolerance in rice (Wanget al., 2005). Interestingly, the expression of pathogenesis relatedprotein-1 (PR1) was 22 fold increased within 24 h of droughtinduction; however root specific pathogenesis related protein-10(PR10) was down regulated in response to drought. PR1 geneexpression is salicylic-acid responsive and induced in responseto a variety of pathogens (Elvira et al., 2008) while PR-10protein expression is induced through activation of jasmonicacid signaling pathway (Hashimoto et al., 2004). Thus, it mightbe suggested that drought stress might responsible in increaseexpression of SA regulated PR proteins and decreases in JA-regulated defense.
Protein Biogenesis and Storage (PBS)A total of 15 proteins were found as differential spots, out ofwhich nuclear pre-mRNA domain-containing protein 1B-likeisoform X2 (OsC-2220), alpha-tubulin (OsC-3516), putativeelongation factor 2 (OsC-7712), putative thiamine biosynthesisprotein (OsC-3318), putative beta-alanine synthase (OsC-7409),and cysteine synthase (OsC-4207, OsC-4207) was upregulated;
tricin synthase 1 (OsC-4607), putative triose phosphateisomerase (OsC-4202), proteasome subunit alpha type-1 (OsC-5205) was downregulated; while three isozymes of 5-methyltetra hydropteroyl triglutamate-homocysteine methyltransferase(OsC-7609, OsC-7614, OsC-8703), elongation factor 1 beta(OsC-1211), proteasome subunit beta type 3 (OsC-3111) andubiquitin-activating enzyme E1 2, putative (OsC-4815) showedmiscellaneous pattern of expression.
Cysteine synthase is a key enzyme for mediating abiotictolerance owing to the production of antioxidants andmetal chelators, such as glutathione, metallothionein, andphytochelatin. In an earlier study cysteine synthase was foundto be increased in the roots of wheat during stress (Yang et al.,2007). We also observed significantly upregulation of cysteinesynthase (OsC-4207) in tolerant variety used for analysis.Tubulin class of proteins is known to play a critical role in celldivision and elongation and alpha-tubulin (OsC-3516) has beenspecifically reported in roots of drought tolerant rice verities(Rabello et al., 2008). Elongation factor 2 (OsC-7712) whichshowed 24 fold increases in abundance is reportedly a stressresponsive protein which is mostly downregulated in sensitiveverities of rice (Chen et al., 2015). There was a decrease inabundance of proteasome subunit beta type 3 protein (OsC-5205) in response to drought while ubiquitin activating enzyme(OsC-4815), required for ubiquitination, showed a gradualincrease in abundance and increased 2-fold by the end of seventhday (Table 1). The protein ubiquitination plays a central rolein regulating the transcriptional changes required for adaptionand survival strategies of plants to different environmentalstresses (Dametto et al., 2015). Glutamate dehydrogenaseexpedites amino acid synthesis in drought stress conditionswhere nitrogen assimilation is decreased and also known tomaintain homeostasis during stress (Masclaux-Daubresse et al.,2010). Present study showed a 4 fold increase in glutamatedehydrogenase (OsC-8412) expression on seventh day (Table 1).Putative thiamine biosynthesis protein (OsC-3318) showedan 18 fold increase on first day and 27 fold increase on theseventh day. Recent reports suggest that vitamin B1 (thiamine)participates in the processes underlying plant adaptations tovarious types stress conditions including cold, heat, drought,and other oxidative stress (Rapala-Kozik et al., 2012). A gradualincrease in putative beta-alanine synthase (OsC-7409) suggestssignificant enrichment of beta-alanine which functions as anosmoprotectant under drought conditions in various plantsystems (Ranjan et al., 2012).
Cell SignalingThe response of roots to water limiting conditions seems to becrucial to trigger drought tolerance mechanisms, since roots areone of the primary sites for stress signal perception in whicha signaling mechanism initiates a cascade of gene expressionresponses to drought. These changes in protein profile canresult in successful adaptations leading to stress tolerance byregulating protein expression and signal transduction in thestress response (regulatory proteins) or directly protecting theplant against environmental stress (functional proteins). Sevencell signaling related proteins were identified as differential
Frontiers in Plant Science | www.frontiersin.org 12 September 2016 | Volume 7 | Article 1466
Agrawal et al. Drought Affected Rice Root Proteome
spots, out of which actin (OsC-2218) and actin-1-like isoformX2 (OsC-2221) were significantly upregulated. Whilst a 20-foldincrease in expression of actin-depolymerizing factor 3 (OsC-2214) was observed on third day (Table 1), a consistent decreasein the expression of actin-depolymerizing factor 4 (OsC-5105)was noticed during the study. Other proteins like the 14-3-3-like protein (OsC-2201), putative membrane protein (OsC-4203), and prolyl endopeptidase-like protein (OsC-5711) showedmiscellaneous pattern of expression in response to drought.Actin depolymerizing factors (ADFs) are small actin-bindingproteins and are known to confer abiotic drought tolerance(Huang et al., 2012). There was an increase in abundance ofactin depolymerizing factor 3 (OsC-2214) by 8 and 20 fold,
respectively, on first and third day of drought induction, howeverthe expression came down as similar to control at seventh day,thus suggesting its regulatory role in drought tolerance. Duringsignal transduction, 14-3-3 proteins mediate several protein-protein interactions through post-translational modificationlike phosphorylation which ultimately affects multiple plantfunctions (Chen et al., 2006). Surprisingly, a non-coherentexpression of 14-3-3 protein (OsC-2201) was observed underdrought stress. A decrease in prolyl endopeptidase (OsC-5711) expression was observed on first of drought induction;however an increase in its abundance was noted at thirdday and the abundance increased by 3.7-fold at the end ofseventh day.
8 proteins
7 proteins 3 proteins 5 proteins 25 proteins
1 protein2 proteins 7 proteins
5 proteins 5 proteins 10 proteins
BEM
PBS
CDR
SIG
MIS
Cluster 1
N=8
Cluster 4
N=7
Cluster 5
N=7
Cluster 7
N=5
Cluster 8
N=25
Cluster 9
N=5
Cluster 10
N=5
Cluster 11
N=10
FIGURE 5 | Expression clustering of 78 differentially expressed proteins showing 8 clusters based on their expression profiles. The gray lines represent
expression profile of protein separately and the mean expression profile is indicated by pink line. Total number of proteins with same expression profile is provided in
inset within the respective cluster. Detailed information on proteins within each cluster present in Figure S1.
Frontiers in Plant Science | www.frontiersin.org 13 September 2016 | Volume 7 | Article 1466
Agrawal et al. Drought Affected Rice Root Proteome
FIGURE 6 | Quantitative validation of relative expression of selected transcripts related to different functional classes of rice root. Heat map of some
differentially regulated known and hypothetical protein showing their relative expression by qRT-PCR along with their protein expression data obtained from 2DE. The
signal intensity of each transcript was normalized using ubiquitin house keeping gene. The top most bar represents the fold change expression values of selected
genes.
Control
Day-1Day-3
Day-7
1
2
5
3 7
6
110
1
11
4
4
20
3
FIGURE 7 | Venn diagram showing the distribution of differentially
expressed proteins in time-specific and overlapping manner during
drought stress. The areas shown in the diagram are not proportional to the
number of proteins in the groups.
MiscellaneousThe proteins with functions other than the defined categorieswere ascribed as miscellaneous proteins and constituted22% of the differentially expressed protein. Amongst thedifferentially expressed miscellaneous proteins, Sgt1 (OsC-2516),submergence induced, nickel-binding protein 2A (OsC-3115),Oryzain alpha chain; Flags: Precursor [O. sativa Japonica group](OsC-1215), hypothetical protein (OsC-3413), hypotheticalprotein OsI_15081 (OsC-3518), and hypothetical protein
SORBIDRAFT_03g034200 (OsC-8717) were upregulated;hypothetical protein OsI_37864 (OsC-3312) was downregulated,while 3′-N-debenzoyl-2′-deoxytaxol N-benzoyltransferase (OsC-3509) and uncharacterized protein LOC100818188 (OsC-6607)showed varied expression.
SGT1 (OsC-2516) is a highly conserved protein amongeukaryotes that binds specifically to the molecular chaperone,HSP90 and helps in regulation of resistance extended by manyresistance proteins (Azevedo et al., 2006; Wang et al., 2010). It isinteresting to note here that there was an upregulation of HSP90(protein abundance increased to 4.6-fold) with SGT1 in responseto drought in Heena roots might involve in drought tolerance.The nickel-binding protein 2A is the member of dehydration-responsive element-binding protein 2A family which acts as atranscriptional activator that binds specifically to the cis-actingdehydration-responsive element (DRE) to regulates high salinityand dehydration-inducible transcription (Dubouzet et al., 2003).A 23-fold increase in protein abundance suggests an importantrole of nickel-binding protein 2A (OsC-3115) in conferringdrought tolerance to Heena, however its exact function asa drought resistant protein needs to be further elucidated.Another hypothetical protein SORBIDRAFT_03g034200(OsC-8717) showed a 63-fold increase in abundance onseventh day which indicates its important role in conferringdrought tolerance in later stage, however it needs furtherelucidation.
Cluster Analysis of Drought ResponsiveProteinThe unbiased hierarchical clustering method, SOTA was usedto study the correlated expression pattern of the droughtresponsive proteins. This method allows integration of themultiple proteins showing similar expression profiles whichprovide comprehensive overview of the rice root protein networkregulation in coordination. The data were log transformed tothe base of 2 to reduce the noise and level the scale of fold-expression for SOTA clustering analysis. The analysis grouped 78proteins into 11 distinct clusters (Figure 5), allowing a maximumdiversity of 0.8 within a single cluster. Figure S1 represents the
Frontiers in Plant Science | www.frontiersin.org 14 September 2016 | Volume 7 | Article 1466
Agrawal et al. Drought Affected Rice Root Proteome
detailed information on proteins within each cluster. Only thoseclusters which is having n ≥ 5 (n represents number of identifiedproteins) were considered for coordinated expression profile ofthe drought responsive proteins. Analysis revealed that cluster 8contains maximum number of proteins with similar expressionfollowed by cluster 11. The proteins of cluster 8, displayeddecreased accumulation during initial stage of drought andsubsequent increase at late stages. While in cluster 1, maximumnumber of the proteins displayed increased accumulation duringinitial stage of drought and then decrease at late stage. Most ofthe clusters except cluster 8 showed selective representation forspecific functional classes of proteins. The proteins in cluster 9,10, and 11 were found to be over-expressed throughout droughtstress in one or all stages, whereas in cluster 5 and 7 showed onlydown-regulated proteins. The cluster 4 which contain 7 proteinswere expressed immediately in day 1 then disappear throughoutsuggested their role in helping in plant to tolerate drought shock.The drought responsive proteins with miscellaneous functionwere found to be distributed in almost all the clusters. SOTAclustering analysis may provide useful information for theirputative function in drought tolerance based on the abundancein rice root.
Validation of Proteomics Data of Selected Proteins
with qRT PCR
The list of putatively regulated proteins depicted in Table 1 is asnapshot of proteins from rice root. We have selected 5 geneswith miscellaneous function for comparing their expression withproteome data. Among them, a highly conserved eukaryoticprotein, SGT1 (OsC-2516) shows its binding specifically to themolecular chaperone to regulation of resistance extended bymany resistance proteins (Wang et al., 2010). An upregulatedTCA cycle enzymes succinyl co-A ligase (OsC-3218) catalyzesthe reversible inter conversion of succinyl-CoA to succinate(Cavalcanti et al., 2014). Triose phosphate isomerase (OsC-4202) is required to produce ATP during glycolysis and ithas been reported to have highly induced by salt or droughtin rice leaves (Lee et al., 2011). Benzoyltransferase (OsC-3509) have role in final step of acylation of taxol biosynthesispathway but its role in drought stress is hitherto unknown.Apart from these, some other unknown genes like hypotheticalprotein (OsC-6201), PREDICTED: uncharacterized protein(OsC-6607), hypothetical protein SORBIDRAFT_03g034200(OsC-8717) were also taken to check their abundance. However,expression profile of these proteins were somewhat similar
TCA
cycle
Glucose-6-P
Glucose
Fructose-6-P
Fructose-1,6-bis-P
GAP DHAP
G6P1
ALD
TPI
1,3-bis-PGA
GAPDH
3-PGA
PGK
2-PGA
PEP
Pyruvate
Citrate
Isocitrate
α-ketoglutarate
Succinyl-CoASuccinate
Fumarate
Malate
Oxaloacetate
AH
IDHFR
PGM
PK
Acetyl-Co-A
L-methionine
SAM
Homocys!ne
SAMS
H2O2
H2O
H2O + O2
O2
HSPs
Cystathione
Cystein
O-acetyl Serine
Serine
SAT
CGL
CBS
SOD
MDAMDA
AA
APx
AA
MDAR
Induced drought
tolerance
Methylglyoxal
Hemiacetal
Lac!c acid
GSSH
GSH
NADPH
NADP+
GR
PFK3314
GLO14315
SD1007 8616
MD5107 7524
PDH3417 4803
SucCoA3218
ENO5504
CoASH8316
CS4207
MDAR5409 7524
POX4416 8506
HSP903703
DHAR7111
MS7609 7614 8703
Drought induced proteins
ACT2218 2221 2221 2516
SGT1
NiBP3115
OAC1215
HYPO3413 3518 4213 4510 6607
STI18601
PR12313
r40C17314
PEP5711
Lignification and
suberization
of cell wall
FIGURE 8 | Illustration of role of differentially regulated proteins involved in different pathways in rice for sustaining during drought stress. Proteins
engaged in carbon breakdown, protein synthesis and ROS pathway are displayed on the corresponding metabolic pathways. Graphs are the representatives of
expression profile of individual protein and number given below in each graph indicates the protein identification number.
Frontiers in Plant Science | www.frontiersin.org 15 September 2016 | Volume 7 | Article 1466
Agrawal et al. Drought Affected Rice Root Proteome
as determined by the 2-DE analysis (Table 1), nevertheless,there was difference in fold-induction in protein expression(Figure 6, Table 1). Since all the differentially expressed proteinswere not identified, it can be assumed that some otherisoforms of the same protein might be part of the rootproteome. The post-translational modifications of some of thedifferentially expressed proteins might also affect the actualexpression levels determined by two different techniques.Further study may explore such type of correlation in abetter way.
CONCLUSION
Crop plants including rice exhibit several adaptive andacclimatization strategies to combat environmental conditionssuch as drought. Such strategies include from visible phenotypicchanges to complex physicochemical traits. At molecularlevel, these traits, often represented by differentially regulatedproteins during stress and may serve as important stresstolerance markers. We report here a systematic proteomicanalysis of the rice root proteins under PEG-simulated droughtstress conditions. In conclusion, the root-specific comparativeproteomes of rice identified a number of proteins that areputatively associated with stage specific drought tolerant. Ofthe 78 differentially expressed proteins, 10 were found to bedifferentially regulated in all the four stages during droughtstress. Three proteins exclusively expressed in the control, weredisappear after drought stresses. Total 8 proteins found to benewly synthesized during early or later stage of drought stress,implying their possible role in drought tolerance (Figure 7).Functional classification revealed that maximum number ofproteins fall in the category of bioenergy and metabolismfollowed by those involved in cell defense (Figure 4). The highernumber of metabolic proteins offers a unique opportunity topredict more activity toward carbon assimilation, respiration,and storage products like sucrose and starch to combat droughtcondition. Based on results from present study, a putative modelelucidating the role of the differentially expressed proteins indifferent pathways like glycolysis, TCA, amino acid metabolism,
ROS, and other possible mechanism(s) underlying dehydrationtolerance is in Figure 8. A large number of proteins withmiscellaneous functions are matter of further investigation fortheir putative role in drought tolerance. However, the presentproteomics study could represent only a small part of the riceproteome, further investigation to assign their putative biologicalfunctions may be useful for a better understanding of complexbiological traits, such as drought tolerance. Many other drought-responsive proteins still need to be identified with advancementto technology which may help in better understanding of thedrought response in rice.
AUTHOR CONTRIBUTIONS
CN, DC, and PC designed the experiment. LA, SG, and SMperformed the experiments. LA, GP, and SK analyze the data andwrote the paper.
FUNDING
The study was supported by New Initiative (as a Cross FlowTechnology project) “Root Biology and its Correlation toSustainable Plant Development and Soil Fertility (BSC0204)”from Council of Scientific and Industrial Research (CSIR), NewDelhi, India and JC Bose Fellowship, Science and EngineeringResearch Board, Department of Science and Technology,Government of India, awarded to CN.
SUPPLEMENTARY MATERIAL
The Supplementary Material for this article can be foundonline at: http://journal.frontiersin.org/article/10.3389/fpls.2016.01466
Figure S1 | Heat map revealing the expression of differentially regulated
proteins of rice root present within each cluster.
Table S1 | List of primers used for Real Time PCR.
Table S2 | Reproducibility of 2-D gels.
Table S3 | List of rice root proteins identified by MS/MS analysis.
REFERENCES
Acevedo, R. M., Maiale, S. J., Pessino, S. C., Bottini, R., Ruiz, O. A., andSansberro, P. A. (2013). A succinate dehydrogenase flavoprotein subunit-like transcript is up-regulated in Ilex paraguariensis leaves in response towater deficit and abscisic acid. Plant Physiol. Biochem. 65, 48–54. doi:10.1016/j.plaphy.2012.12.016
Agrawal, L., Chakraborty, S., Jaiswal, D. K., Gupta, S., Datta, A., and Chakraborty,N. (2008). Comparative proteomics of tuber induction, development andmaturation reveal the complexity of tuberization process in potato (Solanumtuberosum L.). J. Proteome Res. 7, 3803–3817. doi: 10.1021/pr8000755
Agrawal, L., Narula, K., Basu, S., Shekhar, S., Ghosh, S., Datta, A., et al. (2013).Comparative proteomics reveals a role for seed storage protein AmA1 incellular growth, development, and nutrient accumulation. J. Proteome Res. 12,4904–4930. doi: 10.1021/pr4007987
Ahsan, N., Lee, D. G., Alam, I., Kim, P. J., Lee, J. J., Ahn, Y. O., et al. (2008).Comparative proteomic study of arsenic-induced differentially expressed
proteins in rice roots reveals glutathione plays a central role during As stress.Proteomics 8, 3561–3576. doi: 10.1002/pmic.200701189
Ali, G. M., and Komatsu, S. (2006). Proteomic analysis of rice leaf sheath duringdrought stress. J. Proteome Res. 5, 396–403. doi: 10.1021/pr050291g
Almagro, L., Gómez Ros, L. V., Belchi-Navarro, S., Bru, R., Ros Barceló, A., andPedreño, M. A. (2009). Class III peroxidases in plant defence reactions. J. Exp.Bot. 60, 377–390. doi: 10.1093/jxb/ern277
Atkinson, N. J., and Urwin, P. E. (2012). The interaction of plant biotic andabiotic stresses: from genes to the field. J. Exp. Bot. 63, 3523–3543. doi:10.1093/jxb/ers100
Azevedo, C., Betsuyaku, S., Peart, J., Takahashi, A., Noël, L., Sadanandom, A., et al.(2006). Role of SGT1 in resistance protein accumulation in plant immunity.EMBO J. 25, 2007–2016. doi: 10.1038/sj.emboj.7601084
Batelli, G., Verslues, P. E., Agius, F., Qiu, Q., Fujii, H., Pan, S., et al. (2007). SOS2promotes salt tolerance in part by interacting with the vacuolar H+-ATPaseand upregulating its transport activity. Mol. Cell. Biol. 27, 7781–7790. doi:10.1128/MCB.00430-07
Frontiers in Plant Science | www.frontiersin.org 16 September 2016 | Volume 7 | Article 1466
Agrawal et al. Drought Affected Rice Root Proteome
Cavalcanti, J. H. F., Esteves-Ferreira, A. A., Quinhones, C. G. S., Pereira-Lima,I. A., Nunes-Nesi, A., Fernie, A. R., et al. (2014). Evolution and functionalimplications of the Tricarboxylic Acid Cycle as revealed by phylogeneticanalysis. Genome Biol. Evol. 6, 2830–2848. doi: 10.1093/gbe/evu221
Chen, C., Song, Y., Zhuang, K., Li, L., Xia, Y., and Shen, Z. (2015). Proteomicanalysis of copper-binding proteins in excess copper-stressed roots of two rice(Oryza sativa L.) varieties with different Cu tolerances. PLoS ONE 10:e0125367.doi: 10.1371/journal.pone.0125367
Chen, F., Li, Q., Sun, L., and He, Z. (2006). The Rice 14-3-3 gene family and itsinvolvement in responses to biotic and abiotic stress. DNA Res. 13, 53–63. doi:10.1093/dnares/dsl001
Dametto, A., Buffon, G., Dos Reis Blasi, É. A., and Sperotto, R. A. (2015).Ubiquitination pathway as a target to develop abiotic stress tolerance in rice.Plant Signal. Behav. 10:e1057369. doi: 10.1080/15592324.2015.1057369
Davies, W. J., and Zhang, J. (1991). Root signals and the regulation of growth anddevelopment of plant in dry soil. Ann. Rev. Plant Physiol. Plant Mol. Biol. 42,55–76. doi: 10.1146/annurev.pp.42.060191.000415
Drummond, R. D., Guimaraes, C. T., Felix, J., Ninamango-Cardenas, F. E.,Carneiro, N. P., Paiva, E., et al. (2001). Prospecting sugarcane genes involvedin aluminum tolerance. Genet. Mol. Biol. 24, 221–230. doi: 10.1590/S1415-47572001000100029
Dubouzet, J. G., Sakuma, Y., Ito, Y., Kasuga, M., Dubouzet, E. G., Miura, S., et al.(2003). OsDREB genes in rice, Oryza sativa L., encode transcription activatorsthat function in drought-, high-salt- and cold-responsive gene expression. PlantJ. 33, 751–763. doi: 10.1046/j.1365-313X.2003.01661.x
Elmore, J. M., and Coaker, G. (2011). The role of the plasma membraneH+-ATPase in plant-microbe interactions. Mol. Plant 4, 416–427. doi:10.1093/mp/ssq083
El-Shabrawi, H., Kumar, B., Kaul, T., Reddy, M. K., Singla-Pareek, S. L.,and Sopory, S. K. (2010). Redox homeostasis, antioxidant defense, andmethylglyoxal detoxification as markers for salt tolerance in Pokkali rice.Protoplasma 245, 85–96. doi: 10.1007/s00709-010-0144-6
Elvira, M. I., Galdeano, M. M., Gilardi, P., García-Luque, I., and Serra, M. T.(2008). Proteomic analysis of pathogenesis-related proteins (PRs) inducedby compatible and incompatible interactions of pepper mild mottle virus(PMMoV) in Capsicum chinense L3 plants. J. Exp. Bot. 59, 1253–1265. doi:10.1093/jxb/ern032
Hashimoto, M., Kisseleva, L., Sawa, S., Furukawa, T., Komatsu, S., and Koshiba,T. (2004). A novel rice PR10 protein, RSOsPR10, specifically induced in rootsby biotic and abiotic stresses, possibly via the jasmonic acid signaling pathway.Plant Cell Physiol. 45, 550–559. doi: 10.1093/pcp/pch063
Hasibeder, R., Fuchslueger, L., Richter, A., and Bahn, M. (2015). Summer droughtalters carbon allocation to roots and root respiration in mountain grassland.New Phytol. 205, 1117–1127. doi: 10.1111/nph.13146
Herrero, J., Valencia, A., and Dopazo, J. (2001). A hierarchical unsupervisedgrowing neural network for clustering gene expression patterns. Bioinformatics
17, 126–136. doi: 10.1093/bioinformatics/17.2.126Hesse, H., and Willmitzer, L. (1996). Expression analysis of a sucrose synthase
gene from sugar beet (Beta vulgaris L.). Plant Mol. Biol. 30, 863–872. doi:10.1007/BF00020799
Huang, S., and Millar, A. H. (2013). Succinate dehydrogenase: the complexroles of a simple enzyme. Curr. Opin. Plant Biol. 16, 344–349. doi:10.1016/j.pbi.2013.02.007
Huang, S., Taylor, N. L., Ströher, E., Fenske, R. A., and Millar, R. A. (2012).Succinate dehydrogenase assembly factor 2 is needed for assembly and activityof mitochondrial complex II and for normal root elongation in Arabidopsis.Plant J. 73, 429–441. doi: 10.1111/tpj.12041
Jardim-Messeder, D., Caverzan, A., Rauber, R., de Souza, F. E., Margis-Pinheiro,M., and Galina, A. (2015). Succinate dehydrogenase (mitochondrial complexII) is a source of reactive oxygen species in plants and regulates developmentand stress responses. New Phytol. 208, 776–789. doi: 10.1111/nph.13515
Kawasaki, S., Borchert, C., Deyholos, M., Wang, H., Brazille, S., Kawai, K., et al.(2001). Gene expression profiles during the initial phase of salt stress in rice.Plant Cell 13, 889–905. doi: 10.1105/tpc.13.4.889
Khan, M., and Komatsu, S. (2004). Rice proteomics: recent developmentsand analysis of nuclear proteins. Phytochemistry 65, 1671–1681. doi:10.1016/j.phytochem.2004.04.012
Komatsu, S. (2005). Rice Proteome Database: a step toward functional analysis ofthe rice genome. Plant Mol. Biol. 59, 179–190. doi: 10.1007/s11103-005-2160-z
Kumar, A., Bimolata, W., Kannan, M., Kirti, P. B., Qureshi, I. A., and Ghazi, I. A.(2015). Comparative proteomics reveals differential induction of both bioticand abiotic stress response associated proteins in rice during Xanthomonas
oryzae pv. oryzae infection. Funct. Integr. Genomics 15, 425–437. doi:10.1007/s10142-014-0431-y
Lafitte, H. R., Yongsheng, G., Yan, S., and Li, Z. K. (2007). Whole plant responses,key processes, and adaptation to drought stress: the case of rice. J. Exp. Bot. 58,169–175. doi: 10.1093/jxb/erl101
Lee, D. G., Park, K. W., An, J. Y., Sohn, Y. G., Ha, J. K., Kim, H. Y., et al. (2011).Proteomics analysis of salt-induced leaf proteins in two rice germplasms withdifferent salt sensitivity. Can. J. Plant Sci. 91, 337–349. doi: 10.4141/CJPS10022
Lin, M., and Oliver, D. J. (2008). The role of acetyl-coenzyme A synthetase inArabidopsis. Plant Physiol. 147, 1822–1829. doi: 10.1104/pp.108.121269
Masclaux-Daubresse, C., Daniel-Vedele, F., Dechorgnat, J., Chardon, F.,Gaufichon, L., and Suzuki, A. (2010). Nitrogen uptake, assimilation andremobilization in plants: challenges for sustainable and productive agriculture.Ann. Bot. 105, 1141–1157. doi: 10.1093/aob/mcq028
Moumeni, A., Satoh, K., Kondoh, H., Asano, T., Hosaka, A., Venuprasad, R.,et al. (2011). Comparative analysis of root transcriptome profiles of two pairsof drought-tolerant and susceptible rice near-isogenic lines under differentdrought stress. BMC Plant Biol. 11:174. doi: 10.1186/1471-2229-11-174
Muscolo, A., Sidari, M., Anastasi, U., Santonoceto, C., and Maggio, A. (2014).Effect of PEG-induced drought stress on seed germination of four lentilgenotypes. J. Plant Interact. 9, 354–363. doi: 10.1080/17429145.2013.835880
Narsai, R., Devenish, J., Castleden, I., Narsai, K., Xu, L., Shou, H., et al. (2013).Rice DB: an Oryza information portal linking annotation, subcellular location,function, expression, regulation, and evolutionary information for rice andArabidopsis. Plant J. 76, 1057–1073. doi: 10.1111/tpj.12357
Odanaka, S., Bennett, A. B., and Kanayama, Y. (2002). Distinct physiologicalroles of fructokinase isozymes revealed by gene-specific suppression ofFrk1 and Frk2 expression in tomato. Plant Physiol. 129, 1119–1126. doi:10.1104/pp.000703
Palmgren, M. G., and Harper, J. F. (1999). Pumping with plant P-type ATPases. J.Exp. Bot. 50, 883–893. doi: 10.1093/jxb/50.Special_Issue.883
Pandey, S., and Bhandari, H. (2008). “Drought: economics costs and researchimplications,” in Drought Frontiers in Rice: Crop Improvement for Increased
Rainfed Production, eds R. Serraj, J. Bennett, and B. Hardy (Singapore: IRRI;World Scientific Publishing), 3–17.
Pastore, D., Trono, D., Laus, M. N., di Fonzo, N., and Flagella, Z. (2007).Possible plant mitochondria involvement in cell adaptation to drought stress.A case study: durum wheat mitochondria. J. Exp. Bot. 58, 195–210. doi:10.1093/jxb/erl273
Périn, C., Rebouillat, J., Brasileiro, A. C. M., Diévart, A., Gantet, P., Breitler,J. C., et al. (2007). “Novel insights into the genomics of rice root adaptivedevelopment,” in Rice Genetics, eds D. S. Brar, D. J. Mackill, B. Hardy (Londres;Singapore: World Scientific; IRRI), 117–141.
Popova, V. N., Eprintseva, A. T., Fedorina, D. N., and Igamberdiev, A. U. (2010).Succinate dehydrogenase in Arabidopsis thaliana is regulated by light viaphytochrome A. FEBS Lett. 584, 199–202. doi: 10.1016/j.febslet.2009.11.057
Price, A. H., Cairns, J. E., Horton, P., Jones, H. G., and Griffiths, H. (2002). Linkingdrought-resistance mechanisms to drought avoidance in upland rice usinga QTL approach: progress and new opportunities to integrate stomatal andmesophyll responses. J. Exp. Bot. 53, 989–1004. doi: 10.1093/jexbot/53.371.989
Rabello, A. R., Guimarães, C. M., Rangel, P. H. N., da Silva, F. R., Seixas, D., deSouza, E., et al. (2008). Identification of drought-responsive genes in roots ofupland rice (Oryza sativa L). BMC Genomics 9:485. doi: 10.1186/1471-2164-9-485
Ranjan, A., Pandey, N., Lakhwani, D., Dubey, N. K., Pathre, U. V., and Sawant, S. V.(2012). Comparative transcriptomic analysis of roots of contrasting Gossypiumherbaceum genotypes revealing adaptation to drought. BMC Genomics 13:680.doi: 10.1186/1471-2164-13-680
Rapala-Kozik,M.,Wolak, N., Kudja,M., and Banas, A. K. (2012). The upregulationof thiamine (vitamin B1) biosynthesis in Arabidopsis thaliana seedlings undersalt and osmotic stress conditions is mediated by abscisic acid at the early stagesof this stress response. BMC Plant Biol. 12:2. doi: 10.1186/1471-2229-12-2
Frontiers in Plant Science | www.frontiersin.org 17 September 2016 | Volume 7 | Article 1466
Agrawal et al. Drought Affected Rice Root Proteome
Romijin, E. P., Christis, C., Wieffer, M., Gouw, W. J., Fullaondo, A., Sluijs, P.,et al. (2005). Expression clustering reveals detailed co-expression patterns offunctionally related proteins during b cell differentiation. A proteomic studyusing a combination of one-dimensional gel electrophoresis, LC-MS/MS, andstable isotope labeling by amino acids in cell culture (SILAC). Mol. Cell
Proteomics 4, 1297–1310. doi: 10.1074/mcp.M500123-MCP200Ryu, H. S., Han, M., Lee, S. K., Cho, J. I., Ryoo, N., Heu, S., et al. (2006).
A comprehensive expression analysis of the WRKY gene superfamily inrice plants during defense response. Plant Cell Rep. 25, 836–847. doi:10.1007/s00299-006-0138-1
Schwer, B., Bunkenborg, J., Verdin, R., Andersen, J., and Verdin, E. (2006).Reversible lysine acetylation controls the activity of the mitochondrial enzymeacetyl-CoA synthetase. Proc. Natl. Acad. Sci. U.S.A. 103, 10224–10229. doi:10.1073/pnas.0603968103
Sengupta, D., Kannan, M., and Reddy, A. R. (2011). A root proteomics-basedinsight reveals dynamic regulation of root proteins under progressive droughtstress and recovery in Vigna radiata (L.) Wilczek. Planta 233, 1111–1127. doi:10.1007/s00425-011-1365-4
Simova-Stoilova, L. P., Romero-Rodríguez, M. C., Sánchez-Lucas, R., Navarro-Cerrillo, R. M. J., Medina-Aunon, J. A., and Jorrín-Novo, J. V. (2015). 2-DEproteomics analysis of drought treated seedlings ofQuercus ilex supports a rootactive strategy for metabolic adaptation in response to water shortage. Front.Plant Sci. 6:627. doi: 10.3389/fpls.2015.00627
Singer, T. P., Oestreicher, G., and Hogue, P. (1973). Regulation of succinatedehyrogenase in higher plants: I. Some general characteristics of themembrane-bound enzyme. Plant Physiol. 52, 616–621. doi: 10.1104/pp.52.6.616
Sreenivasulu, N., Sopory, S. K., and Kavi Kishor, P. B. (2007). Decipheringthe regulatory mechanisms of abiotic stress tolerance in plants by genomicapproaches. Gene 388, 1–13. doi: 10.1016/j.gene.2006.10.009
Subba, P., Kumar, R., Gayali, S., Shekhar, S., Parveen, S., Pandey, A., et al. (2013).Characterisation of the nuclear proteome of a dehydration-sensitive cultivarof chickpea and comparative proteomic analysis with a tolerant cultivar.Proteomics 13, 1973–1992. doi: 10.1002/pmic.201200380
Tuberosa, R., and Salvi, S. (2006). Genomics-based approaches to improvedrought tolerance of crops. Trends Plant Sci. 11, 405–412. doi:10.1016/j.tplants.2006.06.003
vanWijk, K. J. (2001). Challenges and prospects of plant proteomics. Plant Physiol.126, 501–508. doi: 10.1104/pp.126.2.501
Vander, J. D. (1989). “The glyoxalase system,” in Glutathione: Chemical,
Biochemical, and Medical Aspects, eds D. Dolphin, R. Poulson, and O.Avramovic (New York, NY: John Wiley & Sons Ltd.), 597–641.
Wang, F. Z., Wang, Q. B., Kwon, S. Y., Kwak, S. S., and Su, W. A.(2005). Enhanced drought tolerance of transgenic rice plants expressing apea manganese superoxide dismutase. J. Plant Physiol. 162, 465–472. doi:10.1016/j.jplph.2004.09.009
Wang, K., Uppalapati, S. R., Zhu, X., Dinesh-Kumar, S. P., andMysore, K. S. (2010).SGT1 positively regulates the process of plant cell death during both compatibleand incompatible plant-pathogen interactions.Mol. Plant Pathol. 11, 597–611.doi: 10.1111/j.1364-3703.2010.00631.x
Yang, Q., Wang, Y., Zhang, J., Shi, W., Qian, C., and Peng, X. (2007). Identificationof aluminum-responsive proteins in rice roots by a proteomic approach:cysteine synthase as a key player in Al response. Proteomics 7, 737–749. doi:10.1002/pmic.200600703
Yoshimura, K., Masuda, A., Kuwano, M., Yokota, A., and Akashi, K. (2008).Programmed proteome response for drought avoidance/tolerance in the rootof a C3 xerophyte (wild watermelon) under water deficits. Plant Cell Physiol.49, 226–241. doi: 10.1093/pcp/pcm180
Zivy, M., and de Vienne, D. (2000). Proteomics: a link between genomics, geneticsand physiology. Plant Mol. Biol. 44, 575–580. doi: 10.1023/A:1026525406953
Conflict of Interest Statement: The authors declare that the research wasconducted in the absence of any commercial or financial relationships that couldbe construed as a potential conflict of interest.
Copyright © 2016 Agrawal, Gupta, Mishra, Pandey, Kumar, Chauhan,
Chakrabarty and Nautiyal. This is an open-access article distributed under
the terms of the Creative Commons Attribution License (CC BY). The use,
distribution or reproduction in other forums is permitted, provided the original
author(s) or licensor are credited and that the original publication in this
journal is cited, in accordance with accepted academic practice. No use,
distribution or reproduction is permitted which does not comply with these
terms.
Frontiers in Plant Science | www.frontiersin.org 18 September 2016 | Volume 7 | Article 1466