Impact Factor 4.106 | CiteScore 4.47
More on impact ›

Original Research ARTICLE

Front. Plant Sci., 26 August 2015 |

Comparative proteomic analysis of a membrane-enriched fraction from flag leaves reveals responses to chemical hybridization agent SQ-1 in wheat

Qilu Song, Shuping Wang, Gaisheng Zhang*, Ying Li, Zheng Li, Jialin Guo, Na Niu, Junwei Wang and Shoucai Ma
  • College of Agronomy, Northwest Agriculture and Forestry University, National Yangling Agricultural Biotechnology and Breeding Center, Yangling Branch of State Wheat Improvement Centre, Wheat Breeding Engineering Research Center, Ministry of Education, Key Laboratory of Crop Heterosis of Shaanxi Province, Yangling, China

The induction of wheat male fertile lines by using the chemical hybridizing agent SQ-1 (CHA-SQ-1) is an effective approach in the utilization of heterosis; however, the molecular basis of male fertility remains unknown. Wheat flag leaves are the initial receptors of CHA-SQ-1 and their membrane structure plays a vital role in response to CHA-SQ-1 stress. To investigate the response of wheat flag leaves to CHA-SQ-1 stress, we compared their quantitative proteomic profiles in the absence and presence of CHA-SQ-1. Our results indicated that wheat flag leaves suffered oxidative stress during CHA-SQ-1 treatments. Leaf O2-, H2O2, and malonaldehyde levels were significantly increased within 10 h after CHA-SQ-1 treatment, while the activities of major antioxidant enzymes such as superoxide dismutase, catalase, and guaiacol peroxidase were significantly reduced. Proteome profiles of membrane-enriched fraction showed a change in the abundance of a battery of membrane proteins involved in multiple biological processes. These variable proteins mainly impaired photosynthesis, ATP synthesis protein mechanisms and were involved in the response to stress. These results provide an explanation of the relationships between membrane proteomes and anther abortion and the practical application of CHA for hybrid breeding.


Membranes are highly organized structures and important components of plant cells. The membrane system of plant cells includes the plasma membrane (PM) and organelle membranes (nuclear membrane, endoplasmic reticulum (ER) membrane, Golgi membrane, mitochondrial membrane, chloroplast membrane, and lysosomal membrane). Membranes not only form cellular compartments for performing multiple differential metabolic processes and maintaining organelle homeostasis, but they also play a critical role in the exchange of substances and signals (Jaiswal et al., 2012). Lipids and proteins are major components of membranes. Membrane proteins take part in multiple biological reactions such as metabolite and ion exchange, signal transduction, biosynthesis, photosynthesis, and energy generation (Kota and Goshe, 2011; Liu et al., 2011; Takahashi et al., 2013).

Plants may suffer various biotic and abiotic stresses during growth, which can cause damage to cellular functions. However, within a certain threshold, plants have a battery of protective mechanisms to maintain the stability of the normal mechanical processes of cellular homeostasis. These protective mechanisms include changes in lipid composition and regulation of protein expression. Some reports show that changes in the ratio of lipid composition and an increase in the level of unsaturation of membrane lipids in some species under various stressors are effective in maintaining a stable membrane conformation (Welti et al., 2002; Gigon et al., 2004; López-Pérez et al., 2009). Membrane proteins play a vital role in plant stress response. Nouri and Komatsu (2010) suggested that transcription and expression of H+-ATPase in Arabidopsis thaliana are increased, which then accelerates the transport of ions between the inner and outer PM to maintain cellular homeostasis during hyperosmotic stress. Some studies have indicated that proteins of photosystem II (PS II), including chlorophyll a/b binding protein, ATP generating proteins, and NAD(P)H-quinone oxidoreductase, are up-regulated in Aeluropus lagopoides, barley, and mangrove plants under salt stress, which can facilitate photon capture, provide sufficient energy, detoxify quinines and maintain normal photosynthesis process (Sobhanian et al., 2010; Rasoulnia et al., 2011; Wang et al., 2013). Tan et al. (2012) suggest the existence of changes in mitochondrial membrane proteomes in response to cold and chemical stresses for normal mitochondrial function. Compared with WT rice, salt tolerance was significantly enhanced in SaVHAc1 (encode c1 subunit of vacuole H+-ATPase)-expressing rice (Baisakh et al., 2012). Some reports show that the ER membrane protein bZIP is released into the Golgi and is then degraded when plants suffer stress; subsequently, degradation products move into cell nuclei and eventually lead to the up-regulation of stress proteins encoded by nuclear genes (Liu and Howell, 2010; Zhang and Wang, 2012; Howell, 2013). Nowadays, researchers increasingly focus on the study of membrane proteins in plants, including phosphoproteomics (Nuhse et al., 2007), regulation of related metabolism in C4 plant (Manandhar-Shrestha et al., 2013), response to stress (Nouri and Komatsu, 2010; Hopff et al., 2013) and signaling (Nie et al., 2015); meanwhile, the technology of mass spectrometry offers convenience in the identification and quantification of membrane proteins.

Wheat (Triticum aestivum L.) is a valuable agricultural crop and an important food source for humans. A continual increase in wheat consumption has led to a demand for greater wheat yields to guarantee world food security (Curtis and Halford, 2014). The utilization of heterosis in wheat is still limited due to its complex hereditary basis and self-pollination; however, the chemical hybridizing agent SQ-1 (CHA-SQ-1) offers a new approach to the application of heterosis. Normal fertile male wheat can be made sterile after the spraying of appropriate doses of CHA-SQ-1 at a specific developmental period. The proportion of male sterile wheat can reach 95–100%, and the outcrossing rate can exceed 85%. However, how does CHA-SQ-1 cause male sterility in wheat? In previous studies, attention has been focused more on the anther and pollen in CHA-SQ-1-induced male sterility (Ba et al., 2013, 2014a,b; Wang et al., 2015; Zhu et al., 2015a). Despite some progress in reactive oxygen metabolism (Ba et al., 2013), aliphatic metabolism (Ba et al., 2014a), DNA methylation (Ba et al., 2014b), cell morphology (Wang et al., 2015), and transcriptome (Zhu et al., 2015a), there is a dearth of information on the interceptive mechanism of CHA-SQ-1-induced male sterility. More importantly, the involvement of leaves as the recipient tissues has not yet been reported.

Flag leaves provide sufficient energy and sucrose for meiosis to take place in the pollen mother cells and to allow pollen development and the accumulation of pollen starch during the reproductive stage (Liu et al., 2006; Zhang et al., 2010). SQ-1 is a pyridazine chemical hybridizing agent which can induce male sterility in wheat. During the induction of wheat male sterility, the flag leaves are the initial receptor of CHA-SQ-1, and are stressed by its application. However, after spraying CHA-SQ-1 on flag leaves, the dynamic characteristics of the membrane proteins of the flag leaves remain unknown and the relationship between the dynamic characteristics of membrane proteins and cellular metabolic processes also remain unclear. In order to explore how the membrane proteomes responses to this abiotic stress during CHA-SQ-1-induced male sterility, we profiled the dynamic characteristics of proteomes for flag leaf membrane-enriched fraction in the presence and absence of CHA-SQ-1.

Materials and Methods

Plant Materials

The wheat cultivar cv. Xinong 1376 was grown in an experimental field of the Northwest Agriculture and Forestry University, Yangling, China (108° E, 34° 15′ N). The wheat flag leaves were treated with CHA-SQ-1 (5 kg/hectare) when wheat reached the Feeks’ 8.5 stage (the internal morphology of wheat development where the connectivum is formed in female and male primordia). Wheat flag leaves treated with water were used as a control. Flag leaves were collected from control and treated with CHA-SQ-1 (at 2, 4, 6, 10, and 24 h after treatment, respectively) and analyzed for physiological indices. For proteomics, flag leaves treated with CHA-SQ-1 after 2 and 6 h were collected, respectively. Three independent experiments were performed as biological replicates for all experiments.

Determination of O2- Formation Rate and H2O2 Content

Determination of O2-, and H2O2 content were performed according to Ba et al. (2013). Estimation of Superoxide anion, flag leaves (0.5 g) were homogenized in 65 mM phosphate buffer (pH 7.8) and centrifuged at 4°C, 5000 × g for 10 min. The supernatants were incubated in 65 mM phosphate buffer (pH 7.8) and 10 mM hydroxylamine chlorhydrate at 25°C for 20 min, and then added 17 mM sulfanilamide and 7 mM α-naphthylamine for another 20 min. The absorbance at 530 nm was measured with a Nicolet Evolution 300 spectrophotometer (Thermo, USA) and the formation rate of O2- was calculated from a standard curve of NaNO2.

For determination of H2O2, flag leaves (0.5 g) were homogenized in cold (-20°C) acetone and centrifuged at 4°C, 3000 × g for 10 min. The supernatant was added to ammonia and 95% (v/v) hydrochloric acid containing 20% (v/v) titanium tetrachloride. After centrifugation at 10 000 × g for 10 min, the sediment was repeatedly washed with cold (-20°C) acetone and centrifuged at 14 000 × g for 4°C and then dissolved in 1 M H2SO4. The absorbance was measured at 410 nm and the content of H2O2 in the leaves was calculated with an H2O2 solution-derived standard curve.

