Optimizing plant density and nitrogen application to manipulate tiller growth and increase grain yield and nitrogen-use efficiency in winter wheat

The growth of wheat tillers and plant nitrogen-use efficiency (NUE) will gradually deteriorate in response to high plant density and over-application of N. Therefore, in this study, a 2-year field study was conducted with three levels of plant densities (75 ×104plants ha−1, D1; 300 ×104plants ha−1, D2; 525 ×104plants ha−1, D3) and three levels of N application rates (120 kg N ha−1, N1; 240 kg N ha−1, N2; 360 kg N ha−1, N3) to determine how to optimize plant density and N application to regulate tiller growth and to assess the contribution of such measures to enhancing grain yield (GY) and NUE. The results indicated that an increase in plant density significantly increased the number of superior tillers and the number of spikes per m2(SN), resulting in a higher GY and higher partial factor productivity of applied N (PFPN). However, there was no significant difference in GY and PFPN between plant densities D2 and D3. Increasing the N application rate significantly increased the vascular bundle number (NVB) and area (AVB), however, excess N application (N3) did not significantly improve these parameters. N application significantly increased GY, whereas there was a significant decrease in PFPN in response to an increase in N application rate. The two years results suggested that increasing the plant density (from 75 ×104plants ha−1to 336 ×104plants ha−1) in conjunction with the application of 290 kg N ha−1N will maximize GY, and also increase PFPN(39.7 kg kg−1), compared with the application of 360 kg N ha−1N. Therefore, an appropriate combination of increased planting density with reduced N application could regulate tiller number and favor the superior tiller group, to produce wheat populations with enhanced yield and NUE.


