Next Article in Journal
Amino Acid Transport Defects in Human Inherited Metabolic Disorders
Previous Article in Journal
Characterization of the Inhibitory Effect of Gastrodigenin and Gastrodin on M-type K+ Currents in Pituitary Cells and Hippocampal Neurons
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Proteomic Analysis Reveals the Regulatory Effects of H2S on Salt Tolerance of Mangrove Plant Kandelia obovata

Key Laboratory of the Ministry of Education for Coastal and Wetland Ecosystems, College of the Environment and Ecology, Xiamen University, Xiamen 361102, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(1), 118; https://doi.org/10.3390/ijms21010118
Submission received: 22 November 2019 / Revised: 17 December 2019 / Accepted: 19 December 2019 / Published: 23 December 2019
(This article belongs to the Section Molecular Biology)

Abstract

:
As a dominant mangrove species, Kandelia obovata is distributed in an intertidal marsh with an active H2S release. Whether H2S participates in the salt tolerance of mangrove plants is still ambiguous, although increasing evidence has demonstrated that H2S functions in plant responses to multiple abiotic stresses. In this study, NaHS was used as an H2S donor to investigate the regulatory mechanism of H2S on the salt tolerance of K. obovata seedlings by using a combined physiological and proteomic analysis. The results showed that the reduction in photosynthesis (Pn) caused by 400 mM of NaCl was recovered by the addition of NaHS (200 μM). Furthermore, the application of H2S enhanced the quantum efficiency of photosystem II (PSII) and the membrane lipid stability, implying that H2S is beneficial to the survival of K. obovata seedlings under high salinity. We further identified 37 differentially expressed proteins by proteomic approaches under salinity and NaHS treatments. Among them, the proteins that are related to photosynthesis, primary metabolism, stress response and hormone biosynthesis were primarily enriched. The physiological and proteomic results highlighted that exogenous H2S up-regulated photosynthesis and energy metabolism to help K. obovata to cope with high salinity. Specifically, H2S increased photosynthetic electron transfer, chlorophyll biosynthesis and carbon fixation in K. obovata leaves under salt stress. Furthermore, the abundances of other proteins related to the metabolic pathway, such as antioxidation (ascorbic acid peroxidase (APX), copper/zinc superoxide dismutase (CSD2), and pancreatic and duodenal homeobox 1 (PDX1)), protein synthesis (heat-shock protein (HSP), chaperonin family protein (Cpn) 20), nitrogen metabolism (glutamine synthetase 1 and 2 (GS2), GS1:1), glycolysis (phosphoglycerate kinase (PGK) and triosephosphate isomerase (TPI)), and the ascorbate–glutathione (AsA–GSH) cycle were increased by H2S under high salinity. These findings provide new insights into the roles of H2S in the adaptations of the K. obovata mangrove plant to high salinity environments.

Graphical Abstract

1. Introduction

Mangroves comprise a complex and unique ecosystem that distributes along tropical and subtropical coastal tidal zone [1]. Mangroves are facultative halophytes and potential stress adaptors due to their special morphological, anatomical, physiological, and biochemical features [1]. Different mangrove species have different salinity preferences and achieve optimal growth at different salinity levels. Jayatissa et al. reported that there is an optimal growth for Sonneratia caseolaris with the lowest salt tolerance at 3–5 ppt salinity, and there is an optimal growth for Avicennia marina with the highest salt tolerance at 25–27 ppt salinity [2]. As a dominant species of mangrove plants, Kandelia obovata can distribute in areas with salinities up to 27.58 ppt seawater level [3]. The salinity between 5 and 15 ppt is suitable for indoor-cultured K. obovata growth, while salinity up to 20 ppt inhibits the photosynthesis and growth of K. obovata [4].
Tidal inundation, which originates from coastal waters, influences sediment salinity and pore-water sulfide concentration [5]. Soil salinities and porewater sulfide concentrations are further enhanced with high inundation [6]. Hydrogen sulfide (H2S) is the dominant sulfide-containing gas that is emitted from intertidal sediment. A mangrove forest was found to serve as a perpetual source of H2S, and the annual mean emission of H2S found in this study was 768 ± 240 μg·S·m−2·d−1 [7]. Mangrove plants affect the sulfide concentration and H2S emission flux of sediments [8], and this research focuses on whether H2S influences the adaptation of mangrove plants to harsh environmental conditions such as high salinity.
As seen in huge body of previous studies, plants gain their tolerance to abiotic stress by actively synthesizing H2S. Hou et al. reported that drought stress leads to stomatal closure through an increase in the expression and activity of d/l-cysteine desulfhydrase (d/l-CD), a key enzyme of H2S biosynthesis in the leaves of Vicia faba [9]. The H2S pretreatment was found to significantly increase the regrowth ability of tobacco suspension-cultured cells after heat stress by alleviating a reduction in cell viability [10]. Salinity can detrimentally impact plant biomass production, which is associated with a salinity-induced a reduction in the amount of chlorophyll [11]. Previous studies have mostly examined the role of H2S in the mitigation of oxidative stress in salt-stressed plants [12,13,14]. The activities of superoxide dismutase (SOD), catalase (CAT), ascorbic acid peroxidase (APX), glutathione reductase (GR), glutathione peroxidase (GPX) and dehydroascorbate reductase (DHAR) were found to be increase in salt-stressed cucumber seedlings with the addition of NaHS, which is a widely used H2S donor, while, under the same experimental conditions, lipid peroxides and hydrogen peroxide (H2O2) levels were found to decrease [12]. The improved germination rate of salt-stressed cucumber seeds by H2S might have been due to the breakdown of starch by H2S-mediated α- and β-amylase induction in endosperm, which ultimately led to improved hypocotyl growth [13]. To salt-stressed strawberry, root pre-treatment with H2S has been found to induce an increase in stomatal conductance, photosynthesis, and leaf water content [14]. In salt-stressed rice, NaHS has been found to increase chlorophyll and protein content [15]. However, in a study by Koch et al., H2S negatively affected the anoxic production of energy and the energy-dependent N uptake in roots of the freshwater marsh species Spartina alterniflora [16].
In most true mangroves, the adaptations of seedlings to tolerate saline conditions have been widely studied [17]. However, the role of H2S in salinity tolerance improvement in mangroves is still mostly unknown. In this study, an appropriate concentration of NaHS, a donor of H2S, had a positive effect on K. obovata under high salinity considering both physiological and proteomic aspects. The K. obovata seedlings were exposed to 400 mM of NaCl with or without 200 μM of NaHS for seven days, and proteomics-based methodology (Two-dimensional (2-DE) accompanied by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI TOF MS)) was adopted to investigate the specific metabolic pathways and regulatory mechanisms of H2S on the salt tolerance of K. obovata seedlings.

2. Results

2.1. Effects of H2S and NaCl on Growth and Photosynthesis of K. Obovata Seedlings

As expected, exogenous NaHS treatments (varying from 50 to 350 μM) progressively alleviated the NaCl-induced decrease in inhibition on the leaf photosynthesis of K. obovata seedlings. On the other side, a high concentration of NaHS (500 μM) exhibited no beneficial effects (Figure 1). A powerful relationship was found among the observed parameters through two-tailed Pearson correlation and path analysis (Supplementary Table S2). Most of the physiological indices were positively correlated with the net photosynthetic rate (Pn) except intercellular carbon dioxide (Ci). The Pn had negative but significant associations (−0.655*) with Ci and displayed significant interactions with leaf dry weight (DW), chlorophyll content (Chl), stomatal conductance (Gs) and transpiration rate (Tr). Among different concentrations of NaHS, 200 μM was the most effective in relieving NaCl-induced adverse effects in K. obovata seedlings. The recovery rate of Pn at 200 μM was the highest compared with other concentrations (Supplementary Figure S1 and Supplementary Table S1). We used 200 µM of NaHS in the next experiments.

2.2. H2S Rather than other Derivatives from Nahs Alleviates NaCl-Induced Reduction in Chlorophyll Content

It should be noticed that the NaHS solution contained not only H2S/HS but also other sulfur-containing components. The research results showed that the above chemicals (Na2SO4, NaHSO3, NaHSO4, and CH3COONa), failed to rescue the NaCl-induced reduction of chlorophyll content (Supplementary Figure S2A). Meanwhile, compared with the experimental control group (CK), the endogenous H2S content kept a high level with the NaHS treatment but maintained a stable low level under other sulfur-containing compounds or sodium-containing compounds (Supplementary Figure S2B). HS and/or H2S has an essential role in relieving the reduction in chlorophyll content in K. obovata seedlings under salt treatment but not other sulfur-containing compounds or sodium-containing compounds (Na2SO4, NaHSO3, NaHSO4, and CH3COONa).

2.3. Effects of H2S and NaCl on the Characteristics of Chlorophyll Fluorescence

The variable fluorescent/maximal fluorescent (Fv/Fm) ratio reflects the potential quantum efficiency of photosystem II (PSII), widely used as a sensitive indicator of stress-induced damage to the photosynthetic apparatus [18]. After seven days of treatment with NaCl, K. obovata seedlings’ leaves were used to measure the chlorophyll fluorescence. The results showed that the NaCl treatment significantly decreased the ratio of Fv/Fm by 27.18% compared to the control. However, when exogenous H2S was applied to the NaCl-treated seedlings, the Fv/Fm was increased by 16.8% (Supplementary Figure S3).

2.4. Effects of H2S and NaCl on Oxidative Stress and Activity of the Antioxidant System

Cytomembranes are damaged by reactive oxygen species (ROS) that were induced by salinity. The content of H2O2 was analyzed to evaluate the effect of H2S on the oxidative burst of NaCl-treated K. obovata seedlings. When exogenous H2S was applied to NaCl-treated leaves, the content of H2O2 was decreased, but it did not significantly differ from that of the NaCl treatment (Figure 2A). Leaf electrolyte leakage percentage (REL) and the thiobarbituric acid reactive substances (TBARS) accumulation were measured to estimate the protective role of H2S on the stability of the membrane. The TBARS levels were decreased by the exogenous H2S addition, but these levels did not significantly differ from that of the NaCl treatment (Figure 2B). Meanwhile, high salinity significantly increased the REL of K. obovata leaves, being 2.0-fold higher than that of the control. The exogenous H2S significantly alleviated the NaCl-induced increase in REL (Figure 2C).
Additionally, the content of glutathione (GSH) significantly increased by around 4-fold after the NaCl and H2S treatment compared to the NaCl treatment (Figure 2D). Plant cells have evolved the antioxidant enzymatic system, such as SOD and APX, to cope with stress-induced ROS generation. When exogenous H2S was applied to NaCl-treated seedlings, the activity of APX and SOD were increased but did not significantly differ from that of the NaCl treatment (Figure 2E,F).