Determination of Activities of Antioxidant Enzymes

The activities of SOD, CAT, and POD were performed according to Ba et al. (2013). To extract antioxidant enzymes, 0.5 g fresh flag leaves were ground in 50 mM cool phosphate buffer [containing 1% (w/v) polyvinylpyrrolidone, pH 7.0] and centrifuged at 4°C, 15 000 × g for 20 min. The supernatant was used for assays of enzyme activity.

The activity of SOD was determined by measuring its ability to inhibit the photoreduction of nitro blue tetrazolium (NBT). The reaction solution contained 50 μM NBT, 1.3 μM riboflavin, 13 mM methionine, 75 nM EDTA, 50 mM phosphate buffer (pH 7.8), and enzyme extract. The photo-induced reaction was performed under a light bank at 78 μM m-2 s-1 for 15 min. The absorbance of the irradiated and non-irradiated solution at 560 nm was determined with a Nicolet Evolution 300 spectrophotometer (Thermo, USA). One unit of SOD activity was defined as the amount of enzyme that would inhibit 50% of NBT photo reduction.

The CAT reaction solution contained 50 mM phosphate buffer (pH 7.0), 15 mM H2O2, and enzyme extract. Reaction was initiated by adding enzyme extract. Changes in absorbance of the reaction solution at 240 nm were read every 20 s. One unit CAT activity was defined as an absorbance change of 0.01 units per min.

The POD reaction solution contained 50 mM sodium acetate buffer (pH 5.0), 20 mM guaiacol, 40 mM H2O2, and enzyme extract. Changes in absorbance of the reaction solution at 470 nm were determined every 20 s. One unit POD activity was defined as an absorbance change of 0.01 units per min.

Determination of MDA Content

The level of lipid peroxidation in samples was determined by estimating the MDA content according to Ba et al. (2013). 0.5 g leaf sample was homogenized in 20% (v/v) trichloroacetic acid and 0.5% (v/v) thiobarbituric acid, and centrifuged at 10 000 × g for 10 min. The amount of MDA in the supernatant was estimated by the thiobarbituric acid reaction.

Measurement of Photosynthesis

Net photosynthetic rates (Pn) of wheat flag leaves were measured with a portable photosynthesis measurement system (Li-6400, Li-Cor, USA) at a photosynthetic photo flux density (PPDF) of 1200 μmol m-2 s-1 (provided by a red–blue LED light source) and an ambient CO2 concentration of 400 μmol mol-1.

Measurement of Total Soluble Sugars and Starch Content

Flag leaves were ground in liquid nitrogen and homogenized with 80% ethanol. After centrifugation at 16 000 × g for 10 min, the supernatant was removed to a fresh tube and the pellet was extracted another two times. Starch in the pellet was hydrolyzed with 30% perchloric acid overnight at room temperature, then incubated at 60°C for 10 min and centrifuged at 16 000 × g for 10 min. Total soluble sugar and starch content were determined using the anthrone method as described by Zhu et al. (2015b).

Isolation of Membrane-Enriched Fraction

The membrane fraction was isolated as previously described (Martinec et al., 2000) with few modifications. Approximately, 50 g of tissue was chopped and ground in a Waring blender. The tissue was homogenized in 500 ml homogenizing solution [500 mM sucrose, 50 mM Tris-MES, 50 mM EDTA-Na2, 20 mM NaF, 10% glycerol (v/v), 10 mM ascorbic acid, 0.6% PVP (w/v), 1 mM PMSF, 0.5% BSA (w/v), 1 mM DTT; pH 7.8]. The homogenates were filtered through four layers of cheese-cloth and centrifuged at 6000 × g for 15 min. Membrane vesicles were pelleted from the resulting supernatant by centrifugation at 150 000 × g for 45 min. The resulting pellet containing the membrane fraction was suspended in suspension buffer (50 mM PBS, 330 mM sucrose, 10 mM NaF, and 2 mM DTT; pH 7.8).

ATPase Activity Measurement

Enrichment assessment of various subcellular membranous components was performed by measurement of membrane-specific H+-ATPase activity compared to total ATPase activity. The hydrolytic activity of H+-ATPase was determined according to the procedures of Yan et al. (2002) and Shen et al. (2006). Hydrolytic ATPase activity was determined in 0.5 ml Tris/MES buffer containing 50 mM Tris-MES, 1 mM Na2MoO4, 5 mM MgSO4, 50 mM KCl, 0.02% Brij 58 (w/v), 5 mM ATP-Na2, pH 6.5 for P-ATPase activity or 8.0 for V-ATPase activity, and F-ATPase activities. Subsequently 1 mM NaN3, 50 mM KNO3, and 0.5 mM Na3VO4 were contained respectively for determination of F-ATPase activity, V-ATPase activity, and P-ATPase activity. The reaction was initiated by the addition of the membrane-enriched fraction (equivalent to 3 μg of membrane protein), proceeded for 30 min at 30°C, and stopped with 0.5 ml of stopping reagent (10% SDS). After 2 min, color reagent containing 2% concentrated H2SO4 (v/v) and 0.5% (NH4)2MoO4 (w/v) was added, followed immediately by 1 ml Millipore-Q water and 10 μl of 10% (w/v) ascorbic acid. Color development was completed after 30 min. To select the most suitable wavelength for the activity of H+-ATPase, we employed the full wavelength range for the first optimization using a UV-spectrophotometer. Then, ATPase activity was calculated as phosphate liberated in excess of the boiled membrane fraction control.

2-DE and Gel Image Analysis

2D-PAGE was performed following the method of Ye et al. (2013) with minor modifications. The fractions’ membrane proteins were solubilized in lysis solution containing 7 M urea, 2 M thiourea, 4% (w/v) CHAPS, 65 mM DTT, 0.5% (v/v) Bio-Lyte (pH 4–7), and 0.001% (w/v) bromphenol blue. The protein concentration was determined by using a standard curve (R2 = 0.998) of Bio-Rad Protein Assay Kit II (Bio-Rad, USA), according to the manufacturer’s instructions. About 520 μg of membrane protein was loaded onto a commercially available precast IPG strip (Bio-Rad, USA) with a 17 cm linear pH 4–7 gradient and actively rehydrated at 50 V for 12 h at 20°C. Subsequently, focusing was performed on the IPGphor apparatus (PROTEAN IEF Cell; Bio-Rad, USA) under the following conditions: 250 V for 1 h, 500 V for 1 h, 1000 V for 1 h, 400 V for 1 h, 8000 V for 4 h, and 8000 V to achieve 80000 V-h. Prior to SDS-PAGE, the strips were equilibrated for 15 min in 10 ml of reducing equilibration buffer [6 M urea, 2% SDS (w/v), 0.375 M Tris-HCl at pH 8.8, 20% glycerol (v/v), 2% DTT (w/v)] and then for another 15 min in alkylating equilibration buffer containing 2.5% (w/v) iodoacetamide instead of 2% DTT. The strips were placed on the top of vertical 11% SDS-PAGE. Electrophoresis was carried out at 15°C and 10 mA for 1 h and then at 20 mA until the dye front reached the bottom of the gel. The gels were then stained with Coomassie Brilliant Blue (CBB) G250. Each sample was run in three independent experiments (biological replicates).

Gels were visualized using a PowerLook 2100XL scanner (UMAX, Taiwan, China) at a resolution of 600 dpi. Images were analyzed using the analytical software PDQuest 2-DE 8.0.1 (Bio-Rad, USA) for spot detection, gel matching, and statistical analysis of spots. The selection of protein spots of interest for analysis by MS was based on a fold change ≥1.5 (p < 0.05).

In-Gel Digestion and MALDI-TOF/TOF MS Analysis

The protein spots were excised manually, washed twice with Millipore-Q water, and destained with fresh solution (50% ACN, 40 mM NH4HCO3). The reaction was stopped with Millipore-Q water when the blue color disappeared after about 2 min. Next, 200 μl of 40 mM ammonium bicarbonate (NH4HCO3) was added to cover the gel for 2 min, followed by repeated dehydration with changes of 150 μl 100% acetonitrile (ACN) until the gel pieces turned opaque white; the gel was then dried in a vacuum centrifuge. Enzymatic digestion was performed with 10 μg/ml trypsin (enzyme/sample ratio 1:20) during incubation for 45 min in an ice bath, and the supernatant was subsequently removed. Then, 40 μl 40 mM of NH4HCO3 solution containing 10% ACN was added to cover the gel, followed by incubation for 16 h at 37°C. Following enzymatic digestion, the supernatant was collected and the resultant peptides were extracted during 20 min incubation in 20 μl extraction solution containing 0.1% trifluoroacetic acid (TFA) and 60% ACN. The supernatant was then collected and the peptides were re-dissolved in 20 μl extraction solution (50% ACN, 0.1% TFA) for 10 min. Subsequently, the supernatant was collected and mixed with the previous two supernatants. Mass spectra were collected using a 5800 MALDI Time of Flight (TOF)/TOFTM analyzer with α-cyano-4-hydroxycinnamic acid (CHCA) as matrix and analyzed using TOF/TOFTM Series ExplorerTM Software V4.1.0 (AB Sciex, Foster City, CA, USA).

