Dissecting the proteome dynamics of the salt stress induced changes in the leaf of diploid and autotetraploid Paulownia fortunei

Abstract

Exposure to high salinity can trigger acclimation in many plants. Such an adaptative response is greatly advantageous for plants and involves extensive reprogramming at the molecular level. Acclimation allows plants to survive in environments that are prone to increasing salinity. In this study, diploid and autotetraploid Paulownia fortunei seedlings were used to detect alterations in leaf proteins in plants under salt stress. Up to 152 differentially abundant proteins were identified by Multiplex run iTRAQ-based quantitative proteomic and LC-MS/MS methods. Bioinformatics analysis suggested that P. fortunei leaves reacted to salt stress through a combination of common responses, such as induced metabolism, signal transduction, and regulation of transcription. This study offers a better understanding of the mechanisms of salt tolerance in P. fortunei and provides a list of potential target genes that could be engineered for salt acclimation in plants, especially trees.

About the Authors:

Minjie Deng

Contributed equally to this work with: Minjie Deng, Yanpeng Dong

Roles Supervision, Writing – original draft

Affiliations College of Forestry, Henan Agricultural University, Zhengzhou, Henan, P.R. China, Institute of Paulownia, Henan Agricultural University, Zhengzhou, Henan, P.R. China

Yanpeng Dong

Contributed equally to this work with: Minjie Deng, Yanpeng Dong

Roles Writing – original draft

Affiliations College of Forestry, Henan Agricultural University, Zhengzhou, Henan, P.R. China, Institute of Paulownia, Henan Agricultural University, Zhengzhou, Henan, P.R. China

Zhenli Zhao

Roles Methodology, Software

Affiliations College of Forestry, Henan Agricultural University, Zhengzhou, Henan, P.R. China, Institute of Paulownia, Henan Agricultural University, Zhengzhou, Henan, P.R. China

Yongsheng Li

Roles Methodology, Software

Affiliations College of Forestry, Henan Agricultural University, Zhengzhou, Henan, P.R. China, Institute of Paulownia, Henan Agricultural University, Zhengzhou, Henan, P.R. China

Guoqiang Fan

Roles Conceptualization, Funding acquisition, Supervision, Writing – review & editing

* E-mail: [email protected]

Affiliations College of Forestry, Henan Agricultural University, Zhengzhou, Henan, P.R. China, Institute of Paulownia, Henan Agricultural University, Zhengzhou, Henan, P.R. China

ORCID http://orcid.org/0000-0003-0385-2578

Introduction

Salinization is one of the most adverse environmental stresses affecting plant growth and crop productivity in many areas of the globe, especially in already susceptible regions [1]. Plants have adapted dynamic responses to salt stress at the physiological, biochemical, and molecular levels to help in their survival [2].

In some plants, saltiness triggers the expression of a variety of genes that function in the plant’s salt response as well as in gene expression regulation and signal transduction to increase salt tolerance. Understanding how plants behave under such conditions and applying the knowledge to improve their performance in these situations is of interest to plant scientists. In recent decades, several advances have been made towards identifying potential salt stress-related genes that are capable of increasing the tolerance of plants to salt stresses [3, 4]. It was found that in this process different biochemical and molecular mechanisms act together in a complex regulatory network. These interactions are controlled by hormone signaling pathways that can induce or antagonize one another; in particular, the pathways that involve abscisic acid [5]. Transcription factors [6], kinase cascades [7], reactive oxygen species (ROS) [8], heat shock factors [9], and small RNAs [10] are key components of the molecular crosstalk in these pathways.

High-throughput genomic and transcriptomic sequencing techniques are powerful tools that have been used to gain comprehensive genomic information of the novelty and plasticity underlying the adaptation response of plants under salt stress. Many studies have focused on the transcriptional reprogramming of the related biological processes [11–13]. However, because most biological functions are carried out by proteins, and posttranscriptional regulation and posttranslational modifications can influence a system, quantitation of the proteins in a system is mandatory to understand the biological process and functions involved [14, 15]. Proteomics is a powerful tool for the identification of proteins and mechanisms involved in the salt stress response and tolerance [16–19]. Salt stress responsive proteins include those involved in basic metabolic pathways and lignification. Researches on drought stress in rice leaves suggested that 42 proteins, which were induced in mild drought stress, were completely reversed within 10 days of rewatering. Plants that have adapted to prolonged drought stress are capable of recovering more completely than unadapted plants after rewatering [20]. The concentration of proline, a compatible solute and also an osmoprotectant, was found to increase by up to 20-fold in some plants under severe drought conditions [21, 22]. Proline may protect or induce stress protective proteins and thus confer a protective effect. It was reported that severe salt stress inhibited the antioxidative enzymes catalase and peroxidase, and that their activity was maintained at significantly higher levels in the presence of proline in Pancratium maritimum [23].