INTRODUCTION
Wheat is one of the most important food crops worldwide. In China, where wheat has accounted for more than 20% of the total sowing area in recent years, production has substantially increased during the past few decades (Chen et al., 2017). Nonetheless, owing to continual population growth, the demand for wheat will continue to increase. Moreover, arable areas have undergone a substantial decrease in extent due to urban expansion and environmental degradation (Wei et al., 2016), and wheat production currently faces severe constraints and difficult challenges (Li et al., 2016a). Therefore, it is essential to ensure food security through increasing yield per unit area, to enhance the total production on a diminishing area on cultivated farmland.
Tillering is an important agronomic trait in graminaceous crops such as wheat and rice, as it determines final spike number per m 2 and plays an important role in determining grain yield (Naruoka et al., 2011;Wang et al., 2016b). When fewer tillers are produced, a smaller population is formed, resulting in yield reductions due to a smaller population of spikes per m 2 (Mitchell et al., 2012;Moeller, Evers & Rebetzke, 2014). Conversely, an excess of surviving tillers can lead to a larger population, resulting in an increased risk of lodging risk due to reduced culm quality (Wang et al., 2009;Zheng et al., 2017). In addition, individual tillers exhibit heterogeneity. In this regard, Cai et al. (2014) reported decreases in grain number and kernel weight per spike as the tiller position shifts from low to high, and the flag leaves of low-position tillers have been shown to have a higher photosynthetic rate than high-position tillers (Xu et al., 2015a). There is also competition between lowand high-position tillers for limiting resources (Elhani et al., 2007;Ma et al., 2008). Earlyemerging tillers generally intercept more radiation and shade the late-emerging tillers, and consequently dominate the asymmetric competition from the early growth stage (Moeller, Evers & Rebetzke, 2014). Late-emerging tillers generally make a minimal contribution to grain yield formation (Otteson et al., 2008), and indeed, these unproductive tillers can negatively affect grain yield by competing for resources such as nitrogen (N) and solar radiation (Richards et al., 2002;Berry et al., 2003). Hence, management of tiller emergence and heading, and utilization of early-emerging tillers, are key factors for the establishment of high-yield and high-efficiency populations.
The occurrence of tillers generally involves two developmental stages: the formation of axillary meristems in the leaf axil and the growth of tillers (Li et al., 2003;Bennett & Leyser, 2006). Tiller development is a result of complex interactions between endogenous signals and environmental factors. Internal factors regulating tiller emergence include numerous genetic factors (Xie, Mayes & Sparkes, 2016;Hyles et al., 2017), and endogenous hormones (Lin et al., 2016). External factors including light, plant density, and fertilizer application, have a significant effect on the growth of tillers (Assuero & Tognetti, 2010). The ratio of red to far-red light irradiance influences the outgrowth of tiller buds (Evers et al., 2006), whereas increasing planting density from 135 to 405 plants per m 2 has been shown to significantly increase grain yield (Dai et al., 2014), although further increases in density do not affect grain yield. Although seeding rate has been shown to have no influence on main leaf development, it does have a considerable effect on tiller development (Chen et al., 2008). Increasing the seeding rate has been found to increase the proportion of yield obtained from the main spike and decreases the proportion from high-position tiller spikes (Otteson et al., 2008). Although low plant density produces a higher grain number and grain weight per spike (Li et al., 2016b), this is generally not sufficient to compensate for the lower spike density per m 2 generated by a lower tiller density (Clerget et al., 2016). Therefore, an appropriate increase in plant density to balance yield component factors would appear to be an apposite agronomic management strategy for enhancing grain yield.
Nitrogen is an important and essential nutrient for all plants. N supply has considerable effects on plant growth in terms of the amount of biomass produced, the size and proportion of organs and their structure, and the progress of plant development (Lawlor, Lemaire & Gastal, 2001). N regulates rice tiller bud growth through regulating endogenous hormones and N metabolism (Liu et al., 2011). Nitrogen deficiency is associated with a reduced rate of leaf appearance and reduces the maximum rate of emergence for each tiller (Prystupa, Slafer & Savin, 2003), whereas an increase in N application rate increases tiller density and reduces tiller mortality (Ahmed et al., 2016). Increasing N application increases yield through an increase in grain yield in early-emerging tillers (Wang et al., 2016a). Furthermore, grain number per unit area is significantly increased in response to N application, mainly due to an increase in the contribution of grains produced by tillers (Terrile, Miralles & González, 2017). In China, farmers excessively apply N fertilizers because of their hope to sustain further grain yield increases, but grain yield does not keep synchronous increase with excessive N application (Meng et al., 2016). N application during the wheat growing season generally exceeds 320∼350 kg N ha −1 ; however, some farmers uses rates as high as 750 kg N ha −1 (Cui, Chen & Zhang, 2010;Lu et al., 2015;Zhang et al., 2015). Over-application of N contributes little to enhancing grain yield but can reduce nitrogen-use efficiency (NUE) and enhances the risk of environmental pollution (Cui, Chen & Zhang, 2010;Zhang et al., 2015;Mon et al., 2016;Tian et al., 2017). Therefore, achieving both high yield and high NUE simultaneously is a major challenge. Lu et al. (2014) have reported that optimal N management can optimize population quantity and quality to achieve high grain yield, and therefore a high-yielding wheat population can be obtained by amending N supply to regulate tiller development. Accordingly, N application rates should be calculated to produce the required quantity of tillers and modified over time to generate wheat crops with the desired NUE.
Although previous studies have either reported the effects of plant density on tiller development (e.g., tiller emergence, tiller mortality, and tiller production) or focused on N demand and supply, there have been few studies that have systematically investigated the relationship between these aspects, for example, the mechanisms underlying the differences in grain yield among different tiller spikes and how to manipulate tillers to establish high-yield and high-efficiency groups. Moreover, little is known regarding how plant density and N application interactively regulate tiller growth and wheat population development to produce higher grain yield and NUE. Therefore, the objectives of this study were as follows: (i) to quantify the variation between tillers regarding yield formation at different plant densities and N levels; and (ii) to build a model to determine the optimal planting density and N application strategies for realizing a high-yielding and high N efficient population structure.

