ORIGINAL RESEARCH article

Front. Plant Sci., 11 September 2015

Sec. Plant Proteomics and Protein Structural Biology

Volume 6 - 2015 | https://doi.org/10.3389/fpls.2015.00732

Identification of early salt stress responsive proteins in seedling roots of upland cotton (Gossypium hirsutum L.) employing iTRAQ-based proteomic technique

  • 1. College of Life Sciences, Henan University Kaifeng, China

  • 2. Economic Crop Research Institute, Henan Academy of Agricultural Sciences Zhengzhou, China

Abstract

Soil salinity is a major abiotic stress that limits plant growth and agricultural productivity. Upland cotton (Gossypium hirsutum L.) is highly tolerant to salinity; however, large-scale proteomic data of cotton in response to salt stress are still scant. Here, an isobaric tag for relative and absolute quantitation (iTRAQ)-based proteomic technique was employed to identify the early differentially expressed proteins (DEPs) from salt-treated cotton roots. One hundred and twenty-eight DEPs were identified, 76 of which displayed increased abundance and 52 decreased under salt stress conditions. The majority of the proteins have functions related to carbohydrate and energy metabolism, transcription, protein metabolism, cell wall and cytoskeleton metabolism, membrane and transport, signal transduction, in addition to stress and defense. It is worth emphasizing that some novel salt-responsive proteins were identified, which are involved in cell cytoskeleton metabolism (actin-related protein2, ARP2, and fasciclin-like arabinogalactan proteins, FLAs), membrane transport (tonoplast intrinsic proteins, TIPs, and plasma membrane intrinsic proteins, PIPs), signal transduction (leucine-rich repeat receptor-like kinase encoding genes, LRR-RLKs) and stress responses (thaumatin-like protein, TLP, universal stress protein, USP, dirigent-like protein, DIR, desiccation-related protein PCC13-62). High positive correlation between the abundance of some altered proteins (superoxide dismutase, SOD, peroxidase, POD, glutathione S-transferase, GST, monodehydroascorbate reductase, MDAR, and malate dehydrogenase, MDH) and their enzyme activity was evaluated. The results demonstrate that the iTRAQ-based proteomic technique is reliable for identifying and quantifying a large number of cotton root proteins. qRT-PCR was used to study the gene expression levels of the five above-mentioned proteins; four patterns are consistent with those of induced protein. These results showed that the proteome of cotton roots under NaCl stress is complex. The comparative protein profiles of roots under salinity vs control improves the understanding of the molecular mechanisms involved in the tolerance of plants to salt stress. This work provides a good basis for further functional elucidation of these DEPs using genetic and/or other approaches, and, consequently, candidate genes for genetic engineering to improve crop salt tolerance.

Introduction

Soil salinity is one of the most important environmental factors limiting plant growth and productivity throughout the world (Munns, 2002). Excessive Na+ in the soil inhibits the absorption of mineral nutrients and moisture leading to the accumulation of toxic ions in plants. Plants employ several strategies to cope with salt stress. These include regulating the expression of specific proteins for the reestablishment of proper cellular ion and osmotic homeostasis with other concomitant processes of repair and detoxification (Chinnusamy et al., 2005). The salt signal is primarily perceived through roots, which rapidly respond to maintain function and transmit signals to the shoot for appropriate changes in function (Zhao et al., 2013). Salt-tolerance studies in plants provide insights into the molecular and biochemical basis of plant stress tolerance, which ultimately lead to crop improvement.

Upland Cotton (Gossypium hirsutum L.) is one of the most important textile fiber crops. Although cotton has a higher basal level of tolerance to NaCl compared with other major crops (Maas and Hoffman, 1977), its growth, yield and fiber quality are adversely affected, especially at germination and at the young seedling stage (Ahmad et al., 2002). Understanding the molecular mechanism of salt tolerance can provide many candidate genes for genetic engineering. Many salt-resistant genes have been identified in model plants but only a few salt stress-inducible genes, e.g. Na+/H+antiporter (Wu et al., 2004), DREB (Gao et al., 2009), ERF (Champion et al., 2009; Jin et al., 2010), NAC (Meng et al., 2009), GhMT3a (Xue et al., 2009), MPK (Zhang et al., 2011), MKK (Lu et al., 2013), and ZFP (Guo et al., 2009), have been documented in cotton. Recently, with advances in transcriptome mapping (high-throughput sequencing), some salt-responsive genes and molecular regulatory pathways have been identified in cotton (Yao et al., 2011; Wang et al., 2012; Xu et al., 2013). These studies provide relevant information about the stress-responsive genes, but the transcriptome data may not correlate with results from proteomic analysis due to post-transcriptional and post-translational modifications. Therefore, it is necessary to investigate the change of proteins under salt stress conditions to be able to understand the adaptive mechanism of salt tolerance in cotton.

Proteomic analysis is a tool that facilitates the study of global protein expression and provides a large amount of information about the individual proteins involved in specific biological responses. It has been used to analyze salt stress induced alterations in the root proteome of plant species, such as Arabidopsis (Jiang et al., 2007), rice (Chitteti and Peng, 2007; Cheng et al., 2009), barley (Witzel et al., 2009), wheat (Peng et al., 2009; Guo et al., 2012a), maize (Zörb et al., 2010), soybean (Aghaei et al., 2009), tomato (Manaa et al., 2011; Gong et al., 2014), cucumber (Du et al., 2010), and salt cress (Zhou et al., 2010). Over 850 DEPs of salt-stressed roots have been identified in the above-mentioned studies. Many previous studies relied upon 2D gel electrophoresis data; however, it is difficult to identify low abundant proteins, proteins with low (< 15 kDa) or high (>150 kDa) molecular weights, proteins that are excessively acidic or basic as well as hydrophobic proteins (Zieske, 2006). Non-gel-based quantitative proteomics techniques established in recent years have overcome some of the drawbacks of the above-mentioned method. iTRAQ is a mass spectrometry-based proteomics technique that can be used to evaluate cell metabolic differences. Zhu et al. (2009) employed iTRAQ to reveal functional differentiation of Brassica napus guard cells and mesophyll cells. iTRAQ can also be used to investigate plant responses to deficient or excess mineral nutrients. For example, Yang et al. (2013) successfully analyzed the protein profile of Citrus sinensis roots in response to long-term boron-deficiency with iTRAQ. In addition, Fukao et al. (2011) used iTRAQ analysis to reveal mechanisms of growth defects due to excess zinc in Arabidopsis. It can retain information on the post-translational modification (PTM), simultaneously analyze multiple samples, help to quantify proteins not amenable to the 2D gel approach (Wang et al., 2014) and relatively quantify peptides at a global level (Ghosh et al., 2013). Gong et al. (2014) used iTRAQ to identify a set of DEPs in tomato roots exposed to salt and alkali stress. However, large-scale proteomic data of cotton roots in response to salt stress has not been reported in previous studies. In this present study, an iTRAQ-based proteomic technique was used to identify the early DEPs in order to elucidate the effects of salt stress in cotton seedling roots treated with NaCl for 24 h.

Materials and methods

Plant culture and salt treatments

Seeds of ZMS23, a salt-tolerant variety, were obtained from the Institute for Cotton Research of the Chinese Academy of Agriculture Science. Cotton seeds were sterilized with 10% H2O2 for 30 min and rinsed with distilled water. The sterilized seeds were germinated on filter paper soaked in distilled water in Petri dishes at 26°C. After 7 days, 54 uniform germinated seedlings were transferred to six plastic containers (48 × 36 × 15 cm) and each contained nine seedlings, which were filled with Hoagland's solution (5 mM Ca(NO3), 3 mM KNO3, 2 mM MgSO4, 0.5 mM KH2PO4, 2.5 μM FeNa2(EDTA), 2.5 μM H3BO3, 5 μM MnC12, 0.5 μM ZnSO4, 0.3 μM CuSO4, and 0.05 μM (NH4)6MoO24), and renewed every 2 days. The experiment was performed in a climate chamber under the following conditions: 28/23°C day/night temperature, relative humidity of 70–80% and a 14 h light period/day at an intensity of 600 μmol m−2 s−1. When the plants grew to the trefoil stage, three plastic containers (including 27 seedlings) were renewed with Hoagland's solution and 200 mM NaCl was added, but no NaCl was added to the other three containers used as a control for 24 h, respectively. After treatment, the 1–5 cm portions of root tips were harvested and frozen at −80°C. In the same way, another biological repeat was carried out.

Protein extraction

Cotton roots (approximately 1 g) were immersed in liquid nitrogen and ground to a fine powder. Four milliliter of lysis buffer (5 mM Tris-HCl, pH 7.4, 1 mM PMSF, 2 mM EDTA, 10 mM DTT, and 1%TritonX-100) was added to the powder and subjected to ultrasonic vibrations for 15 min. The supernatant was transferred to a 50 mL tube after centrifugation at 25,000 g for 20 min; then, five volumes of cold acetone was added and incubated at −20°C for 2 h. Thereafter, the tubes were centrifuged at 16,000 g for 20 min and the supernatants discarded. The pellets were resuspended in the lysis buffer and centrifuged as described above. Finally, the protein pellets were washed twice with 30 ml of ice-cold acetone, lyophilized and stored at −80°C.

Protein digestion, iTRAQ labeling and strong cation exchange

iTRAQ analysis was performed at the Beijing Genomics Institute (BGI, Shenzhen, China). Protein samples (100 μg of each protein) were reduced with 10 mM DTT at 56°C for 2 h, alkylated with 55 mM iodoacetamide at room temperature in the dark for 45 min, digested with trypsin at 20:1 mass ratio at 37°C for 12 h, then labeled using the iTRAQ Reagents 8-plex kit according to the manufacturer's instructions (AB Sciex Inc., MA, USA). The salt-treated samples' replicates were labeled with iTRAQ tags 113, 114, and the untreated labeled with tags 115, 116, respectively. After labeling, the samples were mixed and lyophilized before dissolving in 4 mL of strong cation exchange (SCX) buffer A (25 mM NaH2PO4 in 25% acetonitrile pH2.7). The peptides were fractionated on Ultremex SCX column (4.6 × 250 mm) using Shimadzu LC-20AB HPLC. The subsequent experiment was performed as described in Zhu et al. (2009).

Tandem mass spectrometry analysis

The fractionated samples were lyophilized to remove acetonitrile and resuspended in Solvent A (5% acetonitril, 0.1% formic acid). Peptide samples, 5 μL (2.5 μg) each were loaded onto a C18LC Packings PepMap trap column and separated on a PepMapC18 75 μm inner diameter (LC Packings) column at a flow rate of 300 nl/min using Shimadzu LC-20AD HPLC. Peptides were eluted from the HPLC column by a linear gradient from 2% buffer B (95% acetonitrile, 1% formic acid) to 35% for 40 min, followed by ramping up to 80% buffer B for 5 min, and then held on 80% buffer B for 4 min. Peptides separated by liquid chromatography were sprayed into the orifice of the Q-Exactive mass spectrometer (Thermo Fisher Scientific, San Jose, CA, USA) and subsequently analyzed according to previously described methods (Qiao et al., 2012).

Data analysis and protein identification

The MS data were performed using Proteome Discoverer 1.3 software (Thermo Fisher Scientific, San Jose, CA, USA). Relative abundance quantitation and protein identification were processed using Mascot 2.3.02 (Matrix Science, London, United Kingdom). The analysis was carried out with cotton AD genome annotation database (81147 sequences) and the National Center for Biotechnology Information (NCBI) non-redundant fasta database (6833826 sequences). The search parameters were set as follows: Type of search: MS/MS Ion search; Enzyme: Trypsin with one missed cleavage; Monoisotopic mass; Tragment Mass Tolerance: 0.02Da; Peptide Mass Tolerance: 15 ppm; oxidation of methionine and tyrosine labeled by iTRAQ 8-plex as variable modifications, while carbamidomethylation on cysteine, iTRAQ 8-plex labeled N-term and lysine as fixed modifications. False discovery rate (FDR) correction was adopted with a threshold of 0.01 to reduce the false identification of peptide, and a Mascot probability of 95% was set for the identification and quantification of protein. Protein identification was considered if at least one unique peptide was identified for each protein.

Go and KEGG analysis

Differentially expressed proteins were classified according to Gene Ontology (http://www.geneontology.org). Kyoto Encyclopedia of Genes and Genomes (KEGG) (http://www.genome.jp/kegg/ or http://www.kegg.jp/) was used to predict molecular function, biological processes and significant pathways involved in response to salt stress.

Measurement of enzyme activities

The activities of SOD, POD and MDAR were assayed according to Chen et al. (2008). GST and MDH were extracted and assayed according to Gronwald et al. (1987) and Chen et al. (2009), respectively.

qRT-PCR analysis

Total RNA was extracted from salt-treated and control cotton roots by Trizol reagent (TaKaRa), and cDNA was reverse transcribed from 1 μg of to total RNA using a First Strand cDNA Synthesis Kit (Invitrogen). Gene-specific primers (GSPs) used for qRT-PCR were designed using primer3 (http://primer3.ut.ee/) according to cDNA sequences obtained from the cotton (Table S1). The cotton 18s-rRNA gene was used as an endogenous control for normalization. The PCR reaction was carried out in a 20 uL volume containing 10 μL 2 × SYBR Green Master Mix reagent (TaKaRa), 1 μL template cDNA and 0.5 μL of each GSPs with the following reaction conditions: 95°C for 30 s; followed by 40 cycles of 95°C for 10 s; 55°C for 10 s and 72°C for 15 s. Relative gene expression was calculated using the ddCt alogorithm (Zhang et al., 2003).

Results

Primary data analysis and protein detection

A total of 458,751 spectra were generated from the iTRAQ experiment using the proteins of salt-treated and untreated roots as materials. The data were analyzed using Mascot software (version 2.3.02). Mascot detected a total of 11,191 spectra matched to known spectra, 8022 spectra matched to unique spectra, 5603 peptides, 4339 unique peptides, and 1649 proteins (Figure 1). The distribution of the number of peptides defining each protein is shown in Figure 2 and over 64.7% of the proteins included at least two peptides. These proteins were involved in multiple metabolic, regulatory and defense pathways (Figure 3A, Table S2).

Figure 1

Figure 2

Figure 3

Identification and functional classification of DEPs

DEPs were selected based on the following criteria: (i) proteins in which the mean ratio {corresponding to the protein reporter ion intensity originating from salt-treated protein samples (113 and 114) with respect to fully control protein samples (115 and 116)} had a 1.5 fold change; (ii) a p < 0.05. Based on these criteria, 128 DEPs were identified in cotton roots, 76 (59.4%) of which displayed increased, and 52 (40.6%) decreased abundance under salt stress conditions. The main biological functions for the 128 DEPs were: carbohydrate and energy metabolism (13.3%), transcription related (4.7%), protein metabolism (15.6%), cell wall and cytoskeleton metabolism (12.5%), membrane, and transport (10.9%), signal transduction (3.1%), and stress and defense (23.4%). In addition, six proteins were involved in other metabolic processes (4.7%) and 15 in unknown biological processes (11.7%). Detailed information can be found in Figure 3B, Table 1, Figure S1 and Table S3.

Table 1

AccessionProteinsSpeciesPercent coverageNo. of unique peptideMean ratioaUp/downb
CARBOHYDRATE AND ENERGY METABOLISM
1gi|377824753Pectin methylesteraseGossypium hirsutum11.640.666↓
2gi|211906490Malate dehydrogenaseGossypium hirsutum28.920.664↓
3gi|122216326Perakine reductaseRauwolffia serpentina28.030.666↓
4gi|339265919Phosphogluconate dehydrogenaseLotus grandiflorus17.931.521↑
5gi|225455555EnolaseGlycine max38.540.679↓
6gi|55584187Quinone oxidoreductase-like proteinArabidopsis thaliana20.840.623↓
7gi|117940179Dihydrolipoyllysine-residue Acetyltransferase component 1 of pyruvate dehydrogenase complex, mitochondrialArabidopsis thaliana10.511.625↑
8gi|225465847NADPH: quinone oxidoreductaseVitis vinifera25.930.570↓
9gi|75262456ATP-citrate synthase beta chain protein 2Arabidopsis thaliana21.930.671↓
10gi|242129048ATP synthase delta subunit 2Gossypium hirsutum21.920.566↓
11gi|91981275Pectin methylesteraseCitrus bergamia2.110.571↓
12gi|21431800NADP-dependent alkenal double bond reductase P2Arabidopsis thaliana10.720.576↓
13gi|75268018Probable fructose-bisphosphate aldolase 3Arabidopsis thaliana5.621.485↑
14gi|470127114Aldose 1-epimerase-likeFragaria vesca subsp. vesca13.141.517↑
15gi|37193998Phosphoenolpyruvate carboxykinaseMitella japonica6.630.466↓
16gi|224057577Glutathione reductasePopulus trichocarpa7.531.992↑
17gi|356532527Dihydrolipoyl dehydrogenase-likeGlycine max8.320.503↓
TRANSCRIPTION RELATED
18gi|75262442Nuclear transcription factor Y subunit B-2Arabidopsis thaliana15.351.511↑
19gi|341958560CASP-like proteinPopulus trichocarpa7.510.605↓
20gi|356521678Putative DNA repair protein RAD23-1-like isoform 1Glycine max7.811.750↑
21gi|55976204Transcription factor HY5Solanum lycopersicum10.112.103↑
22gi|224133758Histone H1Populus trichocarpa3.811.846↑
23gi|75321585Zinc finger CCCH domain-containing protein 40Oryza sativa subsp. japonica2.011.681↑
PROTEIN TRANSLATION, PROCESSING, AND DEGRADATION
24gi|7526634240S ribosomal protein S20-2Arabidopsis thaliana9.812.103↑
25gi|2209637940S ribosomal protein S30Arabidopsis thaliana7.410.674↓
26gi|1786556660S ribosomal protein L36-3Arabidopsis thaliana21.821.743↑
27gi|2447379660s acidic ribosomal proteinPrunus dulcis57.020.645↓
28gi|6015064Elongation factor 1-deltaPimpinella brachycarpa26.022.089↑
29gi|23503072Eukaryotic translation initiation factor 3 subunit IArabidopsis thaliana15.330.666↓
30gi|18803Polyubiquitin proteinHelianthus annuus16.240.675↓
31gi|95116512Ubiquitin activating enzymeTheobroma cacao3.830.564↓
32gi|356560787Ubiquitin-conjugating enzyme E2 5-likeGlycine max10.321.543↑
33gi|117949833T-complex protein 1 subunit gammaArabidopsis thaliana8.441.578↑
34gi|75115360Protein disulfide isomerase-like 1-6Arabidopsis thaliana3.221.911↑
35gi|255541132Structural constituent of nuclear poreRicinus communis4.831.680↑
36gi|255575861Glycolipid transfer proteinRicinus communis4.111.972↑
37gi|82581521Proteasome subunit beta type-4Arabidopsis thaliana5.421.526↑
38gi|211906494Heat shock protein 70Gossypium hirsutum32.810.654↓
39gi|211906504Heat shock protein 70Gossypium hirsutum37.210.675↓
40gi|289064666S-adenosylmethionine synthase-like proteinEperua falcata34.610.536↓
41gi|255543963Aspartic proteinase nepenthesin-1 precursorRicinus communis19.851.819↑
42gi|229830633L-idonate 5-dehydrogenaseVitis vinifera3.810.601↓
43gi|308743337Asparagine synthetase 1Solanum tuberosum6.630.537↓
CELL WALL AND CYTOSKELETON METABOLISM
44gi|89212812Actin depolymerizing factor 2Gossypium hirsutum44.441.524↑
45gi|117553550Actin-binding protein ABP29Vitis vinifera15.911.663↑
46gi|182627650Actin-related protein4Oryza sativa subsp. indica6.031.503↑
47gi|157273642Fasciclin-like arabinogalactan protein 4Gossypium hirsutum28.310.603↑
48gi|150416583Fasciclin-like arabinogalactan protein 11Gossypium hirsutum28.310.652↑
49gi|157273666Fasciclin-like arabinogalactan protein 16Gossypium hirsutum13.031.563↓
50gi|157273646Fasciclin-like arabinogalactan protein 6Gossypium hirsutum15.931.711↑
51gi|157273640Fasciclin-like arabinogalactan protein 3Gossypium hirsutum12.021.574↑
52gi|157273638Fasciclin-like arabinogalactan protein 2Gossypium hirsutum14.033.059↑
53gi|157273660Fasciclin-like arabinogalactan protein 13Gossypium hirsutum5.622.124↑
54gi|253509569Caffeic acid O-methyltransferase 2Gossypium hirsutum23.050.682↓
55gi|224552010Hybrid proline-rich proteinGossypium hirsutum5.213.507↑
56gi|255547195Glycine-rich RNA-binding proteinRicinus communis27.311.502↑
MEMBRANE AND TRANSPORT
57gi|1336803Vacuolar H(+)-ATPase subunit AGossypium hirsutum48.521.635↑
58gi|2493146V-type proton ATPase 16 kDa proteolipid subunitGossypium hirsutum10.922.155↑
59gi|75273758Cysteine-rich repeat secretory protein 38Arabidopsis thaliana22.730.628↓
60gi|7105717Plasma membrane proton ATPaseKosteletzkya virginica17.541.732↑
61gi|224130846Multidrug/pheromone exporter, MDR family, ABC transporter familyPopulus trichocarpa6.631.699↑
62gi|292653531Aquaporin TIP1;7Gossypium hirsutum6.310.641↓
63gi|292653547Aquaporin TIP2;5Gossypium hirsutum4.410.658↓
64gi|461929Probable aquaporin TIP-typeAntirrhinum majus7.010.503↓
65gi|300793602TPA: TPA_inf: aquaporin TIP1;4Gossypium hirsutum6.810.451↓
66gi|292653535Aquaporin TIP1;10Gossypium hirsutum6.210.506↓
67gi|300793598Aquaporin PIP2;10Gossypium hirsutum21.520.291↓
68gi|256568429PIP proteinGossypium hirsutum20.720.541↓
69gi|164668308PIP2 proteinGossypium hirsutum10.220.522↓
70gi|118132686PIP1 proteinGossypium hirsutum3.510.270↓
SIGNAL TRANSDUCTION
71gi|363807628Probable leucine-rich repeat receptor-like protein kinaseGlycine max9.931.454↑
72gi|1702983Auxin-repressed 12.5 kDa proteinFragaria ananassa69.441.586↑
73gi|349504495Leucine rich repeat-containing proteinCorchorus capsularis1.521.572↑
74gi|1346675Nucleoside diphosphate kinase BFlaveria bidentis16.921.523↑
STRESS AND DEFENSE
75gi|74229677Cytoplasmic Cu/ZnSODGossypium hirsutum16.512.940↑
76gi|357470271PeroxidaseMedicago truncatula29.781.751↑
77gi|115345276PeroxidasePopulus alba6.511.667↑
78gi|73913500PeroxidasePhaseolus lunatus19.951.673↑
79gi|255551599Peroxidase 26 precursorRicinus communis12.221.981↑
80gi|32351452Class III peroxidaseGossypium hirsutum10.621.805↑
81gi|255581003Peroxidase 2 precursorRicinus communis19.841.643↑
82gi|25453205Peroxidase 12Arabidopsis thaliana12.531.700↑
83gi|225447324Peroxidase 27Vitis vinifera20.210.674↓
84gi|220967704Monodehydroascorbate reductaseSolanum lycopersicum14.230.667↓
85gi|195973264Glutathione S-transferaseGossypium hirsutum18.731.648↑
86gi|354620267pCPR10-16Gossypium barbadense49.711.684↑
87gi|15811629Ribonuclease-like PR-10Gossypium arboreum30.821.797↑
88gi|33338347Osmotin-like pathogenesis-related proteinGossypium hirsutum15.331.546↑
89gi|255537367Osmotin precursorRicinus communis8.511.780↑
90gi|383932370Nodulin-like proteinGossypium hirsutum27.420.612↓
91gi|38258655Monocopper oxidase-like protein SKU5Arabidopsis thaliana10.541.444↑
92gi|354620271MLPGossypium barbadense19.110.625↓
93gi|194321204LaccaseGossypium hirsutum1.311.680↑
94gi|65998365Dirigent-like proteinGossypium barbadense12.411.788↑
95gi|66276977Dirigent-like proteinGossypium barbadense12.611.848↑
96gi|118926Desiccation-related protein PCC13-62Craterostigma plantagineum20.542.237↑
97gi|359480830L-ascorbate oxidase-likeVitis vinifera11.641.565↑
98gi|166203457Universal stress protein 1Gossypium arboreum11.011.947↑
99gi|94717590GDP-mannose 3,5-epimerase 2Oryza sativa subsp. Japonica12.531.529↑
100gi|225455388Germin-like protein 11-1Vitis vinifera14.710.603↓
101gi|470122858Plant cadmium resistance 2-like isoform 1Fragaria vesca subsp. vesca8.331.919↑
102gi|75099392Subtilisin-like proteaseArabidopsis thaliana4.611.832↑
103gi|68064400Thaumatin-like proteinPhaseolus vulgaris8.921.714↑
104gi|319433441Copper binding protein 3Gossypium hirsutum13.010.551↓
105gi|259016223Glucan endo-1,3-beta-glucosidase 7Arabidopsis thaliana5.723.108↑
106gi|255546283Glucan endo-1,3-beta-glucosidase precursorRicinus communis7.532.234↑
107gi|255573702Glucan endo-1,3-beta-glucosidase precursorRicinus communis2.311.934↑
OTHER METABOLISM
108gi|89258498Short chain alcohol dehydrogenaseGossypium hirsutum39.110.664↓
109gi|7546402Chain A, Structures of adenylosuccinate Synthetase from Triticum aestivum and Arabidopsis thalianaArabidopsis Thaliana19.461.553↑
110gi|74273629Gibberellin 20-oxidase 1Gossypium hirsutum13.331.631↑
111gi|395406786Putative inactive methylesterase 20Arabidopsis thaliana6.420.656↓
112gi|255554698Homogentisate 1,2-dioxygenaseRicinus communis9.530.629↓
113gi|3183454Uncharacterized oxidoreductase ykwCBacillus subtilis2.912.067↑
UNKNOWN
114gi|224137260Predicted proteinPopulus trichocarpa23.822.309↑
115gi|297736988Unnamed protein productVitis vinifera19.561.889↑
116gi|224106732Predicted proteinPopulus trichocarpa3.421.775↑
117gi|297736988Unnamed protein productVitis vinifera25.681.717↑
118gi|225458697Uncharacterized proteinVitis vinifera12.751.660↑
119gi|147767808Hypothetical protein VITISV_032830Vitis vinifera8.021.596↑
120gi|224092318Predicted proteinPopulus trichocarpa12.651.521↑
121gi|147820236Hypothetical protein VITISV_010210Vitis vinifera12.681.494↑
122gi|217073300UnknownMedicago truncatula14.730.667↓
123gi|359496362Uncharacterized protein LOC100854560Vitis vinifera5.610.664↓
124gi|388499178UnknownLotus japonicus13.220.645↓
125gi|223943077UnknownZea mays17.850.644↓
126gi|388500070Unknown knownLotus japonicus29.420.582↓
127gi|84453208Putative cytosolic factorTrifolium pratense13.530.547↓
128gi|147782603Hypothetical protein VITISV_010455Vitis vinifera19.430.447↓

Differentially expressed proteins in cotton roots subject to salt stress (200 mM NaCl).

a

Mean ratio corresponds to the protein reporter ion intensity originating from salt-treated protein samples (113 and 114) relative to fully control protein samples (115 and 116) with a 1.5 fold-changes and a p < 0.05.

b

Proteins increased in abundance (↑) or decreased in abundance (↓).

Analysis of differentially expressed enzymes

Under the same conditions, the level of activity is positively correlated with the enzyme protein abundance (Yang et al., 2013). To validate the DEPs, five enzymes involved in ROS scavenging and organic acid metabolism were selected for activity analysis. The activities of SOD, POD, and GST were higher in NaCl-treated roots than in the control, whereas the activities of MDH and MDAR were lower (Figure 4). These results agree with the protein profiles of the iTRAQ analysis.

Figure 4

Transcriptional analysis of genes for some DEPs

In order to assess the correlation of expression levels between mRNA and protein, qRT-PCR was applied to five DEP genes (POD, SOD, GST, MDAR and MDH) as shown in Figure 5. The expression of the former four genes (POD, SOD, GST and MDAR) is consistent with the corresponding DEPs, indicating that the expression of these proteins is regulated at the transcriptional level, but this was not the case for MDH (Table 1).

Figure 5

Discussion

Carbohydrate and energy metabolism

Salt stress alters the abundance of many proteins involved in carbon and energy metabolism, including glycolysis, the tricarboxylic acid cycle (TCA), and the pentose phosphate pathway (PPP) in cotton roots. It was found that an FBP3 protein (gi|75268018) was increased. FBP aldolase, a key enzyme in the glycolytic pathways, plays an important role in the production of water-soluble carbohydrates, triose phosphates metabolism and signal transduction (Schaeffer et al., 1997). It was reported that overexpression of FBP could enhance salt tolerance in tobacco by increasing proline content (Konishi et al., 2005). In the present study, up-regulation of FBP3 aldolase increased levels of sugars and starch, and may improve the growth of cotton roots under stress. In our experiment, an enolase protein, which catalyzes the formation of a high-energy phosphoenol pyruvate from 2-phosphoglycerate in the glycolytic pathway, displayed a decrease in abundance after 24 h of salt stress treatment. This result disagreed with the expression profile of this same protein in wheat (Guo et al., 2012a) and cucumber (Du et al., 2010) under salt stress.

Three proteins related to TCA were identified. Pyruvate dehydrogenase E2 (gi|117940179) is an enzyme component of the multienzyme pyruvate dehydrogenase complex and is involved in the formation of cellular energy during the TCA cycle. In this study, the abundance of this enzyme increased under salt stress. However, MDH and ATP-citrate synthase beta chain protein 2 decreased. This suggests that the TCA cycle was inhibited in cotton roots after 24 h of salt stress treatment.

Some ROS-scavenging systems need the PPP pathway that produces NADPH under stress conditions. Phosphogluconate dehydrogenase (PGD, gi|339265919)—the key regulatory enzyme of the PPP pathway—was enhanced under salt stress conditions. Manaa et al. (2011) reported that the PGD activity increased under salt conditions in tomato roots.

The abundance of ATP synthase delta subunit 2 (gi|242129048) decreased under stress. This result supports the previous data on the expression profile of this protein in the roots of Arabidopsis (Jiang et al., 2007), rice (Chitteti and Peng, 2007), and cucumber (Du et al., 2010) under salt stress.

Thus, the flexibility of carbohydrate and energy metabolism may help cotton survive under salt stress conditions.

Proteins involved in transcription

Transcriptional regulation of salt-responsive genes is a crucial part of the plant response to various abiotic and biotic stresses (Jiang et al., 2007). Previous studies showed that chromatin-mediated regulation of gene expression plays an important role in the response to abiotic stress and that histone H1 is involved in stress-induced reactions (Kim et al., 2010). In our data, the expression of histone H1 in salt-treated samples was nearly twice as high as that in the control samples, indicating its role in the salt stress response. Moreover, a zinc-finger transcription factor (gi|75321585) showed higher abundance in roots under salt stress conditions. Zinc finger proteins are well characterized in the regulation of stress responses (Chinnusamy et al., 2006), and the overexpression of CCCH-type zinc finger proteins AtSZF1 and AtSZF2 enhanced salt tolerance in Arabidopsis (Sun et al., 2007).

Protein metabolism

Protein turnover, the balance between synthesis and degradation, is one of the many forms of regulation that is employed to achieve a unified cellular response (Reinbothe et al., 2010). Several proteins, involved in protein translation, processing and degradation, were identified in these iTRAQ data. The abundance of two ribosomal proteins (gi|75266342 and gi|17865566) decreased, whereas two (gi|7526634 and gi|17865566) increased in the present study. Rodriguez-Uribe et al. (2011) has also reported that levels of some of the ribosomal proteins decreased while some specific ribosomal components increased under salt stress. Moreover, our data showed lower expression of a eukaryotic translation initiation factor 3 subunit I (eIF3I, gi|23503072) under salinity, which is consistent with a previous report on Arabidopsis (Jiang et al., 2007). In addition, elongation factor gi|6015064 displayed higher abundance under salt stress conditions. The differential regulation of different components of the translation machinery indicates that complicated regulation mechanisms may govern protein synthesis in order to help plants cope with salt stress.

Proper protein folding and processing is important for normal cellular function under salt stress. Here, it was found that the expression of T-complex protein 1 (TCP1, gi|117949833) and protein disulfide isomerase-like 1-6 (PDIL1-6, gi|75115360) were enhanced. TCP1 is involved in proper folding and assembly of proteins to cope with salinity in wheat roots (Wang et al., 2008). PDIs are molecular chaperones that aid the formation of proper disulfide bonds during protein folding (Houston et al., 2005). Two isoforms of PDIs increased in rice roots under salt stress (Nohzadeh Malakshah et al., 2007). Hsp70 s assists in proper folding of newly synthesized polypeptides and import/translocation of precursor proteins. Two hsp70 members (gi|211906494, gi|211906504) showed lower abundance in NaCl-treated roots. This result is consistent with the expression profile of this protein in Arabidopsis (Jiang et al., 2007). The ubiquitin/26S proteasome pathway selectively degrade key regulatory proteins and enzymes under salt stress conditions (Vierstra, 2003). The abundance of some components of ubiquitin/26S proteasome system including ubiquitin-activating enzyme (gi|95116512) and polyubiquitin protein (gi|18803) decreased under salt stress conditions. These findings suggest that decreased protein degradation compensates for decreased protein biosynthesis in roots under salt stress.

Cell wall and cytoskeleton metabolism

The cytoskeleton is rapidly remodeled to allow cell size adjustment in order to maintain normal cell turgor pressure under salt stress conditions (Zhang et al., 2012). In salt-treated roots, it was found that three actin-binding proteins (ABPs), including actin depolymerizing factor (ADF), actin-related protein 2 (ARP2), and actin-binding protein 29 (ABP29), can bind to actin cytoskeletons and effect remodeling. In a previous study, ABP29 from Lilium pollen played an important role in the remodeling of the actin cytoskeleton during pollen germination and pollen tube growth (Xiang et al., 2007). Thus, depolymerization and subsequent reorganization of the actin cytoskeleton enhanced salt tolerance in cotton roots.

Some proteins, including glycine-rich proteins (GRPs), proline-rich protein (PRPs), and arabinogalactan proteins (AGPs), are essential structural protein components of the cell walls of many higher plants. We found that a GRP (gi|255547195) and a hybrid PRP (HyPRP gi|224552010) displayed higher abundance in roots under salt stress. Biosynthesis of GRPs and their accumulation in vascular tissues are part of the plant's defense mechanism (Mousavi and Hotta, 2005). Overexpression of HyPRP (encoding a HyPRP) in Arabidopsis enhanced germination under cold and high salinity stress conditions (Qin et al., 2013). A sub-group of AGPs that include one or two AGP domains and one or two copies of the fasciclin domain are termed fasciclin-like arabinogalactan protein (FLAs). FLAs, which are located in the cell wall/plasma membrane and cell surface, have many developmental roles. Some FLAs are involved in microspore and lateral root/shoot development, maintaining proper cell expansion and/or keeping the integrity and elasticity of cell wall matrix in Arabidopsis (Johnson et al., 2003). In this present study, five FLAs (gi|157273666, gi|157273646, gi|157273640, gi|157273638, gi|157273660) displayed increased in abundance, but two (gi|157273642, gi|150416583) decreased. The diverse expression of FLAs suggests that these proteins may be involved in a wide range of biological process under salt stress conditions.

Membrane and transport

Under salinity conditions, Na+/K+ ratios and Na+ concentration increase in plant roots causing hyperosmotic stress, ion imbalance and toxicity (Zhao et al., 2013). H+-ATPase plays an essential role in the maintenance of ion homeostasis in plant cells. The plasma membrane H+-ATPase in tomato (Kerkeb et al., 2001) and the vacuolar H+-ATPase in the roots of Arabidopsis (Jiang et al., 2007), rice (Cheng et al., 2009), wheat (Guo et al., 2012a), tomato (Manaa et al., 2011), and cucumber (Du et al., 2010) are induced under salt stress conditions. Here, increased abundance of one plasma membrane H+-ATPase (gi|7105717) and two vacuolar H+-ATPases (gi|1336803, gi|2493146) indicates that the increased activities of these enzymes are considered to be a cost-effective strategy for osmotic adjustment, which reduces the Na+ concentration in the cytosol in plants under salt stress conditions.

ABC transporters transport stress-related secondary metabolites such as alkaloids, terpenoids, polyphenols and quinines (Theodoulou, 2000). In Arabidopsis, ABC transporter affected Na+/K+ homeostasis and elicited a salt stress response (Lee et al., 2004). The up-regulation of an ABC transporter (gi|224130846) in cotton roots suggests that it may play an important role in salt-stressed responses.

Aquaporins (AQPs)—channel proteins that facilitate the transport of water and/or small neutral solutes or gasses in the plasma and intracellular cell membranes—are associated with plant stress tolerance (Wang et al., 2011). PIPs and TIPs, two subfamilies of AQP, are most abundant in the plasma membrane and vacuolar membrane, respectively (Danielson and Johanson, 2008). HvPIP2:1 was down-regulated in barley seedlings, and its overexpression enhanced salt sensitivity in transgenic rice under salt stress conditions (Katsuhara et al., 2003). Overexpression of the Panax ginseng TIP2:1 gene in Arabidopsis enhances tolerance to salt stress, but overexpression of GsTIP2:1 depresses salt tolerance and dehydration stress (Wang et al., 2011). Thus, the regulation mechanism of AQPs under salt stress conditions is complicated and requires further study (Peng et al., 2007). Here, four PIPs (gi|300793598, gi|256568429, gi|118132686, gi|118132686) and five TIPs (gi|29265353, gi|292653547, gi|461929, gi|300793602, gi|292653535) showed lower abundance in response to salt stress. This may be attributed to the reduced hydraulic conductivity of membranes to prevent water loss under salt stress conditions (Sutka et al., 2011).

Signal transduction

LRR-RLKs function in a wide variety of signal transduction pathways related to hormone and abiotic stress responses (Hove et al., 2011). Potato LRPK1 functions under diverse stress conditions, such as wounding, and high-, low-temperature, and salinity stress (Wu et al., 2009). In our present data, the up-regulation of LRR-RLKs (gi|363807628) imply it has a role in Na+ and plant interactions, specific recognition, and signal transduction leading to an induced salt-stressed tolerance. Nucleoside diphosphate kinase B (NDPKB, gi|1346675) is an enzyme that converts GTP to ATP, and is involved in the H2O2 mediated mitogen-activated protein kinase signaling pathway. NDPK increased tolerance in response to NaCl in Arabidopsis, creeping bentgrass and rice (Jiang et al., 2007; Seong et al., 2007; Xu et al., 2010).

Stress and defense

Salt stress causes the production of excessive reactive oxygen species (ROS), which oxidize cellular components and irreversibly damage plant cells (Askim et al., 2014). ROS can be scavenged in plants by SOD, POD and GSTs. Ten of these proteins were identified in this study (Table 1). In most cases, higher expression of these proteins was found in salt-treated samples than in the control. Increased accumulation of SOD was noted in the roots of Arabidopsis (Jiang et al., 2007), wheat (Guo et al., 2012a), cucumber (Du et al., 2010), and salt cress (Zhou et al., 2010) under salt stress conditions. The up-regulation in abundance of Cu/ZnSOD (gi|74229677) also indicates that it helps cope with salt stress in cotton. PODs catalyze the reduction of H2O2 using electron donors such as lignin precursors, phenolic compounds, auxins and secondary metabolites (Zhao et al., 2013). In the present study, the levels of seven POD isozymes (gi|357470271, gi|115345276, gi|73913500, gi|255551599, gi|32351452, gi|255581003, gi|25453205) increased in response to salt stress but POD 27 (gi|225447324) did not. The levels of POD increased in the salt-stressed roots of wheat (Peng et al., 2009), barley (Witzel et al., 2009), cucumber (Du et al., 2010), and rice (Cheng et al., 2009) but decreased in creeping bentgrass (Xu et al., 2010). GST increased in the salt-stressed roots of Arabidopsis (Jiang et al., 2007), rice (Chitteti and Peng, 2007), barley (Witzel et al., 2009), and wheat (Peng et al., 2009). Here, higher levels of a GST (gi|195973264) were also observed in salt-stressed cotton roots. GSTs may play a pivotal role in preventing the degradation of organic hydroperoxides to cytotoxic aldehyde derivatives under salt stress conditions in cotton. Thus, it is demonstrated that antioxidant enzymes protect salt-stressed cotton roots from oxidative damage.

MDAR catalyzes the reduction of monodehydroascorbate to ascorbate (ASA) and is essential in order to maintain a reduced pool of ascorbate. Germin-like proteins (GLP) possess both oxalate activity and SOD activity. Here, the decreased expression of MDAR (gi|220967704) and GLP (gi|225455388) was identified in salt- stressed cotton roots. It is suggested that although plants require MDAR and GLP in order to eliminate ROS, the fine tuning of the levels of various antioxidants is also an important consideration in stress responses (Lisenbee et al., 2005).

In addition to the redox related proteins, plants have developed cross-tolerance mechanisms to be able to cope with different stresses (Zhang et al., 2012). Some biotic and abiotic stress-responsive proteins play important roles in salt tolerance (Table 1). Some biotic stress-related proteins were induced under salt stress conditions, such as pCPR10-16 (gi|354620267), ribonuclease-like PR-10 (gi|15811629), osmotin-like pathogenesis-related proteins (gi|3333834), thaumatin-like protein (TLP, gi|68064400), USP (gi|166203457), and glucan endo-1,3-beta-glucosidases (gi|255546283, gi|259016223, gi|255573702). PR10 mediates tolerance to heavy metals (Wang et al., 2014) and pathogen attack (Coumans et al., 2009). TLP, a subgroup of pathogenesis-related proteins, is induced by phytohormones (SA, JA, and ABA) and stress stimuli (wounding, cold temperature and high salinity) (Wang et al., 2010). Overexpression of the GbTLP1 in tobacco enhances resistance to Verticillium dahliae, salinity and drought (Munis et al., 2009). USP helps cotton plants adapt to water stress (Maqbool et al., 2009). Glucan endo-1, 3-beta-glucosidase accumulates in rice in response to ABA and salt stress (Li et al., 2010). Moreover, some abiotic stress-related proteins, e.g. DIR (gi|65998365, gi|66276977), desiccation-related protein PCC13-62 (gi|118926) also respond to salt stress (Bartels et al., 1990; Guo et al., 2012b). DIR is involved in the response to drought, salts and oxidation (Guo et al., 2012b). Desiccation-related protein PCC13-62 promotes the plant's tolerance to extreme desiccation (Bartels et al., 1990). These proteins provide novel insights into the understanding of the cross-tolerance mechanisms in roots in response to biotic and abiotic stress.

The correlation of protein abundance and gene expression

There might be a weak correlation between the transcript levels of genes and their protein abundance (Yang et al., 2013). The discrepancy between protein and mRNA expression may be caused by the various levels of regulation, e.g., post-transcriptional, translational or post-translational regulation (Tian et al., 2004). A discrepancy between transcript levels of MDH and the abundance of the corresponding proteins (Figure 5E, Table 1) highlights the effect of post-transcriptional modifications.

Conclusion

An iTRAQ-based proteomic technique was employed to compare the abundance of proteins in untreated and salt-treated roots for 24 h. One hundred and twenty-eight DEPs were identified, 76 of which displayed increased abundance and 52 decreased under salt stress conditions. These DEPs are mainly involved in the biological processes of carbohydrate and energy metabolism, transcription, protein metabolism, cell wall and cytoskeleton metabolism, membrane and transport, signal transduction and stress and defense. The diverse array of proteins affected by salt stress conditions indicates that there is a remarkable flexibility in cotton root metabolism, which may contribute to its survival in salinity conditions. High positive correlation between the abundance of some altered proteins (SOD, POD, GST, MDAR, and MDH) and their enzyme activity demonstrates that the iTRAQ-based proteomic technique is sufficiently reliable for the identification and quantification of a large number of cotton root proteins. qRT-PCR results suggest that the expression of some proteins (e.g., MDH) can be regulated by post-transcriptional modifications. With this technology, many new salt-responsive proteins, such as ARP2, FLAs, TIPs, PIPs, LRR-RLKs, TLP, USP, DIR and the desiccation-related protein PCC13-62 were identified from cotton roots. These novel proteins provide a good starting point for further research into their functions using genetic or other approaches. These findings significantly improve the understanding of the molecular mechanisms involved in the tolerance of plants to salt stress.

Conflict of interest statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Statements

Acknowledgments

This work was supported by the National Special Key Program of Transgenic Species Breeding (No. 2015ZX08005-001-005) and the earmarked fund of Henan Academy of Agricultural Sciences (N0. 201412016).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2015.00732

    Abbreviations

  • ABPs

    actin-binding proteins

  • ADF

    actin depolymerizing factor

  • ARP2

    actin-related protein 2

  • AGPs

    arabinogalactan proteins

  • AQPs

    Aquaporins

  • DEPs

    differentially expressed proteins

  • DIR

    dirigent-like protein

  • FBP3

    fructose-bisphosphate aldolase 3

  • FLAs

    fasciclin-like arabinogalactan proteins

  • GST

    glutathione S-transferase

  • GRPs

    glycine-rich proteins

  • iTRAQ

    isobaric tag for relative and absolute quantitation

  • LRR-RLKs

    leucine-rich repeat receptor-like kinase-encoding genes

  • MDAR

    monodehydroascorbate reductase

  • MDH

    malate dehydrogenase

  • NDPK

    nucleoside diphosphate kinase

  • PPP

    pentose phosphate pathway

  • PGD

    phosphogluconate dehydrogenase

  • POD

    peroxidase

  • PDIL1-6

    protein disulfide isomerase-like 1-6

  • PRPs

    proline-rich proteins

  • PIPs

    plasma membrane intrinsic proteins

  • SOD

    superoxide dismutase

  • TCA

    tricarboxylic acid cycle

  • TCP1

    T-complex protein 1

  • TIPs

    tonoplast intrinsic proteins

  • TLP

    thaumatin-like protein

  • USP

    universal stress protein.

References

  • 1

    AghaeiK.EhsanpourA. A.ShahA. H.KomatsuS. (2009). Proteome analysis of soybean hypocotyl and root under salt stress. Amino Acids36, 91–98. 10.1007/s00726-008-0036-7

  • 2

    AhmadS.KhanN.IqbalM.HussainA.HassanM. (2002). Salt tolerance of cotton (Gossypium hirsutum L.). Asian J. Plant Sci.2, 715–719. 10.3923/ajps.2002.715.719

  • 3

    AskimH. S.RenginO.BarisU.IsmailT. (2014). Reactive oxygen species scavenging capacities of cotton (Gossypium hirsutum) cultivars under combined drought and heat induced oxidative stress. Environ. Exp.Bot. 99, 141–149. 10.1016/j.envexpbot.2013.11.010

  • 4

    BartelsD.SchneiderK.TerstappenG.PiatkowskiD.SalaminiF. (1990). Molecular cloning of abscisic acid-modulated genes which are induced during desiccation of the resurrection plant Craterostigma plantagineum. Planta181, 27–34. 10.1007/BF00202321

  • 5

    ChampionA.HebrardE.ParraB.BournaudC.MarmeyP.TranchantC.et al. (2009). Molecular diversity and gene expression of cotton ERF transcription factors reveal that group IXa members are responsive to jasmonate, ethylene and Xanthomonas. Mol. Plant Pathol. 10, 471–485. 10.1111/j.1364-3703.2009.00549.x

  • 6

    ChenF. X.LiuX. H.ChenL. S. (2009). Developmental changes in pulp organic acid concentration and activities of acid-metabolising enzymes during the fruit development of two loquat (Eriobotrya japonicaLindl.) cultivars differing in fruit acidit. Food Chem.114, 657–664. 10.1016/j.foodchem.2008.10.003

  • 7

    ChenL. S.LiP.ChengL. (2008). Effects of high temperature coupled with high light on the balance between photooxidation and photoprotection in the sun-exposed peel of apple. Planta228, 745–756. 10.1007/s00425-008-0776-3

  • 8

    ChengY.QiY.ZhuQ.ChenX.WangN.ZhaoX.et al. (2009). New changes in the plasma-membrane-associated proteome of rice roots under salt stress. Proteomics9, 3100–3114. 10.1002/pmic.200800340

  • 9

    ChinnusamyV.JagendorfA.ZhuJ. K. (2005). Understanding and improving salt tolerance in plants. Crop Sci.45, 437–448. 10.2135/cropsci2005.0437

  • 10

    ChinnusamyV.ZhuJ. H.ZhuJ. K. (2006). Gene regulation during cold acclimation in plants. Physiol. Plant.126, 52–61. 10.1111/j.1399-3054.2006.00596.x

  • 11

    ChittetiB. R.PengZ. (2007). Proteome and phosphoproteome differential expression under salinity stress in rice (Oryza sativa) roots. J. Proteome Res.6, 1718–1727. 10.1021/pr060678z

  • 12

    CoumansJ. V.PoljakA.RafteryM. J.BackhouseD.Pereg-GerkL. (2009). Analysis of cotton (Gossypium hirsutum) root proteomes during a compatible interaction with the black root rot fungus Thielaviopsis basicola. Proteomics9, 335–349. 10.1002/pmic.200800251

  • 13

    DanielsonJ. A.JohansonU. (2008). Unexpected complexity of the aquaporin gene family in the moss Physcomitrella patens. BMC Plant Biol.8:45. 10.1186/1471-2229-8-45

  • 14

    DuC. X.FanH. F.GuoS. R.TezukaT.LiJ. (2010). Proteomic analysis of cucumber seedling roots subjected to salt stress. Phytochemistry71, 1450–1459. 10.1016/j.phytochem.2010.05.020

  • 15

    FukaoY.FerjaniA.TomiokaR.NagasakiN.KurataR.NishimoriY.et al. (2011). iTRAQ analysis reveals mechanisms of growth defects due to excess zinc in Arabidopsis. Plant Physiol.155, 1893–1907. 10.1104/pp.110.169730

  • 16

    GaoS. Q.ChenM.XiaL. Q.XiuH. J.XuZ. S.LiL. C.et al. (2009). Cotton (Gossypium hirsutum) DRE-binding transcription factor gene, GhDREB, confers enhanced tolerance to drought, high salt, and freezing stresses in transgenic wheat. Plant Cell Rep.28, 301–311. 10.1007/s00299-008-0623-9

  • 17

    GhoshD.LiZ.TanX. F.LimT. K.MaoY.LinQ. (2013). ITRAQ based quantitative proteomics approach validated the role of calcyclin binding protein (CacyBP) in promoting colorectal cancer metastasis. Mol. Cell. Proteomics12, 1865–1880. 10.1074/mcp.M112.023085

  • 18

    GongB.ZhangC.LiX.WenD.WangS.ShiQ.et al. (2014). Identification of NaCl and NaHCO3 stress-responsive proteins in tomato roots using iTRAQ-based analysis. Biochem. Biophys. Res. Commun.446, 417–422. 10.1016/j.bbrc.2014.03.005

  • 19

    GronwaldJ. W.FuerstE. P.EberleinC. V.EgliM. A. (1987). Effect of herbicide antidotes on glutathione content and glutathione S-transferase activity of sorghum shoots. Pestic. Biochem. Phys.29, 66–76. 10.1016/0048-3575(87)90085-X

  • 20

    GuoG.GeP.MaC.LiX.LvD.WangS.et al. (2012a). Comparative proteomic analysis of salt response proteins in seedling roots of two wheat varieties. J. Proteomics75, 1867–1885. 10.1016/j.jprot.2011.12.032

  • 21

    GuoJ. L.XuL. P.FangJ. P.SuY. C.FuH. Y.QueY. X.et al. (2012b). A novel dirigent protein gene with highly stem-specific expression from sugarcane, response to drought, salt and oxidative stresses. Plant Cell Rep.31, 1801–1812. 10.1007/s00299-012-1293-1

  • 22

    GuoY. H.YuY. P.WangD.WuC. A.YangG. D.HuangJ. G.et al. (2009). GhZFP1, a novel CCCH-type zinc finger protein from cotton, enhances salt stress tolerance and fungal disease resistance in transgenic tobacco by interacting with GZIRD21A and GZIPR5. New Phytol. 183, 62–75. 10.1111/j.1469-8137.2009.02838.x

  • 23

    HoustonN. L.FanC.XiangJ. Q.SchulzeJ. M.JungR.BostonR. S. (2005). Phylogenetic analyses identify 10 classes of the protein disulfide isomerase family in plants, including single-domain protein disulfide isomerase-related proteins. Plant Physiol.137, 762–778. 10.1104/pp.104.056507

  • 24

    HoveC. A.BochdanovitsZ.JansweijerV. M.KoningF. G.BerkeL.Sanchez-PerezG. F.et al. (2011). Probing the roles of LRR RLK genes in Arabidopsis thaliana roots using a custom T-DNA insertion set. Plant Mol. Biol.76, 69–83. 10.1007/s11103-011-9769-x

  • 25

    JiangY.YangB.HarrisN. S.DeyholosM. K. (2007). Comparative proteomic analysis of NaCl stress-responsive proteins in Arabidopsis roots. J. Exp. Bot.58, 3591–3607. 10.1093/jxb/erm207

  • 26

    JinL. G.LiH.LiuJ. Y. (2010). Molecular characterization of three ethylene responsive element binding factor genes from cotton. J. Integr. Plant Biol.52, 485–495. 10.1111/j.1744-7909.2010.00914.x

  • 27

    JohnsonK. L.JonesB. J.BacicA.SchultzC. J. (2003). The fasciclin-like arabinogalactan proteins of Arabidopsis. A multigene family of putative cell adhesion molecules. Plant Physiol.133, 1911–1925. 10.1104/pp.103.031237

  • 28

    KatsuharaM.KoshioK.ShibasakaM.HayashiY.HayakawaT.KasamoK. (2003). Over-expression of a barley aquaporin increased the shoot/root ratio and raised salt sensitivity in transgenic rice plants. Plant Cell Physiol.44, 1378–1383. 10.1093/pcp/pcg167

  • 29

    KerkebL.DonaireJ. P.Rodríguez-RosalesM. P. (2001). Plasma membrane H+-ATPase activity is involved in adaptation of tomato calli to NaCl. Physiol Plant.111, 483–490. 10.1034/j.1399-3054.2001.1110408.x

  • 30

    KimJ. M.ToT. K.NishiokaT.SekiM. (2010). Chromatin regulation functions in plant abiotic stress responses. Plant Cell Environ.33, 604–611. 10.1111/j.1365-3040.2009.02076.x

  • 31

    KonishiH.KitanoH.KomatsuS. (2005). Identification of rice root proteins regulated by gibberellins using proteome analysis. Plant Cell Environ.28, 328–339. 10.1111/j.1365-3040.2005.01269.x

  • 32

    LeeE. K.KwonM.KoJ. H.YiH.HwangM. G.ChangS.et al. (2004). Binding of sulfonylurea by AtMRP5, an Arabidopsis multidrug resistance-related protein that functions in salt tolerance. Plant Physiol.134, 528–538. 10.1104/pp.103.027045

  • 33

    LiX. J.YangM. F.ChenH.QuL. Q.ChenF.ShenS. H. (2010). Abscisic acid pretreatment enhances salt tolerance of rice seedlings: proteomic evidence. Biochim. Biophys. Acta.1804, 929–940. 10.1016/j.bbapap.2010.01.004

  • 34

    LisenbeeC. S.LingardM. J.TreleaseR. N. (2005). Arabidopsis peroxisomes possess functionally redundant membrane and matrix isoforms of monodehydroascorbate reductase. Planta43, 900–914. 10.1111/j.1365-313X.2005.02503.x

  • 35

    LuW.ChuX.LiY.WangC.GuoX. (2013). Cotton GhMKK1 induces the tolerance of salt and drought stress, and mediates defence responses to pathogen infection in transgenic Nicotiana benthamiana. PLoS ONE8:e68503. 10.1371/journal.pone.0068503

  • 36

    MaasE. V.HoffmanG. J. (1977). Crop salt tolerance-current assessment. J. Irrig. Drain. Div. 103, 115–134.

  • 37

    ManaaA.Ben AhmedH.ValotB.BouchetJ. P.Aschi-SmitiS.CausseM.et al. (2011). Salt and genotype impact on plant physiology and root proteome variations in tomato. J. Exp. Bot.62, 2797–2813. 10.1093/jxb/erq460

  • 38

    MaqboolA.ZahurM.HusnainT.RiazuddinS. (2009). GUSP1 and GUSP2, two drought-responsive genes in Gossypium arboretum have homology to universal stress proteins. Plant Mol. Biol. Rep.27, 109–114. 10.1007/s11105-008-0049-0

  • 39

    MengC. M.CaiC. P.ZhangT. Z.GuoW. Z. (2009). Characterization of six novel NAC genes and their responses to abiotic stresses in Gossypium hirsutum L. Plant Sci.176, 352–359. 10.1016/j.plantsci.2008.12.003

  • 40

    MousaviA.HottaY. (2005). Glycine-rich proteins. Appl. Biochem. Biotech.120, 169–174. 10.1385/ABAB:120:3:169

  • 41

    MunisF. H.TuL. L.DengF. L.TanJ. F.XuL.XuS. C.et al. (2009). A thaumatin-like protein gene involved in cotton fiber secondary cell wall development enhances resistance against Verticillium dahliaeand other stresses in transgenic tobacco. Biochem. Biophys. Res. Commun.39, 38–44. 10.1016/j.bbrc.2010.01.069

  • 42

    MunnsR. (2002). Comparative physiology of salt and water stress. Plant Cell Environ.25, 239–250. 10.1046/j.0016-8025.2001.00808.x

  • 43

    Nohzadeh MalakshahS.Habibi RezaeiM.HeidariM.SalekdehG. H. (2007). Proteomics reveals new salt responsive proteins associated with rice plasma membrane. Biosci. Biotechnol. Biochem.71, 2144–2154. 10.1271/bbb.70027

  • 44

    PengY.LinW.CaiW.AroraR. (2007). Overexpression of a Panax ginseng tonoplast aquaporin alters salt tolerance, drought tolerance and cold acclimation ability in transgenic Arabidopsis plants. Planta226, 729–740. 10.1007/s00425-007-0520-4

  • 45

    PengZ. Y.WangM. C.LiF.LvH. J.LiC. L.XiaG. M. (2009). A proteomic study of the response to salinity and drought stress in an introgression strain of bread wheat. Mol. Cell Proteomics8, 2676–2686. 10.1074/mcp.M900052-MCP200

  • 46

    QiaoJ.WangJ.ChenL.TianX.HuangS.RenX.et al. (2012). Quantitative iTRAQ LC–MS/MS proteomics reveals metabolic responses to biofuel ethanol in cyanobacterial Synechocystis sp. PCC 6803. J. Proteome Res.11, 5286–5300. 10.1021/pr300504w

  • 47

    QinL. X.ZhangD. J.HuangG. Q.LiL.LiJ.GongS. Y.et al. (2013). Cotton GhHyPRP3 encoding a hybrid proline-rich protein is stress inducible and its overexpression in Arabidopsis enhances germination under cold temperature and high salinity stress conditions. Acta. Physiol. Plant.35, 1531–1542. 10.1007/s11738-012-1194-5

  • 48

    ReinbotheC.PollmannS.ReinbotheS. (2010). Singlet oxygen signaling links photosynthesis to translation and plant growth. Trends Plant Sci.15, 499–506. 10.1016/j.tplants.2010.05.011

  • 49

    Rodriguez-UribeL.HigbieS. M.StewartJ. M.WilkinsT.LindemannW.Sengupta-GopalanC.et al. (2011). Identification of salt responsive genes using comparative microarray analysis in Upland cotton (Gossypium hirsutum L.). Plant Sci.180, 461–469. 10.1016/j.plantsci.2010.10.009

  • 50

    SchaefferG. W.SharpeF. T.SicherR. C. (1997). Fructose 1,6-bisphosphate aldolase activity in leaves of a rice mutant selected for enhanced lysine. Phytochemistry46, 1335–1338. 10.1016/S0031-9422(97)00470-6

  • 51

    SeongE. S.GuoJ.KimY. H.ChoJ. H.LimC. K.Hyun HurJ.et al. (2007). Regulations of marker genes involved in biotic and abiotic stress by overexpression of the AtNDPK2 gene in rice. Biochem. Biophys. Res. Commun.363, 126–132. 10.1016/j.bbrc.2007.08.147

  • 52

    SunJ.JiangH.XuY.LiH.WuX.XieQ.et al. (2007). The CCCH-type zinc finger proteins AtSZF1 and AtSZF2 regulate salt stress responses in Arabidopsis. Plant Cell Physiol.8, 1148–1158. 10.1093/pcp/pcm088

  • 53

    SutkaM.LiG.BoudetJ.BoursiacY.DoumasP.MaurelC. (2011). Natural variation of root hydraulics in Arabidopsis grown in normal and salt-stressed conditions. Plant Physiol.155, 1264–1276. 10.1104/pp.110.163113

  • 54

    TheodoulouF. L. (2000). Plant ABC transporters. Biochim. Biophys. Acta.1465, 79–103. 10.1016/S0005-2736(00)00132-2

  • 55

    TianQ.StepaniantsS. B.MaoM.WengL.FeethamM. C.DoyleM. J.et al. (2004). Integrated genomic and proteomic analyses of gene expression in mammalian cells. Mol. Cell Proteomics3, 960–969. 10.1074/mcp.M400055-MCP200

  • 56

    VierstraR. D. (2003). The ubiquitin/26S proteasome pathway, the complex last chapter in the life of many plant proteins. Trends Plant Sci.8, 135–142. 10.1016/S1360-1385(03)00014-1

  • 57

    WangG.ZhuQ. G.MengQ. W.WuC. G. (2012). Transcript profiling during salt stress of young cotton (Gossypium hirsutum) seedlings via Solexa sequencing. Acta Physiol. Plant. 34, 107–115. 10.1007/s11738-011-0809-6

  • 58

    WangM. C.PengZ. Y.LiC. L.LiF.LiuC.XiaG. M. (2008). Proteomic analysis on a high salt tolerance introgression strain of Triticum aestivum/Thinopyrum ponticum. Proteomics8, 1470–1489. 10.1002/pmic.200700569

  • 59

    WangX.TangC.DengL.CaiG.LiuX.LiuB.et al. (2010). Characterization of a pathogenesis-related thaumatin-like protein gene TaPR5 from wheat induced by stripe rust fungus. Physiol. Plantarum139, 27–38. 10.1111/j.1399-3054.2009.01338.x

  • 60

    WangX.LiY.JiW.BaiX.CaiH.ZhuD.et al. (2011). A novel Glycine soja tonoplast intrinsic protein gene responds to abiotic stress and depresses salt and dehydration tolerance in transgenic Arabidopsis thaliana. J. Plant Physiol.168, 1241–1248. 10.1016/j.jplph.2011.01.016

  • 61

    WangZ. Q.XuX. Y.GongQ. Q.WeiC. X.LiF. J.ShanY. Q.et al. (2014). Root proteome of rice studied by iTRAQ provides integrated insight into aluminum stress tolerance mechanisms in plants. J. Proteomics98, 189–205. 10.1016/j.jprot.2013.12.023

  • 62

    WitzelK.WeidnerA.SurabhiG. K.BörnerA.MockH. P. (2009). Salt stress-induced alterations in the root proteome of barley genotypes with contrasting response towards salinity. J. Exp. Bot.60, 3545–3557. 10.1093/jxb/erp198

  • 63

    WuC. A.YangG. D.MengQ. W.ZhengC. C. (2004). The cotton GhNHX gene encoding a novel putative tonoplast Na+/H+ antiporter plays an important role in salt stress. Plant Cell Physiol.45, 600–607. 10.1093/pcp/pch071

  • 64

    WuT.TianZ.LiuJ.XieC. (2009). A novel leucine-rich repeat receptor-like kinase gene in potato, StLRPK1, is involved in response to diverse stresses. Mol. Biol. Rep.36, 2365–2374. 10.1007/s11033-009-9459-9

  • 65

    XiangY.HuangX.WangT.ZhangY.LiuQ. W.HusseyP. J.et al. (2007). Actin binding protein29 from lilium pollen plays an important role in dynamic actin remodeling. Plant Cell19, 1930–1946. 10.1105/tpc.106.048413

  • 66

    XuC.SibickyT.HuangB. (2010). Protein profile analysis of salt-responsive proteins in leaves and roots in two cultivars of creeping bentgrass differing in salinity tolerance. Plant Cell Rep.29, 595–615. 10.1007/s00299-010-0847-3

  • 67

    XuP.LiuZ.FanX.GaoJ.ZhangX.ZhangX.et al. (2013). De novo transcriptome sequencing and comparative analysis of differentially expressed genes in Gossypium aridum under salt stress. Gene525, 26–34. 10.1016/j.gene.2013.04.066

  • 68

    XueT.LiX.ZhuW.WuC.YangG.ZhengC. (2009). Cotton metallothionein GhMT3a, a reactive oxygen species scavenger, increased tolerance against abiotic stress in transgenic tobacco and yeast. J. Exp. Biol.60, 339–349. 10.1093/jxb/ern291

  • 69

    YangL. T.QiY. P.LuY. B.GuoP.SangW.FengH.et al. (2013). iTRAQ protein profile analysis of Citrus sinensis roots in response to long-term boron-deficiency. J. Proteomics93, 197–206. 10.1016/j.jprot.2013.04.025

  • 70

    YaoD.ZhangX.ZhaoX.LiuC.WangC.ZhangZ.et al. (2011). Transcriptome analysis reveals salt-stress-regulated biological processes and key pathways in roots of cotton (Gossypium hirsutum L.). Genomics98, 47–55. 10.1016/j.ygeno.2011.04.007

  • 71

    ZhangH.HanB.WangT.ChenS.LiH.ZhangY.et al. (2012). Mechanisms of plant salt response: insights from proteomics. J. Proteome Res.11, 49–67. 10.1021/pr200861w

  • 72

    ZhangL.XiD.LiS.GaoZ.ZhaoS.ShiJ.et al. (2011). A cotton group C MAP kinase gene, GhMPK2, positively regulates salt and drought tolerance in tobacco. Plant Mol. Biol. 77, 17–31. 10.1007/s11103-011-9788-7

  • 73

    ZhangX.LinC.ChenH.WangH.QuZ.ZhangH.et al. (2003). Cloning of a NaCl-induced fructose-1,6-diphosphate aldolase cDNA from Dualiella salinaand its expression in tobacco. Sci. China Life Sci.46, 49–57. 10.1007/BF03182684

  • 74

    ZhaoQ.ZhangH.WangT.ChenS.DaiS. (2013). Proteomics-based investigation of Salt-responsive mechanisms in plant roots. J. Proteomics82, 230–253. 10.1016/j.jprot.2013.01.024

  • 75

    ZhouY. J.GaoF.LiX. F.ZhangJ.ZhangG. F. (2010). Alterations in phosphoproteome under salt stress in Thellungiella roots. Chinese Sci. Bull.32, 3673–3679. 10.1007/s11434-010-4116-1

  • 76

    ZhuM.DaiS.McClungS.YanX.ChenS. (2009). Functional differentiation of brassica napus guard cells and mesophyll cells revealed by comparative proteomics. Mol. Cell Proteomics8, 752–766. 10.1074/mcp.M800343-MCP200

  • 77

    ZieskeL. R. (2006). A perspective on the use of iTRAQ™ reagent technology for protein complex and profiling studies. J. Exp. Bot.57, 1501–1508. 10.1093/jxb/erj168

  • 78

    ZörbC.SchmittS.MühlingK. H. (2010). Proteomic changes in maize roots after short-term adjustment to saline growth conditions. Proteomics10, 4441–4419. 10.1002/pmic.201000231

Summary

Keywords

Gossypium hirsutum, salt stress, iTRAQ, root, proteomics

Citation

Li W, Zhao F, Fang W, Xie D, Hou J, Yang X, Zhao Y, Tang Z, Nie L and Lv S (2015) Identification of early salt stress responsive proteins in seedling roots of upland cotton (Gossypium hirsutum L.) employing iTRAQ-based proteomic technique. Front. Plant Sci. 6:732. doi: 10.3389/fpls.2015.00732

Received

12 May 2015

Accepted

28 August 2015

Published

11 September 2015

Volume

6 - 2015

Edited by

Hans-Peter Mock, Institute of Plant Genetics and Crop Plant Reserach, Germany

Reviewed by

Keyvan Aghaei, The University of Zanjan, Iran; Giridara Kumar Surabhi, Regional Plant Resource Centre, India

Updates

Copyright

*Correspondence: Weiping Fang, Economic Crop Research Institute, Henan Academy of Agricultural Sciences, NO. 115, Huayuan Road, Zhengzhou 450002, China

†These authors have contributed equally to this work.

This article was submitted to Plant Proteomics, a section of the journal Frontiers in Plant Science

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics