Nitrogen Supply Regulates Vascular Bundle Structure and Matter Transport Characteristics of Spring Maize Under High Plant Density

Nitrogen (N) fertilizer application greatly enhances grain yield by improving dry matter accumulation and grain filling in spring maize. However, how N application rates regulate the vascular bundle structure, matter transport and grain filling of spring maize under a high planting density has been poorly understood thus far. In this study, we analyzed the relationship between grain filling, vascular bundle structure and matter transport efficiency (MTE) of spring maize in the field. Zhongdan909 (ZD909) was used as the experimental material in a 2-year field experiment from 2015 to 2016, and it was grown under different N levels (0, 150, and 300 kg N ha^–1) applied to the grain-filling stage of plots with planting densities of 67,500 plants ha^–1 (ND) and 90,000 plants ha^–1 (HD). Nitrogen application significantly optimized the structure of the big and small vascular bundles. In particular, there was an increase in the total number of small vascular bundles in the peduncle and cob of the ear system, i.e., increases of 51.8% and 25.7%, respectively, and the proportions of small vascular bundles to the total number of vascular bundles in the peduncle and cob were significantly increased. The root bleeding sap and MTE of maize were significantly increased by N application under both ND and HD, as indicated by the significant increase in the rate of ¹³C-photosynthate allocation to grain and amount of postsilking dry matter at maturity. Moreover, N application greatly improved the mean grain-filling rate (G_mean) under ND and HD by 30.0% and 36.1%, respectively, and the grain-filling rate increased, leading to a distinct improvement in the grain sink at the grain-filling stage. We concluded that nitrogen application significantly optimized the vascular bundle structure of the ear system, increased the MTE and improved photosynthate distribution to the grain, ultimately enhancing the filling rate and grain yield.

Introduction

Nitrogen (N) is the key nutrient element required for plant growth, and maize (Zea mays L.) production highly depends on the appropriate application of nitrogen fertilizer (Mu et al., 2016; Di Salvo et al., 2018). However, N requirement of maize is changing with various plant density. For instance, it is necessary to properly increase N input for enhancing dry matter accumulation and distribution and decreasing the competition of nutrients among individuals to increase the grain yield under high plant density, compared with normal plant density (Rajcan and Tollenaar, 1999; Wang et al., 2006; Ren et al., 2017). An insufficient N input will result in the issue of nutrient supply for plants, severely reduce the plant photosynthesis, accumulation of matter and grain yield (Tilman et al., 2002; Chen et al., 2011, 2015; Wei et al., 2017). Hence, the reasonable selection of N fertilizer input is beneficial to grain filling and grain weight in different periods after silking through coordinating crop photosynthesis, dry matter accumulation and distribution and others processes (Austin et al., 1997; Borrás et al., 2004). It has been recognized that the duration and rate at the reproductive growth stage determine the final grain weight and maize yield (Borrás et al., 2004; Zhou et al., 2017).

The grain-filling process can be divided into two stages: the first is the lag stage in which the kernel number and the potential kernel weight (size) are determined, and the second stage is the effective grain-filling stage in which the final kernel weight is established (Jones et al., 1996; Gambín et al., 2006; Hisse et al., 2019). The final kernel weight is highly dependent on the accumulation of assimilates during the initial grain-filling stages, but under stressful environmental conditions (such as low nitrogen and high planting density), the limitations of the supply of assimilates seriously influence the initial grain-filling stages and reduce the final kernel weight (Ouattar et al., 1987; Cirilo and Andrade, 1996; Gambín et al., 2006). Previous studies have focused on increasing the grain weight through nitrogen fertilizers by extending the effective filling period (Cao et al., 2008; Shi et al., 2013). However, Xu et al. (2016) pointed out that appropriate increasing nitrogen fertilizer can change the grain-filling rate and can increase the grain weight.

Actually, the N supply and planting density influence the potential assimilate supply and assimilate accumulation of maize grain, and matter transport efficiency (MTE) regulates the process of grain filling (Austin et al., 1997; Tollenaar et al., 2006; Wang et al., 2006). The assimilation of carbohydrates for grain-filled cereals are generally derived from two resources: (i) the current assimilates transferred directly to the grain and (ii) the assimilates redistributed from reserve pools in vegetative tissues either preanthesis or postanthesis (Gan and Amasino, 1997; Gebbing and Schnyder, 1999). However, under conditions associated with low nitrogen or high density, grain filling is limited by matter accumulation and increased transport (D’Andrea et al., 2016). As a result, the final grain weight decreases dramatically due to growth restrictions during the active grain-filling stage, and the accumulation of dry matter and the matter translocation rate have a decisive effect on grain filling (Ntanos and Koutroubas, 2002). Increasing studies have indicated that the main factor affecting maize grain filling under high-density conditions is not matter production but rather matter transport characteristics (Yang et al., 2002; Wang et al., 2019; Yan et al., 2019).

The smooth flow of matter transport depends on the structure and function of the vascular bundle systems above- and belowground (Zhong et al., 1997; Noguchi et al., 2005). It is generally accepted that xylem transports water and phloem transports organic matter (Shane et al., 2000). Organic matter, such as leaf photosynthetic products, is transported to the grain through the phloem, thereby achieving the purpose of increasing the grain weight, and the phloem unloads 60–98.2% of total carbon into the grain (Khan et al., 1998; Wigoda et al., 2014). Higher numbers and areas of vascular bundles can regulate the distribution of dry matter after silking by regulating the metabolism of carbon (C) and N and, to a certain extent, can determine the accumulation of C and N in the grain and the size of the grain sink, which affects grain filling under various conditions (Housley and Peterson, 1982; Yang et al., 2001; Monneveux et al., 2005). A previous study showed that tillage practice and row spacing could optimize the structure of the vascular bundles, improve matter transport and thus increase maize yield (Piao et al., 2016). However, the mechanism of how N application under high-density conditions regulates the grain-filling efficiency of maize by affecting vascular bundle characteristics is unclear.

Therefore, the purpose of this study is to clarify three points: (1) how the N and plant density interaction influences grain filling and matter transport characteristics, (2) how the N and plant density interaction affects the vascular bundle structure, and (3) the relationship between grain filling and vascular bundle structure under various N fertilizer conditions.