Proteomics data for Paulownia species are particularly desirable because some of these species are of economic and forestry importance. Several transcriptome studies of Paulownia have been reported [24–27]. However, up to now, there is little Paulownia protein information in the public databases, despite the availability of large numbers of transcriptome unigenes. To investigate the salt stress response of the proteome in Paulownia, in this study, we conducted a systematic proteomic analysis during the plants’ acclimation to salt stress. Proteins from control and salt-treated Paulownia fortunei were identified and common proteins associated within this process were detected. A total of 152 proteins that responded to salt stress were identified by mass spectrometry. As a result, we established a comprehensive inventory of salinity-responsive proteins, many of which were involved in functions that included photosynthetic processes, ROS scavenging, signaling, and stress defense. We found that polyploidization affected the salt-responsive pattern, especially that of ion channel proteins. The autotetraploid P. fortunei had more salt-responsive proteins than the diploid P. fortunei. To our knowledge, this is the first reported proteome study on salt-treated Paulownia. This analysis provides a deeper insight of the acclimation process, which will help to understand the general stress response to saltiness and to identify pivotal proteins involved in plant stress responses.

Materials and methods

Plant material

All biological materials in this study were from the Institute of Paulownia, Henan Agricultural University, Zhengzhou, Henan Province, China.

Uniformly grown tissue cultured seedlings of diploid (PF2) and autotetraploid P. fortunei (PF4) were cultured in 100-mL triangular flasks on 1/2 MS medium for 30 days. After that the seedlings were clipped from the roots. Samples with the same height and crown size were then transferred into plastic pots containing the same amount of ordinary garden soil, one pot for each plant.

Prior to salt stress treatment, the plants were housed randomly in an outdoor nursery for 50 days. Twelve plants from each accession were selected, including three controls and nine plants for salt treatment. For salt treatment, the PF2 and PF4 seedlings in the three groups were treated with 0.4% NaCl solution (1 L) every two days for 15 days. The salt-treated PF2 and PF4 samples were named PFS2 and PFS4, respectively. The PF2 and PF4 seedlings in the control group were irrigated with tap water every two days for 15 days. Healthy leaves (second leaf from the apex) were harvested at 15 days, and equal numbers of leaves from each plant in each of the replicate groups were pooled to form the four samples. Terminal buds of 1.5 cm in length of the plantlets in the four treatment groups were sheared, immediately frozen in liquid nitrogen, and stored at −80°C.

Protein extraction, digestion, and iTRAQ labeling

The PF2, PF4, PFS2, and PFS4 samples were ground into powder in liquid nitrogen, extracted with lysis buffer containing 1 mM PMSF and 2 mM EDTA. After that, proteins were extracted according to the method of Tang et al. [28]. Total protein was digested and iTRAQ-labelled according to Meng et al. [29]. Two biological replicates were set for each sample. Eight samples were labeled with iTRAQ tags as follows: PF2-1, 113; PF2-2, 117; PF4-1, 114, PF4-2, 118; PFS2-1, 119, PFS2-2, 115; PFS4-1, 116, PFS4-2, 121. The peptides were labeled with the isobaric tags and incubated at room temperature for 2 h. The labeled peptide mixtures were then pooled for strong cation exchange and dried by vacuum centrifugation.

Strong cation exchange

Strong cation exchange chromatography was performed with a LC-20AB HPLC pump system. The iTRAQ-labeled peptide mixtures were reconstituted with 4 mL buffer A and loaded onto a 4.6 × 250 mm Ultremex SCX column containing 5-μm particles. The peptides were eluted at a flow rate of 1 mL min⁻¹ with a gradient of buffer A for 10 min, 5–60% buffer B for 27 min, and 60–100% buffer B for 1 min. The system was then maintained at 100% buffer B for 1 min before equilibrating with buffer A for 10 min prior to the next injection. Elution was monitored by measuring the absorbance at 214 nm, and fractions were collected every 1 min. The eluted peptides were pooled into 20 fractions, desalted with a Strata X C18 column, and vacuum-dried.

Liquid chromatography tandem-mass spectrometry (LC-MS/MS) analysis

Mass spectroscopy analysis was performed using an AB SCIEX TripleTOF™ 5600 mass spectrometer, coupled with an online micro-flow HPLC system as described in section 1.3. The peptides were separated using the method of Qiao et al. [30].

iTRAQ protein identification and quantification

The acquired raw data files were converted into Mascot generic format files according the method of Lin et al. [31]. The Mascot 2.3.02 (Matrix Science, London, United Kingdom) was used for protein identification and quantification according the method of Guo et al. [32]. Protein identification was performed using the Mascot search engine against a Paulownia transcriptome database that contained 82,934 sequences. The sequencing data have been submitted to the Short Reads Archive under accession number SRP034738 [26]. To reduce the probability of false peptide identification, only peptides with significance scores in the 99% confidence interval greater than “identity” in a Mascot probability analysis were counted as identified. Each confident protein identification involved at least one unique peptide.

To reduce the probability of false peptide identification, only peptides with significance scores in the 99% confidence interval greater than “identity” in a Mascot probability analysis were counted as identified. Each confident protein identification involved at least one unique peptide.

For protein quantization, each protein was required to contain at least two unique peptides. The quantitative protein ratios were weighted and normalized by the median ratio in Mascot. The student’s T-test integrated into the Mascot software was used. Only the proteins with ratios>1.2or<0.84 and p-values<0.05 were considered as differentially abundant proteins (DAPs).