Searches for MS/MS mass spectra were performed using the Mascot search engine1 against the NCBInr and Swissprot databases with a taxonomy parameter set to green plants. For the database search, parameters such as one missed cleavage site by trypsin, peptide tolerance of 100 ppm, and an MS/MS tolerance of 0.3 Da, peptide charge of 1+, carbamidomethylation of cysteine, and oxidation of methionine as fixed and variable modifications were used. The MASCOT score, the number of peptide matches, sequence coverage, pI, and molecular weight were used to evaluate the database search results. Sequences of proteins identified as unknown, hypothetical or proteins with an uncharacterized function were used as queries for searching their homologues with BLASTP algorithm.

Bioinformatic Analysis

The prediction of transmembrane domains (TMDs) of the identified proteins was carried out using TMpred2. Grand Average of Hydropathicity (GRAVY) value for each protein was calculated using the Protein GRAVY tool3. Hierarchical cluster analysis and heat-map were performed using MultipleExperiment Viewer 4.9 software based on the Log 2-transformed fold change. All identified proteins were blasted against the Arabidopsis thaliana TAIR10 (The Arabidopsis Information Resource) protein database4 with the intention of obtaining annotated protein entries for protein–protein interaction network (PPI) tools. Results with the highest score and lowest E value were considered as relevant for each identified protein. A PPI was constructed with the online analysis tool STRING 9.1 and biological processes and molecular functions were predicted by BiNGO 3.0.2, a plugin for Cytoscape. Subcellular localization of the identified proteins was analyzed using AT_CHLORO database5, TAIR10 databas4 and SUBA3 database6.

Statistical Analysis

The results presented are the means of three independent experiments. Sample variability is given as the standard deviation of the mean. The significance of differences between control and treatment mean values was determined by Student’s t-test, at the 0.01 significance level, and where applicable at the 0.05 significance level.


Physiological Changes in Wheat Flag Leaves Under CHA-SQ-1 Treatment Conditions

O2- and H2O2, as the major members of the ROS family, play vital roles in oxidative damage to plants. In order to further explore whether the generation and accumulation of ROS caused oxidative damage in CHA-SQ-1 treated wheat flag leaves, ROS content and antioxidant enzyme activities were measured (Figure 1). CHA-SQ-1-treated flag leaves had significantly higher levels of O2-, the precursor of most ROSs, than was found in control leaves from 2 to 10 h and contents increased by 1.64-fold at 2 h, 1.88-fold at 4 h, 1.74-fold at 6 h, and 1.09-fold at 10 h, respectively (Figure 1A), thus first showing an increase and then a decrease with treatment time and eventually returning to normal levels. Meanwhile, excess O2- was catalyzed to form H2O2. The levels of H2O2 in leaves treated at 2, 4, 6, and 10 h after CHA-SQ-1 were also increased by 1.25-fold, 1.12-fold, 1.14-fold, and 1.07-fold, respectively, as compared to the control plants (Figure 1A). Meanwhile, the activities of SOD, POD, and CAT were measured, and the results showed that activities of all three antioxidant enzymes in wheat flag leaves declined significantly after CHA-SQ-1 treatment when excess ROS was generated (Figure 1B) and this further interfered with the oxidative/antioxidative balance. More importantly, excessive ROS can oxidize membrane lipids and generate MDA, which can aggravate damage to membrane structure. Not surprisingly, within 24 h after CHA-SQ-1 treatments, the levels of MDA were increased by 1.25-fold at 2 h, 1.13-fold at 4 h, 1.13-fold at 6 h, 1.07-fold at 10 h, and 1.09-fold at 24 h over control leaves, respectively (Figure 1A).


FIGURE 1. Assessment of ROS content and oxidative stress analysis in wheat flag leaves under chemical hybridizing agent (CHA)-SQ-1 treatment conditions. (A) The generation rate of O2- and content of H2O2 in control wheat flag leaves and CHA-SQ-1 treated flag leaves; content of malonaldehyde (MDA) in control and CHA-SQ-1 treated flag leaves. (B) Activities of superoxide dismutase (SOD), guaiacol peroxidase (POD) and catalase (CAT) in control and CHA-SQ-1 treated flag leaves. Data are means ± SD of three independent experiments (biological replicates). The significant of differences was assessed by Student’s t-test (P < 0.05, ∗∗P < 0.01).

To determine the effects of CHA-SQ-1 treatment on photosynthesis of flag leaves, net photosynthesis rate (Pn) was measured at different times after CHA-SQ-1 treatment. Meanwhile, contents of soluble sugar and starch were also determined. Results showed that the Pn of the wheat flag leaves decreased significantly from 20.57 μmol m-2 s-1 (control plants) to 3.11 (at 10 h after CHA-SQ-1 treatment) and began to increased slightly at 24 h after treatment (Supplementary Figure S1). Contents of soluble sugar and starch were decreased by 1.73- to 4.19-fold and 3.26- to 8.01-fold, respectively (Supplementary Figure S2).

Enrichment Assessment of Various Subcellular Membranous Components

Enrichment of the membrane microsomal fraction of flag leaves was accomplished by differential centrifugation. It was systematically assessed for the enrichment of various subcellular membranous components using standard marker enzyme assays viz., vanadate-, azide-, nitrate-sensitive ATPase PM, mitochondrial membrane, tonoplast, and Golgi membranes, respectively. The results showed that a maximum absorption at 820 nm (Figure 2A), which means 820 nm was selected as the most adequate wavelength for testing all activity of H+-ATPase. Accordingly, the relative changes in percent inhibition of ATPase activities associated with mitochondria, tonoplast, and PM were 1.60-, 1.63-, and 1.73-fold, respectively in microsomal fractions compared to the crude homogenate (Figure 2B). The higher level activities of the marker enzymes in the membrane fraction indicate the enrichment of various subcellular membranous components.


FIGURE 2. Enrichment assessment of various subcellular membranous components. (A) Wavelength scanning of reaction system of H+-ATPase activity determination at a wavelength of 400–900 nm. (B) Enzymatic characterization of the membrane fraction. Relative H+-ATPase activity was determined by measurement of membrane-specific H+-ATPase activity compared to total ATPase activity. MF, membrane fraction; CH, crude homogenate; VM, vacuole membrane; PM, plasma membrane; MM, mitochondria membrane. Data are means ± SD of three independent experiments (biological replicates).

Analysis of Differential Membrane Proteins in Control and CHA-SQ-1 Treated Wheat Flag Leaves

To understand the membrane proteome response to short-term CHA-SQ-1-treatment of wheat flag leaves, and the changes in membrane proteomes of wheat flag leaves from the control, CHA-SQ-1-treated plants (2 and 6 h) were analyzed by 2-DE. The membrane protein maps produced from three independent protein extractions showed a high reproducibility based on analysis using PDQuest software.

Figure 3 shows a representative gel image of proteins extracted from the control and CHA-SQ-1-treated plants. Protein spots [345 (±12), 361 (±15), 370 (±13)] were reproducibly detected using PDQuest software from the control, 2 and 6 h after CHA-SQ-1 treatment, respectively (n = 3). From a spot-to-spot comparison and based on statistical analysis, a total of 150 spots (numbered from 1 to 150) exhibited at least 1.5-fold (p < 0.05) difference in abundance between the control and CHA-SQ-1 treatment (Figure 3, Supplementary Table S1). In wheat leaves at 2 h after treatment, 68 spots had >1.5-fold change in abundance (p < 0.05) and 37 spots showed a >2.0-fold change; meanwhile, in wheat leaves at 6 h after CHA-SQ-1 treatment, 83 spots had a >1.5-fold change in abundance and 41 spots showed a >2.0-fold change (Supplementary Table S1). Among 150 differential proteins, only spot 20 showed qualitative changes and was detected in CHA-SQ-1-treated flag leaves (Figure 3, Supplementary Table S1).


FIGURE 3. 2-DE image analysis of membrane proteomes in control and CHA-SQ-1-treated wheat flag leaves. (A) 2-DE map of control flag leaves; (B) 2-DE map of flag leaves treated with CHA-SQ-1 for 2 h (C) 2-DE map of flag leaves treated with CHA-SQ-1 for 6 h.

To understand the differentially expressed protein (DEP) profile patterns at two time points during CHA-SQ-1 treatment, the distribution of proteins were analyzed. Figures 4A,B present the number of DEPs under CHA-SQ-1 treatments and how these spots overlap using Venn diagram analysis. There were 35 up-regulated proteins and 70 down-regulated proteins in flag leaves at 2 h after CHA-SQ-1 treatment; meanwhile, there were 43 up-regulated proteins and 81 down-regulated proteins in flag leaves at 6 h after CHA-SQ-1 treatment (Figures 4A,B). Among the up-regulated spots, 24 spots were up-regulated at both CHA-SQ-1 time treatments; of all down-regulated spots, 52 spots were down-regulated under both two different CHA-SQ-1 treatments (Figures 4A,B). Approximately 50% (52.5% for down-regulated proteins and 44.4% for up-regulated proteins) of the DEPs exhibited a similar regulatory pattern under the two different CHA-SQ-1 treatments. Additionally, 77 spots showed up or down-regulation under only one treatment. Of these, 29 spots were found to be differentially expressed at 2 h after CHA-SQ-1 treatment, while the other 48 spots showed significant changes in response to 6 h of CHA-SQ-1 treatment, indicating that these spots were specifically responsive to short-term CHA-SQ-1-treatment of wheat flag leaves (Figures 4A,B).


