Yunze Lu1
Ruiqiong Li1
Ruochen Wang1
Xiaoming Wang1
Weijun Zheng1
Qixin Sun1,2
Shaoming Tong3
Shaojun Dai4
Shengbao Xu1*- 1State Key Laboratory of Crop Stress Biology for Arid Areas, College of Agronomy, Northwest A&F University, Yangling, China
- 2Department of Plant Genetics and Breeding, China Agricultural University, Beijing, China
- 3College of Life Sciences, Liaoning Normal University, Dalian, China
- 4Development Center of Plant Germplasm Resources, College of Life and Environmental Sciences, Shanghai Normal University, Shanghai, China
Hexaploid wheat (Triticum aestivum L.) is an important food crop but it is vulnerable to heat. The heat-responsive proteome of wheat remains to be fully elucidated because of previous technical and genomic limitations, and this has hindered our understanding of the mechanisms of wheat heat adaptation and advances in improving thermotolerance. Here, flag leaves of wheat during grain filling stage were subjected to high daytime temperature stress, and 258 heat-responsive proteins (HRPs) were identified with iTRAQ analysis. Enrichment analysis revealed that chlorophyll synthesis, carbon fixation, protein turnover, and redox regulation were the most remarkable heat-responsive processes. The HRPs involved in chlorophyll synthesis and carbon fixation were significantly decreased, together with severe membrane damage, demonstrating the specific effects of heat on photosynthesis of wheat leaves. In addition, the decrease in chlorophyll content may result from the decrease in HRPs involved in chlorophyll precursor synthesis. Further analysis showed that the accumulated effect of heat stress played a critical role in photosynthesis reduction, suggested that improvement in heat tolerance of photosynthesis, and extending heat tolerant period would be major research targets. The significantly accumulation of GSTs and Trxs in response to heat suggested their important roles in redox regulation, and they could be the promising candidates for improving wheat thermotolerance. In summary, our results provide new insight into wheat heat adaption and provide new perspectives on thermotolerance improvement.
Introduction
Bread wheat (Tricticum aestivum L.) is an important staple food for human beings, however, this crop is vulnerable to heat stress and this can lead to substantial yield loss (Lobell and Tebaldi, 2014; Tack et al., 2015; Lesk et al., 2016). During the life cycle of wheat, the grain filling stage is a critical period for yield and quality, but this stage is particularly vulnerable to heat stress that can lead to greater detrimental effects in wheat production and quality than any of the other developmental stages (Zhao et al., 2007; ur Rehman et al., 2009; Farooq et al., 2011; Pradhan and Prasad, 2015). However, the mechanisms involved in the heat response and adaption in wheat, especially at the filling stage, are far from clearly understood.
Photosynthesis converts light energy into carbohydrates, a process that is sensitive to heat stress. Accumulating evidence has shown that heat stress affects multiple processes that tightly related to photosynthesis, including thylakoid membrane arrangement and damage (Prasad et al., 2008; Narayanan et al., 2015), chlorophyll content reduction (Ristic et al., 2007; Prasad et al., 2008), and Rubisco quantity and activity decline. The deactivation of Rubisco is the primary limiting factor for photosynthesis loss (Law and Crafts-Brandner, 1999; Sharkey, 2005). Moreover, high temperature reduces the quantity of chlorophyll and hastens leaf senescence (Al-Khatib and Paulsen, 1984; Prasad et al., 2008), resulting in a shorter filling duration stage and the subsequent loss of grain weight (Zhao et al., 2007; Farooq et al., 2011).
Heat stress also triggers oxidative stress and activates the reactive oxygen species (ROS)-scavenging system (Mittler et al., 2004; Suzuki et al., 2012). Severe or prolonged heat stress reduces the expression levels and enzyme activities of several major antioxidant enzymes, such as superoxide dismutase (SOD), catalase (CAT), and ascorbate peroxidase (APX), leading to ROS accumulation and a negative impact on plant growth and development (Mittler et al., 2004; Wang et al., 2011, 2015; Narayanan et al., 2015).
To counteract the detrimental effects of heat, plants have evolved multiple adaption systems. The transcriptional reprograming of the plant heat response is chartered by (1) the induction of heat shock factors and heat shock proteins (HSPs), a typical feature of the heat response; (2) the inhibition of general protein synthesis; (3) ROS burst and significant antioxidant accumulation; and (4) anabolism inhibition and catabolism activation (Morimoto, 1998; Qin et al., 2008; Mangelsen et al., 2011; Sarkar et al., 2014; Liu et al., 2015). Proteomics analyses have also revealed that HSPs, detoxification, protein synthesis and photosynthesis process are affected by heat stress (Rollins et al., 2013; Liao et al., 2014). However, in wheat, only HSPs have been characterized by proteomic analysis under heat stress in previous studies because of limited genomic information and technical restrictions (Nguyen et al., 1994; Majoul et al., 2003, 2004; Wang et al., 2015).
With updated genomic information and improvement in proteomic technology, more precise, larger-scale protein characterization in wheat has become practicable (Choulet et al., 2014; IWSGC, 2014). In this study, an isobaric tags for relative and absolute quantitation (iTRAQ)-based proteomics approach was used to analyze proteomic changes in wheat flag leaves under heat treatment. We identified 258 heat-responsive proteins (HRPs) in wheat flag leaves that were involved in several biological pathways. We also presented a detailed landscape of the reduction in photosynthesis caused by heat stress. These results enhanced our understanding of the wheat response and adaption to heat stress and provided new perspectives for improving wheat thermotolerance.
Materials and Methods
Plant Materials and Growth Conditions
Wheat plants (T. aestivum cv. Chinese Spring) were grown in a greenhouse and watered daily to avoid drought stress. The main stem ears were labeled when flowering. At 12 days after anthesis, the plants were transferred to growth chambers with a temperature regime of 24/17°C (14/10 h) day/night cycle. At 15 days after anthesis, half of the plants were applied a temperature regime of 37/17°C for 3 days and were harvested on the 4th day after 4 h under heat with light; the remaining plants were kept at 24/17°C as a control and were harvested at the same sampling time. For the experiments at the seedling stage, 10-day-old seedlings were subjected to the same conditions described above except that a different daytime temperature was used. All obtained materials were stored at −80°C for further analysis.
Measurement of the Chlorophyll and Thiobarbituric Acid Reactive Substances (TBARS) Content of Flag Leaves and 3,3′-Diaminobenzidine (DAB) Staining
The chlorophyll content of flag leaves was measured according to the method of Wang et al. (2011), and the TBARS content was determined according to Lee et al. (2007). The middle part of leaf was used to determine the H2O2 content by DAB staining, as described by Tian et al. (2013). All of these determinations were performed on three biological replicates.
Protein Preparation
About 0.5 g of flag leaf was ground into a fine powder in liquid nitrogen, transferred to a 5-mL Eppendorf tube, and sonicated in lysis buffer (8 M urea, 2 mM EDTA, 1% Triton-100, 65 mM DTT, and 1% protease inhibitor cocktail). After centrifugation (20,000 g, 4°C, 10 min), the debris was removed and the protein was precipitated with cold 15% trichloroacetic acid for 2 h at −20°C. After centrifugation, the pellet was washed with cold acetone three times. The protein was resuspended in buffer (8 M urea, 100 mM triethylammonium bicarbonate, pH 8.0) and the protein concentration was determined using the 2-D Quant kit (GE Healthcare) following the manufacturer's protocol. Three independent biological replicates were performed for each sample.
Trypsin Digestion and iTRAQ Labeling
Approximately 100 μg protein of each biological replicate was used for digestion. Firstly, the protein sample was reduced by the addition of 10 mM DTT for 1 h at 37°C, and then alkylated using 20 mM iodoacetamide for 45 min at room temperature in the dark. Subsequently, the protein sample was diluted with 0.1 M triethylammonium bicarbonate to reduce the urea concentration to <2 M. Then, trypsin (Promega) was added at a trypsin/protein ratio (w/w) of 1:50 for the first digestion overnight at 37°C, and the trypsin/protein ratio was adjusted to 1:100 for another 4 h digestion at 37°C to ensure the fully digestion of proteins. The tryptic peptides were desalted through a Strata X C18 SPE column (Phenomenex) and vacuum-dried, and then reconstituted in 500 mM triethylammonium bicarbonate and processed following the instructions of the 8-plex iTRAQ kit (Applied Biosystems). Briefly, one unit of the iTRAQ reagent (defined as the amount of reagent required to label 100 μg of protein) was thawed and reconstituted in 24 μl acetonitrile. The proteins of heat-stressed samples were labeled with 113, 114, and 116, and the proteins of control samples were labeled with 117, 119, and 121. The peptide mixtures of different groups were labeled differently by incubation for 2 h at room temperature and were then pooled, desalted, and dried by vacuum centrifugation.
HPLC Fractionation
The peptides were fractionated by high pH reverse-phase HPLC using an Agilent 300Extend C18 column (5 μm particles, 4.6 mm ID, 250 mm length). Peptides were separated by a gradient of 2–60% acetonitrile in 10 mM ammonium bicarbonate, pH 10, over 80 min into 80 fractions. Then, the peptides were combined into 18 fractions and dried by vacuum centrifugation.
LC-MS/MS Analysis
The collected fractions dissolved in 0.1% formic acid were separated on an EASY-nLC 1000 UPLC system. Firstly, the reconstituted samples were applied to a reversed-phase pre-column (Acclaim PepMap 100, Thermo Fisher Scientific), then separated on a reversed-phase analytical column (Acclaim PepMap RSLC, Thermo Fisher Scientific) with a gradient at a constant flow rate of 300 nL/min. The gradient was as follows: 7–20% solvent B (0.1% formic acid in 98% acetonitrile; 24 min), 20–35% solvent B (8 min), 35–80% solvent B (3 min), and 80% solvent B (4 min).
The resulting peptides were subjected to nanospray source and analyzed by tandem mass spectrometry (MS/MS) in Q Exactive™ plus (Thermo Fisher Scientific) coupled online to the UPLC system. A resolution of 70,000 in the orbitrap was used to detect intact peptides, and peptides were selected for MS/MS using 33% normalized collision energy. MS data was obtained using a data-dependent “top 20” procedure, which alternatively chooses the most abundant precursor ions with a threshold ion count above 2E4 in the MS survey scan (350–1,800 m/z). The dynamic exclusion duration was 10 s. A resolution of 17,500 was set for ion fragments detection. The electrospray voltage applied was 2.0 kV. Automatic gain control was applied to prevent overfilling of the Orbitrap; 5E4 ions were accumulated to generate MS/MS spectra. The fixed first mass was adjusted to 100 m/z.
Database Searching
The resulting MS/MS data were processed by Proteome Discoverer (Thermo scientific™ software) with Mascot search engine (v.2.3.0) and searched against UniprotTriticum_aestivum_4565.fasta (100,981 sequences in total; downloaded from Uniprot on June 5, 2015). Trypsin/P was selected as the enzyme with two maximum missing cleavages. Mass error of precursor ions and fragment ions was set to 10 ppm and 0.02 Da, respectively. Carbamidomethyl on Cys, iTRAQ-8plex (N-term) and iTRAQ-8plex (K) were specified as fixed modification and oxidation on Met was specified as variable modifications. False discovery rate was adjusted to ≤ 1%. The peptides with score ≥20 were kept for protein identification. Due to the complexity and redundancy of wheat genome, a rigorous standards were applied in protein identification, include (1) Mascot score ≥40; (2) number of unique peptides ≥2. To obtain reliable protein quantification, only the proteins that could be quantified at all six samples were assigned as quantitative proteins. To designate the HRPs, a Student's t-test P < 0.05 and fold change > 1.3 or < 0.77 across the heat-stressed and control samples were applied (Yentrapalli et al., 2013; Trevisan et al., 2015). The mass spectrometry proteomics data had been deposited to the ProteomeXchange Consortium via the PRIDE (Vizcaino et al., 2014, 2016) partner repository with the dataset identifier PXD006207 and doi: 10.6019/PXD006207.
Bioinformatics Analysis
The UniProt-GOA database (http://www.ebi.ac.uk/GOA/) and the InterProScan software (http://www.ebi.ac.uk/interpro/, v5.13-52.0) were used for gene ontology (GO) annotation. The unannotated HRPs were searched against the NCBI database further. After annotation, proteins were mapped to the pathways in Kyoto Encyclopedia of Genes and Genomes (KEGG) database (http://www.kegg.jp/). The InterPro database (http://www.ebi.ac.uk/interpro/) were used for domain annotation. Subcellular localizations of proteins were predicted by Wolfpsort (Horton et al., 2007).
Enrichment Analysis
Firstly, the number of (a) HRPs and (b) all the quantitative proteins in each domain, GO biological process and KEGG pathway term were calculated. Then the total number of (c) HRPs and (d) all the quantitative proteins annotated by domain (or the GO biological process/KEGG pathway) were calculated. Finally, these four numbers (a, b, c, d) were used to calculate Fisher's exact test P-values.
Quantitative Real-Time PCR (qRT-PCR) Analysis
Total RNA was extracted from 50 mg of leaf sample using the TRNzol reagent (TianGen, Beijing, China) following the manufacturer's protocol. About 1 μg RNA was reverse transcribed to cDNA using the PrimeScript™ RT reagent Kit with gDNA Eraser (Perfect Real Time) kit (Takara, Dalian, China) according to the kit protocol. qRT-PCR was performed on an Applied Biosystems 7300 Real-Time PCR System using the SYBR® Premix Ex Taq™ II (Tli RNaseH Plus) kit (Takara, Dalian, China) following the corresponding protocol.
Results
Proteome Profiles of Wheat Flag Leaves under Heat Stress
To investigate the proteome alterations of wheat flag leaves under heat stress, iTRAQ-based quantitative proteomics analysis was conducted (Figure 1A). Among the 3499 proteins identified, 1915 could be further quantified (Supplementary Tables 1, 2), and 258 were assigned as HRPs (Figure 1B), including 155 heat-increased proteins and 103 heat-decreased proteins (Figure 1C).
Figure 1. Proteins identification in this experiment. (A) General workflow of the experiment for iTRAQ-based proteome analysis of wheat with/without heat stress. (B) Schematic diagram of proteins characterization. (C) Number of HRPs.
The HRPs were classified into seven categories according to the GO database and the Uniprot knowledgebase (Figure 2A, Supplementary Table 3, including stress response (29.1%), metabolism (26.7%), protein synthesis/folding/degradation (17.8%), photosynthesis (9.3%), and other categories. The categories affected suggest that heat stress has an effect on diverse biological processes in wheat leaves.
Figure 2. Function categories and enrichment analysis of HRPs. Black bars represent decreased proteins, white bars represent increased proteins. (A) Functional classification of HRPs. (B) GO biological process analysis. The bars represent −log2(P) where P represents the Fisher' exact test P-values, and so as it in (C,D). (C) KEGG pathway analysis. (D) Domain enrichment analysis.
Enrichment Analysis of HRPs
To illustrate the major heat-responsive processes in wheat flag leaves, enrichment analysis was performed on the GO biological process, the KEGG pathway and the protein domain (Supplementary Table 4). Among the biological process, as expected, the GO term “response to stress,” and “protein folding” were significantly enriched in increased HRPs. Whereas, “tetrapyrrole biosynthetic processes,” “carbon fixation,” and “translation” were remarkably enriched in decreased HRPs (Figure 2B). These results indicated that the stress response and protein folding were activated, whereas translation and photosynthesis were impeded by heat stress.
In the pathway enrichment analysis, protein processing in endoplasmic reticulum (ER) and glutathione metabolism were significantly enriched, and most proteins in these pathways were heat increased. Whereas, porphyrin and chlorophyll metabolism, and carbon fixation in photosynthetic organisms were enriched, and most related proteins were decreased in heat conditions (Figure 2C). These results further suggested that protein processing and photosynthesis were remarkably affected by heat stress accompanied by redox system alterations.
In the domain enrichment analysis, HSPs (including HSP100, HSP90, HSP70, and HSP20), the glutathione-S-transferase (GST) domain, the thioredoxin (Trx)-like fold, the tetratricopeptide repeat-containing domain, and others were remarkably enriched as heat-induced domains (Figure 2D), indicating these domains may be the central heat responsive domains. Whereas, only the RNA recognition motif domain was enriched as a heat-suppressive domain. The domain enrichment analysis revealed that the stress response, protein folding, and redox regulation were all activated by the heat stress in flag leaves.
Taken together, the enrichment analyses based on the proteomics data suggested that protein synthesis and folding, redox regulation, and photosynthesis were the most heat-responsive processes.
HRPs in Protein Synthesis and Folding
The proteomics data suggested that protein translation was decreased during heat stress (Figure 2B). To clarify the translational changes in wheat leaves under heat stress, de novo protein synthesis was evaluated based on HRPs (Figure 3). The translation components (translation initiation factor IF-3, elongation factors P and G, peptide chain release factor 1, ribosomal proteins) showed decreased under heat stress. These results indicated that translation elongation was coincidentally interfered by heat, in addition to the previously reported heat-inhibition of translation initiation (Spriggs et al., 2010; de Nadal et al., 2011).
Figure 3. Regulation of protein synthesis and turn over activity. Protein names in blue were decreased, while red ones were increased. The underlined proteins were present in the enriched “translation” biological process in Figure 2B. 30S, 30S ribosomal protein; infC, translation initiation factor IF-3; EF-G, elongation factor G; EF-P, elongation factor P; RP-L, large subunit ribosomal protein; prfA, peptide chain release factor 1; CALR, calreticulin; CANX, calnexin; FKBP, FK506-binding protein.
However, the chaperones involved in protein folding were remarkably heat induced (Figure 2B), especially the HSPs (Figure 2D). In this study, 26 of the 29 heat-responsive HSPs were accumulated, suggesting that their critical roles in heat response in wheat leaves. Other chaperones, for example, calnexin and calreticulin were coincident heat accumulated. A set of proteolytic enzymes were also identified, the cysteine proteinase RD21a and serine carboxypeptidase-like clade I were coincident accumulated, while subtilisin-like protease were differentially regulated.
HRPs involved in protein metabolism suggested that wheat leaves preferentially protected existing proteins by chaperones, instead of de novo substrates synthesis in response to heat stress.
GSTs and Trxs Were Significantly Accumulated in Response to Heat Stress
The glutathione-mediated redox system was significantly enriched in HRPs (Figure 2C), suggesting a remarkable alteration in the redox status of cells under heat stress. However, DAB staining showed that the H2O2 content was slightly affected by the applied heat stress (Figure 4A). Then, the TBARS content, which predicts the peroxidation and destruction of lipids with subsequent membrane damage caused by ROS (Liu et al., 2014), showed significant accumulation (Figure 4B), indicating that membrane damage was induced by the oxidative stress caused by heat.
Figure 4. Redox regulation in wheat flag leaf in heat condition. (A) DAB staining of wheat flag leaves. (B) TBARS content of wheat flag leaves. (C) Schematic representation of redox regulation. Red represent increased, blue represent decreased, brown with italic represent differentially regulated. Map was drawn based on Mittler et al. (2004). APX, Ascorbate peroxidase; ASC, Ascorbate; CAT, Catalase; DHA, Dehydroascorbate; DHAR, Dehydroascorbate reductase; GLR, Glutaredoxin; GPX, Glutathione peroxidase; GR, Glutathione reductase; GSSG, glutathione disulphide; GST, Glutathione S-transferase; MDAR, Monodehydroascorbate reductase; Trx, Thioredoxin; Prx, Peroxiredoxin; SOD, Superoxide dismutase; (D) qRT-PCR analysis of Trxs and GSTs.
To reveal the major players that maintained redox homeostasis in wheat leaves, heat-responsive antioxidant proteins were investigated and mapped to the redox regulation system (Figure 4C). Twelve GSTs and three Trxs showed significant accumulation (Figure 2D), with increase of polyphenol oxidase and decrease of three other redox regulatory enzymes (glutaredoxin, glutathione peroxidase, and peroxiredoxin). Our proteomics data suggested that GSTs and Trxs were the major players in regulating redox alterations in heat-stressed wheat leaves.
Unexpectedly, in contrast to the accumulation of proteins, the transcripts of GSTs and Trxs in the same heat-stressed samples were reduced or remained relatively stable, respectively Figure 4D, Supplementary Table 5. These results suggested that the significant increases in GSTs and Trxs may result from post-transcriptional regulation of long-term heat stress.
Carbon Fixation and Chlorophyllide Synthesis Were Remarkably Affected by Heat Stress
In this study, HRPs in chloroplasts showed significantly decrease under heat stress (Figure 5A), which indicated that processes in chloroplast maybe repressed. Photosynthesis is an already known heat-sensitive process (Farooq et al., 2011; Prasad et al., 2011), correspondingly, 19 out of the 24 photosynthesis-related HRPs involved in chlorophyll synthesis, photosynthetic electron transport, and the Calvin cycle, showed a decrease in abundance after heat stress. Furthermore, porphyrin and chlorophyll metabolism, and carbon fixation were significantly enriched among the decreased HRPs (Figure 2C). These results indicated photosynthesis was affected through multiple processes during heat stress.
Figure 5. Regulation of photosynthesis in wheat leaf under heat stress. (A) Subcellular localization of HRPs predicated by Wolfpsort. White bars represent increased proteins, black bars represent decreased proteins. (B) Map of Calvin cycle regulation. Blue represent decreased, brown represent differentially regulated. 1,3-BPG, 1,3-Bisphospho-D-glycerate; F-6-P, D-Fructose 6-phosphate; GAP, D-Glyceraldehyde 3-phosphate; GPA, D-Glycerate 3-phosphate; PGK, phosphoglycerate kinase; RCA, Rubisco activase; rbcL, Rubisco large chain, rbcS; Rubisco small chain; RuBP, D-Ribulose 1,5-bisphosphate; S7P, D-Sedoheptulose 7-phosphate; SBPase, sedoheptulose-bisphosphatase; (C) Chlorophyll content of wheat flag leaves. Asterisk represent significantly different at p < 0.05 level, double asterisks at 0.01 level. (D) Map of chlorophyll biosynthesis regulation. Red represent increased, blue represent decreased. chlE, magnesium-protoporphyrin IX monomethyl ester (oxidative) cyclase; chlH, magnesium chelatase subunit H; chlI, magnesium chelatase subunit I; hemB, porphobilinogen synthase; hemC, hydroxymethylbilane synthase; hemE, uroporphyrinogen decarboxylase (E) qRT-PCR analysis of genes in chlorophyll synthesis and Calvin cycle.
Rubisco is crucial for carboxylation in the Calvin cycle and it is known to be inhibited by heat stress (Lee et al., 2007; Wang et al., 2015). In our results, all the four heat-responsive Rubisco large subunits and Rubisco small subunits were significantly decreased (Supplementary Table 4). In addition, chloroplast phosphoglycerate kinase and sedoheptulose-bisphosphatase were also decreased (Figure 5B), indicating that heat stress has a negative impact on the Calvin cycle.
Chlorophyll metabolism was significantly enriched in decreased HRPs (Figure 2C), indicating that the chlorophyll content of flag leaves may be reduced during heat stress. Chlorophyll quantification confirmed that the chlorophyll content decreased with heat stress (Figure 5C). Furthermore, both chlorophyll a and b were decreased (Figure 5C), indicating that their common precursor was likely reduced during heat stress. Correspondingly, five HRPs (hemE, hemC, chlE, chlH, and chlI) that were critical for the synthesis of protochlorophyllide from porphobilinogen were significantly decreased (Figure 5D). This finding indicated that the inhibition of protochlorophyllide biosynthesis was likely due to the decreased abundance of these proteins.
Further, qRT-PCR quantification showed that the mRNAs of all of the HRPs involved in chlorophyll synthesis and carbon fixation were significantly downregulated under heat stress (Figure 5E). Changes in the mRNA levels were more severe than those at the protein level, indicating that the decreases in chlorophyll synthesis and carbon fixation may be initially regulated at the transcriptional level during heat stress.
Expression Pattern of Photosynthesis-Related Genes during Continuous Heat Stress
The results above demonstrated that HRPs involved in photosynthesis were all encoded by heat sensitive genes. To illustrate their heat response patterns, transcripts analysis were conducted at wheat seedling stage and showed that most of them were inhibited by temperature rises (Figure 6A), suggesting that the process of photosynthesis was rapidly regulated at the transcriptional level following heat stress. However, the rapid downregulation of genes involved in chlorophyll precursor synthesis did not result in substantial chlorophyll loss until 4–5 days later (Figure 6B).
Figure 6. Chlorophyll content dynamics and the photosynthesis process related genes expressions under heat stress at wheat seedling stage. The Y-axis represent fold changes (heat to control) in all figures. (A) Transcripts analysis of genes involved in chlorophyll synthesis and Calvin cycle under 24, 37, 40°C for 4 h under light condition. (B) Chlorophyll content dynamics under daylight heat stress at 40°C. (C–E), time course analysis of mRNAs of chlorophyll synthesis and Calvin cycle related genes. Samples at time point of 2 and 3 d were obtained 4 h after light with heat stress.
Therefore, the time course analysis of the mRNA of these genes under severe heat stress (40°C) was performed. The results showed that the transcripts of these proteins were encoded by short-term heat-responsive genes that were downregulated within 1–2 h following heat stress, as reported by previous transcriptome analysis (Liu et al., 2015). However, in the prolonged heat treatment, these genes displayed different transcriptional profiles. HemB and HemC were recovered at the transcriptional level, indicating a heat acclimation pattern (Figure 6C), whereas the others genes were severely heat-suppressed (Figures 6D,E), with five of the nine genes displaying a further decrease in transcription with recurring heat stress, indicating the accumulated effects of heat stress (Figure 6D).
Discussion
Induced Chaperones Are Key Players in Wheat Leaf Heat Adaptation
Protein folding and ER-associated degradation system are important components of ER quality control, unfold proteins stress is induced when the misfolded proteins are over aggregated when plants are suffered from environment stress (Liu and Howell, 2010; Liu and Li, 2014; Ruggiano et al., 2014). Chaperones and cochaperones protect substrates from aggregation and aid proteins refolding. Furthermore, increases in chaperones have been proposed to enhance heat tolerance in plants and maintain yields under heat stress (Jia et al., 2008; Salvucci, 2008; Chauhan et al., 2012; Zhong et al., 2013). In our results, many chaperones together with some proteasome subunits were accumulated, suggesting that protein folding and ER-associated degradation system are induced to maintain normal plant growth and development during heat stress. Further, studies to characterize substrates of these chaperones will be critical for clarifying their specific roles. In addition, due to complexity of protein quality control and heat stress response, interactions within chaperones/cochaperones and between substrates also need more efforts to illustrate their protection networks (Wang et al., 2004; Kotak et al., 2007; Hahn et al., 2011). Among these chaperones, accumulating evidence demonstrates that HSP induction is the conserved mechanism for plant adaptation, and may be considered to be a potential target for molecular-assisted selection in thermotolerance breeding.
In contrast to the significant accumulation of chaperones, the proteolytic enzymes or components were differentially regulated, possibly due to their different roles. For example, serine carboxypeptidase-like clade I were accumulated, but the clade II was decreased. In Arabidopsis, clade I may act as acyltransferases while clade II were regarded as genuine carboxypeptidases (Fraser et al., 2005), and their roles in plant growth in heat stress need further characterizations.
GST and Trx Play a Crucial Role in Redox Regulation in Leaves under Heat Stress
As a response to the oxidative burst, several antioxidant enzymes such as CAT, SOD, and APX (Lee et al., 2007; Qin et al., 2008; Liu et al., 2014) are often induced to maintain redox homeostasis in plants. However, their abundance was not induced by heat stress in this study, as shown in previous heat stress studies (Shin et al., 2005; Lee et al., 2007; Wang et al., 2011, 2015). By contrast, GSTs and Trxs that were also involved in maintaining redox homeostasis in plants (Dixon et al., 2002; Park et al., 2009; Zhang et al., 2011; Zagorchev et al., 2013; Labrou et al., 2015), played a prominent role in redox regulation in our analysis. Accumulation of GST was thought to contribute to thermotolerance by reducing ROS production (Xu and Huang, 2008). This finding inferred an unusual redox regulatory system comprising GSTs and Trxs was activated in response to redox changes within leaves induced by heat stress.
After 3 days of heat stress, the mRNA level of GSTs and Trxs were reduced or maintained at a relatively stable level, respectively (Figure 4D). Some of these genes showed similar expression patterns to those reported in previous study (Liu et al., 2015), suggesting that the accumulation of these proteins may result from post-transcriptional regulation. The contrast between the mRNA and protein expression level of GSTs may be due to their translation efficiencies, perhaps mediated by alternative translation mechanisms such as those of HSPs (Spriggs et al., 2010). Furthermore, their interactions with HSPs (Robin et al., 2003; Montrichard et al., 2009; Raza, 2011), may confer GST and Trx proteins more thermostability (Gautier et al., 1998; Yamamoto et al., 2005, 2007; Traverso et al., 2007), and guarantee their functions in heat stress. These results suggest that GSTs and Trxs may be promising candidates for improving wheat thermotolerance.
Inhibition of Photosynthesis Following Heat Stress
Photosynthesis is the basis for plant growth and development. Previous studies have demonstrated that heat stress significantly reduced photosynthesis (Ristic et al., 2007; Fu et al., 2008; Prasad et al., 2008; Wang et al., 2011; Narayanan et al., 2015). Our proteomics data showed that most HRPs involved in chlorophyll synthesis, photosynthetic electron transport, and the Calvin cycle were heat-suppressed, providing a global perspective for the underlying mechanism of decreased photosynthesis during heat stress. Moreover, the HRPs involved in photosynthesis that were decreased in abundance showed reduced levels of transcription, indicating that these HRPs may not be the direct heat damage targets but rather may be indirectly regulated at the transcriptional level.
A potential option for improving wheat thermotolerance would be to enable heat adaptation of photosynthesis process, which we showed to suffer the accumulated effects of heat stress. A previous study showed that short-term heat stress conferred less loss of floret fertility and individual grain weight compared with long-term heat stress, which also demonstrated a cumulative heat response (Prasad and Djanaguiraman, 2014). Thus, short episodes of heat stress may result in mild damage to wheat growth, which permits potential recovery under normal conditions. Therefore, extending the heat tolerance period may be a potential breeding target to achieve improved wheat thermotolerance.
Author Contributions
YL, SX, XW, and WZ performed proteomic data analysis. YL, RL, RW, and ST performed sample collection, mRNA quantification and physiological experiments. SX, YL, QS, and SD designed the study, and organized all the data and wrote the manuscript. All authors were involved with the revision of the manuscript and approved the final manuscript.
Funding
This work was supported by the National Natural Science Foundation of China (31571756 and 31501380), grants for Starting Research of Northwest A&F University (Z111021308 and Z111021505).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank Huatao Yu, Ning Yang and Hao Yang for valuable comments to improve this manuscript; Xiaojun Peng for technical help with proteomic analysis.
Supplementary Material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2017.01086/full#supplementary-material
References
Al-Khatib, K., and Paulsen, G. M. (1984). Mode of high temperature injury to wheat during grain development. Physiol. Plant 61, 363–368. doi: 10.1111/j.1399-3054.1984.tb06341.x
Chauhan, H., Khurana, N., Nijhavan, A., Khurana, J. P., and Khurana, P. (2012). The wheat chloroplastic small heat shock protein (sHSP26) is involved in seed maturation and germination and imparts tolerance to heat stress. Plant Cell Environ. 35, 1912–1931. doi: 10.1111/j.1365-3040.2012.02525.x
Choulet, F., Alberti, A., Theil, S., Glover, N., Barbe, V., Daron, J., et al. (2014). Structural and functional partitioning of bread wheat chromosome 3B. Science 345:1249721. doi: 10.1126/science.1249721
de Nadal, E., Ammerer, G., and Posas, F. (2011). Controlling gene expression in response to stress. Nat. Rev. Genetics 12, 833–845. doi: 10.1038/nrg3055
Dixon, D. P., Lapthorn, A., and Edwards, R. (2002). Plant glutathione transferases. Genome Biol. 3, 3001. doi: 10.1186/gb-2002-3-3-reviews3004
Farooq, M., Bramley, H., Palta, J. A., and Siddique, K. H. (2011). Heat stress in wheat during reproductive and grain-filling phases. Crit. Rev. Plant Sci. 30, 491–507. doi: 10.1080/07352689.2011.615687
Fraser, C. M., Rider, L. W., and Chapple, C. (2005). An expression and bioinformatics analysis of the arabidopsis serine carboxypeptidase-like gene family. Plant Physiol. 138, 1136–1148. doi: 10.1104/pp.104.057950
Fu, J., Momčilović, I., Clemente, T. E., Nersesian, N., Trick, H. N., and Ristic, Z. (2008). Heterologous expression of a plastid EF-Tu reduces protein thermal aggregation and enhances CO2 fixation in wheat (Triticum aestivum) following heat stress. Plant Mol. Biol. 68, 277–288. doi: 10.1007/s11103-008-9369-6
Gautier, M.-F., Lullien-Pellerin, V., de Lamotte-Guéry, F., Guirao, A., and Joudrier, P. (1998). Characterization of wheat thioredoxin h cDNA and production of an active Triticum aestivum protein in Escherichia coli. Eur. J. Biochem. 252, 314–324. doi: 10.1046/j.1432-1327.1998.2520314.x
Hahn, A., Bublak, D., Schleiff, E., and Scharf, K.-D. (2011). Crosstalk between Hsp90 and Hsp70 chaperones and heat stress transcription factors in tomato. Plant Cell 23, 741–755. doi: 10.1105/tpc.110.076018
Horton, P., Park, K. J., Obayashi, T., Fujita, N., Harada, H., Adams-Collier, C. J., et al. (2007). WoLF PSORT: protein localization predictor. Nucleic Acids Res. 35, 585–587. doi: 10.1093/nar/gkm,259
IWSGC (2014). A chromosome-based draft sequence of the hexaploid bread wheat (Triticum aestivum) genome. Science 345:1251788. doi: 10.1126/science.1251788
Jia, X.-Y., Xu, C.-Y., Jing, R.-L., Li, R.-Z., Mao, X.-G., Wang, J.-P., et al. (2008). Molecular cloning and characterization of wheat calreticulin (CRT) gene involved in drought-stressed responses. J. Exp. Bot. 59, 739–751. doi: 10.1093/jxb/erm369
Kotak, S., Larkindale, J., Lee, U., von Koskull-Döring, P., Vierling, E., and Scharf, K.-D. (2007). Complexity of the heat stress response in plants. Curr. Opin. Plant Biol. 10, 310–316. doi: 10.1016/j.pbi.2007.04.011
Labrou, N. E., Papageorgiou, A. C., Pavli, O., and Flemetakis, E. (2015). Plant GSTome: structure and functional role in xenome network and plant stress response. Curr. Opin. Biotechnol. 32, 186–194. doi: 10.1016/j.copbio.2014.12.024
Law, R. D., and Crafts-Brandner, S. J. (1999). Inhibition and acclimation of photosynthesis to heat stress is closely correlated with activation of ribulose-1, 5-bisphosphate carboxylase/oxygenase. Plant Physiol. 120, 173–182. doi: 10.1104/pp.120.1.173
Lee, D. G., Ahsan, N., Lee, S. H., Kang, K. Y., Bahk, J. D., Lee, I. J., et al. (2007). A proteomic approach in analyzing heat-responsive proteins in rice leaves. Proteomics 7, 3369–3383. doi: 10.1002/pmic.200700266
Lesk, C., Rowhani, P., and Ramankutty, N. (2016). Influence of extreme weather disasters on global crop production. Nature 529, 84–87. doi: 10.1038/nature16467
Liao, J.-L., Zhou, H.-W., Zhang, H.-Y., Zhong, P.-A., and Huang, Y.-J. (2014). Comparative proteomic analysis of differentially expressed proteins in the early milky stage of rice grains during high temperature stress. J. Exp. Bot. 65, 655–671. doi: 10.1093/jxb/ert435
Liu, G.-T., Ma, L., Duan, W., Wang, B.-C., Li, J.-H., Xu, H.-G., et al. (2014). Differential proteomic analysis of grapevine leaves by iTRAQ reveals responses to heat stress and subsequent recovery. BMC Plant Biol. 14:110. doi: 10.1186/1471-2229-14-110
Liu, J.-X., and Howell, S. H. (2010). Endoplasmic reticulum protein quality control and its relationship to environmental stress responses in plants. Plant Cell 22, 2930–2942. doi: 10.1105/tpc.110.078154
Liu, Y., and Li, J. (2014). Endoplasmic reticulum-mediated protein quality control in Arabidopsis. Front. Plant Sci. 5:162. doi: 10.3389/fpls.2014.00162
Liu, Z., Xin, M., Qin, J., Peng, H., Ni, Z., Yao, Y., et al. (2015). Temporal transcriptome profiling reveals expression partitioning of homeologous genes contributing to heat and drought acclimation in wheat (Triticum aestivum L.). BMC Plant Biol. 15:5118. doi: 10.1186/s12870-015-0511-8
Lobell, D. B., and Tebaldi, C. (2014). Getting caught with our plants down: the risks of a global crop yield slowdown from climate trends in the next two decades. Environ. Res. Lett. 9:074003. doi: 10.1088/1748-9326/9/7/074003
Majoul, T., Bancel, E., Triboi, E., Ben Hamida, J., and Branlard, G. (2003). Proteomic analysis of the effect of heat stress on hexaploid wheat grain: characterization of heat-responsive proteins from total endosperm. Proteomics 3, 175–183. doi: 10.1002/pmic.200390026
Majoul, T., Bancel, E., Triboi, E., Ben Hamida, J., and Branlard, G. (2004). Proteomic analysis of the effect of heat stress on hexaploid wheat grain: characterization of heat-responsive proteins from non-prolamins fraction. Proteomics 4, 505–513. doi: 10.1002/pmic.200300570
Mangelsen, E., Kilian, J., Harter, K., Jansson, C., Wanke, D., and Sundberg, E. (2011). Transcriptome analysis of high-temperature stress in developing barley caryopses: early stress responses and effects on storage compound biosynthesis. Mol. Plant 4, 97–115. doi: 10.1093/mp/ssq058
Mittler, R., Vanderauwera, S., Gollery, M., and Van Breusegem, F. (2004). Reactive oxygen gene network of plants. Trends Plant Sci. 9, 490–498. doi: 10.1016/j.tplants.2004.08.009
Montrichard, F., Alkhalfioui, F., Yano, H., Vensel, W. H., Hurkman, W. J., and Buchanan, B. B. (2009). Thioredoxin targets in plants: the first 30 years. J. Proteomics 72, 452–474. doi: 10.1016/j.jprot.2008.12.002
Morimoto, R. I. (1998). Regulation of the heat shock transcriptional response: cross talk between a family of heat shock factors, molecular chaperones, and negative regulators. Genes Dev. 12, 3788–3796. doi: 10.1101/gad.12.24.3788
Narayanan, S., Prasad, P., Fritz, A., Boyle, D., and Gill, B. (2015). Impact of high night-time and high daytime temperature stress on winter wheat. J. Agron. Crop Sci. 201, 206–218. doi: 10.1111/jac.12101
Nguyen, H. T., Joshi, C. P., Klueva, N., Weng, J., Hendershot, K. L., and Blum, A. (1994). The heat-shock response and expression of heat-shock proteins in wheat under diurnal heat stress and field conditions. Funct. Plant Biol. 21, 857–867. doi: 10.1071/pp9940857
Park, S. K., Jung, Y. J., Lee, J. R., Lee, Y. M., Jang, H. H., Lee, S. S., et al. (2009). Heat-shock and redox-dependent functional switching of an h-type Arabidopsis thioredoxin from a disulfide reductase to a molecular chaperone. Plant Physiol. 150, 552–561. doi: 10.1104/pp.109.135426
Pradhan, G. P., and Prasad, P. V. (2015). Evaluation of wheat chromosome translocation lines for high temperature stress tolerance at grain filling stage. PLoS ONE 10:e0116620. doi: 10.1371/journal.pone.0116620
Prasad, P., Pisipati, S., Ristic, Z., Bukovnik, U., and Fritz, A. (2008). Impact of nighttime temperature on physiology and growth of spring wheat. Crop Sci. 48, 2372–2380. doi: 10.2135/cropsci2007.12.0717
Prasad, P. V., and Djanaguiraman, M. (2014). Response of floret fertility and individual grain weight of wheat to high temperature stress: sensitive stages and thresholds for temperature and duration. Funct. Plant Biol. 41, 1261–1269. doi: 10.1071/FP14061
Prasad, P. V. V., Pisipati, S. R., Momčilović, I., and Ristic, Z. (2011). Independent and combined effects of high temperature and drought stress during grain filling on plant yield and chloroplast EF-Tu expression in spring wheat. J. Agron. Crop Sci. 197, 430–441. doi: 10.1111/j.1439-037X.2011.00477.x
Qin, D., Wu, H., Peng, H., Yao, Y., Ni, Z., Li, Z., et al. (2008). Heat stress-responsive transcriptome analysis in heat susceptible and tolerant wheat (Triticum aestivum L.) by using wheat genome array. BMC Genomics 9:432. doi: 10.1186/1471-2164-9-432
Raza, H. (2011). Dual localization of glutathione S-transferase in the cytosol and mitochondria: implications in oxidative stress, toxicity and disease. FEBS J. 278, 4243–4251. doi: 10.1111/j.1742-4658.2011.08358.x
Ristic, Z., Bukovnik, U., and Prasad, P. V. (2007). Correlation between heat stability of thylakoid membranes and loss of chlorophyll in winter wheat under heat stress. Crop Sci. 47, 2067–2073. doi: 10.2135/cropsci2006.10.0674
Robin, M.-A., Prabu, S. K., Raza, H., Anandatheerthavarada, H. K., and Avadhani, N. G. (2003). Phosphorylation enhances mitochondrial targeting of GSTA4-4 through increased affinity for binding to cytoplasmic Hsp70. J. Biol. Chem. 278, 18960–18970. doi: 10.1074/jbc.M301807200
Rollins, J., Habte, E., Templer, S., Colby, T., Schmidt, J., and Von Korff, M. (2013). Leaf proteome alterations in the context of physiological and morphological responses to drought and heat stress in barley (Hordeum vulgare L.). J. Exp. Bot. 64, 3201–3212. doi: 10.1093/jxb/ert158
Ruggiano, A., Foresti, O., and Carvalho, P. (2014). ER-associated degradation: protein quality control and beyond. J. Cell Biol. 204, 869. doi: 10.1083/jcb.201312042
Salvucci, M. E. (2008). Association of Rubisco activase with chaperonin-60β: a possible mechanism for protecting photosynthesis during heat stress. J. Exp. Bot. 59, 1923–1933. doi: 10.1093/jxb/erm343
Sarkar, N. K., Kim, Y.-K., and Grover, A. (2014). Coexpression network analysis associated with call of rice seedlings for encountering heat stress. Plant Mol. Biol. 84, 125–143. doi: 10.1007/s11103-013-0123-3
Sharkey, T. D. (2005). Effects of moderate heat stress on photosynthesis: importance of thylakoid reactions, rubisco deactivation, reactive oxygen species, and thermotolerance provided by isoprene. Plant Cell Environ. 28, 269–277. doi: 10.1111/j.1365-3040.2005.01324.x
Shin, S.-Y., Lee, H.-S., Kwon, S.-Y., Kwon, S.-T., and Kwak, S.-S. (2005). Molecular characterization of a cDNA encoding copper/zinc superoxide dismutase from cultured cells of Manihot esculenta. Plant Physiol. Biochem. 43, 55–60. doi: 10.1016/j.plaphy.2004.12.005
Spriggs, K. A., Bushell, M., and Willis, A. E. (2010). Translational regulation of gene expression during conditions of cell stress. Mol. Cell 40, 228–237. doi: 10.1016/j.molcel.2010.09.028
Suzuki, N., Koussevitzky, S., Mittler, R., and Miller, G. (2012). ROS and redox signalling in the response of plants to abiotic stress. Plant Cell Environ. 35, 259–270. doi: 10.1111/j.1365-3040.2011.02336.x
Tack, J., Barkley, A., and Nalley, L. L. (2015). Effect of warming temperatures on US wheat yields. Proc. Natl. Acad. Sci. U.S.A. 112, 6931–6936. doi: 10.1073/pnas.1415181112
Tian, F., Gong, J., Zhang, J., Zhang, M., Wang, G., Li, A., et al. (2013). Enhanced stability of thylakoid membrane proteins and antioxidant competence contribute to drought stress resistance in the tasg1 wheat stay-green mutant. J. Exp. Bot. 64, 1509–1520. doi: 10.1093/jxb/ert004
Traverso, J. A., Vignols, F., Cazalis, R., Pulido, A., Sahrawy, M., Cejudo, F. J., et al. (2007). PsTRXh1 and PsTRXh2 are both Pea h-type thioredoxins with antagonistic behavior in redox imbalances. Plant Physiol. 143, 300–311. doi: 10.1104/pp.106.089524
Trevisan, S., Manoli, A., Ravazzolo, L., Botton, A., Pivato, M., Masi, A., et al. (2015). Nitrate sensing by the maize root apex transition zone: a merged transcriptomic and proteomic survey. J. Exp. Bot. 66, 3699–3715. doi: 10.1093/jxb/erv165
ur Rehman, A., Habib, I., Ahmad, N., Hussain, M., Khan, M. A., Farooq, J., et al. (2009). Screening wheat germplasm for heat tolerance at terminal growth stage. Plant Omics 2, 9–19.
Vizcaino, J. A., Csordas, A., Del-Toro, N., Dianes, J. A., Griss, J., Lavidas, I., et al. (2016). 2016 update of the PRIDE database and its related tools. Nucleic Acids Res. 44, 11033. doi: 10.1093/nar/gkw880
Vizcaino, J. A., Deutsch, E. W., Wang, R., Csordas, A., Reisinger, F., Rios, D., et al. (2014). ProteomeXchange provides globally coordinated proteomics data submission and dissemination. Nat. Biotechnol. 32, 223–226. doi: 10.1038/nbt.2839
Wang, W., Vinocur, B., Shoseyov, O., and Altman, A. (2004). Role of plant heat-shock proteins and molecular chaperones in the abiotic stress response. Trends Plant Sci. 9, 244–252. doi: 10.1016/j.tplants.2004.03.006
Wang, X., Cai, J., Jiang, D., Liu, F., Dai, T., and Cao, W. (2011). Pre-anthesis high-temperature acclimation alleviates damage to the flag leaf caused by post-anthesis heat stress in wheat. J. Plant Physiol. 168, 585–593. doi: 10.1016/j.jplph.2010.09.016
Wang, X., Dinler, B. S., Vignjevic, M., Jacobsen, S., and Wollenweber, B. (2015). Physiological and proteome studies of responses to heat stress during grain filling in contrasting wheat cultivars. Plant Sci. 230, 33–50. doi: 10.1016/j.plantsci.2014.10.009
Xu, C., and Huang, B. (2008). Root proteomic responses to heat stress in two Agrostis grass species contrasting in heat tolerance. J. Exp. Bot. 59, 4183–4194. doi: 10.1093/jxb/ern258
Yamamoto, K., Fujii, H., Aso, Y., Banno, Y., and Koga, K. (2007). Expression and characterization of a sigma-class glutathione S-transferase of the fall webworm, Hyphantria cunea. Biosci. Biotechnol. Biochem. 71, 553–560. doi: 10.1271/bbb.60592
Yamamoto, K., Zhang, P., Miake, F., Kashige, N., Aso, Y., Banno, Y., et al. (2005). Cloning, expression and characterization of theta-class glutathione S-transferase from the silkworm, Bombyx mori. Comp. Biochem. Physiol. 141, 340–346. doi: 10.1016/j.cbpc.2005.04.012
Yentrapalli, R., Azimzadeh, O., Barjaktarovic, Z., Sarioglu, H., Wojcik, A., Harms-Ringdahl, M., et al. (2013). Quantitative proteomic analysis reveals induction of premature senescence in human umbilical vein endothelial cells exposed to chronic low-dose rate gamma radiation. Proteomics 13, 1096–1107. doi: 10.1002/pmic.201200463
Zagorchev, L., Seal, C. E., Kranner, I., and Odjakova, M. (2013). A central role for thiols in plant tolerance to abiotic stress. Int. J. Mol. Sci. 14, 7405–7432. doi: 10.3390/ijms14047405
Zhang, C.-J., Zhao, B.-C., Ge, W.-N., Zhang, Y.-F., Song, Y., Sun, D.-Y., et al. (2011). An apoplastic h-type thioredoxin is involved in the stress response through regulation of the apoplastic reactive oxygen species in rice. Plant Physiol. 157, 1884–1899. doi: 10.1104/pp.111.182808
Zhao, H., Dai, T., Jing, Q., Jiang, D., and Cao, W. (2007). Leaf senescence and grain filling affected by post-anthesis high temperatures in two different wheat cultivars. Plant Growth Regul. 51, 149–158. doi: 10.1007/s10725-006-9157-8
Keywords: wheat, iTRAQ, heat, photosynthesis, redox regulation
Citation: Lu Y, Li R, Wang R, Wang X, Zheng W, Sun Q, Tong S, Dai S and Xu S (2017) Comparative Proteomic Analysis of Flag Leaves Reveals New Insight into Wheat Heat Adaptation. Front. Plant Sci. 8:1086. doi: 10.3389/fpls.2017.01086
Received: 06 April 2017; Accepted: 06 June 2017;
Published: 20 June 2017.
Edited by:
Ive De Smet, Flanders Institute for Biotechnology, BelgiumReviewed by:
Lyudmila Petrova Simova-Stoilova, Institute of Plant Physiology and Genetics (BAS), BulgariaTong Zhang, University of Florida, United States
Copyright © 2017 Lu, Li, Wang, Wang, Zheng, Sun, Tong, Dai and Xu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shengbao Xu, xushb@nwsuaf.edu.cn