2.5. Identification and Classification of DEPs

The 2-DE proteomic approach was employed to analyze the protein expression profile to investigate the mechanism of the H2S-mediated salt tolerance of K. obovata. Three biological replicates were performed with similar results. According to the image analysis by the PDQuest software, approximately 450 protein spots were reproducibly resolved on each gel. The protein spots were detected in the isoelectric point (pI) range from 4 to 7, and the range of molecular weight (MW) was 14.4–66.2 kDa. A total of 37 spots with significant changes (higher than 2.0-fold, p < 0.05) were identified as differentially expressed proteins (DEPs) and presented in representative image (Figure 3A). Close-up images from four protein spots were displayed in Figure 3B.
We provide the information for 37 DEPs in Table 1 and Supplementary Table S2. The difference between the two groups (NaCl vs. CK and NaCl and H2S vs. CK) was statistically significant. A comprehensive overview of the expression patterns of 37 identified proteins is shown in Figure 4A. They were divided into five functional categories, including carbohydrate and energy metabolism (29.73%), photosynthesis (18.92%), stress response proteins and cell structure (18.92%), amino acid and protein metabolism (16.22%), and hormone biosynthesis and transcription factor (10.81%) by their biological functions, each of which is in the UniProt database (Figure 4B). There were fifteen proteins that significantly changed under the NaCl treatment; specifically, the abundances of seven proteins were increased, and the abundances of eight proteins were decreased. Under the NaCl and H2S treatment, fourteen proteins were significantly changed, the abundances of nine proteins were increased, and five proteins were decreased. There were six proteins that were affected by both the NaCl and NaCl and H2S treatments (Figure 4A). Compared with the NaCl vs. CK group, there were more up-regulated DEPs in the NaCl and H2S vs. CK group. These up-regulated DEPs were mostly divided into carbohydrate metabolism, energy metabolism, and stress response proteins (Figure 4B). There were less down-regulated DEPs divided into carbohydrate metabolism, energy metabolism, and photosynthesis under the NaCl and H2S treatment (Figure 4B).

2.6. STRING Interaction Networks among DEPs

To predict protein–protein interaction (PPI) networks with the Search Tool of the Retrieval of Interaction Genes/Proteins (STRING) database, an alignment of the amino acid sequence to find the corresponding homologous proteins from Arabidopsis thaliana was carried out. The PPI network visualized by Cytoscape is shown in Supplementary Figure S4. The PPI network of DEPs was found to be composed of 27 nodes and 129 edges. A total of two modules were identified in the PPI network, and these modules include nine and eight proteins (Figure 5, Supplementary Table S3). Module A has nine nodes and 32 interactions. The key node in module A is glutamine synthetase 2 (GS2). Module B has eight nodes and 22 interactions. The key node in module B is ribulose-1,5-bisphosphate carboxylase large subunit (RBCL).
An acyclic graph was constructed and directed by using the Biological Network Gene Ontology tool (BiNGO) to depict the visual interactions of functions based on the enrichment levels of the GO terms. The present study identified eight significantly enriched GO terms in biological processes, and especially in terms of response to abiotic stimulus, response to stress, and the generation of precursor metabolites and energy (Supplementary Figure S5).

2.7. Comparisons of Expression Patterns between Protein and Transcript of Selected DEPs

To further verify the results of the proteome, we randomly selected six DEPs to further analyze their transcriptional expressions under different treatments by using quantitative real-time PCR. In particular, the consistency between the protein and mRNA levels of the following proteins was compared: transcription factor basic helix-loop-helix 145 (bHLH145, spot 13), superoxide dismutase 4 (SOD, spot 18), heat-shock protein (HSP, spot 19), oxygen-evolving enhancer protein 1 (OEE1, spot 14), abscisic acid stress ripening protein (Asr, spot 31) and glutamine synthetase 1;1 (GS1;1, spot 33). The detailed list of the corresponding primer pairs and the detailed PCR procedures used for the real-time PCR are summarized in Supplemental Table S4. As shown in Figure 6, four selected genes showed a similar pattern with their corresponding protein, except for bHLH 145 and GS1;1.

3. Discussion

3.1. H2S but Not Other NaHS Derivatives Contributes to Chlorophyll Content Recovery Decreased by NaCl Treatment

NaHS, which in solution dissociates to Na+ and HS (of which the latter associates with H+ to produce H2S), is widely used to examine the biological effects of H2S [19]. In contrast to the results from NaHS, treatment with other compounds failed to induce endogenous H2S accumulation and to alleviate the reduction in chlorophyll content caused by NaCl exposure (Supplementary Figure S2). Accordingly, this study provides evidence that NaHS contributes to chlorophyll content recovery by HS, rather than Na+. NaHS-associated responses are H2S-specific in K. obovata.

3.2. Exogenous H2S Alleviates Growth and Photosynthesis Inhibition Induced by NaCl in K. obovata Seedlings

Parameters such as relatively low fresh weight, relatively low dry weight, and reduced photosynthesis are potential indicators of salinity stress in plants [11]. In this study, significant variations in the studied physiological traits of K. obovata under different concentrations of NaHS with 400 mM of NaCl were recorded. When the concentration of H2S reached 500 μM, the leaf dry weight and photosynthesis decreased (Figure 1). These findings are in agreement with those of Chen et al. [20], who reported a similar decline in the biomass and photosynthesis of Spinacia oleracea after treatment with a high concentration of NaHS. The change in Ci was analyzed as an indicator of the dominant factor of stomatal or non-stomatal limitations to photosynthesis. Ci showed a reverse trend to Gs, which suggests that non-stomatal closure is the dominant limitation to photosynthesis [21]. In the present research, the decrease in Ci showed a reverse trend to Gs under various concentrations of NaHS, which indicated that non-stomatal limitations dominated.
The previous study indicated that H2S could indirectly enhance the photosynthesis reaching its maximal values, as well as the maximal photochemical efficiency of photosystem II (Fv/Fm) [20]. In the present research, the NaCl and H2S treatment also increased Fv/Fm (Supplementary Figure S3). Additionally, the proteomics analyses revealed six proteins involved in the light reactions of photosynthesis, including electron transfer, light-harvesting, and the light-induced oxidation of water. The level of ferredoxin-NADP reductase (FNR, spot 2), which catalyzes the removal of the electron from photosystem I to NADP+, and otherwise nicotinamide adenine dinucleotide phosphate (NADPH) is oxidized to produce reducing power via ferredoxin (Fd) [22], which is recovered under the NaCl and H2S treatment. The O2 evolving complex 33kD family protein (OEC spot 12) has been termed as a chloroplast manganese-stabilizing protein due to its functions in assembly, stabilization and protection of the manganese cluster required for oxygen evolution in higher plants [23]. Though the accumulation of OEC was down-regulated in two groups, the accumulation of OEC was up-regulated relative to the NaCl treatment, which implied that H2S maintained the photosynthetic electron transport activity during K. obovata exposure to high salt stress. Additionally, NaCl-induced the down-regulation of photosynthetic electron transfer C (PETC, spot 10), which is an important subunit of the Cytb6f complex that participates in linear electron transport [24], was alleviated by the addition of H2S. Nevertheless, plastocyanin (PC, spot 1) was up-regulated under NaCl, which is consistent with the study on K. candel under high salinity did by Wang et al. [25]. On another hand, the abundances of PC remained unchanged under the combined NaCl and NaHS treatment. The results of proteomics analyses are in agreement with physiological measurements that showed that the treatment with NaHS alleviated the NaCl-caused inhibition of photosynthesis.
Photosynthetic carbon assimilation is driven by ribulose-1,5-bisphosphate (RuBP) carboxylase (Rubisco). Ribulose-1,5-bisphosphate carboxylase large subunit (RBCL, spot 8) was down-regulated by the NaCl treatment and unchanged under the NaCl and NaHS treatment. It is noteworthy that ribulose activase small isoform precursor (RBCS, spot 36) showed a reverse trend that was unchanged under the NaCl treatment and down-regulated by the NaCl and NaHS treatment. The capacity of photosynthesis is associated with the level of the RBCL/RBCS ratio and the higher photosynthetic rate exhibited the higher RBCL/RBCS ratio in rice leaves [26]. The up-regulation of the RBCL/RBCS ratio by the NaHS addition under salt treatment in our study indicated that the alleviation of Rubisco inhibition inevitably leads to the improvement of photosynthesis.
Phosphoglycolate phosphatase (PGP, spot 25) is involved in photorespiration, which is required to give carbon from 2-phosphoglycolate back into metabolism [27]. The allocation of photosynthetic electron transport to photorespiration has been found to be enhanced with the decrease of the chlorophyll content [27]. In this study, the NaCl treatment strongly up-regulated PGP, which indicated that the allocation of photosynthetic electron transport to the Calvin cycle was reduced [28]. The accumulation of PGP was not significantly different from the control treatment by the NaHS addition.
The photosynthetic electron transport chain converts light into chemical energy, supplying ATP and NADPH to drive the carbon dioxide reduction and fixation processes [29]. ATP synthase F1 complex (CF1) epsilon subunit (atpC, spot 7) is involved in the energy supply needed for the carbon dioxide reduction and fixation processes. Though the accumulation of atpC was down-regulated in our study, the accumulation of atpC was modified relative to the NaCl treatment, thus implying the NaHS addition might have had a positive effect to supply the energy.

3.3. H2S Rescues the Primary Metabolism Altered by NaCl