FIGURE 4. Venn diagram analysis of the differentially expressed membrane protein spots and functional classification of 149 identified membrane proteins. (A) The numbers of differentially expressed protein spots with down-regulation in wheat flag leaves after CHA-SQ-1 treatment for 2 and 6 h. (B) The numbers of differentially expressed protein spots with up-regulation in wheat flag leaves after CHA-SQ-1 treatments for 2 and 6 h. (C) Distribution of proteins according to their biological functions. Unclassified proteins include those whose functions have not been described.

Identification and Functional Classification of DEPs

A total of 150 DEPs were analyzed by MALDI-TOF/TOF MS. Amongst them, 149 were successfully identified by MS/MS (Supplementary Table S2). Among the 149 identities, 142 have been functionally annotated in the current database, whereas the remaining seven identities were either unnamed proteins (spot 84, spot 92, spot 115, and spot 120) or hypothetical proteins (spot 27, spot 78, and spot 109; Supplementary Table S2). To annotate their identities, their sequences were used as a query to search for homologs using BLASTP (NCBI). The corresponding homologs with the highest similarity are listed in Supplementary Table S3. All 7 proteins shared at least 80% sequence similarity, suggesting that they may have similar function with their homologues. In summary, 149 identities represented 103 unique proteins.

Based on the metabolic and functional features of wheat flag leaves, all of the 149 identities were classified into 12 major categories, including photosynthesis, ATP synthesis and ion transport, protein folding and assembly, unclassified proteins, protein biosynthesis, cell rescue and defense, redox homeostasis, carbohydrate metabolism, protein degradation, signal transduction, protein transport, and chlorophyll biosynthesis (Figure 4C). Eighty percent of these identified proteins were implicated in the first six functional groups, whereas the largest functional groups that were greatly affected by CHA-SQ-1 treatment were proteins involved in photosynthesis (27.5%). Further analysis of the change of abundance in each group revealed that proteins involved in ATP synthesis (13.4%), protein folding and assembly (11.4%), protein biosynthesis (9.4%), and stress-related proteins (9.4%) were overrepresented, either in number or in expression level, suggesting that these processes were susceptible to CHA-SQ-1 treatment. In order to visualize the protein expression patterns of all 12 categories, hierarchical clustering of proteins was analyzed (Figure 5).


FIGURE 5. Hierarchical clustering of membrane proteins of all 12 categories in control and two different CHA-SQ-1 treatment regimes of wheat flag leaves. The hierarchical cluster analysis was conducted using MultiExperiment Viewer 4.9 and the Log-transformed values of -fold change ratios listed in Supplementary Table S4. A, photosynthesis; B, signal transduction; C, protein transport; D, chlorophyll biosynthesis; E, ATP synthesis and ion transport; F, protein biosynthesis; G, electron transport; H, protein degradation; I, protein folding and assembly; J, unclassified proteins; K, cell rescue and defense; L, carbohydrate metabolism.

In general, the apparent Mr predicted by SDS-PAGE has an error of about ±10% compared with the theoretical value. However, amongst all of the identified proteins, a set of 37 identities with known function were found with observed Mr values much smaller than the theoretical values (Supplementary Table S2), suggesting that these proteins appeared to be partially degraded products of their intact proteins. Of these, twelve identities were involved in photosynthesis (spot 26, spot 39, spot 48, spot 67, spot 71, spot 72, spot 76, spot 116, spot 121, spot 133, spot 141, and spot 145); five identities were involved in the process of protein folding and assembly (spot 74, spot 129, spot 139, spot 142, and spot 143); respectively, four identities were involved in the process of ATP synthesis (spot 33, spot 64, spot 66, and spot 111), protein biosynthesis (spot 107, spot 118, spot 132, and spot 137) and cell rescue and defense (spot 68, spot 69, spot 119, and spot 138); three identities were involved in the process of carbohydrate metabolism (spot 61, spot 77, and spot 135); two identities were related to signal transduction (spot 124 and spot 146); and respectively, one identity was suggested to be related to protein degradation (spot 123) and protein transport (spot 49) processes; one identity (spot 27) was unclassified protein.

By contrast, 17 identities with annotated function were found with observed Mr values much larger than theoretical values (Supplementary Table S2), indicating that these proteins may be products of post translation modified proteins.

Physicochemical Characteristics of Membrane Proteins

To evaluate the physicochemical characteristics of flag leaf membrane proteins, all of the identified proteins were analyzed in terms of hydrophobicity (GRAVY values) and the number of TMDs. The GRAVY score takes into account the size and the charge of the whole protein and ranges for instance from -2 to +2, positive values referring to hydrophobic proteins while negative values refer to hydrophilic proteins. The majority of analyzed proteins have a GRAVY between -0.4 and +0.4, which could not discriminate their hydrophobic or hydrophilic nature (Morel et al., 2006). In our case, the GRAVY of membrane proteins of the wheat flag leaves analyzed ranged from -1.122 to +1.125, and most of the proteins (78.5%) had a GRAVY index between -0.4 and +0.4 (Supplementary Table S2). TMDs prediction programs (TMpred) were used to predict putative TMDs in all of 149 identified proteins (Supplementary Table S2). TM candidate proteins formed 92.7% (138/149) of the proteins, displaying at least one TMD. Among them, 47.8% (66/138) exhibited one to three TMDs, 31.2% (43/138) showed four to six TMDs, and 21.0% (29/138) had 7 to 19 TMDs (Figure 6). It has been reported that most of the integral cytoplasmic membrane proteins are hydrophobic, while the majority of integral outer membrane proteins are hydrophilic (Nouwens et al., 2000; Santoni et al., 2000), which can cause the observed ambiguity between GRAVY values and TMDs. A similar phenomena was also observed in our study. Eighty-seven percent (120/138) of proteins with putative TMDs showed negative GRAVY values; while LRR disease resistance protein/transmembrane receptor kinase PS4 with no predicted TMD scored positive GRAVY values (1.125; Supplementary Table S2). According to the localization tools, the results indicated that identified proteins were higher enrichment in membrane proteins, but a number of non-membrane proteins were identified in this study (Supplementary Table S5), which is a common problem in isolating membrane from plant tissue and also observed in many other membrane proteomic studies (Li et al., 2012; Manandhar-Shrestha et al., 2013; Nie et al., 2015).


FIGURE 6. Distribution of the number of predicted transmembrane domains (TMDs) of membrane protein in wheat leaves. TMDs analysis of identified proteins was performed using TMpred program. The algorithm is based on the statistical analysis of TMbase. The number of TMDs for every DEP was listed in Supplementary Table S2.

PPI Analysis of Identified Membrane Proteins

In order to explore the relationship among all identified differential proteins, protein–protein interaction was analyzed. All 149 identified proteins were blasted against the Arabidopsis thaliana TAIR10 protein database (Supplementary Table S5). Identified proteins were grouped into functional classes according to the biological processes in which they are involved. STRING and BiNGO were used to visualize the protein–protein interaction, biological pathways, and molecular functions (Maere et al., 2005; Franceschini et al., 2013).

The STRING analysis revealed the functional links between different proteins in which proteins involved in photosynthesis, ATP synthesis, response to stress, and protein synthesis were major clusters (Figure 7). Actually, these four clusters were not separated and together they formed a related-network in response to CHA-SQ-1 treatment. Proteins overlapped among the four clusters, especially proteins involved in photosynthesis and energy metabolism. Abbreviations of the specific protein names in the network are shown in Supplementary Table S6.


FIGURE 7. Analysis of protein interaction network by STRING 9.1. TAIR homologous proteins from identified proteins were mapped by searching the STRING 9.1 software with a confidence level of 0.4. Colored lines between the proteins indicate the various types of interaction evidence.

To obtain statistically over- or under-represented categories of biological pathways and molecular functions related to CHA-SQ-1 treatment, BiNGO was used to analyze identified differential proteins (Figure 8, Supplementary Tables S7 and S8). Several overrepresented biological pathways were mostly significant (Figure 8A, Supplementary Table S7), including photosynthesis (p = 2.3849e–21), photosynthesis and light reaction (p = 2.8509e–21), generation of precursor metabolites and energy (p = 5.2254e–16), and response to stimulus (p = 6.6954e–13). Meanwhile, a complete list of the enriched Gene Ontology (GO) molecular functions for the proteins was presented in Figure 8B and Supplementary Table S8. Of them, several most highly enriched molecular functions are poly(U) RNA binding (p = 1.5997e–7), poly-pyrimidine tract binding (p = 1.5997e–7), chlorophyll binding (p = 5.8727e–7), single-stranded RNA binding (p = 8.3523e–7), hydrogen ion transporting ATP synthase activity and rotational mechanism (p = 5.2896e–6), ATPase activity (p = 5.4790e–6), NADPH dehydrogenase activity (p = 2.5621e–5), antioxidant activity (p = 4.1425e–5), ATP-dependent peptidase activity (p = 4.1601e–5), metallopeptidase activity (p = 6.8133e–5) and translation elongation factor activity (p = 7.6145e–5).