Materials and Methods

Site Description

The field experiment was conducted in Gongzhuling, Jilin Province, China (43°31′N, 124°48′E) in 2015 and 2016. The experimental soil (brown loam) consisted of organic matter 25.9 g kg^–1, total N 1.6 g kg^–1, available N 140.2 mg kg^–1, available P 62.3 mg kg^–1, and available K 148.40 mg kg^–1. Precipitation (mm) and mean air temperatures measured during the growth stage of spring maize are shown in Supplementary Figure 1. The effective cumulative temperature was calculated using the mean daily air temperature above 10°C; the respective rainfall and effective cumulative temperature during the growth stage of spring maize were 409.6 mm and 1631.3°C in 2015 and 643.7 mm and 1616.0°C in 2016.

Experimental Design and Field Management

The field experiment with a split plot design was performed with two planting densities (67,500 plants ha^–1, ND; 90,000 plants ha^–1, HD) as the main plot and three N levels, 0 kg N ha^–1 (N0), 150 kg N ha^–1 (N150), and 300 N kg ha^–1 (N300), as subplots. All treatments were set with three replicates. The size of each plot was 6 m × 7.5 m with 0.6 m row spacing. The application of N fertilizer (urea) was performed before sowing at the elongation and silking stages of maize at a 5:3:2 ratio. Both Ca₃(PO₄)₂ and KCl were applied before sowing at 100 kg ha^–1 in 2015 and 2016. Pests, weeds and diseases were well controlled, and no irrigation was applied throughout the growing season. Maize seed was planted with hand planters on 30 April 2015 and 30 April 2016, and the respective harvest dates were October 1, 2015 and 2016.

Data Collection

Characterization of Vascular Bundles and Matter Transport Efficiency

The characterization of the vascular bundles of maize was performed according to the description provided by He et al. (2005). Briefly, at the milking stage, five random plants in each plot were manually sampled and dissected to obtain the basal stem (1 cm above the ground), ear internodes, peduncle and cob fractions. All samples collected were fixed using Kano fixative solution containing acetic acid/alcohol = 1:3 (v/v) and then stored in 70% ethanol until imaging. Freehand sectioning, safflower staining, and observation of the vascular bundle structure were performed with a Zeiss Axio A1 upright microscope with a 5×/0.3 numerical aperture (NA) and 10×/0.3 NA Axio HRc camera (Carl Zeiss Inc., Ontario, CA, United States). The numbers of small vascular bundles (average of 18 adjacent vascular bundles) were recorded, and the built-in ZEN analysis system was used to calculate the area occupied by the vascular bundles in both the xylem and phloem.

To determine the MTE, root-bleeding sap was sampled from the above maize plants by using lopping shears to cut 1 cm above the basal stem internode of five representatives per plot. Twelve grams of absorbent cotton was wrapped around the basal stem and secured with rubber bands to collect the root-bleeding sap. The weight of the root-bleeding sap was measured after 12 h of collection (from 17:00 to 05:00). Then, the MTE (mg mm^–2h^–1) was calculated as follows, Eq. 1 (Piao et al., 2016).

$$MTE=RBS/VBA(1)$$

where RBS refers to the root-bleeding sap collected from 17:00 to 05:00 of the next day (mg h^–1) and VBA refers to the vascular bundle area of the basal stem internode (mm²).

Grain-Filling Parameters

During the grain-filling stage, three plants per plot were randomly sampled. In 2015, the N0 treatment was sampled at 10, 25, 32, 43, and 57 days after pollination for the ND and HD plots; the N150 and N300 treatments were sampled at 10, 20, 30, 38, and 53 days after pollination for the ND plot and at 10, 20, 30, 40, and 55 days after pollination for the HD plot. In 2016, the N0 treatment was sampled at 7, 14, 24, 33, 48, and 56 days after pollination for the ND plot and at 7, 14, 24, 31, 46 and 56 days after pollination for the HD plot; and the N150 and N300 treatments were sampled at 7, 14, 25, 32, 47, and 56 days after pollination for the ND and HD plots. We sampled 100 kernels in the middle of the ear according to the grain position; these kernels were dried to a constant weight at 80°C. The dry matter accumulation of the kernel was estimated by fitting a logistic equation (Eq. 2) to the dry kernel weight plotted against the number of days after pollination during the grain-filling stage (Cao et al., 2008).

$$y=a/(1+b\times \text{exp}(-cx))(2)$$

where x is the number of days after pollination, y is the dry weight of one kernel (g), and a, b, and c are estimated parameters representing the ultimate growth mass, primary parameter, and growth rate parameter, respectively. The following equation (Eq. 3), derived by taking the first derivative of Eq. 2, was used to estimate the effective grain-filling duration and kernel growth rate:

$$dy/dx=a\times b\times c\times \text{exp}(-\text{cx})/{(1+b\times \exp \left(-\text{cx}\right))}^2(3)$$

The estimated parameters a, b, and c were used to calculate the grain-filling parameters. The day reaching the rate maximum grain-filling rate is T_max = ln⁡b/c, the increase in the weight of the kernel at the maximum grain-filling rate is W_max=a/2, the maximum grain-filling rate is G_max =a×c/2, the mean grain-filling rate is G_mean =c×W_max/6, the active grain-filling phase is P =  6/c, and the lag grain-filling period is T₁ = (ln⁡b−1.317)/c.

¹³C-Photosynthate Allocation Measurements Using ¹³C Isotope Labeling

To determine the ¹³C-photosynthate distribution among plant organs, 10 individual plants in each plot were selected for ¹³C labeling at the silking (R1) stage of maize. Carbon dioxide-¹³C (SRICI, Shanghai, China) with an enrichment >99 atom% was used for the tracing test. Ear leaves of selected plants were encased within 0.1-mm-thick mylar plastic bags, which allowed sunlight transmission at levels up to 95% of natural intensity. The bags were sealed with plastic tape at the base and then injected with 50 ml ¹³CO₂. The bags were removed after plant photosynthesis proceeded for 1 h. After 24 h, five labeled plants were sampled at ground level, and the remaining plants were collected at maturity (R6) stage. All labeled plants were dissected into stems, sheaths, ear leaves, other leaves, cobs and grains, oven-dried at 80°C to a constant weight, weighed to determine the dry matter, and then finely ground into powder. Subsequently, 5 mg of powder was used to determine the ¹³C isotopic abundance using an Isoprime 100 instrument (Isoprime 100, Cheadle, United Kingdom). The calculated ¹³C atom% excess in the plant organs refers to the distribution of ¹³C in the organs after ¹³CO₂ feeding, as described in previous studies (Wei et al., 2009; Liu et al., 2015).