In this study, triosephosphate isomerase (TPI, spots 16), and phosphoglycerate kinase (PGK, spot 35) were found to be down-regulated or unchanged under high salinity, respectively. PGK plays a crucial role in catalyzing the ATP-dependent phosphorylation of phosphoglycerate in the Calvin cycle [30]. The two glycolytic enzymes were all positivity modulated by the NaHS addition, suggesting much more energy was produced by glycolysis.
All levels of plant function are affected by nitrogen metabolism [31]. In higher plants, glutamine synthetase (GS2 and GS1:1, spots 32 and 33, respectively) is a crucial enzyme of primary ammonium assimilation and nitrogen metabolism [32]. Glutamine synthetase participates in the synthesis of GSH through the glutamate biosynthesis pathway; thus, the overexpression of glutamine synthetase leads to more GSH formation [33,34]. It is interesting to note that the accumulation of glutamine synthetase was positivity modulated by the NaHS addition in our study.
HSP (spot 19) and 20 kDa chaperonin family protein (Cpn 20, spot 26) were involved in protein synthesis, a process that was positively modulated by the NaHS addition. As molecular chaperones, HSP and chaperonin family protein participate in protein transport, protein folding, and protein assembly processes [35]. During abiotic stress, chaperonins help plants combat increasing amount of incorrectly folded proteins [35].
Previous publications have reported that the TCA cycle provides essential precursors for amino acid biosynthesis and energy metabolism. ATP synthase subunit d (ATP5PD, spot 9) and ATP synthase beta subunit (ATP5F1B, spot 29) are located in the mitochondria. The ATP synthase is an important enzyme that catalyzes energy production by synthesizing ATP from ADP [36], providing energy to the TCA cycle. In mammals, increased oxidative stress represents an important factor for isolated disorders of ATP synthase [37], and there is a connection between ATP synthase and oxidative stress. Moreover, it has recently been predicted that a protein network in which ATP synthase is involved can contribute to protect plants against photo-oxidative damage [38]. In our study, the NaCl and H2S treatment significantly up-regulated the accumulation of ATP5PD and ATP5F1B. In addition, nucleoside diphosphate kinase 1 (NDPK1, spot 30), a housekeeping enzyme that maintains levels of CTP, GTP, and UTP in cells, was identified to be down-regulated under high salinity and significantly up-regulated by the NaHS addition. A recent study showed that NDPK grants transgenic plants tolerance to multiple stresses, including salt and extreme temperatures [39].
Remarkably, H2S rescued the primary metabolism altered by NaCl in our study. In a previous study, Li et al. also reported that increment in energy metabolism by H2S might be responsible for abiotic stress alleviation [40].

3.4. H2S Relieves Oxidative Stress Induced by NaCl

The salt-induced disruption of normal plant metabolism results in the accumulation of harmful ROS [41]. Excess ROS not only cause growth inhibition but also programmed cell death (PCD). The accumulation of H2O2 and consequent increase in malondialdehyde and REL are indicators of salt-induced oxidative damage to membranes [42]. In the present study, the NaHS pretreatment reduced the accumulation of REL under salt stress, which is in agreement with a previous study by Shi et al. [43]. In addition, Actin7 isoform 1 (Actin7, spot 23)—which is thought to aid in the stability of organelle membranes in plant cells and to induce PCD with the stabilization of actin depolymerization [44,45]—was significantly up-regulated under the NaCl and H2S treatment in our study. Meanwhile, plants trigger various enzymatic and non-enzymatic antioxidants to detoxify ROS and prevent cellular damage [46]. In this study, compared with the NaCl treatment, the content of GSH was found to be induced by application of H2S (Figure 2). Ascorbate (AsA) and GSH play important roles in the AsA–GSH cycle [47]. As an electron donor, APX reduces H2O2 into water by using AsA, and monodehydroascorbate can react with GSH to produce AsA and GSSG catalyzed by DHAR [47]. In our study, the increase in GSH content and the activity of APX indicate that H2S promoted the AsA–GSH cycle that was responsible for the removal of H2O2 under high salinity. In agreement, APX (spot 37) protein expression was found to be up-regulated under the NaCl and H2S treatment. Moreover, copper/zinc superoxide dismutase CSD2A-1 (CSD2, spot 17) was found to be up-regulated under the NaCl and H2S treatment, but SOD (spot 18) showed the opposite trend. The SODs were divided into four types according to the different metal cofactors in the catalytic site, and CSD2 is mainly found in higher plants [48]. We propose that CSD2 plays a critical role in a direct and rapid mechanism of ROS detoxification under high salinity. The SODs (spot 18) might be iron SOD, which was found in the cytoplasmic of plant or the bacterial cells and played no major role in K. obovata leaves [48]. In short, based on the above results, we can infer that H2S increased non-enzymatic antioxidants to detoxify ROS and to prevent cellular damage under high salinity.
Moreover, some abiotic stress-related proteins, such as cysteine proteinase (CP, spot 4), pyridoxine biosynthesis PDX1-like protein 3 (PDX1, spot11), and alcohol dehydrogenase (ADH, spots 20) were also found to respond to high salinity and the NaCl and H2S treatment. An increase in the CP activity occurs just before the first senescence symptoms became visible in most plants [49]. CP is involved in PCD, and the balance between the levels of CP and phytocystatins is kept in a certain state through the antagonistic activities of abscisic acid (ABA) and gibberellins [50,51]. In our study, the accumulation of CP was significantly down-regulated by the NaHS addition. PDX1 is regarded as pyridoxal phosphate-binding protein [52]. Pyridoxal phosphate-binding proteins protect plants from high ion concentrations by inducing cysteine biosynthesis or the synthesis of cystatin to inhibit PCD [53]. We observed that PDX1 was significantly up-regulated to inhibit PCD under the NaCl and H2S treatment. ADH catalyzes the regeneration of NAD+ from the reduction of acetaldehyde to ethanol [54]. H2S has been shown to inhibit ADH activity in the roots of the freshwater marsh species Spartina alterniflora [16], but the accumulation of ADH was down-regulated under the NaCl treatment in our study and did not change significantly with the NaHS addition. We observed that H2S helped alleviate the reduction of ADH under the NaCl treatment.

3.5. H2S Regulates Hormone Biosynthesis and Transcription Factor

Abscisic acid stress ripening protein (Asr, spot 31) is a downstream protein that participates in ABA signaling pathways [55]. The bHLH transcription factors serve as a negative feedback regulatory loop in ABA signaling in Arabidopsis thaliana [56]. In the present study, Asr and bHLH145 were up-regulated under the salt treatment and did not change significantly under the NaCl and H2S treatment, which indicated the signal transduction of ABA was inhibited by the addition of NaHS. Furthermore, in the promoter regions of the genes encoding CP exist ABA-responsive elements, and the transcription of CP is repressed by ABA through DOF (DNA binding One Zinc Finger) [51]. Additionally, NaHS up-regulated the accumulation of basic region-leucine zipper (bZIP, spot 27), which has been found to positively alter the adaptation of plants to adverse circumstances [57,58].

3.6. Protein–Protein Interaction Networks Analysis

The DEPs in the PPI network exhibited a higher enrichment, indicating a higher degree of modularization; therefore, the DEPs were divided into two modules to investigate their interactions by using the Molecular Complex Detection (MCODE) (Figure 5). In module A, chloroplastic drought-induced stress protein (CDSP), NAD(P)-binding Rossmann-fold superfamily protein isoform 1 (FLDH, spot 24), FNR, GS2, OEE, PC, PETC, and RBCS were enriched in response to stimulus. The key node in this module was GS2. In module B, six nodes (APX, Cpn20, TPI, PGK, 2-Cys peroxiredoxin (Prx), and SOD) were involved in response to abiotic stimulus, and another two nodes (RBCL, ATP5F1B) were related to the Calvin cycle. The two modules implied that response to abiotic stimulus influenced other biological processes.

4. Materials and Methods

4.1. Plant Growth and Treatment

Mature propagules of K. obovata were collected from the mangrove forest in the National Nature Reserve for Mangroves in Zhangjiang River Estuary (23°55′ N, 117°26′ E), Yunxiao County, Fujian Province, PR China, where salinity levels range from 8% to 20%. Healthy propagules of similar size (20 cm in length) were selected and pre-cultivated in a plastic pot (dimension of 12 cm in diameter and 11 cm in depth) with clean sands. The pots were placed in a growth chamber with a temperature of 25/30 °C (night/day), a relative humidity of 70%, and a photoperiod of 8 h dark/16 h light with around 1250 μmol m−2 s−1. A 1/4 strength Hoagland’s nutrient solution was used to cultivate the hypocotyls. The solutions were replaced every week.
The treatments of NaHS and NaCl were set up when the 3rd pair of leaves appeared for physiological assays. The healthy and uniform seedlings were randomly cultivated in three groups, and each group had three replicates.
In the first group, plants were supplied with NaHS and NaCl for 7 days. Different concentrations of NaHS at 0, 50, 100, 200, 350 and 500 µM (The average of H2S was 243.1 ± 234.9 μM in the pore water samples of mangrove forest sediment) [59] were prepared in the 1/4 strength Hoagland’s nutrient solution in the presence or absence of 400 mM of NaCl. After 7 days of treatment, photosynthesis, chlorophyll content, and dry weight were measured to determine the optimal concentration of NaHS for the salinity-induced growth inhibition.
The second group was supplied with a 1/4 strength Hoagland’s nutrient solution containing 400 mM of NaCl in the presence of NaHS (200 μM), Na2SO4 (200 μM), NaHSO3 (200 μM), NaHSO4 (200 μM) and CH3COONa (200 μM). After 7 days of treatments, the contents of chlorophyll and endogenous H2S were measured to distinguish the actual role of H2S/HS from the other compounds.
According to the results from above two experiments, the third group was supplied with a 1/4 strength Hoagland’s nutrient solution containing both 400 mM of NaCl (NaCl) and 400 mM of NaCl and 200 μM of NaHS (NaCl and H2S). The same volume of a 1/4 strength Hoagland’s nutrient solution was used as the control (CK). After chlorophyll fluorescence measurements, the second pair of leaves from the apex of the growing shoots was harvested, rapidly frozen in liquid nitrogen, and stored at −80 °C for a proteomic study and further biochemical analysis.

4.2. Measurements of Leaf Photosynthetic Pigment, Photosynthetic Rate and Chlorophyll Fluorescence Quenching

The contents of chlorophyll were measured according to the method of Lichtenthaler with little modifications [60]. Eighty percent (v/v) acetone was used to extract chlorophyll, and the content was estimated according to the absorbances at 470, 646 and 663 nm.
The net photosynthetic rate (Pn) of fully expanded leaves, as well as transpiration (Tr), stomatal conductance (Gs), and internal CO2 concentration (Ci), were determined by using a portable photosynthesis system (Li-6400, Li-Cor, Lincoln, NE, USA). Those measurements were carried out in the morning from 9:00 to 11:30. The recovery coefficient of Pn was calculated according to the following equation: PnNaCl/PnCK. Here, PnNaCl and PnCK stand for the leaf photosynthetic rate measured in the NaCl and CK treatments.
We applied a fluorometer for measuring chlorophyll fluorescence (Li-6400, Li-Cor, Lincoln, NE, USA). Fv/Fm was calculated by using (FmF0)/Fm, where minimum (dark) fluorescence (F0) was obtained by applying to measure light pulses at low frequency (0.03 μmol m−2 s−1 for 1 s). The maximum fluorescence (Fm) was determined by applying a saturating light pulse (6000 μmol m−2 s−1 for 0.8 s) to a dark-adapted sample [61].