FIGURE 8. Biological pathway (A) and molecular function (B) networks generated by BiNGO. Homologous proteins were used for the GO analysis. The size of the node is related to the number of proteins and the color represents the p-value for the statistical significance of the overrepresented GO term (see the color scale on the right bottom), whereas white nodes are not enriched.


Oxidative Stress of Wheat Flag Leaves Caused by CHA-SQ-1

Various abiotic stresses lead to the overproduction of ROS in plants, which result in oxidative stress and cause damage to multiple cellular components such as proteins (Gill and Tuteja, 2010). It has been reported previously that CHA-SQ-1, in a manner similar to other stress factors, can generates excessive ROS during the anther abortion (Ba et al., 2013). And, indeed, our current results showed that CHA-SQ-1 also induced the overproduction of ROS in flag leaves and further led to degradation of multiple membrane proteins (Supplementary Table S2). Generally, in order to cope with continuous ROS production under stress, major antioxidant enzymes including SOD, POD, and CAT, are activated to scavenge ROS (Gill and Tuteja, 2010). We found an inverse correlation between the level of intracellular ROS and the activities of leaf SOD, CAT, and POD. This meant that excessive ROS cannot be eliminated effectively at the early stage of CHA-SQ-1 treatment, which aggravated damage to cellular components. In plants, overproduction of ROS can destroy membrane stability through the formation of MDA (Gill and Tuteja, 2010). Similarly, in the present study, MDA content increased significantly under CHA-SQ-1 treatment, showing that the peroxidative reaction in membrane lipids had become stronger and the plant injury would be more serious. In this case, membrane proteomes were damaged.

Proteins Involved in Photosynthesis

Photosynthesis is the primary pathway for the production of carbohydrates, which are essential for cell growth and proliferation. It is composed of two steps, photoreaction (light dependent) and dark reaction (carbon fixation). The capability of a plant to maintain a stable photosynthetic rate is significant for sustaining plant growth under stress conditions (Wang et al., 2013). Previously, our results suggest that CHA-SQ-1 affected the major products of photosynthesis during the anther abortion, such as sucrose, starch and aliphatic metabolism (Ba et al., 2014a; Zhu et al., 2015b). Similarly, in the present study, 41 DEPs were found to be associated with the photosynthetic process under CHA-SQ-1 treatment (Supplementary Table S2).

Of these, six DEPs (spot 25, spot 50, spot 67, spot 121, spot 127, and spot 136) belonged to Chlorophyll a/b binding proteins (LHCB; Supplementary Table S2). As a whole, their abundance was down-regulated in CHA-SQ-1 treated flag leaves (Supplementary Table S1). LHCB are the most abundant membrane proteins in plants and play a vital role in maintaining a stable photosynthetic rate. Previous study showed that gene expression of LHCB is down-regulated under ROS stress (Staneloni et al., 2008), and decrease LHCB expression leads to plants being more vulnerable under stress conditions (Andersson et al., 2001; Ganeteg et al., 2004; Kovács et al., 2006).

Oxygen evolving enhancer proteins 1 (OEE 1) and OEE 2 are responsible for the stability of PS II and play a vital role in catalyzing the splitting of water (Mayfield et al., 1987a,b). Study has shown that OEE 2 expression is increased in some species under drought and salt stress (Sugihara et al., 2000; Abbasi and Komatsu, 2004; Gazanchian et al., 2007). However, a few studies have shown that the expression level of OEE decreased during salt stress in potatoes and wheat (Aghaei et al., 2008; Gao et al., 2011). In our study, the abundance of OEE 1 (spot 18 and spot 19) was up-regulated at 2 h after CHA-SQ-1 treatment, and down-regulated at 6 h after CHA-SQ-1 treatment (Supplementary Tables S1 and S2). At 2 h after CHA-SQ-1 treatment, increased abundance of OEE 1 helped repair the damaged PS II. However, decreased abundance of OEE1 at 6 h after treatment caused instability of PS II under peroxide stress caused by CHA-SQ-1. Meanwhile, the abundance of OEE 2 (spot 32) was reduced at 6 h after CHA-SQ-1 treatment in flag leaves (Supplementary Tables S1 and S2), which further slowed the photosynthetic process.

Proteins Involved in ATP Production and Electron Transport

ATP synthase is the universal enzyme that manufactures ATP from ADP and provides energy for a large number of fundamental biological processes (Kang et al., 2012). In the present study, a total of 16 DEPs were involved in ATP synthesis, with 14 DEPs down-regulated and two DEPs (spot 34 and spot 64) up-regulated in response to CHA-SQ-1 treatment of flag leaves (Figure 3, Supplementary Tables S1 and S2). During biotic and abiotic stresses in plants, high energy costs are required at the stage of stress acclimation, such as the increased relative abundance of components of mitochondrial ATP-synthase (Kosová et al., 2014). Here, as a whole, the abundance of different subunits of ATP synthesis was decreased in CHA-SQ-1 treated flag leaves compared to the control plant; this interrupted multiple normal metabolic processes dependent on ATP. In addition, the abundance of vacuolar proton-ATPase subunit A (spot 46, spot 47, and spot 86) was decreased significantly in CHA-SQ-1 treated flag leaves (Figure 3, Supplementary Tables S1 and S2), which would cause an imbalance of homeostasis and abnormal cellular activity. These negative events might trigger an insufficiency of ATP production for maintaining normal metabolic processes under CHA-SQ-1 treatment in flag leaves.

The process of ATP production is always along with electron transport. Ferredoxin-NADP+ oxidoreductase (FNR) is a ubiquitous enzyme encoded by nuclear genes in higher plants. It catalyzes reversible electron transfer between ferredoxin (Fd; or flavodoxin) and NAD(P)H (Mulo, 2011). In the present study, the expression level of FNR (spot 10, spot 12, spot 14, spot 62, and spot 63) was decreased in CHA-SQ-1 treated wheat flag leaves compared to the control plants (Supplementary Tables S1 and S2). As mentioned above, LHCB was down-regulated under CHA-SQ-1 treatment, and this decrease may hinder photon capturing and the transfer the excitation energy to reaction centers to decrease NADP+ to NADPH generation (Komatsu et al., 2014). Subsequently, the marked release of FNR from the thylakoid membrane followed by a reduction in NADP+ photoreduction capacity might maintain the NADP+/NADPH homeostasis of the stressed plants (Palatnik et al., 1997; Lehtimäki et al., 2010). On the other hand, decreased reduction of NADP+ caused the photosynthetic electron transport chain to become over-reduced, therefore accelerating the formation and accumulation of ROS (Komatsu et al., 2014). As a consequence, the photosynthetic process was further impaired and slowed the growth of wheat plants and the development of pollen. In addition, NADH-ubiquinone oxidoreductase (spot 37) was down-regulated in CHA-SQ-1 treated flag leaves in our study (Figure 3, Supplementary Tables S1 and S2). NADH-ubiquinone oxidoreductase is an enzyme that detoxifies quinones and their derivatives (Wang et al., 2013) and its down-regulation in CHA-SQ-1 treated flag leaves might cause further damage to cellular metabolic process.

Proteins Involved in the Stress Response

In the present study, ROS and MDA, with the potential to cause cellular damage, were generated and accumulated. Under abiotic stresses, plants have evolved protective mechanisms to eliminate or reduce ROS and MDA levels and defend against a stressful environment (Gill and Tuteja, 2010). In our study, some anti-stress proteins, such as Temperature stress-induced lipocalin (TIL; spot 134) and 2-Cys peroxiredoxin (spot 138 and spot 140), and Abscisic acid stress ripening (ASR; spot 20 and spot 21) who play vital roles in response to stresses (König et al., 2002; Yang et al., 2005; Abo-Ogiala et al., 2014), were up-regulated to response to CHA-SQ-1 treatment in flag leaves (Figure 3, Supplementary Tables S1 and S2). To some extent, this increase can restrain accumulation of ROS. However, in our results, other antioxidative proteins, such as LRR disease resistance protein (spot 87 and spot 149) were down-regulated in CHA-SQ-1 treated flag leaves (Figure 3, Supplementary Tables S1 and S2). This means that excessive ROS was not been completely scavenged and this aggravated the damage to cellular biological processes.

Protein Metabolism-Related Proteins

Synthesis, assembling, folding, and degradation-related proteins are necessary for maintaining cellular protein homeostasis (Saikawa et al., 2004). In the present study, DEPs involved in protein metabolism are implicated in three functional subgroups: (1) protein synthesis-related proteins; (2) assembling/folding-related proteins; (3) degradation-related proteins.

In the first subgroup, five elongation factor Tu (EF-Tu; spot 7, spot 54, spot 56, spot 103, and spot 118) were up-regulated in flag leaves under CHA-SQ-1 treatment (Figure 3, Supplementary Tables S1 and S2). EF-Tu is a protein that plays a central role in the elongation phase of protein synthesis in plants. In additional, EF-Tu can prevent aggregation of denatured proteins caused by environmental stresses and chemicals (Kudlicki et al., 1997; Fu and Dooner, 2002). In our case, on the one hand, the up-regulation of EF-Tu could accumulate synthesis of proteins and replaced damaged proteins under CHA-SQ-1 treatment; on the other hand, it prevented the aggregation of damaged proteins caused by ROS.