The natural abundance of ¹³C% was calculated as follows (Eq. 4).

$${\mathrm{\delta }}^{13}C\%=(R_{sample}/R_{PBD}-1)\times 100(4)$$

where R_sample = ¹³C/¹²C_sample and R_PDB = ¹³C/¹²C_PDB, which has a value of 0.0112372. The labeled plant ¹³C abundance was calculated as follows (Eq. 5):

$${}^{13}C(\%)=({\mathrm{\delta }}^{13}C+100)\times R_{PBD}/$$

$$\times 100(5)$$

According to Eqs. 4 and 5, we calculated the ¹³C mg excess in each organ as follows (Eq. 6):

$${}^{13}Cmgexcess$$

$$=(dryweightofplanteachorgan\left(mg\right)\times C\%)\times$$

$$\left({}^{13}C\%-{\mathrm{\delta }}^{13}C\%\right)\times 100(6)$$

The distribution of the ¹³C excess ratio in each plant organ 24 h after ¹³CO₂ feeding and at maturity (R6) was calculated as follows (Eq. 7):

$${}^{13}C(\%){=}^{13}Cmgexcessineachplantorgan/$$

$$tota{l}^{13}Cmgexcessintheplant\times 100(7)$$

Dry Matter Accumulation and Grain Yield Measurement

In each plot, five individual plants were cut manually at ground level to determine the dry matter (DM) aboveground at the silking and maturity stages in 2015 and 2016. The harvested plants were dried initially at 105°C for 1 h to deactivate enzymes and then oven-dried at 80°C to a constant weight before weighing. The DM translocation efficiency of vegetative organs (%) (DMTE) was calculated for the postsilking stage following previous studies (Zhu et al., 2011) (Eq. 8).

$$DMTE(\%)=\mathrm{\hspace{0.33em}100}\times (DMallvegetativeorgansatsilking$$

$$-DMofallvegetativeorgansatmilkingormaturitystage)/$$

$$\mathrm{\hspace{0.33em}}DMallvegetativeorganssilking(8)$$

At physiological maturity, maize was harvested manually in a bordered area (5 m × 3.6 m) of each plot to determine grain yield. The grain yield was adjusted to 14% moisture content. The kernel number per ear was recorded, and five samples were randomly selected from each plot and then air-dried to record the 1000-kernel weight.

Statistical Analysis

Curve Expert 1.3 software was used to obtain the grain-filling logistic equation and parameters a, b, and c. Analysis of variance (ANOVA) was performed using SPSS 20.0 software (SPSS Institute Inc., United States), and figures were produced by Origin 9.0 and R 3.4.1.

Results

Characteristics of the Vascular Bundle in Internodes

The microscopic structure of vascular bundles from the basal stem, ear-located, peduncle and cob internodes of maize in 2016 is shown in Figures 1A,B. Generally, the numbers of both big and small vascular bundles gradually decreased from the basal stem to cob fractions (Figure 2). Planting density tended to decrease the numbers of both big and small vascular bundles in all internode fractions of maize. For instance, the numbers of big vascular bundles in cobs was decreased by an average of 28.8%, and small vascular bundles in peduncles was decreased by an average of 21.8% (Figures 2D,G). Compared with N0, increased N input significantly increased the number of vascular bundles, by an average of 20.7%, 7.9%, 18.2%, and 21.9% for big vascular bundles and 16.5%, 48.1%, 51.8%, and 25.7% for small vascular bundles in the basal stem, ear-located, peduncle and cob, respectively (P < 0.05).

FIGURE 1

Figure 1. (A) The micrograph of big vascular bundles structure at basal stem internode (a–f), ear located (g–i), peduncle (m–r), and cob (s–x) internodes. ND and HD indicate planting density of 67,500 and 90,000 plants ha^{− 1}; N0, N150, and N300 indicate N applied at 0, 150, and 300 kg ha^{− 1} levels, respectively. (B) The micrograph of small vascular bundles structure at basal-stem internode (a–f), ear located (g–i), peduncle (m–r), and cob (s–x) internodes. ND and HD indicate planting density of 67,500 and 90,000 plants ha^{− 1}; N0, N150, and N300 indicate N applied at 0, 150, and 300 kg ha^{− 1} levels, respectively.

FIGURE 2

Figure 2. Effect of planting density and nitrogen application on number of big (A–D) and small (E–H) vascular bundles at maize R3 stage in 2016. ND and HD indicate planting density of 67,500 and 90,000 plants ha^{− 1}; N0, N150, and N300 indicate N applied at 0 kg N ha^{− 1}, 150 kg N ha^{− 1}, and 300 kg N ha^{− 1} levels, respectively. Different letters indicate significant difference at 5% level. Values shown are the mean ± SE (n = 3).

Planting density and N fertilization level significantly influenced vascular bundles in the phloem and xylem, including the total area of the vascular bundles (TAVB) and the proportions of big vascular bundles (PBVB) and small vascular bundles (PSVB) relative to the total vascular bundle area (Table 1, P<0.05). The average values of PBVB in both the basal stem (0.36) and ear location (0.45) were lower than those of PSVB (0.64 and 0.55). However, the single average area of small vascular bundles decreased and the single average area of big vascular bundles increased at the cob internode, resulting in PBVB in both the peduncle (0.45) and cob (0.30) being higher than those of PSVB (0.55 and 0.70). Through increasing the area of the small vascular bundles, the PSVB was enhanced by N application increased, especially in the peduncle and cob internodes. In the peduncle and cob, the average single area of small vascular bundles increased by 53.4% and 102.6% while the average single area of vascular bundles increased by only 16.2% and 63.6%, respectively. Moreover, the average single area of small vascular bundles in the phloem increased by 34.8% in the peduncle and 35.0% in the cob while this area in the xylem increased by 27.4% in the peduncle and 18.4% in the cob.