4.3. Determination of Dry Weight and the Content of Endogenous H2S

Endogenous H2S content was determined according to the procedure of Zhang et al. [62], and the method of Chen et al. was employed to determine the dry weight [63].

4.4. Measurements of Oxidative Stress and Antioxidant System Activity

The H2O2 content was measured according to the method of Hung and Kao [64]. Lipid peroxidation was determined by measuring thiobarbituric acid reactive substances (TBARS), according to the method of Yan et al. [65].
For the relative electrolyte leakage (REL) measurement, fresh leaves (0.2 g) were cut, added to deionized water (20 mL), and degassed for 10 min at room temperature. Using an electrical conductivity meter (DDS-11A), the initial reading of conductivity (E1) was recorded. Then, the solution containing plant materials was incubated at 100 °C for 15 min. The final reading of conductivity (E2) of the solution was recorded after cooling. Additionally, the reading of conductivity for deionized water (E0) was measured. Finally, the REL was estimated based on the following equation: REL (%) = (E1 − E0)/(E2 − E0) × 100%.
Superoxide dismutase (SOD) activity was determined according to the method described by Beauchamp and Fridovich [66]. Ascorbate peroxidase (APX) activity was estimated according to the method of Chen and Asada [67].
The total content of glutathione (GSH) was measured by a GSH kit (JBI, Nanjing, China). Leaves (0.3 g) were ground with 0.3 mL of H3PO3 (25%) and 0.9 mL of an NaH2PO4-EDTA buffer (0.1 M, pH 8.0). The homogenates were separated at 10,000× g at 4 °C for 20 min. The supernatant was used for GSH content measurement [68]. The absorbance at 420 nm was measured according to the manufacturer’s instructions.

4.5. Protein Extraction from K. Obovata Leaves and 2-DE PAGE and Image/Data Analysis

Proteins were extracted by a phenol extraction procedure and followed by methanolic ammonium acetate precipitation, according to the method described by Delaplace et al. [69]. We rehydrated 17 cm immobilized pH 4–7 gradient strips (BioRad, CA, USA) with 340 μL of a rehydration buffer (containing a 1 mg protein sample) in the tray overnight. An Ettan immobilized pH gradient (IPG)phor3 system (GE Healthcare Amersham Bioscience, Little Chalfont, UK) was used for isoelectric focusing (IEF), following the conditions described by Shen et al. [70]. In 12.5% acrylamide gels, we observed gel electrophoresis. Each experiment was repeated three times. Each treatment was performed for three biological replications for quantitative analysis.
We stained SDS-PAGE gels with Coomassie Brilliant Blue R-250, and these gels were scanned by Uniscan M3600 (China) at 600 dpi. Gels were analyzed by the PDQuest software (Version 8.0, Bio-Rad) [71]. Afterward, DEPs were obtained by pairwise comparison with a fold change ≥2.0 and a Student’s t-test (p < 0.05).

4.6. In-Gel Digestion, Identification and Classification of Differentially Expressed Proteins

We excised the DEPs with a more than 2.0-fold change from 2-DE gels [70]. Gel slices were incubated with 25 mM of NH4HCO3 and acetonitrile for four times to remove the color of coomassie brilliant blue-R250 (CBB-R250). Trypsin/Lys-C Mix (Promega, Madison, WI, USA) was used to digest the proteins. Trifluoroacetic acid (0.5–1%) was added to terminate digestion. The supernatant was collected at 12,000× g for 10 min. A MALDI-TOF-TOF mass spectrometer was used to analyze the supernatant for mass spectrometry identification (Applied Biosystems, Massachusetts, USA), according to the method of Hu et al. [71]. The search parameters were set as follows: database NCBInr (release date: 2018.12.01); taxonomy Viridiplantae (green plants); peptide mass ranged from 10 to 130 kDa; the coverage of protein sequence had to reach a minimum of 10%; proteins had scores higher than 60 (p < 0.05); results with a confidence interval % (C.I.%) value higher than 95% were considered to be an identification.
The functions of the identified proteins were determined in http://www.uniprot.org/uniprot, and then the proteins were divided into five groups according to their biological functions in the plant. The subcellular localization of identified proteins was obtained by the protein subcellular localization prediction tool (WoLF PSORT, http://wolfpsort.hgc.jp/), as well as from previously published papers when possible.
The STRING database was used to predict the protein–protein interactions (PPIs) [72], and Cytoscape was used to visualize significant protein–protein associations in the PPI network. A combined score of >0.4 was selected as the cut-off value. Then, module analysis was carried out by the Molecular Complex Detection (MCODE) plugin to highlight the biological significance of gene modules that respond to the NaCl and/or NaHS treatment in K. obovata leaves [73]. Subsequently, the BiNGO tool was used to visualize the level of enriched GO terms.

4.7. Quantitative Real-Time PCR Analysis

For the total RNA extraction, the frozen K. obovata leaves (0.1g) were ground in liquid nitrogen and extracted by using a Total RNA Kit (TaKaRa, Dalian, China). We observed the RNA quality and integrity with an ultraviolet spectrophotometer (Cary 50, Varian, USA) and agarose gel electrophoresis. The RNA was used to synthesized cDNAs with moloney murine leukemia virus (M-MLV) reverse transcriptase First-Strand cDNA synthesis kit (TaKaRa, Dalian, China), and the cDNA mixture was used as a template for subsequent PCRs. The primers used for real-time PCR are shown in Supplementary Table S4. A 10 μL real-time PCR mixture contained 2 μL of primers, 2 μL of cDNA, and 6 μL of SYBR Green (Sangon, Shanghai, China). Five independent biological replicates were used to perform gene expression. The relative gene expression was calculated by the 2−ΔΔCT method, and actin was used as an an internal control [74]. The Bio-Rad iQ5 Multicolor Real-Time PCR Detection System (Bio-Rad, Hercules, CA) was used to run qRT-PCR.

4.8. Statistical Analysis

A two-tailed Student’s t-test (p < 0.05) and a Duncan’s multiple range test were used for statistical analysis. Pearson’s correlation among the amendment’s characteristics was evaluated with SPSS Statistics for Windows (Version 22.0, IBM Corp, Armonk, NY). The direct and indirect relationships among leaf drought weight, chlorophyll content, net photosynthetic rate, intercellular carbon dioxide, stomatal conductance and transpiration rate were evaluated through path analysis by using edgeR package. The heatmap of the DEPs was drawn by using the R software “heatmap package.”

5. Conclusions

Physiological and proteomic evidence in this study support beneficial role of H2S on K. obovata exposed to high salinity: (1) H2S increased carbon fixations and electron transfer, maintaining photosynthesis under high salinity; (2) H2S stimulated glycolysis to generate precursor metabolites, energy, and protein synthesis; (3) H2S increased non-enzymatic antioxidants to detoxify ROS and prevent cellular damage under high salinity; and (4) H2S affected the ABA signaling pathway, which led to a sustained plant adaptation to high salinity. As shown in Figure 7, we propose a working model to illustrate the detailed mechanism by which H2S alleviates NaCl-induced inhibition. In conclusion, these findings have important implications for understanding the functional role of H2S in the salt tolerance of K. obovata.

Supplementary Materials

Supplementary Materials can be found at https://www.mdpi.com/1422-0067/21/1/118/s1; Figure S1. Net photosynthetic rate (Pn) recovery coefficient of K. obovata seedling leaf under different concentrations of NaHS; Figure S2. H2S and/or HS, but not other compounds derived from NaHS, contribute to (A) chlorophyll content and (B) endogenous H2S accumulation; Figure S3. Effects of NaHS on the chlorophyll fluorescence characteristics of K. obovata seedling leaves; Figure S4. Protein–protein interactions (PPIs) network analysis of 27 identified proteins from the NaCl- and NaHS-treated K. obovata leaves; Figure S5. Gene ontology analysis of the proteins identified from the leaves of K. obovate; Table S1. Co-relations among different physiological indices in K. obovata leaves; Table S2. The detail information of the peptide sequence, expected and observed m/z values, and sequence coverage for the different expressed proteins identified by MALDI-TOF/TOF MS in the leaves of K. obovata under the salinity and NaHS treatments; Table S3. The modules of the protein–protein interaction networks and nodes biological process (GO) analysis; Table S4. Sequences of forward and reverse primers used in qPCR for the gene expression analysis in K. obovata leaves.

Author Contributions