Bioinformatics analysis

Functional analysis of the identified proteins was conducted using Gene Ontology (GO) annotations and the proteins were assigned terms under the biological process, molecular function, and cellular component categories. The proteins were further analyzed using the Clusters of Orthologous Groups (COG) of proteins database and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database [33].

RNA preparation and quantitative RT-PCR

In this analysis, the salt-treated PF2 samples were named PFS2-1 and PFS2-2, which corresponded to PF2 treated with salt for 15 days and 30 days, respectively. PFS4-1 and PFS4-2 were named in the same way. The RNA from these samples were extracted with Trizol (Invitrogen, CA, USA). The RNA samples from PF2, PF4, PFS2-1, PFS2-2, PFS4-1, and PFS4-2 were extracted with Trizol. The RNA was then precipitated with isopropanol. Purified and concentrated RNA was denatured and first-strand cDNAs for all the samples were synthesized using a PrimeScript RT reagent kit (Takara Biotechnology, Dalian, China). The genes encoding eight DAPs were selected for qRT-PCR analysis. Primers were shown in Table 1, 18S rDNA was selected as the reference gene. The method of qRT-PCR and relative expression calculation was described by Fan et al. [10].

[Figure omitted. See PDF.]

Table 1. Primers of quantitative RT-PCR analysis of candidate DAP genes.

https://doi.org/10.1371/journal.pone.0181937.t001

Results

Proteomics characterization

Proteins were extracted from control and salt-treated diploid and autotetraploid P. fortunei (PF2, PFS2, PF4, and PFS4) seedling leaves and a total of 312,926 spectra were generated. 34,608 spectra that matched known spectra were identified by the Mascot software. Among them, 23,881 unique spectra were matched to 7,040 unique peptides and 2,634 proteins (Fig 1A), approximately 66% of which consisted of at least two of the unique peptides (Fig 1B). The vast majority of these proteins were larger than 10 kDa, although their molecular weights covered a wide range (Fig 1C). Most of the identified proteins had good peptide coverage; 45% had more than 10% sequence coverage, and 22% had 20% sequence coverage (Fig 1D).

[Figure omitted. See PDF.]

Fig 1. Summary of the iTRAQ-based Paulownia fortunei proteome.

(a) Spectra, peptides and proteins that were identified by searching against the Paulownia fortunei transcriptome database; (b) Number of peptides that match proteins using MASCOT; (c) Distribution of the proteins that were identified among different molecular weights; (d) Coverage of the proteins by the identified peptides.

https://doi.org/10.1371/journal.pone.0181937.g001

To avoid identification omissions, less than 2 ppm were confined as the database search strategy for the peptide matching error. The reproducibility of the proteomic analysis was shown in Fig 2. These results indicate that our proteomics analyses were reliable.

[Figure omitted. See PDF.]

Fig 2. Repeatability of the expression of duplicate samples.

(a) PF2; (b) PF4; (c) PFS2; (d) PFS4. X-axis represents the difference of the quantitative ratios between the first and the second biological repeats of the two samples. The right y-axis represents the cumulative percentage between the proteins of a certain range and the quantitative proteins, while the left y-axis represents the number of total protein in a certain range.

https://doi.org/10.1371/journal.pone.0181937.g002

To carry out a functional analysis, we assigned gene ontology (GO) terms to all the quantitated proteins. Under the biological process category, metabolic process (1,305) and cellular process (1,003) were the most represented groups; under the cellular component category, cell (1,040) and cell part (1,040) were the two largest groups; and under molecular function category, catalytic activity (1,153) and binding (900) were most represented (Fig 3). These results indicate that the identified proteins are involved in almost every aspect of P. fortunei metabolism.

[Figure omitted. See PDF.]

Fig 3. Gene Ontology classification of distinct proteins that were detected in Paulownia fortunei leaves.

https://doi.org/10.1371/journal.pone.0181937.g003

To further understand the functions of the 2,634 proteins, they were assigned to 24 categories in the Clusters of Orthologous Groups (COG) database. The most highly represented functional categories were “general function prediction only” and “posttranslational modification, protein turnover, chaperones”, which represented approximately 14% and 12% of the identified proteins, respectively (Fig 4).

[Figure omitted. See PDF.]

Fig 4. Classification of the clusters of orthologous groups function classification of distinct proteins that were detected in Paulownia fortunei leaves.

https://doi.org/10.1371/journal.pone.0181937.g004

The annotated proteins were mapped onto 128 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways (S1 Table). Among them, “metabolic pathways” (830) was significantly more highly represented than the other pathways, followed by “biosynthesis of secondary metabolites” (452), and “carbon metabolism” (192). Based on these analyses, we concluded that most of the mapped proteins may affect cellular component biogenesis, posttranslational modification, protein turnover, and response to stimulus.

Analysis of the differentially abundant proteins related to salt stress in the autotetraploid