TABLE 1

Table 1. Effect of planting density and nitrogen application on vascular bundle area of maize at milking stage in 2016.

Root-Bleeding Sap and Matter Transport Efficiency

Root-bleeding sap was significantly affected by planting density (D) and N fertilizer level (N), and MTE was significantly affected by the N application level and D × N (Figure 3, P < 0.05). Compared with N0, root-bleeding sap and MTE under both ND and HD conditions were significantly increased by N application, the averages increase of 93.0% and 89.3% for root-bleeding sap and 16.3% and 30.4% for MTE, respectively (P < 0.05). In particular, the highest MTE value (23.1 mg mm^–2 h^–1) was measured under HD in the N300 treatment. Compared with ND treatment, HD increased MTE by 11.6% under N300 conditions.

FIGURE 3

Figure 3. Response of root bleeding sap and matter transport efficiency on planting density and nitrogen fertilization level at maize R3 stage in 2016. ND and HD indicate planting density of 67,500 and 90,000 plants ha^{− 1}; N0, N150, and N300 indicate N applied at 0, 150, and 300 kg ha^{− 1} levels, respectively. Different letters indicate significant difference at 5% level. Values shown are the mean ± SE (n = 3). ns, no significance; *, **, and *** indicate significance at 5%, 1%, and 0.1% level, respectively.

Grain-Filling Characteristics

The factors year (Y), planting density (D) and N level (N) significantly affected the grain-filling parameters. including the kernel weight increment achieving maximum grain-filling rate (W_max), the day reaching the maximum grain-filling rate (T_max), the maximum filling rate (G_max), the active grain filling phase (P), the mean grain-filling rate (G_mean), and the lag period (T1) (Table 2). In addition, Tmax and T1 were significantly affected by D × N, and except for G_max, other grain-filling parameters were significantly affected by Y × D × N (P < 0.05). The G_max, P, and G_mean were significantly improved by N application increased, especially under HD conditions. Compared with N0, the G_max, P, and G_mean under HD conditions were average 35.2%, 8.1%, and 36.1% increased by N150 and N300 treatments, and these values were observed with 29.8%, 10.5%, and 30.0% increased under ND condition. Comparing those grain-filling parameters between ND and HD, the G_max, P, and G_mean under HD condition were average 8.1%, 1.8%, and 7.8% lower than those under ND condition in 2015 and 2016. Moreover, the effects of the effects of planting density and N input on the on the G_mean were more obvious.

TABLE 2

Table 2. Effect of planting density and nitrogen application on grain filling parameters of maize in 2015 and 2016.

According to the grain-filling parameters, the grain-filling process of maize was well simulated using logistic equations. In Figure 4, N input significantly increased the kernel weight (sink): in the initial grain-filling stages (the sink size was defined), the kernel weight increased by 42.7–46.3%, and the completion time of the lag period was 3–7 days earlier than that of the N0 treatment. HD resulted in a 1- to 5-day delay when the peak grain-filling rate was reached compared with ND.

FIGURE 4

Figure 4. Changes of grain filling rate and 100-kernel weight of maize affected by planting density and nitrogen fertilization level in 2015 and 2016. The reference line represents the number of days to reach the peak of grain filling. ND and HD indicate planting density of 67,500 and 90,000 plants ha^{− 1}; N0, N150, and N300 indicate N applied at 0, 150, and 300 kg ha^{− 1} levels, respectively. Values shown are the mean ± SE (n = 3).

¹³C-Photosynthate Allocation Between Maize Organs at Silking and Maturity

Generally, after ¹³C labeling for 24 h, the highest ratio of ¹³C-photosynthates in maize shoot at silking stage was recorded in stem, followed by other leaves and sheath, cob and ear leaves were found with very low ¹³C-photosynthates distribution (Figure 5). It was interesting to find that higher ¹³C-photosynthates distribution of N0 treatment in stem than those of N300 treatment under both ND and HD conditions. At maturity stage, the highest ratio of ¹³C-photosynthates allocation to grain from N300 treatment under HD condition were 48.3% and 57.3% in 2015 and 2016. Compared with N0, the ¹³C-photosynthates allocations to grain from N300 treatment were increased by 0.7–10.5% in 2015 and 2016.

FIGURE 5

Figure 5. ¹³C-photosynthates allocation between maize organs at R1 and R6 stages in 2015 and 2016. ND and HD indicate planting density of 67,500 and 90,000 plants ha^{− 1}; N0, N150, and N300 indicate N applied at 0, 150, and 300 kg ha^{− 1} levels, respectively. Different letters indicate significant difference at 5% level. Values shown are the mean ± SE (n = 3).

Yield Components and Dry Matter Accumulation

According to ANOVA results, the factors year (Y), planting density (D), N level (N) and their interactions (D × N and Y × D × N) had significant effects on the 1000-kernel weight (TKW) and grain yield (Table 3). N fertilizer input significantly increased the kernel number per ear (KNP) and TKW under both ND and HD (P < 0.05), by 23.7% and 38.8% for KNP and by 17.1% and 12.8% for TKW in 2015 and 2016, respectively. The highest grain yields were 11.0 Mg ha^–1 and 12.3 Mg ha^–1, obtained from the N300 treatment under HD conditions in 2015 and 2016, respectively.

TABLE 3

Table 3. Effect of planting density and nitrogen application on yield components of maize in 2015 and 2016.

Significant effects on total dry matter (DM) were obtained from both planting density and N level at the silking (R1), milking (R3) and maturity (R6) stages of maize, and the accumulation of DM and its translocation were significantly different under the different levels of N and its interaction with the planting density treatments. Total DM was significantly increased with increasing N application (P < 0.05), especially under HD, by 33.0% at R1, 59.0% at R3 and 87.1% at R6. Similarly, the Postsilking DM was significantly increased by N input, compared with N0, the postsilking DM, i.e., an increase of 218.1% in the accumulation of DM from R1 to R3 and of 161.4% from R3 to R6, and of the 170.3% from R1 to R6. However, N input decreased the dry matter translocation efficiency of vegetative organs (DMTE) from R1 to R3 and from R3 to R6 by 45.1% and 56.7%, respectively (Figure 6).

