Original Research ARTICLE
Plant Density and Nitrogen Supply Affect the Grain-Filling Parameters of Maize Kernels Located in Different Ear Positions
- 1College of Agriculture, Nanjing Agricultural University, Nanjing, China
- 2State Key Laboratory of Crop Biology, College of Agriculture, Shandong Agricultural University, Tai’an, China
Although yield output of maize (Zea mays L.) has improved markedly over the last century, procedures for improving the grain-filling process remain elusive. Our aim in this study was to relate grain-filling variation in maize (including kernels in apical and middle positions in the ears) to plant density and nitrogen (N) application rate using a crossed experimental design. We also investigated changes in zeatin riboside (ZR), indole-3-acetic acid (IAA), abscisic acid (ABA), and gibberellic acid (GA) in the kernels during the grain-filling period. Two high-yield maize varieties cultivated extensively in China were field grown under normal (67,500 pl ha-1) and high (97,500 pl ha-1) densities, and supplied with low, normal and high (0, 180, and 360 kg N ha-1) concentrations of N. Kernel weight (KW), the maximum grain-filling rate (Gmax), the average grain-filling rate (Gave), and the kernel weight increment achieving Gmax (Wmax) were all significantly depressed under high density (HD) conditions, but increased N supply partially offset the losses. The apical kernels were more sensitive to density and N application rate than middle kernels. Correlation analysis indicated that plant density and N rate affected KW mainly by influencing the grain-filling rate. Variation in ZR, IAA, and ABA content tracked the variation in KW, but variation in GA content did not. Furthermore, the grain-filling parameters (closely related to TKW) had strong canonical correlation with the content of all hormones across the filling period and ZR content had the strongest relationship. Based on our study, high N supply is beneficial to optimize grain-filling parameters and improve KW of maize kernels under crowded condition.
A gradual increase in plant density has been a vital contributor to maize yield enhancement worldwide (Tokatlidis and Koutroubas, 2004; Antonietta et al., 2014; Testa et al., 2016). However, kernel number per individual plant (KNP) and kernel weight (KW) were decreased due to the exacerbated competition for light and nutrients among individuals under high densities (Tollenaar et al., 2006; Borrás et al., 2007; Boomsma et al., 2009; Al-Naggar et al., 2015). Maize yield is determined by the kernel number per unit area (KNA) and the kernel weight (KW) (Borrás and Otegui, 2001; Mansfield and Mumm, 2014). Once the number of kernels is established, the final KW becomes a main factor determining maize yield (Wang et al., 2012). Hence, improving KW under crowded conditions remains a challenging issue for further enhancements in maize production.
The final weight of maize kernels is generally determined by genotype characteristics (Reddy and Daynard, 1983), which varies greatly with the position on the ear (cob) (Xu et al., 2013), and is also affected by the growth environment (Gambin et al., 2006; Zhou et al., 2017). Sub-optimal conditions related to population density, nutrient status, temperature, drought, and light intensity may restrict kernel set and development (Borrás et al., 2003; Brekke et al., 2011), as well limit the maximum weight achievable (Ouattar et al., 1987; Borras and Gambin, 2010; Shi et al., 2013). For KW, maize grains are usually classified as inferior grains or superior grains, where the former are generally located on the upper part of the ear, and the latter are mostly located on the middle and lower parts of the ear (Chen Y.J. et al., 2013). KW varies with the rate of kernel growth and the duration of grain-filling (Gambin et al., 2007; Sadras and Egli, 2008), which are directly related to grain weight and yield (Borras et al., 2009; Chen T. et al., 2013). During much of the grain-filling period, grain dry weight increases at an essentially linear rate over time. This period of linear increase is known as the effective filling-period duration (EFPD) (Johnson and Tanner, 1972). Borrás and Otegui (2001) reported that KW was closely related to the rate of kernel growth during the EFPD rather than the duration of this period. However, Bolaños (1995) argued that improvements in maize yield were largely explained by the prolongation of the grain-filling period.
The grain-filling process can be accurately modeled with logistic equations, thereby generating a series of important grain-filling parameters, such as maximum grain-filling rate (Gmax), average grain-filling rate (Gave) and active grain-filling period (Xu et al., 2013). Besides, the whole grain-filling process could be further divided to three periods according to the shape of the logistic curve, i.e., (i) gradual increase period (GIP), (ii) rapid increase period (RIP), and (iii) slight increase period (SIP) (Wang et al., 2014). Previous work showed the period of maximal grain-filling rate was protracted as fertilizer dosage increased, and the whole process of grain-filling was prolonged (Cao et al., 2008). Xu et al. (2016) indicated the Gmax and Gave varied with different sowing dates, but the duration of the active grain-filling period did not. Variations of light intensity also affected grain yield by changing the Gmax (Shi et al., 2013).
Exogenous growth conditions and endogenous plant hormones are both key regulators of kernel development (Lv et al., 2017). Plant hormones play important roles in the development and enrichment of maize, wheat and rice endosperms (Yang et al., 2000; Zhang et al., 2009; Liang et al., 2017). For instance, high zeatin (Z) and zeatin riboside (ZR) content in the initial stage of grain-filling can promote kernel setting and enhance maximum endosperm cell division rate (Morris et al., 1993; Yang et al., 2002). Auxins, gibberellic acid (GA), and abscisic acid (ABA) also regulate grain development (Kende and Zeevaart, 1997; Hansen and Grossmann, 2000; Yang et al., 2006; Zhu et al., 2011). Xu et al. (2013) showed that kernel filling rate was strongly positively correlated with the concentration of Z+ZR, IAA, and ABA in the kernels, but significantly negatively correlated with concentration of GA3. Besides, the relative level of hormone in both superior and inferior florets were shown to determine the development of inferior florets (Wang et al., 2006). All early reports have shown that the grain-filling of kernels is controlled by hormones.
Improvements in grain-filling remain a substantial challenge for modern crop production systems with high yield output (Kato et al., 2007). Previous studies regarding grain-filling, especially the relationships between different pollination time and position, plant growth rate at different growth periods, various resource-sink ratio and kernel number set and KW have been established under diverse stress conditions (Zinselmeier et al., 1999; Andrade et al., 2002; Uribelarrea et al., 2002; Gambin et al., 2006). Nevertheless, many studies were concentrated on the early filling stage or the filling rate over the entire filling period. Besides, studies on grain-filling characteristics and KW were mainly focused on variety comparisons, ecological conditions and modifications of single cultivation practices. However, these previous studies showed no evidence on the response of grain-filling characteristics to contrasting stand densities and N availabilities. Information on the relationship between kernel ear position and grain-filling characteristics is also limited. Thus, we examined the grain-filling parameters of maize kernels (including apical and middle position kernels in the ears) under different plant densities and N supply rates using logistic modeling. Besides, the hormone content in grains were determined to analyze their relationships with grain-filling characteristics.
Materials and Methods
Trial Site and Conditions
Field evaluations were performed at the Shandong Agricultural University Experimental Farm in Shandong, China during 2014 growing season. This area has a semi-humid, warm temperate continental climate with monsoons. The average monthly maximum temperature, minimum temperature, total solar radiation, and precipitation during the experimental period are shown in Table 1. The soil at the site is neutral sandy loam, and the nutrient status of the top 30 cm before seeding consisted of 11.1 g kg-1 organic matter, 0.76 g kg-1 available N, 25.4 mg kg-1 available phosphate and 105.1 mg kg-1 available potassium.
Table 1. Monthly average maximum temperature (Tmax, °C) and minimum temperature (Tmin, °C) and monthly total solar radiation (SR, MJ m-2) and rainfall (mm) during the growing season (June to September).
Two high-yield varieties planted widely in local production were selected as the testing materials, i.e., Denghai618 (521 × DH392) and Denghai605 (DH351 × DH382). The common density was 67,500 pl ha-1 and the N application rate ranged from 180 to 360 kg ha-1 in this location (Jin et al., 2012; Gao et al., 2017). Maize seeds were planted with hand planters on June 13 and harvested on September 29. Treatments consisted of two densities 67,500 pl ha-1 (low density, LD) and 97,500 pl ha-1 (high density, HD) and three N levels, 0 kg ha-1 (low N, N0), 180 kg ha-1 (normal N, N180), and 360 kg ha-1 (high N, N360). The experimental treatments were laid out in a split-split plot design, where the N rate was set as the main plot, and variety and density were set as subplot and sub-subplot with three replicates per treatment. The plot size was 12 m long × 3 m wide contained five rows and the space between rows was 60 cm. Plants within a row were placed 24.7 and 17.1 cm apart to achieve population density of 67,500 and 97,500 pl ha-1, respectively. Phosphorus (P2O5) and potassium (K2O) fertilizer were applied before sowing, at a rate of 90 and 120 kg ha-1, respectively. The N fertilizer was in the form of urea and applied as topdressing in two equal portions—one half at the jointing stage, and the other half at the large bell stage. Irrigation, weeds, diseases, and insect pests were controlled adequately during each growing season so that no factors other than density and N supply limiting growth.
In each plot, at least 50 representative plants with robust and uniform growth were labeled and bagged before silking. After the silks were fully exerted, artificial pollination was conducted uniformly on the same day. The shoot bags were replaced after pollination and removed only after pollen stopped shedding. All data were obtained from three rows of plants in the center of the cultivation plots.
Sampling and Measurements
Five labeled ears were sampled from each plot to determine grain-filling parameters at 10, 20, 30, 40, and 50 days after pollination (DAP). Kernels from the apical 5–25% files were collected and mixed to provide samples of apical kernels (inferior kernels). Kernels from the 40 to 60% files (counted from the apex to bottom) were collected and mixed to provide middle position grain (superior kernel) samples. We sampled 200 kernels from the five ears and randomly divided them into two equal portions. One portion was frozen in liquid N for 10 min and then stored at -80°C. The second portion was held in an oven for 30 min at 105°C and then dried to constant weight at 80°C.
Grain-filling process was fitted by the logistic equation W = A/ (1+Be-Ct). Where W refers to the kernel weight, t refers to days after pollination (the day of pollination set as t0 = 0), A refers to ultimate growth mass, B refers to primary parameter and C refers to growth rate parameter, respectively. The grain-filling parameters of the tested maize varieties were calculated as follows.
Occurrence time of maximal grain-filling rate (Tmax) = lnB/C;
Kernel weight increment achieving maximum grain-filling rate (Wmax) = A/2;
Maximum grain-filling rate (Gmax) = (C × Wmax) [1-(Wmax/A)];
Initial filling power (R0) = C;
Active grain-filling period (P) = 6/C.
The whole process of grain-filling was divided into three periods according to the shape of the logistic curve: (i) gradual increase period, (ii) rapid increase period, and (iii) slight increase period. The first and second order derivatives of the logistic equation were calculated to determine the filling parameters of these three grain-filling periods. The start time of peak grain-filling (t1) = (ln B-1.317)/C, the grain weight w1 = A/(1+Be-Ct1); the end time of peak grain-filling (t2) = (ln B+1.317)/C, the grain weight w2 = A/(1+Be-Ct2); the effective grain-filling period (when the grain weight reached 99% A, t3) = (ln B+4.59512)/C, the grain weight was w3. The grain-filling parameters of the three grain-filling period were calculated as the following equations.
The gradual increase period (T1) = t1-t0;
The rapid increase period (T2) = t2-t1;
The slight increase period (T3) = t3-t2;
The increment in grain weight during gradual increase period (W1) = w1-w0;
The increment in grain weight during rapid increase period (W2) = w2-w1;
The increment in grain weight during slight increase period (W3) = w3-w2;
The average grain-filling rate during gradual increase period (G1) = W1/T1;
The average grain-filling rate during rapid increase period (G2) = W2/T2;
The average grain-filling rate during slight increase period (G3) = W3/T3.
Besides, the average grain-filling rate during the whole grain-filling period was calculated as:
Average grain-filling rate (Gave) = W3/t3.
The endogenous ABA, IAA, GA, and ZR were extracted with 80% (v/v) methanol and quantified by enzyme-linked immunosorbent assay (ELISA) according to previous studies (Liang et al., 2017). The ELISA kits were provided by the Phytohormones Research Institute, China Agricultural University. The recovery rates for ABA, IAA, GA, and ZR were (90.2 ± 4.1), (86.1 ± 4.5), (79.4 ± 5.3), and (93.4 ± 6.5)%, respectively.
Grain Yield, Kernel Weight, and Kernel Number
At physiological maturity, 30 ears from the center three rows of each plot were harvested and air-dried to investigate yield, KNP and 1000-kernel weight (TKW). Grain yield (Mg ha-1) was expressed at 15.5% moisture.
The Curve Expert 1.3 software was used to fit the grain-filling process equation and obtain parameters A, B, and C. Analysis of variance (ANOVA) was performed with SPSS ver. 18.0 (SPSS Institute, Inc.). Duncan’s multiple range test was used to evaluate differences among treatments, and the significance level was set at the 0.05 probability level. Pearson correlation analysis and canonical correlation analysis were also performed. Figures in the article were plotted using SigmaPlot ver. 12.0 and GraphPad Prism ver. 7.
Grain Yield and Yield Components
Density, N rate, maize variety and the interactions between them (the density × variety interaction term was not included) had significant effects on GY (Table 2). The highest GY for variety DH618 was obtained under a combination of high density (HD) and high N rate (360 kg N ha-1, N360). The highest yield of variety DH605 was obtained with a combination of HD and a N rate of 180 kg ha-1 (N180). The GY response of DH618 to N (N180, N360, vs. N0) was greater under HD (GY increases of 38.1 and 46.3% with N180 and N360 fertilization, respectively) than at low density (LD) (GY increases of 31.5 and 32.4%, respectively). In DH605, the respective GY increases under HD were smaller: 22.2% (N180) and 22.5% (N360).
Table 2. Effects of density and N rate on grain yield (155 g kg-1 water content), kernel number and kernel weight.
Density and N rate had significant effects on KNP and 1000-kernel weight (TKW), but the density × N interaction term was not significant (P > 0.05). Increases in stand density were accompanied by decreases in KNP and TKW, but the values of both of these response parameters increased with increasing N rate, especially under HD (KNP values for DH618 increased on average by 27.9 and 17.9% under HD and LD, respectively). The KNP response of DH618 (average 23.0%) exceeded that of DH605 (average 9.9%). The difference in TKW between the two densities as the N rate increased was less marked than the KNP response. The TKW responses of DH618 exceeded those of DH605. For example, in DH618, the TKW value in the N360 × HD treatment increased by 7.9% in comparison with the N0 treatment, whereas the TKW proportional response to the same treatment was 4.8% in DH605.
Grain Weight Dynamic
Density, N rate and variety significantly affected KW in both the apical and middle ear positions (data not shown). The density × N rate interaction term was significant after 10 DAP for middle position kernels, but not for the apical grains. The KW of middle position kernels was higher than that of apical kernels over the whole grain-filling period. In each of the grain-filling periods, increased density reduced KW, while increased N application rates improved KW values (Figure 1), as expected. The responses of kernels in different ear positions were always different. Apical kernels responded more strongly to density and N supply than middle position kernels. Thus, in variety DH618, the KW values of apical and middle position kernels under HD were reduced on average by 8.6 and 5.9%, respectively, in comparison with the LD treatment. The respective average KW values in the high N treatment (N360) were 7.7 and 6.3% above the values in the N0 treatment. The response to N (N180, N360, vs. N0) was greater at HD. For example, under HD, the KW values of middle position kernels of DH618 at 30 DAP were 5.3% (N180) and 7.0% (N360) higher than those in the N0 treatment, whereas the respective increases under LD were 3.4% (N180) and 4.6% (N360).
Figure 1. Effects of density and nitrogen application rate on the grain dry weight and grain-filling rate of (A,B,E,F) upper kernels and (C,D,G,H) middle kernels (A–D) in Denghai618 and (E–H) Denghai605 during the 2014 growing season. LD and HD refer to low density and high density, respectively. N0, N180, and N360 represent nitrogen rates of 0, 180, and 360 kg ha-1, respectively.
Grain-Filling Characteristic Parameters
We divided the grain-filling parameter values under different D × N combinations by the value for the LDN0 treatment. These ratios are listed in Supplementary Table S1 and are presented in Figure 2. At 10 DAP, the different density and N rate combinations showed no significant effects. However, after 20 DAP, the differences in grain-filling rates among treatments became obvious. Progress in grain-filling was well-simulated by the Richards model: all correlation coefficients were > 0.99 (Table 3). The maximum grain-filling rate (Gmax), the kernel weight increment achieving Gmax (Wmax), average grain-filling rate (Gave) and the active grain-filling period (P) all decreased with increased density, but increased with the N supply rate. The initial grain-filling potential (R0) and the day on which the maximum grain-filling rate (Tmax) reached responded differently. Compared with the values for superior kernels, the Gmax, Gave, and P-values of inferior kernels were depressed. However, apical kernels responded more strongly to N supply. Gmax values in treatment N360 (DH618) were elevated by 8.6 and 5.6% in apical and middle position kernels, respectively, above the values in the N0 treatment. The grain-filling rate was higher in DH618 than in DH605, but the active grain-filling period of DH618 was shorter.
Figure 2. The variations of grain-filling characteristic parameters under different density and N application rate combinations. The results are the ratios of grain-filling characteristic parameters values under different D × N combinations to LDN0 for each variety (i.e., ratios of LDN180/LDN0, HDN180/LDN0, et al.), where LD and HD refer to low density and high density, respectively. N0, N180, and N360 represent nitrogen rates of 0, 180, and 360 kg ha-1, respectively. Gave, average grain-filling rate; Tmax, the day reaching the maximum grain-filling rate; Wmax, kernel weight increment achieving maximum grain-filling rate; Gmax, maximum grain-filling rate; R0, initial grain-filling potential; P, active grain-filling period; T1, grain-filling duration of gradual increase period; W1, increased grain weight of gradual increase period; G1, mean grain-filling rate of gradual increase period; T2, grain-filling duration of rapid increase period; W2, increased grain weight of rapid increase period; G2, mean grain-filling rate of rapid increase period; T3, grain-filling duration of slight increase period; W3, increased grain weight of slight increase period; G3, mean grain-filling rate of slight increase period.
Table 3. Effects of density and N rate on grain-filling characteristics parameters of maize kernels at upper and middle position.
Density increases prolonged the duration of the GIP in most cases, but shortened the durations of two other periods on the logistic curve. The GIP, RIP and SIP durations were all prolonged in both inferior and superior kernels of the two varieties by increasing the N supply rate (Table 4). The increases in KW of inferior kernels through all the three filling periods were smaller compared with the increases in superior kernels. The duration of the apical kernel RIP and SIP were shorter, but the duration of the GIP was longer in comparison with middle position kernels. In addition, the duration of the DH618 GIP was shorter than the DH605 GIP, but the grain-filling rate and the increases in grain weight during the GIP were higher in DH618.
Table 4. Effects of density and N rate on characteristics parameters of the three grain-filling phases of maize kernels at upper and middle position.
Variations in Hormone Content
During the grain-filling process, the overall trends in IAA, ZR, and ABA content in the grains showed unimodal curves. The GA content declined continually in all treatments (Figure 3). Plant density and N rate both had significant effects on hormonal content during the whole filling period. The density × N rate interaction was significant for most hormone content (Supplementary Table S2). Density and N rate did not alter the general trends of any of the hormones over time. At each grain-filling stage, IAA and ZR content were higher in middle position kernels than in the apical kernels, but the GA content was higher in the apical kernels. Before 30 DAP, the concentration of ABA was much lower in the apical kernels than in the middle position grains. The hormone content in apical kernels responded more strongly to density and N supply than in middle position kernels. For example, the Z+ZR content at 20 DAP was reduced by average 8.8 and 3.7% (HD vs. LD), but were increased by average 17.3 and 5.9% (N360 vs. N0) for inferior and superior kernels, respectively. The response to N (N180, N360, vs. N0) was also elevated in the HD treatment. Thus, the Z+ZR content at 20 DAP (middle position kernels of DH618) was increased by 2.8% (N180) and 3.5% (N360) under LD, and by 6.2% (N180) and 7.7% (N360) under HD.
Figure 3. Effects of density and nitrogen application rate on the content of zeatin riboside (ZR), indole-3-acetic acid (IAA), abscisic acid (ABA), and gibberellic acid (GA) in (A,C,E,G,I,K,M,O) upper kernels and (B,D,F,H,J,L,N,P) middle kernels of the two varieties during 2014 growing season. LD and HD refer to low density and high density, respectively. N0, N180, and N360 represent nitrogen rates of 0, 180, and 360 kg ha-1, respectively.
Correlation Between 1000-Grain Weight and Grain-Filling Parameters
We examined the relationships between grain-filling parameters in different ear positions and the TKW (Table 5). Correlation analysis showed that Gave, Wmax, Gmax, the grain-filling rate and KW increased during the GIP, RIP, and SIP were significantly positively correlated with TKW at physiological maturity in kernels at both ear positions. We also found a significant negative relationship between KW and the duration of the GIP in middle position kernels.
Table 5. Pearson correlation of 1000-kernel weight (TKW) and filling parameters of maize kernels at different ear positions.
Correlation Between Grain-Filling Parameters and Hormone Content
Canonical correlation analysis was used to examine the relationships between hormone content (X) and grain-filling parameters (Y). Hormone content at 10, 20, 30, 40, and 50 DAP were set as variables x1, x2, x3, x4, and x5, respectively. Gave, Wmax, Gmax, G1, G2, and G3 were set as y1, y2, y3, y4, y5, and y6, respectively. As shown in Table 6, the first canonical correlation coefficients were all statistically significant, indicating that the correlations between the grain-filling parameters and hormone content can be described by the first pair of typical variables (except the ABA content of superior kernels). The linear combinations of the first pair of canonical correlation models are shown in Table 7. From the coefficients for each variable, we were able to show that the grain-filling parameters of apical kernels were mainly determined by the content of ZR at 40 DAP, IAA at 20 DAP, ABA at 10 DAP and GA at 40 DAP, whereas content of ZR at 50 DAP, IAA and ABA at 20 DAP, and GA at 10 DAP were the main factors determining the grain-filling parameters of the middle position kernels. Most information of X and Y were concentrated in the first pair of canonical variables (L1 and M1, Table 8). Among the hormones examined in our study, ZR had the strongest correlation with grain-filling parameters of both inferior and superior kernels (86.0 and 79.6% of the ZR variation may affect variation in grain-filling parameters in inferior and superior kernels, respectively; 90.2 and 91.9% of the variation in grain-filling parameters were related to ZR content changes in the inferior and superior kernels, respectively).
Table 6. Canonical correlation coefficient between grain filling parameters and hormone content and the significance testing.
Table 8. The proportion of variable variation can be explained by canonical variable (redundancy analysis).
Although maize GY is determined mainly by KNA, variations in KW are also influential (Zhou et al., 2016). Grain-filling is a critical physiological process that determines the KW and yield of corn kernels. Poor grain-filling limits the potential GY (Yang and Zhang, 2010) and is among the major scientific and practical issues that should be addressed to improve agricultural production. Appropriate plant densities and N fertilization rates are important elements of efficient management that have contributed to maize yield improvement (Wei et al., 2017). We showed in this study that density, N rate and their interaction affected KW during the whole filling period, while the impact of density tended to be stronger. The variation in KW between different floret positions of the maize ear was substantial and mostly related to changes in the grain-filling rate. Furthermore, the low KW values in apical kernels were exacerbated under HD, though the yield potential under HD was nevertheless 9.8% higher than that under LD. Taken together with the strong response of apical kernels to N rate (especially under HD), we concluded that apical grains were more sensitive to environmental factors than middle grains. Hence, promoting grain-filling of grains at apical position will be key for improving the KW of maize.
The KW is determined by the grain-filling rate, the effective filling duration, or both (Zhou et al., 2017). In our study, Gave, Gmax and the filling rate during GIP, RIP, and SIP were all positively correlated with KW, but we found no obvious correlation between KW and the active grain-filling period in neither superior kernels nor inferior kernels (Table 6). Thus, changes in the filling rate were mainly a result of variation in KW under different density and N rate combinations. Li et al. (2000) showed that difference in KW within the same maize variety was determined by the grain-filling rate, whereas KW variation among varieties was determined by grain-filling duration. We attributed KW differences between the two maize varieties that we tested to differences in filling rate, rather than filling duration. Unlike apical kernels, the duration of the GIP in middle position kernels was strongly negatively correlated to KW. A mechanism for shortening the duration of the GIP whilst maintaining the filling rate could prolong the duration of the subsequent RIP relative to the SIP. This procedure could increase KW.
The processes of cell division, differentiation, enlargement, and enrichment of inclusions during the development of maize kernels are controlled by a range of hormones operating during different filling periods (Yang et al., 2001; Artlip et al., 2010). In the current study, density and N level did not change the regulation period of the hormones in grains, however, hormone content varied significantly through kernel development. HD decreased the highest values, and accelerated reductions in ZR and IAA content, whereas increased N application offset these declines, especially under HD. The strong canonical correlation between grain-filling parameters and ZR content indicated that density and N rate combinations affected grain-filling progress largely by adjusting the content of ZR, which regulated the rates of endosperm cell division and grain filling (Ahmad et al., 2018). Additionally, we also showed that the ZR and IAA content of apical kernels were significantly lower than the content in middle position kernels. Hence, we attributed the lower KW in the apical kernels to their lower ZR and IAA content, these low hormone levels constrained the proliferation of endosperm cells (Wang et al., 2006), reduced the sink capacity, and ultimately decreased KW. ZR and IAA are considered to be the main regulators of endosperm cell development, whereas ABA likely has a vital regulatory role in the process of photosynthate transport and accumulation in the grain (Hansen and Grossmann, 2000; Travaglia et al., 2007; Qin et al., 2013). We found that the ABA concentration in the kernels was reduced by plant crowding, but enhanced by increased N availability, in agreement with the findings of Gawronska et al. (2003). Furthermore, we determined that about 75% of the variation in grain-filling parameters was related to variation in ABA content in both apical and middle position kernels. The ABA content during the early filling period (10–20 DAP) was most strongly positively correlated with the grain-filling parameters, as shown in previous studies on maize, wheat, and rice (Xu et al., 2007, 2013; Zhang et al., 2009; Lv et al., 2017). The combination of elevated ABA concentration after 30 DAP and the low KW in the HD × N0 treatment also suggested that maintaining a high level of ABA during the late grain-filling stage may negatively impact the progress of grain-filling.
Unlike the other hormones, the GA levels in both apical and middle position kernels were highest at 10 DAP, when the embryo was enlarging rapidly. In constract to the ZR content, GA content increased when density was increased, but declined with increased N application. The concentration of GA in the apical kernels exceeded that in the middle position kernels. Earlier studies have indicated that GA in the kernels played an important role in the rapid amplification of endosperm cells. However, we found that the GA content in the HD × N0 was maintained at a high level over the whole filling period, suggesting that elevated GA content adversely affected the improvement of KW. This outcome is congruent with earlier work (Xu et al., 2013), showing that GA can enhance the activities of α-amylase and other hydrolytic enzyme that promote the decomposition of starch (Yang et al., 2001). Thus, changes in ZR, IAA, ABA, and GA content in grains were all significantly correlated with the grain-filling parameters. Further work is needed to explore the mechanisms controlling (i) hormone content variation under different density × N rate combinations, and (ii) hormone distribution within maize ears over the whole grain-filling period.
In short, variation in grain-filling rate was the main factor affecting TKW under different density × N rate combinations. The apical kernels were more sensitive to density and N application rate than middle position kernels. High N input weakened the reduction in KW, Gmax, Gave, Wmax, RIP, and SIP durations under high density. The grain-filling parameters (closely related to TKW) had strong canonical correlations with the content of all hormones across the filling period, and ZR content had the strongest relationship. Therefore, we conclude from our results that high N supply is beneficial to optimize grain-filling parameters and improve KW of maize kernels under crowded condition.
SW designed the work, carried out all experiments and data analysis and drafted the manuscript. XW performed field experiments and sampling and helped drafting the work. SD conceived the study and planned the experiments. GL and YQ helped draft the manuscript. DJ helped draft and revise the manuscript. All authors read and approved the final manuscript.
This research was funded by the Natural Science Foundation of Jiangsu Province (Grant No. BK20170720), the China Postdoctoral Science Foundation funded project (Grant No. 2017M611832), National Natural Science Foundation of China (Grant No. 3117 1497), the National Basic Research Program of China (973 Program, Grant No. 2011CB100105), the National Food Science and Technology of High-yield Program of China (Grant No. 2011BAD16B09), the Special Fund for Agro-scientific Research in the Public Interest of China (Grant No. 20120306), and supported by the 111 Project (Grant No. B16026).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors gratefully acknowledge Liyuan Hou, Meifeng Zeng, and Ting Liang for helping us to perform the study. They would also like to thank Mr. Mingcai Xu and all the workers in their farm for managing the field experiments.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2019.00180/full#supplementary-material
TABLE S1 | The data used in Figure 2.
TABLE S2 | Analysis of variance of hormone contents in maize grains under different density and N rate conditions.
G1, mean grain-filling rate of gradual increase period; G2, mean grain-filling rate of rapid increase period; G3, mean grain-filling rate of slight increase period; Gave, average grain-filling rate; Gmax, maximum grain-filling rate; KW, kernel weight; P, active grain-filling period; R0, initial grain-filling potential; T1, grain-filling duration of gradual increase period; T2, grain-filling duration of rapid increase period; T3, grain-filling duration of slight increase period; TKW, 1000-kernel weight; Tmax, the day reaching the maximum grain-filling rate; W1, increased grain weight of gradual increase period; W2, increased grain weight of rapid increase period; W3, increased grain weight of slight increase period; Wmax, kernel weight increment achieving maximum grain-filling rate.
Ahmad, I., Kamran, M., Ali, S., Cai, T., Bilegjargal, B., Liu, T., et al. (2018). Seed filling in maize and hormones crosstalk regulated by exogenous application of uniconazole in semiarid regions. Environ. Sci. Pollut. Res. Int. 25, 33225–33239. doi: 10.1007/s11356-018-3235-0
Al-Naggar, A. M. M., Shabana, R. A., Atta, M. M. M., and Al-Khalil, T. H. (2015). Maize response to elevated plant density combined with lowered N-fertilizer rate is genotype-dependent. Crop J. 3, 96–109. doi: 10.1016/j.cj.2015.01.002
Andrade, F. H., Echarte, L., Rizzalli, R., Della Maggiora, A., and Casanovas, M. (2002). Kernel number prediction in maize under nitrogen or water stress. Cropence 42, 1173–1179. doi: 10.2135/cropsci2002.1173
Antonietta, M., Fanello, D. D., Acciaresi, H. A., and Guiamet, J. J. (2014). Senescence and yield responses to plant density in stay green and earlier-senescing maize hybrids from Argentina. Field Crops Res. 155, 111–119. doi: 10.1016/j.fcr.2013.09.016
Artlip, T. S., Madison, J. T., and Setter, T. L. (2010). Water deficit in developing endosperm of maize - cell division and nuclear DNA endoreduplication. Plant Cell Environ. 18, 1034–1040. doi: 10.1111/j.1365-3040.1995.tb00614.x
Boomsma, C. R., Santini, J. B., Tollenaar, M., and Vyn, T. J. (2009). Maize morphophysiological responses to intense crowding and low nitrogen availability: an analysis and review. Agron. J. 101, 1426–1452. doi: 10.2134/agronj2009.0082
Borras, L., and Gambin, B. L. (2010). Trait dissection of maize kernel weight: towards integrating hierarchical scales using a plant growth approach. Field Crops Res. 118, 1–12. doi: 10.1016/j.fcr.2010.04.010
Borras, L., Zinselmeier, C., Lynn, M., Westgate, M. E., and Muszynski, M. G. (2009). Characterization of grain-filling patterns in diverse maize germplasm. Crop Sci. 49, 999–1009. doi: 10.2135/cropsci2008.08.0475
Brekke, B., Edwards, J., and Knapp, A. (2011). Selection and adaptation to high plant density in the iowa stiff stalk synthetic maize (Zea mays L.) population. Crop Sci. 51, 1965–1972. doi: 10.2135/cropsci2010.09.0563
Chen, T., Xu, Y., Wang, J., Wang, Z., Yang, J., and Zhang, J. (2013). Polyamines and ethylene interact in rice grains in response to soil drying during grain filling. J. Exp. Bot. 64, 2523–2538. doi: 10.1093/jxb/ert115
Gao, J., Zhao, B., Dong, S., Liu, P., Ren, B., and Zhang, J. (2017). Response of summer maize photosynthate accumulation and distribution to shading stress assessed by using 13CO2 stable isotope tracer in the field. Front Plant Sci. 8:1821. doi: 10.3389/fpls.2017.01821
Gawronska, H., Deji, A., Sakakibara, H., and Sugiyama, T. (2003). Hormone-mediated nitrogen signaling in plants: implication of participation of abscissic acid in negative regulation of cytokinin-inducible expression of maize response regulator. Plant Physiol. Biochem. 41, 605–610. doi: 10.1016/S0981-9428(03)00059-7
Jin, L., Cui, H., Li, B., Zhang, J., Dong, S., and Liu, P. (2012). Effects of integrated agronomic management practices on yield and nitrogen efficiency of summer maize in North China. Field Crops Res. 134, 30–35. doi: 10.1016/j.fcr.2012.04.008
Johnson, D. R., and Tanner, J. W. (1972). Comparisons of Corn (Zea mays L.) inbreds and hybrids grown at equal leaf area index, light penetration, and population 1. Crop Sci. 12, 482–485. doi: 10.2135/cropsci1972.0011183X001200040027x
Kato, T., Dai, S., and Taniguchi, A. (2007). Activities of enzymes for sucrose-starch conversion in developing endosperm of rice and their association with grain filling in extra-heavy panicle types. Plant Prod Sci. 10, 442–450. doi: 10.1626/pps.10.442
Liang, W., Zhang, Z., Wen, X., Liao, Y., and Liu, Y. (2017). Effect of non-structural carbohydrate accumulation in the stem pre-anthesis on grain filling of wheat inferior grain. Field Crops Res. 211, 66–76. doi: 10.1016/j.fcr.2017.06.016
Lv, X., Han, J., Liao, Y., and Liu, Y. (2017). Effect of phosphorus and potassium foliage application post-anthesis on grain filling and hormonal changes of wheat. Field Crops Res. 214, 83–93. doi: 10.1016/j.fcr.2017.09.001
Morris, R. O., Blevins, D. G., Dietrich, J. T., Durley, R. C., Gelvin, S. B., Gray, J., et al. (1993). Cytokinins in plant-pathogenic bacteria and developing cereal-grains. Aust. J. Plant Physiol. 20, 621–637. doi: 10.1071/PP9930621
Ouattar, S., Jones, R. J., and Crookston, R. K. (1987). Effect of water deficit during grain filling on the pattern of maize kernel growth and development. Crop Sci. 27, 726–730. doi: 10.1111/gcb.12381
Qin, S., Zhang, Z., Ning, T., Ren, S., Su, L., and Li, Z. (2013). Abscisic acid and aldehyde oxidase activity in maize ear leaf and grain relative topost-flowering photosynthetic capacity and grain-filling rate under different water/nitrogen treatments. Plant Physiol. Biochem. 70, 69–80. doi: 10.1016/j.plaphy.2013.04.024
Shi, J., Cui, H., Zhao, B., Dong, S., Liu, P., and Zhang, J. (2013). Effect of light on yield and characteristics of grain-filling of summer maize from flowering to maturity. Sci. Agric. Sin. 46, 4427–4434.
Testa, G., Reyneri, A., and Blandino, M. (2016). Maize grain yield enhancement through high plant density cultivation with different inter-row and intra-row spacings. Eur. J. Agron. 72, 28–37. doi: 10.1016/j.eja.2015.09.006
Tokatlidis, I. S., and Koutroubas, S. D. (2004). A review of maize hybrids’ dependence on high plant populations and its implications for crop yield stability. Field Crops Res. 88, 103–114. doi: 10.1016/j.fcr.2003.11.013
Travaglia, C., Cohen, A. C., Reinoso, H., Castillo, C., and Bottini, R. (2007). Exogenous abscisic acid increases carbohydrate accumulation and redistribution to the grains in wheat grown under field conditions of soil water restriction. J. Plant Growth Regul. 26, 285–289. doi: 10.1007/s00344-007-9018-3
Wei, S., Wang, X., Zhu, Q., Jiang, D., and Dong, S. (2017). Optimising yield and resource utilisation of summer maize under the conditions of increasing density and reducing nitrogen fertilization. Sci. Nat. 104, 86. doi: 10.1007/s00114-017-1509-x
Xu, G., Zhang, J., Lam, H. M., Wang, Z., and Yang, J. (2007). Hormonal changes are related to the poor grain filling in the inferior spikelets of rice cultivated under non-flooded and mulched condition. Field Crops Res. 101, 53–61. doi: 10.1016/j.fcr.2006.09.008
Xu, Y., Gu, D., Zhang, B., Zhang, H., Wang, Z., and Yang, J. (2013). Hormone contents in kernels at different positions on an ear and their relationship with endosperm development and kernel filling in maize. Acta Agron. Sin. 39, 1452–1461. doi: 10.3724/SP.J.1006.2013.01452
Yang, J., Peng, S., Visperas, R. M., Sanico, A. L., Zhu, Q., and Gu, S. (2000). Grain filling pattern and cytokinin content in the grains and roots of rice plants. Plant Growth Regul. 30, 261–270. doi: 10.1023/A:1006356125418
Yang, J., Zhang, J., Huang, Z., Wang, Z., Zhu, Q., and Liu, L. (2002). Correlation of cytokinin levels in the endosperms and roots with cell number and cell division activity during endosperm development in rice. Ann. Bot. 90, 369–377. doi: 10.1093/aob/mcf198
Yang, J., Zhang, J., Liu, K., Wang, Z., and Liu, L. (2006). Abscisic acid and ethylene interact in wheat grains in response to soil drying during grain filling. New Phytol. 171, 293–303. doi: 10.1111/j.1469-8137.2006.01753.x
Yang, J., Zhang, J., Wang, Z., Zhu, Q., and Wang, W. (2001). Hormonal changes in the grains of rice subjected to water stress during grain filling. Plant Physiol. 127, 315–323. doi: 10.1104/pp.127.1.315
Zhang, H., Tan, G. L., Yang, L. N., Yang, J. C., Zhang, J. H., and Zhao, B. H. (2009). Hormones in the grains and roots in relation to post-anthesis development of inferior and superior spikelets in japonica/indica hybrid rice. Plant Physiol. Biochem. 47, 195–204. doi: 10.1016/j.plaphy.2008.11.012
Zhou, B., Yue, Y., Sun, X., Ding, Z., Ma, W., and Zhao, M. (2017). Maize kernel weight responses to sowing date-associated variation in weather conditions. Crop J. 5, 43–51. doi: 10.1016/j.cj.2016.07.002
Zhou, B. Y., Yue, Y., Sun, X. F., Wang, X. B., Wang, Z. M., Ma, W., et al. (2016). Maize grain yield and dry matter production responses to variations in weather conditions. Agron. J. 108, 196–204. doi: 10.2134/agronj2015.0196
Zhu, G., Ye, N., Yang, J., Peng, X., and Zhang, J. (2011). Regulation of expression of starch synthesis genes by ethylene and ABA in relation to the development of rice inferior and superior spikelets. J. Exp. Bot. 62, 3907–3916. doi: 10.1093/jxb/err088
Keywords: canonical correlation analysis, grain-filling characteristic, kernel position, hormone, logistic curve, maize
Citation: Wei S, Wang X, Li G, Qin Y, Jiang D and Dong S (2019) Plant Density and Nitrogen Supply Affect the Grain-Filling Parameters of Maize Kernels Located in Different Ear Positions. Front. Plant Sci. 10:180. doi: 10.3389/fpls.2019.00180
Received: 18 June 2018; Accepted: 05 February 2019;
Published: 01 March 2019.
Edited by:Abraham J. Escobar-Gutiérrez, Institut National de la Recherche Agronomique (INRA), France
Reviewed by:Axel Tiessen, Centro de Investigación y de Estudios Avanzados (CINVESTAV), Mexico
Michael Gerard Muszynski, University of Hawai’i at Mānoa, United States
Copyright © 2019 Wei, Wang, Li, Qin, Jiang and Dong. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.