Differentially abundant proteins (DAPs) were defined as the proteins that showed a >1.2-fold or <0.84-fold change in relative abundance and a p-value <0.05. A total of 767 DAPs were identified in the PF4 vs. PF2 comparison; 360 were more abundant and 407 were less abundant. A total of 916 DAPs were identified in the PFS2 vs. PF2 comparison; 488 were more abundant and 428 were less abundant. A total of 712 DAPs were identified in the PFS4 vs. PF4 comparison; 390 were more abundant and 322 were less abundant. A total of 693 DAPs were identified in the PFS4 vs. PFS2 comparison; 399 were more abundant and 294 were less abundant. We identified the common DAPs between the PFS2 vs. PF2 and PFS4 vs. PF4 comparisons to detect DAPs that were involved in salt stress in both diploid and autotetraploid plants. The list of these common DAPs was then aligned with the list of DAPs from the PF4 vs. PF2 comparison and the DAPs that were represented in the two lists were identified as DAPs that were involved in salt stress only in the autotetraploid plants (S2 Table).

The autotetraploid-specific salt stress DAPs were then analyzed using the GO, COG, and KEGG databases. A total of 152 DAPs were assigned to 26 functional groups under the three main GO categories (Fig 5 and S3 Table).

[Figure omitted. See PDF.]

Fig 5. Gene Ontology classification of differentially abundant proteins that were detected in Paulownia fortunei leaves.

https://doi.org/10.1371/journal.pone.0181937.g005

Under biological process, metabolic process and cellular process were the most represented groups; under cellular component, cell and cell part were the two largest groups; and under molecular function, catalytic activity and binding were the most represented. Interestingly, the GO classifications for the DAPs were similar to the classifications for all the identified proteins. The most enriched GO terms under biological process, cellular component, and molecular function were the same for all the proteins and for the autotetraploid-specific salt stress DAPs subset.

To predict and classify the possible functions of the DAPs, we assigned them to COG categories. Based on their sequence homology with known proteins in GenBank, 28 DAPs (1.06% of all the proteins) were divided into 11 categories (Table 2). The “posttranslational modification, protein turnover, chaperones” category contained five proteins and was the largest, followed by “energy production and conversion” (4), “translation, ribosomal structure and biogenesis” (4), “carbohydrate transport and metabolism” (4), “general function prediction only” (3) (S1 Fig).

[Figure omitted. See PDF.]

Table 2. Statistics of the functional categories of the DAPs that were detected in P. fortunei leaves.

https://doi.org/10.1371/journal.pone.0181937.t002

The DAPs were mapped to 22 KEGG metabolic pathways (S4 Table); among them, “metabolic pathways” were highly enriched, followed by “biosynthesis of secondary metabolites”, “starch and sucrose metabolism”, and “ribosome”.

Confirmation of DAPs by qRT-PCR

To confirm the DAPs identified by the iTRAQ analysis, qRT-PCR was performed to detect the transcript expression levels of the corresponding genes. The qRT-PCR results showed that in the PFS2 vs. PF2 and PFS4 vs. PF4 comparisons, the expression level of the gene corresponding to one of the DAPs was consistent with the iTRAQ LC-MS/MS analysis results, while the expression levels of the genes corresponding to three DAPs showed opposite trends (Fig 6). Interestingly, in the PFS2 vs. PF2 comparison, the expression levels of the genes corresponding to six DAPs displayed the same trends in the qRT-PCR and iTRAQ analysis and only two displayed opposite trends, while in the PFS4 vs. PF4 comparison, the expression levels of the genes corresponding to seven DAPs displayed opposite trends and only one displayed the same trend. These discrepancy may be attributed to posttranscriptional and posttranslational regulatory processes [34].

[Figure omitted. See PDF.]

Fig 6. Quantitative RT-PCR analysis of Paulownia fortunei candidate proteins.

The 18S rRNA of Paulownia was chosen as an internal reference gene for normalization. CL2476.Contig2: strand cell wall invertase; CL1435.Contig2: chaperonin-60 alpha subunit; Unigene38438: UDP-glucuronate decarboxylase 1; Unigene35182: RSI3; CL6922.Contig2: peroxidase; Unigene14040: stromal 70 kDa; Unigene31115: Photosystem I reaction center subunit II; CL10694.Contig2: 2-oxoglutarate dehydrogenase.

https://doi.org/10.1371/journal.pone.0181937.g006

Discussion

Salt stress is expected to become more severe and has been recognized as a major threat to global plants in the future. Although several approaches have been advocated to address this problem, due to the genetic and physiological complexity of the mechanisms involved, limited success has been achieved in attempts to breed salt-tolerant plants [35]. Thus, understanding the molecular mechanisms of salt tolerance of plants has been a focus of much research. An enhanced understanding of the associated genes and molecular mechanisms involved, should make it easier to breed salt-tolerant plants [36]. Altered abundance levels of proteins due to changes in the expressions of many genes occur when plants are subjected to salt stress. These changes may confer susceptibility or protection to salt stress. These changes may confer susceptibility or protection to salt stress. To identify the functions of individual genes that control traits associated with tolerance to salt stress, iTRAQ techniques are used to identify candidate protein for salt stress tolerance in many plants [37, 38]. In this study, to detect variations in the responses of diploid and autotetraploid P. fortunei to salt stress, we employed a proteomic approach on salt-treated and normal diploid and autotetraploid P. fortunei seedlings.