FIGURE 6

Figure 6. Effect of planting density and nitrogen application on the accumulation and translocation of dry matter during maize R1, R3, and R6 stages in 2015 and 2016. ND and HD indicate planting density of 67,500 and 90,000 plants ha^{− 1}; N0, N150, and N300 indicate N applied at 0 kg N ha^{− 1}, 150 kg N ha^{− 1}, and 300 kg N ha^{− 1} levels, respectively. Different letters indicate significant difference at 5% level. Values shown are the mean ± SE (n = 3); ns, no significance; *, **, and *** indicate significance at 5%, 1%, and 0.1% level, respectively.

Correlation Analysis

As shown in Figure 7, correlation analysis between the vascular bundle structure, MTE, the accumulation of postsilking dry matter (PDM) and mean grain-filling rate (G_mean) revealed that MTE was positively correlated with PDM and G_mean (R_MD = 0.86, R_MG = 0.82), and both PDM and G_mean were positively correlated with KNP and TKW (R_DKN = 0.85, R_DKW = 0.77; R_GKN = 0.90, R_GKW = 0.88). Correlation analysis showed that the proportion of the small vascular bundle area to the total vascular bundle area in the peduncle (PSVBP) and cob (PSVBC) and the corresponding number of vascular bundles in the peduncle (ASP) and cob (ASC) were positively correlated with MTE (R_PM = 0.70, R_CM = 0.78; R_APM = 0.67; R_ACM = 0.73). Additionally, in the phloem of the peduncle, the area of small vascular bundles (PASP) was positively correlated with TKW (R_SM = 0.82).

FIGURE 7

Figure 7. Correlation analysis of post-silking dry matter, yield components, vascular bundle structure and grain filling parameters under different planting density and nitrogen fertilization levels in 2016. PBVBP and PSVBC, the proportion of small vascular bundle area to total vascular bundle area in peduncle and cob; ASP and ASC, the number of small vascular bundle in peduncle and cob; PASP, the phloem area of small vascular bundle in peduncle; MTE, matter transport efficiency of basal-stem internode; G_mean, mean grain-filling rate; PDM, the accumulation of postsilking dry matter (R1–R6); TKW, 1000-kernel weight; KNP, kernel number per ear.

Discussion

Grain filling is an important process in determining grain yield that can explain yield changes, which influenced by agricultural management, such as N input and planting density (Xu et al., 2013; Wei et al., 2017, 2019). An appropriate increase in N application could improve the photosynthate supply and matter transport from vegetative to reproductive organs, contributing to the grain-filling rate, grain-filling phase and grain yield under high-density planting conditions (He et al., 2005; Ciampitti and Vyn, 2012; Hou et al., 2012). Grain filling is determined by the grain-filling rate and the active grain-filling phase (Zhou et al., 2017). Our results showed that the increase in plant density decreased the grain-filling rate without a significant change in the active filling phase. Increasing N fertilizer input increased both the grain-filling rate and active filling phase. In particular, the increase in the maize grain-filling rate tended to be stronger with increased N input under high planting density (HD). The first grain-filling period (the lag grain-filling period) determines the potential kernel weight size (Hisse et al., 2019). Our study showed that in this period, the effect of increased N input on the kernel weight (i.e., the demand of the sink) significantly increased and advanced the completion of the lag grain-filling period, especially under HD conditions (Figure 4). The grain sink demand also adjusted the ¹³C-photosynthate transport ratio of vegetative organs and translocation from leaves to grains (Figure 5). Therefore, the effect of increased N fertilizer input under high planting density might occur through a change in the grain-filling rate and an increase in the potential to increase grain weight, not an extension of the active period proposed in a previous study (Poneleit and Egli, 1979). Increased N input significantly accelerated the activation of grain filling and the sink demand so that the kernel weight significantly increased during the initial grain-filling stages.

Additionally, previous studies pointed out that the grain-filling rate is chiefly restricted by insufficient photosynthate supply and sink capacity limitation (Fischer and Wilson, 1975; Evans, 1993; Gebbing and Schnyder, 1999). However, other studies have clarified that the source and sink were not the dominant limiting factors for the grain filling rate of maize under appropriate N input and without abiotic stresses, but the dominant limiting factors was the dry matter (DM) accumulation and how much of assimilation distribute to each organ (Pan et al., 2011; Zhu et al., 2011; Gasura et al., 2013; Shen et al., 2017). Additionally, previous study indicated that the assimilation source of grain filling mainly came from the current assimilation post-silking and the vegetative organ stored in the pre-silking (Fu et al., 2011; Tan et al., 2020). In this study, although translocation efficiency of vegetative organs at the at the silking (R1) stage to the milking (R3) stages was higher than that at the R3 to maturity (R6) stages, the grain translocation efficiency was only 2.3% (R1 to R3) to 17.2% (R3 to R6) during the whole grain-filling period of maize, and with increased N input, this efficiency decreased. Moreover, compared with N0, increased N input significantly enhanced the accumulation of postsilking DM (Figure 6) as well as the MTE (Figure 3), particularly under HD, increased 170.3% for the postsilking DM (R1 to R6) and increased 30.4% for the MTE. These result indication that the higher grain filling rate and grain yield may were attributed to the accumulation of DM and the MTE (i.e., current assimilation post-silking transport) rate rather than the vegetative organ stored in pre-silking, because of the low vegetative organ translocation efficiency cannot meet the assimilation required for grain filling rate and grain yield formation.

As important channels of matter transport in crops, vascular bundles can increase the accumulation of postsilking DM, adjusting the distribution of photosynthetic products in various organs (Shane et al., 2000; Noguchi et al., 2005). Previous studies showed that N application increased the number and area of small vascular bundles in the phloem in certain positions (panicle neck or cob), which was beneficial for matter transport and grain yield improvement (He et al., 2005; Terao et al., 2010). However, increased N input how to change the structure of the vascular bundle and the proportion of the phloem in the structure of the vascular bundle, and the relationship between the changed structure and the grain filling is still unclear. The results from this study suggest that the increase in N application significantly increased the number and area of small vascular bundles, which led to an increase in the proportion of the small vascular bundle area to the total vascular bundle area (PSVB), especially in peduncles and cobs (Table 1). Furthermore, N application coordinated the relationship between the proportion of and small vascular bundle areas to the total vascular bundle area and the MTE, which stimulated ¹³C-photosynthate allocation to grains to meet the sink (kernel) demand (Figure 5). Additionally, the accumulation of postsilking DM and grain-filling rate were improved by increasing N application, which could be explained by the change in the MTE. A positive relationship was found between MTE and vascular bundles in this study (Figure 7). These results suggest that PSVB in peduncles and cobs improved through an increase in the N supply, which was advantageous for MTE, contributing to the accumulation of DM, the grain-filling rate and the grain yield of maize under intensive planting.

