- Open Access
Proteomics reveals the effects of drought stress on the kernel development and starch formation of waxy maize
BMC Plant Biology volume 21, Article number: 434 (2021)
Kernel development and starch formation are the primary determinants of maize yield and quality, which are considerably influenced by drought stress. To clarify the response of maize kernel to drought stress, we established well-watered (WW) and water-stressed (WS) conditions at 1–30 days after pollination (dap) on waxy maize (Zea mays L. sinensis Kulesh).
Kernel development, starch accumulation, and activities of starch biosynthetic enzymes were significantly reduced by drought stress. The morphology of starch granules changed, whereas the grain filling rate was accelerated. A comparative proteomics approach was applied to analyze the proteome change in kernels under two treatments at 10 dap and 25 dap. Under the WS conditions, 487 and 465 differentially accumulated proteins (DAPs) were identified at 10 dap and 25 dap, respectively. Drought induced the downregulation of proteins involved in the oxidation–reduction process and oxidoreductase, peroxidase, catalase, glutamine synthetase, abscisic acid stress ripening 1, and lipoxygenase, which might be an important reason for the effect of drought stress on kernel development. Notably, several proteins involved in waxy maize endosperm and starch biosynthesis were upregulated at early-kernel stage under WS conditions, which might have accelerated endosperm development and starch synthesis. Additionally, 17 and 11 common DAPs were sustained in the upregulated and downregulated DAP groups, respectively, at 10 dap and 25 dap. Among these 28 proteins, four maize homologs (i.e., A0A1D6H543, B4FTP0, B6SLJ0, and A0A1D6H5J5) were considered as candidate proteins that affected kernel development and drought stress response by comparing with the rice genome.
The proteomic changes caused by drought were highly correlated with kernel development and starch accumulation, which were closely related to the final yield and quality of waxy maize. Our results provided a foundation for the enhanced understanding of kernel development and starch formation in response to drought stress in waxy maize.
Drought is the most serious environmental stress that obstructs crop growth and production worldwide and become more frequent due to climate change . Maize (Zea mays L.) is one of the most important cereal crops, and maize production is frequently compromised because of the increasing frequency and intensity of drought . Drought stress (DS) has drastically reduced maize production in many regions around the world, such as China and the United States [2, 3]. DS leads to abnormal embryo development and decrease kernel numbers, grain yield and quality, especially during the critical period of maize kernel formation [4, 5]. Thus, investigating the mechanism of drought response at the maize grain-filling stage is helpful for improving maize yield and quality under DS.
The development of endosperm (i.e., early development, differentiation, and maturation) has a remarkable influence on grain weight and quality and is highly sensitive to DS . DS reduces cell viability and accelerates nuclear deformation and programmed cell death in wheat endosperm cells . DS occurs during early endosperm development and inhibits kernel growth by decreasing endosperm cell division . In the endosperm of crops, the primary function of endosperm cells is to synthesize starch and store proteins . Maize endosperm, which is the foremost kernel component that determines grain weight and quality, contains about 75% starch. However, the biosynthesis and accumulation of starch are seriously threatened by DS . For instance, DS has a negative effect on the development of endosperm starch granules and the composition and physicochemical properties of starches . Studies reported that DS during the grain-formation stage reduced the number of endosperm cells and starch granules, but increased the starch granule size [5, 12]. Water shortages accelerate grain-filling and ripening periods and reduce the total starch accumulation, which is directly correlated with grain productivity [13, 14]. During grain development, DS impedes starch synthesis by reducing the enzymatic activities related to starch biosynthesis, which ultimately leading to a decline in grain yield [15, 16]. These enzymes mainly include: (1) Sucrose phosphate synthase (SPS) is a key enzyme in the regulation of sucrose biosynthesis in source tissues . (2) Sucrose synthase (SUS) catalyzes the cleavage of sucrose, which is the first step in the conversion of sucrose to starch. (3) Adenosine diphosphate-glucose pyrophosphorylase (AGPase) catalyzes ADP-Glucose, which is the rate-limiting enzyme in starch biosynthesis. (4) Soluble starch synthase (SSS) catalyzes the glucosyl moiety of ADP-Glucose to synthesize amylose and amylopectin. (5) Starch branching enzyme (SBE) catalyzes α-(1,6)-linkages within the polymer to form the branched structure of amylopectin . The study of how endosperm and starch development respond to DS is crucial to ensure the high yield and quality of maize under adverse conditions. However, our understanding of drought response mechanisms during the grain-filling stage remains unclear.
Recent advances in high-throughput technologies enable the quantitation of abundance changes in transcriptomics, proteomics, and other omics under DS . Generally, many studies on the response of maize to DS focus on large-scale transcriptomic analyses [20,21,22,23]. However, transcriptome profiling has limitations because the changes in gene expression levels do not correspond directly to protein expression levels due to post-transcriptional and post-translational modifications . Proteomics, as a research of post-genome era, is a crucial link between transcriptomics and metabolomics . High-throughput proteomics has become a powerful tool that can provide new insights into how plants respond to DS at the protein level . Proteins are important for plant DS response because of their direct effect on metabolism and other cellular processes . For example, alterations in proteins related to the metabolism of carbohydrates and antioxidation and defense systems are found to be characteristic features between drought-tolerant and drought-sensitive wheat . In maize, the changes in seeding leaf protein abundance are consistent with the observed phenotypes of YE8112 (drought-tolerant line) and Mo17 (drought-sensitive line) under DS . Moreover, comparative proteomics analysis is an effective strategy to identify pivotal candidate functional proteins and regulatory networks under DS . In maize, studies on the proteome response to DS on kernel development and starch formation are limited.
Thus, in the present study, we performed the comparative proteomic analysis of waxy maize hybrid kernels under DS during the grain filling and analyzed the dynamic changes of kernel development, starch accumulation, and related enzymatic activities. The objectives of the present study were to identify candidate proteins that might play important roles in the kernel development, starch accumulation, and drought stress response. The results may enhance our understanding of the regulatory networks of endosperm development and starch formation in response to DS. Furthermore, drought responsive proteins identified herein could be harnessed for improving drought tolerance breeding in waxy maize yield and quality.
Dynamic changes in kernel development under drought stress
Ears were harvested and analyzed at six stages after pollination (from 5 to 30 dap with five days apart). Kernel fresh and dry weights gradually increased with kernel development under well-watered (WW) and water-stressed (WS) conditions. However, kernel development was significantly inhibited under WS conditions (P < 0.01, Fig. 1A and B). From 25 to 30 dap, the increase rates of kernel fresh and dry weights under WW conditions were significantly higher than those under WS conditions. The kernel water content gradually decreased with kernel development, and the decline rate of WS was higher than that of WW during the whole grain-filling stage (Fig. 1C). These results indicated that WS conditions accelerated grain maturity, shortened grain filling, and reduced grain weight.
Dynamic changes in the starch accumulation and enzymatic activities under drought stress
With the progress of grain filling, the accumulation of starch in endosperm cells increased gradually. The starch content increased sharply from 10 to 20 dap under both water levels. Our results showed no difference in the starch content under two treatments at 10 dap, but the starch content under WS treatment was significantly higher than that under WW treatment at 15 dap (Fig. 1D). Compared with WW, WS caused a significant decrease in starch accumulation from 20 to 30 dap (Fig. 1D). The starch in endosperm cells at 10–30 dap was observed and revealed by scanning electron microscopy to evaluate the morphological features and accumulation of starch granules under WW and WS conditions (Fig. 1E). The starch granule was irregularly shaped due to crowding under WW treatment at 25 and 30 dap but mostly spherical under WS (Fig. 1E). In addition, the activities of starch synthetic enzymes (i.e., AGPase, SUS, SBE, SDBE, SPS, and SSS) increased initially and declined thereafter with kernel development (Fig. 2). Furthermore, all enzymatic activities were significantly reduced under WS treatment during grain filling (P < 0.01). Therefore, DS significantly inhibited the kernel development, starch synthetic activities, and starch accumulation during grain filling.
Quantitative proteomic analysis
A total of 6,769 proteins were identified from waxy maize kernels by using DIA mass spectroscopy (Table S1). The principal component analysis (PCA) was performed to characterize all proteins in the three replicates at 10 dap and 25 dap under the WW and WS treatments. In the PCA plot, a clear separation was observed between WW and WS treatments at 10 dap and 25 dap, suggesting that the proteome was distinct between the two treatments (Figure S1). To identify the differentially accumulated proteins (DAPs) among treatments, two fold change and P < 0.05 were chosen as the criteria. Compared with the WW treatment, the numbers of DAPs for WS treatment were 487 (320 upregulated and 167 downregulated) at 10 dap, and 465 (242 upregulated and 223 downregulated) at 25 dap (Fig. 3A, Table S2). At 10 dap, 19 and 2 transcription factors (TFs) were identified in upregulated and downregulated proteins, respectively. Moreover, 10 and 15 TFs were found in upregulated and downregulated proteins, respectively, at 25 dap (Fig. 3B). These TFs belonged to 31 families, and members of the bHLH, bZIP, and C2H2 TF families were highly represented among the DAPs (Fig. 3B). In addition, we randomly chose ten genes for expression validation by using qRT-PCR. The results showed that seven genes expression level displayed the same trend with the abundance of the corresponding protein species, three genes showed opposite trends to the abundance of their corresponding proteins (Figure S2). Thus, changes in mRNA expression might not always adequately reflect corresponding protein levels.
Among all DAPs, 54 (38 upregulated and 16 downregulated) and 45 (24 upregulated and 21 downregulated) DAPs with experimentally confirmed functions [32, 33] were identified at 10 dap and 25 dap, respectively (Figure S3). These DAPs that responded to DS had different functions. For instance, A0A1D6JMR5 and A0A1D6Q9W2, which related to catalase and glutamine synthetase, respectively, were downregulated at 10 dap. Abscisic acid stress ripening 1 (B4FKG5) and lipoxygenase (A0A1D6DVH5) were downregulated at 25 dap. In addition, several DAPs, which involved in endosperm development (i.e., A0A1D6H9T3, A0A1D6GRB4, B4F832, and A0A1D6PZG8) and starch accumulation (i.e., A0A1D6EBS7, A0A1D6NUU8, A0A1D6PXQ9, and A0A1D6LVT4), were upregulated at 10 dap.
KEGG pathway and GO annotation analyses of DAPs
KEGG pathways were classified into hierarchical categories according to the website (http://www.kegg.jp/kegg/kegg1.html). In present study, KEGG pathway analysis showed that DAPs at 10 dap were significantly enriched in the 11 KEGG pathway (such as “Plant hormone signal transduction”, “Starch and sucrose metabolism”, “Biosynthesis of secondary metabolites”, etc.), but only one significantly enriched pathway, i.e., “Ribosome” at 25 dap (P < 0.01; Figure S4). In addition, a high proportion of the upregulated DAPs at 10 dap were enriched in “Phenylalanine, tyrosine, and tryptophan biosynthesis” and “Plant hormone signal transduction”. A large number of downregulated DAPs at 10 dap were enriched in the “Metabolic pathways”, “Phenylalanine metabolism”, “Flavonoid biosynthesis”, “Thiamine metabolism”, and “Phenylpropanoid biosynthesis”. However, only one pathway, i.e., “Ribosome”, was enriched among the upregulated proteins at 25 dap (Fig. 3C).
We also analyzed the GO terms represented by all DAPs at two stages. The enriched GO terms at 10 dap included mainly “cellular polysaccharide metabolic process”, “polysaccharide metabolic process”, and “cellular carbohydrate metabolic process”, etc.; and at 25 dap included mainly “organonitrogen compound metabolic process”, “cellular protein metabolic process”, and “cellular nitrogen compound metabolic process”, etc. (Figure S5). These results demonstrated that the DAPs at 10 dap and 25 dap were involved in different biological processes. Additionally, in order to clarify the biological processes involved in up- and down-regulated DAPs, we performed GO analysis for these DAPs. For the upregulated DAPs at 10 dap, three significantly enriched GO terms, including “chloroplast”, “vitamin binding”, and “lyase activity”, had the most DAPs (Fig. 4A). In addition, the “oxidation–reduction process”, “oxidoreductase activity”, and “catalytic activity” had the most DAPs in downregulated DAPs at 10 dap (Fig. 4B). In the GO enrichment analysis of upregulated proteins at 25 dap, significant GO terms, i.e., the “organonitrogen compound metabolic process”, “cytoplasm”, and “cell”, contained more DAPs (Fig. 4C). Significant GO terms, i.e., “oxidation–reduction process” and “oxidoreductase activity”, had most DAPs in downregulated DAPs at 25 dap (Fig. 4D). This result suggested that the oxidation–reduction process and oxidoreductase activity were enriched in downregulated DAPs at the two stages and might negatively affect the drought tolerance of kernels under DS.
DAPs between 10 and 25 dap
Comparison analysis was conducted for the upregulated and downregulated DAPs to discover common DAPs in response to DS among different stages. Seventeen common DAPs were sustained in the upregulated group, and 11 common DAPs were consistent in the downregulated group (Table 1). We found that the peroxidase (B4FBC8) was downregulated in the two stages. Unfortunately, these DAPs have not been reported in maize. Thus, we performed the homology comparison of these proteins through the protein–protein BLAST method and found four homologous genes with known functions in rice, including A0A1D6H543 (homolog of PDIL1-1), B4FTP0 (homolog of OsADF2), B6SLJ0 (homolog of MODD), and A0A1D6H5J5 (homolog of BS1). These genes have been shown to involved in cell wall (BS1), endosperm development (PDIL1-1), and DS response (OsADF2 and MODD, Table 1) in rice. Thus, these homologous genes may have similar functions in maize.
Effects of DS on the accumulation level of proteins involved in maize endosperm and starch biosynthesis
The main component of maize kernels is the endosperm, which includes starchy endosperm (SE), aleurone (AL), basal endosperm transfer layer (BETL), and embryo-surrounding region (ESR) cells. Moreover, the endosperm adjacent to scutellum (EAS) cells appears around 9 dap and persists throughout embryo growth . Considering that DS affects the development of endosperm, we specifically explored the accumulation patterns of proteins related to endosperm development under WS conditions. Therefore, we collected 33 genes related to different cells of the endosperm from the report of Dai et al.  (Fig. 5). In the present study, A0A1D6EAS6 (O11) and A0A1D6N2E7 (ZmICEa), as central hub factors regulating endosperm development, were upregulated under WS treatment at 10 dap (Fig. 5). The proteins of Q84LS4, A0A1D6I8V4, A0A060D1R6, A0A1D6PZG8, A0A1D6M414, and K7TFF5, which involved in AL, EAS, SE, and BETL, respectively, were upregulated under WS treatment at 10 dap, but these proteins were reversed at 25 dap (Fig. 5). Proteins of A0A1D6MEW5, A0A0B4J2Y9, B4FQK1, B7ZXT3, A0A1D6M2N3, A0A1D6DXK1, and B4FQK1 were downregulated under WS treatment at 25 dap (Fig. 5).
In addition, we identified 21 proteins involved in starch biosynthesis and analyzed the accumulation patterns through proteomic data. Among these 21 proteins, 8 were differentially accumulated at 10 dap and 25 dap. The accumulation levels of A0A1D6NT44 (sh1), A0A1D6PXQ9 (su1), A0A1D6NHZ0 (sh2), Q947B9 (agpsl1), A0A1D6NUU8 (ss1), and A0A1D6EBS7 (sbe4), which associated with SUS, AGPase, SSS, and SBE, respectively, were upregulated under WS treatment at 10 dap. The proteins A0A1D6HWJ3 (agpsl2) and A0A1D6JEU3 (ss2), which associated with AGPase and SSS, were upregulated by DS at 25 dap. Only one protein, i.e., A0A1D6IAT2 (which related to UGPase), was downregulated under WS treatment at 25 dap (Figure S6). These results indicated that the changes in the proteome might accelerate endosperm development and kernel maturation and increase the accumulation of starch in the early stages of kernel development under DS.
Dynamic response of kernel development to drought stress during grain filling
In the context of climate change, drought stress will have a profound impact on maize productivity in the future [40, 41]. During the whole growth period of maize, the grain-filling stage is important for kernel development and quality formation, and DS at this stage results in kernel abortion and yield reduction . In the present study, DS after pollination inhibited kernel development (Fig. 1). From 25 to 30 dap, the kernel increase rate under WW treatment was significantly higher than that under WS treatment, indicating that kernel is close to maturity under WS. In addition, the changes in the kernel water content proved this point, and the kernel water content under WS treatment had a fast decline rate (Fig. 1C). Endosperm development affects kernel weight and quality. DS may have affected the development of endosperm cells and reduced the number of endosperm nuclei, resulting in obstructed grain filling and decreased yield [6, 7].
Thus, the molecular regulation mechanism of endosperm development is crucial for the formation of grain yield under DS. In this study, we collected 33 proteins related to endosperm development and analyzed their protein accumulation patterns under WW and WS treatments through proteomics. The upregulation of nine proteins that related to AL, EAS, SE, and BETL under WS treatment at 10 dap might jointly regulate the rapid development of kernels under WS treatment. For example, Feng et al.  reported that ZmICEa (which encodes A0A1D6N2E7) interacts with O11 (which encodes A0A1D6EAS6), and they coregulate endosperm development and stress response. ZmNAC128 (which encodes A0A060D1R6) positively regulates the accumulation of starch and protein in grains . The TF ZmbZIP22 (which encodes A0A1D6I8V4) regulates the 27 kD γ-zein gene expression and affects protein body initiation during endosperm development . Moreover, the accumulation of seven other proteins were downregulated under DS at 25 dap, which may be consistent with the significant reduction in the fresh and dry weights of grains (P < 0.01). Feng et al.  reported that O11 is a central regulator linking maize endosperm development, which regulates the expression of ZmYODA (which encodes A0A0B4J2Y9) specifically expressed in ESR. The maize VKS1 (which encodes A0A1D6DXK1) regulates mitosis and cytokinesis during early-endosperm development, and its loss of function will lead to reduced cell proliferation and grain size . These proteins may affect the rate of grain filling and shorten the grain-filling and ripening periods by affecting the endosperm development. Overall, these proteins may coregulate kernel development under WS conditions.
Effects of drought stress on starch formation during the grain-filling stage
The most important nutrient stored in the endosperm of maize is starch, accounting for about 75% of the dry weight of grain. The morphological and physicochemical properties of starch granules affect the grain quality . Numerous studies showed that the biosynthesis and accumulation of starch are remarkably influenced by DS [16, 47]. DS inhibits the activities of starch synthetic enzymes (such as AGPase, SSS, SBE, and GBSS), which suppress grain development and ultimately reduce grain weight [16, 48]. In present experiment, we found that the activities of AGPase, SUS, SBE, SDEB, and SSS were reduced by WS treatment during the grain-filling stage (Fig. 2). In addition, WS treatment decreased the starch accumulation from 20 to 30 dap, and the spherical shapes of starch granules were unchanged after the WS treatment (Fig. 1E). Li et al.  reported a similar phenomenon for starch granules under WS treatment during grain filling. In the endosperm development, the expression of starch biosynthesis proteins is also reported to be affected by DS [50, 51]. In this study, 21 proteins related to starch biosynthesis were found using proteomics (Figure S1). Compared with WW conditions, six DAPs (i.e., A0A1D6NT44, A0A1D6PXQ9, A0A1D6NHZ0, Q947B9, A0A1D6NUU8, and A0A1D6EBS7) were upregulated under WS treatment at 10 dap, and this finding might be related to the upregulation of proteins that regulated SE in kernel (Fig. 5). Notably, no difference in starch content was observed between WW and WS treatments at 10 dap, but the starch content under WS treatment was significantly higher than that under WW treatment at 15 dap (Fig. 1D). Our results showed that DS accelerated the accumulation of starch at early grain development stage. Therefore, these enzymes and proteins had coordinated participation in the complex starch biosynthesis in maize kernel under DS.
Proteomic analysis of drought stress response during the grain filling
Understanding how kernels respond to drought stress at the molecular level is important for improving the maize yield and quality. In this study, proteomic changes at 10 dap and 25 dap were surveyed to gain insight into the protein accumulation patterns in kernel in response to DS (Fig. 3). We identified 487 and 465 DAPs at 10 dap and 25 dap, respectively. Among these DAPs, 21 and 25 TFs were identified at 10 dap and 25 dap, respectively (Fig. 3). In previous studies, TFs play important roles in DS response and plant development . For example, the rice gene OsICE1 (OsbHLH002), which is orthologous to A0A1D6N2E7 (upregulated at 10 dap), has been reported to improve tolerance to abiotic stresses in tobacco . The accumulation of A0A1D6IHZ6 (an ortholog of OsWRKY78), a stem elongation and seed development regulator in rice , was downregulated at 25 dap. In our data set, four DAPs (catalase, glutamine synthetase, abscisic acid stress ripening 1, and lipoxygenase) were downregulated by DS (Figure S3). Catalase plays a key role not only in the defense against oxidative stress but also in governing cellular redox status and regulating cellular signaling . The ectopic expression of the maize Cat2 gene (A0A1D6JMR5) in tobacco plants changes the levels of catalase and increases the resistance to oxidative stress . The glutamine synthetase (GS) activity and expression are reported to be modulated in various plants to respond to DS [28, 57]. James et al.  previously reported that OsGS1;1 (a homolog of A0A1D6Q9W2), which overexpresses transgenic rice, enhances tolerance to osmotic stress and reduces ROS production. A previous study suggested that abscisic acid-, stress-, and ripening-induced (ASR) proteins are involved in the response of plant cell to various stressors . The overexpression of OsASR1 (a homolog of B4FKG5) in rice plants enhances salinity and drought tolerance and improves crop yields . As a multifunctional enzyme, lipoxygenases play an important role when plants are subjected to biotic and abiotic stresses . Consequently, catalase, glutamine synthetase, abscisic acid stress ripening 1, and lipoxygenase were downregulated under DS, which indicated that DS had a serious influence on kernel development.
The GO functional classification analysis revealed that upregulated and downregulated DAPs at the two stages were associated with a wide range of functions. The GO terms oxidation–reduction process (GO: 0,055,114) and oxidoreductase activity (GO: 0,016,672) were the most highly enriched among the downregulated DAPs at both stages (Fig. 4). These proteins were involved in the oxidation–reduction process, which also provide further information about the response of kernel development and starch formation responses to drought. In addition, 17 and 11 common DAPs were identified in upregulated and downregulated groups at 10 dap and 25 dap, respectively (Table 1). Among these data set, four proteins had homologous genes with known functions in rice (Table 1). Kim et al.  reported that PDIL1-1 (a homolog of A0A1D6H543) controls endosperm development through the regulation of the amount and composition of seed proteins in rice. A recent study found that OsADF2-OE (a homolog of B4FTP0) maintains a relatively high relative water content, plant height, membrane stability index, and dry weight under DS . MODD (a homolog of B6SLJ0) negatively regulates ABA signaling and drought tolerance in rice . The phenotypic analysis of the bs1 mutant revealed that BS1 (a homolog of A0A1D6H5J5) plays an important role in the maintenance of proper acetylation level on the xylan backbone . Consequently, these proteins may be important candidates for the endosperm development and DS response of the waxy maize kernel during grain filling.
In summary, drought stress during grain filling seriously affected the physiological and biochemical parameters related to waxy maize kernel development and starch formation. Comparative proteome analysis identified 487 and 465 DAPs at 10 dap and 25 dap under drought stress, respectively. Among them, some proteins associated with endosperm development and starch biosynthesis were upregulated at 10 dap, which could be responsible for the shortened grain filling, and increased starch accumulation in the early stages of kernel development. Furthermore, some of the DAPs that were positively related to stress responses were downregulated in the waxy maize kernels and might negatively affect the drought tolerance of kernels. Notably, we determined four promising candidate proteins for kernel development and drought tolerance for subsequent functional validation. Our findings provided novel insights into the proteomic level and physiological response mechanism of waxy maize kernel development and starch formation under drought stress.
Plant materials and DS treatment
A pot trial was performed at Yangzhou University experimental farm (Yangzhou, China) in 2018–2019. The seeds of waxy maize ‘Suyunuo5’ (provided by Jiangsu Yanjiang Institute of Agricultural Sciences, Nantong, China), a widely planted waxy maize variety in Southern China, was selected as experimental material. Seeds with the same size were sown in seedling trays. Three 7-day-old healthy seedlings were transplanted into pots (diameter = 38 cm, height = 43 cm) containing 30 kg loam soil obtained from the field. At the jointing stage, a single plant was left in the pot, and the plants in all pots grew to the same state. Plants were provided a basal dressing of 10 g per pot (commercial fertilizer, N:P2O5:K2O = 15%:15%:15%) at the transplanting date and a top dressing of 6.6 g per pot (commercial urea, N = 46%) at the jointing stage.
All plants were grown under natural conditions until pollination. After pollination, these plants were moved to a mobilizable transparent waterproof shed and treated with water until maturity. The soil relative moisture contents under well-watered (WW) and water-stressed (WS) conditions were 75%–80% and 50%–55%, respectively. The amount of water loss was calculated and resupplied by weighing at 09:00 every day. These treatments were applied from 1 to 30 days after pollination (dap). After 30 dap, they were rewatered normally until maturity.
Grain harvest and phenotypic identification
The ears of five plants were harvested at 5, 10, 15, 20, 25, and 30 dap. The kernels in the middle of the ear were used for analysis. The kernel dry weight was weighed after the kernel was dried in an oven at 80 °C for 24 h. The starch content was determined using the anthrone–sulfuric acid method with three biological replicates of each treatment . The kernel water content (%) was calculated as: Kernel water content (%) = (fresh weight – dry weight)/fresh weight×100.
Scanning electron micrographs of endosperm cells
The resin slicing method was used to cut the kernels at 10, 15, 20, 25, and 30 dap into slices with thickness of about 300 mm (n = 5 per stage per treatment). After the slices were plated with ion-sputtered gold, endosperm cells and starch granules were observed using the GeminiSEM 300 (Carl Zeiss, Jena, Germany) environmental scanning electron microscope.
Enzymatic activity assays for starch biosynthesis
Waxy maize kernels (middle position of ears) were stripped from ears at 5, 10, 15, 20, 25, and 30 dap and stored at a refrigerator maintained at − 80 °C. We assayed the activities of five starch biosynthesis enzymes, including SUS, AGPase, SSS, SBE, SPS and starch-debranching enzyme (SDBE), respectively, using the reagent kits ml10561, ml20341, ml10983, ml10724, ml056613 and ml10726 obtained from Shanghai Enzyme-linked Biotechnology Co., Ltd. All enzymatic activities were measured by enzyme-linked immunoassay, with three biological replicates. The OD value of SUS, AGPase, SSS, SBE, SPS and SDBE was measured with a microplate reader (Labsystems Multiskan, MS 352, Finland) at 450 nm and quantified using a standard curve.
Total proteins were extracted from the kernels under WW and WS treatments at 10 dap (10 days after treatment) and 25 dap (25 days after treatment) with three biological replicates in accordance, and the methods described by Zhu et al.  with minor modifications. Approximately 500 mg maize kernels were ground in liquid nitrogen, then dissolved with the lysis buffer and centrifuged at 25,000 g at 4 °C for 15 min. The supernatant was then transferred into a new tube and precipitated with 4 volumes of prechilled 10% TCA-acetone at –20 °C overnight and centrifuged at 25,000 g at 4 °C for 30 min. Pellets were washed by prechilled acetone three times. The protein dissolved with the lysis buffer and protein concentration was determined by the Bradford assay. Protein digestion was performed using the filter-aided sample preparation method with trypsin (Promega, Madison, WI) in 50 mM NH4HCO3 (Sigma, MO, USA) . The absorbance was determined at A280 nm by using the NanoDrop 2000 instrument (Thermo Fisher Scientific, Mass, USA) to measure the peptide concentration. Resulting peptides were then separated and analyzed using the Orbitrap Fusion Tribrid mass spectrometer (Thermo Fisher Scientific, Mass, USA) equipped with the EASY-nLC 1200 system (Thermo Fisher Scientific, Mass, USA).
Protein identification and bioinformatics analysis
The detailed protocol for data-dependent acquisition of library construction and data-independent acquisition (DIA) of peptide quantification were performed on the Orbitrap Fusion Tribrid mass spectrometer (Thermo Fisher Scientific, Mass, USA) equipped with the EASY-nLC 1200 system (Thermo Fisher Scientific, Mass, USA) in accordance with a previous method  with minor modifications. Samples (2 μg) were analyzed on the C18 trap column (150 μm I.D. × 25 cm, C18, 1.9 μm, 120 Å, Dr. Maisch GmbH). For the DIA analysis, the full scan was set at a resolution of 60,000 at m/z 350–1500. A total of 45 variable DIA windows were set for DIA acquisition ranging from m/z 200 to m/z 2000 with resolution of 30,000. Protein identification and quantification were finished using the Spectronaut pulsar X 12.0 (Biognosys) at default setting. In general, the raw spectral data were searched against the Uniprot (https://www.uniprot.org/, 2019.1.30) B73 complete proteome database (99,254 sequences). All identified proteins were filtered at 1% false discovery rates (FDRs). A DAP was defined with |log2 (fold change) |≥ 1 and FDR < 0.05. Functional annotations of DAPs were performed using the Gene Ontology (GO) annotation (http://www.geneontology.org/) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of DAPs were based on KEGG database (https://www.kegg.jp/kegg/kegg1.html) . The homologous genes of maize in rice (identity > 40% and coverage > 60%) were identified using the NCBI protein BLAST (https://blast.ncbi.nlm.nih.gov/Blast.cgi) .
We selected ten identified DAPs and analyzed the mRNA expression levels of corresponding genes by using quantitative real-time PCR (qRT-PCR). The qRT-PCR was performed on ABI ViiA™ 7 Real-Time PCR System (Applied Biosystems, CA, USA). Maize GAPDH was used as the internal reference gene . The relative expression levels were calculated by the 2−∆∆Ct method . The Primer Premier 6 software was used to design primers on the basis of cDNA sequences in the maizeGDB (https://www.maizegdb.org/, Table S3).
The statistical analyses of physiological data were performed with SPSS v19.0 using ANOVA, followed by LSD test to evaluate the significant differences at P < 0.01.
Adenosine diphosphate-glucose pyrophosphorylase
Basal endosperm transfer layer
Days after pollination
Differentially accumulated protein
Endosperm adjacent to scutellum
Kyoto Encyclopedia of Genes and Genomes
Principal component analysis
Quantitative real-time PCR
Starch branching enzyme
Sucrose phosphate synthase
Soluble starch synthase
Hu HH, Xiong LZ. Genetic engineering and breeding of drought-resistant crops. Annu Rev Plant Biol. 2014;65:715–41.
Lobell DB, Roberts MJ, Schlenker W, Braun N, Little BB, Rejesus RM, Hammer GL. Greater sensitivity to drought accompanies maize yield increase in the U.S. Midwest. Science. 2014;344(6183):516–5198.
Liu SL, Wu WB, Yang XG, Yang P, Sun J. Exploring drought dynamics and its impacts on maize yield in the Huang-Huai-Hai farming region of China. Clim Change. 2020;163(1):415–30.
Zinselmeier C, Jeong BR, Boyer JS. Starch and the control of kernel number in maize at low water potentials. Plant Physiol. 1999;121(1):25–35.
Wang LF, Yan YN, Lu WP, Lu DL. Application of exogenous phytohormones at silking stage improve grain quality under post-silking drought stress in waxy maize. Plants-Basel. 2021;10(1):48.
Begcy K, Walia H. Drought stress delays endosperm development and misregulates genes associated with cytoskeleton organization and grain quality proteins in developing wheat seeds. Plant Sci. 2015;240:109–19.
Li C, Li C, Wang BB, Zhang RQ, Fu KY, Gale WJ, Li CY. Programmed cell death in wheat (Triticum aestivum L.) endosperm cells is affected by drought stress. Protoplasma. 2018;255(4):1039–52.
Ober ES, Setter TL, Madison JT, Thompson JF, Shapiro PS. Influence of water deficit on maize endosperm development: enzyme activities and RNA transcripts of starch and zein synthesis, abscisic acid, and cell division. Plant Physiol. 1991;97(1):154–64.
Zheng YK, Wang Z. The cereal starch endosperm development and its relationship with other endosperm tissues and embryo. Protoplasma. 2015;252(1):33–40.
Yu XR, Li B, Wang LL, Chen XY, Wang WJ, Gu YJ, Wang Z, Xiong F. Effect of drought stress on the development of endosperm starch granules and the composition and physicochemical properties of starches from soft and hard wheat. J Sci Food Agr. 2016;96(8):2746–54.
Lu DL, Cai XM, Lu WP. Effects of water deficit during grain filling on the physicochemical properties of waxy maize starch. Starch-Stärke. 2015;67(7–8):692–700.
Singh S, Singh G, Singh P, Singh N. Effect of water stress at different stages of grain development on the characteristics of starch and protein of different wheat varieties. Food Chem. 2008;108(1):130–9.
Setter TL, Flannigan BA, Melkonian J. Loss of kernel set due to water deficit and shade in maize: carbohydrate supplies, abscisic acid, and cytokinins. Crop Sci. 2001;41(5):1530–40.
Yang JC, Zhang JH, Wang ZQ, Zhu QS. Activities of starch hydrolytic enzymes and sucrose-phosphate synthase in the stems of rice subjected to water stress during grain filling. J Exp Bot. 2001;52(364):2169–79.
Yang JC, Zhang JH, Wang ZQ, Xu GW, Zhu QS. Activities of key enzymes in sucrose-to-starch conversion in wheat grains subjected to water deficit during grain filling. Plant Physiol. 2004;135(3):1621–9.
Yang H, Gu XT, Ding MQ, Lu WP, Lu DL. Activities of starch synthetic enzymes and contents of endogenous hormones in waxy maize grains subjected to post-silking water deficit. Sci Rep-Uk. 2019;9(1):7059.
Geigenberger P, Reimholz R, Deiting U, Sonnewald U, Stitt M. Decreased expression of sucrose phosphate synthase strongly inhibits the water stress-induced synthesis of sucrose in growing potato tubers. Plant J. 1999;19(2):119–29.
Thitisaksakul M, Jiménez RC, Arias MC, Beckles MD. Effects of environmental factors on cereal starch biosynthesis and composition. J Cereal Sci. 2012;56(1):67–80.
Shanker AK, Maheswari M, Yadav SK, Desai S, Bhanu D, Attal NB, Venkateswarlu B. Drought stress responses in crops. Funct Integr Genomics. 2014;14(1):11–22.
Zhang QB, Liu H, Wu XL, Wang W. Identification of drought tolerant mechanisms in a drought-tolerant maize mutant based on physiological, biochemical and transcriptomic analyses. BMC Plant Biol. 2020;20(1):315.
Zhang XJ, Liu XY, Zhang DF, Tang HJ, Sun BC, Li CH, Hao LY, Liu C, Li YX, Shi YS, et al. Genome-wide identification of gene expression in contrasting maize inbred lines under field drought conditions reveals the significance of transcription factors in drought tolerance. PloS One. 2017;12(7):e0179477.
Guo J, Li CH, Zhang XQ, Li YX, Zhang DF, Shi YS, Song YC, Li Y, Yang DG, Wang TY. Transcriptome and GWAS analyses reveal candidate gene for seminal root length of maize seedlings under drought stress. Plant Sci. 2020;292:110380.
Wang BM, Liu C, Zhang DF, He CM, Zhang JR, Li ZX. Effects of maize organ-specific drought stress response on yields from transcriptome analysis. BMC Plant Biol. 2019;19(1):335.
Zhao Q, Zhang H, Wang T, Chen SX, Dai SJ. Proteomics-based investigation of salt-responsive mechanisms in plant roots. J Proteomics. 2013;82:230–53.
Tan BC, Lim YS, Lau SE. Proteomics in commercial crops. An overview. J Proteomics. 2017;169:176–88.
Su XY, Fan XC, Shao RX, Guo JM, Wang YC, Yang JP, Yang QH, Guo L. Physiological and iTRAQ-based proteomic analyses reveal that melatonin alleviates oxidative damage in maize leaves exposed to drought stress. Plant Physiol Bioch. 2019;142:263–74.
Wu S, Ning F, Zhang QB, Wu XL, Wang W. Enhancing omics research of crop responses to drought under field conditions. Front Plant Sci. 2017;8:174.
Cheng LX, Wang YP, He Q, Li HJ, Zhang XJ, Zhang F. Comparative proteomics illustrates the complexity of drought resistance mechanisms in two wheat (Triticum aestivum L.) cultivars under dehydration and rehydration. BMC Plant Biol. 2016;16(1):188.
Zenda TS, Liu ST, Wang XA, Jin HY, Liu G, Duan HJ. Comparative proteomic and physiological analyses of two divergent maize inbred lines provide more insights into drought-stress tolerance mechanisms. Int J Mol Sci. 2018;19(10):3225.
Deng X, Liu Y, Xu XX, Liu DM, Zhu GR, Yan X, Wang ZM, Yan YM. Comparative proteome analysis of wheat flag leaves and developing grains under water deficit. Front Plant Sci. 2018;9:425.
Kanehisa M, Goto S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000;28(1):27–30.
Schaeffer ML, Harper LC, Gardiner JM, Andorf CM, Campbell DA, Cannon EK, Sen TZ, Lawrence CJ. MaizeGDB: curation and outreach go hand-in-hand. Database (Oxford). 2011; 022.
Schnable JC, Freeling M. Genes identified by visible mutant phenotypes show increased bias toward one of two subgenomes of maize. PloS One. 2011;6(3):e17855.
Kim YJ, Yeu SY, Park BS, Koh HJ, Song JT, Seo HS. Protein disulfide isomerase-like protein 1–1 controls endosperm development through regulation of the amount and composition of seed proteins in rice. PloS One. 2012;7(9):e44493.
Sengupta S, Mangu V, Sanchez L, Bedre R, Joshi R, Rajasekaran K, Baisakh N. An actin-depolymerizing factor from the halophyte smooth cordgrass, Spartina alterniflora (SaADF2), is superior to its rice homolog (OsADF2) in conferring drought and salt tolerance when constitutively overexpressed in rice. Plant Biotechnol J. 2019;17(1):188–205.
Tang N, Ma SQ, Zong W, Yang N, Lv Y, Yan C, Guo ZL, Li J, Li X, Xiang Y, et al. MODD mediates deactivation and degradation of OsbZIP46 to negatively regulate ABA signaling and drought resistance in rice. Plant Cell. 2016;28(9):2161–77.
Zhang BC, Zhang LJ, Li F, Zhang DM, Liu XL, Wang H, Xu ZP, Chu CC, Zhou YH. Control of secondary cell wall patterning involves xylan deacetylation by a GDSL esterase. Nat Plants. 2017;3:17017.
Doll NM, Just J, Brunaud V, Caius J, Grimault A, Depege-Fargeix N, Esteban E, Pasha A, Provart NJ, Ingram GC, et al. Transcriptomics at maize embryo/endosperm interfaces identifies a transcriptionally distinct endosperm subdomain adjacent to the embryo scutellum. Plant Cell. 2020;32(4):833–52.
Dai DW, Ma ZY, Song RT. Maize endosperm development. J Integr Plant Biol. 2021;63(4):613–27.
Edmeades GO. Progress in achieving and delivering drought tolerance in maize-an update. ISAA: Ithaca; 2013. p. 1–39.
Feller U, Vaseva II. Extreme climatic events: Impacts of drought and high temperature on physiological processes in agronomically important plants. Front Environ Sci. 2014;2:39.
Sah RP, Chakraborty M, Prasad K, Pandit M, Tudu VK, Chakravarty MK, Narayan SC, Rana M, Moharana D. Impact of water deficit stress in maize: Phenology and yield components. Sci Rep-Uk. 2020;10(1):2944.
Feng F, Qi WW, Lv YD, Yan SM, Xu LM, Yang AY, Yuan Y, Chen YH, Zhao H, Song RT. OPAQUE11 is a central hub of the regulatory network for maize endosperm development and nutrient metabolism. Plant Cell. 2018;30(2):375–96.
Zhang ZY, Dong JQ, Ji C, Wu YR, Messing J. NAC-type transcription factors regulate accumulation of starch and protein in maize seeds. Proc Natl Acad Sci USA. 2019;116(23):11223–8.
Li CB, Yue YH, Chen HJ, Qi WW, Song RT. The ZmbZIP22 transcription factor regulates 27-kD gamma-Zein gene transcription during maize endosperm development. Plant Cell. 2018;30(10):2402–24.
Huang YC, Wang HH, Huang X, Wang Q, Wang JC, An D, Li JQ, Wang WQ, Wu YR. Maize VKS1 regulates mitosis and cytokinesis during early endosperm development. Plant Cell. 2019;31(6):1238–56.
Prathap V, Ali K, Singh A, Vishwakarma C, Krishnan V, Chinnusamy V, Tyagi A. Starch accumulation in rice grains subjected to drought during grain filling stage. Plant Physiol Bioch. 2019;142:440–51.
Lu HF, Hu YY, Wang CY, Liu WX, Ma G, Han QX, Ma DY. Effects of high temperature and drought stress on the expression of gene encoding enzymes and the activity of key enzymes involved in starch biosynthesis in wheat grains. Front Plant Sci. 2019;10:1414.
Li C, Li CY, Zhang RQ, Liang W, Kang XL, Jia Y, Liao YC. Effects of drought on the morphological and physicochemical characteristics of starch granules in different elite wheat varieties. J Cereal Sci. 2015;66:66–73.
Kimbembe RER, Li GY, Fu GF, Feng BH, Fu WM, Tao LX, Chen TT. Proteomic analysis of salicylic acid regulation of grain filling of two near-isogenic rice (Oryza sativa L.) varieties under soil drying condition. Plant Physiol Bioch. 2020;151:659–72.
Boehlein SK, Liu P, Webster A, Ribeiro C, Suzuki M, Wu S, Guan JC, Stewart JD, Tracy WF, Settles AM, et al. Effects of long-term exposure to elevated temperature on Zea mays endosperm development during grain fill. Plant J. 2019;99(1):23–40.
Leng PF, Zhao J. Transcription factors as molecular switches to regulate drought adaptation in maize. Theor Appl Genet. 2020;133(5):1455–65.
Budhagatapalli N, Narasimhan R, Rajaraman J, Viswanathan C, Nataraja KN. Ectopic expression of AtICE1 and OsICE1 transcription factor delays stress-induced senescence and improves tolerance to abiotic stresses in tobacco. J Plant Biochem Biot. 2016;25(3):285–93.
Zhang CQ, Xu Y, Lu Y, Yu HX, Gu MH, Liu QQ. The WRKY transcription factor OsWRKY78 regulates stem elongation and seed development in rice. Planta. 2011;234(3):541–54.
Xing Y, Jia WS, Zhang JH. AtMEK1 mediates stress-induced gene expression of CAT1 catalase by triggering H2O2 production in Arabidopsis. J Exp Bot. 2007;58(11):2969–81.
Polidoros AN, Mylona PV, Scandalios JG. Transgenic tobacco plants expressing the maize Cat2 gene have altered catalase levels that affect plant-pathogen interactions and resistance to oxidative stress. Transgenic Res. 2001;10(6):555–69.
Singh KK, Ghosh S. Regulation of glutamine synthetase isoforms in two differentially drought-tolerant rice (Oryza sativa L.) cultivars under water deficit conditions. Plant Cell Rep. 2013;32(2):183–93.
James D, Borphukan B, Fartyal D, Ram B, Singh J, Manna M, Sheri V, Panditi V, Yadav R, Achary VMM, et al. Concurrent overexpression of OsGS1;1 and OsGS2 genes in transgenic rice (Oryza sativa L.):impact on tolerance to abiotic stresses. Front Plant Sci. 2018;9:786.
Maskin L, Gudesblat GE, Moreno JE, Carrari FO, Frankel N, Sambade A, Rossi M, Iusem ND. Differential expression of the members of the Asr gene family in tomato (Lycopersicon esculentum). Plant Sci. 2001;161(4):739–46.
Park SI, Kim JJ, Shin SY, Kim YS, Yoon HS. ASR enhances environmental stress tolerance and improves grain yield by modulating stomatal closure in rice. Front Plant Sci. 2019;10:1752.
Eckardt NA. Oxylipin signaling in plant stress responses. Plant Cell. 2008;20(3):495–7.
Lu DL, Sun XL, Yan FB, Wang X, Xu RC, Lu WP. Effects of heat stress at different grain-filling phases on the grain yield and quality of waxy maize. Cereal Chem. 2014;91(2):189–94.
Zhu L, Wang DP, Sun JS, Mu YY, Pu WJ, Ma B, Ren FL, Yan WX, Zhang ZG, Li GY, Li YB, Pan YH. Phenotypic and proteomic characteristics of sorghum (Sorghum bicolor) albino lethal mutant sbe6-a1. Plant Physiol Bioch. 2019;139:400–10.
Wiśniewski JR, Zougman A, Nagaraj N, Mann M. Universal sample preparation method for proteome analysis. Nat Methods. 2009;6(5):359-U360.
Song YM, Zhong LJ, Zhou JT, Lu M, Xing TY, Ma LL, Shen J. Data-independent acquisition-based quantitative proteomic analysis reveals potential biomarkers of kidney cancer. Proteom Clin Appl. 2017;11:11–2.
Guo J, Chen L, Li YX, Shi YS, Song YC, Zhang DF, Li Y, Wang TY, Yang DG, Li CH. Meta-QTL analysis and identification of candidate genes related to root traits in maize. Euphytica. 2018;214(12):1–15.
Ma J, Chen T, Wu SF, Yang CY, Bai MZ, Shu KX, Li KL, Zhang GQ, Jin Z, He FC, Hermjakob H, Zhu YP. iProX: an integrated proteome resource. Nucleic Acids Res. 2019;47(D1):1211–7.
Research involving plants
The authors confirm that all the experimental methods and plants complied with relevant institutional, national, and international guidelines and legislation.
This work, including the design of the study and collection, analysis, and interpretation of data, and write manuscript, was supported by the National Natural Science Foundation of China (31771709, 32071958), the National Key Research and Development Program of China (2016YFD0300109, 2018YFD0200703), Jiangsu Agriculture Science and Technology Innovation Fund (CX3147), Jiangsu Agricultural Industry Technology System of China (JATS444), Priority Academic Program Development of Jiangsu Higher Education Institutions, and High-end Talent Support Program of Yangzhou University.
Ethics approval and consent to participate
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. Total proteins identified in waxy maize kernels.
: Table S2. Total differentially accumulated proteins identified at 10 dap and 25 dap.
: Table S3. Primer sequences of DAPs encoding genes used for qRT-PCR. Figure S1. Principal component analysis (PCA) of all proteins at 10 dap (A) and 25 dap (B) under WW and WS treatments. Figure S2. mRNA expression level analysis (qRT-PCR) of ten DAPs. Log2 (fold change) was calculated for WS conditions relative to WW. Figure S3. Classical maize genes among drought stress-responsive proteins. Log2 (fold change) of the proteins values was calculated for WS treatment relative to WW. Figure S4. KEGG pathway enrichment analysis of all DAPs at 10 dap and 25 dap . Figure S5. GO functional classification of all DAPs at 10 dap and 25 dap. Figure S6. Relative proteins level from proteomic involved in starch biosynthesis. The Log2 (fold change) was calculated for WW conditions relative to WS. Asterisks designate statistical (using Student’s t-test) significance P < 0.01.
About this article
Cite this article
Guo, J., Qu, L., Hu, Y. et al. Proteomics reveals the effects of drought stress on the kernel development and starch formation of waxy maize. BMC Plant Biol 21, 434 (2021). https://doi.org/10.1186/s12870-021-03214-z
- Waxy maize
- Drought stress
- Kernel development