Cell expansion and growth are the first processes to decline when a plant is subjected to salt stress. After that, with more severe salt stress, photosynthesis is adversely affected. Then, a reduction in hydration and increase in ROS can damage cellular membranes and proteins [35]. Plants have developed various mechanisms to overcome such damage, for example, osmotic pressure is often changed during salt stress in plants due to the accumulation of solutes [39] and alterations in the structure of aquaporins or ion channel proteins. Among this process, proline and glycine betaine are the most commonly accumulated compounds in plants in response to environmental stresses [40]. At the same time, soluble sugars can act as moderators of osmotic pressure, in general, salt stress can alter the concentrations of soluble sugars, such as glucose, sucrose, and fructose [41]. However, in severe dehydration conditions, soluble sugars act opposite with their former alteration [42]. Soluble sugars can act not only as moderators of osmotic pressure, but also as signal transductors in salt stress [41] when they may interact with plant hormones to regulate the expression of other genes, such as those involved in photosynthetic metabolism [43]. In our study, sugar transporter as the DAP was up-regulated in the salt stress in the leaf of diploid and autotetraploid P. fortunei. It seems that the accumulation of sugar was closely related to accumulation of sugar transporter genes. In the salt stress wheat, two sugar transporter genes SuT4 and SuT5 were up-regulated, which had greater inhibitory effects on plant growth [36]. In addition, changes of osmotic pressure also as signal transductors may interact with plant hormones to regulate the expression of other genes, such as those involved in photosynthetic process, which further affected plant growth [43, 44]. Therefore, changes of osmotic pressure in leaves may be one of the important regulation pathway in the diploid and autotetraploid P. fortunei under salt stress.

Leaf Na⁺ and the K⁺/Na⁺ discrimination ratio, were proved to be important in salt tolerance in plants[45]. A series of K⁺/Na⁺ channel proteins, Na⁺/H⁺ and Ca²⁺/H⁺ antiporters, and enzymes involved in the transport of ions were found to be altered in autotetraploid and salt-treated Paulownia, indicating that the proton gradient was important in the salt tolerance process [27]. Similarly, a gene encoding vacuolar ATPase subunit B in a salt-sensitive barley cultivar was identified by differential display-reverse transcriptase PCR (DDRT-PCR) [46]. In roots, the ATPase can produce a proton gradient that can drive the uptake of Na⁺ and Ca²⁺ by Na⁺/H⁺ and Ca²⁺/H⁺ antiporters, respectively. Thus, it is supposed to be involved in detaining of Ca²⁺ and Na⁺ ions into the vacuole [47]. Genetic manipulation to add expression of K⁺/Na⁺ channel proteins, Na⁺/H⁺ and Ca²⁺/H⁺ antiporters, and enzymes involved in the transport of ions so as to increase osmoprotectant flux could be a good method to acclimate plants to salt stress [48]. Overall, the mechanisms of the Paulownia response to salt stress can be divided into two main types: minimizing the entry of salt into Paulownia plants, and minimizing the concentration of salt in the cytoplasm. The DAP, ATP synthase delta chain was up-regulated in the salt stress paulownia, play important role in stabilization of thylakoid protein composition, it was reported that knock-out of the genes coding for ATP synthase delta chain in Arabidopsis could affect the integrity of the thylakoid membrane, which further affect electron and proton transport across the thylakoid membrane, resulting in the reduction of the synthesis of NADPH and ATP [49].on the contrary, the up-regulated ATP synthase delta chain in Paulownia may enhance the ATP production and increase the activity of transport of ions to response salt stress.

In a previous study, we examined the differences between drought-treated diploid and autotetraploid P. fortunei using RNA sequencing [50]. The differentially expressed genes that we found in that study do not overlap well with the DAPs detected in this study. There may be several reasons for this discrepancy. Firstly, the plant response and acclimate bioprocesses under salt stress and under drought stress will not be identical. For example, Oztur et al. [51] reported that in barley, the expression of 5% of the genes were changed under salinity stress, while the expression of 15% of the genes were changed under drought stress. Secondly, posttranscriptional regulation and posttranslational modifications influence translational efficiency, which could lead to discordance between mRNA expression levels and protein abundance. Such discordance had been widely reported; for example, by Lan et al. in Arabidopsis roots [52].

Conclusion

Dynamic changes of proteins abundance in the leaf of diploid and autotetraploid P. fortunei were investigated by iTRAQ analysis. A total of 2,634 proteins were identified. After analysis of the differentially abundant proteins, 152 autotetraploid-specific DAPs were obtained. The key metabolic pathways of DAP enriched pathways were detected. they were assigned to 26 GO functional groups and mapped to 22 KEGG metabolic pathways. The up-regulated sugar transporter might promote the accumulation of soluble sugars which might increase changes of osmotic pressure, and the up-regulated ATP synthase delta in ATP production increased the activity of transport of ions which may for proton gradient to response the salt stress. Currently, based on the fact of discordance between mRNA expression levels and protein abundance, we have only a limited understanding the roles of these proteins in the leaf of diploid and autotetraploid P. fortunei under salt stress, further research may focus on the functional genomic of the candidate proteins identified by our proteomic investigation.