Conceptualization, Z.-J.S.; Data curation, Y.-L.L., H.L. and D.-N.M.; Formal analysis, Y.-L.L.; Project administration, M.S.; Writing—original draft, X.-Y.Z.; Writing—review & editing, Y.-L.L. and H.-L.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This study was jointly supported by the funds from the National Key Research and Development Program of China (2017YFC0506102), and the Natural Science Foundation of China (NSFC 31570586 and 31870581).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Parida, A.K.; Jha, B. Salt tolerance mechanisms in mangroves: A review. Trees-Struct. Funct. 2010, 24, 199–217. [Google Scholar]
  2. Jayatissa, L.P.; Wickramasinghe, W.; Dahdouh-Guebas, F.; Huxham, M. Interspecific variations in responses of mangrove seedlings to two contrasting salinities. Int. Rev. Hydrobiol. 2008, 93, 700–710. [Google Scholar] [CrossRef] [Green Version]
  3. Yang, S.C.; Shih, S.S.; Hwang, G.W.; Adams, J.B.; Lee, H.Y.; Chen, C.P. The salinity gradient influences on the inundation tolerance thresholds of mangrove forests. Ecol. Eng. 2013, 51, 59–65. [Google Scholar] [CrossRef]
  4. Li, N.; Chen, S.L.; Zhou, X.Y.; Li, C.Y.; Shao, J.; Wang, R.G.; Fritz, E.; Hüttermann, A.; Polle, A. Effect of NaCl on photosynthesis, salt accumulation and ion compartmentation in two mangrove species, Kandelia candel and Bruguiera gymnorhiza. Aquat. Bot. 2008, 88, 303–310. [Google Scholar] [CrossRef]
  5. Watson, E.B.; Szura, K.; Wigand, C.; Raposa, K.B.; Blount, K.; Cencer, M. Sea level rise, drought and the decline of Spartina patens in New England marshes. Biol. Conserv. 2016, 196, 173–181. [Google Scholar] [CrossRef]
  6. Luo, M.; Huang, J.F.; Zhu, W.F.; Tong, C. Impacts of increasing salinity and inundation on rates and pathways of organic carbon mineralization in tidal wetlands: A review. Hydrobiologia 2019, 827, 31–49. [Google Scholar] [CrossRef]
  7. Ganguly, D.; Ray, R.; Majumdar, N.; Chowdhury, C.; Jana, T.K. Biogenic hydrogen sulphide emissions and non-sea sulfate aerosols over the Indian Sundarban mangrove forest. J. Atmos. Chem. 2018, 75, 319–333. [Google Scholar]
  8. Lyimo, T.J.; Pol, A.; den Camp, H.J.O. Sulfate reduction and methanogenesis in sediments of Mtoni mangrove forest, Tanzania. Ambio 2002, 31, 614–617. [Google Scholar] [CrossRef]
  9. Hou, Z.H.; Liu, J.; Hou, L.X.; Li, X.D.; Liu, X. H2S may function downstream of H2O2 in jasmonic acid-induced stomatal closure in Vicia faba. Chin. Bull. Bot. 2011, 46, 396–406. [Google Scholar]
  10. Li, Z.G.; Gong, M.; Xie, H.; Yang, L.; Li, J. Hydrogen sulfide donor sodium hydrosulfide-induced heat tolerance in tobacco (Nicotiana tabacum L.) suspension cultured cells and involvement of Ca2+ and calmodulin. Plant Sci. 2012, 185, 185–189. [Google Scholar] [CrossRef]
  11. Yarsi, G.; Sivaci, A.; Dasgan, H.Y.; Altuntas, O.; Binzet, R.; Akhoundnejad, Y. Effects of salinity stress on chlorophyll and carotenoid contents and stomata size of grafted and ungrafted galia C8 melon cultivar. Pak. J. Bot. 2017, 49, 421–426. [Google Scholar]
  12. Yu, L.X.; Zhang, C.J.; Shang, H.Q.; Wang, X.F.; Min, W.E.I.; Yang, F.J.; Shi, Q.H. Exogenous hydrogen sulfide enhanced antioxidant capacity, amylase activities and salt tolerance of cucumber hypocotyls and radicles. J. Integr. Agric. 2013, 12, 445–456. [Google Scholar] [CrossRef]
  13. Sun, Y.D.; Luo, W.R. Effects of exogenous hydrogen sulphide on seed germination and seedling growth of cucumber (Cucumis sativus) under sodium bicarbonate stress. Seed Sci. Technol. 2014, 42, 126–131. [Google Scholar] [CrossRef]
  14. Christou, A.; Manganaris, G.A.; Papadopoulos, I.; Fotopoulos, V. Hydrogen sulfide induces systemic tolerance to salinity and non-ionic osmotic stress in strawberry plants through modification of reactive species biosynthesis and transcriptional regulation of multiple defence pathways. J. Exp. Bot. 2013, 64, 1953–1966. [Google Scholar] [CrossRef]
  15. Mostofa, M.G.; Saegusa, D.; Fujita, M.; Tran, L.S.P. Hydrogen sulfide regulates salt tolerance in rice by maintaining Na+/K+ balance, mineral homeostasis and oxidative metabolism under excessive salt stress. Front. Plant Sci. 2015, 6, 1055. [Google Scholar] [CrossRef] [Green Version]
  16. Koch, M.S.; Mendelssohn, I.A.; Mckee, K.L. Mechanism for the hydrogen sulfide-induced growth limitation in wetland macrophytes. Limnol. Oceanogr. 1990, 35, 399–408. [Google Scholar] [CrossRef] [Green Version]
  17. Krishnamurthy, P.; Tan, X.F.; Lim, T.K.; Lim, T.M.; Kumar, P.P.; Loh, C.S.; Lin, Q. Proteomic analysis of plasma membrane and tonoplast from the leaves of mangrove plant Avicennia officinalis. Proteomics 2014, 14, 2545–2557. [Google Scholar] [CrossRef]
  18. Maxwell, K.; Johnson, G.N. Chlorophyll fluorescence-a practical guide. J. Exp. Bot. 2000, 51, 659–668. [Google Scholar] [CrossRef]
  19. Hosoki, R.; Matsuki, N.; Kimura, H. The possible role of hydrogen sulfide as an endogenous smooth muscle relaxant in synergy with nitric oxide. Biochem. Biophys. Res. Co. 1997, 237, 527–531. [Google Scholar] [CrossRef]
  20. Chen, J.; Wu, F.H.; Wang, W.H.; Zheng, C.J.; Lin, G.H.; Dong, X.J.; He, J.X.; Pei, Z.M.; Zheng, H.L. Hydrogen sulphide enhances photosynthesis through promoting chloroplast biogenesis, photosynthetic enzyme expression, and thiol redox modification in Spinacia oleracea seedlings. J. Exp. Bot. 2011, 62, 4481–4493. [Google Scholar] [CrossRef] [Green Version]
  21. Flexas, J.; Medrano, H. Drought-inhibition of photosynthesis in C3 plants: Stomatal and non-stomatal limitations revisited. Ann. Bot-London 2002, 89, 183–189. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Mulo, P.; Medina, M. Interaction and electron transfer between ferredoxin-NADP+ oxidoreductase and its partners: Structural, functional, and physiological implications. Photosynth. Res. 2017, 134, 265–280. [Google Scholar] [CrossRef] [PubMed]
  23. Sugihara, K.; Hanagata, N.; Dubinsky, Z.; Baba, S.; Karube, I. Molecular characterization of cDNA encoding oxygen evolving enhancer protein 1 increased by salt treatment in the mangrove Bruguiera gymnorrhiza. Plant Cell Physiol. 2000, 41, 1279–1285. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Sanda, S.; Yoshida, K.; Kuwano, M.; Kawamura, T.; Munekage, Y.N.; Akashi, K.; Yokota, A. Responses of the photosynthetic electron transport system to excess light energy caused by water deficit in wild watermelon. Physiol. Plant 2011, 142, 247–264. [Google Scholar] [CrossRef] [PubMed]
  25. Wang, L.; Pan, D.Z.; Li, J.; Tan, F.; Hoffmann-Benning, S.; Liang, W.; Chen, W. Proteomic analysis of changes in the Kandelia candel chloroplast proteins reveals pathways associated with salt tolerance. Plant Sci. 2015, 231, 159–172. [Google Scholar] [CrossRef] [PubMed]
  26. Wang, D.; Lu, Q.; Li, X.F.; Jiang, Q.S.; Wu, J.X.; Jiang, D.A. Relationship between Rubisco activase isoform levels and photosynthetic rate in different leaf positions of rice plant. Photosynthetica 2009, 47, 621–629. [Google Scholar] [CrossRef]
  27. Bauwe, H.; Hagemann, M.; Kern, R.; Timm, S. Photorespiration has a dual origin and manifold links to central metabolism. Curr. Opin. Plant Biol. 2012, 15, 269–275. [Google Scholar] [CrossRef]
  28. Wingler, A.; Lea, P.J.; Quick, W.P.; Leegood, R.C. Photorespiration: Metabolic pathways and their role in stress protection. Philos. Trans. R. Soc. London Ser. B 2000, 355, 1517–1529. [Google Scholar] [CrossRef] [Green Version]
  29. Matuszyńska, A.; Saadat, N.P.; Ebenhoeh, O. Balancing energy supply during photosynthesis–a theoretical perspective. Physiol. Plant 2019, 166, 392–402. [Google Scholar] [CrossRef] [Green Version]
  30. Zhao, Q.; Suo, J.W.; Chen, S.X.; Jin, Y.D.; Ma, X.L.; Yin, Z.P.; Zhang, Y.H.; Wang, T.; Luo, J.; Jin, W.H.; et al. Na2CO3–responsive mechanisms in halophyte Puccinellia tenuiflora roots revealed by physiological and proteomic analyses. Sci. Rep. 2016, 6, 32717. [Google Scholar] [CrossRef] [Green Version]
  31. Schultz, C.J.; Hsu, M.; Miesak, B.; Coruzzi, G.M. Arabidopsis mutants: Define an in vivo role for isoenzymes of aspartate aminotransferase in plant nitrogen assimilation. Genetics 1998, 149, 491–499. [Google Scholar] [PubMed]
  32. Lea, P.J.; Miflin, B.J. Glutamate synthase and the synthesis of glutamate in plants. Plant Physiol. Bioch. 2003, 41, 555–564. [Google Scholar] [CrossRef]
  33. Sarry, J.E.; Kuhn, L.; Ducruix, C.; Lafaye, A.; Junot, C.; Hugouvieux, V.; Jourdain, A.; Bastien, O.; Fievet, J.B.; Vailhen, D.; et al. The early responses of Arabidopsis thaliana cells to cadmium exposure explored by protein and metabolite profiling analyses. Proteomics 2006, 6, 2180–2198. [Google Scholar] [CrossRef] [PubMed]
  34. Semane, B.; Dupae, J.; Cuyper, A.; Noben, J.P.; Tuomainen, M.; Tervahauta, A.; Kärenlampi, S.; Van Belleghem, F.; Smeets, K.; Vangronsveld, J. Leaf proteome responses of Arabidopsis thaliana exposed to mild cadmium stress. J. Plant Physiol. 2010, 167, 247–254. [Google Scholar] [CrossRef]
  35. Omar, S.A.; Fu, Q.T.; Chen, M.S.; Wang, G.J.; Song, S.Q.; Elsheery, N.I.; Xu, Z.F. Identification and expression analysis of two small heat shock protein cDNAs from developing seeds of biodiesel feedstock plant Jatropha curcas. Plant Sci. 2011, 181, 632–637. [Google Scholar] [CrossRef]
  36. Yasuda, R.; Noji, H.; Yoshida, M.; Kinosita, J.K.; Itoh, H. Resolution of distinct rotational substeps by submillisecond kinetic analysis of F-1-ATPase. Nature 2001, 410, 898–904. [Google Scholar] [CrossRef]
  37. Houštěk, J.; Pícková, A.; Vojtíšková, A.; Mráček, T.; Pecina, P.; Ješina, P. Mitochondrial diseases and genetic defects of ATP synthase. BBA-Bioenergetics 2006, 1757, 1400–1405. [Google Scholar] [CrossRef] [Green Version]
  38. Bertini, L.; Palazzi, L.; Proietti, S.; Arrigoni, G.; Polverino de Laureto, P.; Caruso, C. Proteomic Analysis of MeJa-Induced Defense Responses in Rice against Wounding. Int. J. Mol. Sci. 2019, 20, 2525. [Google Scholar] [CrossRef] [Green Version]
  39. Dorion, S.; Rivoal, J. Clues to the functions of plant NDPK isoforms. N-S Arch. Pharmacol. 2015, 388, 119–132. [Google Scholar] [CrossRef]
  40. Li, D.; Limwachiranon, J.; Li, L.; Du, R.; Luo, Z. Involvement of energy metabolism to chilling tolerance induced by hydrogen sulfide in cold-stored banana fruit. Food Chem. 2016, 208, 272–278. [Google Scholar] [CrossRef]
  41. Hernandez, J.A.; Almansa, M.S. Short-term effects of salt stress on antioxidant systems and leaf water relations of pea leaves. Physiol. Plant 2002, 115, 251–257. [Google Scholar] [CrossRef] [PubMed]
  42. Demiral, T.; Turkan, I. Comparative lipid peroxidation, antioxidant defense systems and proline content in roots of two rice cultivars differing in salt tolerance. Environ. Exp. Bot. 2005, 53, 247–257. [Google Scholar] [CrossRef]
  43. Shi, H.; Ye, T.; Chan, Z. Exogenous application of hydrogen sulfide donor sodium hydrosulfide enhanced multiple abiotic stress tolerance in bermudagrass (Cynodon dactylon (L). Pers.). Plant Physiol. Bioch. 2013, 71, 226–234. [Google Scholar]
  44. Geitmann, A.; Nebenfuhr, A. Navigating the plant cell: Intracellular transport logistics in the green kingdom. Mol. Biol. Cell 2015, 26, 3373–3378. [Google Scholar] [CrossRef] [PubMed]
  45. Thomas, S.G.; Huang, S.J.; Li, S.T.; Staiger, C.J.; Franklin-Tong, V.E. Actin depolymerization is sufficient to induce programmed cell death in self-incompatible pollen. J. Cell Biol. 2006, 174, 221–229. [Google Scholar] [CrossRef] [Green Version]
  46. Yin, Z.P.; Zhang, H.; Zhao, Q.; Yoo, M.J.; Zhu, N.; Yu, J.L.; Yu, J.J.; Guo, S.Y.; Miao, Y.C.; Chen, S.X.; et al. Physiological and comparative proteomic analyses of saline-alkali NaHCO3-responses in leaves of halophyte Puccinellia tenuiflora. Plant Soil. 2019, 437, 137–158. [Google Scholar] [CrossRef]
  47. Suo, J.W.; Zhao, Q.; Zhang, Z.X.; Chen, S.X.; Cao, J.G.; Liu, G.J.; Wei, X.; Wang, T.; Chuan, P.; Dai, S.J. Cytological and proteomic analyses of osmunda cinnamomea germinating spores reveal characteristics of fern spore germination and rhizoid tip growth. Mol. Cell Proteomics 2015, 14, 2510–2534. [Google Scholar] [CrossRef] [Green Version]
  48. Abreu, I.A.; Cabelli, D.E. Superoxide dismutases-a review of the metal-associated mechanistic variations. BBA-Proteins Proteom 2010, 1804, 263–274. [Google Scholar] [CrossRef]
  49. Buono, R.A.; Hudecek, R.; Nowack, M.K. Plant proteases during developmental programmed cell death. J. Exp. Bot. 2019, 70, 2097–2112. [Google Scholar] [CrossRef]
  50. Prins, A.; Van Heerden, P.D.; Olmos, E.; Kunert, K.J.; Foyer, C.H. Cysteine proteinases regulate chloroplast protein content and composition in tobacco leaves: A model for dynamic interactions with ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) vesicular bodies. J. Exp. Bot. 2008, 59, 1935–1950. [Google Scholar] [CrossRef] [Green Version]
  51. Szewińska, J.; Simińska, J.; Bielawski, W. The roles of cysteine proteases and phytocystatins in development and germination of cereal seeds. J. Plant Physiol. 2016, 207, 10–21. [Google Scholar] [CrossRef]
  52. Komatsu, S.; Kamal, A.H.; Hossain, Z. Wheat proteomics: Proteome modulation and abiotic stress acclimation. Front Plant Sci. 2014, 5, 684. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Youssefian, S.; Nakamura, M.; Sano, H. Tobacco plants transformed with the O-acetylserine (thiol) lyase gene of wheat are resistant to toxic levels of hydrogen sulphide gas. Plant J. 1993, 4, 759–769. [Google Scholar] [CrossRef] [PubMed]
  54. Huang, Y.; Picha, D.H.; Kilili, A.W. Atmospheric oxygen level influences alcohol dehydrogenase and pyruvate decarboxylase activities in sweet potato roots. J. Plant Physiol. 2002, 159, 129–136. [Google Scholar] [CrossRef]
  55. Shkolnik, D.; Bar-Zvi, D. Tomato ASR1 abrogates the response to abscisic acid and glucose in Arabidopsis by competing with AB14 for DNA binding. Plant Biotechnol. J. 2008, 6, 368–378. [Google Scholar] [CrossRef]
  56. Zheng, K.; Wang, Y.; Wang, S. The non-DNA binding bHLH transcription factor Paclobutrazol Resistances are involved in the regulation of ABA and salt responses in Arabidopsis. Plant Physiol. Bioch. 2019, 139, 239–245. [Google Scholar] [CrossRef]
  57. Joo, J.; Lee, Y.H.; Kim, Y.K.; Nahm, B.H.; Song, S.I. Abiotic stress responsive rice ASR1 and ASR3 exhibit different tissue-dependent sugar and hormone-sensitivities. Mol. Cells 2013, 35, 421–435. [Google Scholar] [CrossRef] [Green Version]
  58. Sen, S.; Chakraborty, J.; Ghosh, P.; Basu, D.; Das, S. Chickpea WRKY70 regulates the expression of a homeodomain-leucine zipper (HD-Zip) I rranscription factor CaHDZ12, which confers abiotic stress tolerance in transgenic tobacco and chickpea. Plant Cell Physiol. 2017, 58, 1934–1952. [Google Scholar] [CrossRef]
  59. Pérez, A.; Gutiérrez, D.; Saldarriaga, M.S.; Sanders, C.J. Tidally driven sulfidic conditions in Peruvian mangrove sediments. Geo-Mar. Lett. 2018, 38, 457–465. [Google Scholar]
  60. Lichtenthaler, H.K. Chlorophylls and carotenoids: Pigments of photosynthetic biomembranes. Method Enzymol. 1987, 148, 350–382. [Google Scholar]
  61. Neves, L.H.; Santos, R.I.N.; dos Santos Teixeira, G.I.; de Araujo, D.G.; Silvestre, W.V.D.; Pinheiro, H.A. Leaf gas exchange, photochemical responses and oxidative damages in assai (Euterpe oleracea Mart.) seedlings subjected to high temperature stress. Sci. Hortic. 2019, 257, 108733. [Google Scholar] [CrossRef]
  62. Zhang, H.; Hu, L.Y.; Hu, K.D.; He, Y.D.; Wang, S.H.; Luo, J.P. Hydrogen sulfide promotes wheat seed germination and alleviates oxidative damage against copper stress. J. Integr. Plant Biol. 2008, 50, 1518–1529. [Google Scholar] [CrossRef] [PubMed]
  63. Chen, J.; Liu, T.W.; Hu, W.J.; Simon, M.; Wang, W.H.; Chen, J.; Liu, X.; Zheng, H.L. Comparative proteomic analysis of differentially expressed proteins induced by hydrogen sulfide in Spinacia oleracea leaves. PLoS ONE 2014, 9, e105400. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Hung, K.T.; Kao, C.H. Hydrogen peroxide is necessary for abscisic acid-induced senescence of rice leaves. J. Plant Physiol. 2004, 161, 1347–1357. [Google Scholar] [CrossRef] [PubMed]
  65. Yan, K.; Chen, W.; Zhang, G.; Xu, S.; Liu, Z.; He, X.; Wang, L. Elevated CO2 ameliorated oxidative stress induced by elevated O3 in Quercus mongolica. Acta Physiol. Plant 2010, 32, 375–385. [Google Scholar] [CrossRef]
  66. Beauchamp, C.; Fridovich, I. Superoxide dismutase: Improved assays and an assay applicable to acrylamide gels. Anal. Biochem. 1971, 44, 276–287. [Google Scholar] [CrossRef]
  67. Chen, G.X.; Asada, K. Ascorbate peroxidase in tea leaves: Occurrence of two isozymes and the differences in their enzymatic and molecular properties. Plant Cell Physiol. 1989, 30, 987–998. [Google Scholar]
  68. Devi, S.R.; Prasad, M.N.V. Copper toxicity in Ceratophyllum demersum L. (Coontail), a free floating macrophyte: Response of antioxidant enzymes and antioxidants. Plant Sci. 1998, 138, 157–165. [Google Scholar] [CrossRef]
  69. Delaplace, P.; van der Wal, F.; Dierick, J.F.; Cordewener, J.H.G.; Fauconnier, M.L.; du Jardin, P.; America, A.H.R. Potato tuber proteomics: Comparison of two complementary extraction methods designed for 2-DE of acidic proteins. Proteomics 2006, 6, 6494–6497. [Google Scholar] [CrossRef]
  70. Shen, Z.J.; Chen, J.; Ghoto, K.; Hu, W.J.; Gao, G.F.; Luo, M.R.; Li, Z.; Simon, M.; Zhu, X.Y.; Zheng, H.L. Proteomic analysis on mangrove plant Avicennia marina leaves reveals nitric oxide enhances the salt tolerance by up-regulating photosynthetic and energy metabolic protein expression. Tree Physiol. 2018, 38, 1605–1622. [Google Scholar] [CrossRef]
  71. Hu, W.J.; Chen, J.; Liu, T.W.; Wu, Q.; Wang, W.H.; Liu, X.; Shen, Z.J.; Simon, M.; Chen, J.; Wu, F.H.; et al. Proteome and calcium-related gene expression in Pinus massoniana needles in response to acid rain under different calcium levels. Plant Soil 2014, 380, 285–303. [Google Scholar] [CrossRef]
  72. Szklarczyk, D.; Franceschini, A.; Wyder, S.; Forslund, K.; Heller, D.; Huerta-Cepas, J.; Simonovic, M.; Roth, A.; Santos, A.; Tsafou, K.P.; et al. STRING v10: Protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015, 43, D447–D452. [Google Scholar] [CrossRef] [PubMed]
  73. Bader, G.D.; Hogue, C.W. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinform. 2003, 4, 2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Effects of various concentrations of NaHS on the leaf photosynthetic characteristics of Kandelia obovata seedlings treated by high salinity. (A) leaf dry weight, (B) total chlorophyll contents, (C) net photosynthetic rate (Pn), (D) intercellular CO2 concentration (Ci), (E) stomatal conductance (Gs) and (F) transpiration rate (Tr). Seedlings were treated with different concentrations of NaHS (0, 50, 100, 200, 350 and 500 μM) together with 400 mM of NaCl for seven days. The sample without 400 mM of NaCl was regarded as the experimental control group (CK), and the sample with 400 mM of NaCl was regarded as the NaCl group. Values are means ± SE of three independent experiments with at least three replicates for each. Bars with different letters are significantly different at p < 0.05 according to Duncan’s multiple range test.