Conclusion

The increase in N application significantly improved the mean grain-filling rate, postsilking dry matter of maize and ¹³C-photosynthate allocation to grains at maturity and optimized the transport system of vascular bundles between organs as well, further coordinating the relationship between the big and small vascular bundles in peduncles and cobs, which greatly enhanced the MTE after maize silking (16.3–30.4%). These results suggest that increasing N input mainly optimized the vascular bundle structure of the ear system which improvement the accumulation of postsilking dry matter and MTE to enhance the mean grain filling and grain yield in maize under high planting density. These results might provide pivotal physiological evidence for how N management changes the vascular bundle transport system and clarifying the effect of the matter transport system on grain filling and grain yield under intensive planting conditions.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author Contributions

HR collected the samples, analyzed the samples, and wrote the manuscript. YJ and MZ made a contribution to writing and editing of the manuscript. HQ and CL made a contribution on to design of the work, analysis, and revised the manuscript. All authors read and approved the article.

Funding

This study was financially supported by the National Key Research and Development Program of China (2016YFD0300103), the National Natural Science Foundation of China (31971852), and the National Modern Agriculture Industry Technology System (CARS-02-12).

Conflict of Interest

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.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2020.602739/full#supplementary-material

References

Austin, R. B., Edrich, J. A., Ford, M. A., and Blackwell, R. D. (1997). The fate of the dry matter, carbohydrates and ¹⁴C lost from the leaves and stems of wheat during grain filling. Ann. Bot. 41, 1309–1321. doi: 10.1093/oxfordjournals.aob.a085419

CrossRef Full Text | Google Scholar

Borrás, L., Slafer, G. A., and Otegui, M. E. (2004). Seed dry weight response to source–sink manipulations in wheat, maize and soybean: a quantitative reappraisal. Field Crops Res. 86, 131–146. doi: 10.1016/j.fcr.2003.08.002

CrossRef Full Text | Google Scholar

Cao, C. Y., Li, K. J., Cui, Y. H., Ma, J. Y., and Zheng, C. L. (2008). Modeling the effect of long-term fertilization on grain filling of summer maize. Plant Nutr. Fert. Sci. 14, 48–53. doi: 10.3321/j.issn:1008-505X.2008.01.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X. P., Cui, Z. L., Vitousek, P. M., Cassman, K. G., Matson, P. A., Bai, J. S., et al. (2011). Integr-ated soil–crop system management for food security. Proc. Natl. Acad. Sci.U.S.A. 108, 6399–6404. doi: 10.1073/pnas.1101419108

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y. L., Xiao, C. X., Wu, D., Xia, T. T., Chen, Q. W., Chen, F. J., et al. (2015). Effects of nitrogen application rate on grain yield and grain nitrogen concentration in two maize hybrids with contrasting nitrogen remobilization efficiency. Eur. J. Agron. 62, 79–89. doi: 10.1016/j.eja.2014.09.008

CrossRef Full Text | Google Scholar

Ciampitti, I. A., and Vyn, T. J. (2012). Physiological perspectives of changes over time in maize yield dependency on nitrogen uptake and associated nitrogen efficiencies: a review. Field Crops Res. 133, 48–67. doi: 10.1016/j.fcr.2012.03.008

CrossRef Full Text | Google Scholar

Cirilo, A. G., and Andrade, F. H. (1996). Sowing date and kernel weight in maize. Crop Sci. 36:325. doi: 10.2135/cropsci1996.0011183X003600020019x

CrossRef Full Text | Google Scholar

D’Andrea, K. E., Piedra, C. V., Mandolino, C. I., Bender, R., Cerri, A. M., Cirilo, A. G., et al. (2016). Contribution of reserves to kernel weight and grain yield determination in maize: phenotypic and genotypic variation. Crop Sci 56:697. doi: 10.2135/cropsci2015.05.0295

CrossRef Full Text | Google Scholar

Di Salvo, L. P., Cellucci, G. C., Carlino, M. E., and García de Salamone, I. E. (2018). Plant growth promoting rhizobacteria inoculation and nitrogen fertilization increase maize (Zea mays L.) grain yield and modified rhizosphere microbial communities. Appl. Soil Ecol. 126, 113–120. doi: 10.1016/j.apsoil.2018.02.010

CrossRef Full Text | Google Scholar

Evans, L. T. (1993). Crop Evolution, Adaptation and Yield. Cambridge: Cambridge University Press, 75.

Google Scholar

Fischer, K. S., and Wilson, G. L. (1975). Studies of grain production in Sorghum bicolor (L. Moench). The relative importance of assimilate supply, grain growth capacity and transport system. Aust. J. Agric. Res 26, 11–23. doi: 10.1071/ar9750031

CrossRef Full Text | Google Scholar

Fu, J., Huang, Z. H., Wang, Z. Q., Yang, J. C., and Zhang, J. H. (2011). Pre-anthesis non-structural ca-rbohydrate reserve in the stem enhances the sink strength of inferior spikelets during grain filling of rice. Field Crops Res. 123, 170–182. doi: 10.1016/j.fcr.2011.05.015

CrossRef Full Text | Google Scholar

Gambín, B. L., Borrás, L., and Otegui, M. E. (2006). Source–sink relations and kernel weight differences in maize temperate hybrids. Field Crops Res. 95, 316–326. doi: 10.1016/j.fcr.2005.04.002

CrossRef Full Text | Google Scholar

Gan, S., and Amasino, R. M. (1997). Making sense of senescence molecular genetic regulation and manipulation of leaf senescence. Plant Physiol. 113, 313–319. doi: 10.1104/pp.113.2.313

