Figures
Abstract
Low temperature at the seedling stage is a major damaging factor for rice production in southern China. To better understand the cold response of cultivated and wild rice, cold-sensitive cultivar 93–11 (Oryza sativa L. ssp. Indica) and cold-resistant hybrid wild rice DC907 with a 93–11 genetic background were used for a quantitative proteomic analysis with tandem mass tags (TMT) in parallel. Rice seedlings grown for four weeks at a normal temperature (25°C) were treated at 8–10°C for 24, 72 and 120 h. The number of differentially expressed proteins increased gradually over time in the cold-exposed rice in comparison with the untreated rice. A total of 366 unique proteins involved in ATP synthesis, photosystem, reactive oxygen species, stress response, cell growth and integrity were identified as responding to cold stress in DC907. While both DC907 and 93–11 underwent similar alterations in proteomic profiles in response to cold stress, DC907 responded in a prompter manner in terms of expressing cold-responding proteins, maintained a higher level of photosynthesis to power the cells, and possessed a stable and higher level of DIR proteins to prevent the plant from obtaining irreversible cell structure damage. The observations made in this study may lay a new foundation for further investigation of cold sensitivity or tolerance mechanisms in rice.
Citation: Wang J, Wang J, Wang X, Li R, Chen B (2018) Proteomic response of hybrid wild rice to cold stress at the seedling stage. PLoS ONE 13(6): e0198675. https://doi.org/10.1371/journal.pone.0198675
Editor: Niranjan Baisakh, Louisiana State University, UNITED STATES
Received: February 7, 2018; Accepted: May 23, 2018; Published: June 7, 2018
Copyright: © 2018 Wang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: This work was supported in part by the grants from Guangxi Department of Science and Technology (GKG 1123001-3B, GKZ 14121001-1-5), and Guangxi University Key Program of Funds (XDZ 110082) to RL. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Low temperature is a major environmental stress affecting plant growth. Chilling stress causes water reduction and osmotic changes in the cellular milieu and suppresses the activities of cellular macromolecules, resulting in reduced growth and extensive losses in agricultural production [1]. Rice (Oryza sativa), a monocot plant and widely grown as food crop in tropical and subtropical areas, is particularly sensitive to cold stress at the seedling and flowering stages [2, 3]. Molecular genetic studies have already identified components of cold tolerance, such as CTB4a, which confers cold resistance by mediating ATP supply [4], and the WRKY gene superfamily in rice [5]. COLD1, as one of the best-characterized rice genes, is considered as a regulator of G-protein signaling (RGS) that regulates Ca2+ signaling in cells and confers chilling tolerance in rice [2, 6]. The dehydrin gene OsDhn1 has been identified as being highly expressed in developing seeds under low temperatures and protects rice floral organs against abiotic stress [7]. qCTS-9, found in hybrid rice under different cold environments, was confirmed to be a functional gene associated with cold tolerance at rice seedling stage [8]. In a genome-wide association mapping of cold tolerance in cultivated rice from rice diversity panel 1 (RDP1), 87 cold tolerance-related quantitative trait loci (QTLs) with significant enrichment for genes related to lipid metabolism, response to stress and oxygen binding were identified [9, 10].
Recently, proteomic technologies have been used to monitor and characterize protein profiles in rice [11]. For example, two-dimensional gel electrophoresis (2-DE) and isobaric tags labelling approaches were used to monitor the proteomic response of rice to the cold treatment and proteins involved in energy metabolism, transport, photosynthesis, precursor metabolites generation, histones and vitamin B biosynthesis, which were found to be differentially expressed by cold stress [12–14]. However, to date, only a limited number of proteins in the cold-response pathway have been identified.
Early season rice in South China generally suffers from cold weather characterized by a temperature drop to approximately 10°C or lower that results in seedling rot one in every three years in mid- to late March [15]. However, wild rice (Oryza rufipogon Griff.) in the same region survives the cold stress. Efforts to introduce the cold tolerance trait of wild rice into cultivar rice have been carried out by crossing cultivar rice with wild rice. To better understand the cold resistance mechanism of wild rice, a cold-tolerant hybrid wild rice DC907 with cultivar 93–11 genetic background and cold-sensitive 93–11 were investigated in parallel in this study by a comparative proteomics approach. Our results show that cold-tolerant DC907 was different from cold-sensitive 93–11 in its protein expression pattern. While a small portion of the differentially expressed proteins match those previously reported, a large proportion of cold stress-induced proteins were reported for the first time.
Materials and methods
Rice growth and cold treatment conditions
Seeds of the indica rice cultivar 93–11 and hybrid wild rice DC907, derived from crossing of Guangxi wild rice (Oryza rufipogon Griff.) with 93–11 and sequential back cross with 93–11 for 4 rounds, were germinated in soil and grown in a phytotron with a 12-h day/night cycle, at 25°C in the day and 18°C at night. Seedlings at the four-leaf stage were subject to cold treatment at 10°C in the day and 8°C at night to simulate natural cold conditions for varied time durations. Seedlings were separated into groups for varied cold treatment conditions. Each group contained a total of 100 individual plants. After cold treatment for a fixed amount of time, plants were then transferred to an environment of 25°C in the day and 18°C at night for 5 days for survival rate determination, defined as the ratio of surviving plants to total plants. The light intensity was set at 30000 lux.
Preparation of protein samples for TMT analysis
Whole rice seedlings were ground into powder with liquid nitrogen and then five volumes of pre-cold acetone containing 10% trichloroacetic acid (TCA) and 0.07% β-mercaptoethanol was added. The mixture was kept at -20°C overnight and centrifuged at 18,000 g for 30 min. The crude protein pellet was washed with pre-cold acetone containing 0.07% β-mercaptoethanol three times by centrifugation at 18,000 g. After vacuum drying, lysis buffer (7 M urea, 2% 3-((3-cholamidopropyl) dimethylammonio)-1-propanesulfonate (CHAPS), 40 mM Tris, 1 mM phenylmethanesulfonyl fluoride (PMSF)) was added to dissolve the protein pellet that was then centrifuged at 18,000 g for 30 min to remove debris. Five volumes of pre-chilled acetone was added to the supernatant and kept at -20°C overnight. The mixture was centrifuged at 18,000 g for 30 min, and the pellet was dissolved in triethylamine borane (TEAB, 100 mM). The protein concentration was determined using the Bradford method [16].
An amount of 100 μg of protein was mixed with 5 μl of 200 mM tris (2-carboxyethyl) phosphine (TCEP) and incubated at 55°C for 60 min. Then, 5 μl of 375 mM iodoacetamide was added and incubated for another 30 min in the dark. Six volumes of pre-cold acetone was added to precipitate proteins overnight. The protein mixture was centrifuged at 18,000 g for 30 min, and the pellet was dried at room temperature and dissolved in 100 μl of triethylamonium bicarbonat (TEAB). Trypsin solution (2.5 μg/100 μg protein) was used to digest protein samples at 37°C overnight. The peptides were labeled with 41 μl of the TMTsixplex label reagent set (Thermo Fisher Scientific, CAT 90061) and incubated for 60 min at room temperature. Then, 8 μl of 5% hydroxylamine was used to quench the labeling reaction for 15 min, and the labeled peptides were stored at -80°C.
Strong cation exchange chromatography
A PolyLC polysulfoethyl aspartamide column (100 mm X 2.1 mm, 5 μm, 300Å pore size) was used on a Waters high-performance liquid chromatography system (HPLC, Waters, series 2695) for off-line strong-cation exchange (SCX) chromatography fractionation. A 40 min gradient elution of 100% solvent A (10 mM monopotassium phosphate, 15% acetonitrile) to 100% solvent B (500 mM potassium chloride in solvent A) at 200 μl/min flow rate was performed. The SCX elution was monitored under a Waters 2998 PDA detector module (220 nm) and collected into 25 fractions for further mass spectrometry (MS) analysis.
Nanoflow LC-MS/MS analysis and data processing
The SCX fractions were loaded onto the trap column (nanoViper C18, 75 μm X 2 cm) and then eluted using capillary analytical column (nanoViper C18, 50 μm X 15 cm) at a 300 nl/min flow rate using Easy-nLC 1000 nanoflow liquid chromatography system (Thermo Fisher Scientific). The linear gradient for peptide elution was from 95% solvent A (0.1% formic acid) and 5% solvent B (0.1% formic acid, 98% acetonitrile) to 40% solvent B for a 60 min program. The peptides from the untreated control and cold treated samples were labeled, mixed and fractionated by SCX chromatography and then analyzed using LTQ-Orbitrap Elite hybrid mass spectrometer system. Three biological replicates for each sample were performed.
The scan range of mass spectrometric analysis was set at 350–1800 m/z in a data-dependent mode. The survey scan was set at 400 m/z with a mass resolution of 60,000. Tandem mass spectrometry (MS/MS or MS2) was preceded with ten of the most intense precursor ions in the collision-induced dissociation (CID) mode with 35% normalized collision energy. MS2 spectrum was acquired in the ion trap analyzer at normal speed. The software Proteome Discoverer 1.3 was used to search the mass spectrometric data against rice genome database v7.0 (http://rice.plantbiology.msu.edu/). Search parameters were set as a standard method: 2 missed cleavages using trypsin as endoprotease, lysine residues as fixed modification, peptide N-termini as variable modification, 10 ppm precursor ion mass tolerance, 0.8 Da fragment mass tolerance, and 1% maximum false discovery rate (FDR). The identified proteins were filtered with high peptide confidence. Proteins with a 1.5-fold change (p<0.05) were considered to be differentially expressed.
Western blot analysis
A Western blot analysis was carried out according to a previous study [17]. Samples of 40 μg of total protein were loaded on sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) gel and transferred to polyvinylidene fluoride (PVDF) membranes using a semi-dry transfer unit after electrophoretic separation. Specific antibodies against ribulose bisphosphate carboxylase oxygenase (Rubisco) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) were purchased from Abcam (California, USA) and used to detect and verify the protein expression level in different samples by enhanced chemiluminescence (SuperSignal West Pico Substrate, Thermo Scientific).
Bioinformatics
Identified differentially expressed proteins were classified according to Gene Ontology (GO) and KEGG. The protein information from the Database of Rice Genome Annotation Project [18] was converted to corresponding access numbers of the UniProt protein database and classified using QuickGo online annotation tool (http://www.ebi.ac.uk/QuickGO/GMultiTerm) [19]. The protein network of differentially expressed proteins was performed using agriGO v2.0 [20].
Results and discussion
Cold-response proteins and their time course
Cold stress that is harmful to the seedlings of early season rice in South China (mid- to late March) is typically at approximately 10°C [15]. Thus, temperatures of 8–10°C were selected for cold treatment in this study. Based on our previous observations that 93–11 primarily survived cold stress for 24 h but died completely after 120 h (survival rate was counted at day 5 after the stress was relieved) and wild rice survived at both conditions, we opted to use 24, 72, and 120 h as the times for proteomic analysis. As seen in Fig 1, the survival rates of 93–11 after cold stress for 24 h, 72 h, and 120 h were 90%, 20%, and 0%, respectively, while the survival rates of DC907 were 100% at all time points. A total of 1781 unique proteins were identified by TMT labeling from cold-stressed DC907 and 93–11 (S1 Table). By comparing protein abundance between the two accessions, 99 proteins (75 up- and 24 down-regulated) with changes greater than 1.5-fold were found after the first 24 h of cold treatment (Table 1). As the cold treatment time was extended, the number of differentially expressed proteins slightly increased: 120 (37 up- and 83 down-regulated) at 72 h (Table 2) and 105 (46 up- and 59 down-regulated) at 120 h (Table 3). These proteins fall into different functional groups with obvious time course characteristics; ATP synthesis, photosystems and other functional groups increased while DNA binding and transcription, cell growth and integrity, and structural protein decreased as the cold treatment time was extended (Fig 2A). The structural proteins, DNA binding proteins, and transcriptional factors were among the most differentially expressed at the time point of 24 h, and the most significant changes found at 72 h were the proteins involved in stress response. Accumulations of protein related to photosynthesis and ATP synthesis were affected more severely in 93–11 as the cold treatment proceeded. A total of 366 unique proteins from DC907 were found to significantly change in terms of accumulation under cold treatment for all three time points when using the untreated sample as a reference (S2 Table). More interestingly, we found that most of the differentially expressed proteins at 24 h returned to the level at time zero at 72 h, but changed to an even higher level as the cold treatment continued to 120 h in DC907 (Fig 2B). However, no such change was found in 93–11. According to these observations, we assume that self-regulation in the hybrid wild rice is important for cold resistance to ensure that the expression level of key proteins would not overload and threaten plant survival. Compared with previous studies of cold-response proteomics of cultivar rice, a proportion of proteins were identified that included rubisco and GAPDH whose expression level were both decreased under cold stress [12, 13]. As shown in Fig 3, the expression level of rubisco and GAPDH in DC907 detected by Western blotting also decreased, in good accordance with TMT quantification.
Rice seedlings were cultured in a phytotron under light and temperature controls (the day/night cycle: 12 h with 25°C and 12 h with 18°C). The samples labeled as 24 h, 72 h and 120 h were treated under corresponding artificial low temperatures and recovered for 5 days under normal culture conditions. A, comparative phenotypes of 120 h low-temperature treated hybrid wild and cultivar rice after 5 days recovery. B, cultivar rice phenotypes of different low-temperature treatment times after 5 days of recovery.
A, numbers of differentially expressed proteins between hybrid wild rice DC-907 and cultivar rice 93–11 at the same cold treatment time points. B, numbers of differentially expressed proteins in hybrid wild rice DC-907 at different cold treatment time points.
An amount of 40 μg of total protein from different samples was used for the Western blot analysis by the enhanced chemiluminescence (ECL) method. The specific antibodies against rubisco (1:1000) and GAPDH (1:1000) were used to detect the corresponding protein expressions. The change trends of these two proteins was consistent with the observations from TMT labeling.
In terms of time frames, the differentially expressed proteins in DC907 contained several functional groups: ATP synthesis, photosystem, reactive oxygen species (ROS), stress response, transcription factors, structural proteins, and cell growth and integrity (S2 Table). Under cold stress, the proteome pattern of hybrid wild rice showed a more sensitive and faster change than cultivar rice, e.g., the vigorous change in protein functional classification was seen at 24 h (Fig 4A), but a similar change occurred in cultivar rice at 72 h (Fig 4B). The delay of cold response in cultivar rice could be a crucial reason for low survival rates. The protein networks of differentially expressed proteins in DC907 show that the centrality of protein change is mainly related to cell structure, stress response, ATP synthesis and photosynthesis (S1 Fig). These networks of functional proteins were generally consistent with those from cultivated rice previously reported using the 2-DE method [13, 14], suggesting that the timely response to cold stress and self-regulation of wild rice is more important than a change in single proteins during cold stress.
GO_BP, Gene Ontology Biological Processes. GO_CC: Gene Ontology Cellular Component. 93–11, C0 represents the control group without cold treatment; C1 represents cold treatment for 24 h; C2 represents cold treatment for 72 h; C3 represents cold treatment for 120 h. DC907, T0 represents the control group without cold treatment; T1 represents cold treatment for 24 h; T2 represents cold treatment for 72 h; T3 represents cold treatment for 120 h.
Functional grouping of cold-responding proteins
Differentially expressed proteins between DC907 and 93–11 with assigned functions at all three cold treatment times are summarized in Fig 5. Some of these proteins have been shown or implicated to have functions against cold stress.
Non-specific lipid-transfer proteins.
Non-specific lipid-transfer proteins (nsLTPs) are a group of small lipophilic proteins that accumulate between the plant epidermis and cell wall [21]. Previous reports have revealed that nsLTPs participate in the process of plant biotic and abiotic stress resistance [22–24]. As seen from the comparative proteomic results between 93–11 and DC907, approximately 10% of the cold-responsive proteins were related to stress response. Although the nsLTPs were both up-regulated in DC907 and 93–11 after cold stress, the comparative expression level of nsLTPs was higher in hybrid wild rice DC907 in the first 24 h and 120 h (Tables 1 and 2), implying that nsLTPs may play an important part in cold tolerance in DC907.
Heat shock proteins.
Heat shock proteins (HSPs) have been shown to facilitate plant adaptation to environmental changes [25]. A higher level of HSPs, acting as molecular chaperones, may help plants adapt to abnormal temperature, light, drought and salt [26]. As shown in S2 Table, many small HSPs were found to be up-regulated by cold stress in DC907, similar to those found for the cold-tolerance response in japonica rice [13].
ROS-related proteins
Reactive oxygen species (ROS) function to oxidize the harmful substrates that may produce and accumulate in the cell. It is now known that ROS-derived signals regulate plant growth, development and stress adaption [30, 31] and are crucial for removing harmful substances from cells and avoiding plant frostbite. Under stress conditions, the ROS level could be increased by 3–10 folds to help the plant adapt to the harsh environment [32]. The ROS related proteins were mainly down-regulated under cold stress in 93–11. The death of 93–11 indicated that harmful substrates may not have been efficiently scavenged in a timely manner and damage of the cells may occur eventually, whereas the expression level of this kind of proteins was relatively stable in DC907 (S1 Table). Thus, insufficiency in ROS may result in plant cell dysfunction and cell death [33, 34].
Cell structure proteins
A total of 60 differentially expressed proteins related to cell growth, integrity and structure were found (S2 Table). A large amount of ribosomal proteins was found to increase in expression during cold stress both in hybrid DC907 and cultivar rice 93–11, but the time points were different; for 93–11, the highest peak was at 72 h and returned to a normal level at 120 h, and for hybrid DC907, the highest peak was at 24 h and returned to a normal level at 72 h, thus a much faster response in the cold-tolerant hybrid DC907. In Escherichia coli, a 70 kDa ribosomal-associated protein (CsdA) was induced when the temperature was shifted from 37 to 15°C, and this protein was further demonstrated to be involved in derepression of HSPs and cell growth at low temperatures [35]. In soybean, three low-temperature inducible ribosomal proteins were found to be increasingly expressed after cold treatment [36].
Dirigent (DIR) proteins were identified as being induced at cold conditions in cultivated rice for the first time in this study. This kind of protein was reported to be involved in lignification and to respond to pathogen infection and abiotic stress in plants [37, 38]. In response to a temperature shift, the expression level of DIRs showed a relatively stable pattern in DC907, similar to hardy plants [39]. In contrast, three DIR proteins were all down-regulated and decreased substantially following the prolonged cold treatment in cultivar rice 93–11. It was reported that a decrease in DIR proteins weakened the process of lignification, which is crucial for the structural integrity of the plant cell wall and cell wall apposition (CWA)-mediated defense [40]. Thus, it is speculated that the decreased expression of DIR proteins in 93–11 results in the vulnerability of indica rice to the cold stress since lignification may prevent cells from collapsing and responding to abiotic stress [41].
Proteins related to photosynthesis and energy metabolism
In the current study, most of the differentially expressed proteins involved in photosynthesis are from photosystem II (PSII), indicating that PSII is more sensitive than photosystem I (PSI) to cold stress [42]. When the cold treatment started, the expression level of chlorophyll a-b binding proteins from the light-harvesting complex (LHC) as a light receptor decreased very fast and was significantly lower in DC907 than in 93–11 (Tables 1–3). This phenomenon is comparable to the previous observation that photoinhibition was a protection mechanism for cold tolerant plants but more significant and protective for cultivar rice in low temperature after long time exposure [43, 44].
Stress-induced inhibition of plant photosynthesis is always coupled with a loss of ATP and ATP synthases [45, 46]. Thus, most of the proteins related to ATP synthesis were found to decrease gradually in DC907 (S2 Table). The decreased ATP synthase resulted in a shortage of ATP, implying a weakened biological activity. The reduction in ATP and ATP synthases directly relates to the reduction in photosynthetic energy captured at low temperatures, and a shortage of energy supply would result in the restriction of normal metabolic processes of plant cells. It is believed that this is a protective mechanism for plant cells under abnormal cold stress [47, 48].
Interaction networks among cold induced proteins
As show in S1 Fig, protein interaction networks constructed with differentially expressed proteins from hybrid wild rice contain biological processes, molecular functions, and cellular components. In biological processes, proteins that function in photosynthesis, metabolite and energy generation, protein translation, and stress response are the highest positively regulated in DC907; in terms of molecular function, proteins that function in structural activities are most positively regulated in DC907. However, the largest group of regulated proteins at the highest level was in the cellular component domain, including membranes, macromolecule complexes, vacuoles, ribosomes, mitochondria and other organelles. An important observation is that hybrid wild rice responded to cold stress in a more timely manner by mobilizing its signal transduction and self-regulation mechanisms, similar to other cold tolerant plants studied [49].
Conclusions
In this work, the proteomes of cold-resistant DC907 with 93–11 genetic background and the cold-sensitive 93–11 were compared. The protein expression level of several important functional categories, including photosynthesis, energy generation, ROS, cell growth and development, were found to be changed under cold stress. While both DC907 and 93–11 underwent similar alterations in proteomic profiles to cold stress, DC907 responded in a prompter manner in expressing cold-response proteins, maintained a higher level of photosynthesis to power the cells, and possessed a stable and higher level of DIR proteins to prevent the plant from obtaining irreversible cell structure damage induced by ROS activity (Fig 6). Since DC907 carries chromosome fragments of wild rice, future studies should focus on the genetic elements of wild rice that confer the cold tolerant trait in DC907.
Supporting information
S1 Fig. The interaction protein networks of hybrid wild rice DC907.
A, biological process; B, molecular function; C, cellular component.
https://doi.org/10.1371/journal.pone.0198675.s001
(TIF)
S1 Table. Identified information of rice proteins by TMT labeling.
https://doi.org/10.1371/journal.pone.0198675.s002
(XLSX)
S2 Table. List of differentially expressed proteins of hybrid wild rice DC907 during the cold treatment.
https://doi.org/10.1371/journal.pone.0198675.s003
(DOC)
Acknowledgments
This work was supported in part by grants from the Guangxi Department of Science and Technology (GKG 1123001-3B, GKZ 14121001-1-5) and Guangxi University Key Program of Funds (XDZ 110082).
References
- 1. Xiong L, Schumaker KS, Zhu JK. Cell signaling during cold, drought, and salt stress. The Plant cell. 2002;14 Suppl:S165–83. pmid:12045276; PubMed Central PMCID: PMC151254.
- 2. Ma Y, Dai X, Xu Y, Luo W, Zheng X, Zeng D, et al. COLD1 confers chilling tolerance in rice. Cell. 2015;160(6):1209–21. pmid:25728666.
- 3. Saito K, Miura K, Nagano K, Hayano-Saito Y, Araki H, Kato A. Identification of two closely linked quantitative trait loci for cold tolerance on chromosome 4 of rice and their association with anther length. Theoretical and Applied Genetics. 2001;103(6–7):862–8.
- 4. Zhang Z, Li J, Pan Y, Li J, Zhou L, Shi H, et al. Natural variation in CTB4a enhances rice adaptation to cold habitats. Nature communications. 2017;8:14788. pmid:28332574; PubMed Central PMCID: PMC5376651.
- 5. Ryu HS, Han M, Lee SK, Cho JI, Ryoo N, Heu S, et al. A comprehensive expression analysis of the WRKY gene superfamily in rice plants during defense response. Plant cell reports. 2006;25(8):836–47. pmid:16528562.
- 6. Manishankar P, Kudla J. Cold tolerance encoded in one SNP. Cell. 2015;160(6):1045–6. pmid:25768901.
- 7. Lee SC, Lee MY, Kim SJ, Jun SH, An G, Kim SR. Characterization of an abiotic stress-inducible dehydrin gene, OsDhn1, in rice (Oryza sativa L.). Molecules and cells. 2005;19(2):212–8. pmid:15879704.
- 8. Zhao J, Zhang S, Dong J, Yang T, Mao X, Liu Q, et al. A novel functional gene associated with cold tolerance at the seedling stage in rice. Plant biotechnology journal. 2017;15(9):1141–8. pmid:28173633; PubMed Central PMCID: PMC5552475.
- 9. Shakiba E, Edwards JD, Jodari F, Duke SE, Baldo AM, Korniliev P, et al. Genetic architecture of cold tolerance in rice (Oryza sativa) determined through high resolution genome-wide analysis. PloS one. 2017;12(3):e0172133. pmid:28282385; PubMed Central PMCID: PMC5345765.
- 10. Wang D, Liu J, Li C, Kang H, Wang Y, Tan X, et al. Genome-wide Association Mapping of Cold Tolerance Genes at the Seedling Stage in Rice. Rice. 2016;9(1):61. pmid:27848161; PubMed Central PMCID: PMC5110459.
- 11. Agrawal GK, Rakwal R. Rice proteomics: A move toward expanded proteome coverage to comparative and functional proteomics uncovers the mysteries of rice and plant biology. Proteomics. 2011;11(9):1630–49. pmid:21462347.
- 12. Neilson KA, Mariani M, Haynes PA. Quantitative proteomic analysis of cold-responsive proteins in rice. Proteomics. 2011;11(9):1696–706. pmid:21433000.
- 13. Cui S, Huang F, Wang J, Ma X, Cheng Y, Liu J. A proteomic analysis of cold stress responses in rice seedlings. Proteomics. 2005;5(12):3162–72. pmid:16078185.
- 14. Ji L, Zhou P, Zhu Y, Liu F, Li R, Qiu Y. Proteomic Analysis of Rice Seedlings Under Cold Stress. The protein journal. 2017;36(4):299–307. pmid:28555319.
- 15. Huang M, Zhang WX, Jiang LG, Zou YB. Impact of temperature changes on early-rice productivity in a subtropical environment of China. Field Crops Research. 2013;146:10–5.
- 16. Bradford MM. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Analytical biochemistry. 1976;72:248–54. pmid:942051.
- 17. Li X, Bai H, Wang X, Li L, Cao Y, Wei J, et al. Identification and validation of rice reference proteins for western blotting. Journal of experimental botany. 2011;62(14):4763–72. pmid:21705388; PubMed Central PMCID: PMC3192993.
- 18. Kawahara Y, de la Bastide M, Hamilton JP, Kanamori H, McCombie WR, Ouyang S, et al. Improvement of the Oryza sativa Nipponbare reference genome using next generation sequence and optical map data. Rice (N Y). 2013;6(1):4. Epub 2013/11/28. pmid:24280374.
- 19. Binns D, Dimmer E, Huntley R, Barrell D, O'Donovan C, Apweiler R. QuickGO: a web-based tool for Gene Ontology searching. Bioinformatics. 2009;25(22):3045–6. Epub 2009/09/12. pmid:19744993; PubMed Central PMCID: PMC2773257.
- 20. Tian T, Liu Y, Yan H, You Q, Yi X, Du Z, et al. agriGO v2.0: a GO analysis toolkit for the agricultural community, 2017 update. Nucleic acids research. 2017;45(W1):W122–W9. pmid:28472432; PubMed Central PMCID: PMC5793732.
- 21. Liu F, Zhang X, Lu C, Zeng X, Li Y, Fu D, et al. Non-specific lipid transfer proteins in plants: presenting new advances and an integrated functional analysis. Journal of experimental botany. 2015;66(19):5663–81. pmid:26139823.
- 22. Wei K, Zhong X. Non-specific lipid transfer proteins in maize. BMC plant biology. 2014;14:281. pmid:25348423; PubMed Central PMCID: PMC4226865.
- 23. Yu G, Hou W, Du X, Wang L, Wu H, Zhao L, et al. Identification of wheat non-specific lipid transfer proteins involved in chilling tolerance. Plant cell reports. 2014;33(10):1757–66. pmid:25037996.
- 24. Salcedo G, Sanchez-Monge R, Barber D, Diaz-Perales A. Plant non-specific lipid transfer proteins: an interface between plant defence and human allergy. Biochimica et biophysica acta. 2007;1771(6):781–91. pmid:17349819.
- 25. Timperio AM, Egidi MG, Zolla L. Proteomics applied on plant abiotic stresses: role of heat shock proteins (HSP). Journal of proteomics. 2008;71(4):391–411. Epub 2008/08/23. pmid:18718564.
- 26. Sun W, Van Montagu M, Verbruggen N. Small heat shock proteins and stress tolerance in plants. Biochimica et biophysica acta. 2002;1577(1):1–9. Epub 2002/08/02. pmid:12151089.
- 27. Boucher V, Buitink J, Lin X, Boudet J, Hoekstra FA, Hundertmark M, et al. MtPM25 is an atypical hydrophobic late embryogenesis-abundant protein that dissociates cold and desiccation-aggregated proteins. Plant, cell & environment. 2010;33(3):418–30. pmid:20002332.
- 28. C ND, Danyluk J, Wilson KE, Pocock T, Huner NP, Sarhan F. Cold-regulated cereal chloroplast late embryogenesis abundant-like proteins. Molecular characterization and functional analyses. Plant physiology. 2002;129(3):1368–81. pmid:12114590; PubMed Central PMCID: PMC166530.
- 29. Ukaji N, Kuwabara C, Takezawa D, Arakawa K, Fujikawa S. Cold acclimation-induced WAP27 localized in endoplasmic reticulum in cortical parenchyma cells of mulberry tree was homologous to group 3 late-embryogenesis abundant proteins. Plant physiology. 2001;126(4):1588–97. pmid:11500557; PubMed Central PMCID: PMC117158.
- 30. Gechev TS, Van Breusegem F, Stone JM, Denev I, Laloi C. Reactive oxygen species as signals that modulate plant stress responses and programmed cell death. BioEssays: news and reviews in molecular, cellular and developmental biology. 2006;28(11):1091–101. pmid:17041898.
- 31. Van Breusegem F, Dat JF. Reactive oxygen species in plant cell death. Plant physiology. 2006;141(2):384–90. pmid:16760492; PubMed Central PMCID: PMC1475453.
- 32. Polle A. Dissecting the superoxide dismutase-ascorbate-glutathione-pathway in chloroplasts by metabolic modeling. Computer simulations as a step towards flux analysis. Plant physiology. 2001;126(1):445–62. pmid:11351106; PubMed Central PMCID: PMC102317.
- 33. Wu J, Sun Y, Zhao Y, Zhang J, Luo L, Li M, et al. Deficient plastidic fatty acid synthesis triggers cell death by modulating mitochondrial reactive oxygen species. Cell research. 2015;25(5):621–33. pmid:25906995; PubMed Central PMCID: PMC4423084.
- 34. Zhang L, Xing D. Methyl jasmonate induces production of reactive oxygen species and alterations in mitochondrial dynamics that precede photosynthetic dysfunction and subsequent cell death. Plant & cell physiology. 2008;49(7):1092–111. pmid:18535010.
- 35. Jones PG, Mitta M, Kim Y, Jiang W, Inouye M. Cold shock induces a major ribosomal-associated protein that unwinds double-stranded RNA in Escherichia coli. Proceedings of the National Academy of Sciences of the United States of America. 1996;93(1):76–80. pmid:8552679; PubMed Central PMCID: PMC40181.
- 36. Kim KY, Park SW, Chung YS, Chung CH, Kim JI, Lee JH. Molecular cloning of low-temperature-inducible ribosomal proteins from soybean. Journal of experimental botany. 2004;55(399):1153–5. Epub 2004/03/17. pmid:15020631.
- 37. Zulak KG, Khan MF, Alcantara J, Schriemer DC, Facchini PJ. Plant defense responses in opium poppy cell cultures revealed by liquid chromatography-tandem mass spectrometry proteomics. Molecular & cellular proteomics: MCP. 2009;8(1):86–98. pmid:18682378.
- 38. Jin-Long G, Li-Ping X, Jing-Ping F, Ya-Chun S, Hua-Ying F, You-Xiong Q, et al. A novel dirigent protein gene with highly stem-specific expression from sugarcane, response to drought, salt and oxidative stresses. Plant cell reports. 2012;31(10):1801–12. pmid:22696141.
- 39. Ralph S, Park JY, Bohlmann J, Mansfield SD. Dirigent proteins in conifer defense: gene discovery, phylogeny, and differential wound- and insect-induced expression of a family of DIR and DIR-like genes in spruce (Picea spp.). Plant molecular biology. 2006;60(1):21–40. pmid:16463097.
- 40. Bhuiyan NH, Selvaraj G, Wei Y, King J. Role of lignification in plant defense. Plant signaling & behavior. 2009;4(2):158–9. Epub 2009/08/04. pmid:19649200; PubMed Central PMCID: PMC2637510.
- 41. Le Gall H, Philippe F, Domon JM, Gillet F, Pelloux J, Rayon C. Cell Wall Metabolism in Response to Abiotic Stress. Plants. 2015;4(1):112–66. pmid:27135320; PubMed Central PMCID: PMC4844334.
- 42. Lei Y-b, Zheng Y-l, Dai K-j, Duan B-l, Cai Z-q. Different responses of photosystem I and photosystem II in three tropical oilseed crops exposed to chilling stress and subsequent recovery. Trees. 2014;28(3):923–33.
- 43. Allen DJ, Ort DR. Impacts of chilling temperatures on photosynthesis in warm-climate plants. Trends in plant science. 2001;6(1):36–42. Epub 2001/02/13. pmid:11164376.
- 44. Huner NP, Oquist G, Hurry VM, Krol M, Falk S, Griffith M. Photosynthesis, photoinhibition and low temperature acclimation in cold tolerant plants. Photosynthesis research. 1993;37(1):19–39. pmid:24317651.
- 45. Sobczyk EA, Marszalek A, Kacperska A. ATP involvement in plant tissue responses to low temperature. Physiologia Plantarum. 1985;63(4):399–405.
- 46. Tezara W, Mitchell VJ, Driscoll SD, Lawlor DW. Water stress inhibits plant photosynthesis by decreasing coupling factor and ATP. Nature. 1999;401:914–7.
- 47. Lawlor DW. Limitation to photosynthesis in water-stressed leaves: stomata vs. metabolism and the role of ATP. Annals of botany. 2002;89 Spec No:871–85. Epub 2002/07/10. pmid:12102513.
- 48. Graham D, Patterson BD. Responses of plants to low, nonfreezing temperatures: proteins, metabolism, and acclimation. Annual Review of Plant Physiology. 1982;33:347–72.
- 49. Yadav SK. Cold stress tolerance mechanisms in plants. A review. Agronomy for Sustainable Development. 2010;30(3):515–27.