Figure 1. Effects of various concentrations of NaHS on the leaf photosynthetic characteristics of Kandelia obovata seedlings treated by high salinity. (A) leaf dry weight, (B) total chlorophyll contents, (C) net photosynthetic rate (Pn), (D) intercellular CO2 concentration (Ci), (E) stomatal conductance (Gs) and (F) transpiration rate (Tr). Seedlings were treated with different concentrations of NaHS (0, 50, 100, 200, 350 and 500 μM) together with 400 mM of NaCl for seven days. The sample without 400 mM of NaCl was regarded as the experimental control group (CK), and the sample with 400 mM of NaCl was regarded as the NaCl group. Values are means ± SE of three independent experiments with at least three replicates for each. Bars with different letters are significantly different at p < 0.05 according to Duncan’s multiple range test.
Ijms 21 00118 g001
Figure 2. Effects of NaHS on oxidative stress and antioxidant system activity of the leaves of K. obovata seedlings treated by high salinity. (A) H2O2 content, (B) thiobarbituric acid reactive substances (TBARS) accumulation, (C) electrolyte leakage percentage (REL), (D) glutathione (GSH) content, (E) ascorbate acid peroxidase (APX) activity and (F) superoxide dismutase (SOD) activity. CK stands for the control treated only by a 1/4 strength Hoagland’s nutrient solution. NaCl stands for the 400 mM of the NaCl treatment. NaCl and H2S stands for the treatment with 400 mM of NaCl and 200 μM of NaHS. Error bars are SE (n = 3). The columns labeled with different letters are significantly different at p < 0.05 according to Duncan’s multiple range test.
Figure 2. Effects of NaHS on oxidative stress and antioxidant system activity of the leaves of K. obovata seedlings treated by high salinity. (A) H2O2 content, (B) thiobarbituric acid reactive substances (TBARS) accumulation, (C) electrolyte leakage percentage (REL), (D) glutathione (GSH) content, (E) ascorbate acid peroxidase (APX) activity and (F) superoxide dismutase (SOD) activity. CK stands for the control treated only by a 1/4 strength Hoagland’s nutrient solution. NaCl stands for the 400 mM of the NaCl treatment. NaCl and H2S stands for the treatment with 400 mM of NaCl and 200 μM of NaHS. Error bars are SE (n = 3). The columns labeled with different letters are significantly different at p < 0.05 according to Duncan’s multiple range test.
Ijms 21 00118 g002
Figure 3. Two-dimensional (2-DE) analysis of proteins extracted from the leaves of K. obovata seedlings. The numbers assigned to the protein spots correspond to those listed in Table 1. (A) Representative coomassie brilliant blue-R250 (CBB-R250) stained 2D gel of total proteins extracted from K. obovata leaves under the control condition. Proteins (1.5 mg) were loaded onto the pH 4–7 gradient immobilized pH gradient (IPG) strip for the first-dimensional isoelectric focusing (IEF) and then separated in the second dimension on a 15% SDS-PAGE gel. The isoelectric point (pI) and molecular weight (MW) in kilodaltons are indicated on the top and left of the gel, respectively. Arrows indicate 37 spots with at least 2.0-fold changes (p < 0.05) that were analyzed by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI TOF MS). (B) The enlarged windows of the representative protein spots with different expressions under both the NaCl treatment and the NaCl and H2S treatment. NaCl stands for the 400 mM of NaCl treatment, NaCl and H2S stands for the treatment with 400 mM of NaCl and 200 μM NaHS, and CK stands for the control treated only by a 1/4 strength Hoagland solution.
Figure 3. Two-dimensional (2-DE) analysis of proteins extracted from the leaves of K. obovata seedlings. The numbers assigned to the protein spots correspond to those listed in Table 1. (A) Representative coomassie brilliant blue-R250 (CBB-R250) stained 2D gel of total proteins extracted from K. obovata leaves under the control condition. Proteins (1.5 mg) were loaded onto the pH 4–7 gradient immobilized pH gradient (IPG) strip for the first-dimensional isoelectric focusing (IEF) and then separated in the second dimension on a 15% SDS-PAGE gel. The isoelectric point (pI) and molecular weight (MW) in kilodaltons are indicated on the top and left of the gel, respectively. Arrows indicate 37 spots with at least 2.0-fold changes (p < 0.05) that were analyzed by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI TOF MS). (B) The enlarged windows of the representative protein spots with different expressions under both the NaCl treatment and the NaCl and H2S treatment. NaCl stands for the 400 mM of NaCl treatment, NaCl and H2S stands for the treatment with 400 mM of NaCl and 200 μM NaHS, and CK stands for the control treated only by a 1/4 strength Hoagland solution.
Ijms 21 00118 g003
Figure 4. (A) Functional classification and hierarchical clustering analysis. (B) Functional categories for the 37 differentially expressed proteins in K. obovata seedling leaves under the NaCl treatment and the NaCl and H2S treatment. The rows represent the individual proteins. The protein cluster is on the left side, and the treatment cluster is on the top. The up- or down-regulated proteins are indicated in red or green, and white represents no change. The intensity of the color increases with increasing expression differences, as shown in the bar at the bottom of the figure.
Figure 4. (A) Functional classification and hierarchical clustering analysis. (B) Functional categories for the 37 differentially expressed proteins in K. obovata seedling leaves under the NaCl treatment and the NaCl and H2S treatment. The rows represent the individual proteins. The protein cluster is on the left side, and the treatment cluster is on the top. The up- or down-regulated proteins are indicated in red or green, and white represents no change. The intensity of the color increases with increasing expression differences, as shown in the bar at the bottom of the figure.
Ijms 21 00118 g004
Figure 5. The two modules obtained from the protein–protein interaction network of differentially expressed proteins (DEPs). (A) Proteins in cluster 1 and (B) proteins in cluster 2.
Figure 5. The two modules obtained from the protein–protein interaction network of differentially expressed proteins (DEPs). (A) Proteins in cluster 1 and (B) proteins in cluster 2.
Ijms 21 00118 g005
Figure 6. Comparison of expression changes at (A) mRNA and (B) protein levels for the selected six DEPs. They are basic helix-loop-helix 145 (bHLH 145), superoxide dismutase (SOD), heat-shock protein (HSP), oxygen-evolving enhancer protein 1 (OEE1), abscisic acid stress ripening protein (Asr), and glutamine synthetase 1;1 (GS1;1). The columns labeled with different letters are significantly different at p < 0.05 according to a Student’s t-test.
Figure 6. Comparison of expression changes at (A) mRNA and (B) protein levels for the selected six DEPs. They are basic helix-loop-helix 145 (bHLH 145), superoxide dismutase (SOD), heat-shock protein (HSP), oxygen-evolving enhancer protein 1 (OEE1), abscisic acid stress ripening protein (Asr), and glutamine synthetase 1;1 (GS1;1). The columns labeled with different letters are significantly different at p < 0.05 according to a Student’s t-test.
Ijms 21 00118 g006
Figure 7. The proposed regulatory networks of H2S on the salt tolerance of K. obovata seedling leaves. The left and right dots stand for the NaCl and NaCl and H2S treatments, respectively. The red color of dot indicates the up-regulated change, the green color of dot indicates the down-regulated change, and the white color of dot indicates no change in comparison to the control. The left and right arrows stand for the NaCl treatment and the NaCl and H2S treatment, respectively. The red color arrow indicates physiological indexes of increased change, and the green color arrow indicates physiological indexes of decreased change. ADH: alcohol dehydrogenase; APX: cytosolic ascorbate peroxidase; Asr: abscisic acid stress ripening protein; atpC: ATP synthase CF1 epsilon subunit; ATP5F1B: ATP synthase beta subunit; ATP5PD: ATP synthase subunit d, mitochondrial; CP: cysteine proteinase; Cpn 20: 20 kDa chaperonin family protein; Cpn 60: TCP-1/cpn60 chaperonin family protein; CSD2: copper/zinc superoxide dismutase; CDSP: chloroplastic drought-induced stress protein; FLDH: NAD(P)-binding Rossmann-fold superfamily protein isoform 1; FNR: ferredoxin-NADP reductase; GS2: glutamine synthetase 2; GS1:1: glutamine synthetase 1;1; HSP: heat-shock protein; MAN: Glycosyl hydrolase superfamily protein; MDH: malate dehydrogenase; NDPK: nucleoside diphosphate kinase; OEE1: oxygen-evolving enhancer protein; OEC: O2 evolving complex 33kD family protein; PC: plastocyanin; TPI: triosephosphate isomerase; PCD: programmed cell death; PDX1: pyridoxin biosynthesis; PETC: photosynthetic electron transfer C; PGK: phosphoglycerate kinase; PGP: phosphoglycolate phosphatase; Prx: 2-Cys peroxiredoxin; RBCL: ribulose-1,5-bisphosphate carboxylase large subunit; RBCS: RuBisCO activase small isoform precursor; and SOD: superoxide dismutase. The detailed information for each spot is shown in Table 1.
Figure 7. The proposed regulatory networks of H2S on the salt tolerance of K. obovata seedling leaves. The left and right dots stand for the NaCl and NaCl and H2S treatments, respectively. The red color of dot indicates the up-regulated change, the green color of dot indicates the down-regulated change, and the white color of dot indicates no change in comparison to the control. The left and right arrows stand for the NaCl treatment and the NaCl and H2S treatment, respectively. The red color arrow indicates physiological indexes of increased change, and the green color arrow indicates physiological indexes of decreased change. ADH: alcohol dehydrogenase; APX: cytosolic ascorbate peroxidase; Asr: abscisic acid stress ripening protein; atpC: ATP synthase CF1 epsilon subunit; ATP5F1B: ATP synthase beta subunit; ATP5PD: ATP synthase subunit d, mitochondrial; CP: cysteine proteinase; Cpn 20: 20 kDa chaperonin family protein; Cpn 60: TCP-1/cpn60 chaperonin family protein; CSD2: copper/zinc superoxide dismutase; CDSP: chloroplastic drought-induced stress protein; FLDH: NAD(P)-binding Rossmann-fold superfamily protein isoform 1; FNR: ferredoxin-NADP reductase; GS2: glutamine synthetase 2; GS1:1: glutamine synthetase 1;1; HSP: heat-shock protein; MAN: Glycosyl hydrolase superfamily protein; MDH: malate dehydrogenase; NDPK: nucleoside diphosphate kinase; OEE1: oxygen-evolving enhancer protein; OEC: O2 evolving complex 33kD family protein; PC: plastocyanin; TPI: triosephosphate isomerase; PCD: programmed cell death; PDX1: pyridoxin biosynthesis; PETC: photosynthetic electron transfer C; PGK: phosphoglycerate kinase; PGP: phosphoglycolate phosphatase; Prx: 2-Cys peroxiredoxin; RBCL: ribulose-1,5-bisphosphate carboxylase large subunit; RBCS: RuBisCO activase small isoform precursor; and SOD: superoxide dismutase. The detailed information for each spot is shown in Table 1.
Ijms 21 00118 g007
Table 1. Identification of differentially expressed proteins with greater than 2.0-fold change in the leaves of K. obovata seedlings under the salinity and NaHS treatments.
Table 1. Identification of differentially expressed proteins with greater than 2.0-fold change in the leaves of K. obovata seedlings under the salinity and NaHS treatments.
Ratio 9
Spot 1Accession (gb) 2Protein Name 3Theoretical (Mr/pI) 4Observed (Mr/pI) 5Score 6MP 7Species 8NaCl vs. CKNaCl and H2S vs. CK
Photosynthesis
1gi|502131189PREDICTED: plastocyanin, chloroplastic-like17.13/5.0414.05/4.371181Cicer arietinum2.257 ± 0.1930.841 ± 0.068
2gi|226497434Ferredoxin-NADP reductase, leaf isozyme40.98/8.5337.32/6.743114Zea mays0.287 ± 0.0550.686 ± 0.153
8gi|363981020Ribulose-1,5-bisphosphate carboxylase/oxygenase large subunit, partial (chloroplast)50.67/5.8719.87/5.191167Cercidiphyllum japonicum0.086 ± 0.0140.736 ± 0.084
10gi|508707371Photosynthetic electron transfer C24.67/8.1519.14/4.791035Theobroma cacao0.100 ± 0.0050.483 ± 0.068
12gi|222853091O2 evolving complex 33kD family protein35.41/5.8536.80/5.224054Populus trichocarpa0.067 ± 0.0060.518 ± 0.065
14gi|527190719Oxygen-evolving enhancer protein 1, chloroplastic34.80/6.4834.78/5.474938Genlisea aurea1.213 ± 0.0610.068 ± 0.012
36gi|62733297RuBisCO activase small isoform precursor52.39/5.5942.29/5.032767Oryza sativa Japonica Group0.741 ± 0.1570.079 ± 0.033
Carbohydrate and energy metabolism
6gi|508785499Glycosyl hydrolase superfamily protein, putative37.45/8.3934.33/4.66681Theobroma cacao2.072 ± 0.4311.279 ± 0.188
7gi|336041766ATP synthase CF1 epsilon subunit13.69/5.4314.89/5.03972Justicia americana0.323 ± 0.0980.528 ± 0.060
9gi|460411739PREDICTED: ATP synthase subunit d, mitochondrial-like19.78/5.3318.95/5.041194Solanum lycopersicum0.704 ± 0.0593.630 ± 0.341
16gi|508786769Triosephosphate isomerase27.50/5.5423.96/6.351647Theobroma cacao0.142 ± 0.0380.605 ± 0.061
20gi|223535342Alcohol dehydrogenase, putative41.61/8.6137.23/5.885136Ricinus communis0.643 ± 0.1580.311 ± 0.064
22gi|226503019Malate dehydrogenase, cytoplasmic35.85/5.7640.48/6.782507Zea mays3.659 ± 0.5202.997 ± 0.469
24gi|508715598NAD(P)-binding Rossmann-fold superfamily protein isoform 1 36.50/9.2928.24/5.532684Theobroma cacao0.663 ± 0.1310.469 ± 0.079
25gi|502105712Phosphoglycolate phosphatase40.78/6.8937.32/6.7422610Cicer arietinum2.400 ± 0.3731.927 ± 0.395
29gi|335059237ATP synthase beta subunit52.29/5.0957.59/4.9745311Clerodendrum trichotomum1.591 ± 0.0622.706 ± 0.178
30gi|475549973Nucleoside diphosphate kinase 117.03/6.8516.16/6.401075Aegilops tauschii0.419 ± 0.0409.366 ± 1.258
35gi|223547261Phosphoglycerate kinase, putative50.11/8.7451.05/5.992946Ricinus communis1.115 ± 0.2563.040 ± 0.592
Amino acid and protein metabolism
4gi|355492134Cysteine proteinase53.52/8.1240.03/4.48965Medicago truncatula1.512 ± 0.2060.162 ± 0.012
19gi|223544592Heat-shock protein, putative17.81/5.9316.81/5.691296Ricinus communis0.734 ± 0.0712.172 ± 0.799
26gi|55033629220 kDa chaperonin family protein26.89/8.7523.19/5.622455Populus trichocarpa0.196 ± 0.0922.934 ± 0.618
28gi|508779629TCP-1/cpn60 chaperonin family protein65.05/5.5756.09/5.2771415Theobroma cacao2.301 ± 0.3191.562 ± 0.134
32gi|508713595Glutamine synthetase 261.64/8.3844.25/5.98633Theobroma cacao0.208 ± 0.0321.068 ± 0.160
33gi|332006826Glutamine synthetase 1;139.32/5.2840.86/5.712417Arabidopsis thaliana1.294 ± 0.2801.807 ± 0.246
Hormone biosynthesis and transcription factor
13gi|332008500Transcription factor bHLH14535.26/5.0840.10/5.526411Arabidopsis thaliana4.254 ± 0.4981.032 ± 0.079
21gi|508712975Mitochondrial transcription termination factor family protein, putative isoform 151.01/9.3244.71/6.385611Theobroma cacao0.252 ± 0.0320.393 ± 0.063
27gi|355514936BZIP transcription factor24.12/8.6523.37/5.88609Medicago truncatula2.478 ± 0.1583.198 ± 0.157
31gi|355481146Abscisic acid stress ripening protein27.33/5.1539.71/5.92072Medicago truncatula2.024 ± 0.2071.567 ± 0.182
Stress response proteins
5gi|502112102PREDICTED: 2-Cys peroxiredoxin BAS1-like, chloroplastic-like isoform X229.14/6.1221.46/4.963616Cucumis sativus2.446 ± 0.6162.257 ± 0.502
11gi|222867611Pyridoxin biosynthesis PDX1-like protein 333.31/6.5515.98/5.29608Populus trichocarpa0.886 ± 0.2662.230 ± 0.521
17gi|409900374Copper/zinc superoxide dismutase CSD2A-123.29/6.1219.07/5.743934Musa acuminata1.184 ± 0.1672.006 ± 0.334
18gi|414866828Superoxide dismutase 415.65/5.1015.02/5.80773Zea mays2.510 ± 0.6880.988 ± 0.100
23gi|508776520Actin 7 isoform 141.80/5.3143.65/5.523429Theobroma cacao0.827 ± 0.2103.400 ± 0.725
34gi|508723241Chloroplastic drought-induced stress protein of 32 kD40.91/7.6627.45/5.772136Theobroma cacao0.237 ± 0.0780.210 ± 0.058
37gi|498923199Cytosolic ascorbate peroxidase27.09/5.5232.52/6.662343Arachis hypogaea0.734 ± 0.2352.172 ± 0.334
Unknown protein
3gi|548848586Hypothetical protein AMTR_s00157p0006468054.99/6.9128.30/4.86681Amborella trichopoda0.120 ± 0.0110.759 ± 0.362
15gi|462410037Hypothetical protein PRUPE_ppa019045mg25.73/5.6529.45/6.011505Prunus persica0.488 ± 0.1090.669 ± 0.116
1 The spot number corresponds to the number listed in the table. 2 Database accession numbers (gb) according to NCBInr. 3 The names of proteins were identified by LC-MALDI-TOF/TOF. 4 Theoretical mass (kDa) and pI of identified proteins. Theoretical values were retrieved from the NCBInr database. 5 Experimental mass (kDa) and pI of identified proteins. Experimental values were calculated by using the PDQuest software (Version 8.0, Bio-Rad, Hercules, CA, USA) and standard molecular mass. 6 The Mascot searched score against the database NCBInr. 7 Number of matched peptide fragments. 8 The species which has the high homology of the identified protein. 9 Ratio between the different treatments. NaCl vs. CK means 400 mM of NaCl vs. control; NaCl and H2S vs. CK means 400 mM of NaCl and 200 μM of NaHS vs. control.