Supporting information

[Figure omitted. See PDF.]

S1 Table. KEGG pathway analysis of all identified proteins for Paulownia fortunei.

https://doi.org/10.1371/journal.pone.0181937.s001

(XLSX)

S2 Table. Differentially abundant proteins identified in salt treated Paulownia fortunei.

https://doi.org/10.1371/journal.pone.0181937.s002

(XLSX)

S3 Table. Gene Ontology analysis of differentially abundant proteins in salt treated Paulownia fortunei.

https://doi.org/10.1371/journal.pone.0181937.s003

(XLS)

S4 Table. KEGG analysis of the identified differentially abundant proteins in salt treated Paulownia fortunei.

https://doi.org/10.1371/journal.pone.0181937.s004

(XLSX)

S1 Fig. Classification of the clusters of orthologous groups for the Itraq-based proteotome of Paulownia fortunei.

https://doi.org/10.1371/journal.pone.0181937.s005

(TIF)

Acknowledgments

This work was supported by the Fund of the Transformation Project of the National Agricultural Scientific and Technological Achievement of China (2012GB2D000271), the Central Financial Forestry Science Promotion Project (GTH [2012]01), the Fund of the Technology Innovation Team Project of Zhengzhou (121PCXTD515), the Joint Funds of the National Natural Science Foundation of China (NSFC) (Grant No. U1204309), the Fund of the Research Cooperation Project of Henan Province (Grant No. 152107000097) and the Fund of Zhongyuan Scholarship Foundation of Henan Province (122101110700). we thank Z. W. for his interest and invaluable assistance during this work.

Citation: Deng M, Dong Y, Zhao Z, Li Y, Fan G (2017) Dissecting the proteome dynamics of the salt stress induced changes in the leaf of diploid and autotetraploid Paulownia fortunei. PLoS ONE 12(7): e0181937. https://doi.org/10.1371/journal.pone.0181937

References

1. Gong Q, Li P, Ma S, Rupassara SI, Bohnert HJ. Salinity stress adaptation competence in the extremophile Thellungiella halophila in comparison with its relative Arabidopsis thaliana. Plant J. 2005;44(5):826–39. pmid:16297073

2. Munns R, Tester M. Mechanisms of Salinity Tolerance. Annu Rev Plant Biol. 2008;59(1):651–81.

3. Sunkar R, Bartels D, Kirch HH. Overexpression of a stress-inducible aldehyde dehydrogenase gene from Arabidopsis thaliana in transgenic plants improves stress tolerance. Plant J. 2003;35(4):452–64. pmid:12904208

4. Zhou L, Wang NN, Gong SY, Lu R, Li Y, Li XB. Overexpression of a cotton (Gossypium hirsutum) WRKY gene, GhWRKY34, in Arabidopsis enhances salt-tolerance of the transgenic plants. Plant Physiol Biochem. 2015;96:311–20. pmid:26332661

5. Raghavendra AS, Gonugunta VK, Christmann A, Grill E. ABA perception and signalling. Trends Plant Sci. 2010;15(7):395–401. pmid:20493758

6. Fu M, Kang HK, Son SH, Kim SK, Nam KH. A subset of Arabidopsis RAV transcription factors modulates drought and salt stress responses independent of ABA. Plant Cell Physiol. 2014;55(11):1892–904. pmid:25189341

7. Morris PC. Integrating lipid signalling, mitogen-activated protein kinase cascades and salt tolerance. New Phytol. 2010;188(3):640–3. pmid:20976867

8. Vaidyanathan H, Sivakumar P, Chakrabarty R, Thomas G. Scavenging of reactive oxygen species in NaCl-stressed rice (Oryza sativa L.)—differential response in salt-tolerant and sensitive varieties. Plant Sci. 2003;165(6):1411–8.

9. Xiang J, Ran J, Zou J, Zhou X, Liu A, Zhang X, et al. Heat shock factor OsHsfB2b negatively regulates drought and salt tolerance in rice. Plant Cell Rep. 2013;32(11):1795–806. pmid:23949687

10. Fan G, Li X, Deng M, Zhao Z, Yang L. Comparative analysis and identification of miRNAs and their target genes responsive to salt stress in diploid and tetraploid Paulownia fortunei seedlings. PloS one. 2016;11(2):e0149617. pmid:26894691; PubMed Central PMCID: PMC4764520.

11. Fan H, Xu Y, Du C, Xue W. Phloem sap proteome studied by iTRAQ provides integrated insight into salinity response mechanisms in cucumber plants. J Proteomics. 2015;125:54–67. pmid:25958826

12. Ngom M, Gray K, Diagne N, Oshone R, Fardoux J, Gherbi H, et al. Symbiotic performance of diverse Frankia strains on salt-stressed Casuarina glauca and Casuarina equisetifolia plants. Front Plant Sci. 2016;7:1331. pmid:27630656; PubMed Central PMCID: PMCPMC5006599.

