- Research article
- Open Access
Physiological and proteome studies of maize (Zea mays L.) in response to leaf removal under high plant density
BMC Plant Biologyvolume 18, Article number: 378 (2018)
Under high plant density, intensifying competition among individual plants led to overconsumption of energy and nutrients and resulted in an almost dark condition in the lower strata of the canopy, which suppressed the photosynthetic potential of the shaded leaves. Leaf removal could help to ameliorate this problem and increase crop yields. To reveal the mechanism of leaf removal in maize, tandem mass tags label-based quantitative analysis coupled with liquid chromatography–tandem mass spectrometry were used to capture the differential protein expression profiles of maize subjected to the removal of the two uppermost leaves (S2), the four uppermost leaves (S4), and with no leaf removal as control (S0).
Excising leaves strengthened the light transmission rate of the canopy and increased the content of malondialdehyde, whereas decreased the activities of superoxide dismutase and peroxidase. Two leaves removal increased the photosynthetic capacity of ear leaves and the grain yield significantly, whereas S4 decreased the yield markedly. Besides, 239 up-accumulated proteins and 99 down-accumulated proteins were identified between S2 and S0, which were strongly enriched into 30 and 23 functional groups; 71 increased proteins and 42 decreased proteins were identified between S4 and S0, which were strongly enriched into 22 and 23 functional groups, for increased and decreased proteins, respectively.
Different defoliation levels had contrastive effects on maize. The canopy light transmission rate was strengthened and proteins related to photosynthetic electron-transfer reaction were up-regulated significantly for treatment S2, which improved the leaf photosynthetic capacity, and obtained a higher grain yield consequently. In contrast, S4 decreased the grain yield and increased the expressions of proteins and genes associated with fatty acid metabolism. Besides, both S2 and S4 exaggerated the defensive response of maize in physiological and proteomic level. Although further studies are required, the results in our study provide new insights to the further improvement in maize grain yield by leaf removal.
Maize (Zea mays L.) yield has advanced through breeding complemented with evolving management technologies . Increasing the maize plant population is an effective practice that has undergone a constant evolution over the years . Under high plant densities, however, intensifying competition occurred among individual plants, and led to overconsumption of energy and nutrients including stronger root systems or bigger leaf area . Meanwhile, the close distances between plants in the group resulted in an almost dark condition in the lower strata of the canopy, which suppressed the photosynthetic potential of the leaves [4, 5]. Nevertheless, leaves in the middle canopy are the main source of corn grain yield, and the photosynthetic intensity is closely related to yield production [6, 7]. Besides, shading condition could accelerate the reduction of chlorophyll (Chl) content and leaf area of leaves at lower canopy status .
Excising vegetative organs partially is an effective method to modify the canopy structure, which is benefit to improve the light environment within the canopy, and ultimately alter crop yield [9,10,11]. Nevertheless, the response of yield to leaf removal levels differs greatly . When plants are injured after artificial defoliation, eaten by animals or pests, the leaf area decrease thereafter [13, 14], however residual organs have a compensating effect when the photosynthetic organs injured above a certain threshold level [15, 16]. The effect (negative, positive, or zero) of source-reducing on plants growth depends on the frequency and intensity of defoliation . Liu et al.  has demonstrated that defoliation above the cob decreased leaf area index significantly, whereas it markedly improved light condition within the canopy. Besides, removing the uppermost two or four leaves in maize appeared to stimulate an increase in net photosynthetic rate (Pn), stomatal conductance, and Chl content of the ear leaf. Hao et al.  also evidenced that an increase in Pn of the remaining ear leaf came up by excising 1/4 and 1/2 of maize leaves over the whole plant. Increased intensity of leaf removal, however, do not conduce to maintain the photosynthetic ability of remaining leaves during late filling stage . The photosynthesis extent of leaves during grain filling can be affected by canopy structure  and the corresponding variations in light conditions may lead to changes in the expression levels of proteins, which invariably leads to changes in plant metabolism .
Leaf removal has also been reported affecting antioxidant metabolism of plants , for instance, altered the activities of superoxide dismutase and peroxidase as well as the content of malondialdehyde [19, 23]. To date, though several physiology variations induced by leaf removal have been studied in maize and other plants, there is still little published information at the proteomic level regarding the effects of leaf removal on maize characteristics under high plant density. Therefore, we employed a quantitative proteomic analysis based on tandem mass tag (TMT) labels, coupled with liquid chromatography-tandem mass spectrometry (LC-MS/MS), to capture the differential protein expression profiles of maize subjected to defoliation. This research compared changes in physiology and proteins induced by different leaf removal treatments using a high-yield and density-tolerance variety under a optimized density, hoping to elucidate the physiological mechanism of leaf removal on maize production and provide a theoretical basis for further improvement in maize grain yield.
The experiment was conducted at the Corn Research Center of Shandong Agricultural University, Tai’an, Shandong Province, China (36° 10′ N, 117° 09′ E) from June 18 to October 8 during 2015 growing season. This area has a semi-humid, warm temperate continental climate with monsoons. The average content of organic matter in the tillage layer was 18.6 g kg− 1 and the total nitrogen (N), rapidly available phosphorous (P), and rapidly available potassium (K) were 1.03 g kg− 1, 43.05 mg kg− 1 and 78.91 mg kg− 1, respectively.
The summer maize hybrid Denghai 618 (a high-yield and density-tolerance variety grown extensively in North China) was selected as the material for testing. Maize seeds were planted with hand planters at a uniform density of 9.75 plants m− 2, which was a optimum density for Denghai 618 selected during 2013 to 2014 growing seasons (a relatively high density for the growing conditions of the North China Plain). Pre-sowing, phosphorus (P2O5) and potassium (K2O) fertilizer were applied at a rate of 90 kg·ha− 1 and 120 kg ha− 1 per plot, respectively. Urea (N 46%) was applied by furrow at six and twelve leaves unfolded stage respectively, at a rate a 180 kg ha− 1 each time.
Three treatments were set up in our study, including the uppermost two leaves removal (S2), the uppermost four leaves removal (S4) and the control with no leaf removal (S0). Plants were grown until silking stage, when leaf removal treatments were imposed. Each treatment had three replicate plots, with each plot area measuring 3 m × 15 m, and the spacing between rows was 0.6 m. Besides, irrigation, weeds, diseases and insect pests were controlled adequately during the whole growing season so that no factors other than leaf removal affect plants’ growth.
For plant sampling, the uniform and healthy maize plants were marked at silking stage. The middle portion of five marked ear leaves from five individual plant of each plot was collected and mixed as one replicate at three days and seven days after leaf removal and plunged directly into liquid nitrogen, then stored at − 80 °C prior to analysis. The remaining marked plants were used to determine photosynthetic parameters. At physiological maturity, 20 ears from three center rows of each plot were harvest to measure yield (adjusted to a moisture content of 15.5%), kernel number per ear (KN) and 1000-kernel weight (TKW). Harvest index (HI) was calculated as the ratio of grain yield to the total above-ground biomass.
The plant canopy digital image analyzer (CI-100, CID Bio-Science, Inc. USA) was used to calculate the light transmission rate, and the hemispheric gray images of ear and bottom layers were also taken. In each plot, the photosynthetic effective radiation (PAR) at the top, the ear and bottom (four leaves below the ear leaf) layers were taken. The PAR values at ear and bottom layers for each plot were the average of five measurements. Light transmission rate (%) was calculated as the following equation.
Light transmission rate (%) = PAR of the determined layer (ear or bottom) / PAR of the top canopy layer × 100%.
Gas exchange parameters, which including net photosynthetic rate (Pn), stomatal conductance (gs) and intercellular CO2 concentration (Ci), were measured using a portable photosynthesis system (CIRAS-II, UK). The artificial light was set at 1600 μmol m− 2 s− 1 and CO2 concentration in the leaf chamber was maintained at 360 μmol mol− 1 using a CO2 injector. The measurements were conducted between 09: 00 AM and 11: 00 AM and each treatment had three replications.
Three representative plants were selected to determine the green leaf area (GLA) nondestructively and leaf area index (LAI) was then calculated. The equations were as GLA = ∑ (leaf length × maximum width × 0.75); LAI = GLA × n / S, where n is the number of plants within a unit area of land and S is the unit area of land.
Leaf chlorophyll content was determined using spectrometry, following standard methods . Nitroblue tetrazolium and guaiacol colorimetry methods  were used to measure the activities of superoxide dismutase (SOD) and peroxidase (POD), respectively. The content of malondialdehyde (MDA) was measured with thio-barbituric acid method .
TMT-based quantitative proteomics analysis
Samples collected at three days after leaf removal were used for the TMT-based proteomics analysis (3 biological replicates × 3 treatments). Total proteins from each sample were extracted using the trichloroacetic acid (TCA)-acetone precipitation method. First, the samples were ground in liquid nitrogen and transferred to 5-mL centrifuge tubes. Then, lysis buffer (8 M urea, 1% Triton-100, 65 mM DTT and 0.1% Protease Inhibitor Cocktail) was added to the tubes, which were sonicated three times on ice using a high-intensity ultrasonic processor (Scientz). Next, the remaining debris was removed by centrifugation at 20,000 g at 4 °C for 10 min. Finally, the protein was precipitated with cold 15% TCA for 2 h at − 20 °C. The supernatant was discarded after centrifuging at 4 °C for 10 min. The remaining precipitate was washed with cold acetone three times. The protein was resuspended in buffer (8 M urea, 100 mM TEAB, pH 8.0), and the protein concentration was determined with a 2-D Quant kit according to the manufacturer’s instructions. Then, trypsin digestion was carried out with 10 mM DTT for 1 h at 37 °C, and 20 mM IAA was added to alkylate the proteins for 45 min at room temperature in the dark. This protein sample was diluted by adding 100 mM TEAB to a urea concentration of less than 2 M. Finally, trypsin was added in a 1:50 w/w ratio of trypsin to protein for a first overnight digestion and 1: 100 w/w ratio of trypsin to protein for a second 4-h digestion. Approximately 100 μg protein for each sample was digested with trypsin for the following experiments.
After trypsin digestion, six-plex TMT labelling (Thermo Scientific) was performed following the manufacturer’s protocol. Briefly, one unit of TMT reagent (defined as the amount of reagent required to label 50 μg of protein) was thawed and reconstituted in 24 μL ACN. The peptide was reconstituted in 0.2 M TEAB, mixed with the TMT reagent, and incubated for 2 h at room temperature. The samples were desalted in a Strata X C18 SPE column (Phenomenex) and vacuum-dried. Each dried peptide sample was fractionated using high-pH reverse-phase HPLC with an Agilent 300 Extend C18 column (5-μm particles, 4.6-mm ID, 250-mm length). Eighteen fractions were collected.
The peptides were dissolved in 0.1% formic acid (FA) and loaded directly onto a reversed-phase pre-column (Acclaim PepMap 100, Thermo Scientific). The peptides were separated using a reversed-phase analytical column (Acclaim PepMap RSLC, Thermo Scientific). The peptide samples were subsequently eluted with a four-step linear gradient of solvent B (0.1% FA in 98% ACN): 6–22%, 26 min; 22–35%, 8 min; 80%, 3 min; and 80%, hold, 5 min. A constant flow rate was maintained as 300 mL/min with an EASY-nLC 1000 ultra-performance liquid chromatography (UPLC) system. The resulting peptides were processed using a Q Exactive™ Plus hybrid quadrupole-Orbitrap mass spectrometer (Thermo Fisher Scientific) coupled online to the UPLC. The MS was processed with a data-dependent procedure that alternated between single-MS and MS/MS scans. Intact peptides were detected in the Orbitrap (350–1800 m/z, 70,000 resolution) and subjected to 20 MS/MS scans using an NCE setting of 31. The top 20 precursor ions above a threshold ion count of 1E4 in the MS survey scan were identified with 30.0-s dynamic exclusion. Ion fragments were detected in the Orbitrap at 17,500 resolution. Automatic gain control (AGC) was set as 5E4 ions to prevent overfilling of the ion trap.
The Mascot search engine (v.2.3.0) was used to search the resulting MS/MS data against the UniProt Zea mays database (58,493 sequences). The cleavage enzyme was specified as trypsin/P, and two missing cleavages were allowable. The mass error was set to 10 ppm for precursor ions and to 0.02 Da for fragment ions. Carbamidomethyl on Cys, TMT-6plex (N-term), and TMT-6plex (K) were specified as fixed modifications, and oxidation of Met was specified as a variable modification. FDR was adjusted to ≤1%, and the peptide ion score was set at ≥20.
The UniProt-GOA database (http://www.ebi.ac.uk/GOA) was used to obtain the Gene Ontology (GO) proteome annotation. First, the identified protein ID was converted to a UniProt ID, and then these were mapped to GO IDs using the protein ID. If some identified proteins were not annotated by the UniProt-GOA database, InterProScan was used to annotate the GO function of the protein based on a protein sequence-alignment method. Then, the proteins were classified into three categories using the GO annotation: biological process, cellular component, and molecular function. For each category, a two-tailed Fisher’s exact test was used to test the enrichment of the differentially expressed protein against all identified proteins. Correction for multiple hypothesis testing was carried out using standard false discovery rate control methods. A GO with a corrected P-value ≤0.05 was considered significant.
Quantitative real-time polymerase chain reaction (qRT-PCR)
The analysis of qRT- PCR was performed following the method of Wang et al. . Total RNA was extracted from 0.05–0.1 g maize ear leaf by the use of RNAiso Plus reagent (Takara Bio, Japan). The cDNA were synthesized by using HiScript II Q RT SuperMix for qPCR (+gDNA wiper) (Vazyme Bio, China). Specific primers for each gene tested in our study are listed in Additional file 1: Table S1. The relative expression level of each gene was calculated according to the 2−ΔΔCt method, using Zmactin as the reference gene. The equation is:
SPSS 18.0 (SPSS Institute Inc.) was used to perform analyses of variation (ANOVAs) for physiological parameter. The results for each parameter are presented as the means of the three replicates (except for kernel numbers). Differences were judged by the least significant differences (LSD) test, and the significance level was set at the 0.05 probability level. Figures were plotted using SigmaPlot 12.0 (Systat Software Inc.).
Yield and physiological indices
The grain yield and yield components were different between leaf removal treatments and the control (Table 1). Relative to S0, S2 plants obtained significantly greater (P ≤ 0.05) 1000-kernel weight, total dry matter and harvest index, which resulted in an increase in final yield of 5.2%. In contrast, S4 plants had significantly lower (P ≤ 0.05) kernel numbers, 1000-kernel weight, total dry matter and harvest index, which resulted in a decrease in final yield of 11%.
Leaf area index (LAI) were decreased significantly after leaf removal (P ≤ 0.05, Table 2). Compared to S0, the LAI after defoliation in S2 and S4 decreased for 5.8 and 19.5%, respectively. Light transmission rate of canopy was significantly enhanced (P ≤ 0.05) at the level of the ear leaf strata and the bottom leaf strata with the increased levels of leaf removal (Fig. 1). In addition, by watching the hemispheric gray images (Fig. 2), we also intuitively found that there was an increased light transmittance, especially in the ear layer, induced by leaf removal.
Leaf removal treatments also affected gas exchange parameters of ear leaves (Fig. 3). Net photosynthesis rate (Pn) and stomatal conductance (gs) were significantly increased in S2 compared to S0, whereas intercellular CO2 concentration (Ci) was significantly decreased after leaf removal. In S4 plants, Pn was not significantly changed compared to the control at three days after leaf excising, but decreased significantly at seven days after leaf excising. Besides, chlorophyll content had the same trend as Pn in response to different leaf removal levels (Fig. 3d).
Table 2 showed the dynamic activities of SOD, POD and content of MDA in maize ear leaves after defoliation. The activities of SOD and POD decreased obviously after leaf removal in both stages compared to the control, and these indices decreased faster (P ≤ 0.05) for S4 treatment. In contrast, MDA content was increased in both treatments compared to the control.
Identification of differentially accumulated proteins
The mass error of all the identified peptides met the requirements (centered at 0 and set within 10 ppm). Besides, most peptides were distributed in 8–20 amino acid residues (sample preparation reached the standard). After merging data from three biological replicates, a total of 3586 proteins were identified, and the repeatability of the three replicates were tested using Person’s correlation coefficient (Additional file 2: Figure S1). We considered a ratio of > 1.3 to indicate up-regulation and a ratio of < 0.77 (1/1.3) to indicate down-regulation (P ≤ 0.05). Using these two criteria, we identified differentially abundant proteins in leaves subjected to leaf removal. We identified 239 increased proteins and 99 decreased proteins between the S2 and S0 treatments, and 71 increased proteins and 42 decreased proteins between the S4 and S0 treatments (Additional file 3: Table S2 and Additional file 4: Table S3).
Bioinformatic analysis of differentially abundant proteins between the S2 and S0 treatments
To identify the significantly enriched GO functional groups of differentially expressed proteins, GO annotation was conducted. The up-accumulated proteins with S2 treatment were strongly enriched into 30 functional groups compared with S0 (Additional file 5: Figure S2A), of which biological processes, cellular components and molecular functions accounted for 14, 8, and 8 GO terms, respectively. We found that 21 proteins that were up-accumulated in S2 compared with S0 played roles in photosynthesis (Fig. 4a), including “protein-chromophore linkage”, “photosynthesis, light harvesting”, “photosynthesis, light reaction”, “photosynthesis”, “response to light stimulus”, “response to radiation”, “response to red or far red light”, and “chlorophyll biosynthetic process”. Proteins that were down-accumulated in S2 compared with S0 were strongly enriched into 23 functional groups (Additional file 5: Figure S2B), of which biological processes, cellular components, and molecular functions accounted for 14, 1, and 8 GO terms, respectively. We found that 13 proteins of the down-accumulated proteins were involved in disease/defence categories (Fig. 4a), including “phenylpropanoid biosynthetic process”, “phenylpropanoid metabolic process”, and “response to stress”.
Bioinformatic analysis of differentially abundant proteins between S4 and S0 treatments
Proteins that were up-accumulated in S4 compared with S0 were strongly enriched into 22 functional groups (Additional file 5: Figure S2C), of which biological processes and molecular functions accounted for 14 and 8 GO terms, respectively. We found that 8 up-accumulated proteins were involved in fatty acid metabolism (Fig. 4b), including “fatty acid biosynthetic process”, “fatty acid metabolic process”, “monocarboxylic acid biosynthetic process”, “lipid biosynthetic process”, “cellular lipid metabolic process,” and “regulation of macromolecule metabolic process”. The proteins that were down-accumulated in S4 compared with S0 were strongly enriched into 23 functional groups (Additional file 5: Figure S2D), of which biological processes, cellular components and molecular functions accounted for 14, 1 and 8 GO terms. Besides, 10 down-accumulated proteins, which similar to the two leaves removal treatment, were enriched in disease/defence (Fig. 4b), including “response to stress”, “response to stimulus”, “defence response”, and “response to biotic stimulus”.
qRT-PCR analysis of the expression of genes between treatments and control
We next assayed whether leaf removal treatment had effects on the relative expression of genes encoding proteins changed in S2 and S4. Quantitative RT-PCR was employed to determine seventeen genes' (gpm571, LOC100273752, LOC100284847, LOC100281879, LOC100282512, lhcb6, pco103778a, psbB, psbC, Lhcb5–1, ACC1, LOC100281026, LOC100383323, LOC103634525, gpm853, Zlp, and LOC100280979) relative expression in S2 and S4. In accordance with the protein results, the relative expression of genes involved in the photosynthesis pathways were mostly increased in S2 treatment (Additional file 6: Figure S3B), and the expression of genes involved in the fatty acid metabolism were increased in S4 treatment (Additional file 6: Figure S3D). Moreover, the relative expressions of defense-related genes were decreased in both S2 and S4 (Additional file 7: Figure S4B).
Leaf removal affected the grain yield and physiological parameters of maize
The maize grain yield mainly depends on photosynthesis production by leaves after silking, and the subsequent biomass allocation to kernels . An optimized canopy structure can enhance the light utilization of plants, inhibit protein degradation in leaf and maximize grain yield [29, 30]. In our study, S2 enhanced the light transmission rate of both ear and bottom layers (Fig. 1), which enabled the leaves in lower canopy to obtain more light energy and achieve a higher grain yield ultimately . Although LAI decreased (Table 2), Chl content and net photosynthetic rate of ear leaf were enhanced after two leaves removal , which may account to the positive compensatory effect of plants . On the contrary, a higher amount of leaf removal (S4) resulted in a significant decrease in grain yield compared to the control, which might be due to the insufficient sources to favor the formation of assimilates after four leaves removal. Reactive oxygen species (ROS) are toxic molecules which can cause early senescence and ultimately cell death, and antioxidant enzymes play important roles in detoxifying ROS . In our research, the activities of antioxidant enzymes reduced, whereas the MDA content increased significantly in S2 and S4 compared to control (Table 2), which demonstrated the extend of peroxidation of membrane lipid was aggravated due to leaf removal after silking. These different physiological reactions between two and four leaves removal in this research indicated that the degree of defoliation affected maize production involved diverse processes. In order to obtain deeper insight into the nature of leaf removal, we focused on a number of proteins involved in notable function categories.
Two leaves removal enhanced the expression of photosynthesis-related proteins
Chlorophyll molecules are important photoreceptor pigments that absorb light energy and transfer electrons into the photosynthesis reaction centre . Magnesium chelatase catalyses the magnesium-insertion process in the synthesis of Chl. The mutant gene GUN5 encodes the Mg-chelatase H (Chl H) subunit of Mg-chelatase, which determines the pale phenotypes of this mutant . Magnesium-protoporphyrin IX monomethyl ester (oxidative) cyclase (Chl 27), which is involved in chlorophyll biosynthesis, catalyses the formation of protochlorophyllide. In our study, one Chl H (K7U7W9) and one Chl 27 (K7USR3) were observed to be up-accumulated with S2 treatment compared to S0 treatment, which may account for the increase of Chl content in ear leaves with S2 treatment (Fig. 3d).
Photosynthesis comprises two sets of reactions: photosynthetic electron-transfer reaction and carbon-fixation reaction. In the current research, we found that removing two leaves affected a series of proteins involved in this process (Additional file 3: Table S2). Photosynthetic electron-transfer reaction involves three key events. Firstly, the antenna complexes capture photons and produce high-energy electron. Next, photosystem II (PSII) catalyzes light-driven oxidation of water, releasing oxygen and electrons in this process. Then, the electrons take part in ATP synthesis via the electron transport chain. Finally, the electrons are transferred to photosystem I (PSI) to produce NADPH . We identified a series of proteins involved with this process, including chlorophyll a-b binding proteins (CABs; A0A096RF43, A0A096RM67, A0A096UJK9, A0A096S5Z5, B4FV94, B4FXB0, B6SZT9, B6T892, K7TXI5, and Q41746), PSII reaction centre protein (P05641, P24993, and P48187), oxygen-evolving enhancer protein (A0A096U686), cytochrome oxidase (A0A096U038 and K7UZJ0), cytochrome (A0A096Q1T0 and B6UBZ9), plastocyanin (B6SSB9), PSI assembly protein (A0A096TR75), PSI reaction centre subunit (B4G1K9), ferredoxin (B6TVC7), ATP synthase (K7VI25, K7VN08, P00835, P17344, and P48186) and F1F0-ATPase inhibitor protein (B6T5U0). In plants, CABs can capture light and transfer the excitation energy to PSI and PSII, which plays a central role in the light-harvesting complex (LHC). However, numerous environmental stressors can affect the expression of CABs [36, 37]. In the current study, ten CABs were identified and were all up-accumulated with S2 treatment. As the accumulation of CABs, PSI, and PSII can be regulated by the Chl content , the up-regulation of CABs in the ear leaves may relate to the increase in Chl content with S2 treatment in this study.
The plant PSII core complex has about 20 subunits, consisting of individual proteins and protein complexes [39, 40]. The present study identified four PSII reaction centre proteins as up-accumulated with S2 treatment, including one PSII CP47 reaction centre protein (PsbB), one PSII CP43 reaction centre protein (PsbC), one PSII reaction centre protein H (PsbH), and one oxygen-evolving enhancer protein 3–1 (OEE3). PsbB and PsbC, which comprise the PSII reaction centre, play important roles in water splitting [39, 41]. The down-regulation of these two proteins can completely destroy the oxygen-forming capacity of plants . PsbB protein also binds several small transmembrane subunits, except pigments [40, 42]. PsbH is a single transmembrane helix subunit that binds within the PsbB protein as a small subunit, and it plays a key role in the proper functioning of PSII and its stable assembly. In this work, the expression of PsbH was increased with S2 treatment compared with the control. Based on previous research , we postulated that increased accumulation of PsbH can stabilise the PSII core and promote the PSII electron transfer between the quinone acceptors QA and QB, to some extent. OEE proteins help to increase the efficiency of the oxygen-evolving complex  and several proteomics studies have shown that the abundance of OEE protein in plants is affected by stresses, nevertheless the underlying mechanism remains unknown [44, 45]. In our study, this protein was up-accumulated with S2 treatment. The function of OEE in the oxygen-evolving complex is considered supplementary, therefore, the up-regulation of this protein may represent a mechanism for the optimisation of oxygen-evolving complex. These results could demonstrate that removing two leaves could affect a series of proteins involved in the system of PSII significantly.
Cytochrome c oxidase and cytochrome bc1 complex are located within the inner mitochondrial membrane, and they function as part of the electron transfer complexes. Meanwhile, cytochrome b is located within the cytochrome b6f and bc1 complexes as part of the electron transport chain . Otherwise, plastocyanin is a copper-containing protein that can receive the electrons from the reduced cytochrome b6f and transmit them to PSI complexes in the photosynthetic electron transfer chain. There are few reports on the change in the abundances of these proteins involved in photosynthesis in response to leaf removal. In our study, the abundances of cytochrome c oxidase and cytochrome bc1 complex were increased but the abundances of cytochrome b and plastocyanin were decreased with S2 treatment. These conflicting results may attribute to that cytochrome b is one part of cytochrome bc1, and further studies are required to explore the underlying reasons.
PSI comprises pigment protein super-complexes in higher plants that have approximately 19 subunits [47, 48]. Similar to PSII mentioned above, PSI complexes are also related to a light harvest antenna (LHCI). In our study, the PSI reaction centre proteins were up-accumulated with S2 treatment. Moreover, ferredoxin (Fdx), a reducing agent that catalyses the formation of NADPH using NADP+  was more abundant in S2 treatment than in the control. Previous studies also showed that PSI can use the light energy collected in LHCI to generate reduced ferredoxin using plastocyanin oxidation [47, 48, 50]. This result may account for the decreased abundance of plastocyanin mentioned above.
ATP synthase, which involved in photosynthetic electron-transfer reaction, plays a key role in the non-photochemical quenching of photosynthesis through reactive oxygen species (ROS)-promoted photo inhibition . All differentially expressed ATP synthase was up-accumulated with S2 treatment. Moreover, we also found that one ATPase inhibitor protein was down-accumulated with S2 treatment. The high levels of ATPase and low levels of ATPase inhibitor protein that we observed in maize ear leaves with the two-leaf removal treatment may result in increased ATP synthesis, thereby decreasing ROS-promoted photo inhibition.
All the proteins discussed above were involved in photosynthetic electron-transfer reaction, and were more abundance in S2 treatment compared to S0. In concordance with these results, leaf net photosynthetic rate increased apparently in S2 treatment in our research (Fig. 3a). Moreover, qRT-PCR results showed that a series of genes which participated in photosynthesis pathway were also more expressed in S2 compared to S0 (Additional file 6: Figure S3B). The combined results in this study supported the conclusion that two leaves removal enhanced the expression of photosynthesis-related proteins and hence increased the capacity of leaf photosynthesis.
Four leaves removal increased activities of key enzymes associated with fatty acid metabolism
Unlike with the two-leaf removal treatment, the up-accumulated proteins from the removal of four leaves were strongly enriched into metabolism and secondary metabolism categories. Fatty acids are important in plants because, when metabolised, they produce many ATP molecules. In the fatty acid elongation process, the initial and rate-limiting step is catalysed by the membrane-bound 3-ketoacyl-CoA synthase, which was first identified in Arabidopsis thaliana . Acetyl-CoA carboxylases (ACCs) catalyse the formation of malonyl-CoA using acetyl-CoA, and malonyl-CoA is an important substrate in de novo lipogenesis . Acyl carrier protein (ACP) is an independent protein in dissociative type II fatty acid synthase (FAS) found in plants and other organisms [54, 55]. ACP plays a central role in the FAS system by shuttling acyl chain intermediates in its hydrophobic cavity to various enzymes. We identified one acetyl-CoA carboxylase 2, one acyl carrier protein, and one 3-ketoacyl-CoA synthase with S4 treatment; their abundances were all increased relative to the control. Besides, qRT-PCR results also showed that two genes (ACC1, LOC100281026) related to fatty acid metabolism were more expressed in S4 treatment compared to S0 (Additional file 6: Fig. S3D). These results indicated that removing four leaves can promote the activity of fatty acid metabolism in maize ear leaves.
Both leaf removal treatments decreased the expression of defense-related proteins
Plants lack of animals’ immune system. Instead, plants have evolved a set of defence mechanisms to protect themselves when attacked by pathogens under natural conditions . For example, phenylpropanoids play central roles in many aspects of the plant responses to abiotic and biotic stimuli and are becoming important indicators in plant’s stress responses to light changing . As previous research shows, decreasing the phenylpropanoid biosynthesis rate can significantly lower plant resistance [58, 59]. In our study, we identified three proteins (B4FQP4, K7VC35, and Q6VWJ0) related to the phenylpropanoid biosynthetic process, which abundances were decreased with S2 treatment. Besides, we also found that a series of proteins (A0A096RTN1, A0A096T686, A0A0B4J3G7, B4FA32, B4FR89, B6SIF0, B6SQM0, C0HGH7, K7VC35, K7VH58, and P33679) involved with the defence categories were down-accumulated with S2 treatment. Similar results have been reported in rice that received a stress .
Moreover, a series of proteins that play roles in plant defences were identified in S4 treatment. Many host proteins are induced in plants during pathogen attacks, and the majority are pathogen-related (PR) proteins. PR proteins are categorised into 17 families (PR1 to PR 17) by their structures and biological activities . Among these families, PR10 proteins play vital roles in resisting biotic and abiotic stresses . In this research, we identified one PR protein 10b down-accumulated with S4 treatment compared to the control, which suggested that leaf removal may also affect the capability of maize to resist stresses.
The level of ROS in plants always increases rapidly in response to abiotic or biotic stresses [63, 64]. Peroxidases catalyse the reduction of peroxide or hydroperoxides using an oxidised donor substrate (typically a thiol), thereby regulating H2O2 levels. In our study, we identified two peroxidases (A0A0B4J3G7 and B4FA32) that were differentially expressed with S4 treatment. The down-regulation of these two peroxidases indicates that changes in ROS levels occur in ear leaves with S4 treatment. As investigations in plants under different abiotic stresses [63,64,65,66], we supposed that peroxidase was damaged after four leaves removing and therefore cannot detoxify ROS-induced lipid peroxidation products.
Asr proteins function in response to abiotic stresses in plants [67,68,69,70,71]. Therefore, the down-regulation of Asr protein (A8IK79) in ear leaves with S4 treatment indicated that removing four leaves weakened the tolerance of maize to abiotic stress. Zeamatin was first identified in corn seeds with high amino acid homology to PR-5 proteins. Besides, it has potent antifungal activity against a number of plant pathogens . We identified one zeamatin, the expression of which was decreased significantly in ear leaves with S4 treatment. Otherwise, the activities of antioxidant enzymes were found significantly decreased in both treatments (Table 2). In addition, some genes encoded the defence-related proteins were also less expressed in both S2 and S4 treatments compared to S0 (Additional file 7: Figure S4B). Therefore, to a certain extent, both two and four leaves removal may affect the defensive system of maize in both physiological and molecular level.
Plant yield production and defense response are trade off regarding the energy distribution. Due to the contrasting yield production between treatment S2 and S4, it would be really interesting to test the effects of different defoliation numbers (one, three or more leaves removal) upper the ear leaves. That would be conducive to fully elucidate the effects of leaf removal in maize.
Based on our study, we demonstrated that different defoliation levels had contrastive effects on maize. The canopy light transmission rate was strengthened and proteins related to photosynthetic electron-transfer reaction were up-regulated significantly for treatment S2, which improved the leaf photosynthetic capacity, and obtained a higher grain yield consequently. In contrast, S4 decreased the grain yield and increased the expressions of proteins and genes associated with fatty acid metabolism. Besides, both S2 and S4 treatments exaggerated the defensive response of maize in physiological and proteomic level. Although further studies of leaf removal are required, the results in our study provide new insights to the effects of leaf removal in maize.
acyl carrier protein
chlorophyll a-b binding protein
- Ci :
intercellular CO2 concentration
days after silking
fatty acid synthetase
- gs :
light harvest antenna
oxygen-evolving enhancer protein
- Pn :
net photosynthetic rate
reactive oxygen species
tandem mass tags
Ciampitti IA, Vyn TJ. Maize nutrient accumulation and partitioning in response to plant density and nitrogen rate: II. Calcium, magnesium, and micronutrients. Agron J. 2013;105:783–95.
Testa G, Reyneri A, Blandino M. Maize grain yield enhancement through high plant density cultivation with different inter-row and intra-row spacings. Eur J Agron. 2016;72:28–37.
Deng J, Wang G, Morris E, Wei X, Li D, Chen B, Zhao C, Liu J, Wang Y. Plant mass-density relationship along a moisture gradient in north-West China. J Ecol. 2006;94:953–8.
Jin LB, Zhang JW, Bo LI, Cui HY, Dong ST, Liu P, Zhao B. Canopy structure and photosynthetic characteristics of high yield and high nitrogen efficiency summer maize. Scientia Agri Sinica. 2013;46:2430–9.
Borrás L, Maddonni GA, Otegui ME. Leaf senescence in maize hybrids: plant population, row spacing and kernel set effects. Field Crops Res. 2003;82:13–26.
Liu T, Huang R, Cai T, Han Q, Dong S. Optimum leaf removal increases nitrogen accumulation in kernels of maize grown at high density. Sci Rep. 2017;7:39601.
Hua HL. Zhao Q, Zhou Y, Li GS, Wang ZH, Bian Yl. Effect of photosynthetic characters of three ear leaves on stalk sugar content in maize. J Maize Sci. 2016;24:92–8.
Wei L, Xiong YC, Ma C, Zhang HQ, Shao Y, Li PF, Cheng ZG, Wang TC. Photosynthetic characterization and yield of summer corn (Zea mays L.) during grain filling stage under different planting pattern and population densities. Acta. Ecol Sinica. 2011;31:2524–31.
Ma SC, Xu BC, Li MF, Huang ZB. Ecological significance of redundancy in tillers of winter wheat (Triticum aestivum) and effect of reducing redundancy on water use efficiency. Acta Ecol Sinica. 2008;28:321–6.
Zhu GX, Midmore DJ, Radford BJ, Yule DF. Effect of timing of defoliation on wheat (Triticum aestivum ) in Central Queensland: 1. Crop response and yield Field Crops Res. 2004;88:211–26.
Gambín BL, Borrás L, Otegui ME. Source-sink relations and kernel weight differences in maize temperate hybrids. Field Crops Res. 2006;95:316–26.
Liu TN, Gu LM, Xu CL, Dong ST. Responses of group and individual leaf photosynthetic characteristics of two summer maize (Zea mays L.) to leaf removal under high plant density. Can J Plant Sci. 2014;94:1449–59.
Chen Z, Kolb TE, Clancy KM. Mechanisms of Douglas-fir resistance to western spruce budworm defoliation: bud burst phenology, photosynthetic compensation and growth rate. Tree Physiol. 2001;21:1159–69.
Vanderklein DW, Reich PB. European larch and eastern white pine respond similarly during three years of partial defoliation. Tree Physiol. 2000;20:283–7.
Van DHF, Stock WD. Regrowth of a semiarid shrub following simulated browsing: the role of reserve carbon. Funct Ecol. 1996;10:647–53.
Du JY, Zhou XC, Yang LR. Effect of damage of the photosynthesis organs in different wheat cultivars on grain yield per spike and their compensative effect. J Tritic. Crops. 2004;24:35–9.
Ruizr N, Ward D, Saltz D. Leaf compensatory growth as a tolerance strategy to resist herbivory in Pancratium sickenbergeri. Plant Ecol. 2008;198:19–26.
Hao MB, Wang KJ, Dong ST, Zhang JW, Li DH, Liu P, Yang JS, Liu JG. Leaf redundancy of high-yielding maize (Zea may L.) and its effects on maize yield and photosynthesis. Chin J Appl Ecol. 2010;21:344–50.
Liu TN, Xu CL, Gu LM, Dong ST. Effects of leaf removal on canopy apparent photosynthesis and individual leaf photosynthetic characteristics in summer maize under high plant density. Acta Agrono Sinica. 2014;40:143–53.
Wei SS, Wang XY, Dong ST. Effects of row spacing on canopy structure and grain-filling characteristics of high-yield summer maize. Chin J Appl Ecol. 2014;25:441–50.
Wei S, Wang X, Zhang J, Peng L, Zhao B, Geng L, Dong S. The role of nitrogen in leaf senescence of summer maize and analysis of underlying mechanisms using comparative proteomics. Plant Sci. 2015;233:72–81.
Dreccer MF, Grashoff C, Rabbinge R. Source-sink ratio in barley (Hordeum vulgare L.) during grain filling: effects on senescence and grain protein concentration. Field Crops Res. 1997;49:269–77.
Tollenaar M, Daynard TB. Effect of source-sink ratio on dry matter accumulation and leaf senescence of maize. Can J Plant Sci. 1982;62:855–60.
Fadeel AA. Location and properties of chloroplasts and pigment determination in roots. Physiol Plant. 2010;15:130–46.
Kaufman PB. Laboratory experiments in plant physiology. Kew Bull. 1975;12:49–50.
Wei S, Wang X, Zhang J, Liu P, Zhao B, Li G, Dong S. The role of nitrogen in leaf senescence of summer maize and analysis of underlying mechanisms using comparative proteomics. Plant Sci. 2015;233:72–81.
Wang X, Zhang X, Chen J, Wang X, Cai J, Zhou Q, Dai T, Cao W, Jiang D. Parental drought-priming enhances tolerance to post-anthesis drought in offspring of wheat. Front Plant Sci. 2018;9.
Eifediyi K, Remison SU. Effects of different organic materials and NPK fertilizer on the performance of maize. Poljoprivreda. 2014;20:16–21.
Wei S, Wang X, Zhu Q, Jiang D, Dong S. Optimising yield and resource utilisation of summer maize under the conditions of increasing density and reducing nitrogen fertilization. Sci Nat. 2017;104:86.
Weaver LM, Amasino RM. Senescence is induced in individually darkened Arabidopsis leaves, but inhibited in whole darkened plants. Plant Physiol. 2001;127:876–86.
Liu T, Gu L, Dong S, Zhang J, Liu P, Zhao B. Optimum leaf removal increases canopy apparent photosynthesis, 13C-photosynthate distribution and grain yield of maize crops grown at high density. Field Crops Res. 2015;170:32–9.
Wang X, Huang M, Zhou Q, Cai J, Dai T, Cao W, Jiang D. Physiological and proteomic mechanisms of waterlogging priming improves tolerance to waterlogging stress in wheat (Triticum aestivum L.). Environ Exp Bot. 2016;132:175–82.
Fromme P, Melkozernov A, Jordan P, Krauss N. Structure and function of photosystem I: interaction with its soluble electron carriers and external antenna systems. FEBS Lett. 2003;555:40–4.
Mochizuki N, Brusslan JA, Larkin R, Nagatani A, Chory J. Arabidopsis genomes uncoupled 5 (GUN5) mutant reveals the involvement of mg-chelatase H subunit in plastid-to-nucleus signal transduction. P Natl Acad Sci USA. 2001;98:2053–8.
Yokono M, Takabayashi A, Akimoto S, Tanaka A. A megacomplex composed of both photosystem reaction centres in higher plants. Nat Commun. 2015;6:6675.
Yamori W, Shikanai T. Physiological functions of cyclic electron transport around photosystem I in sustaining photosynthesis and plant growth. Annu Rev Plant Biol. 2016;67:81–106.
Staneloni RJ, Rodriguez-Batiller MJ, Casal JJ. 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. 2008;1:75–83.
Hansson A, Jensen PE. Chlorophyll limitation in plants remodels and balances the photosynthetic apparatus by changing the accumulation of photosystems I and II through two different approaches. Physiol Plant. 2009;135:214–28.
Amerongen HV, Croce R. Light harvesting in photosystem II. Photosynth Res. 2013;116:251–63.
Umena Y, Kawakami K, Shen JR, Kamiya N. Crystal structure of oxygen-evolving photosystem II at a resolution of 1.9 Å. Nature. 2011;473:55–60.
Bricker TM, Roose JL, Fagerlund RD, Frankel LK, Eaton-Rye JJ. The extrinsic proteins of photosystem II. Biochim Biophys Acta. 2012;1817:121–42.
Boehm M, Romero E, Reisinger V, Yu J, Komenda J, Eichacker LA, Dekker JP, Nixon PJ. Investigating the early stages of photosystem II assembly in Synechocystis sp. PCC 6803: isolation of CP47 and CP43 complexes. J Biol Chem. 2011;286:14812–9.
Mayes SR, Dubbs JM, Vass I, Hideg E, Nagy L, Barber J. Further characterization of the psbH locus of Synechocystis sp. PCC 6803: inactivation of psbH impairs QA to QB electron transport in photosystem 2. Biochemistry. 1993;32:1454–65.
Wang L, Pan D, Li J, Tan F, Hoffmann-Benning S, Liang W, Chen W. Proteomic analysis of changes in the Kandelia candel chloroplast proteins reveals pathways associated with salt tolerance. Plant Sci. 2015;231:159–72.
Pang Q, Chen S, Dai S, Chen Y, Wang Y, Yan X. Comparative proteomics of salt tolerance in Arabidopsis thaliana and Thellungiella halophila. J Proteome Res. 2010;9:2584–99.
Mishra S. “Cytochrome B”- analysis of hydrophobicity, surface accessibility, antigenicity and prediction of MHC I and MHC II binders from Dracunculiasis. Med Chem. 2016;6:041–6.
Roman K, Ondřej S, Lukáš N, René L, Ivo C, Egbert JB, Marek Š, Petr I. Structural characterization of a plant photosystem I and NAD(P)H dehydrogenase supercomplex. Plant J. 2014;77:568–76.
Amunts A, Drory O, Nelson N. The structure of a plant photosystem I supercomplex at 3.4 Å resolution. Nature. 2007;447:58–63.
Silveira JA, Carvalho FE. Proteomics, photosynthesis and salt resistance in crops: an integrative view. J Proteome. 2016;143:24–35.
Hiyama T. Isolation of photosystem I particles from spinach. Methods Mol Biol. 2004;274:11–7.
Nishiyama Y, Allakhverdiev S, Murata N. A new paradigm for the action of reactive oxygen species in the photoinhibition of photosystem II. Biochim Biophys Acta. 2006;1757:742–9.
James DW, Lim E, Keller J, Plooy I, Ralston E, Dooner HK. Directed tagging of the Arabidopsis fatty acld elongatlon1 (FAE1) gene with the maize transposon activator. Plant Cell. 1995;7:309–19.
Lee WM, Elliott JE, Brownsey RW. Inhibition of acetyl-CoA carboxylase isoforms by pyridoxal phosphate. J Biol Chem. 2005;280:41835–43.
Yin J, Straight PD, Hrvatin S, Dorrestein PC, Bumpus SB, Jao C, Kelleher N, Kolter R, Walsh CT. Genome-wide high-throughput mining of natural-product biosynthetic gene clusters by phage display. Cell Chem Bio. 2007;14:303–12.
Van Lanen S, Dorrestein P, Christenson S, Liu W, Ju J, Kelleher N, Shen B. Biosynthesis of the beta-amino acid moiety of the enediyne antitumor antibiotic C-1027 featuring beta-amino acyl-S-carrier protein intermediates. J Am Chem Soc. 2005;127:11594–5.
Ruiz J, Rivero R, Romero L, Pc. Response of phenolic metabolism to the application of carbendazim plus boron in tobacco. Physiol Plant 1999; 106: 151–157.
Camera SL, Gouzerh G, Dhondt S, Hoffmann L, Fritig B, Legrand M, Heitz T. Metabolic reprogramming in plant innate immunity: the contributions of phenylpropanoid and oxylipin pathways. Immunol Rev. 2004;198:267–84.
Maher EA, Lamb CJ. Increased disease susceptibility of transgenic tobacco plants with suppressed levels of preformed phenylpropanoid products. P Natl Acad Sci USA. 1994;91:7802–6.
Yao K, De LV, Brisson N. Creation of a metabolic sink for tryptophan alters the phenylpropanoid pathway and the susceptibility of potato to phytophthora infestans. Plant Cell. 1995;7:1787–99.
Wang YD, Wang X, Wong YS. Proteomics analysis reveals multiple regulatory mechanisms in response to selenium in rice. J Proteome. 2012;75:1849–66.
Christensen AB, Cho BH, Næsby M, Gregersen PL, Brandt J, Madriz-Ordeñana K, Collinge DB, Thordal-Christensen H. The molecular characterization of two barley proteins establishes the novel PR-17 family of pathogenesis-related proteins. Mol Plant Pathol. 2002;3:135–44.
Liu JJ, Ekramoddoullah AKM. The family 10 of plant pathogenesis-related proteins: their structure, regulation, and function in response to biotic and abiotic stresses. Physiol Mol Plant Pathol. 2006;68:3–13.
Suzuki N, Koussevitzky S, Mittler R, Miller G. ROS and redox signalling in the response of plants to abiotic stress. Plant Cell Environ. 2011;35:259–70.
Torres MA. ROS in biotic interactions. Physiol Plant. 2009;138:414–29.
Durand TC, Sergeant K, Renaut J, Planchon S, Hoffmann L, Carpin S, Label P, Morabito D, Hausman JF. Poplar under drought: comparison of leaf and cambial proteomic responses. J Proteome. 2011;74:1396–410.
Shin SY, Lee HS, Kwon SY, Kwon ST, Kwak SS. Molecular characterization of a cDNA encoding copper/zinc superoxide dismutase from cultured cells of Manihot esculenta. Plant Physiol Biochem. 2005;43:55–60.
Wang HB, Liu HY, Dai JR, Feng DR, Bing L, Wang JF. Characterization of a novel plantain Asr gene, MpAsr, that is regulated in response to infection of fusarium oxysporum f. Sp. Cubense and abiotic stresses. J Integr Plant Biol. 2010;52:315–23.
Schneider A, Gebhardt C. Expression patterns and promoter activity of the cold-regulated gene Ci21a of potato. Plant Physiol. 1997;113:335–45.
Kalifa Y, Perlson E, Gilad A, Konrad Z, Scolnik PA, Bar-Zvi D. Over-expression of the water and salt stress-regulated Asr1 gene confers an increased salt tolerance. Plant Cell Environ. 2004;27:1459–68.
Yang CY, Wang CSA. Lily ASR protein involves abscisic acid signaling and confers drought and salt resistance in Arabidopsis. Plant Physiol. 2005;139:836–46.
Dai JR, Liu B, Feng DR, Liu HY, He YM, Qi KB, Wang HB, Wang JF. MpAsr encodes an intrinsically unstructured protein and enhances osmotic tolerance in transgenic Arabidopsis. Plant Cell Rep. 2011;30:1219–30.
Trudel J, Grenier J, Potvin C, Asselin A. Several thaumatin-like proteins bind to beta-1,3-glucans. Plant Physiol. 1998;118:1431–8.
The authors gratefully acknowledge Deyang Shi, Guanghao Li, Xiulin Wang, Meifeng Zeng and Ting Liang for helping us to perform the study and acknowledge Qicen Zhu for helping draft part of the manuscript. We would also like to thank Mr. Mingcai Xu and all the workers in our farm for managing the field experiments.
This work was funded by the national key research and development program of China (no. 2016YFD0300308), the Natural Science Foundation of Jiangsu Province (no. BK20170720), the China Postdoctoral Science Foundation funded project (no. 2017 M611832), the Postdoctoral Science Foundation funded project of Jiangsu Province (no. 1701040A), the National Natural Science Foundation of China (no. 31171497), the National Basic Research Program of China (973 Program, no. 2011CB100105), the National Food Science and Technology of High-yield Program of China (no. 2011BAD16B09), the Special Fund for Agro-scientific Research in the Public Interest of China (no. 20120306) and supported by 111 Project (B16026).
Availability of data and materials
The resulting MS/MS data were processed using Mascot search engine (v.2.3.0). Tandem mass spectra were searched against uniprot Zea mays database (http://www.uniprot.org). All other additional data generated or analyzed during this study are included in this published article and its supplementary information files.
Ethics approval and consent to participate
The variety used in this study (Denghai 618), is a high-yield and density-tolerance variety grown extensively in North China. The maize seeds were brought from Shandong Denghai Seeds Co., Ltd. As our experiment involves neither transgenic materials nor technology, it does not require ethical approval. The experimental research on plants performed in this study complies with institutional, national and international guidelines. The field study was conducted in accordance with local legislation.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Primers used in quantitative RT-PCR in this study. (DOCX 31 kb)
Figure S1. The Pearson correlation analysis of the three replicates of each treatment. S0 refers to control (no leaf removal); S2 and S4 refer to the removal of two and four uppermost leaves, respectively. (TIF 1038 kb)
Table S2. Differences in protein abundances between the two-leaf removal treatment (S2) and the control (S0). (DOCX 85 kb)
Table S3. Differences in protein abundances between the four-leaf removal treatment (S4) and the control (S0). (DOCX 48 kb)
Figure S2. Gene ontology (GO) classification of differentially accumulated proteins. (A) Up-regulated proteins with S2 treatment compared to S0 treatment; (B) down-regulated proteins with S2 treatment compared to S0 treatment; (C) up-regulated proteins with S4 treatment compared to S0 treatment; and (D) down-regulated proteins with S4 treatment compared to S0 treatment. S0 refers to control (no leaf removal); S2 and S4 refer to the removal of two and four uppermost leaves, respectively. (PDF 199 kb)
Figure S3. Effects of leaf removal on relative expression of photosynthesis related proteins (A), the corresponding encoding genes (B) in S2 and relative expression of fatty acid metabolism related proteins (C) and the encoding genes (D) in S4, compared to S0 respectively. The gene candidates are selected by proteins which accumulate in photosynthesis and fatty acid biosynthetic process terms in Fig. 4. S0 refers to no leaf removal (control); S2 and S4 refer to the removal of two or four uppermost leaves, respectively. Data are means ± SE (n = 3). * indicates the significant difference at P ≤ 0.05 level. (PDF 57 kb)
Figure S4. Effects of leaf removal on relative expression of defense related proteins (A) and the encoding genes (B) in S2 and S4 compared to S0. S0 refers to no leaf removal (CK); S2 and S4 refer to the removal of two or four leaves, respectively. Data are means ± SE (n = 3). Different lowercase letters indicate the significant difference at P ≤ 0.05 level. (PDF 16 kb)