In the second subgroup, eight DEPs (spot 43, spot 44, spot 55, spot 79, spot 81, spot 82, spot 83, and spot 122) were identified as heat shock proteins (HSPs) and chaperones and their expression levels were increased except for the 20 kD chaperonin (spot 122) whose abundance was decrease in flag leaves under CHA-SQ-1 treatment (Figure 3, Supplementary Tables S1 and S2). HSPs and chaperones play crucial roles in protecting plants against stress and they are involved in a wide range of crucial cellular processes (Wang et al., 2004; Al-Whaibi, 2011). Our results indicated that the accumulation of ROS caused instability of proteins and membrane structure under CHA-SQ-1 treatments. An increased abundance of HSPs was a protective mechanism in response to oxidative stress. However, HCF136 (spot 65) which is related to the stability of PS II, was down-regulated in CHA-SQ-1-treated wheat leaves (Figure 3, Supplementary Tables S1 and S2) and this might cause the damage of PS II.

In the third subgroup, four DEPs (spot 3, spot 4, spot 45, and spot 98) identified as cell division protease FtsH-like proteins (Supplementary Table S2). FtsH proteins are involved in photosynthesis and control of the cell cycle in eukaryotic cells. In vivo, studies have shown that FtsH proteins are responsible for the degradation of unassembled proteins (Ostersetzer and Adam, 1997) and D1 proteins (Lindahl et al., 2000; Malnoë et al., 2014). In the present study, down-regulation of FtsH protein in CHA-SQ-1 treated wheat flag leaves interrupted the process of photosynthesis and impacted the formation of thylakoid membranes (Supplementary Table S1). Meanwhile, variegation in wheat flag leaves after CHA-SQ-1 treatment may be related to a decrease in the abundance of FtsH protein (Chen et al., 2000). These events resulted in there being a low photosynthesis rate (Supplementary Figure S1) and insufficient source of energy (Supplementary Figure S2) to response to CHA-SQ-1 treatment in wheat flag leaves.

Other Proteins

In addition to the above described DEPs, 14-3-3 protein (spot 28) related to signal transduction was up-regulated in flag leaves after CHA-SQ-1 treatment (Figure 3, Supplementary Tables S1 and S2), which suggested signal transduction was enhanced in response to CHA-SQ-1 treatment; Toc75 (spot 80) is the protein translocation channel located in the outer envelope membrane of plastids (Baldwin et al., 2005) and its abundance was decreased in flag leaves after CHA-SQ-1 treatment (Figure 3, Supplementary Tables S1 and S2), which resulted in precursor proteins not being transported in the chloroplast stroma and this subsequently impaired chloroplast function.


During pollen development, wheat flag leaves are an important component of the source-sink unit as they provide sucrose and energy, and initial recipient tissue of CHA-SQ-1. The results of the present work have shown that ROS (O2- and H2O2) contents accumulated rapidly in wheat flag leaves after CHA-SQ-1 treatment and exceeded the cell threshold. An increase in ROS content was accompanied by an inhibition of SOD, CAT and POD activities. Excessive ROS were not effectively removed by the antioxidative system, which induced oxidative damage in membrane structures and membrane proteins. Comparative membrane proteomic analysis revealed that four biological processes (photosynthesis, ATP production, response to stress, and proteins metabolism) were interrupted after CHA-SQ-1 treatment in flag leaves. These abnormal biological processes caused complete dysfunction of the flag leaves. And then dysfunctional cells of flag leaf affected photosynthesis, energy supply and carbohydrate production. These results provide the basic insight needed to further investigate the mechanism of anther abortion induced by CHA-SQ-1.

Author Contributions

QS, SW, and GZ conceived the study. QS and SW performed the experiments. SW performed the MS/MS analysis. QS carried out the analysis of the data, made the identification of the proteins and drafted the manuscript. QS and SW contributed in the preparation of the final draft of the manuscript. QS, SW, YL, ZL, JG, SM, JW, and NN provided reagents, materials and analysis tools. All authors read and approved the final manuscript.

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.


The study was sponsored by the National High Technology Research and Development Program of China (No. 2011AA10A106), the National Support Program of China (No. 2015BAD27B01), the National Natural Science Foundation of China (No. 31171611, 31371697), the Technological Innovation and Over Planning Projects of Shaanxi Province (No. 2014KTZB02-01-02).

Supplementary Material

The Supplementary Material for this article can be found online at:


  1. ^
  2. ^
  3. ^
  4. ^
  5. ^
  6. ^


Abbasi, F. M., and Komatsu, S. (2004). A proteomic approach to analyze salt-responsive proteins in rice leaf sheath. Proteomics 4, 2072–2081. doi: 10.1002/pmic.200300741

PubMed Abstract | CrossRef Full Text | Google Scholar

Abo-Ogiala, A., Carsjens, C., Diekmann, H., Fayyaz, P., Herrfurth, C., Feussner, I., et al. (2014). Temperature-induced lipocalin (TIL) is translocated under salt stress and protects chloroplasts from ion toxicity. J. Plant Physiol. 171, 250–259. doi: 10.1016/j.jplph.2013.08.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Aghaei, K., Ehsanpour, A. A., and Komatsu, S. (2008). Proteome analysis of potato under salt stress. J. Proteome Res. 7, 4858–4868. doi: 10.1021/pr800460y

PubMed Abstract | CrossRef Full Text | Google Scholar

Al-Whaibi, M. H. (2011). Plant heat-shock proteins: a mini review. J. King Saud Univ. Sci. 23, 139–150. doi: 10.1016/j.jksus.2010.06.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Andersson, J., Walters, R. G., Horton, P., and Jansson, S. (2001). Antisense inhibition of the photosynthetic antenna proteins CP29 and CP26: implications for the mechanism of protective energy dissipation. Plant Cell 13, 1193–1204. doi: 10.1105/tpc.13.5.1193

PubMed Abstract | CrossRef Full Text | Google Scholar

Ba, Q. S., Zhang, G. S., Che, H. X., Liu, H. Z., Ng, T. B., Zhang, L., et al. (2014a). Aliphatic metabolism during anther development interfered by chemical hybridizing agent in wheat. Crop Sci. 54, 1458–1467. doi: 10.2135/cropsci2013.04.0277

CrossRef Full Text | Google Scholar

Ba, Q. S., Zhang, G. S., Wang, J. S., Niu, N., Ma, S. C., and Wang, J. W. (2014b). Gene expression and DNA methylation alterations in chemically induced male sterility anthers in wheat (Triticum aestivum L.). Acta Physiol. Plant. 36, 503–512. doi: 10.1007/s11738-013-1431-6

CrossRef Full Text | Google Scholar

Ba, Q. S., Zhang, G. S., Wang, J. S., Che, H. X., Liu, H. Z., Niu, N., et al. (2013). Relationship between metabolism of reactive oxygen species and chemically induced male sterility in wheat (Triticum aestivum L.). Can. J. Plant Sci. 93, 675–681. doi: 10.4141/cjps2012-280

CrossRef Full Text | Google Scholar

Baisakh, N., Ramanarao, M. V., Rajasekaran, K., Subudhi, P., Janda, J., Galbraith, D., et al. (2012). Enhanced salt stress tolerance of rice plants expressing a vacuolar H+-ATPase subunit c1 (SaVHAc1) gene from the halophyte grass Spartina alterniflora Löisel. Plant Biotechnol. J. 10, 453–464. doi: 10.1111/j.1467-7652.2012.00678.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Baldwin, A., Wardle, A., Patel, R., Dudley, P., Park, S. K., Twell, D., et al. (2005). A molecular-genetic study of the Arabidopsis Toc75 gene family. Plant Physiol. 138, 715–733. doi: 10.1104/pp.105.063289

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, M., Choi, Y., Voytas, D. F., and Rodermel, S. (2000). Mutations in the Arabidopsis VAR2 locus cause leaf variegation due to the loss of a chloroplast FtsH protease. Plant J. 22, 303–313. doi: 10.1046/j.1365-313x.2000.00738.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Curtis, T., and Halford, N. G. (2014). Food security: the challenge of increasing wheat yield and the importance of not compromising food safety. Ann. Appl. Biol. 164, 354–372. doi: 10.1111/aab.12108

PubMed Abstract | CrossRef Full Text | Google Scholar

Franceschini, A., Szklarczyk, D., Frankild, S., Kuhn, M., Simonovic, M., Roth, A., et al. (2013). STRING v9.1: protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res. 41, D808–D815. doi: 10.1093/nar/gks1094

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, H., and Dooner, H. K. (2002). Intraspecific violation of genetic colinearity and its implications in maize. Proc. Natl. Acad. Sci. U.S.A. 99, 9573–9578. doi: 10.1073/pnas.132259199

PubMed Abstract | CrossRef Full Text | Google Scholar

