Simultaneous gains in grain yield and nitrogen efficiency over 70 years of maize genetic improvement

The competing demands of increasing grain yields to feed a growing population and decreasing nitrogen (N) fertilizer use and loss to the environment poses a grand challenge to farmers and society, and necessitates achieving improved N use efficiency (NUE) in cereal crops. Although selection for increased yield in maize has improved NUE over time, the present understanding of the physiological determinants of NUE and its key components hampers the design of more effective breeding strategies conducive to accelerating genetic gain for this trait. We show that maize NUE gains have been supported by more efficient allocation of N among plant organs during the grain filling period. Comparing seven maize hybrids commercialized between 1946 and 2015 from a single seed company in multiple N fertilizer treatments, we demonstrate that modern hybrids produced more grain per unit of accumulated N by more efficiently remobilizing N stored in stems than in leaves to support kernel growth. Increases in N fertilizer recovery and N harvest index at maturity were mirrored by a steady decrease in stem N allocation in this era study. These insights can inform future breeding strategies for continued NUE gains through improved conversion efficiency of accumulated plant N into grain yield.

A physiological framework for the analysis of NUE (Fig. 1A) considers two components: how much N is taken up by the crop (N recovery efficiency, NRE) and how efficiently N is transformed into grain yield once it has been accumulated in the crop (N internal efficiency, NIE). NIE is further broken down into the ratio between: • N content in the grain expressed as percent of the grain mass (grain N concentration) and the nitrogen harvest index (NHI), which is the fraction of the total plant N at maturity that is present in the grain 22 • the ratio between grain mass and total plant mass (harvest index, HI) and the N present in the plant at maturity expressed as percent of total dry matter (total N concentration) 23 Nitrogen internal efficiency is genetically controlled and synthesis analyses across diverse environments, managements, and genotypes, have estimated that about 56 kg of grain are produced for every kg of above-ground plant N uptake 23,24 . Nitrogen recovery efficiency is typically below 50%, meaning that much of N fertilizer applied is not used by the recipient crop. Nitrogen fertilizers not taken up by the crop, or retained in soil, are lost to the environment 4,6,23 . Improvements in NRE through agronomic management has been a research focus for several decades and has resulted in best management practices for N fertilizer including optimized fertilizer source and N application timing 22,24 . In contrast, research seeking to identify opportunities for genetic improvement in NRE and NIE is only recent 11,12,25 . Understanding how NUE and its components have changed after decades of maize breeding provides a first assessment into the opportunities for intentional breeding-driven gains in NUE. Several studies have found that NUE has increased with hybrid improvement [11][12][13]26 . However, what has not been adequately demonstrated is the source of hybrid improvement in NUE.
As pointed out in reviews by Egli (2015) 17 and Duvick (2005a) 14 , N fertilizer rates for maize production in the United States have plateaued since the 1980's. However, over this same period, grain yields have continued to increase. In order to best evaluate changes in NUE, a 0 N control is needed to account for soil N accumulation and response to N fertilizer. Two experiments that have included 0 N treatments, Haegele et al. 11 and Woli et al. 26 , found contrasting results for NUE change over time with the former noting an increase of 0.16 kg kg N −1 year −1 while the latter found no trend in NUE. The lack of difference in NUE with YOC reported by Woli et al. 26 may have been partially due to the different plant populations used for each hybrid. In the absence of a 0 N control, partial factor productivity can serve as a substitute for understanding the ability of hybrids to produce grain yield per unit of N applied. Chen and Vyn (2017) 12 and DeBruin et al. 13 both reported a positive relationship between partial factor productivity and YOC. Furthermore, the role of NIE has had mixed results showing both improvement with YOC 11 and no response. To address the knowledge gap of source of hybrid improvement in NUE over time, the objective of this study was to evaluate the mechanisms of NUE gains, through the assement of both NIE and NRE, over the past 70 years of commercial hybrid selection.

Methods
Experimental design and site description. The methods of this experiment are discussed in-depth in Mueller and Vyn (2018) 27 . In brief, a two-year experiment was carried out at the Purdue Agriculture Center for Research and Education in West Lafayette, Indiana (40.471, −86.992) in 2016 and 2017. This site consists of silty clay loam soil (fine-sily, mixed, mesic Typic Haplaquolls). The site was rainfed and managed in a maize-soybean rotation. A split-plot design was used with N treatment as the main plot and hybrid as the sub-plot with three replications. Nitrogen treatments included a 0 N control and four treatments which all received a total of 220 kg N ha −1 but varied in timing of N application. Non-0N treatments included the following where the number before the underscore indicates the rate applied at V4 and the number after the underscore represents the rate applied at R1: 220_0, 55_165, 165_55, and 0_220. These N treatments were intended to create a range of N stress at R1. Seven hybrids, representing a subset of the Pioneer ERA hybrids 28 and a more recent hybrid were used in this experiment. The hybrids ranged in year of commercialization (YOC) from 1946-2015. Hybrid names and YOC are listed in Table 1. These ERA hybrids were recommended for this study by Pioneer scientists because of their popularity and wide use during their commercial lifespan. Hybrids ranged in crop relative maturity from 111 to 114 days and the established plant population was 78,500 plants ha −1 . Although plant populations changed during the decades represented by these hybrids, previous era studies have also utilized a single "modern-era" plant population 11,19,20 . Several authors have attempted to separate the confounding effect of increasing plant density relative to genetic contributions to grain yield with conflicting results. While some researchers comparing hybrids at multiple plant densities do find that modern hybrids out-perform older hybrids at higher densities 13,29-31 , others find no YOC by plant density interaction [32][33][34] . All hybrids in this experiment were similar in silking date and the mean anthesis-silking interval did not exceed 2.1 days for any hybrid ( Table 1).

Plant measurements and nitrogen use efficiency.
Whole-plant biomass samples were collected at flowering (R1) and at physiological maturity (R6). At each of these sampling times, 10 consecutive plants were removed from the field and partitioned into stems (including tassels), leaves (including husks), and ears. At R6 www.nature.com/scientificreports www.nature.com/scientificreports/ ears were further divided into cobs and kernels. All biomass sampling zones were pre-selected shortly after seedling emergence. All biomass samples were dried to a constant weight, ground to 1 mm, and analyzed for N concentration using the combustion method by Pioneer (Johnston, IA). Grain yield, kernel number per ear, and kernel weight were determined from the 10 plants harvested at R6.
Remobilization of N content and dry matter from both leaves and stems were calculated as the differences between R1 and R6. We use remobilization to refer only to the N changes in the vegetative organs (i.e. leaves and stems).
Nitrogen use efficiency (NUE), N internal efficiency (NIE), and N recovery efficiency (NRE) were calculated as: Where GY fert and GY 0N refer to grain dry matter (0% moisture) under the fertilized and unfertilized treatments, respectively, and R6Nc fert and R6Nc 0N refer to the total R6 N content (kg N ha −1 ) in the fertilized and unfertilized treatments, respectively. Statistical analysis. Statistical analysis was conducted using SAS 9.3 35 . For the majority of measured variables, F-tests based on the mean square between the two years (2016 and 2017) had Pr (F > F 0 ) > 0.01 36 indicating homogeneity of variance between years. For this reason, the two years were pooled together and the means presented are the average of both years. Analysis of variance was conducted using PROC MIXED and means separation was determined using the LSMEANS statement in SAS. The N treatment (whole plot) and hybrid (sub-plot) were considered fixed effects. Block was treated as nested within year (block(year)) and year, block(year) and year × block(year) × N treatment were considered random effects. Standard error calculations were derived from three replications, five N treatments, and two years.
Rate of gain was calculated using the lm function in RStudio 37 with YOC as the explanatory variable and mean (including all five N treatments) grain yield, NUE, NRE, NIE, grain N concentration, NHI, HI, stover N content, stover dry weight, total N concentration, total N content, and total dry weight as the response variables 28 .
All data are represented as the mean of all N treatements, including the 0N control plots.The interaction of YOC × N treatment was rarely significant (p < 0.05) and, therefore, discussion will be based on the average of all N treatments unless otherwise specified.

Results and Discussion
Grain yield improved linearly with hybrid improvement. Over the 70 years of genetic improvement in the U.S. Corn Belt represented by the seven hybrids evaluated in this experiment, grain yield increased at a rate 0.10 Mg ha −1 year −1 , resulting in an 89% increase from 1946 to 2015 when averaged across all N fertilizer treatments (Fig. 1B, Table 2) There was no evidence that modern hybrids responded differently to the non-0 N treatments and, in all cases, grain yield increased linearly with YOC. The rate of yield gain was much greater when 220 kg N ha −1 was applied (0.12 Mg ha −1 year −1 , average of four treatments) compared to when there was no fertilizer added (0.05 Mg ha −1 year −1 ) (data not shown).
The increase in grain yield (averaged over all N treatments) was driven by a combination of both higher kernel number per ear and kernel weight ( Table 1). The increase in kernel number primarily occurred between 1946 and 1976 while the gains in kernel weight were most notable from 1976 until 2015.
The rate of yield gain in this experiment was similar to previous era studies. The average increase in the non-0N treatments of 0.12 Mg ha −1 year −1 is similar to the gains of 0.14 Mg ha −1 year −1 reported by DeBruin et al. 13 and 0.09 Mg ha −1 year −1 found by Haegele et al. 11 under similar plant densities (79,000 and 80,000 plants ha −1 , respectively) and non-limiting N rates. Our finding of a yield increase of 0.05 Mg ha −1 year −1 when no N was applied was also similar to the 0 N 11 or low N 13 gains reported by those authors of 0.06 and 0.06 Mg ha −1 year −1 , respectively.

Nitrogen use efficiency increased over time.
In parallel with grain yield, NUE increased at a rate of 0.21 kg ha −1 year −1 ( Table 2), resulting in a overall NRE increase of 73% (Fig. 1C) across the evaluated hybrids in this study. This increase was due to improvements in both NRE and NIE (Fig. 1D, E). NRE made rapid gains of 0.26 percent year −1 ( Table 2) and this was mirrored by the strong increase in total plant N uptake at R6, which increased from 160 kg ha −1 in 1946 to 192 kg ha −1 in 2015 (Fig. 2E).
In addition to the rapid gain in NRE, modern hybrids produced more grain yield per unit of N uptake (NIE, Fig. 1E). Nitrogen internal efficiency increased nearly 50% between 1946 and 2015 (average of all N treatments), improving at a rate of 0.28 kg kg −1 year −1 . This increase in NIE occurred in conjunction with a reduction in grain N concentration (numerator in NIE calculation, Fig. 1A), from 1.34 (1958) to 1.07% (2015) (Fig. 1F).
The reduction in grain protein concentration has often been implicated as having a role in improved grain yields 28,38 . However, the simultaneous increase in NHI (Fig. 1G) shows grain yields increased faster than the decrease in grain N concentration. Harvest index also increased with YOC (Fig. 2A). The change in both of www.nature.com/scientificreports www.nature.com/scientificreports/ these indices (NHI and HI) was driven entirely by increased grain N and grain dry matter accumulation because there was no change in stover (stems and leaves) N content (p = 0.221) or dry matter (p = 0.956) over the same period ( Table 2, Fig. 2B, C). In contrast, but consistent with prior studies 13,24,25 , total N concentration at maturity decreased with YOC from 0.90 to 0.83% (p = 0.042, Fig. 2D) indicating that the increase in total N content was less than the increase in total dry matter (Fig. 2E, F). These changes suggest that improved maize grain yield has not simply been due to greater total biomass accumulation, but rather to a shift in the proportional allocation of accumulated N within the plant to the grain and away from the stover in modern maize hybrids.

Leaf versus stem contributions to nitrogen efficiency changes over time.
Because of the different biological function of leaves and stems for yield determination, an observed increase in N allocation to the grain could be explained by changes in N allocation between these two plant organs in response to selection for yield. While leaf N determines rates of photosynthesis that support both vegetative and reproductive growth 38 , stem N is primarily structural and serves as transient storage for N accumulated in excess of that required to support plant growth 39 . The dynamic relationship between the stems and leaves during the growing season can be visualized through the stem to leaf ratio (S:L) of dry matter and N content (Fig. 3).
While experimental results show only small changes in stover dry matter or N content over the 70 years of hybrid improvement (Fig. 2B, C), major shifts occurred in S:L ratios for both dry matter and N content (Fig. 3). At flowering, stem dry matter was nearly twice that of leaf dry matter, but there was little pattern with YOC (Fig. 3A). This is in contrast with Chen et al. (2015a) 40 who found higher S:L DM at R1 in Dekalb hybrids released in 1967 compared to 2005 (2.5 and 2.3, respectively). In the present experiment, by R6, S:L DM generally decreased with YOC; however, the change with hybrid improvement was much starker for S:L N content. At flowering, S:L N content was greater in hybrids released prior to 1976 (Fig. 3B). However, these differences were small compared to the prominent hybrid gradient at maturity when 1946 realized a S:L N content double that of 2015 (0.8 compared to 0.4, Fig. 3B). This demonstrates the change among hybrids in the apparent remobilization during grain filling of N content present at flowering in the stems versus the leaves. To provide the N needed for the growing kernels, maize depends on a combination of N remobilized from the stems and leaves, as well as continued N uptake from the soil during grain fill. Modern hybrids remobilized a significantly higher percentage of their stem N content at flowering (58.2% compared to 31.6% for 2015 and 1946, respectively, p < 0.001, Table 1). Meanwhile, there was no pattern between YOC and the percentage of leaf N present at flowering that was remobilized, as apparent leaf N remobilization ranged from 47.7% (1976) to 60.6% (1958) (p < 0.004). These data, except 2015, agree with previous research showing that most of the total remobilized N originates from the leaves of maize hybrids compared to the stems [41][42][43] .
The greater stem N remobilization in modern hybrids led to genotypic differences in N allocation among the plant organs at maturity (Table 3) and allowed for improved NIE. As previously discussed, total dry matter at maturity increased in modern hybrids due to the greater grain dry weight while stover dry matter and N content remained unchanged (Fig. 2, Table 2). However, there were large differences in N allocation among the plant organs. As shown by the NHI results (Fig. 1G), the proportional allocation of N content to the grain increased steadily with YOC at a rate of 0.16 percent year −1 ( Table 2). The increase in NHI was mirrored by a steady decrease in stem N allocation at R6 (Table 3) falling significantly from 17.2% in 1946 to 7.3% in 2015 (p < 0.001). There was little difference in cob N content among hybrids (representing 3.0-5.9% of the total N content, p < 0.001, Table 3) and the variation in leaf N allocation was also minimal (range of 17.2-20.9% of total N content, p < 0.001, Table 3). This analysis suggests that the improved NIE in modern hybrids can be attributed www.nature.com/scientificreports www.nature.com/scientificreports/  www.nature.com/scientificreports www.nature.com/scientificreports/ to more efficient remobilization of stem N content to the grain to support kernel growth while the proportional leaf N allocation remained relatively unchanged.
In conclusion, this research addressed an important aspect of the quandary for producers and society alike of simultaneous imperatives to increase grain yields and decrease N fertilizer use and environmental losses. We have shown that modern hybrids are meeting this challenge as genetic improvement in the U.S. over the past 70 years has resulted in both a 89% increase in grain yields and a 73% increase in NUE. Modern hybrids are better at accumulating a greater percent of applied N fertilizer (NRE), which limits environmentally damaging field nutrient losses while more efficiently producing grain yield per unit of accumulated N (NIE). Although previous research has suggested that maize grain yield and NUE gains over time have arisen primarily from greater total N accumulation and dilution of grain N concentration, we provide novel evidence that increased stem N remobilization, and retention of leaf N during reproductive growth, played key roles in achieving NIE gains. Historically, green leaf area development, architecture, and duration have received the focus in maize yield improvement programs, but results from this study suggest there are additional prospects for genetic improvement in NIE. Understanding the physiological underpinnings of NIE such as the dynamics and form of stem N storage and remobilization, and transporter activity during the grain filling period, are warranted to help design breeding strategies to select genotypes with increased NRE and NIE. Table 3. Proportional N allocation among plant organs at physiological maturity. Percent of total N at physiological maturity present in the grain, cob, leaf, and stem. All means are presented as the average of two years and five N treatments. Values in parentheses represent standard error. Means denoted with different letters are significantly different from each other at p < 0.05.