Plant materials and growth conditions
The field experiments conducted in this study were carried out over two winter wheat growing seasons (2013-2014 and 2014-2015), at the experimental station of Shandong Agricultural University, Tai 'an,China (36 • 09 N,117 • 09 E, 128 m above sea level). This region has a warm and semi-humid continental monsoon climate, with an average total annual solar irradiance of 5.08×10 6 kJ cm −2 , an average annual temperature of 13.7 • C, and an average annual rainfall of 631.5 mm. The rainfall and mean temperature distributions of winter wheat growing stages are shown in Fig. S1. The soil in the study area is classified as a Eutric Cambisol according to the World Reference Base for Soil Resources (2014). The top 30 cm of the soil contained 14.7 g kg −1 organic matter, 1.24 g kg −1 total N, 87.2 mg kg −1 available N, 9.6 mg kg −1 available P, and 85.3 mg kg −1 available K. The winter wheat cultivar used for the purposes of this study was Jimai 22 (JM22), which was grown in the experimental plots. Seeds were sown on October 10, 2013 and October 8, 2014, and plants were respectively harvested on June 9, 2014 and June 8, 2015. Disease, pests, and weeds in each treatment were well controlled by managers.

Treatments and experimental design
The experiments were laid out in a two-factor completely randomized design with three replicates. The main plots were assigned to three plant densities (75 × 10 4 plants ha −1 , D 1 ; 300 × 10 4 plants ha −1 , D 2 ; and 525 × 10 4 plants ha −1 , D 3 ). Subplots were assigned to three urea (N) fertilizer application rates (120 kg N ha −1 , N 1 ; 240 kg N ha −1 , N 2 ; and 360 kg N ha −1 , N 3 ). Each replicate plot (3 m × 3 m) consisted of 10 rows of wheat and two ridges. Row spacing, the width of rides, and a diagram of plant densities are shown in Fig. S2. Half the amount of N, 75 kg ha −1 P 2 O 5 , and 150 kg ha −1 K 2 O were mixed into the soil before planting the wheat. The remaining amount of N was applied at the jointing stage in each treatment.

Observations of wheat population dynamics
Uniform plants in a 1 m 2 area per plot were selected and tagged with labels for observations of population dynamics at the following stages: three-leaf, over-wintering, stem elongation, heading time (50% of plants headed), and maturity. Twenty plants were randomly selected and tagged with labels in each treatment for observation of tiller number at over-wintering, stem elongation, heading time, and maturity. The main stem was denoted as 0; primary tillers on the main stem in emergence order were referred to as I, II, III, IV, and V. Secondary tillers on the primary tillers in emergence order were referred to as I 1 and II 1 (Fig. S3). The effective tiller rate (ETR) was calculated using the equation: ETR = the number of tiller spikes per plant at maturity/the maximal number of tillers per plant ×100%. The ineffective tiller rate (ITR) = 100-ETR.

Assessment of the grain-filling process
From 4 days after anthesis (DAA), 30 plants from each treatment in the field experiments were sampled at 4-day intervals until 36 DAA. Spikes were dried at 70 • C to constant weight, dehulled, and weighed. These data were used to characterize the grain-filling process using Richards' growth equation, as described by Yang et al. (2004): Grain filling rate (G) was calculated as the derivative of Eq. (1): Integration of Eq.
(2) gave the mean grain-filling rate (Eq. (3)) and the maximum grainfilling rate (Eq. (4)): where W is the kernel weight (g ), A is the final kernel weight (g ), t is the time after anthesis (d), and B,k, and N are the coefficients determined by regression. The active grain-filling period (T ) was defined as the period during which W constituted from 5% (t 1 ) to 95% (t 2 ) of A.

Assessment of first internode microstructure
The microstructure of the first internode was assessed according to the methods described by Zheng et al. (2017). Fifteen days after anthesis, three plants from each treatment were sampled for microstructural observations. The middle section of the first internode (2 cm in length) of plants from each treatment was fixed in formalin:acetic acid:ethanol for 24 h, dehydrated with ethanol, and embedded in paraffin. Approximately 4-µm-thick sections were obtained using a microtome (Leica, Wetzlar, Germany), and these were sequentially stained with 1% safranin and 0.5% fast green. The stained cross-sections were observed and photographed using a Nikon DS-V3 microscope (Nikon, Tokyo, Japan). The number and area of vascular bundles at the internode were measured and calculated using Image-Pro Plus (Version 6.0; Media Cybernetics, Rockville, MD, USA).