13. Feng K, Yu J, Cheng Y, Ruan M, Wang R, Ye Q, et al. The SOD gene family in tomato: identification, phylogenetic relationships, and expression patterns. 2016;7:1279. pmid:27625661

14. Jin K, Chen S, Zhu N, Yu F, Soltis PS, Soltis DE. Comparative proteomics of the recently and recurrently formed natural allopolyploid Tragopogon mirus (Asteraceae) and its parents. New Phytol. 2012;196(1):292–305. pmid:22861377

15. Alam I, Sharmin SA, Kim KH, Yang JK, Choi MS, Lee BH. Proteome analysis of soybean roots subjected to short-term drought stress. Plant Soil. 2010;333(1):491–505.

16. Riccardi F, Gazeau P, de Vienne D, Zivy M. Protein changes in response to progressive water deficit in maize. Quantitative variation and polypeptide identification. Plant Physiol. 1998;117(4):1253–63. pmid:9701581; PubMed Central PMCID: PMC34889.

17. Yoshimura K, Masuda A, Kuwano M, Yokota A, Akashi K. Programmed proteome response for drought avoidance/tolerance in the root of a C(3) xerophyte (wild watermelon) under water deficits. Plant Cell Physiol. 2008;49(2):226–41. pmid:18178965.

18. Salekdeh GH, Siopongco J, Wade LJ, Ghareyazie B, Bennett J. Proteomic analysis of rice leaves during drought stress and recovery. Proteomics. 2002;2(9):1131–45. pmid:12362332.

19. Hajheidari M, Abdollahian-Noghabi M, Askari H, Heidari M, Sadeghian SY, Ober ES, et al. Proteome analysis of sugar beet leaves under drought stress. Proteomics. 2005;5(4):950–60. pmid:15712235.

20. Parent B, Suard B, Serraj R, Tardieu F. Rice leaf growth and water potential are resilient to evaporative demand and soil water deficit once the effects of root system are neutralized. Plant Cell Environ. 2010;33(8):1256–67. pmid:20302604.

21. Abraham E, Hourton-Cabassa C, Erdei L, Szabados L. Methods for determination of proline in plants. Methods Mol Biol. 2010;639:317–31. pmid:20387056.

22. Hoque MA, Okuma E, Banu MN, Nakamura Y, Shimoishi Y, Murata Y. Exogenous proline mitigates the detrimental effects of salt stress more than exogenous betaine by increasing antioxidant enzyme activities. J Plant Physiol. 2007;164(5):553–61. pmid:16650912.

23. Khedr AH, Abbas MA, Wahid AA, Quick WP, Abogadallah GM. Proline induces the expression of salt-stress-responsive proteins and may improve the adaptation of Pancratium maritimum L. to salt-stress. J Exp Bot. 2003;54(392):2553–62. pmid:14512386.

24. Chaires M, Gupta D, Joshee N, Cooper KK, Basu C. RNA-seq analysis of the salt stress-induced transcripts in fast-growing bioenergy tree, Paulownia elongata. J Plant Interact, 2017;12(1):128–36.

25. Zhang X, Deng M, Fan G. Differential Transcriptome analysis between Paulownia fortunei and its synthesized autopolyploid. Int J Mol Sci. 2014;15(3):5079–93. pmid:24663058

26. Li Y, Fan G, Dong Y, Zhao Z, Deng M, Cao X, et al. Identification of genes related to the phenotypic variations of a synthesized Paulownia (Paulownia tomentosa x Paulownia fortunei) autotetraploid. Gene. 2014;553(2):75–83. pmid:25300252.

27. Fan G, Wang L, Deng M, Zhao Z, Dong Y, Zhang X, et al. Changes in Transcript Related to Osmosis and Intracellular Ion Homeostasis in Paulownia tomentosa under Salt Stress. Front Plant Sci. 2016;7:384. pmid:27066034

28. Tang Z, Du W, Du XL, Ban YY, Cheng JL. iTRAQ Protein Profiling of Adventitious Root Formation in Mulberry Hardwood Cuttings. J Plant Growth Regul. 2016;35(3):618–31.

29. Meng Q, Hou L, Zhao Y, Huang X, Huang Y, Xia S, et al. iTRAQ-based proteomic study of the effects of Spiroplasma eriocheiris on Chinese mitten crab Eriocheir sinensis hemocytes. Fish Shellfish immunol. 2014;40(1):182–9. pmid:25017370.

30. Qiao J, Wang J, Chen L, Tian X, Huang S, Ren X, et al. Quantitative iTRAQ LC-MS/MS proteomics reveals metabolic responses to biofuel ethanol in Cyanobacterial Synechocystis sp. PCC 6803. J Proteome Res. 2012;11(11):5286–300. pmid:23062023.

31. Lin XQ, Liang SL, Han SY, Zheng SP, Ye YR, Lin Y. Quantitative iTRAQ LC-MS/MS proteomics reveals the cellular response to heterologous protein overexpression and the regulation of HAC1 in Pichia pastoris. J Proteomics. 2013;91:58–72. pmid:23851310.

