- Open Access
Proteome profiling reveals insights into cold-tolerant growth in sea buckthorn
Proteome Science volume 14, Article number: 14 (2016)
Low temperature is one of the crucial environmental factors limiting the productivity and distribution of plants. Sea buckthorn (Hippophae rhamnoides L.), a well recognized multipurpose plant species, live successfully in in cold desert regions. But their molecular mechanisms underlying cold tolerance are not well understood.
Physiological and biochemical responses to low-temperature stress were studied in seedlings of sea buckthorn. Differentially expressed protein spots were analyzed using multiplexing fluorescent two-dimensional fluorescence difference gel electrophoresis (2D-DIGE) coupled with matrix-assisted laser desorption/ionization (MALDI) time-of-flight/time-of-flight (TOF/TOF) mass spectrometry (MS), the concentration of phytohormone was measured using enzyme-linked immunosorbent assay, and a spectrophotometric assay was used to measure enzymatic reactions.
With the increase of cold stress intensity, the photosynthesis rate, transpiration rate, stomatal conductance in leaves and contents of abscisic acid (ABA) and indole acetic acid (IAA) in roots decreased significantly; however, water-use efficiency, ABA and zeatin riboside in leaves increased significantly, while cell membrane permeability, malondialdehyde and IAA in leaves increased at 7 d and then decreased at 14 d. DIGE and MS/MS analysis identified 32 of 39 differentially expressed protein spots under low-temperature stress, and their functions were mainly involved in metabolism, photosynthesis, signal transduction, antioxidative systems and post-translational modification.
The changed protein abundance and corresponding physiological–biochemical response shed light on the molecular mechanisms related to cold tolerance in cold-tolerant plants and provide key candidate proteins for genetic improvement of plants.
Plants being sessile have developed intricate adaptation mechanisms to stresses. Among abiotic environmental stresses, low temperature is one of the crucial factors limiting plant productivity and distribution [1–3]. Therefore, identifying the reasons behind cold adaptation in cold-tolerant plants could be of central importance to gaining a deeper understanding of the mechanisms involved, and be beneficial for plant growth and production [4, 5].
Sea buckthorn (Hippophae rhamnoides) is well adapted to extreme temperatures ranging from −43 to 40 °C , and it is distributed in cold regions of Russia, China, Finland, Sweden and many other countries of Asia, Europe and North America , which provided an ideal material to study mechanisms related to low-temperature (LT) stress tolerance . Although there is a considerable wealth of studies on origin, distribution, ecology and nutrient and medicinal values in Hippophae [8–15], the physiological and biochemical basis of tolerance and mechanisms of abiotic stress response, especially low-temperature or cold response, are not well understood [2, 16, 17].
Due to the direct roles of proteins in plant stress responses, profound changes in proteome composition can be observed during plant acclimation to stress. Mass spectrometry (MS)-based proteomics has become an essential tool in unraveling possible relationships between protein abundance and plant stress acclimation [18, 19]. The present study discusses the proteome-wide protein responses to low-temperature stress of H. rhamnoides cv. ‘Chuisk’, a widely cultivated hybrid of H. rhamnoides subsp. mongolica and H. rhamnoides subsp. rhamnoides in northeast China, which has excellent cold resistance and good characteristics, including large berries, high content of oil and high production. Using physiological, biochemical and comparative proteomic analyses, we hope to provide insights into cold adaptation mechanisms in this cold-tolerant species.
Physiological and biochemical responses during LT stress
Under the LT treatments 7 d (T1) and 14 d (T2), the values of superoxide dismutase (SOD), glutathione reductase (GR) and zeatin riboside in root [ZR(R)] showed no significant changes, but leaf area (LA) and gibberellins (GA3) did. The decreases in stomatal conductance (Cond) and abscisic acid in root [ABA(R)] were significant with the extension of LT treatment time (p-value ≤ 0.05). There were significant increases in water-use efficiency (WUE), ZR (L) (L indicated leaf) and ABA (L) for T1. Gradual decreases in photosynthesis rate (Pn), transpiration rate (Tr) and indole acetic acid [IAA(R)] were observed among both controls and treatments. Cell membrane permeability (CMP), malondialdehyde (MDA) and IAA(L) showed non-consistent changes (increase at T1 and decrease at T2) under LT stress (Fig. 1).
Variations in LT responsive proteins
The DeCyder image analysis of fluorescent images detected 1466 ± 35 protein spots (Fig. 2, Additional file 1: Figure S1). One-way analysis of variance (ANOVA) showed that 39 different protein spots were significantly affected by LT stress (p-value ≤ 0.05) (Fig. 2), of which 32 proteins were identified using the MS/MS method: eight exhibited a gradual increase (A-type), five exhibited a gradual decrease (B-type), 11 exhibited an increase at T1 and decrease at T2 (C-type) and eight exhibited a decrease at T1 and increase at T2 (D-type) (Additional file 2: Table S1).
Gene ontology (GO) and pathway enrichment analysis of differentially expressed protein spots
Ninety-seven GO enrichment terms were obtained using GOEAST . The enriched biological process ontology included the metabolic process (nitrate assimilation and nucleotide metabolic, oxidoreduction coenzyme metabolic and carotenoid biosynthetic processes), biological regulation (regulation of protein dephosphorylation), response to stimulus (responses to cold, glucose and fructose) and localization (mitochondrial transport) (Fig. 3 and Additional file 2: Table S2).
KEGG pathway analysis  showed that 18 terms including genetic information processing (folding, sorting and degradation, and translation) and metabolism (amino acid, carbohydrate and energy metabolisms) were enriched under LT stress (Additional file 2: Table S3). Of the 32 identified protein spots, 14 (43.75 %) proteins were classified to metabolism pathways including energy metabolism (14 protein spots), carbohydrate metabolism (12 spots) and amino acid metabolism (spots 486 and 693, C-type; and spot 437, D-type). Energy metabolism was enriched in carbon fixation in photosynthetic organisms including ribulose-bisphosphate carboxylase large chain (rbcL; 10 homologous/isoform protein spots), photosynthesis involving oxygen-evolving enhancer protein 1 (spot 770, A-type), photosystem I reaction center subunit II (spot 966, C-type) and nitrogen metabolism (spot 486, C-type, glutamine synthetase, glnA). Carbohydrate metabolism was enriched in glyoxylate and dicarboxylate metabolism (11 protein spots), and citrate cycle and pyruvate metabolism (spot 693, C-type, Nodule-enhanced malate dehydrogenase). Four of 32 proteins belonged to the pathway of genetic information processing (spots 298 and 345, A-type, TCP-1/cpn60 chaperonin family protein; spot 1117, D-type, 18.5-kDa class I heat shock protein; and spot 722, D-type, unnamed protein) (Additional file 2: Tables S1 and S3).
mRNA expression of differentially expressed protein spots between different treatments
To confirm the reliability of the MALDI-TOF MS/MS and 2D-DIGE in detecting H. rhamnoides, we used RealTime PCR to compare the expression levels of identified protein genes. Ten genes were randomly selected from the differentially expressed protein spots in the three groups. Similar expression patterns between protein and mRNA were observed in six of the ten detected genes (Fig. 4). Of which, two were up-regulated (Spots 298 and 770), one was down-regulated (Spot 1285), one was up-regulated under D1 stress and down-regulated under D2 stress (Spot 966), and two was down-regulated under D1 stress and up-regulated under D2 stress (Spots 1094 and 1117).
Endogenous hormones and antioxidant enzymes play roles in physiological and biochemical responses of sea buckthorn under LT stress
Compared with GA3 and ZR, ABA and IAA responded more actively. The LT stress induced the decline of ABA content in sea buckthorn roots, and then accumulation in leaves, which led to stomatal closure and restrained aboveground Pn, Tr and Cond. Meanwhile, the complementary/equilibrium of growth-promoting hormones such as ZR and IAA possibly contributed to maintaining plant growth, while the increase of ABA in leaves also showed growth inhibition under LT stress . The changed protein abundance of probable LRR receptor-like serine/threonine-protein kinase RKF3-like (spot 1426, C-type) involved in hormone-mediated (ABA and GA) signaling pathways, and hypothetical protein (spot 830, C-type) participating in IAA biosynthetic/metabolic process, further showed that ABA, GA and IAA played an important role in responding to LT stress in sea buckthorn, as found in many other cold-stress studies [23–27]. These confirmed that many hormonal proteins contributed to the changes in ABA, GA and IAA content under LT stress and so induced the hormone signaling pathways responding to LT stress in sea buckthorn.
In addition, the CMP and MDA increase under LT stress suggested that they induced cell damage and membrane lipid peroxidation. Some protein spots of antioxidant enzymes, such as zinc-binding dehydrogenase family oxidoreductase (spot 544, A-type), peroxiredoxin family protein (Prx, spot 888, D-type), dehydroascorbate reductase (DHAR, spot 900, C-type) and glnA (spot 486, C-type), changed dynamically in protein abundance under LT stress, suggesting that they had taken part in changes of enzymatic activity or antioxidants (i.e. SOD and GR). Many proteins involved in antioxidative systems were also induced by LT stress in H. rhamnoides  and Arabidopsis . This study hint at cooperation between hormone and redox signaling pathways in responding to LT stress, and this confers enormous regulatory potential because it allows plants to adapt to changing and often challenging stress conditions .
Proteins related to metabolism, signal transduction and post-translational modification contribute to LT stress responses in sea buckthorn
In forest tree species, many genes taking part in metabolic processes were differentially expressed under most abiotic stress conditions [30–34]. Pn significantly decreased under LT stress, and expression of the rbcLs and oxygen-evolving enhancer proteins in sea buckthorn greatly changed consistent with a number of previous reports [35–37]. This suggested negative effects on photosynthetic efficiency and photosynthesis rate. Many studies have shown that LT triggers a cellular signal transduction pathway leading to molecular and metabolic changes [26, 38–40]. This study found glutamine synthetase (spot 486, C-type) was differently expressed during LT stress, suggesting that cold stress stimuli initiated the signal transmission processes in H. rhamnoides. During cold acclimation, plants accumulate specific types of solutes, such as soluble sugars (i.e. glucose and fructose)  and carotenoids [42, 43] to enhance freezing tolerance. The different expression of Prx (spot 888, D-type) and glnA (spot 486, C-type) under LT stress suggested the accumulation of soluble sugars in H. rhamnoides. Meanwhile, two proteins (spots 544 and 594, A-type) involving carotenoid biosynthetic process were enriched in the up-regulated differently expressed protein spots, suggesting an important role of carotenoid in responding to LT stress in sea buckthorn.
It is noteworthy that two identified isoforms had multiple spots located at different positions on the gel with different molecular weight, pI or both. For example, spots 298 and 345 were both identified as TCP-1/cpn60 chaperonin family protein and spots 369, 550, 764, 811, 865, 1012, 1285, 1290, 1316 and 1326 were rbcL homologs. These isoforms might represent post-translationally modified forms of the same protein, such as phosphorylated and glycosylated forms under LT stress, hinting at the post-translational regulation of cold acclimation responsein H. rhamnoides [44, 45].
According to the findings, we proposed some scenarios involved in LT stress tolerance of H. rhamnoides. During the growth season, LT abiotic stress triggered signal transduction, stimulated genetic information processing (such as protein degradation and synthesis, and changes in enzymatic activity and hormone content), promoted cellular processes (cell damage, membrane lipid peroxidation and stomatal opening/closing) and regulated necessary metabolism pathways (especially the energy and carbohydrate metabolism mostly from photosynthesis). This then induced physiological (such as Pn, Tr, Cond, CMP and WUE) and biochemical changes (such as activity of SOD and GR, and contents of MDA and endogenous hormones), eventually leading to growth retardation of H. rhamnoides (Fig. 5). Additionally, post-translational modification (PTM), such as dephosphorylation (Fig. 3) might also play important roles in cold stress responses and adaptation in sea buckthorn .
The responses of sea buckthorn to cold stress were complex and involved numerous physiological, molecular and cellular adaptations. The proteins involved in photosynthesis, plant hormone signal transduction, antioxidative systems and PTM played vital roles. Furthermore, the 2D-PAGE profile of sea buckthorn under LT stress provided useful candidate proteins for genetic improvement and understanding the general ability of plants to respond to a wide range of external environmental stresses.
Plant materials and treatments
Root turions of H. rhamnoides cv. ‘Chuisk’ were planted in 18 cm × 22 cm × 26 cm plastic pots containing 6 kg of a mixture of turfy soil/clay soil/sand (5:3:2). The plants were grown in a greenhouse under 25/20 °C of day/night temperatures, relative humidity of 60–70 % and CO2 concentration of 380 ± 10 μmol mol−1. After three months, the plants were divided into three groups/treatments (CK, T1 and T2), each with three biological replicates, applied as follows. CK was an unstressed control, T1 and T2 were 4 °C LT treatment, in which soil water content was maintained at about 30 % (i.e. 75–85 % of field moisture capacity) by regular watering during the whole experiment (i.e. from 0 to 14 d). Samples taken at 7 (T1) and 14 d (T2) were snap-frozen in liquid nitrogen, and stored at −80 °C for further studies. Soil water content was checked every 2 d using a moisture probe type HH2 Meter (Delta-T Devices Ltd., Cambridge, UK). For this, the plants were cultured separately in a 3.5 m × 2.2 m × 3.2 m AGC-2 Growth Cabinet (Zhejiang University Electric Equipment Factory, Hangzhou, China), whose parameters could be regulated strictly to be consistent with the greenhouse conditions.
Measurement of physiological and biochemical parameters
Three fully expanded leaves in the middle portion of each stem were selected for measurement with a LI-6400 Portable Photosynthesis System (LI-COR Inc., Lincoln, NE, USA). Every measurement was repeated five times. The changes in Pn (μmol · m−2 · s−1), Tr (μmol · m−2 · s−1) and Cond (mol · m−2 · s−1) are shown in Fig. 1. The activities of major antioxidant enzymes (SOD and GR) (U · g−1), MDA (μmol · g−1) content were analyzed according to the method of Gao . Molar conductances, used for CMP (%) were measured using a digital conductivity meter DDS-12DW Microprocessor (LIDA Instrument Factory, Shanghai, China). Fresh LA (cm2) was measured using a leaf area meter (CI-202, CID Inc., WA, USA). The content (ng · g−1) of four kinds of phytohormones (ABA, IAA, GA3 and ZR) were respectively measured in leaves and roots using enzyme-linked immunosorbent assay method. Each measure had three biological replicates.
Protein extraction and quantification
Of frozen leaves, 2 g was homogenized in a pre-chilled mortar in 25 ml of cold acetone (−20 °C) containing 10 % (V/V) trichloroacetic acid and 65 mM dithiothreitol. The mixture was kept at −20 °C for 1 h, and then centrifuged at 12,000 g for 45 min. The supernatant was discarded and the pellet precipitated with 25 ml of cold acetone. After 1 h at −20 °C, the mixture was centrifuged at 12,000 g for 45 min, and the white pellet was lyophilized and stored at −80 °C. The protein pellet was dissolved in DIGE lysis buffer (7 M urea, 2 M thiourea, 4 % 3-[(3-chol-amidopropyl) dimethylammonio]-1-propane sulfonate and 0.2 % IPG buffer). The mixture was shaken and centrifuged at 12,000 g for 1 h. The supernatant was collected and protein concentrations were determined using the Bradford method with the protein assay reagent (Bio-Rad Laboratories, Hercules, CA, USA).
Two-dimensional difference gel electrophoresis and imaging analysis
Nine protein samples (50 μg, pH 8.0–9.0) labeled with Cy3 and Cy5 dyes (labeled with Cy2 dye for use as internal standard for normalization) were mixed with a strip rehydration and then separated by 2D-polyacrylamide gel electrophoresis (PAGE). Total volume was adjusted to 50 μL, mixed 1:1 with DIGE lysis buffer and incubated for 30 min on ice in the dark. The reaction was then quenched with 10 mM lysine and additionally incubated for 10 min [48, 49]. Isoelectrofocusing was carried out with the IPGphor system (Amersham Pharmacia Biotech, Uppsala, Sweden) using Immobiline DryStrip gels (13 cm) with nonlinear pH gradients (pH 3–10). The IPGphor system was then programmed as follows: 12 h at 30 V, 1 h at 500 V, 1 h at 1000 V, 8 h at 8000 V and 4 h at 500 V. After IEF, the gel strips were incubated with equilibration buffer, and the strips were then transferred onto vertical 12.5 % (w/v) sodium dodecyl sulfate (SDS)-PAGE gels. The SDS-PAGE was performed on a Hoefer SE600 system (Amersham Pharmacia Biotech) with 15 mA/gel for 30 min and then 30 mA/gel until the bromophenol blue dye reached the bottom (about 0.5 cm) of the gels. The gels were scanned using an UMax Powerlook 2110XL laser scanner (Amersham Pharmacia Biotech), and the images analyzed using DeCyderTM 2D software (Amersham Pharmacia Biotech). Intra-gel spot detection, quantification and inter-gel matching and quantification were performed using Differential In-gel Analysis and Biological Variation Analysis modules respectively. To test for significant differences in expression of protein spots among pairwise groups, one-way analysis of variance (ANOVA) was performed at p-value ≤ 0.05. The differentially expressed protein spots were filtered based on an average volume ratio of twofold with significance at p-value ≤ 0.05. Each treatment had three biological replicates.
In-gel trypsin digestion and MS analysis
These gels were fixed and stained with colloidal Coomassie brilliant blue (CBB). Proteins of interest, as defined by the DeCyder analysis TM 2D software, were excised from the CBB-stained gels for an in-gel trypsin digestion procedure . The resolved peptides were analyzed on matrix-assisted laser desorption/ionization (MALDI) time-of-flight/time-of-flight (TOF/TOF) based on an MALDI TOF/TOF system (ultrafleXtreme, Bruker Daltonics, MA, USA). Protein identification was performed using the Mascot search engine version 2.1.03 (Matrix Science, London, UK) with a mass tolerance of 0.05 Da permitted for intact peptide masses and fragmented ions, and allowance for one missed cleavage in trypsin digests. Gln- > pyro-Glu (N-term Q), Oxidation (M) and Deamidated (NQ) as the potential variable modifications, Carbamidomethyl (C) as fixed modifications. The search that was performed using MASCOT v2.1.03 software (Matrix Science, London, UK) and the following settings: Uniprot_plant_database, one missed cleavage, fixed modifications of carbamidomethyl, variable modifications of oxidation, peptide tolerance 100 ppm, fragment mass tolerance 0.5 Da, peptide charge 1+. Only peptides with MS/MS ion scores significantly (P < 0.05) exceeding the MASCOT identity or extensive homology threshold were reported. Proteins were annotated using Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. GO and pathway enrichment were analyzed using Gene Ontology Enrichment Analysis Software Toolkit (GOEAST)  and KEGG Orthology Based Annotation System (KOBAS)  Software Toolkit, respectively.
Relative quantification of mRNA expression by RealTime PCR
Total RNAs from H. rhamnoides leaves in three groups were used for quantitative PCR analysis. Briefly, the first cDNA strands were obtained using a Thermo First cDNA Synthesis Kit, and were then subjected to quantification with 18S rRNA as an internal control using a standard SYBR Green PCR kit (Bio-Rad, Hercules, CA) on the StepOnePLUS RealTime PCR Detection System (Applied Biosystems, Foster City, CA). Quantitative PCR was then performed using the following conditions: 95 °C for 10 min; 40 cycles of 95 °C for 20 s and a 60 °C annealing temperature for 30 s; 95 °C for 15 s; 60 °C for 30 s; and 95 °C for 15 s. The primers for all 10 genes are listed in Online Additional file 2: Table S4. All reactions were performed in triplicate for each sample. Gene expression was quantified relative to 18S rRNA expression using the comparative cycle threshold (ΔCT) method. Differences in gene expression between the CK and LT stresses (T1 and T2) were detected by using the t-test.
Theocharis A, Clément C, Barka EA. Physiological and molecular changes in plants grown at low temperatures. Planta. 2012;235:1091–105.
Gupta R, Deswal R. Low temperature stress modulated secretome analysis and purification of antifreeze protein from Hippophae rhamnoides, a Himalayan wonder plant. J Proteome Res. 2012;11:2684–96.
Chinnusamy V, Zhu J, Zhu J-K. Cold stress regulation of gene expression in plants. Trends Plant Sci. 2007;12:444–51.
Lawlor DW. Genetic engineering to improve plant performance under drought: physiological evaluation of achievements, limitations, and possibilities. J Exp Bot. 2013;64:83–108.
Mickelbart MV, Hasegawa PM, Bailey-Serres J. Genetic mechanisms of abiotic stress tolerance that translate to crop yield stability. Nature Rev Genet. 2015;16:237–51.
Ruan CJ, Teixeira da Silva J, Jin H, Li H, Li DQ. Research and biotechnology in sea buckthorn (Hippophae spp.). Med Aromat Plant Sci Biotechnol. 2007;1:47–60.
Chaudhary S, Sharma PC. DeepSAGE based differential gene expression analysis under cold and freeze stress in seabuckthorn (Hippophae rhamnoides L.). PLoS One. 2015;10:e0121982.
Arimboor R, Arumughan C. HPLC-DAD-MS/MS profiling of antioxidant flavonoid glycosides in sea buckthorn (Hippophae rhamnoides L.) seeds. Int J Food Sci Nutr. 2012;63:730–8.
Michel T, Destandau E, Floch GL, Lucchesi ME, Elfakir C. Antimicrobial, antioxidant and phytochemical investigations of sea buckthorn (Hippophaë rhamnoides L.) leaf, stem, root and seed. Planta Med. 2012;131:754–60.
Dalija S, Daina K, Silvija R, Inta K. The effect of processing on the composition of sea buckthorn juice. J Fruit Ornam Plant Res. 2006;14:257-64.
Fatima T, Snyder CL, Schroeder WR, Cram D, Datla R, Wishart D, Weselake RJ, Krishna P. Fatty acid composition of developing sea buckthorn (Hippophae rhamnoides L.) berry and the transcriptome of the mature seed. PLoS One. 2012;7:e34099.
Bartish, Jeppsson, Nybom. Population genetic structure in the dioecious pioneer plant species Hippophae rhamnoides investigated by RAPD markers. Mol Ecol. 1999;8:791–802.
Sun K, Chen X, Ma R, Li C, Wang Q, Ge S. Molecular phylogenetics of Hippophae L. (Elaeagnaceae) based on the internal transcribed spacer (ITS) sequences of nrDNA. Plant Syst Evol. 2002;235:121–34.
Kallio H, Yang B, Peippo P. Effects of different origins and harvesting time on vitamin C, tocopherols, and tocotrienols in sea buckthorn (Hippophae rhamnoides) berries. J Agric Food Chem. 2002;50:6136–42.
Jia DR, Abbott RJ, Liu TL, Mao KS, Bartish IV, Liu JQ. Out of the Qinghai–Tibet Plateau: evidence for the origin and dispersal of Eurasian temperate plants from a phylogeographic study of Hippophaë rhamnoides (Elaeagnaceae). New Phytologist. 2012;194:1123–33.
Ghangal R, Raghuvanshi S, Sharma PC. Expressed sequence tag based identification and expression analysis of some cold inducible elements in seabuckthorn (Hippophae rhamnoides L.). Plant Physiol Biochem. 2012;51:123–8.
Li C, Yang Y, Junttila O, Palva ET. Sexual differences in cold acclimation and freezing tolerance development in sea buckthorn (Hippophae rhamnoides L.) ecotypes. Plant Sci. 2005;168:1365–70.
Xu G, Li C, Yao Y. Proteomics analysis of drought stress-responsive proteins in Hippophae rhamnoides L. Plant Mol Biol Report. 2009;27:153–61.
Zhang M, Lv D, Ge P, Bian Y, Chen G, Zhu G, Li X, Yan Y. Phosphoproteome analysis reveals new drought response and defense mechanisms of seedling leaves in bread wheat (Triticum aestivum L.). J Proteomics. 2014;109:290–308.
Jin ZP, Wen XF, Zhang JK, Gu YK, Jin SY. Analysis for affect factors of drought stress on production of different kinds of Seabuckthorn sources. The Global Seabuckthorn Research and Development. 2008;2:33–8.
Xie C, Mao X, Huang J, Ding Y, Wu J, Dong S, Kong L, Gao G, Li C-Y, Wei L. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011;39:W316–22.
Huang X, Chen MH, Yang LT, Li YR, Wu JM. Effects of exogenous abscisic acid on cell membrane and endogenous hormone contents in leaves of sugarcane seedlings under cold stress. Sugar Tech. 2015;17:59–64.
Sreenivasulu N, Harshavardhan VT, Govind G, Seiler C, Kohli A. Contrapuntal role of ABA: does it mediate stress tolerance or plant growth retardation under long-term drought stress? Gene. 2012;506:265–73.
Miao ZY, Xu W, Li DF, Hu XN, Liu JX, Zhang RX, Tong ZY, Dong JL, Su Z, Zhang LW. De novo transcriptome analysis of Medicago falcata reveals novel insights about the mechanisms underlying abiotic stress-responsive pathway. BMC Genomics. 2015;16:818.
Maruyama K, Urano K, Yoshiwara K, Morishita Y, Sakurai N, Suzuki H, Kojima M, Sakakibara H, Shibata D, Saito K. Integrated analysis of the effects of cold and dehydration on rice metabolites, phytohormones, and gene transcripts. Plant Physiol. 2014;164:1759–71.
Colebrook EH, Thomas SG, Phillips AL, Hedden P. The role of gibberellin signalling in plant responses to abiotic stress. J Exp Biol. 2014;217:67–75.
Shibasaki K, Uemura M, Tsurumi S, Rahman A. Auxin response in Arabidopsis under cold stress: underlying molecular mechanisms. Plant Cell. 2009;21:3823–38.
Goulas E, Schubert M, Kieselbach T, Kleczkowski LA, Gardeström P, Schröder W, Hurry V. The chloroplast lumen and stromal proteomes of Arabidopsis thaliana show differential sensitivity to short‐and long‐term exposure to low temperature. Plant J. 2006;47:720–34.
Bartoli CG, Casalongué CA, Simontacchi M, Marquez-Garcia B, Foyer CH. Interactions between hormone and redox signalling pathways in the control of growth and cross tolerance to stress. Environ Exp Bot. 2013;94:73–88.
Kieffer P, Schröder P, Dommes J, Hoffmann L, Renaut J, Hausman J-F. Proteomic and enzymatic response of poplar to cadmium stress. J Proteomics. 2009;72:379–96.
He CY, Zhang JG, Duan AG, Zheng SX, Sun HG, Fu LH. Proteins responding to drought and high-temperature stress in Populus × euramericana cv. ‘74/76’. Trees. 2008;22:803–13.
Ferreira S, Hjernø K, Larsen M, Wingsle G, Larsen P, Fey S, Roepstorff P, Pais MS. Proteome profiling of Populus euphratica Oliv. upon heat stress. Ann Bot. 2006;98:361–77.
Costa P, Bahrman N, Frigerio J-M, Kremer A, Plomion C. Water-deficit-responsive proteins in maritime pine. Plant Mol Biol. 1998;38:587–96.
Echevarría-Zomeño S, Ariza D, Jorge I, Lenz C, Del Campo A, Jorrín JV, Navarro RM. Changes in the protein profile of Quercus ilex leaves in response to drought stress and recovery. J Plant Physiol. 2009;166:233–45.
Evers D, Legay S, Lamoureux D, Hausman J, Hoffmann L, Renaut J. Towards a synthetic view of potato cold and salt stress response by transcriptomic and proteomic analyses. Plant Mol Biol. 2012;78:503–14.
Gharechahi J, Alizadeh H, Naghavi MR, Sharifi G. A proteomic analysis to identify cold acclimation associated proteins in wild wheat (Triticum urartu L.). Mol Biol Rep. 2014;41:3897–905.
Rocco M, Arena S, Renzone G, Scippa GS, Lomaglio T, Verrillo F, Scaloni A, Marra M. Proteomic analysis of temperature stress-responsive proteins in Arabidopsis thaliana rosette leaves. Mol Biosyst. 2013;9:1257–67.
Wang XC, Zhao QY, Ma CL, Zhang ZH, Cao HL, Kong YM, Yue C, Hao XY, Chen L, Ma JQ. Global transcriptome profiles of Camellia sinensis during cold acclimation. BMC Genomics. 2013;14:415.
Yang QS, Gao J, He WD, Dou TX, Ding LJ, Wu JH, Li CY, Peng XX, Zhang S, Yi GJ. Comparative transcriptomics analysis reveals difference of key gene expression between banana and plantain in response to cold stress. BMC Genomics. 2015;16:1.
Sergeant K, Kieffer P, Dommes J, Hausman J-F, Renaut J. Proteomic changes in leaves of poplar exposed to both cadmium and low-temperature. Environ Exp Bot. 2014;106:112–23.
Xu J, Li Y, Sun J, Du L, Zhang Y, Yu Q, Liu X. Comparative physiological and proteomic response to abrupt low temperature stress between two winter wheat cultivars differing in low temperature tolerance. Plant Biology. 2013;15:292–303.
Bano S, Aslam M, Saleem M, Basra S, Aziz K. Evaluation of maize accessions under low temperature stress at early growth stages. J Anim Plant Sci. 2015;25:392–400.
Wu YQ, Wei W, Pang XY, Wang XF, Zhang HL, Dong B, Xing YP, Li XG, Wang MY. Comparative transcriptome profiling of a desert evergreen shrub, Ammopiptanthus mongolicus, in response to drought and cold stresses. BMC Genomics. 2014;15:1.
Kashyap P, Deswal R. CBF-dependent cold stress signaling relevant post translational modifications. In: Stress Signaling in Plants: Genomics and Proteomics Perspective, Volume 1. New York: Springer; 2013. p. 105–122.
Barrero-Gil J, Salinas J. Post-translational regulation of cold acclimation response. Plant Sci. 2013;205:48–54.
Ichimura K, Mizoguchi T, Yoshida R, Yuasa T, Shinozaki K. Protein phosphorylation and dephosphorylation in environmental stress responses in plants. Adv Bot Res. 2000;32:355–77.
Gao JF. Plant physiology laboratory technology. Xi’an: World Publishing Corporation; 2000. p. 135–49.
García-Sevillano MA, Abril N, Fernández-Cisnal R, García-Barrera T, Pueyo C, López-Barea J, Gómez-Ariza JL. Functional genomics and metabolomics reveal the toxicological effects of cadmium in Mus musculus mice. Metabolomics. 2015;11:1432-50.
Yohannes E, Chang J, Christ GJ, Davies KP, Chance MR. Proteomics analysis identifies molecular targets related to diabetes mellitus-associated bladder dysfunction. Mol Cell Proteomics. 2008;7:1270–85.
Li BY, Cheng M, Gao HQ, Ma YB, Xu L, Li XH, Li XL, You BA. Back-regulation of six oxidative stress proteins with grape seed proanthocyanidin extracts in rat diabetic nephropathy. J Cell Biochem. 2008;104:668–79.
We acknowledge Shanghai Applied Protein Technology and Beijing Genomics Institute at Shenzhen for its assistance in original data processing. We would like to thank many other people for helping with the sampling. This study was supported by the Special Fund for Forest Scientific Research in the Public Welfare (201504103) and National Natural Science Foundation of China (31470616).
This study was supported by the Special Fund for Forest Scientific Research in the Public Welfare (201504103) and National Natural Science Foundation of China (31470616).
Availability of data and materials
Proteomics raw data have been submitted to PRIDE. Data are available via ProteomeXchange with identifier PXD004609.
CH carried out the experiments, analyzed the data and drafted the manuscript. AD. and HL. collected samples. GG participated in the RT-PCR analysis. JZ. conceived of the study. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
The 1D- and 2D- DIGE maps of H. rhamnoides leaves under control and low-temperature stress. (JPG 635 kb)
The identified protein summary of H. rhamnoides cv. ‘Chuisk’ responding to low-temperature stress using MALDI-TOF-TOF. Table S1B. The identified protein detail of H. rhamnoides cv. ‘Chuisk’ responding to low-temperature stress using MALDI-TOF-TOF. Table S2. GO enrichment terms under low-temperature stress with Arabidopsis thaliana homologues using GOEAST. Table S3. KEGG pathways of protein spots under low-temperature stress. Table S4. Primers used for detecting the gene expression by reverse real-time PCR. (XLS 171 kb)