Measurements of single-stem biological yield and grain yield and its components
At maturity, the 20 tagged plants in each treatment were selected and divided into different spikes according to tillering order. The grains from individual spikes were oven dried at 70 • C to constant weight for determinations of grain number per spike and yield per spike. Euclidean distances were used to identify dissimilarities between the different tillers based on grain number and yield per spike. Spikes in an area of 1 m 2 from which no spikes had been sampled were harvested by hand to determine grain yield, number of kernels per spike, and 1,000-grain weight. Each measurement was performed on plants from three different plots.

Polynomial regression with response surface analysis and evaluation of nitrogen-use efficiency
Trend surface simulation was used to analyze the effect of plant densities and N application rates on grain yield. The quadratic polynomial trend surface equation is as follows: where z is the grain yield (kg ha −1 ), x is plant density (plant ha −1 ), y is the N application rate (kg N ha −1 ), and a 0 , a 1 , a 2 , a 3 , a 4 and a 5 are the coefficients determined by regression. These coefficients were calculated by the least square method.
The partial factor productivity for the applied N (PFP N ) was used to indicate NUE. PFP N was calculated using the following equation: where Y N is grain yield in the N-applied plot and N is the nitrogen application rate.

Statistical analysis and processing
Statistical analyses were carried out using Data Processing System software version 7.05 (DPS, Hangzhou, China). A multivariate ANOVA was conducted to determine the mean squares, degrees of freedom, and significance levels. Means were compared using LSD test and differences were considered significant at P < 0.05. The data for changes in the wheat population, effective tiller rate, grain number and yield per spike and the partial factor productivity for the applied N were averaged from the 2 years of the study. The trend surface simulation was analyzed using SPSS Statistics 24. Graphs were plotted using SigmaPlot 10.

Results of variance analysis of the year (Yr), plant density (D), nitrogen application rate (N), and their interactions
Variance analysis was conducted using DPS7.05 to find the mean squares and significance (Tables 1 and 2). Gain yield (GY), spike number per m 2 (SN), 100-grain weight (TGW), the effective tiller rate (ETR), partial factor productivity of applied N (PFP), yield per spike, grain number per spike, the maximum grain-filling rate (Gmax); the mean grain-filling rate (Gmean), the active grain-filling period (T), and the associated vascular bundle parameters showed no significant effects of year (Yr), the interaction between Yr, plant density (D) and nitrogen application rate (N). In contrast, GY, SN, TGW, ETR, PFP, and the associated grain-filling parameters were significantly influenced by D, N and D × N interactions (Tables 3 and 4).

Grain yield and its components
Grain yield (GY), spike number per m 2 (SN), grain number (GN), and 1,000-grain weight (TGW) were found to be significantly influenced by planting density and N application rate (Table 5). Although increasing plant density significantly enhanced GY and SN, both GN and TGW were significantly decreased at higher plant densities. As N application rate was increased from 120 to 240 kg ha −1 , the 2-year average GY was significantly enhanced by 43.5%, 34.4%, and 35.2% under planting densities D1, D2, and D3, respectively. Compared with the D2N1 treatment, SN was significantly increased, by 5% and 12% in treatments D2N2 and D2N3, respectively. No significant differences in GY were observed between the N2 and N3 treatments across densities D1 and D2. However, compared with D3N2, the D3N3 treatment significantly decreased GY.
In addition, we found that that SN was only significantly positively correlated with GY (Table S1), whereas it was significantly negatively correlated with TGW and GN. Furthermore, we identified positive direct path coefficients for SN and TGW to GY, respectively, whereas the indirect path coefficients were negative for SN and TGW to GY (Tables S2 and S3).   Notes. D, plant density; N, nitrogen application rate; GY, grain yield; SN, spike number per m 2 ; TGW, 100-grain weight; ETR, the effective tiller rate; PFP, partial factor productivity of applied N ; WP, the wheat population at the jointing stage. Data represent mean square. ** and * represent significance at the 0.01 and 0.05 probability level, respectively.

Changes in the wheat population and effective tiller rate
We found that the wheat populations rapidly increased from the three-leaf stage to the jointing stage, reaching a maximum before subsequently decreasing. Wheat populations increased significantly with increasing planting density and N application rate. At the jointing stage, compared with D1N1 treatment, the populations in treatments D2N1 and D3N1 increased by 125% and 141%, respectively. Similarly, compared with the D2N1 treatment, the populations obtained in the D2N2 and D2N3 treatments increased by 26.7% and 45.4%, respectively (Fig. 1). Effective tiller rate (ETR) tended to be higher than ineffective tiller rate (ITR) in the low plant density treatments (D1) (Fig. 2). ITR increased significantly with increases in plant density, with that in treatments D2N1 and D3N1 being more than doubling in both growing seasons compared with D1N1. Although application of N significantly increased ETR (P < 0.05), with that of D2N2 increasing by 24.1% compared with D2N1, the increase was less pronounced (20.8%) at the highest N application level (N3), and in both growing seasons, ETR in the D2N3 treatment was lower than that in the D2N2 treatment.

Grain number and yield per spike of different tillers
Grain number and yield per spike appeared to vary among different tiller orders (Fig. 3). The main stem had the highest GN and yield per spike. With the sequence of tiller emergence, GN and yield per spike of tillers gradually declined. For all tillers, increasing plant density resulted in a decrease in GN and yield per spike. For example, compared with treatment D1N2, yield per spike of I order tillers for D2N2 and D3N2 decreased by 29.5% and 56.5%, respectively. An increase in N application resulted in an increase in GN and yield per spike. In comparison with treatment D1N1, the yield per spike of order 0, I, and II tillers in Notes. D, plant density; N, nitrogen application rate; Gmax, the maximum grain-filling rate; Gmean, the mean grain-filling rate; T, the active grain-filling period; NBVB, the number of big vascular bundles; NSVB, the number of small vascular bundle; NTVB, the number of total vascular bundle; ABVB, the area of big vascular bundle; ASVB, the area of small vascular bundle; S, superior tiller group; I, inferior tiller group. Data represent mean square. ** and * represent significance at the 0.01 and 0.05 probability level, respectively.   response to treatment D1N2 was 7.8%, 13.0%, and 6.9%, respectively (Fig. 3A). Compared with the D1N2 treatment, the GN per spike of order I tillers for D2N2 and D3N2 decreased by 24.7% and 57.9%, respectively, whereas compared with the D1N1 treatment, the GN of order 0, I, and II tillers increased by 5.2%, 15.9%, and 8.3%, respectively, when N was applied at 240 kg N ha −1 (Fig. 3B). In addition, with an increase in plant density, there was a decrease in the number of fertile tillers. Under higher plant densities, only order I and II tillers were ear-bearing, and in the D3 treatment, only one order I tiller produced ears.

Grain-filling process and microstructure of the first internode of superior and inferior tiller groups
The main stem and tillers in density D1 plants were categorized into two groups, and defined as the superior tiller group, including order 0, I, II, and III tillers, and the inferior tiller group, including order IV, I1, and V tillers, whereas the superior tiller group for the density D2 plants included order 0 and I tillers and the inferior tiller group included order II tillers (Fig. 4). Plants at a density of D3 had just a single tiller; however, dissimilarities between GN and yield per spike meant that a distinction could be made between the main stem, as the superior group, and tiller I as the inferior tiller group. The grain-filling processes and the microstructure of the first internode of the two tiller groups were then investigated to identify the regulatory mechanisms whereby plant density and N application influence yield formation. The maximum grain-filling rate (Gmax) and the mean grain-filling rate (Gmean) of superior and inferior tiller groups decreased with an increase in plant density, but were enhanced by increasing N application (Table 6). For example, compared with treatment D2N1, the Gmax of the superior tiller group increased by 17.6% and 15.7% for treatments D2N2 and D2N3, respectively, whereas compared with treatment D1N1, the Gmean of the inferior tiller group increased by 15.0% and 19.2% for D1N2 and D1N3, respectively. Furthermore, we found that the active grain-filling period (T) of the inferior tiller group increased with increasing N application.
On the basis of the results of cluster analysis, we examined the microstructures of the first internode in superior and inferior tiller groups. Our observations revealed that the number of large vascular bundles (NBVB), small vascular bundles (NSVB), and total number of vascular bundles (NTVB), as well as the area of the large vascular bundles (ABVB) and small vascular bundles (ASVB) of the superior tiller group, were all higher than those of the inferior tiller group (Table 6). An increase in plant density significantly decreased NBVB, NSVB, NTVB, ABVB, and ASVB, whereas an increase in N application significantly increased NBVB, NSVB, NTVB, ABVB, and ASVB. For example, compared with treatment

Notes.
Gmax, the maximum grain-filling rate; Gmean, the mean grain-filling rate; T, the active grain-filling period; NBVB, the number of big vascular bundles; NSVB, the number of small vascular bundle; NTVB, the number of total vascular bundle; ABVB, the area of big vascular bundle; ASVB, the area of small vascular bundle. Values followed by different letters within the columns are significantly different at the 0.05 probability level.
D2N1, the NTVB of the superior and inferior tiller groups in the D2N2 treatment increased by 6.5 and 6.4 per culm, respectively. Compared with the D2N1 treatment, the ABVB of the superior tiller group increased significantly by 99% and 98% in the D2N2 and D2N3 treatments, and that of the inferior tiller group increased by 46.6% and 62.8%, respectively.

Relationships between individual yield per spike and vascular bundles
We observed significant correlations between individual yield and vascular bundles (Fig. 5). The grain number and yield per spike of the superior tiller group were significantly positively correlated with the NTVB (r = 0.86 * * , P <0.01; r = 0.79*, P <0.05, respectively), as were the grain number and yield per spike of the inferior tiller group (r = 0.77 *, P <0.05; r = 0.87 **, P <0.01, respectively). Significant correlations were likewise observed between vascular bundle areas and grain-filling rates (Table 7).

Trend surface analysis and evaluation of nitrogen-use efficiency
The effects of plant density and N application on GY are shown by the coefficients of second order polynomials (Fig. 6). These coefficients were calculated by the least square method, as given in the equations (Table 8). For example, the total determination coefficient R 2 = 0.992 implies that variations of 99.2% for GY are attributable to plant density and N application rate. If the partial derivative of the equation is zero, two equations can be

Notes.
Gmax, the maximum grain-filling rate; Gmean, the mean grain-filling rate; ABVB, the area of big vascular bundle; ASVB, the area of small vascular bundle. Correlation coefficients (r) are calculated and asterisks (**) represent significance at the 0.01 probability level.
constructed as follows: 13.9 − 0.012y − 0.032x = 0 and 66 − 0.012x − 0.214y = 0, and 11.3 − 0.011y − 0.024x = 0 and 68 − 0.011x − 0.222y = 0, respectively. By solving the equation groups, the following results can be obtained: x = 326 (×10 4 plant ha −1 ), y = 290 (kg N ha −1 ) in the first growing season, and x = 338 (×10 4 plant ha −1 ), y = 289 (kg N ha −1 ) in the second growing season, respectively. Similarly, the average results in two years can be obtained: x = 336 (×10 4 plant ha −1 ), y = 290 (kg N ha −1 ). Via inputting the data of these two variables into the equation, the maximum GY and its partial factor productivity of applied N could be calculated: Z MaxGY = 11,524 (kg ha −1 ) and Max GY PFP = 39.7 (kg kg −1 ). The results showed that an increase in plant density (D1 to D2) significantly increased PFP, but there was no significant difference in PFP between densities D2 and D3 (Fig. 7). In contrast, an increase in N application significantly decreased PFP. Compared with treatment D2N3, the average PFP N in the 2 years of the study was increased by 50.4% in response to treatment D2N2. Furthermore, we found that PFP was significantly negatively correlated with the rate of N application (Fig. 8A). Taking GY and NUE into consideration, treatments could be classified into four groups: (I) lower GY with higher NUE (D1N1, D2N1, and D3N1), (II) high yield with high NUE (D1N2), (III) higher yield with lower  Notes. R 2 , the coefficient of determination; z, represents the grain yield (kg ha −1 ); x, represents plant density (plant ha −1 ); y, represents the N application rate (kg N ha −1 ); PFP, partial factor productivity of applied N (kg kg −1 ). ** represents significance at the 0.01 probability level.
NUE (D1N3, D2N3, and D3N3), and (IV) higher yield with high NUE (D2N2 and D3N2) (Fig. 8B). These results indicate that increasing plant density and decreasing N application could be a useful strategy for achieving high GY and NUE.

Effects of plant density and nitrogen application on grain yield and its components
Yield variability is largely controlled by plant density and N fertilizer input (Sukumaran et al., 2015;Liu et al., 2016). In the current study, we found that plant density, N application rate, and their interactions significantly affected GY. Grain yield increased with an increase in plant density and N application, which can be explained in terms of an increase in SN, whereas, both GN and TGW were reduced by increasing plant density. Previously, Li et al. (2016b) observed that SN was significantly increased by an increase in plant density, resulting in higher grain yield, whereas TGW and GN were decreased, which is partially consistent with our results. Large changes in yield can only be determined by variation in the grain number per m 2 , which is primarily associated with SN (Slafer, Savin & Sadras, 2014).
In the present study, path analysis revealed that SN appears to be the most important factor The relationship between PFP and nitrogen application rate and grain yield. PFP represents the partial factor productivity of applied N . Max GY N represents the nitrogen application rate, which achieve the maximum grain yield calculated by the second order polynomials trend surface. Max GY PFP represents the partial factor productivity of applied N , which achieve the maximum grain yield. YG represents yield gap between the maximum grain yield and D2N2 treatment. EG represents the PFP gap between the maximum grain yield and D2N2 treatment. Full-size DOI: 10.7717/peerj.6484/ fig-8 determining GY, thereby indicating that preferentially optimizing the SN is an important measure for increasing grain yield. Spike number includes the number main stem and tiller spikes and is determined by seeding rate and tiller generation and survival (Lloveras et al., 2004;Dreccer et al., 2013).
Tillers arise from tiller buds initiated from the axillary meristems in the axils of leaves on the main shoot and tiller shoots, and the development of wheat tillers is regulated by numerous factors. Regarding genetic characteristics, multi-tiller or inhibition genes have been shown to control the development of tillers (Xue et al., 2013;Hendriks et al., 2016). Environmental factors and agronomic management, such as planting density and N application, also affect the initiation and cessation of tillering (Evers et al., 2011;Huang et al., 2013). The results of the present study indicate that an increase in plant density increases the effective tiller rate and decreases the number of high-position tiller spikes. In contrast with the low plant density (D1) treatments, we found that only the main stem and order I and II tillers survive to constitute the final population under density D2 (300 plants m −2 ).
Plant density-regulated tiller development may be related mainly to the variation in light within the canopy structure. In previous studies, light interception has been shown to increase with an increase in plant density, leading to less photosynthetically active radiation penetrating to the lower vegetation layer (Mao et al., 2014;Xue et al., 2015). Given that the development of tillers is asynchronous, the low-position tillers develop earlier than the high-position tillers. The former occupy the upper space and preempt the uppermost light source, shading the late-emerging tillers (Wang et al., 2016a). The cessation of tillering is induced when the fraction of light intercepted by the canopy exceeds a threshold (Evers et al., 2006). Plant density affects not only light quantity but also light quality, and the ratio of red to far-red light at the base of the canopy is reduced by higher plant densities (Rondanini et al., 2017), resulting in the inhibition of bud growth (Holalu & Finlayson, 2017). In addition, leaf area and leaf number both decrease as the tiller position shifts from low to high (Lafarge, Broad & Hammer, 2002). These factors may contribute to decreases in the photosynthetic activity of leaves on the late emerging tillers, due to the lack of photosynthetically active radiation and photosynthetic area, thereby resulting in lower sucrose levels, which would suppress tiller bud outgrowth (Mason et al., 2014;Kebrom & Mullet, 2015).
Nitrogen is one of the essential nutrient elements for crop growth and yield formation. In the present study, we found that the number of stems per m 2 was significantly increased by an increase in N application, and the effective tiller rate was also increased. Therefore, the final SN was increased by a high N supply (N2 and N3), compared with low N application (N1). N deficiency inhibits tiller bud elongation (Luo et al., 2017), whereas N application has been shown to promote tillering via the accumulation of microRNA393 (miR393) in response to upregulation of OsmiR393 (Li et al., 2016c). In addition to genetic factors, an interaction between N and certain hormones is known to play a role in tiller development (Domagalska & Leyser, 2011;Xu et al., 2015b). Nitrogen application promotes cytokinin biosynthesis, and inhibits its degradation, thereby inducing tiller bud development (Liu et al., 2011), and also affects auxin transport and strigolactone synthesis to regulate axillary bud activation (Jong et al., 2014). However, supra-optimal N application may result in the vigorous vegetative growth (Nasim et al., 2016;Ata-Ul-Karim et al., 2017). In the present study, we found that our highest N application rate (N3) decreased the effective tiller rate compared with N2 treatments, and, notably, a high N supply combined with a high plant density (D3N3) resulted in lodging (Fig. S4).
Plant density and N application have significant effects on culm development. Vascular bundle number (NVB) and area (AVB) are the major anatomical features influencing stem-breaking strength and transportation (Tian et al., 2015). We found that whereas increasing plant density significantly decreased NVB and AVB at the first internode, an increase in N application increased NVB and AVB, and GN and yield per spike also increased with excessive N application, Moreover, the results of correlation analysis indicated that the number vascular bundle was significantly positively correlated with yield per spike, and that vascular bundle area was significantly positively correlated with grain filling rate. These results indicate that an appropriate combination of increased planting density and N application could regulate culm quality to facilitate the transport of photoassimilates to the grain sink, resulting in an acceleration of the grain-filling rate and enhancement of yield per spike.
Optimized plant density and nitrogen application to enhance grain yield and nitrogen-use efficiency Nutrient competition, mainly with respect to N, occurs between tillers (Assuero et al., 2012), and is intensified by increased plant density. Therefore, it should be possible to optimize plant density and N application to regulate tiller development and spike number, thereby achieving high yield and high NUE via a trade-off between the three yield components. Cluster analysis showed that the low-and high-position tiller spikes can be classified into superior and inferior tiller groups, respectively. The yield per spike, grain-filling rate and NVB and AVB in the superior tiller group were all higher than those in the inferior tiller group. In the present study, we found that increasing plant density reduced the number of high-position tillers, and that under a high plant density (300 × 10 4 plant ha −1 ) the spikes consisted of the superior group. Increasing plant density increased the partial factor productivity of applied N (PFP), because the density of roots in soil was increased by increasing plant density, which enhances N uptake (Dai et al., 2014). Although increasing N application increased GY, the PFP was significantly decreased. We found that PFP was significantly negatively correlated with the N application rate, which indicates that increasing NUE should properly reduce N application. Second order polynomial trend surface analysis revealed that a plant density of 336 ×10 4 plants ha −1 and an N application rate of 290 kg N ha −1 could produce the maximum GY (Max GY ). Although the yield gap between Max GY and that obtained with treatment D2N2 was 453 kg ha −1 , the corresponding PFP (39.7 kg kg −1 ) was obviously high than that for N application at 360 kg N ha −1 .

CONCLUSION
Wheat spikes were classified into superior (including the main stem and low-position tillers) and inferior (high-position tillers) groups. The superior groups had higher grain number and yield per spike owing to a larger number of vascular bundles and faster grain-filling rate. We proposed and verified a technical approach to achieve a higher grain yield and high plant nitrogen-use efficiency through optimizing plant density and N application to regulate tiller growth, which appropriately increasing plant densities (from 75 to 336 plants m −2 ) and reducing nitrogen application (from 360 to 290 kg N ha −1 ), the