PubMed Abstract | CrossRef Full Text | Google Scholar

Gasura, E., Setimela, P., Edema, R., Gibson, P. T., Okori, P., and Tarekegne, A. (2013). Exploiting grain-filling rate and effective grain-filling duration to improve grain yield of early-maturing maize. Crop Sci. 53:2295–2303. doi: 10.2135/cropsci2013.01.0032

CrossRef Full Text | Google Scholar

Gebbing, T., and Schnyder, H. (1999). Pre-Anthesis reserve utilization for protein and carbohydrate synthesis in grains of wheat. Plant Physiol. 121, 871–878. doi: 10.1104/pp.121.3.871

PubMed Abstract | CrossRef Full Text | Google Scholar

He, Q. P., Dong, S. T., and Gao, R. Q. (2005). Relationship between development of spike vascular bundle and sink capacity of ear and kernel in maize (Zea mays L.). Acta Agron. Sin. 31, 995–1000. (in Chinese with English abstract),Google Scholar

Hisse, I. R., D’Andrea, K. E., and Otegui, M. E. (2019). Source-sink relations and kernel weight in maize inbred lines and hybrids: responses to contrasting nitrogen supply levels. Field Crops Res. 230, 151–159. doi: 10.1016/j.fcr.2018.10.011

CrossRef Full Text | Google Scholar

Hou, P., Gao, Q., Xie, R. Z., Li, S. K., Meng, Q. F., Kirkby, E. A., et al. (2012). Grain yields in relation to N requirement: optimizing nitrogen management for spring maize grown in China. Field Crops Res. 129, 1–6. doi: 10.1016/j.fcr.2012.01.006

CrossRef Full Text | Google Scholar

Housley, T. L., and Peterson, D. M. (1982). Oat stem vascular size in relation to kernel number and w-eight. I. Controlled environment1. Crop Sci. 22, 259–263. doi: 10.2135/cropsci1982.0011183X002200020014x

CrossRef Full Text | Google Scholar

Jones, R. J., Schreiber, B. M. N., and Roessler, J. A. (1996). Kernel sink capacity in maize: genotypic and maternal regulation. Crop Sci. 36, 301–306. doi: 10.2135/cropsci1996.0011183X003600020015x

CrossRef Full Text | Google Scholar

Khan, N. A., Murayama, S., Ishimine, Y., Tsuzuki, E., and Nakamura, I. (1998). Physio-morphological studies of F1 hybrids in rice (Oryza sativa L.): photosynthetic ability and yield. Plant Prod. Sci. 1, 233–239. doi: 10.1626/pps.1.233

CrossRef Full Text | Google Scholar

Liu, T. N., Gu, L. M., Dong, S. T., Zhang, J. W., Liu, P., and Zhao, B. (2015). Optimum leaf removal increases canopy apparent photosynthesis, ¹³C-photosynthate distribution and grain yield of maize crops grown at high density. Field Crops Res. 170, 32–39. doi: 10.1016/j.fcr.2014.09.015

CrossRef Full Text | Google Scholar

Monneveux, P., Zaidi, P. H., and Sanchez, C. (2005). Population density and low nitrogen affects yield associated traits in tropical maize. Crop Sci. 45, 535–545. doi: 10.2135/cropsci2005.0535

CrossRef Full Text | Google Scholar

Mu, X. H., Chen, Q. W., Chen, F. J., Yuan, L. X., and Mi, G. H. (2016). Within-leaf nitrogen allocation in adaptation to low nitrogen supply in maize during grain-filling stage. Front. Plant Sci. 7:699. doi: 10.3389/fpls.2016.00699

PubMed Abstract | CrossRef Full Text | Google Scholar

Noguchi, A., Kageyama, M., Shinmachi, F., Schmidhalter, U., and Hasegawa, I. (2005). Potential for using plant xylem sap to evaluate inorganic nutrient availability in soil. Soil Sci. Plant Nutr. 51, 333–341. doi: 10.1111/j.1747-0765.2005.tb00038.x

CrossRef Full Text | Google Scholar

Ntanos, D., and Koutroubas, S. (2002). Dry matter and N accumulation and translocation for Indica and Japonica rice under mediterranean conditions. Field Crops Res. 74, 93–101. doi: 10.1016/s0378-4290(01)00203-9

CrossRef Full Text | Google Scholar

Ouattar, S., Jones, R. J., Crookston, R. K., and Kajeiou, M. (1987). Effect of drought on water relations of developing maize kernels1. Crop Sci. 27, 730–735. doi: 10.2135/cropsci1987.0011183x002700040

CrossRef Full Text | Google Scholar

Pan, J. F., Cui, K. H., Wei, D., Huang, J. L., Xiang, J., and Nie, L. X. (2011). Relationships of non-structural carbohydrates accumulation and translocation with yield formation in rice recombinant inbred lines under two nitrogen levels. Physiol. Plant 141, 321–331. doi: 10.1111/j.1399-3054.2010.01441.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Piao, L., Qi, H., Li, C. F., and Zhao, M. (2016). Optimized tillage practices and row spacing to improve grain yield and matter transport efficiency in intensive spring maize. Field Crops Res. 198, 258–268. doi: 10.1016/j.fcr.2016.08.012

CrossRef Full Text | Google Scholar

Poneleit, C. G., and Egli, D. B. (1979). Kernel growth rate and duration in maize as affected by plant density and genotype. Crop Sci 19, 385–388. doi: 10.2135/cropsci1979.0011183x001900030027x

CrossRef Full Text | Google Scholar

Rajcan, I., and Tollenaar, M. (1999). Source: sink ratio and leaf senescence in maize:I. Dry matter accumulation and partitioning during grain filling. Field Crops Res. 60, 245–253. doi: 10.1016/s0378-4290(98)00142-7

CrossRef Full Text | Google Scholar

Ren, B. Z., Dong, S. T., Zhao, B., Liu, P., and Zhang, J. W. (2017). Responses of nitrogen metabolism, uptake and translocation of maize to waterlogging at different growth stages. Front. Plant Sci. 8:1216. doi: 10.3389/fpls.2017.01216