Ganeteg, U., Kulheim, C., Andersson, J., and Jansson, S. (2004). Is each light-harvesting complex protein important for plant fitness? Plant Physiol. 134, 502–509. doi: 10.1104/pp.103.033324

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, L., Yan, X., Li, X., Guo, G., Hu, Y., Ma, W., et al. (2011). Proteome analysis of wheat leaf under salt stress by two-dimensional difference gel electrophoresis (2D-DIGE). Phytochemistry 72, 1180–1191. doi: 10.1016/j.phytochem.2010.12.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Gazanchian, A., Hajheidari, M., Sima, N. K., and Salekdeh, G. H. (2007). Proteome response of Elymus elongatum to severe water stress and recovery. J. Exp. Bot. 58, 291–300. doi: 10.1093/jxb/erl226

PubMed Abstract | CrossRef Full Text | Google Scholar

Gigon, A., Matos, A. R., Laffray, D., Zuily-Fodil, Y., and Pham-Thi, A. T. (2004). Effect of drought stress on lipid metabolism in the leaves of Arabidopsis thaliana (ecotype Columbia). Ann. Bot. 94, 345–351. doi: 10.1093/aob/mch150

PubMed Abstract | CrossRef Full Text | Google Scholar

Gill, S. S., and Tuteja, N. (2010). Reactive oxygen species and antioxidant machinery in abiotic stress tolerance in crop plants. Plant Physiol. Biochem. 48, 909–930. doi: 10.1016/j.plaphy.2010.08.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Hopff, D., Wienkoop, S., and Luthje, S. (2013). The plasma membrane proteome of maize roots grown under low and high iron conditions. J. Proteomics 91, 605–618. doi: 10.1016/j.jprot.2013.01.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Howell, S. H. (2013). Endoplasmic reticulum stress responses in plants. Annu. Rev. Plant Biol. 64, 477–499. doi: 10.1146/annurev-arplant-050312-120053

PubMed Abstract | CrossRef Full Text | Google Scholar

Jaiswal, D. K., Ray, D., Subba, P., Mishra, P., Gayali, S., Datta, A., et al. (2012). Proteomic analysis reveals the diversity and complexity of membrane proteins in chickpea (Cicer arietinum L.). Proteome Sci. 10, 59. doi: 10.1186/1477-5956-10-59

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, G., Li, G., Xu, W., Peng, X., Han, Q., Zhu, Y., et al. (2012). Proteomics reveals the effects of salicylic acid on growth and tolerance to subsequent drought stress in wheat. J. Proteome Res. 11, 6066–6079. doi: 10.1021/pr300728y

PubMed Abstract | CrossRef Full Text | Google Scholar

Komatsu, S., Kamal, A. H., and Hossain, Z. (2014). Wheat proteomics: proteome modulation and abiotic stress acclimation. Front. Plant Sci. 5:684. doi: 10.3389/fpls.2014.00684

PubMed Abstract | CrossRef Full Text | Google Scholar

König, J., Baier, M., Horling, F., Kahmann, U., Harris, G., Schürmann, P., et al. (2002). The plant-specific function of 2-Cys peroxiredoxin-mediated detoxification of peroxides in the redox-hierarchy of photosynthetic electron flux. Proc. Natl. Acad. Sci. U.S.A. 99, 5738–5743. doi: 10.1073/pnas.072644999

PubMed Abstract | CrossRef Full Text | Google Scholar

Kosová, K., Vitámvás, P., and Prášil, I. T. (2014). Proteomics of stress responses in wheat and barley-search for potential protein markers of stress tolerance. Front. Plant Sci. 5:711. doi: 10.3389/fpls.2014.00711

PubMed Abstract | CrossRef Full Text | Google Scholar

Kota, U., and Goshe, M. B. (2011). Advances in qualitative and quantitative plant membrane proteomics. Phytochemistry 72, 1040–1060. doi: 10.1016/j.phytochem.2011.01.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Kovács, L., Damkjaer, J., Kereïche, S., Ilioaia, C., Ruban, A. V., Boekema, E. J., et al. (2006). Lack of the light-harvesting complex CP24 affects the structure and function of the grana membranes of higher plant chloroplasts. Plant Cell 18, 3106–3120. doi: 10.1105/tpc.106.045641

PubMed Abstract | CrossRef Full Text | Google Scholar

Kudlicki, W., Coffman, A., Kramer, G., and Hardesty, B. (1997). Renaturation of rhodanese by translational elongation factor (EF) Tu: protein refolding by EF-Tu flexing. J. Biol. Chem. 272, 32206–32210. doi: 10.1074/jbc.272.51.32206

PubMed Abstract | CrossRef Full Text | Google Scholar

Lehtimäki, N., Lintala, M., Allahverdiyeva, Y., Aro, E. M., and Mulo, P. (2010). Drought stress-induced upregulation of components involved in ferredoxin-dependent cyclic electron transfer. J. Plant Physiol. 167, 1018–1022. doi: 10.1016/j.jplph.2010.02.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, B., Takahashi, D., Kawamura, Y., and Uemura, M. (2012). Comparison of plasma membrane proteomic changes of arabidopsis suspension-cultured cells (T87 line) after cold and ABA treatment in association with freezing tolerance development. Plant Cell Physiol. 53, 543–554. doi: 10.1093/Pcp/Pcs010

PubMed Abstract | CrossRef Full Text | Google Scholar

Lindahl, M., Spetea, C., Hundal, T., Oppenheim, A. B., Adam, Z., and Andersson, B. (2000). The thylakoid FtsH protease plays a role in the light-induced turnover of the photosystem II D1 protein. Plant Cell 12, 419–431. doi: 10.1105/tpc.12.3.419

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J. X., Liao, D. Q., Oane, R., Estenor, L., Yang, X. E., Li, Z. C., et al. (2006). Genetic variation in the sensitivity of anther dehiscence to drought stress in rice. Field Crops Res. 97, 87–100. doi: 10.1016/j.fcr.2005.08.019

CrossRef Full Text | Google Scholar

Liu, X. Y., Wu, Y. D., Shen, Z. Y., Shen, Z., Li, H. H., Yu, X. M., et al. (2011). Shotgun proteomics analysis on maize chloroplast thylakoid membrane. Front. Biosci. (Elite Ed.) 3:250–255. doi: 10.2741/e240

PubMed Abstract | CrossRef Full Text | Google Scholar

López-Pérez, L., Martínez-Ballesta Mdel, C., Maurel, C., and Carvajal, M. (2009). Changes in plasma membrane lipids, aquaporins and proton pump of broccoli roots, as an adaptation mechanism to salinity. Phytochemistry 70, 492–500. doi: 10.1016/j.phytochem.2009.01.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Maere, S., Heymans, K., and Kuiper, M. (2005). BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics 21, 3448–3449. doi: 10.1093/bioinformatics/bti551

PubMed Abstract | CrossRef Full Text | Google Scholar

Malnoë, A., Wang, F., Girard-Bascou, J., Wollman, F. A., and De Vitry, C. (2014). Thylakoid FtsH protease contributes to photosystem II and cytochrome b6f remodeling in Chlamydomonas reinhardtii under stress conditions. Plant Cell 26, 373–390. doi: 10.1105/tpc.113.120113

PubMed Abstract | CrossRef Full Text | Google Scholar

Manandhar-Shrestha, K., Tamot, B., Pratt, E. P., Saitie, S., Brautigam, A., Weber, A. P., et al. (2013). Comparative proteomics of chloroplasts envelopes from bundle sheath and mesophyll chloroplasts reveals novel membrane proteins with a possible role in C4-related metabolite fluxes and development. Front. Plant Sci. 4:65. doi: 10.3389/fpls.2013.00065

PubMed Abstract | CrossRef Full Text | Google Scholar

Martinec, J., Feltl, T., Scanlon, C. H., Lumsden, P. J., and Macháčková, I. (2000). Subcellular localization of a high affinity binding site for D-myo-inositol 1,4,5-trisphosphate from Chenopodium rubrum. Plant Physiol. 124, 475–483. doi: 10.1104/pp.124.1.475

PubMed Abstract | CrossRef Full Text | Google Scholar

Mayfield, S. P., Bennoun, P., and Rochaix, J. D. (1987a). Expression of the nuclear encoded OEE1 protein is required for oxygen evolution and stability of photosystem II particles in Chlamydomonas reinhardtii. EMBO J. 6, 313–318.

PubMed Abstract | Google Scholar

Mayfield, S. P., Rahire, M., Frank, G., Zuber, H., and Rochaix, J. D. (1987b). Expression of the nuclear gene encoding oxygen-evolving enhancer protein 2 is required for high levels of photosynthetic oxygen evolution in Chlamydomonas reinhardtii. Proc. Natl. Acad. Sci. U.S.A. 84, 749–753. doi: 10.1073/pnas.84.3.749

PubMed Abstract | CrossRef Full Text | Google Scholar

Morel, J., Claverol, S., Mongrand, S., Furt, F., Fromentin, J., Bessoule, J. J., et al. (2006). Proteomics of plant detergent-resistant membranes. Mol. Cell. Proteomics 5, 1396–1411. doi: 10.1074/mcp.M600044-MCP200

PubMed Abstract | CrossRef Full Text | Google Scholar

Mulo, P. (2011). Chloroplast-targeted ferredoxin-NADP+ oxidoreductase (FNR): structure, function and location. Biochim. Biophys. Acta 1807, 927–934. doi: 10.1016/j.bbabio.2010.10.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Nie, L., Feng, J., Fan, P., Chen, X., Guo, J., Lv, S., et al. (2015). Comparative proteomics of root plasma membrane proteins reveals the involvement of calcium signalling in NaCl-facilitated nitrate uptake in Salicornia europaea. J. Exp. Bot. 66, 4497–4510. doi: 10.1093/jxb/erv216