32. Guo H, Chen C, Lee DJ, Wang A, Gao D, Ren N. Coupled carbon, sulfur and nitrogen cycles of mixotrophic growth of Pseudomonas sp. C27 under denitrifying sulfide removal conditions. Bioresource Technology. 2014;171:120–6. pmid:25189517

33. Fan G, Dong Y, Deng M, Zhao Z, Niu S, Xu E. Plant-pathogen interaction, circadian rhythm, and hormone-related gene expression provide indicators of phytoplasma infection in Paulownia fortunei. Int J Mol Sci. 2014;15(12):23141–62. pmid:25514414; PubMed Central PMCID: PMC4284757.

34. 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. J Proteome Res. 2015;14(1):268. pmid:25367710

35. Hu L, Zhang P, Jiang Y, Fu J. Metabolomic analysis revealed differential adaptation to salinity and alkalinity stress in Kentucky Bluegrass (Poa pratensis). Plant Mol Biol Rep. 2014;33(1):56–68.

36. Charkazi F, Ramezanpour SS, Soltanloo H. Expression pattern of two sugar transporter genes (SuT4 and SuT5) under salt stress in wheat. Plant Omics. 2010;3(6):194–8.

37. Rabey HAE, Al-Malki AL, Abulnaja KO. Proteome analysis of Date Palm (Phoenix dactylifera L.) under severe drought and salt stress. Int J Genomics. 2016;2016:7840759. pmid:27840818

38. Jiang Q, Li X, Niu F, Sun X, Hu Z, Zhang H. iTRAQ-based quantitative proteomic analysis of wheat roots in response to salt stress. Proteomics. 2017;17(8):1600265. pmid:28191739

39. Orcutt DM, Nilsen ET. 2000. The physiology of plants under stress: soil and biotic factors. Wiley, New York, pp 30

40. Hasegawa PM, Bressan RA, Zhu JK, Bohnert HJ. Plant cellular and molecular responses to high salinity. Annu Rev Plant Physiol Plant Mol Biol. 2000;51:463–99. pmid:15012199

41. Chaves MM, Oliveira MM. Mechanisms underlying plant resilience to water deficits: prospects for water-saving agriculture. J Exp Bot. 2004;55(407):2365–84. pmid:15475377.

42. Pinheiro C, Chaves MM, Ricardo CP. Alterations in carbon and nitrogen metabolism induced by water deficit in the stems and leaves of Lupinus albus L. J Exp Bot. 2001;52(358):1063–70. pmid:11432922.

43. Chaves MM, Flexas J, Pinheiro C. Photosynthesis under drought and salt stress: regulation mechanisms from whole plant to cell. Ann Bot. 2009;103(4):551–60. pmid:18662937; PubMed Central PMCID: PMC2707345.

44. Chaves MM, Oliveira MM. Mechanisms underlying plant resilience to water deficits: prospects for water-saving agriculture. J Exp Bot. 2004;55(407):2365–84. pmid:15475377

45. Saleh B. Effect of salt stress (NaCl) on biomass and K+/Na+ ratio in cotton. J Stress Physiol Biochem. 2011;7(4):05–14.

46. Wei W, Bilsborrow P, Hooley P, Fincham D, Forster B. Variation between two near isogenic barley (Hordeum vulgare) cultivars in expression of the B subunit of the vacuolar ATPase in response to salinity. Hereditas. 2000;135(2–3):227–31.

47. Taiz L. The plant vacuole. J Exp Biol. 1992;172(1):113–22.

48. Huang J, Hirji R, Adam L, Rozwadowski KL, Hammerlindl JK, Keller WA, et al. Genetic engineering of glycinebetaine production toward enhancing stress tolerance in plants: metabolic limitations. Plant Physiol. 2000;122(3):747–56. pmid:10712538; PubMed Central PMCID: PMC58910.

49. Maiwald D, Dietzmann A, Jahns P, Pesaresi P, Joliot P, Joliot A, et al. Knock-out of the genes coding for the Rieske protein and the ATP-synthase δ-subunit of Arabidopsis. Effects on photosynthesis, thylakoid protein composition, and nuclear chloroplast gene expression. Plant Physiol. 2003;133(1):191–202. pmid:12970486

50. Dong Y, Fan G, Zhao Z, Deng M. Compatible solute, transporter protein, transcription factor, and hormone-related gene expression provides an indicator of drought stress in Paulownia fortunei. Funct Integr Genomics. 2014;14(3):479–91. pmid:24801596; PubMed Central PMCID: PMC4137158.

51. Oztur ZN, Talamé V, Deyholos M, Michalowski CB, Galbraith DW, Gozukirmizi N, et al. Monitoring large-scale changes in transcript abundance in drought- and salt-stressed barley. Plant Mol Biol. 2002;48(5):551–73.

52. Lan P, Li W, Schmidt W. Complementary proteome and transcriptome profiling in phosphate-deficient Arabidopsis roots reveals multiple levels of gene regulation. Mol Cell Proteomics. 2012;11(11):1156–66. pmid:22843991; PubMed Central PMCID: PMC3494196.