PubMed Abstract | CrossRef Full Text | Google Scholar

Shane, M. W., McCully, M. E., and Canny, M. J. (2000). The vascular system of maize stems revisited: implications for water transport and xylem safety. Ann. Bot. 86, 245–258. doi: 10.1006/anbo.2000.1171

CrossRef Full Text | Google Scholar

Shen, L. X., Huang, Y. K., and Li, T. (2017). Top-grain filling characteristics at an early stage of maize (Zea mays L.) with different nitrogen use efficiencies. J. Integr. Agric. 16, 626–639. doi: 10.1016/s2095-3119(16)61457-0

CrossRef Full Text | Google Scholar

Shi, J. G., Cui, H. Y., Zhao, B., Dong, S. T., Liu, P., and Zhang, J. W. (2013). Effect of light on yield and characteristics of grain-filling of summer maize from flowering to maturity. Sci. Agric. Sin. 46, 4427–4434. (in Chinese with English abstract),Google Scholar

Tan, D., Guo, L. L., Liu, J. M., Fan, Y. L., and Li, Q. Q. (2020). Response of dry matter translocation and grain yield of summer maize to biodegradable film in the North China Plain. Int. Agrophys. 34, 87–94. doi: 10.31545/intagr/112269

CrossRef Full Text | Google Scholar

Terao, T., Nagata, K., Morino, K., and Hirose, T. (2010). A gene controlling the number of primary rachis branches also controls the vascular bundle formation and hence is responsible to increase the harvest index and grain yield in rice. Theor. Appl. Genet. 120, 875–893. doi: 10.1007/s00122-009-1218-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Tilman, D., Cassman, K. G., Matson, P. A., Naylor, R., and Polasky, S. (2002). Agricultural sustainability and intensive production practices. Nature 418, 671–677. doi: 10.1038/nature01014

PubMed Abstract | CrossRef Full Text | Google Scholar

Tollenaar, M., Deen, W., Echarte, L., and Liu, W. (2006). Effect of crowding stress on dry matter accumulation and harvest index in maize. Agron. J. 98, 930–937. doi: 10.2134/agronj2005.0336

CrossRef Full Text | Google Scholar

Wang, F., Cheng, F. M., and Zhang, G. P. (2006). The relationship between grain filling and hormone content as affected by genotype and source–sink relation. Plant Growth Regul. 49, 1–8. doi: 10.1007/s10725-006-0017-3

CrossRef Full Text | Google Scholar

Wang, P., Wang, Z. K., Sun, X. C., Mu, X. H., Chen, H., Chen, F. J., et al. (2019). Interaction effect of nitrogen form and planting density on plant growth and nutrient uptake in maize seedlings. J. Integr. Agric 18, 1120–1129. doi: 10.1016/s2095-3119(18)61977-x

CrossRef Full Text | Google Scholar

Wei, S. S., Wang, X. Y., Li, G. H., Qin, Y. Y., Jiang, D., and Dong, S. T. (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

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, S. S., Wang, X. Y., Zhu, Q., Jiang, D. C., and Dong, S. T. (2017). Optimizing yield and resource utilization of summer maize under the conditions of increasing density and reducing nitrogen fertilization. Sci. Nat. 104, 86. doi: 10.1007/s00114-017-1509-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, X. W., Wei, L., Hao, H. L., Ying, X. C., Medha, D., and Janice, T. (2009). Use of ¹³C labeling to assess carbon partitioning in transgenic and nontransgenic (parental) rice and their rhizosphere soil microbial communities. FEMS Microbiol. Ecol. 67, 93–102. doi: 10.1111/j.1574-6941.2008.00599

CrossRef Full Text | Google Scholar

Wigoda, N., Moshelion, M., and Moran, N. (2014). Is the leaf bundle sheath a “smart flux valve” for K⁺ nutrition? J. Plant Physiol. 171, 715–722. doi: 10.1016/j.jplph.2013.12.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, T. F., Lv, T., Zhao, J. R., Wang, R. H., Chen, C. Y., Liu, Y. E., et al. (2016). Response of grain filling characteristics of maize to sowing date. Chin. J. Appl. Ecol. 27, 2513–2519. (in Chinese with English abstract),Google Scholar

Xu, Y. J., Gu, D. J., Zhang, B. B., Zhang, H., Wang, Z. Q., and Yang, J. C. (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

CrossRef Full Text (in Chinese with English abstract),Google Scholar

Yan, S. C., Wu, Y., Fan, J. L., Zhang, F. C., Qiang, S. C., Zheng, J., et al. (2019). Effects of water and fertilizer management on grain filling characteristics, grain weight and productivity of drip-fertigated winter wheat. Agric. Water Manage. 213, 983–995. doi: 10.1016/j.agwat.2018.12.019

CrossRef Full Text | Google Scholar

Yang, J. C., Zhang, J. H., Huang, Z. L., Wang, Z. Q., Zhu, Q. S., and Liu, L. J. (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

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J. C., Zhang, J. H., Wang, Z. Q., Zhu, Q. S., and Liu, L. J. (2001). Water deficit–induced senescence and its relationship to the remobilization of pre-stored carbon in wheat during grain filling. Agron. J. 93, 196–206. doi: 10.2134/agronj2001.931196x

CrossRef Full Text | Google Scholar

Zhong, R., Taylor, J. J., and Ye, Z. H. (1997). Disruption of interfascicular fiber differentiation in an Arabidopsis mutant. Plant Cell 9, 2159–2159. doi: 10.2307/3870576

CrossRef Full Text | Google Scholar

Zhou, B. Y., Yue, Y., Sun, X. F., Ding, Z. S., Ma, W., and Zhao, M. (2017). Maize kernel weight responses to sowing date-associated variation in weather conditions. Crops J. 5, 43–51. doi: 10.1016/j.cj.2016.07.002

CrossRef Full Text | Google Scholar

Zhu, M., Ge, Y. X., Li, F. H., Wang, Z. B., Wang, H. W., and Shi, Z. S. (2011). Accumulation and translocation of dry matter and nitrogen in different purple corn hybrids (Zea mays L.). Afr. J Agric. Res. 6, 2820–2827. doi: 10.5897/AJAR10.1164

CrossRef Full Text | Google Scholar