PubMed Abstract | CrossRef Full Text | Google Scholar

Nouri, M. Z., and Komatsu, S. (2010). Comparative analysis of soybean plasma membrane proteins under osmotic stress using gel-based and LC MS/MS-based proteomics approaches. Proteomics 10, 1930–1945. doi: 10.1002/pmic.200900632

PubMed Abstract | CrossRef Full Text | Google Scholar

Nouwens, A. S., Cordwell, S. J., Larsen, M. R., Molloy, M. P., Gillings, M., Willcox, M. D., et al. (2000). Complementing genomics with proteomics: the membrane subproteome of Pseudomonas aeruginosa PAO1. Electrophoresis 21, 3797–3809. doi: 10.1002/1522-2683(200011)21

PubMed Abstract | CrossRef Full Text | Google Scholar

Nuhse, T. S., Bottrill, A. R., Jones, A. M., and Peck, S. C. (2007). Quantitative phosphoproteomic analysis of plasma membrane proteins reveals regulatory mechanisms of plant innate immune responses. Plant J. 51, 931–940. doi: 10.1111/j.1365-313X.2007.03192.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ostersetzer, O., and Adam, Z. (1997). Light-stimulated degradation of an unassembled Rieske FeS protein by a thylakoid-bound protease: the possible role of the FtsH protease. Plant Cell 9, 957–965. doi: 10.1105/tpc.9.6.957

PubMed Abstract | CrossRef Full Text | Google Scholar

Palatnik, J. F., Valle, E. M., and Carrillo, N. (1997). Oxidative stress causes ferredoxin-NADP+ reductase solubilization from the thylakoid membranes in methyl viologen-treated plants. Plant Physiol. 115, 1721–1727. doi: 10.1104/pp.115.4.1721

PubMed Abstract | CrossRef Full Text | Google Scholar

Rasoulnia, A., Bihamta, M. R., Peyghambari, S. A., Alizadeh, H., and Rahnama, A. (2011). Proteomic response of barley leaves to salinity. Mol. Biol. Rep. 38, 5055–5063. doi: 10.1007/s11033-010-0651-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Saikawa, N., Akiyama, Y., and Ito, K. (2004). FtsH exists as an exceptionally large complex containing HflKC in the plasma membrane of Escherichia coli. J. Struct. Biol. 146, 123–129. doi: 10.1016/j.jsb.2003.09.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Santoni, V., Molloy, M., and Rabilloud, T. (2000). Membrane proteins and proteomics: un amour impossible? Electrophoresis 21, 1054–1070. doi: 10.1002/(SICI)1522-2683(20000401)21

CrossRef Full Text | Google Scholar

Shen, H., Chen, J., Wang, Z., Yang, C., Sasaki, T., Yamamoto, Y., et al. (2006). Root plasma membrane H+-ATPase is involved in the adaptation of soybean to phosphorus starvation. J. Exp. Bot. 57, 1353–1362. doi: 10.1093/jxb/erj111

PubMed Abstract | CrossRef Full Text | Google Scholar

Sobhanian, H., Motamed, N., Jazii, F. R., Nakamura, T., and Komatsu, S. (2010). Salt stress induced differential proteome and metabolome response in the shoots of Aeluropus lagopoides (Poaceae), a halophyte C(4) plant. J. Proteome Res. 9, 2882–2897. doi: 10.1021/pr900974k

PubMed Abstract | CrossRef Full Text | Google Scholar

Staneloni, R. J., Rodriguez-Batiller, M. J., and Casal, J. J. (2008). Abscisic acid, high-light, and oxidative stress down-regulate a photosynthetic gene via a promoter motif not involved in phytochrome-mediated transcriptional regulation. Mol. Plant 1, 75–83. doi: 10.1093/mp/ssm007

PubMed Abstract | CrossRef Full Text | Google Scholar

Sugihara, K., Hanagata, N., Dubinsky, Z., Baba, S., and Karube, I. (2000). Molecular characterization of cDNA encoding oxygen evolving enhancer protein 1 increased by salt treatment in the mangrove Bruguiera gymnorrhiza. Plant Cell Physiol. 41, 1279–1285. doi: 10.1093/pcp/pcd061

PubMed Abstract | CrossRef Full Text | Google Scholar

Takahashi, D., Li, B., Nakayama, T., Kawamura, Y., and Uemura, M. (2013). Plant plasma membrane proteomics for improving cold tolerance. Front. Plant Sci. 4:90. doi: 10.3389/fpls.2013.00090

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, Y. F., Millar, A. H., and Taylor, N. L. (2012). Components of mitochondrial oxidative phosphorylation vary in abundance following exposure to cold and chemical stresses. J. Proteome Res. 11, 3860–3879. doi: 10.1021/pr3003535

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Liang, W., Xing, J., Tan, F., Chen, Y., Huang, L., et al. (2013). Dynamics of chloroplast proteome in salt-stressed mangrove Kandelia candel (L.) Druce. J. Proteome Res. 12, 5124–5136. doi: 10.1021/pr4006469

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S. P., Zhang, G. S., Song, Q. L., Zhang, Y. X., Li, Z., Guo, J. L., et al. (2015). Abnormal development of tapetum and microspores induced by chemical hybridization agent SQ-1 in wheat. PLoS ONE 10:e0119557. doi: 10.1371/journal.pone.0119557

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Welti, R., Li, W., Li, M., Sang, Y., Biesiada, H., Zhou, H. E., et al. (2002). Profiling membrane lipids in plant stress responses. Role of phospholipase D alpha in freezing-induced lipid changes in Arabidopsis. J. Biol. Chem. 277, 31994–32002. doi: 10.1074/jbc.M205375200

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, F., Zhu, Y., Müller, C., Zörb, C., and Schubert, S. (2002). Adaptation of H+-pumping and plasma membrane H+-ATPase activity in proteoid roots of white lupin under phosphate deficiency. Plant Physiol. 129, 50–63. doi: 10.1104/pp.010869

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, C. Y., Chen, Y. C., Jauh, G. Y., and Wang, C. S. (2005). A Lily ASR protein involves abscisic acid signaling and confers drought and salt resistance in Arabidopsis. Plant Physiol. 139, 836–846. doi: 10.1104/pp.105.065458

PubMed Abstract | CrossRef Full Text | Google Scholar

Ye, J. X., Wang, S., Zhang, F., Xie, D., Yao, Y., and Zhang, G. H. (2013). Proteomic analysis of leaves of different wheat genotypes subjected to PEG 6000 stress and rewatering. Plant Omics J. 6, 286–294.

Google Scholar

Zhang, H., Liang, W. Q., Yang, X. J., Luo, X., Jiang, N., Ma, H., et al. (2010). Carbon starved anther encodes a myb domain protein that regulates sugar partitioning required for rice pollen development. Plant Cell 22, 672–689. doi: 10.1105/tpc.109.073668

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, L., and Wang, A. (2012). Virus-induced ER stress and the unfolded protein response. Front. Plant Sci. 3:293. doi: 10.3389/fpls.2012.00293

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, Q., Song, Y., Zhang, G., Ju, L., Zhang, J., Yu, Y., et al. (2015a). De Novo assembly and transcriptome analysis of wheat with male sterility induced by the chemical hybridizing Aagent SQ-1. PLoS ONE 10:e0123556. doi: 10.1371/journal.pone.0123556

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, W. W., Ma, S. C., Zhang, G. S., Liu, H. Z., Ba, Q. S., Li, Z. K., et al. (2015b). Carbohydrate metabolism and gene regulation during anther development disturbed by chemical hybridizing agent in wheat. Crop Sci. 55, 868–876. doi: 10.2135/cropsci2014.04.0333

CrossRef Full Text | Google Scholar

Keywords: wheat, flag leaves, membrane proteome, CHA-SQ-1, ROS

Citation: Song Q, Wang S, Zhang G, Li Y, Li Z, Guo J, Niu N, Wang J and Ma S (2015) Comparative proteomic analysis of a membrane-enriched fraction from flag leaves reveals responses to chemical hybridization agent SQ-1 in wheat. Front. Plant Sci. 6:669. doi: 10.3389/fpls.2015.00669

Received: 20 May 2015; Accepted: 13 August 2015;
Published: 26 August 2015.

Edited by:

Hans-Peter Mock, Leibniz Institute of Plant Genetics and Crop Plant Research, Germany

Reviewed by:

Ben M. Long, Australian National University, Australia
Katja Witzel, Leibniz Institute of Vegetable and Ornamental Crops, Germany

Copyright © 2015 Song, Wang, Zhang, Li, Li, Guo, Niu, Wang and Ma. 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: Gaisheng Zhang, College of Agronomy, Northwest Agriculture and Forestry University, National Yangling Agricultural Biotechnology and Breeding Center, Yangling Branch of State Wheat Improvement Centre, Wheat Breeding Engineering Research Center, Ministry of Education, Key Laboratory of Crop Heterosis of Shaanxi Province, Yangling 712100, China,

These authors have contributed equally to this work.