Share and Cite

MDPI and ACS Style

Liu, Y.-L.; Shen, Z.-J.; Simon, M.; Li, H.; Ma, D.-N.; Zhu, X.-Y.; Zheng, H.-L. Comparative Proteomic Analysis Reveals the Regulatory Effects of H2S on Salt Tolerance of Mangrove Plant Kandelia obovata. Int. J. Mol. Sci. 2020, 21, 118. https://doi.org/10.3390/ijms21010118

AMA Style

Liu Y-L, Shen Z-J, Simon M, Li H, Ma D-N, Zhu X-Y, Zheng H-L. Comparative Proteomic Analysis Reveals the Regulatory Effects of H2S on Salt Tolerance of Mangrove Plant Kandelia obovata. International Journal of Molecular Sciences. 2020; 21(1):118. https://doi.org/10.3390/ijms21010118

Chicago/Turabian Style

Liu, Yi-Ling, Zhi-Jun Shen, Martin Simon, Huan Li, Dong-Na Ma, Xue-Yi Zhu, and Hai-Lei Zheng. 2020. "Comparative Proteomic Analysis Reveals the Regulatory Effects of H2S on Salt Tolerance of Mangrove Plant Kandelia obovata" International Journal of Molecular Sciences 21, no. 1: 118. https://doi.org/10.3390/ijms21010118

APA Style

Liu, Y. -L., Shen, Z. -J., Simon, M., Li, H., Ma, D. -N., Zhu, X. -Y., & Zheng, H. -L. (2020). Comparative Proteomic Analysis Reveals the Regulatory Effects of H2S on Salt Tolerance of Mangrove Plant Kandelia obovata. International Journal of Molecular Sciences, 21(1), 118. https://doi.org/10.3390/ijms21010118

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop