Figures
Abstract
Salinity is a major abiotic stress that affects plant growth and development. In this study, we performed a proteomic analysis of cotton roots and leaf tissue following exposure to saline stress. 611 and 1477 proteins were differentially expressed in the roots and leaves, respectively. In the roots, 259 (42%) proteins were up-regulated and 352 (58%) were down-regulated. In the leaves, 748 (51%) proteins were up-regulated and 729 (49%) were down-regulated. On the basis of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, we concluded that the phenylalanine metabolism and starch and sucrose metabolism were active for energy homeostasis to cope with salt stress in cotton roots. Moreover, photosynthesis, pyruvate metabolism, glycolysis / gluconeogenesis, carbon fixation in photosynthetic organisms and phenylalanine metabolism were inhabited to reduce energy consumption. Characterization of the signaling pathways will help elucidate the mechanism activated by cotton in response to salt stress.
Citation: Chen T, Zhang L, Shang H, Liu S, Peng J, Gong W, et al. (2016) iTRAQ-Based Quantitative Proteomic Analysis of Cotton Roots and Leaves Reveals Pathways Associated with Salt Stress. PLoS ONE 11(2): e0148487. https://doi.org/10.1371/journal.pone.0148487
Editor: Wujun Ma, Murdoch University, AUSTRALIA
Received: August 13, 2015; Accepted: January 19, 2016; Published: February 3, 2016
Copyright: © 2016 Chen et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: This work was supported by the National Natural Science Foundation of China (31301262 and 31501343).
Competing interests: The authors have declared that no competing interests exist.
Introduction
Salinity is one of the major abiotic stresses that severely limit plant growth and development worldwide. It is estimated that more than 6% of the world’s total landmass and approximately 20% of irrigated landmass is affected by salinity[1]. Salt stress causes water deficit, ion toxicity, nutrient imbalance, and oxidative stress, leading to cellular damage and growth reduction, and even plant death [2–4].
Plant salinity tolerance is a complex phenomenon that involves physiological, biochemical, and molecular processes [1,5]. At the cellular level, plant roots function as the primary organ for sensing salinity so that they can respond rapidly to maintain functionality. Toxic ions are transported into the plant root along with the water stream, which moves from the soil to the vascular system of the root via several routes, including the symplastic and apoplastic pathways [2]. The ability of plants to maintain low cytosolic salt concentrations is due to their employment of selective ion uptake, ion exclusion[6] and compartmentalization; this is exemplified by the management of toxic levels of Na+ in vacuoles [7,8]. At the molecular level, several root proteins, such as plasma membrane and vacuolar H+-ATPase and Na+/H+ antiporters, play an essential role in ion uptake and transport [9,10]. In addition, antioxidant enzymes and soluble proteins in the cytoplasm have an important function in protecting cells from salt-induced oxidative stress[11,12].
Previously reported plant proteomics studies were limited to 2D gel electrophoresis analysis of roots [13,14]. However, the disadvantages of the 2D gel technique limited its application for comprehensive analysis of proteome changes [15–19]. Although most proteomics studies in crops have used 2D gel approaches [20], alternative methods are now available. Isobaric tags for relative and absolute quantitation (iTRAQ) is currently one of the most robust techniques that allows quantification of proteins on the basis of peptide labeling, as well as the identification and accurate quantification of proteins from multiple samples within broad dynamic ranges of protein abundance[21–24]. Moreover, some important post-translational modification information, such as proteolytic cleavage, glycosylation, or phosphorylation, may be retained. In addition, the proteome generated by iTRAQ typically consists of millions of peptides, from which substantial information on the proteins can be extracted with the use of bioinformatics tools.
Cotton (Gossypium hirsutum L.) produces an essential commodity, namely fibers for use in textiles, and cottonseed is a source of oil. Although it is classified as a salt-tolerant crop, this tolerance is actually limited, and varies according to the growth and developmental stages of the plant [25]. The effect of salinity on the germination, vegetative growth, and yield of cotton has been reported [26–28]. Breeders have sought to make cotton more tolerant to salt through various methods, including traditional plant breeding and biotechnological approaches such as creating transgenic cotton.
In this study, we used iTRAQ to assess proteome changes and identify proteins that were differentially expressed in cotton roots and leaves in response to NaCl concentration. On the basis of enrichment analysis of Gene Ontology (GO) annotations and Kyoto Encyclopedia of Genes and Genomes (KEGG), these differentially expressed proteins with functions in several biologically important pathways (phenylalanine metabolism and starch and sucrose metabolism for roots; photosynthesis, pyruvate metabolism, glycolysis / Gluconeogenesis, carbon fixation in photosynthetic organisms, and phenylalanine metabolism for leaves) we were likely important in the response to salt stress. Further characterization of these differentially expressed proteins will help to elucidate the signaling pathways activated by cotton in response to salt stress.
Materials and Method
Plant materials
Saline-tolerant seeds of Gossypium hirsutum cv. CCRI-79 were obtained from the National Medium-Term Gene Bank of Cotton in China and soaked in sterile deionized water at 28°C for 6 h. The method of seed planting was performed as described previously[29].
Salt treatment
Plants were cultured under non-saline conditions for 10 d to ensure full establishment before starting the salinity treatments. Salt stress treatment was initiated by providing the plants with full-strength Hoagland’s solutions containing 0 or 240 mM NaCl. To avoid osmotic shock, salt concentrations were increased daily by 40 mM NaCl, until reaching the required concentration. Three experimental replicates were performed. After 1 week, the roots and leaves of cotton were frozen immediately in liquid nitrogen and stored at −80°C.
Protein preparation and iTRAQ labeling
iTRAQ analysis was implemented at BGI(Shenzhen, China). Total proteins were extracted from the roots and leaves tissue of CCRI-79 plants using a phenol extraction procedure [30,31]. Two biological replicates were carried out for each sample.
Protein concentrations were determined using the Bradford method [32]. An equal amount of protein was prepared for each biological replicate. Protein samples from the roots and leaves were reduced with 10 mM DTT, alkylated with 55 mM iodoacetamide, digested using sequencing grade trypsin (Promega) at a ratio of 1:10 (w:w) for12 h at 37°C, and labeled using iTRAQ 8-plex kits (AB Sciex Inc., MA, USA) according to manufacturer’s manual. The root samples were labeled with iTRAQ tags 113 and114 (CK), 117 and 118 (salt treatment). The leaf samples were labeled with tags 115 and116 (CK), 119 and 120 (salt treatment).
Strong cation exchange
After labeling and quenching, the samples of roots and leaves were combined and lyophilized, respectively. The peptide mixture was dissolved in 4 mL strong cation exchange (SCX) buffer A (25% v/v acetonitrile, 25 mMNaH2PO4, pH 2.7). The peptides were fractionated on a Shimadzu LC-20ABHPLC system (Shimadzu, Kyoto, Japan) using an Ultremex SCX column (4.6 × 250 mm). Peptides were eluted at a flow rate of 1 mL/min with elution buffer B (25% v/v acetonitrile, 25 mM NaH2PO4, 1 M KCl, pH 2.7). The absorbance at 214 nm was monitored and 20 fractions were collected. Samples of each fraction were dried and desalted before LC-ESI MS/MS analysis [33,34].
LC-MS/MS analysis
The LC-MS/MS analysis was performed as described previously [33,34]. Peptides of each fraction (5-μL injections) were resolved insolvent A (5% acetonitrile, 0.1% formic acid) and centrifuged at 20,000 g for 10 min. The supernatant was separated using a Shimadzu LC-20AD Nano-HPLC system with a flow rate of 300 nL/min. Peptides were eluted by application of a linear gradient from 5% solvent B (95% acetonitrile v/v, 0.1% formic acid) to 35% solvent B over 35 min, followed by ramping up to 60% solvent B over 5 min, up to 80% in 2 min and maintained for 1 min; chromatographic conditions (5%) were restored in 1 min and equilibrated in solvent A for 10 min.
For the Triple TOF analysis, the AB SCIEX Triple TOF 5600 System (Concord, USA) was applied. Data were acquired using an ion spray voltage of 2.5 kV, curtain gas of 30 PSI, nebulizer gas of 15 PSI, and an interface heater temperature of 150°C. The MS was operated with a RP greater than or equal to 30,000 FWHM for TOFMS scans. For information-dependent data acquisition (IDA), survey scans were acquired in 250 ms and as many as 30 product ion scans were collected if exceeding a threshold of 120 counts per second (counts/s) and with a 2+ to 5+ charge-state. Total cycle time was fixed to 3.3 s. The Q2 transmission window was 100 Da for 100%. Four time bins were summed for each scan at a pulser frequency value of 11 kHz through monitoring of the 40 GHz multichannel TDC detector with four-anode/channel detection. A sweeping collision energy setting of 35 ± 5 eV, coupled with iTRAQ adjust rolling collision energy was applied to all precursor ions for collision-induced dissociation. Dynamic exclusion was set for 1/2 of peak width (18 s), and then the precursor was refreshed off the exclusion list. 5600 MS Converter was used to convert raw data files acquired from the Orbitrap into MGF files.
iTRAQ protein identification and quantification
Protein identification and quantification was performed using the Mascot 2.3.02 search engine (Matrix Science, Boston, MA). The protein mass were predicted by website http://www.expasy.ch/tools/ based on the protein sequences. The following search settings were used as described previously [35,36]. Searches were made against our cotton_AD_nr database, including protein sequences from the G. hirsutum L. genome (AADD) [37]. The search parameters were as follows: threshold set-off at 0.05 in the ion-score cutoff (with 95% confidence); MS/MS fragment ion mass tolerance of ±0.1 Da; enzyme specificity was set to trypsin with one missed cleavage; peptide tolerance was set at 10 ppm; fixed modifications of carbamido methylation at Cys and iTRAQ 8 plex at Lys and the N-terminal amino group of peptides; variable modifications of oxidation at Met, iTRAQ 8 plex at Tyr, and glutamine as pyroglutamic acid; peptide change was set at Mr and monoisotopic mass was chosen; charge states of peptides were set to +2 and +3. Only peptides with significance scores greater than “identity_ score” were counted as identified. Considering that multiple MS/MS spectra match to one peptide, normalization of the signal intensities of each MS/MS spectra was performed to find the most likely expression ratio for a given peptide. To demonstrate the reproducibility of the replicates, protein abundances between various biological replicates were compared, and ratios for each protein in each comparison were normalized to 1. The difference was plotted against the percentage of the proteins quantified. For quantitative changes, a 1.2-fold cutoff was set to determine up-regulated and down-regulated proteins, with a p-value < 0.05 present in at least two replicates[34].
Bioinformatic analysis of proteins
Functional category analysis was performed with Blast2GO software (http://www.geneontology.org)[38]. To take advantage of the current knowledge of biochemical pathways and other types of molecular interactions, we used KEGG databases (http://www.genome.jp/kegg/pathway.html)[39,40]. GO and pathway enrichment analysis were performed to determine which functional subcategories and metabolic pathways were overrepresented by the differentially accumulated proteins.
Results
Quantitative identification of cotton roots and leaves proteins using iTRAQ
In previous study, it was found that salt stress significantly reduced the growth rates, surface area, volume, and average diameter of the cotton roots and the dry weights of roots and leaves in cotton variety CCRI-79. To monitor the dynamic changes, cotton variety CCRI-79 was treated with 0 and 240 mM NaCl, respectively. A salt-induced proteomic analysis by iTRAQ was performed in roots and leaves. After merging data from the two replicates, a total of 404201 and 395831 spectra were generated from roots and leaves. Ultimately, we obtained 43828 unique spectra, 36382 identified peptides including 23299 unique peptides and 8974 identified proteins from roots. Meanwhile, we obtained 38679 unique spectra, 24395 identified peptides including 16920 unique peptides and 7008 identified proteins from leaves (Fig 1A). The protein mass was normally distributed. In roots, the number of 10–50, 50–100 and above 100 kDa protein is 4758 (53%), 3206 (36%) and 1010 (11%), respectively. In leaves, the number of 10–50, 50–100 and above 100 kDa protein is 4247 (60%), 2225 (32%) and 536 (8%), respectively (Fig 1B). The distribution of peptide number is shown in Fig 1C. In roots, the proteins with a single peptide, 2–5 peptides, 6–10 peptides and above 11 peptides comprise 4191, 3865, 663 and 156, respectively. In leaves, the proteins with a single peptide, 2–5 peptides, 6–10 peptides and above 11 peptides comprise 3560, 2942, 416 and 90, respectively (Fig 1C).
(A) Total spectra, spectra, unique spectra, peptides, unique peptide, and proteins identified from iTRAQ proteomic analysis. (B) Identified proteins were grouped according to the protein mass. (C) Number of peptides that match to proteins as indicated by MASCOT 2.3.02.
Functional categories of differentially expressed proteins
According to previous studies [15,41], differentially regulated proteins are defined based on a 1.2–1.5-fold change threshold. In this research, any proteins changed with ≥ 1.2-fold difference and a P value <0.05 would be designated as significant differently expression proteins (DEPs). Using these criteria, 611 and 1477 DEPs were identified in the roots and leaves, respectively (S1 Table). In the roots, 259(42%) proteins were up-regulated and 352 (58%) were down-regulated compared with the control. In the leaves, 748 (51%) proteins were up-regulated and 729 (49%) were down-regulated (Fig 2).
To gain insight into the functional categories that were altered following treatment with 240 mM NaCl, GO enrichment analysis of DEPs was conducted. The DEPs were classified into two groups, biological process and cellular component on the basis of GO enrichment analysis (Fig 3). In the roots, the main biological functional categories represented were metabolic process (356), organic substance metabolic process (284), primary metabolic process (268), and cellular metabolic process (251). For cellular components, proteins were predominantly distributed in the cell (306), cell part (306), intracellular part (271), intracellular organelle (241) and organelle (241)-related components (Fig 3A, S2 Table). In the leaves, the main biological functional categories were cellular metabolic process (527), single-organism metabolic process (404) and response to stimulus (274). According to the cellular components properties, these proteins were mainly classified into cell (716), cell part (716), intracellular (660) and intracellular part (649) (Fig 3B, S2 Table).
(A) GO classification of the DEPs in roots; (B) GO classification of the DEPs in leaves. The number of each category is displayed based on cellular components, molecular function, or biological process.
In order to analyze the metabolic pathways that responded to salt stress, 484 and 1100 DEPs were further investigated using the KEGG database in roots and leaves, respectively. These DEPs from roots were found to been enriched in biosynthesis of secondary metabolites (22.1%), ribosome (20.9%), phenylpropanoid biosynthesis (5.8%), phenylalanine metabolism (4.1%) and starch and sucrose metabolism (4.1%). The DEPs from leaves were found to been enriched in metabolic pathways (37.4%), biosynthesis of secondary metabolites (21.6%), ribosome (6.7%), photosynthesis (4.9), pyruvate metabolism (3.7%), glycolysis/Gluconeogenesis (3.7%), carbon fixation in photosynthetic organisms (3.4%), glutathione metabolism (3.0%) and phenylalanine metabolism (2.3%) (Table 1).
Carbon metabolism involved in the response to salt stress in cotton leaves.
According to our proteomic profile, carbon metabolism has important roles in response to salinity, particularly with regard to the contributions made by respiration and photosynthesis. The DEPs were conducted with an enrichment analysis and annotated in certain GO categories and KEGG pathways, respectively. The 453 annotated proteins were strongly enriched in GO categories “photosynthesis” (P = 3.06E-10), “photosynthesis, light reaction” (P = 7.09E-05), “electron transport chain” (P = 5.47E-07), “photosynthesis, light harvesting” (P = 3.20E-05), “oxidation-reduction process” (P = 0.00030) and “photosynthetic electron transport chain” (P = 0.00031). The 170 annotated proteins were significantly enriched in pathways of “photosynthesis” (P = 4.39E-13), “Pyruvate metabolism” (P = 0.0018), “Glycolysis / Gluconeogenesis” (P = 0.0116) and “Carbon fixation in photosynthetic organisms” (P = 0.0024). The chlorophyll and carotenoid contents were reduced. Simultaneously, the net photosynthetic rate (Pn) and stomatal conductance (gs) were inhibited by salinity due to NaCl in cotton variety CCRI-79 [29]. Consistent with this finding, 72 (7.7%) DEPs were enriched in the GO category “photosynthesis”. Among these proteins, 23 were down-regulated and enriched in the same pathway category, including chloroplast proteins in PSII (PsbA-E) and PSI (PsaA, PsaB, PsaF, PsaG, PsaL and PsaN), rieske fes protein (PetA and PetD) in cytochrome b6/f complex, F-type H+-transporting ATPase subunit b. The protein level of two light-harvesting proteins PSI PsaB and PSII PsbE were decreased by 0.4 to 1.2-fold in leaves (Fig 4, S3 Table). A number of proteins involved in “carbon fixation in photosynthetic organism” were down-regulated in leaves, including phosphoenolpyruvate carboxylase (PEPC), malate dehydrogenase, triosephosphate isomerase (TIM), fructose-1, 6-bisphosphatase I (FBP), transketolase (tak) and phosphoglycerate kinase. The protein level were decreased by 0.3 to 1.1-fold (Fig 4, S3 Table). This finding indicated that the photosynthetic capacity were inhibited under salt stress through the involvement of photosynthesis and carbon fixation.
Proteins involved in phenylalanine and starch and sucrose metabolism were mostly up-regulated in roots; Proteins involved in photosynthesis, carbon fixation, glycolysis, pyruvate metabolism and phenylalanine were mostly down-regulated in leaves. Red and green numbers represent expression level of up- and down-regulated proteins, respectively. PAL: phenylalanine ammonialyase; 4CL: 4-coumarate—CoA ligase; POD: peroxidase; UDPG: uridine diphosphate glucose; PGAL: glyceraldehyde 3-phosphate; PGA: 3-glyceric acid phosphate.
In the process of carbon metabolism, glycolysis and pyruvate metabolism were important as the part of respiration. Most of DEPs were enriched in the pathway category “glycolysis / gluconeogenesis”, including 2, 3-bisphosphoglycerate-dependent phosphoglycerate mutase, phosphoglucomutase and fructose-bisphosphate aldolase. Three protein level was decreased by 0.4 to 0.5-fold. Pyruvate kinase was enriched in the pathway category “pyruvate metabolism”, which was decreased by 1.4-fold (Fig 4, S3 Table). It was concluded that respiration were inhabited to reduce energy consumption in cotton leaves.
Starch and sucrose metabolism involved in the response to salt stress in roots.
Many proteins related to starch and sucrose metabolism were altered in response to increased salinity in cotton roots. Consistent with this finding, a number of proteins involved in starch and sucrose metabolism were up-regulated in roots (Fig 4, S3 Table). Such proteins included beta-fructofuranosidase, phosphoglucomutase (pgm), pectinesterase, polygalacturonase, trehalose 6-phosphate synthase (otsA), trehalose 6-phosphate phosphatase (otsB), beta-amylase, alpha-amylase and sucrose synthase. The protein level was increased by 0.2- to 0.7-fold (Fig 4, S3 Table). This finding indicated that the carbohydrate biosynthesis and metabolism was possibly developed for energy homeostasis to cope with salt stress in cotton roots.
The phenylalanine metabolism involved in the response to salt stress both in roots and leaves.
The phenylpropanoid pathway, which includes phenylalanine metabolism and phenylpropanoid biosynthesis, is involved in the stress response of plant cells. On the basis of the proteomic profile, in roots, PAL and 4-coumarate-CoA ligase (4CL) were up-regulated by 1.2 and 0.5-fold, respectively. In leaves, PAL and 4CL were down-regulated by 0.7 and 0.8-fold, respectively (S3 Table). Thus we concluded that phenylalanine metabolism was inhabited and active under salt stress in roots and leaves, respectively.
Discussion
To cope with salt stress, cotton plants have evolved complex salt-responsive signaling and metabolic processes at the cellular, organ, and whole-plant levels. Concomitant analysis of proteomics results from both the seedling roots and leaves of soybean has facilitated our understanding of how this plant responds to salt stress at the whole-plant level[42]. Several transcriptomic and proteomic analyses of the plant response to salinity have been performed [35,36,43]. In our study, we examined the molecular responses of cotton root and leaves to salt stress by iTRAQ technology. Our proteomics approach has provided a systematic comparison between the salt-responsive metabolic pathways found in cotton seedling leaves and roots under salt stress.
The growth rates of cotton plant roots and leaves decreases with increasing salt concentration [44,45], which may be result in osmotic injury or specific ion toxicity [46]. Compared with the no-salt control treatment, salt stress significantly reduces the growth rates, surface area, volume, and average diameter of the cotton roots, and the dry weights of roots and leaves in cotton variety CCRI-79 [45]. This is accompanied by strong changes in carbohydrate metabolism owing to severe impairments in the photosynthetic and respiration apparatus [47]. In the present study, we found evidence that leaf photosynthetic capacity was weakened under salt conditions, which indicates that this process is sacrificed during salt adaptation in order to conserve energy that then can be redirected to maintaining leaf growth [48]. This was underscored by the down-regulation of chloroplast proteins, an important enzyme involved in the first major photosynthetic step of photosynthesis.
Consistent with this observation, the DEPs were significantly enriched in the photosynthesis pathway. In light reaction. PSI, PSII, cytochrome b6/f and F-type H+-transporting ATP synthase subunits are involved. Photosynthesis takes place in the chloroplasts and the Chl content is important in the regulation and accumulation of PSI and PSII [49,50].
In our previous study, Zhang et al. found that the contents of Chl a, Chl b,and Chl (a+b) were decreased significantly during salinity stress in cotton CCRI-79 [45], which might directly reduce the abundance of these proteins and affect the accumulation of the chloroplast-targeted proteins. In our research, most of PSI and PSII proteins were down-regulated by 0.4- to 1.5-fold. It is reported that photosynthetic proteins were dynamically down-accumulation in mutant defected in Chl biosynthesis [51].
Respiration is an important metabolic pathway that is involved in carbohydrate metabolism in virtually all living organisms. The central role of glycolysis in plants is to provide energy in the form of ATP and to generate precursors such as fatty acids and amino acids for anabolism [48,52]. In this study, Glycolysis and pyruvate pathway were inhabited in leaves, accompanying 2, 3-bisphosphoglycerate-dependent phosphoglycerate mutase, phosphoglucomutase, fructose-bisphosphate aldolase and pyruvate kinase were down-regulated. Our results showed was consistent with the suggestion that increased respiration favors growth maintenance and adaptation to salt stress[48]. Together, these observations led us to conclude that increased energy production via glycolysis and pyruvate metabolism, together with a decreased photosynthetic capacity, maintains energy homeostasis. This adaptive balance permits normal growth of cotton plants even in the presence of salt stress.
Engagement of the phenylpropanoid pathway is one of the critical adaptations that allow plants to survive under stress conditions [53]. Phenylpropanoids are precursors to a wide range of phenolic compounds with many functions in plants. Lignin is a complex phenolic polymer that imparts strength, rigidity, and hydrophobicity to plant secondary cell walls[54]; its formation and subsequent deposition on the cell wall can affect the mechanical strength of supportive and water-conducting tissues[55]. In order to with stand various biotic and abiotic stresses, plants can trigger the production of lignin at specific sites [54,56]. Salt treatment induces root lignification, and increases the number of lignified vessels [57,58]. Lignin monomers are synthesized from phenylalanine[59]. PAL plays an essential role in the phenylpropanoid pathway and is responsive to both biotic and abiotic stress, including pathogen attack, wounding, cold, and UV light [60,61]. PAL catalyzes the first committed step of phenylpropanoid biosynthesis (the deamination of phenylalanine to cinnamic acid) and, along with 4CL, is essential for the synthesis of all phenylpropanoids [59].
Although LC-MS/MS techniques are instrumental in characterizing the peptides, quantitative peptides based biomarker discovery still remains challenging due to the several technological limitations. Among the emerging quantitative technologies, iTRAQ allows the concurrent protein sequence identification and relative quantification of those peptides with known protein sequences in up to eight different biological samples in single experiments [62]. However, due to its limited throughput and current cost, iTRAQ is not feasible to simultaneously compare large sample sizes of salinity subjects to achieve the discovery of differential features of sufficient statistical power [63]. In addition, the success of iTRAQ efforts depends on the peptide sequence determination. Despite our increased understanding of the peptide composition, peptide sequence characterization by the cotton AD genome annotation database (81147 sequences) and the National Center for Biotechnology Information (NCBI) non-redundant fasta database (6833826 sequences) combined, could only determine a small portion peptide features revealed by the HPLC or CE coupled mass spectrometric analysis, which leads to undersampling and incomplete analytic coverage of the peptide. Moreover, it is difficult to meet with the following objectives: to reduce background interference and to increase the signal-to-noise ratio [64]. Therefore, it is necessary to address clearly the limits of iTRAQ and the sources of departure from the theoretical values before using, in order to take full advantage of iTRAQ’s capability in discovering important biomarkers for salt tolerant cotton breeding in the future.
Conclusions
A proteomic analysis by iTRAQ technology was performed to analyze differential protein expression under salt stress in cotton roots and leaves. In total, 611 and 1477 differentially expressed proteins were identified in roots and leaves, respectively. On the basis of GO and pathway enrichment analysis, we concluded that the “phenylalanine metabolism” and “starch and sucrose metabolism” were active for energy homeostasis to cope with salt stress in cotton roots. Moreover, “photosynthesis”, “pyruvate metabolism”, “glycolysis / Gluconeogenesis”, “carbon fixation in photosynthetic organisms”, and “phenylalanine metabolism” were inhabited to reduce energy consumption in cotton leaves. Further studies for gene function analysis were needed to further clarify the molecular mechanism under salt stress in cotton roots and leaves.
Ethical standards
The authors note that this research is performed and reported in accordance with ethical standards of the scientific conduct.
Supporting Information
S1 Table. Differentially expressed proteins identified in roots and leaves.
https://doi.org/10.1371/journal.pone.0148487.s001
(XLSX)
S2 Table. The lists of DEPs in GO enrichment analysis.
https://doi.org/10.1371/journal.pone.0148487.s002
(XLSX)
S3 Table. The expression level of DEPs in metabolism pathway.
https://doi.org/10.1371/journal.pone.0148487.s003
(XLSX)
Acknowledgments
This work was supported by the National Natural Science Foundation of China (31301262 and 31501343).
Author Contributions
Conceived and designed the experiments: TTC LZ XHZ YLY. Performed the experiments: JWL JWG QG AYL. Analyzed the data: HHS YZS WKG JP. Contributed reagents/materials/analysis tools: SDL SPZ HJM. Wrote the paper: TTC LZ XHZ YLY.
References
- 1. Munns R, Tester M. Mechanisms of salinity tolerance. Annual Review of Plant Biology. 2008; 59: 651–681. pmid:18444910
- 2. Tester M, Davenport R. Na+ tolerance and Na+ transport in higher plants. Annals of Botany. 2003; 91: 503–527. pmid:12646496
- 3. Nohzadeh Malakshah S, Habibi Rezaei M, Heidari M, Hosseini Salekdeh G. Proteomics reveals new salt responsive proteins associated with rice plasma membrane. Bioscience, Biotechnology, and Biochemistry. 2007; 71: 2144–2154. pmid:17827676
- 4. Munns R. Genes and salt tolerance: bringing them together. New Phytologist. 2005; 167: 645–663. pmid:16101905
- 5. Greenway H, Munns R. Mechanisnn of salt tolerance in non-halophytes. Annu Rev Plant Physiology. 1980: 149–190.
- 6. Shi H, Ishitani M, Kim C, Zhu J-K. The Arabidopsis thaliana salt tolerance gene SOS1 encodes a putative Na+/H+ antiporter. Proceedings of the National Academy of Sciences. 2000; 97: 6896–6901.
- 7. Blumwald E, Poole RJ. Salt tolerance in suspension cultures of sugar beet: induction of Na+/H+ antiport activity at the tonoplast by growth in salt. Plant Physiology. 1987; 83: 884–887. pmid:16665356
- 8. Apse MP, Aharon GS, Snedden WA, Blumwald E. Salt tolerance conferred by overexpression of a vacuolar Na+/H+ antiport in Arabidopsis. Science. 1999; 285: 1256–1258. pmid:10455050
- 9. Munns R. Comparative physiology of salt and water stress. Plant, Cell & Environment. 2002; 25: 239–250.
- 10. Aharon G, Apse M, Duan S, Hua X, Blumwald E. Characterization of a family of vacuolar Na+/H+ antiporters in Arabidopsis thaliana. Plant and Soil. 2003; 253: 245–256.
- 11. Gossett DR, Millhollon EP, Lucas MC. Antioxidant response to NaCl stress in salt-tolerant and salt-sensitive cultivars of cotton. Crop Science. 1994; 34: 706–714.
- 12. Marcum KB, Murdoch CL. Salinity tolerance mechanisms of six C4 turfgrasses. Journal of the American Society for Horticultural Science. 1994; 119: 779–784.
- 13. Gao L, Yan X, Li X, Guo G, Hu Y, Ma W, et al. Proteome analysis of wheat leaf under salt stress by two-dimensional difference gel electrophoresis (2D-DIGE). Phytochemistry. 2011; 72: 1180–1191. pmid:21257186
- 14. Song Y, Zhang C, Ge W, Zhang Y, Burlingame AL, Guo Y. Identification of NaCl stress-responsive apoplastic proteins in rice shoot stems by 2D-DIGE. Journal of Proteomics. 2011; 74: 1045–1067. pmid:21420516
- 15. Lan P, Li W, Wen T-N, Shiau J-Y, Wu Y-C, Lin W, et al. iTRAQ protein profile analysis of Arabidopsis roots reveals new aspects critical for iron homeostasis. Plant Physiology. 2011; 155: 821–834. pmid:21173025
- 16. Chandramouli K, Qian P-Y. Proteomics: challenges, techniques and possibilities to overcome biological sample complexity. Human Genomics and Proteomics. 2009; 2009:239–242.
- 17. Owiti J, Grossmann J, Gehrig P, Dessimoz C, Laloi C, Hansen MB, et al. iTRAQ-based analysis of changes in the cassava root proteome reveals pathways associated with post-harvest physiological deterioration. The Plant Journal. 2011; 67: 145–156. pmid:21435052
- 18. Wienkoop S, Weckwerth W. Relative and absolute quantitative shotgun proteomics: targeting low-abundance proteins in Arabidopsis thaliana. Journal of Experimental Botany. 2006; 57: 1529–1535. pmid:16574747
- 19. Schulze WX, Usadel B. Quantitation in mass-spectrometry-based proteomics. Annual Review of Plant Biology. 2010; 61: 491–516. pmid:20192741
- 20. Salekdeh GH, Komatsu S. Crop proteomics: aim at sustainable agriculture of tomorrow. PROTEOMICS. 2007; 7: 2976–2996. pmid:17639607
- 21. Casado-Vela J, Martínez-Esteso MJ, Rodriguez E, Borrás E, Elortza F, Bru-Martínez R. iTRAQ-based quantitative analysis of protein mixtures with large fold change and dynamic range. PROTEOMICS. 2010; 10: 343–347. pmid:20029838
- 22. Nogueira FCS, Palmisano G, Schwämmle V, Campos FAP, Larsen MR, Domont GB, et al. Performance of isobaric and isotopic labeling in quantitative plant proteomics. Journal of Proteome Research. 2012; 11: 3046–3052. pmid:22452248
- 23. Jorrin-Novo JV, Komatsu S, Weckwerth W, Wienkoop S, Martínez-Esteso M, Casado-Vela J, et al. Differential plant proteome analysis by isobaric tags for relative and absolute quantitation (iTRAQ). Plant Proteomics: Humana Press. 2013.
- 24. Ross PL, Huang YN, Marchese JN, Williamson B, Parker K, Hattan S, et al. Multiplexed protein quantitation in Saccharomyces cerevisiae using amine-reactive isobaric tagging reagents. Molecular & Cellular Proteomics. 2004; 3: 1154–1169.
- 25. Qidar M, Shams M. Some agronomic and physiological aspects of salt tolerance in cotton (Gossypium hirsutum L.). Journal of Agronomy and Crop Science. 1997; 179: 101–106.
- 26. Guo WX, Mass SJ, Bronson KF. Relationship between cotton yield and soil eletrical conductivity, topography, and Landsat imagery. Precision Agronomy. 2012; 13: 678–692.
- 27. Zhang L, Zhang GW, Wang YH, Zhou ZG, Meng YL, Chen BL. Effect of soil salinity on physiological characteristics of functional leaves of cotton plants. Journal of Plant Research. 2013; 126: 293–304. pmid:23114969
- 28. Ahmads S, Khan N, Iqbal MZ. Salt tolerance of cotton (Gossypium hirsutum L.). Asian Journal of Plant Science. 2002; 1: 715–719.
- 29. Zhang L, Ma H, Chen T, Pen J, Yu S, Zhao X. Morphological and physiological responses of cotton (Gossypium hirsutum L.) plants to salinity. PLoS One. 2014; 9: e112807. pmid:25391141
- 30. Saravanan RS, Rose JKC. A critical evaluation of sample extraction techniques for enhanced proteomic analysis of recalcitrant plant tissues. PROTEOMICS. 2004; 4: 2522–2532. pmid:15352226
- 31. Yao Y, Yang Y-W, Liu J-Y. An efficient protein preparation for proteomic analysis of developing cotton fibers by 2-DE. ELECTROPHORESIS. 2006; 27: 4559–4569. pmid:17058302
- 32. Kruger N. The Bradford Method for Protein Quantitation. Basic Protein and Peptide Protocols: Humana Press. 1994.
- 33. Liu J, Pang C, Wei H, Song M, Meng Y, Ma J, et al. iTRAQ-facilitated proteomic profiling of anthers from a photosensitive male sterile mutant and wild-type cotton (Gossypium hirsutum L.). Journal of Proteomics. 2015; 126: 68–81. pmid:26047712
- 34. Chu P, Yan GX, Yang Q, Zhai LN, Zhang C, Zhang FQ, et al. iTRAQ-based quantitative proteomics analysis of Brassica napus leaves reveals pathways associated with chlorophyll deficiency. Journal of Proteomics. 2015; 113: 244–259. pmid:25317966
- 35. Ge X, Zhang C, Wang Q, Yang Z, Wang Y, Zhang X, et al. iTRAQ protein profile differential analysis between somatic globular and cotyledonary embryos reveals stress, hormone, and respiration involved in increasing plantlet regeneration of Gossypium hirsutum L. Journal of Proteome Research. 2015; 14: 268–278. pmid:25367710
- 36. Yang L-T, Qi Y-P, Lu Y-B, Guo P, Sang W, Feng H, et al. iTRAQ protein profile analysis of Citrus sinensis roots in response to long-term boron-deficiency. Journal of Proteomics. 2013; 93: 179–206. pmid:23628855
- 37. Li F, Fan G, Lu C, Xiao G, Zou C, Kohel RJ, et al. Genome sequence of cultivated Upland cotton (Gossypium hirsutum TM-1) provides insights into genome evolution. Nat Biotech. 2015; 33: 524–530.
- 38. Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005; 21: 3674–3676. pmid:16081474
- 39. Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Research. 2012.
- 40. Kanehisa M, Goto S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Research. 2000; 28: 27–30. pmid:10592173
- 41. Zhu M, Dai S, McClung S, Yan X, Chen S. Functional differentiation of Brassica napus guard cells and mesophyll cells revealed by comparative proteomics. Molecular & Cellular Proteomics. 2009; 8: 752–766.
- 42. Zhao Q, Zhang H, Wang T, Chen S, Dai S. Proteomics-based investigation of salt-responsive mechanisms in plant roots. Journal of Proteomics. 2013; 82: 230–253. pmid:23385356
- 43. Xie F, Stewart CN, Taki FA, He Q, Liu H, Zhang B. High-throughput deep sequencing shows that microRNAs play important roles in switchgrass responses to drought and salinity stress. Plant Biotechnology Journal. 2014; 12: 354–366. pmid:24283289
- 44. Garratt LC, Janagoudar BS, Lowe KC, Anthony P, Power JB, Davey MR. Salinity tolerance and antioxidant status in cotton cultures. Free Radical Biology and Medicine. 2002; 33: 502–511. pmid:12160932
- 45. Zhang L, Ma H, Chen T, Pen J, Yu S, Zhao X. Morphological and physiological responses of cotton (Gossypium hirsutum L.) plants to salinity. PLoS One. 2014; 9: e112807. pmid:25391141
- 46. Meloni DA, Oliva MA, Ruiz HA, Martinez CA. Contribution of proline and inorganic solutes to osmotic adjustment in cotton under salt stress. Journal of Plant Nutrition. 2001; 24: 599–612.
- 47. Singh A, Dubey R. Changes in chlorophyll a and b contents and activities of photosystems I and II in rice seedlings induced by NaCl. Photosynthetica 1995; 31: 489–499.
- 48. Shi H, Ye T, Chan Z. Comparative proteomic responses of two bermudagrass (Cynodon dactylon (L). Pers.) varieties contrasting in drought stress resistance. Plant Physiology and Biochemistry. 2014; 82: 218–228. pmid:24992888
- 49. Leister D. Chloroplast research in the genomic age. Trends in Genetics. 2003; 19: 47–56. pmid:12493248
- 50. Hansson A, Jensen PE. Chlorophyll limitation in plants remodels and balances the photosynthetic apparatus by changing the accumulation of photosystems I and II through two different approaches. Physiol Plant 2009; 135: 214–228. pmid:19055541
- 51. Zhou Y, Gong ZY, Yang ZF, Yuan Y, Zhu JY, Wang M, et al. Mutation of the Light-Induced Yellow Leaf 1 Gene, Which Encodes a Geranylgeranyl Reductase, Affects Chlorophyll Biosynthesis and Light Sensitivity in Rice. PloS One. 2013; 8: e75299. 75210.71371/journal.pone.0075299. pmid:24058671
- 52. Plaxton WC. The organization and regulation of plant glycolysis. Annual Review of Plant Physiology and Plant Molecular Biology. 1996; 47: 185–214. pmid:15012287
- 53. Ferrer JL, Austin MB, Stewart C Jr, Noel JP. Structure and function of enzymes involved in the biosynthesis of phenylpropanoids. Plant Physiology and Biochemistry. 2008; 46: 356–370. pmid:18272377
- 54. Rogers LA, Campbell MM. The genetic control of lignin deposition during plant growth and development. New Phytologist. 2004; 164: 17–30.
- 55. Degenhardt B, Gimmler H. Cell wall adaptations to multiple environmental stresses in maize roots. Journal of Experimental Botany. 2000; 51: 595–603. pmid:10938816
- 56. Lange BM, Lapierre C, Sandermann H Jr. Elicitor-induced spruce stress lignin (structural similarity to early developmental lignins). Plant Physiology. 1995; 108: 1277–1287. pmid:12228544
- 57. Sánchez-Aguayo I, Rodríguez-Galán JM, García R, Torreblanca J, Pardo JM. Salt stress enhances xylem development and expression of S-adenosyl-l-methionine synthase in lignifying tissues of tomato plants. Planta. 2004; 220: 278–285. pmid:15322882
- 58. Neves GYS, Marchiosi R, Ferrarese MLL, Siqueira-Soares RC, Ferrarese-Filho O. Root growth inhibition and lignification induced by salt Stress in soybean. Journal of Agronomy and Crop Science. 2010; 196: 467–473.
- 59. Bonawitz ND, Chapple C. The genetics of lignin biosynthesis: connecting genotype to phenotype. Annual Review of Genetics. 2010; 44: 337–363. pmid:20809799
- 60. Xu L, Zhu L, Tu L, Liu L, Yuan D, Jin L, et al. Lignin metabolism has a central role in the resistance of cotton to the wilt fungus Verticillium dahliae as revealed by RNA-Seq-dependent transcriptional analysis and histochemistry. Journal of Experimental Botany. 2011; 62: 5607–5621. pmid:21862479
- 61. Huang J, Gu M, Lai Z, Fan B, Shi K, Zhou Y-H, et al. Functional analysis of the Arabidopsis PAL gene family in plant growth, development, and response to environmental stress. Plant Physiology. 2010; 153: 1526–1538. pmid:20566705
- 62. Ong SE, Mann M. Mass spectrometry-based proteomics turns quantitative. Nature Chemical Biology. 2005: 252–262. pmid:16408053
- 63.
Makowski GS. Advances in clinical chemistry. Academic Press. 2010: 188.
- 64. Ow SY, Salim M, Noirel J, Evans C, Rehman I, Wright PC. iTRAQ underestimation in simple and complex mixtures: “The Good, the Bad and the Ugly. Journal of Proteome Research. 2009; 8: 5347–5355. pmid:19754192