
==== Front
Theor Appl Genet
Theor Appl Genet
TAG. Theoretical and Applied Genetics. Theoretische Und Angewandte Genetik
0040-5752
1432-2242
Springer Berlin Heidelberg Berlin/Heidelberg

38329519
4521
10.1007/s00122-023-04521-9
Original Article
Breeding progress of nitrogen use efficiency of cereal crops, winter oilseed rape and peas in long-term variety trials
http://orcid.org/0000-0002-7431-549X
Laidig F. Friedrich.Laidig@uni-hohenheim.de

1
Feike T. 2
Lichthardt C. 3
Schierholt A. 4
Piepho H. P. 1
1 https://ror.org/00b1c9541 grid.9464.f 0000 0001 2290 1502 Institute of Crop Science, Biostatistics Unit, University of Hohenheim, Fruwirthstrasse 23, 70599 Stuttgart, Germany
2 https://ror.org/022d5qt08 grid.13946.39 0000 0001 1089 3517 Julius Kühn Institute – Federal Research Centre for Cultivated Plants, Institute for Strategies and Technology Assessment, Stahnsdorfer Damm 81, 14532 Kleinmachnow, Germany
3 https://ror.org/04f7aqa58 0000 0004 7591 3592 Bundessortenamt, Osterfelddamm 60, 30627 Hannover, Germany
4 https://ror.org/01y9bpm73 grid.7450.6 0000 0001 2364 4210 Plant Breeding Methodology, Georg-August-University Göttingen, Carl-Sprengel-Weg 1, 37075 Göttingen, Germany
Communicated by Kai Peter Voss-Fels.

8 2 2024
8 2 2024
2024
137 2 4526 6 2023
7 12 2023
© The Author(s) 2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Key message

Grain yield and NUE increased over time while nitrogen yield did not drop significantly despite reduced nitrogen input. Selection for grain and nitrogen yield is equivalent to selection for NUE.

Abstract

Breeding and registration of improved varieties with high yield, processing quality, disease resistance and nitrogen use efficiency (NUE) are of utmost importance for sustainable crop production to minimize adverse environmental impact and contribute to food security. Based on long-term variety trials of cereals, winter oilseed rape and grain peas tested across a wide range of environmental conditions in Germany, we quantified long-term breeding progress for NUE and related traits. We estimated the genotypic, environmental and genotype-by-environment interaction variation and correlation between traits and derived heritability coefficients. Nitrogen fertilizer application was considerably reduced between 1995 and 2021 in the range of 5.4% for winter wheat and 28.9% for spring wheat while for spring barley it was increased by 20.9%. Despite the apparent nitrogen reduction for most crops, grain yield (GYLD) and nitrogen accumulation in grain (NYLD) was increased or did not significantly decrease. NUE for GYLD increased significantly for all crops between 12.8% and 35.2% and for NYLD between 8% and 20.7%. We further showed that the genotypic rank of varieties for GYLD and NYLD was about equivalent to the genotypic rank of the corresponding traits of NUE, if all varieties in a trial were treated with the same nitrogen rate. Heritability of nitrogen yield was about the same as that of grain yield, suggesting that nitrogen yield should be considered as an additional criterion for variety testing to increase NUE and reduce negative environmental impact.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00122-023-04521-9.

http://dx.doi.org/10.13039/501100001659 Deutsche Forschungsgemeinschaft PI 377/20-2 PI 377/20-2 Laidig F. Piepho H. P. BMELNot available Feike T. Universität Hohenheim (3153)Open Access funding enabled and organized by Projekt DEAL.

issue-copyright-statement© Springer-Verlag GmbH Germany, part of Springer Nature 2024
==== Body
pmcIntroduction

Nitrogen is an essential nutrient in crop production and human nutrition. Prior to the industrial production of synthetic nitrogen through the Haber–Bosch process, nitrogen was the limiting factor in agro-ecosystems and for food production as a whole (Rütting et al. 2018). The great success of the Green Revolution was largely driven by the availability of mineral fertilizers, especially nitrogen (Erisman et al. 2008), and by the parallel development of high-yielding fertilizer-responsive varieties (Borlaug 1968). Those were able to use higher nitrogen fertilizer rates more efficiently for the production of higher yields. The tremendous increase of global food production in the past was possible mainly by increasing the input of nitrogen fertilizer (Stewart et al. 2005). In Western and Central Europe, synthetic fertilizer N application in kg ha−1 increased until 1990 and then steadily decreased while crop N removal in kg ha−1 increased (Cassman and Dobermann 2022; Einarsson et al. 2021). Today 50% of the nitrogen produced synthetically by the Haber–Bosch process is used for the three major cereals, i.e., maize (16%), rice (16%) and wheat (18%). Those grain crops cover the majority of human food calories and proteins consumed either directly as grain or indirectly through livestock (Cassman and Dobermann 2022).

In Germany, the reported total nitrogen input per hectare (ha) which consists of mineral and organic N fertilizer, N deposition and biological N fixation as well as N from seeds and planting material was quite stable since the 1990s, but recently shows a decreasing trend from 2015 onwards (Fig. 1). At the same time, the total output of N in the harvested plant material slightly increased between 1991 and 2021, while the N surplus budget and the area of utilized agricultural land decreased slightly. The average N input in 2017–2021 was 201 kg N ha−1, of which the N input from mineral and organic N was 174 kg N ha−1. For the same period, the average annual N output was 139  kg N ha−1 and the N surplus budget 62 kg N ha−1, which corresponds to an input–output efficiency of 69% and a surplus budget of 31% of total N input (BMEL 2022a, b), respectively. We should note that these figures are aggregated over all German regions. The highest surplus was reported in the northwestern and the southeastern German regions. Häußermann et al. (2020) reported slightly differing numbers for N budgets in Germany, i.e., that 46% of total N input stem from mineral fertilizer, 42% from organic fertilizer, 6% from biological N fixation and 6% from atmospheric deposition in 1995–2017. Unfortunately, reliable statistical data on mineral fertilizer sales are published only for Germany as a whole. From the fact that nitrogen statistics in Fig. 1 are the average over arable and grassland, it can be assumed that significantly more nitrogen was applied in arable land than reflected by the average. In addition, there are so far no publicly accessible statistics available regarding on-farm N fertilization rates for different crops.Fig. 1 Development of total nitrogen input, N from mineral and organic fertilizer, N output and N surplus in kg ha−1 for Germany 1991–2021, and area of utilized agricultural land in Mill. ha. Source https://www.bmel-statistik.de/landwirtschaft/tabellen-zur-landwirtschaft#c8262, SJT-3070400-0000.xlsx Landwirtschaftlich genutzte Fläche nach Kulturarten. https://www.bmel-statistik.de/landwirtschaft/tabellen-zur-landwirtschaft#c8273, MBT-0111260-0000 Nährstoffbilanz insgesamt von 1990 bis 2021 in kg N/ha

The nitrogen surplus lost to the air, surface water and groundwater causes serious environmental damage. Lassaletta et al. (2014) reported that on global scale more than half of the N added to cropland is lost to the environment. While high-resolution data on the nitrogen use efficiency (NUE) of crop production, i.e., the ratio of grain produced per unit of available nitrogen in the soil (Moll et al. 1982), are lacking for most countries worldwide, few studies aimed to provide estimates of global NUE. As such, Raun and Johnson (1999) reported that NUE was approximately 33% in worldwide cereal production. In comparison, Cassman and Dobermann (2022) recently reported a global NUE of 40–50%.

Improving NUE and reducing N surplus is crucial for reducing the negative impacts on water quality as well as riparian and aquatic ecosystems in Germany and globally (Cassman and Dobermann 2022). With nitrogen being the major source of greenhouse gas emissions of crop production, an increased NUE would also help to reduce nitrogenous gases and their contribution to climate change (Riedesel et al. 2022). The reduction of N surplus needs to be realized without comprising global food security, where N fertilizer is an indispensable input. According to Connor (2008), only half of the current world population could be fed without synthetic mineral fertilizer.

Despite the above-mentioned rather low NUE in crop production, NUE increased in several countries across relevant crops (e.g., Ladha et al. 2016; Lassaletta et al. 2014). Furthermore, different studies report considerable breeding progress in NUE in winter wheat (Cormier et al. 2016; Sieling and Kage 2021; Ivic et al. 2021; Hawkesford and Riche 2020), triticale (Neuweiler et al. 2022), spring barley (Bingham et al. 2012) and oilseed rape (Kessel et al. 2012; Bouchet et al. 2016; Stahl et al. 2019, 2017).

NUE is a complex trait with many contributing processes. Apart from the type, amount and timing of fertilization, it is the result of the interaction between soil, weather, genotype and management measures. In this study, we refer to the NUE definition given by Hawkesford and Riche (2020): “NUE is the yield of grain produced per unit of N available to the crop. It is expressed as kg N in grain per kg N available.” For a more detailed description of NUE and its components, see, e.g., Moll et al. (1982), Good et al. (2004) and Hawkesford and Riche (2020).

Diverse factors and their interactions influence soil available nitrogen. Yan et al. (2019) reported based on published experiments with 15N-labeled fertilizer that most of the N in small grain crops (63%) came from sources other than the current year’s fertilizer. Part of the fertilizer N applied is assimilated into the soil. In cereals, Ladha et al. (2016) estimated that between 10% and 40% of applied N is fixed in the soil through microbial biomass and crop residues during the season of application. In a near steady-state situation, this quantity is approximately balanced by the N released from soil organic matter through mineralization. Accordingly, plant available N not only depends on the applied nitrogen fertilizer rate but also on the available soil mineral N, which is affected by the N in crop residues, atmospheric N deposition, soil quality, moisture and temperature, as well as microbial activity (e.g., Capriel 2014; Hawkesford 2014; Pituello et al. 2016; Cormier et al. 2016).

To reduce nitrogen surplus, a number of policies have been launched at the European Union and national level over the last decades to monitor and minimize N pollution. The EU Farm-to-Fork-Strategy and national agricultural policy have the goal to reduce nitrogen surplus by at least 50% and nitrogen fertilizer use by at least 20% by 2030 (BMEL 2019; EU 2020). At the same time, global demand for food and non-food agricultural products is increasing continuously. The current world population of 8 billion people is expected to grow to nearly 10 billion by 2050 resulting in an increased demand for food and respective carbohydrates and protein to be produced on limited global cropland resources.

Against the background of the required reduction of nitrogen surplus and respective nitrogen fertilizer rates, and the need to feed a growing world population, crop varieties with improved NUE are needed in combination with adapted crop management (e.g., Stahl et al. 2019; de Oliveira Silva et al. 2020). Several studies reported breeding progress in terms of yield (e.g., Laidig et al. 2014; Mackay et al. 2011; Voss-Fels et al. 2019) and disease resistance (e.g., Laidig et al. 2021, 2022; Zetzsche et al. 2020), while Cassman and Dobermann (2022) question whether new crop varieties with an apparent improvement in specific traits for NUE have actually emerged over the last decades. They attribute the improved NUE in high-fertilizer-use regions largely to a more judicious use of N fertilizers due to policies and regulations to reduce N use, rather than to the benefits of increasing crop yields.

In our study, we evaluate crops, whose protein concentrations have been routinely assessed in trials, i.e., for four cereal crops, winter oil seed rape and grain peas. Hence, the overall goals of this study are i) to quantify the breeding progress for NUE and related traits, ii) to investigate the question, whether high grain yield is closely linked to high grain nitrogen, (iii) to quantify the relationship between the varieties’ grain and nitrogen yield and the NUE of the grain and nitrogen yield, and (iv) to evaluate the potential for a direct selection for high grain nitrogen yield in trials. In particular, we first look at the relative frequency of pre-crops, the distribution of N rates and the soil-mineralized nitrogen at the beginning of the vegetation period and respective differences between crops. Secondly, we estimate the different crops’ long-term trends for breeding progress and determine changes in trends between 1995 and 2021. Thirdly, we compare the estimates of genotypic and environmental variation, correlation and heritability for crops and their traits.

Materials and methods

Variety trial and data

This study is based on data from official variety trials conducted by the Federal Plant Variety Office (Bundessortenamt, Hannover) for field crops at multiple locations during 1983 to 2021 to assess their value for cultivation and use. The investigated crops were winter wheat (WW), winter wheat under an organic testing regimen (WWORG), winter rye hybrid (WR Hyb) and population (WR Pop) varieties, spring wheat (SW), spring barley (SB), winter oil seed rape (WOSR) and grain peas (PEAS) (Table 1). On average, these crops accounted for 42.9% of the total arable land in Germany between 2019 and 2021. WW was the most important crop (24.7%), followed by WOSR (8.7%), WR (5.2%), SB (3.1%), PEAS (0.7%) and SW (0.5%) (BMEL 2022a). The trials were integrated in crop-specific continuous crop rotation regimen. The regular testing period for a newly applied variety was three years for WW, WWORG, WR, SW and SB, while for WOSR the regular testing period was two years between 1995 and 2010 and three years from 2011 onwards. Grain peas were tested for three years until 2006 and from then on for two years. Depending on the crop, time period and number of applied varieties, up to three parallel trial series (S1–S3) were run at each location in each year. In WW and SB, S1, S2 and S3 contain the varieties that were tested in their first, second and third year, respectively. In WR, varieties in their first and second testing year were tested together in one trial series and the varieties in their third year were tested separately. For PEAS, the varieties in the first and second and in SW and WWORG also those in the third test year were tested all together in one assortment.Table 1 Overview on number of observations in varieties and trials

Crop	Code	First year	Observations	No. of trials	No. of varieties	
Winter wheat	WW	1983	25,290	897	852	
Winter wheat organica	WWORG	2013	842	69	31	
Winter rye hybrid varieties	WR Hyb	1989	7712	569	244	
Winter rye population varieties	WR Pop	1989	2352	569	52	
Spring wheat	SW	1983	9546	640	155	
Spring barley	SB	1983	18,619	836	738	
Winter oil seed rape	WOSR	1995	25,655	696	797	
Grain peas	PEAS	1985	8769	533	277	
aWinter wheat under organic testing regimen

For all crops, except WWORG and PEAS, up to three different treatment intensities with different N rates were applied until 1991, while from 1992 onwards, only two intensities with different N rates were applied. From 2005 on, both intensities received identical N rates. WWORG, WOSR and PEAS were tested under only one treatment. At least three reference varieties were included in each series. The references were included in each trial series and updated on a regular basis, ensuring at least partial overlap of sets of references used in successive years.

In WOSR and PEAS, grain protein concentration (GPC) and grain oil concentration (GOC) were assessed routinely at each location and each trial. For the cereals, GPC was assessed only at a subset of six to eight locations within a trial series for which grain quality samples were taken for testing baking and malting quality of varieties. GPC was derived from N assessed by near-infrared spectroscopy using the protein equivalent factor c = 5.7 for WW, WWORG and SW, and c = 6.25 for the other crops (DIN EN 15948, 2012). As only GPC was available in the dataset, nitrogen yield (NYLD) was calculated as N in kg ha−1 accumulated in grain by NYLD = 100 × GYLD × (DM/100)  × (GPC/100)/c), where DM is the percent dry matter content in grain and c is the protein equivalent factor.

Trials for all crops were managed by the regimen of good local agronomic practice, including the application of fertilizer and growth regulators as well as the control of pests and diseases. The recorded N fertilization rates (henceforth referred to as “N rate”) for each individual trial were accumulated as total nitrogen applied in kg ha−1. If organic fertilizer was applied, the N equivalent was taken into account and added to the applied mineral N quantity accordingly. In Fig. 2a, the distribution and the average N rates (magenta color) and predicted Nmin values (green color) 1995–2021 are shown, indicating considerable differences between crops and a large variability between trials within crops. Unfortunately, no data on plant available mineralized nitrogen in the soil (Nmin) was available until 2018. Hence, we predicted these missing data by using the available Nmin data in 2019–2021 (see section “Prediction of soil-mineralized nitrogen (Nmin)”). The total available nitrogen per trial is then the sum of N rate and Nmin. We evaluated three measures of NUE: NUE for grain yield (GYLDNUE), NUE for oil yield (OYLDNUE) and NUE for nitrogen yield (NYLDNUE) in grain, expressed as kg grain, kg oil and kg N per kg available N, respectively. NUE for PEAS was not considered in this study, because this crop usually received no N fertilizer, except at a very low rate as starter in spring; therefore, we did not evaluate NUE for GYLD and NYLD.Fig. 2 Frequency distribution based on years 1995–2021 of (a) nitrogen fertilization rates (N rate) for trials and predicted soil-mineralized nitrogen (Nmin) by Eq. (9) before beginning of vegetation period in harvest year and of (b) soil fertility points (SLF). WW winter wheat, WWORG winter wheat under organic testing regimen, WR winter rye, SW spring wheat, SB spring barley, WSOR winter oil seed rape, PEAS grain peas

Trials were laid out as split-plot designs with main plots arranged in complete blocks. Trials with only one intensity were designed as complete blocks. Subplots within main plots were either laid out as randomized complete blocks or according to an alpha-lattice design. The harvested plot size was about 10 m2 on average. WR hybrid and WR population varieties were grown in the same trial and treated identically. Nevertheless, we analyzed both variety types separately. In WOSR, line and hybrid varieties were analyzed together, as only 30% of all tested varieties were line varieties and no single line variety has been registered since 2014. Compared to the large yield difference in winter rye between hybrids and population varieties, where heterosis can be utilized very effectively due to two divergent gene pools, the differences between lines and hybrids are relatively small in WOSR as shown, e.g., in Stahl et al. (2017). Accordingly, we decided that a separation is not meaningful in the context of this study.

We used only data from varieties tested for at least two years to achieve a good representation of the trial conditions. Data included in this study are shown in Table 1. The dataset was highly non-orthogonal with respect to variety × year combinations, whereas the variety × location combinations were orthogonal within year and trial series, i.e., all varieties were grown together at all locations within the same year and trial series. The data were checked for recording errors and outliers by calculating standardized residuals based on Eq. (1). We excluded observations with standardized residuals greater than ± 5.0 from further analysis.

Soil fertility and pre-crops

Variety trials were conducted in the crops’ typical growing regions across Germany. For each trial, soil fertility points (In German: Ackerzahl) were recorded. Soil fertility describes the site-specific productiveness of arable land in Germany (Weiser et al. 2018) and it is the most important factor describing the natural yield potential of a trial site (Laidig et al. 2022). Soil fertility and respective yield potential, however, is not only dependent on soil type, geological age of the parent rock and soil development stage, but also influenced by factors like climate, temperature, precipitation and topography. In a field rating, soil fertility is assessed by taking into account natural environmental conditions of a specific area of arable land (BodSchätzG 2007; Blume et al. 2015, Chap. 11.2, p. 564 ff). In Germany, soil fertility is graded on a scale from 1 to 120 points, where 1 means very poor and 120 very good soil fertility. In Fig. 2b, we show the distribution of soil fertility points for the different crops assessed in our study.

The specific pre-crops affect growth and yields of the succeeding crop. Legumes are generally more beneficial than, e.g., cereals or foliage crops. We categorized pre-crops into three groups: cereals, foliage crops (e.g., sugar beet, oil seed rape and maize) and legumes (e.g., beans, peas and clover). The percentage share of groups is shown in Fig. 3. In WW, foliage crops were the most frequent pre-crop with more than 60%, while in WOSR cereals predominated with 80%. In WWORG, legumes were grown as pre-crops in more than 80% and foliage crops in about 15% of trials.Fig. 3 Pre-crops as percent of total number of trials. Observations are based on individual trials (year × location × trial series combinations) from 1995 to 2021. Pre-crops were categorized as cereals, foliage (e.g., sugar beet, oil seed rape and maize) and legumes (e.g., beans, peas and clover). WW winter wheat, WWORG winter wheat under organic testing regimen, WR winter rye, SW spring wheat, SB spring barley, WOSR winter oil seed rape, PEAS grain peas

Statistical analysis

Basic Model

For a given observation (average over replications), we used a model with factors genotype, location, trial series and year and considering linear genetic and quadratic non-genetic long-term trends given by1 yijkl=μ+βri+γ1tk+γ2tk2+Gi+Lj+Yk+LYTjkl+GLij+GYik+ϵijkl,

where yijkl is the mean yield of the ith genotype in the jth location and kth year within the lth trial series, μ is the overall mean, β is a fixed regression coefficient for the genetic trend, ri is the first year in trial of the ith genotype, γ1 and γ2 are fixed linear and quadratic regression coefficients for the non-genetic trend, tk is the covariate for the kth calendar year, Gi is the main effect of the ith genotype, Lj is the main effect of the jth location, Yk is the main effect of the kth year, T indicates the trial series (S1, S2, S3) and (LYT)jkl is the effect of the lth trial series within the jkth location × year combination, (GL)ij is the ijth genotype × location interaction effect, (GY)ik is the ikth genotype × year interaction effect and ϵijkl is a residual comprising the genotype × location × year interaction GLYijk, the genotype × location × year × trial series interaction GLYTijkl and the error of a mean arising from sampling the replications. We confounded GLYijk and GLYTijkl with the residual error, because they were only based on the few reference varieties and were of about the same magnitude as the residual without these interactions (Hartung et al. 2023). All effects, except μ, β and γ, are assumed to be random and independent with constant variance for each effect. We modeled the genetic trend by linear and the non-genetic trend by quadratic regression terms, because inspection of graphical representation of trends indicated that genetic trends were approximately linear and non-genetic trends of quadratic shape. We estimated variance components for the random effects in Eq. (1) to get insight on the relative impact of genotypes and environmental factors on total variation by restricted maximum likelihood method (REML). We included a genetic and a non-genetic trend in Eq. (1) to avoid inflated variance components for the genotypic and year effects in case that these effects were subject to time trends. Hence, both effects can be interpreted as deviations from their respective trend functions.

Estimation of heritability

Traits assessed in registration trials should be useful for evaluating the value for the cultivation and use of a variety. Further, they should also have a reasonable predictive power to proof their performance also on farmers’ fields. Therefore, we estimated broad sense heritability for cycle means of varieties. A testing cycle includes all trials of a variety's two- or three-year testing cycle. We used equation2a H2=σG2/σG2+vd2,

where σG2 is the genotypic variance component and vd is the average variance of a difference between the means of two varieties tested in the same cycle (Piepho and Möhring 2007). To obtain vd, a dummy dataset was generated having the exact same structure in terms of the number of trials in each of the three or two years of a cycle, including the overlap of locations used in more than one year (Table 2). Dummy observations were generated for the responses of two varieties tested for three years. The dataset was analyzed using the GLIMMIX procedure of SAS, plugging on the variance component estimates of all random effects in Eq. (1) and fixing these during analysis. As the variances were fixed, the values of the dummy response were immaterial (Piepho et al. 2022). The variety main effect was taken as fixed, so adjusted means could be computed, as well as the variance of the difference between the two means (vd).Table 2 Heritability H2 of traits based a crop’s testing cycle (Eq. 2a) and heritability based on an individual trials q2 (Eq. 2b), where n corresponds to total number of trials per testing cycle, n1, n2, n3 to trials in the first second and third testing year, n12, n13, n23 to number of overlapping locations

	Number of trails per cycle	NYLD	GPC/GOC	GYLD/OYLD	NYLDNUE	GYLD/OYLDNUE	
Crop	n	n1	n2	n3	n12	n13	n23	H2	q2	H2	q2	H2	q2	H2	q2	H2	q2	
Winter wheat	24	8	8	8	0	0	4	0.83	0.22	0.96	0.65	0.90	0.41	0.82	0.21	0.90	0.39	
Winter wheat organic	22	7	7	8	5	6	6	0.61	0.14	0.95	0.68	0.83	0.37	0.56	0.12	0.77	0.29	
Winter rye Hyb	24	8	8	8	0	0	6	0.61	0.15	0.89	0.48	0.65	0.19	0.54	0.11	0.63	0.17	
Winter rye Pop	24	8	8	8	0	0	6	0.72	0.18	0.85	0.35	0.82	0.28	0.67	0.16	0.82	0.30	
Spring wheat	21	7	7	7	5	6	5	0.74	0.30	0.90	0.58	0.62	0.21	0.69	0.26	0.59	0.19	
Spring barley	23	8	8	7	0	0	3	0.57	0.08	0.78	0.20	0.72	0.21	0.56	0.07	0.72	0.16	
Winter oil seed rape	36	9	13	14	0	1	7	0.85	0.25	0.94	0.49	0.80	0.15	0.84	0.23	0.79	0.20	
Winter oil seed rapea										0.95	0.58	0.77	0.20			0.79	0.20	
Grain peas	20	10	10		7			0.88	0.27	0.96	0.57	0.88	0.28					
Mean	24.3	8.1	8.6	8.6	2.1	1.9	5.3	0.73	0.20	0.91	0.51	0.78	0.26	0.67	0.17	0.75	0.24	
All results are based on years 1995–2021. Grain peas were tested for only two years. NUE for grain peas were not listed as they received no regular nitrogen fertilizer

Organic Winter wheat under organic testing regimen with no mineral N and no pesticides, Hyb hybrid varieties, Pop population varieties, NYLD nitrogen yield, GPC/GOC grain/oil protein concentration, GYLD/OYLD grain/oil yield, NYLDNUE, GYLD/OYLDNUE nitrogen use efficiency for NYLD and GYLD/OYLD, respectively

aGrain oil concentration (GOC), Oil yield (OYLD) and NUE for oil yield (OYLDNUE) in winter oilseed rape

Heritability H2 estimated by Eq. (2a) depends on the relative magnitude of variance components but also on the number of years and locations typical for a specific crop’s testing cycle. To get a measure of heritability independent of the number of locations and years and hence comparable between crops, we estimated a trial-specific measure, given by2b q2=σG2/σG2+σGY2+σGL2+σRes2,

where σG2, σGY2,σGL2 and σRes2 are variance components for genotypes, the interaction of genotypes × year, genotypes × location and the residual error.

Model for overall trend

The overall trend was modeled by confounding year and genotypes within years, i.e., genotypes are nested within years (Laidig et al. 2014). Thus, compared with the basic model (Eq. 1), for this analysis we dropped effects involving genotypes that are not nested within years, i.e., the effects Gi and (GL)ij. Then the reduced model is given by3 yijkl=μ+Lj+Yk+LYTjkl+GYik+GLYTijkl.

Further, we assumed that Yk is subject to an overall long-term time trend, confounding genetic and non-genetic trends. We extended Eq. (3) by fixed linear and quadratic regression coefficients. Then, the model is given by4 yijkl=μ+α1tk+α2tk2+Lj+Yk+LYTjkl+GYik+GLYTijkl.

where α1 and α2 are fixed linear and quadratic regression coefficients for the overall trend, tk is the continuous covariate for the calendar year. The expected value under this model is given by5 Eyijkl=μ+α1tk+α2tk2.

Estimation of breeding progress

As WOSR data were only available from 1995 onwards, we estimated the change achieved between tk= 1995 and 2021 based on the overall trend given by Eq. (4) to ensure that estimation of breeding progress between crops and traits was based on the same time period. Hence, for all crops we calculated the change as the difference of the predicted values for year 2021 and 1995 byDiff=Eyijkltk=2021-Eyijkltk=1995=

6 α12021-1995+α220212-19952usingEq.(4).

Genotypic, environmental and G × E correlation

As in long-term trials, the strength of the association between pairs of traits can be influenced by several effects. Simple correlation coefficients over all observations, for example the Pearson sample correlation coefficient, are not always appropriate to allow valid inferences as the structure of the trial series is not considered. The simple sample correlation coefficient does not indicate which effect, genotype or environment, was dominating. To get insight, we therefore decomposed the correlation between traits by their individual random effects as given in Eq. (1). We estimated correlation coefficients (Piepho 2018) between traits based on variety × year × location × trial series observations (Eq. 1).

The correlations between random effects of Eq. (1) were calculated assuming a multivariate model with traits as independent variables. We choose a univariate approach from which correlations for pairs of traits can be inferred (Piepho et al. 2014):We calculated variance components of random effects according to the model of Eq. (1) for trait U and V and for the difference U-V between both traits.

We computed covariances between the random effects of trait U and V from variance components obtained from univariate models by using the equation7 var(U-V)=var(U)+var(V)-2cov(U,V)⟺

8 Cov(U,V)=12(varU+varV-varU-V)

We used variances of random effects from Eq. (1) and their covariance from Eq. (7) to calculate the correlation coefficients.

The marginal correlation coefficient was derived by the marginal variances and covariances which are the sum over individual random effects of var (U), var (V) and cov (U,V). Compared to the simple Pearson correlation coefficient (rP), the marginal correlation (rM) is the correlation on the level of observations (variety × year × location × trial series combinations), which takes into account the model structure of the trial series (Piepho 2018). We aggregated the variances and covariances of the random effects Lj,Yk,LYTjkl to the environmental (E) effect and GLij,GYik to the genotype × environment interaction (G × E) effect. Finally, we obtained the correlation coefficients for the marginal (rM), genotypic (rG), the G × E (rG×E), the environmental (rE) and the residual effects (rRes). The magnitude of the marginal correlation depends mostly on the magnitude of the random effect with the largest variance and covariance components.

Prediction of soil-mineralized nitrogen (Nmin)

Soil-mineralized nitrogen was assessed up to 60 cm soil depth for each trial in spring before the start of vegetation. Nmin was considered by fixing the target N rate for a given trial such that available N is the sum of applied N rate and the assessed Nmin (DUEV 2017). However, Nmin data were only available for 2019–2021. We utilized those available data (across crops, n = 259) to predict Nmin data for trials where no data were available. As the distribution of Nmin values showed a right-skewed shape, we transformed the data by a logarithmic function to achieve a more symmetric distribution. The model is given by9 yijkl=μ+(CP)im+δmajkl+ηicjkl+Lj+Yk+LYTjkl,

where yijkl is the log-transformed Nmin assessed for the ith crop at the lth trial series within the jth location and kth year, μ is the overall mean, (CP)im is a categorial effect of the ith crop and the mth pre-crop, δm the fixed pre-crop-specific regression coefficient of the linear trend for SLF, ajkl is the covariate represented by the SLF point of the jklth trial and ηi the fixed crop-specific regression coefficient of the linear trend for the N rate and cjkl is the covariate represented by the N rate in kg ha−1 of the jklth trial, Lj is the main effect of the jth location, Yk is the main effect of the kth year and LYTijkl is the residual error. We assumed that the effects Lj, Yk and (LYT)jkl are random and independent with constant variance, while all other effects are considered as fixed. In the model selection procedure we started with a basic model which was given by the random effects Lj, Yk and LYTjkl, only, and added stepwise the fixed effects as given in Eq. (9). As selection criterion we used the coefficient of determination R2 (Piepho 2019). We stopped the selection of model terms until it reached R2 = 41.5% and could not be improved further. The resulting model is given by Eq. (9). The back-transformed best linear unbiased predictors (BLUP) for Nmin of Eq. (9) were used for all trials between 1983 and 2021, assuming that the Nmin was not subject to a time trend in 1983–2021.

Results

Overall trends for breeding progress of NUE and related traits

Generally, we should note that breeding progress estimated for all traits was subject to two confounded processes, the introduction of continously improved new varieties and decreasing N rates. Overall breeding progress was estimated by a mixed linear model where year effects were assumed to follow a quadratic time trend (Eqs. 4 and 5). The difference of the estimated trends in 2021 and 1995 was considered as the breeding progress (Eq. 6). In Table 3 the estimated levels in 1995 and 2021, plus the absolute and relative differences are shown for NYLD, N rate, GYLD/OYLD, GPC/GOC, and for NYLDNUE and GYLD/OYLDNUE. In Fig. 4, trends and Nmin levels are displayed.Table 3 Performance levels of overall trends at years 1995 and 2021 and difference between levels 2021 and 1995 expressed in absolute (Diff) and relative (%) values based on level 1995 (Eq. 6). Table for (a) nitrogen yield in grain kg ha−1, N fertilizer kg ha−1, grain protein concentration % and grain/oil yield dt ha−1 and (b) for nitrogen use efficiency for nitrogen yield kg kg−1 and grain yield/oil kg kg−1. NUE for grain peas were not listed as they received no regular nitrogen fertilizer

(a)					
	Nitrogen yield kg ha−1	Nitrogen fertilizer kg ha−1	Grain/Oil protein concentration %	Grain/Oil yield dt ha−1	
Crop	1995	2021	Diff	%	Sign	1995	2021	Diff	%	Sign	1995	2021	Diff	%	Sign	1995	2021	Diff	%	Sign	
Winter wheat	187.4	194.0	6.6	3.5	ns	180.1	170.3	−9.8	−5.4	ns	13.1	12.5	−0.6	−4.4	**	95.4	103.3	7.9	8.3	**	
Winter wheat organica	100.6	94.9	−5.7	−5.7	ns	8.9	23.1	14.2	158.8	ns	11.7	11.1	−0.6	−5.3	ns	55.9	57.2	1.3	2.3	ns	
Winter rye Hyb	134.1	135.1	1.1	0.8	ns	130.4	107.8	−22.6	−17.3	***	10.5	9.3	−1.2	−11.2	***	84.9	96.6	11.6	13.7	**	
Winter rye Pop	122.0	117.9	−4.1	−3.3	ns	130.4	107.8	−22.6	−17.3	***	10.8	10.0	−0.9	−8.1	*	74.4	78.4	3.9	5.3	ns	
Spring wheat	160.9	155.0	−5.9	−3.6	ns	163.1	115.9	−47.2	−28.9	***	14.0	13.9	0.0	−0.3	ns	76.8	74.4	−2.3	−3.0	ns	
Spring barley	94.3	102.3	8.0	8.5	*	68.6	83.0	14.3	20.9	***	10.8	10.3	−0.5	−4.5	ns	64.0	72.4	8.4	13.2	**	
Winter oil seed rape	51.3	51.4	0.0	0.1	ns	176.9	143.5	−33.4	−18.9	***	18.8	17.4	−1.4	−7.6	**	44.8	48.6	3.8	8.5	*	
Winter oil seed rapeb											41.6	42.7	1.2	2.8	*	17.9	20.3	2.4	13.5	*	
Grain peas	176.0	164.6	−11.3	−6.4	ns	5.5	3.0	−2.5	−45.1	ns	20.1	20.3	0.1	0.7	ns	54.7	50.1	−4.5	−8.3	ns	
(b)			
	NUE Nitrogen yield kg kg−1	NUE Grain/Oil yield kg kg−1	
Crop	1995	2021	Diff	%	Sign	1995	2021	Diff	%	Sign	
Winter wheat	0.84	0.91	0.07	8.0	*	43.0	48.5	5.5	12.8	***	
Winter wheat organica	1.76	1.31	−0.45	−25.5	ns	96.8	79.6	−17.2	−17.8	ns	
Winter rye Hyb	0.89	1.01	0.12	13.0	*	56.7	71.9	15.2	26.8	***	
Winter rye Pop	0.81	0.88	0.07	9.1	ns	49.4	58.4	9.0	18.2	***	
Spring wheat	0.76	0.88	0.12	15.9	*	36.2	42.4	6.2	17.2	*	
Spring barley	0.76	0.83	0.06	8.0	ns	51.9	58.9	7.0	13.5	*	
Winter oil seed rape	0.25	0.30	0.05	20.7	**	20.9	27.3	6.4	30.8	***	
Winter oil seed rapeb						8.7	11.8	3.1	35.2	***	
aYears 2013–2021; bOil

Hyb hybrid varieties, Pop population varieties, Diff absolute difference between levels at years 2021 and 1995 based on overall quadratic regression estimates (Eq. 6), % relative difference as percent of yield level 1995, NUE nitrogen use efficiency expressed as nitrogen yield or grain/oil yield per kg applied N fertilizer and soil-mineralized nitrogen, Sign significance level; *Significant at 5% level; **Significant at 1% level; ***Significant at 0.1% level

Fig. 4 Adjusted overall year means (Eq. 3) (blue circles), quadratic regression lines (Eq. 5) (blue lines) based on years 1983–2021 (a) for nitrogen yield in grain (kg N ha−1), nitrogen fertilizer application rate (kg ha−1) (magenta circles and lines) and average predicted soil-mineralized nitrogen (green horizontal lines), grain/oil yield (dt ha−1), grain/oil protein concentration (%), (b) NUE for grain yield and nitrogen yield. NUE for grain peas were not displayed as they received no regular nitrogen fertilizer. WW winter wheat; WWORG winter wheat under organic testing regimen; WR winter rye, Hyb hybrid and Pop population varieties, SW spring wheat, SB spring barley, WSOR winter oil seed rape, PEAS grain peas, NYLD nitrogen yield in grain, GPC grain protein concentration, GOC grain oil concentration, GYLD grain yield, OYLD oil yield, NUE GYLD, NUE NYLD, NUE OYLD NUE for grain yield, nitrogen yield in grain, oil yield in grain

The N rate decreased in most crops, especially in WR, SW and WOSR. SB was the only crop where the N rate increased by 14.3 kg ha−1, corresponding to 20.9%. Highest N rates were applied in WW with 180.1 kg ha−1 in 1995 and 170.3 kg ha−1 in 2021 followed by WOSR with 176.9 kg ha−1 in 1995 and 143.5 kg ha−1 in 2021. WWORG and PEAS received only negligible mineral fertilizer so that their results are not directly comparable with the other crops.

In WWORG, no significant trends were estimated (Table 3) due to the large year-to-year variation relative to the short 9-year period. This is another reason for the limited comparability with the other crops which build on data from 27 years. Despite the low number of years, a rough comparison with WW indicated a slight nonsignificant increase in GYLD (2.3%) but on a GYLD and NYLD level of more than 50% below the corresponding levels in WW. This may be partially attributed to the lower recovery rate of organic N fertilizer compared to the mineral N fertilizer which was supported by Yan et al. (2019) reporting a higher recovery from mineral N (37%) than from organic N (27%). The lower recovery of organic N fertilizer may be partially compensated by the observed higher Nmin in WWORG (59.6 kg ha−1) compared to WW (45.5 kg ha−1) as about 80% of WWORG trials had legumes as pre-crops (Fig. 2).

The change for NYLD between 1995 and 2021 was in the range of −5.7% (WWORG) to 8.5% (SB), but only the increase in SB was significant, yet at a low NYLD level in 1995 (94.3 kg ha−1) and in 2021 (102.3 kg ha−1) We found the highest NYLD in WW and PEAS whereas WOSR had a rather low NYLD (about 51 kg ha−1).

In most crops, GPC decreased, especially in WR hybrid varieties (−11.2%). Noticeably, GPC in WOSR decreased (−7.6%) while GOC in WOSR increased slightly (2.8%). GYLD increased in the range of 2.3% (WWORG) and 13.7% (WR hybrid varieties). A nonsignificant decrease in GYLD was found in SW (−3.0%) and PEAS (−8.3%).

In Fig. 4b, the trends for NYLDNUE and GYLD/OYLDNUE are displayed, and in Table 3b, the corresponding changes between 1995 and 2021 are given. Results for NUE in PEAS are not shown, because they received only limited N fertilization and were therefore not comparable with other crops in NUE traits. NYLDNUE increased in all crops, but only significant at the 5%-level in WW (8.0%), WR hybrid varieties (13.0%), SW (15.9%) and at the 1%-level in WOSR (20.7%).

Genotypic, environmental and G × E variation

Varieties grown over many years and locations are exposed to a wide range of environmental conditions as shown in Figs. 2 and 3. Variance components for the random effects given in the basic model (Eq. 1) were estimated. The graphical representation of genotypic, G × E interaction, environmental and residual components and the marginal variance is shown in Fig. 5 and more details are given in Supplementary Material Table S1. First, we consider the relative magnitude of variation for individual traits across crops. Environmental variation was the dominating component. For NYLD and GYLD/OYLD it ranged roughly between 80% and 90%, only for GPC/GOC the environmental variation was smaller (60–90%). Nearly the complete variation (around 95%) for NYLDNUE and GYLD/OYLDNUE was caused by the influence of environmental conditions. In contrast to the environmental variation, the genotypic variation was small. For NYLDNUE it was about 1% on average, for GYLD/OYLDNUE about 1.5%, for NYLD about 2.5% and GYLD/OYLD about 4%. Among all traits, the genotypic variation of GPC/GOC was highest, especially for WW, WOSR and PEAS. The G × E component was of about the same magnitude as the genotypic variation. The residual error showed the second largest variation. Across all traits, WW showed the largest and WR hybrid varieties the lowest genotypic variation.Fig. 5 Variance components as percent of marginal variance (total sum of variance components for random effects given by Eq. 1) considering linear genetic trend in the genotype effects and quadratic non-genetic trends in the year effects based on 1995–2021. Traits for OYLD, GOC and OYLDNUE in WOSR are shown in (b). Y-axis was truncated at the 40% level. NUE for grain peas were not listed as they received no regular nitrogen fertilizer. WW winter wheat, WWORG winter wheat under organic treatment regimen, WR winter rye, Hyb hybrid, Pop population varieties, SW spring wheat, SB spring barley, WSOR winter oil seed rape, PEAS grain peas, NYLD nitrogen yield in grain, GPC grain protein concentration, GOC grain oil concentration, GYLD grain yield, OYLD oil yield, NYLD_NUE (GYLDNUE) NUE of nitrogen yield in grain (GYLDNUE = NYLD/available N), GYLD_NUE (GYLDNUE) NUE of grain yield (GYLDNUE = GYLD/available N), OYLD_NUE (OYLDNUE) NUE of oil yield (OYLDNUE = OYLD/available N), available N = N rate + Nmin, G genotype, G × E Genotype ×  environment interaction, Res residual, E environment

The strikingly low genotypic but high environmental variation for NUE traits compared to the other traits requires an explanation. We think this is due to the fact that NUE is derived as the ratio of two random variables which were approximately independent, e.g., GYLD and available N, which increases the variation according to the law of error propagation (average GYLD per trial and available N were nearly uncorrelated, data not shown). Importantly, GYLD varies between genotypes within one trial, whereas the amount of applied N is the same for all genotypes in the same trial and hence only varies between environments. Further, looking at this ratio in terms of model effects (Eq. 1), we find that the ranking of genotypic effects within a trial is unaffected by the constant value of available N in the denominator, while the environmental effects show an increased variance. This results in a higher percentage of environmental variance in the total variance and hence in lower percentages of the other components as Fig. 5 and Supplementary Material Table S1 show.

Despite the remarkable very large environmental variation, compared to the genotypic component, considerable breeding progress has been achieved as shown in the previous section, because each year new variation due to new varieties enter trials where this variation can be used to achieve progress.

Genotypic, environmental and G × E correlation

The overall correlation between traits based on observations of variety × year × location × trial series combinations were estimated by marginal correlation coefficients using a univariate approach (Eqs. 7 and 8). We further dissected the marginal correlation coefficient rM by the genotypic rG, genotype × environment interaction rG×E, residual rRes and environmental correlation rE coefficients to show the strength of association of genotypic, environmental and residual effects. To check the plausibility of the marginal correlation, we estimated the Pearson sample correlation coefficient rP over all observations, which was approximately of the same magnitude as the marginal one. The strength of correlation between traits was denoted by the following categories: |r|< 0.15 very weak, 0.15 ≤|r|< 0.35 weak, 0.35 ≤|r|< 0.55 moderate, 0.55 ≤|r|< 0.75 strong, 0.75 ≤|r| very strong.

The marginal correlation coefficients rM in Table 4 show that NYLD was moderate to strong correlated with GYLD/OYLD while NYLD was only weak to moderate correlated with GPC. As expected, GYLD was negative but weak correlated with GPC/GOC. The environmental correlation rE between NYLD and GYLD/OYLD showed strong to very strong while rE between NYLD and GPC was less strong in the range of weak to moderate. The genotypic correlation was positive for all crops for NYLD with GYLD/OYLD (rG = 0.37 on average) and GPC (rG = 0.45 on average) varying in a wider range than environmental correlation coefficients. This indicates that first, NYLD was stronger linked to GYLD than to GPC/GOC and second that the selection for varieties with high genotypic value for NYLD does not counteract with the selection for varieties with high genotypic value for GYLD/OYLD and GPC. In WOSR, the genotypic correlation for GYLD with GOC (rG = 0.49) was positive, contrary to the correlation with GPC (rG = −0.21). However, GPC was negative correlated with GOC (rG = −0.47).Table 4 Decomposition of marginal correlation (rM) by genotypic (rG), genotype × environment interaction (rG×E), environmental (rE) and residual (rRes) effects (Eqs. 1 and 7), where rP is the Pearson sample correlation coefficient, sign its significance level, n number ot total observations and nG number of genotypes. Correlation coefficients are based on years 1995–2021

WW winter wheat, WWORG winter wheat under organic testing regimen, WR winter rye, Hyb hybrid varieties, Pop population varieties, WOSR winter oil seed rape, SW spring wheat, SB spring barley, PEAS grain peas, NYLD nitrogen yield in grain, GPC grain protein concentration, GYLD grain yield, OYLD oil yield, GOC grain oil concentration, G genotype, Y year, L location, Y × L × T interaction of trials within Y × L, G  × E genotype × environment interaction (G × E = G × Y + G × L), E environment (E = Y + L + Y × L × T), Res residual

To evaluate the question of how strong NYLD was correlated with NYLDNUE and GYLD/OYLD with GYLD/OYLDNUE, we also dissected the marginal correlation as described above. For these traits the genotypic, G × E and residual correlations were approximately one (rG ≈ 1, rG×E ≈ 1, rE ≈ 1) while the environmental correlation rE was lower in the range of 0.33 and 0.83 (Supplementary Material Table S2). This is a general result which holds only if all varieties in an individual trial received the same N rate. The numerical results in Table S2 have been confirmed by a mathematical proof shown in Appendix. In the case of our study, the genotypic value of a variety is about equivalent to the least square estimate of a variety’s 3-year cycle mean using Eq. (1). Hence, the ranks of the estimated variety means for GYLD and GYLDNUE, and NYLD and NYLDNUE are about the same. This is demonstrated graphically by a few examples given in Supplementary Material Fig. S1. In consequence, this result indicates that the selection of varieties for high GYLD or high NYLD implies the selection of varieties for high GYLDNUE and NYLDNUE.

Heritability

Traits assessed in registration trials should not only be relevant for evaluating the value for cultivation and use but results obtained during the two- or three-year testing cycle of a variety should also be repeatable at farmers’ fields. As a measure of repeatability, we show heritability coefficients H2 (Eq. 2a) in Table 2 for individual traits based on the crop-specific testing systems. Further, we show the trial-specific heritability coefficients q2 (Eq. 2b), which are independent of the number of years and locations for a given crop’s testing cycle.

The total number of trials for a crop’s present testing cycle was in the range of 36 in WOSR and 21 in SW. In the second and third testing year, more overlapping locations were available compared to the other years. PEAS was the only crop with a testing period of 2 years. Across crops, cycle-based heritability was on average across crops highest for GPC/GOC (H2 = 0.91) followed by GYLD/OYLD (H2 = 0.78) and NYLD (H2 = 0.73) while NYLDNUE reached only H2 = 0.69. Noticeably, the average heritability for GYLD/OYLDNUE (H2 = 0.75) was higher than for NYLDNUE (H2 = 0.67). The magnitude of the average trial-based heritability coefficients was of about the same rank order as for the cycle-based coefficients in the range of 0.17 ≤ q2 ≤ 0.26. Table 2 indicates that across traits, the largest values for cycle-based heritability were achieved in WW, WOSR and PEAS, while heritability in WR hybrid varieties and SB was relatively low. Overall, the average heritability coefficient of NYLD was of about the same magnitude as for GYLD/OYLD, which is generally the most important registration criterion.

Discussion

Soil-mineralized nitrogen (Nmin)

When considering the predicted Nmin values as shown in Fig. 2a, it should be noted that the predictive power of the BLUP values was based on a moderate coefficient of determination of R2 = 41.5% and therefore predicted Nmin values only approximately. Further, the variation between trials is likely lower than the variation for the actual Nmin values as BLUP shrinks toward the mean. We assumed that no noteworthy time trends in Nmin during 1983 and 2021 were present, because registration trials were integrated in the crops typical crop rotation sequence. This is consistent with the fact that in the course of long-term experiments only small net changes in soil-mineralized N were observed (Johnston et al. 2009; Sylvester-Bradley and Kindred 2009; Ladha et al. 2016), because a near steady state may be achieved in fields under continuous crop rotation regimen. In a long-term experiment with winter wheat in the UK, van Grinsven et al. (2022) reported an available nitrogen from soil (nitrogen deposition and natural biological nitrogen fixation from free-living bacteria) between 4 and 64 kg N ha−1 with a mean of 30 kg N ha−1 showing no trend in time. In an another UK wheat study (Hawkesford and Riche 2020), soil mineral N ranged from 25.6 to 115.7 kg N ha−1 between 2006 and 2017. These results indicate a large variation of Nmin, which is in accordance with the distribution of Nmin predicted in this study as shown in Fig. 2b.

Overall trends for breeding progress of NUE and related traits

Table 3a has shown that N rates decreased considerably. Cassman and Dobermann (2022) reported a reduction of N input in Western and Central Europe during the last decades. They assumed that this was probably due to general agricultural policy and the gradual tightening of fertilizer regulation requirements to reduce nitrate leaching and environmental pollution. We assume that the observed reduction in N rates in registration trials conducted according to good local agronomic practice is also due to the same reasons. SB was the only crop where a significant increase of N rate (20.9%) was found. However, the increasing N rate in SB may be explained by the interaction of several factors. In this crop, N was applied at the rather low malting barley level to balance grain yield and malting quality criteria to fulfill the strict requirements of malting industry (Barmeier et al. 2021). For example, malting barley varieties’ GPC must be below 12% as required by the Federal Plant Variety Office (BSL 2023). We believe that increasing N rates in SB was possible due to the introduction of shorter rh-genotypes without compromising malting quality due to a lower lodging risk. Further, higher N rate increased yield without increasing GPC due to the trade-off mechanism between GYLD and GPC (de Oliveira Silva et al. 2020). This is supported by the fact that GPC did not increase, but decreased by −4.5% as shown in Table 3a.

Despite the strongly reduced N rates, GYLD/OYLD mostly increased significantly, which suggests a breeding progress through new genotypes. Numerous other studies confirmed breeding progress with regard to yield (e.g., Laidig et al. 2014; Mackay et al. 2011; Voss-Fels et al. 2019). Nevertheless, this study shows that breeding progress was even possible under reduced N rates. The reduced GPC in all crops can be attributed to two effects; first to the reduced N fertilizer rates and second to the trade-off effect between yield and protein due to the well-known negative association between both traits. The strong reduction of N fertilizer and consequently the decay in GPC were not reflected by the same magnitude in the changes for NYLD. This can be explained by the stronger link of GYLD with NYLD than of GPC with NYLD, shown by the correlation coefficients in Table 4.

The highest breeding progress for NUE was achieved in WOSR compared with other crops. Despite this progress WOSR still had the lowest NUE levels (e.g., for NYLDNUE, 0.25 kg kg−1 in 1995 and 0.30 kg kg−1 in 2021) among all crops (Table 3, Fig. 4). Our results are in line with Sylvester-Bradley and Kindred (2009) who compared NUE of major agricultural crops in the UK, indicating that harvested WOSR had the lowest NUE, lower than for cereals. WOSR depends on higher N fertilization than other crops due to its low NUE which may partially be explained by the negative correlation between oil and protein (r = −0.76). The breeding focus on high oil concentrations - as one target trait of an oil crop - could conflict with a higher mobilization of N from source to sink. In other studies, reported NUE in WOSR was often not exceeding 60% due to its low ability to remobilize plant stored N (Bouchet et al. 2016; Stahl et al. 2017, 2019).

PEAS are recovering N mostly from air by N2 fixation of symbiotic bacteria. Kelstrup et al. (1996) estimated a bacteria fixed N amount of 122 kg N ha−1 and a soil accumulated N of 57 kg ha−1 achieving a grain yield of about 43 dt ha−1 and Ruisi et al. (2012) an average N of 25 kg ha−1. We found a higher yield level for PEAS of 50 dt ha−1, which may likely be due to the breeding progress of newer varieties investigated in our study. Yang et al. (2017) reported lower GYLD of 21.9–51.5 dt ha−1 from six pea varieties in western Canada. Table 3 and Fig. 4 show that PEAS accumulated the highest amount of N in grain besides WW, however, with nearly zero N fertilizer, indicating the effect of a legume crop.

Our study showed a very strong increase in NUE for GYLD/OYLD and a lower one for NYLD. We have to look at this increase by taking two aspects into account. Firstly, the increase due to new improved varieties and secondly, the increase due to the reduction of N fertilizer. Regarding the first aspect, Lassaletta et al. (2014) reported on 50-year NUE trends (% harvested N in protein /N input onto cropland) across cropping systems in Western Europe and found considerable increase in NUE, e.g., in France from about 40% in 1980 to nearly 80% in 2010. Higher NUE of newer varieties in WW was also reported by Guarda et al. (2004), Ladha et al. (2016), Guttieri et al. (2017), Ivic et al. (2021), Sieling and Kage (2021), in WOSR by Stahl et al. (2017, 2019), and in winter triticale by Neuweiler et al. (2022). The second aspect was confirmed by studies testing varieties under different N rates, which found a general agreement that NUE increases with decreasing N rates or, conversely, NUE decreases with increasing available N (Cormier et al. 2016; Sieling and Kage 2021, 2022). However, Cassman and Dobermann (2022) question whether new crop varieties with an apparent improvement of specific traits for NUE were actually generated. They largely attribute the improved NUE in high-fertilizer-use regions to a more judicious use of N fertilizers as a result of policies and regulations aiming to reduce N use, rather than to the benefits of increasing crop yields.

Disentangling the confounded effects of breeding progress by new varieties and the increase of NUE by reduction of N rates was statistically not possible. However, we elaborate in the following that breeding progress for NUE of GYLD was at least as high as for GYLD. Let us assume a constant N level from 1995 to 2021, e.g., in WW (180.1 kg). Then NUE in 1995 is equal to the ratio of GYLD and N level in 1995 (95.4 dt ha−1/180.1 kg ha−1 ≙ 53.0 kg kg−1) and NUE in 2021 is equal to the ratio of GYLD in 2021 and N level in 1995 (103.3 dt ha−1/180.1 kg ha−1 ≙ 57.4 kg kg−1). The derived difference of NUE for GYLD in 2021 and in 1995 (4.4 kg kg−1) relative to the NUE level in 1995 (53.0 kg kg−1) then equals the relative change of GYLD in WW (4.4 kg kg−1/53.0 kg kg−1 ≙ 8.3%), both in the hypothetical example and in the observed data (see Table 3a). This calculated 8.3% correspond to the lower limit of breeding progress for NUE, as in reality N levels went down from 1995 to 2021 comprising the actual breeding progress for GYLD and especially NYLD, which would be higher if N levels were maintained at 1995-level. These considerations allow concluding that the breeding progress for NUE of NYLD and GYLD (Table 3a) is actually somewhere in between the breeding progress for NYLD and GYLD (Table 3b) if N rate would have been unchanged at level 1995. Accordingly, breeding progress for NUE was actually achieved by new varieties (under considerable N fertilizer reduction), and improved NUE was not just a result of reduced N fertilizer use over time as Cassman and Dobermann (2022) assumed. Furthermore, the absolute reduction in N rate far outweighed the small reduction of grain accumulated N, indicating considerable mitigation of adverse environmental impacts.

Genotypic, environmental and G × E variation

Our study built on trial data with a single N rate following good agronomic local practice by providing N according to crop demand. Figure 2 showed that N rates differed strongly between environments, which was most likely due to differences in trial-specific Nmin supply and differences in actual N demand driven by differences in yield potential between trials. Hence, we could not estimate variance components for genotype by N rate interaction. This raises the often-asked question whether trials with only one N level, as in this study, will be efficient enough to select genotypes with high NUE also under lower N fertilizer rates. Most studies on crops’ NUE are conducted with two or more N rates. Cormier et al. (2016) stated in a review paper that numerous studies on wheat (e.g., Ortiz-Monasterio et al. 1997; Le Gouis et al. 2000; Laperche et al. 2006; Barraclough et al. 2010; Cormier et al. 2013), detected significant genotype × N rate interactions for agronomic traits, meaning that the genotypic values of varieties differ between N levels and that selection for varieties in low N target regions may be efficient if the magnitude of genotype × N rate interaction was large compared to the genotypic variation. Also, several studies with contrasting N rates found genotype × N rate interactions for NUE and related traits, but of very low magnitude compared to the genotypic variation, as reported by Voss-Fels et al. (2019), Ivic et al. (2021) and Brasier et al. (2020) in WW, Anbessa et al. (2009) in SB, and Bouchet et al. (2016), Kessel et al. (2012) and Stahl et al. (2017, 2019) in WOSR. Accordingly, results from a UK germplasm diversity trial with recent winter wheat varieties grown under five different N fertilization levels indicated that the ranking of the varieties at each of the N rates is almost identical (Hawkesford and Riche 2020). Further, Hasegawa (2003) and Büchi et al. (2016) argued that crop varieties selected under reduced N rates are not necessarily better adapted to low-input conditions and that breeders should instead devote a majority of their resources to multi-environment testing. These often-observed low genotype × N rate interactions in previous assessments of genotypic variation for NUE indicate that it is rather promising to increase the number of testing environments instead of increasing the number of different N rates when selecting for NUE (Brasier et al. 2020). Registration trials with multiple N rates per trial would presumably provide more insights although one would not expect any differences in the ranking of the varieties regarding genotype-specific N responses. Besides, it would be much more expansive and harder to manage considering the generally large number of varieties and environments in variety testing systems. Accordingly, we expect no different results regarding the variety ranking and respective approval for variety release. From these results, we conclude that a trial system, as evaluated in this study, is at least as efficient to select varieties with high NUE as a system with more N levels, but fewer environments.

The dominating share of environmental variation compared to the genotypic and G × E interaction variation shown in Fig. 5 was in line with results reported by Brasier et al. (2020) in WW, Anbessa et al. (2009) in barley and Stahl et al. (2019) in WOSR. In contrast to these results, Ivic et al. (2021) found that a genotypic variation from a Croatian study in WW varieties developed between 1936 and 2016 for NUE of GYLD which was 42.1% compared to the environmental of only 3.4%. This is a strongly biased result which can be explained by the large trend in breeding progress for NUE between 1936 and 2016, which inflated the genotypic variation. In our study, the genetic trend was taken into account such that we estimated unbiased genotypic variance components (conf. Eq. 1).

Among the cereal crops, WW had the largest and SB the lowest genotypic variation and, vice versa, WW the smallest and SB the largest environmental variation as shown in Fig. 5. This can be ascribed to the fact that in WW varieties with a large spread of baking quality, from fodder to elite types, were included with very different GPC and GYLD (Laidig et al. 2017a). However, it should be pointed out that this large genotypic variation for WW could not be fully used for selection, because genotypes with low GPC have higher yield and are more likely of fodder and not of baking quality whereas the reverse applies for genotypes with lower yield but with higher GPC restricting the range for selection. SB was grown under N rates corresponding to malting barley level to achieve high malting quality. This requires a rather low and well-balanced N fertilization, which is below optimum grain yield. Further, it is known that SB genotypes are not very different, because most of the present malting barley varieties are descendants of “Hana-type” varieties from Moravia in the early 1900s and of the semidwarf variety “Trumpf” in the early 1970s which was likely the reason for this low genotypic variability (Laidig et al. 2017b).

The generally low genotypic variation for NUE and NYLD shown in Fig. 5 raises the question if a further increase of NUE could be achieved by broadening the genotypic variation from a wider germplasm pool outside. The low genotypic variation might have different causes, such as the existence or fixation of unfavorable alleles for NUE in the elite breeding material, or the presence of unfavorable allele combinations which are affecting the highly quantitative trait. Both assumptions could be the result of a low selection pressure on the trait NUE and NYLD in the last decades due to higher fertilization rates. In case of the first assumption that a wider germplasm pool should be considered for future breeding efforts, by broadening the genetic basis of NUE and NYLD. However, the analysis of recent varieties (released and grown before mineral fertilizers were commonly used) or genetic resources like resynthesized lines as source for NUE genetic variation in WOSR was not considered as promising (Kessel et al. 2012). The latter assumption that best combinations of favorable alleles for NUE which already exist in the current elite material have not yet been either produced or found, implicates the need of larger population sizes in breeding programs for NUE to increase the number of meiosis and the chance for selecting best progenies.

This study revealed a strong influence of environmental conditions on total variation of NUE for GYLD/OYLDNUE and NYLDNUE of more than 90%. While trial-specific soil and climate conditions cannot be influenced, improved management including crop rotation, soil tillage and pesticide application may help to improve NUE. Most importantly, optimizing fertilization with regard to distribution, timing, type and amount, under consideration of pre-crop, seasonal weather course and actual crop demand provides substantial potential for increasing NUE in crop production.

Genotypic, environmental and G × E correlation

The very strong marginal correlation coefficients rM between NYLD and GYLD/OYLD, shown in Table 4, were in line with results reported by Ivic et al. (2021) and Guttieri et al. (2017) for WW, Anbessa et al. (2009) and Sinebo et al. (2004) for SB and Stahl et al. (2019) for WOSR. This strong marginal correlation between NYLD and GYLD/OYLD confirmed that breeding progress achieved in GYLD/OYLD also resulted in a higher NYLD, while the lower correlation with GPC showed that this trait is of lower influence on NYLD.

The fact that the genotypic correlation coefficient rG for GYLD/OYLD and NYLD with the corresponding traits for NUE is approximately one, as shown in Supplementary Material Table S2 and confirmed by the proof in Appendix, implies that the rank order of their genotypic values, or approximately equivalent between their least square estimates of variety means, is about the same. In numerous studies strong to very strong correlations coefficients between variety means for GYLDNUE with GYLD were reported indicating this very strong association, however, without mentioning that this is due to the functional relation between trait and its NUE (e.g., Ivic et al. 2021; Muurinen et al. 2006; Stahl et al 2019).

Heritability

Many studies on NUE were published, however, only a few reported heritability coefficients, but of very different magnitude in the range between 0.25 ≤ H2 ≤ 0.94, e.g., in WW by Hitz et al. (2017), Ivic et al. (2021), Guttieri et al. (2017) and in SB by Sinebo et al. (2004) and Anbessa et al. (2009). Compared to reported heritability coefficients, we estimated larger values for NUE as shown in Table 2. This spread in reported heritability coefficients may likely be explained, as mentioned previously, by the different sets of genotypes and environments from which they were estimated, which makes it difficult to compare them between studies. So far, NYLD of submitted varieties is not a direct criterion for registration, except in PEAS and in WOSR where protein yield is considered. The fact that the average heritability across crops for NYLD (H2 = 0.73) was not much lower than for GYLD/OYLD (H2 = 0.78), which is a very important criterion for registration, indicates that NYLD is about equally reliant as GYLD. Further, weak positive genotypic correlation coefficients between NYLD with GYLD and GPC, as shown in Table 4, give evidence that selection of varieties for high NYLD does not counteract the selection for high GYLD or GPC. For these reasons, we suggest to use NYLD as further criterion in registration trials, or, in the case of WOSR and PEAS where NYLD is already assessed, NYLD should be given a higher weight. Oberforster and Werteker (2005) already advocated for the use of NYLD in Austrian winter wheat registration trials as additional criterion. Reliability of NYLD as registration criterion could be improved if the number of locations at which GPC and NYLD is assessed is increased in WW, WR, SW and SB, beyond the subset of locations where quality samples are taken so far.

Conclusions

This study assessed breeding progress for NUE and related traits in important crops. We showed that nitrogen fertilization rates in variety trials of cereal crops and winter oilseed rape, except spring barley, were considerably reduced between 1995 and 2021. Despite this reduction, grain and oil yield increased while grain nitrogen yield did not decrease significantly in all crops. NUE for grain yield, oil and nitrogen yield increased strongly, which indicates a large breeding progress for NUE due to improved varieties. However, some of the increase may be attributed to the known effect that NUE increases when nitrogen fertilizer use is reduced. Genotypic variation was low compared to the environmental variation for grain, oil yield and for nitrogen yield, while the environmental variation of NUE for grain and nitrogen yield was even higher. The correlation coefficients showed that grain yield had a stronger influence on nitrogen yield than grain protein concentration. Furthermore, the low positive genotypic correlation of nitrogen yield with grain yield and grain protein concentration suggest that the selection for high nitrogen yield does not counteract with grain yield and grain protein concentration. The result that the ranking of the genotypic values of varieties for grain and nitrogen yield, which is approximately the same as the ranking of the least square estimates for variety means, was the same as for the genotypic values of their corresponding NUE traits, indicates that NUE for grain yield was already taken into account by grain yield as an important registration criterion. Heritability of cycle means for nitrogen yield was only little lower than for grain yield, which means that nitrogen yield is approximately as reliable as grain yield when using it as an additional trait for registration. Therefore, nitrogen yield should be given a higher weight in breeding and variety registration to increase NUE and reduce adverse environmental impact. This study’s results highlight that despite considerable reduction in nitrogen fertilizer inputs in cereal crops and oilseed rape breeding progress was achieved without comprising land use efficiency.

Electronic supplementary material

Below is the link to the electronic supplementary material.Supplementary file1 (PDF 67 KB)SM3 Plots between least square means of winter wheat varieties grown in the same cycle

Supplementary file2 (PDF 28 KB)SM1 Variance components for NUE traits and related traits

Supplementary file3 (PDF 27 KB)SM2 Correlation coefficients of NUE traits with related traits

Appendix

Why genotypic and G × E correlation coefficients for grain yield with NUE for grain yield and nitrogen yield with NUE for nitrogen yield are approximately equal to one

Consider two random variables Wij and Uij = Wij/Nj, where Wij = grain yield of ith genotype in jth environment and Nj = nitrogen input in jth environment. For Wij, we assume the random effects model10 Wij=μ+Gi+Lj+GLij

with independent random effects having zero mean and variances varGi=σG2, varLj=σL2 and varGLij=σGL2. Here, we will investigate what can be said about the covariance of effects for Wij and Uij assuming model (10) for Wij. For simplicity, we will assume that Nj is independent of all effects in the model for Wij. This could be relaxed by assuming a covariance with Lj, but this is not expected to have any substantive bearing on the subsequent derivation, as will be explained briefly at the end. The model for Uij can be written as11 Uij=μ+Gi+Lj+GLijNj=μj′+Gij′+Lj′+GLij′

where μj′=μNj, Gij′=GiNj, Lj′=LjNj and GLij′=GLijNj.

A challenge in the subsequent derivation is that all effects in (11) involve the random variable Nj, hence we need to check if there are correlations among the effects. It will be useful to make use of the laws of total variance and total covariance (Rudary 2009). Let X, Y and Z be three random variables. Then the law of total variance states thatvarY=EvarYX+varEYX.

The law of total covariance states thatcovX,Y=EcovX,YZ+covEXZ,EYZ.

If we identify Z with Nj, and X and Y with any of the two effects in (11), then it emerges from the law of total covariance that the effects in (11) must all be uncorrelated from one another. For example, we find thatcovGij′,Lj′=EcovGij′,Lj′Nj+covEGij′Nj,ELj′Nj=

E1Nj2covGi,Lj+cov0,0=0

Before we consider covariances among effects for Wij and Uij, we need to transition (11) into a model of the same form as (10). To do so, we may defineμj′=μ″+Fj″,

where μ″=Eμj′ and Fj″=μj′-μ″. Similarly, we defineGij′=Gi″+Hij″,

where Gi″=EGij′Gi and Hij″=Gij′-EGij′Gi. With these definitions, we can rewrite (11) as12 Wij=μ″+Gi″+Lj″+GLij″,

where Lj″=Lj′+Fj″ and GLij″=GLij′+Hij″. With this reparameterization, we can now study covariances among effects between models (10) and (12), using the law of total covariance and the delta method (Johnson et al. 1993).(i) covGi,Gi″=EcovGi,Gi″Nj+covEGiNj,EGi″Nj≈σG2φN

(ii) covLj,Lj″=EcovLj,Lj″Nj+covELjNj,ELj″Nj≈σL2φN

(iii) covGLij,GLij″=EcovGLij,GLij″Nj+covEGLijNj,EGLij″Nj≈σGL2φN

where φN=ENj. Similarly, the variances of the random effects in (12) can be approximated as follows:(i) varGi″≈EGi2φN4σN2+σG2φN2=σG2φN2

(ii) varLj″=varLj′+Fj″≈σL2φN2+μ2σN2φN4

(iii) varGLij″=varGLij′+Gij′-EGij′Gi≈σGL2φN2+varGi1Nj-1φN≈σGL2φN2

where σN2=varNj. From these results, the correlations can be approximated as(i) corrGi,Gi″≈1

(ii) corrLj,Lj″≈11+μ2σN2φN2σL2<1

(iii) corrGLij,GLij″≈1

As noted before, we have assumed here that Lj and Nj are uncorrelated. If a covariance is allowed, this will only affect the correlation corrLj,Lj″, but it will still be the case that corrLj,Lj″<1.

Abbreviations

GOC Grain oil concentration

GPC Grain protein concentration

GYLD Grain yield

GYLDNUE NUE of grain yield

Nmin Soil-mineralized nitrogen

NUE Nitrogen use efficiency

NYLD Nitrogen yield in grain

NYLDNUE NUE of nitrogen yield

OYLD Oil yield

PEAS Grain peas

SB Spring barley

SLF Soil fertility

SW Spring wheat

WOSR Winter oil seed rape

WR Winter rye

WW Winter wheat

WWORG Winter wheat organic regimen

Acknowledgements

FL and HPP were supported by DFG project PI 377/20-2. TF was supported by the project KlimaN funded by the Federal Ministry of Food and Agriculture (BMEL) as part of the German Climate Protection Program 2022. We thank Andrea Richter, Institute of Crop Science, University of Hohenheim, for carefully recording trial data from printed variety reports, the Bundessortenamt for providing the data and Thomas Drobek, Bundessortenamt, for preparing the data and making them available for this study.

Author Contribution

FL conceived the study, carried out the analyses, prepared the figures and tables and wrote the manuscript. TF contributed in refining the study goals, interpretation and presentation of results and writing of the manuscript. HPP supported with statistical advice and participated in editing the manuscript. CL and AS have read and commented on the manuscript.

Funding

Open Access funding enabled and organized by Projekt DEAL.

Data availability

Data were provided by the Federal Plant Variety Office for exclusive use in this study and are in general not publicly available. Reasonable requests may be addressed to the Federal Plant Variety Office, Hannover, Germany.

Declarations

Conflict of interest

The authors have no relevant financial or non-financial interests to disclose.

Ethical approval

The authors declare that the experiments comply with the current laws of Germany.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

Anbessa Y Juskiw P Good A Nyachiro J Helm J Genetic variability in nitrogen use efficiency of spring barley Crop Sci 2009 49 1259 1269 10.2135/cropsci2008.09.0566
Barmeier G Hu Y Schmidhalter U Partitioning and translocation of dry matter and nitrogen during grain filling in spring barley varieties and their roles in determining malting quality Front Plant Sci 2021 12 722871 10.3389/fpls.2021.722871 34497628
Barraclough PB Howarth JR Jones J Lopez-Bellido R Parmar S Shepherd CE Hawkesford MJ Nitrogen efficiency of wheat: genotypic and environmental variation and prospects for improvement Eur J Agron 2010 33 1 11 10.1016/j.eja.2010.01.005
Bingham IJ Karley AJ White PJ Thomas WTB Russell JR Analysis of improvements in nitrogen use efficiency associated with 75 years of spring barley breeding Europ J Agron 2012 42 49 58 10.1016/j.eja.2011.10.003
Blume HP et al (2015) Scheffer/Schachtschabel, soil science. Springer Berlin/Heidelberg. ProQuest Ebook Central
BMEL (2019) Discussion paper 2035. Arable farming strategy. Prospects for productive and diverse crop farming. Federal Ministry of Food and Agriculture (BMEL). Published Dec 2019. https://www.bmel.de/SharedDocs/Downloads/EN/Publications/ackerbaustrategie-en.pdf?__blob=publicationFile&v=6. Accessed 15 Feb 2022
BMEL (2022a) Tabellen aus dem Statistischen Jahrbuch über Ernährung, Landwirtschaft und Forsten und Tabellen des Monatsberichts zum Thema Landwirtschaft. https://www.bmel-statistik.de/landwirtschaft/tabellen-zur-landwirtschaft#c8262; SJT-3070400-0000.xlsx Landwirtschaftlich genutzte Fläche nach Kulturarten. Accessed 6 Oct 2023
BMEL (2022b) Tabellen aus dem Statistischen Jahrbuch über Ernährung, Landwirtschaft und Forsten und Tabellen des Monatsberichts zum Thema Landwirtschaft. https://www.bmel-statistik.de/landwirtschaft/tabellen-zur-landwirtschaft#c8273; MBT-0111260-0000 Nährstoffbilanz insgesamt von 1990 bis 2021 - in kg N/ha. Accessed 6 Oct 2023
BodSchätzG (2007) Gesetz zur Schätzung des landwirtschaftlichen Kulturbodens. https://www.gesetze-im-internet.de/bodsch_tzg_2008/BodSch%C3%A4tzG.pdf. Accessed 15 Feb 2022
Borlaug NE (1968) Wheat breeding and its impact on world food supply. In: Proceedings 3rd international wheat genetics symposium. Canberra. Aust. Acad. Sci., Canberra, pp 1–36
Bouchet AS Laperche A Bissuel-Belaygue C Snowdon RJ Nesi N Stahl A Nitrogen use efficiency in rapeseed. A review Agron Sustain Dev 2016 36 2029 10.1007/s13593-016-0371-0
Brasier K Oakes J Balota M Thomason W Griffey C Greater biomass accumulation at anthesis increases nitrogen use efficiency in winter wheat Agron J 2020 111 1 11 10.2134/agronj2018.10.0693
BSL (2023) Descriptive variety list. Cereal, maize, large grained pulse crops, root crops (except potato) (in German), Hannover, Bundessortenamt. https://www.bundessortenamt.de/internet30/filedmin/Files/PDF/bslgetreide_2023.pdf. Accessed 18 Oct 2023
Büchi L Charles R Schneider D Sinaj S Maltas A Fossati D Mascher F Performance of eleven winter wheat varieties in a long term experiment on mineral nitrogen and organic fertilisation Field Crops Res 2016 191 111 122 10.1016/j.fcr.2016.02.022
Capriel P Trends in organic carbon and nitrogen contents in agricultural soils in Bavaria (south Germany) between 1986 and 2007 Eur J Soil Sci 2014 64 445 454 10.1111/ejss.12054
Cassman KG Dobermann A Nitrogen and the future of agriculture: 20 years on Ambio 2022 51 17 24 10.1007/s13280 33715091
Connor DJ Organic agriculture cannot feed the world Field Crops Res 2008 106 187 190 10.1016/j.fcr.2007.11.010
Cormier F Faure S Dubreuil P Heumez E Beauchene K Lafarge S Praud S Le Gouis J A multi-environmental study of recent breeding progress on nitrogen use efficiency in wheat (Triticum aestivum L.) Theor Appl Genet 2013 126 3035 3048 10.1007/s00122-013-2191-9 24057081
Cormier F Foulkes J Hirel B Gouache D Moënne-Loccoz Y Le Gouis J Breeding for increased nitrogen-use efficiency: a review for wheat (T. aestivum L.) Plant Breed 2016 135 255 278 10.1111/pbr.12371
de Oliveira Silva A, Ciampitti IA, Slafer GA, Lollato RP (2020) Nitrogen utilization efficiency in wheat: a global perspective. Eur J Agron 114. 10.1016/j.eja.2020.126008
DIN EN 15948 (2012) Cereals—determination of moisture and protein—method using near-infrared-spectroscopy in whole kernels; German version EN 15948:2012 http://www.nal.din.de. Accessed 24 Feb 2023 Titel (englisch)
DUEV (2017) Verordnung über die Anwendung von Düngemitteln, Bodenhilfsstoffen, Kultursubstraten und Pflanzenhilfsmitteln nach den Grundsätzen der guten fachlichen Praxis beim Düngen (Düngeverordnung - DüV). http://www.gesetze-im-internet.de/dv_2017/. Accessed 15 Feb 2022
Einarsson R Sanz-Cobena A Aguilera E Billen G Garnier J van Grinsven HJM Lassaletta L Crop production and nitrogen use in European cropland and grassland 1961–2019 Sci Data 2021 8 288 10.1038/s41597-021-01061-z 34716345
Erisman J Sutton M Galloway J Klimont Z Winiwarter W How a century of ammonia synthesis changed the world Nature Geosci 2008 1 636 639 10.1038/ngeo325
EU (2020) European Commission. COM (2020) 381 final. A farm to fork strategy for a fair, healthy and environmentally-friendly food system. European Union. “Farm to Fork Strategy: For a fair, healthy and environmentally-friendly food system.” https://ec.europa.eu/food/farm2fork_en (2020). Accessed 24 Jan 2023
Good AG Shrawat AK Muench DG Can less yield more? Is reducing nutrient input into the environment compatible with maintaining crop production? Trends Plant Sci 2004 9 597 605 10.1016/j.tplants.2004.10.008 15564127
Guarda G Padovan S Delogu G Grain yield, nitrogen-use efficiency and baking quality of old and modern Italian bread-wheat cultivars grown at different nitrogen levels Eur J Agron 2004 21 181 192 10.1016/j.eja.2003.08.001
Guttieri MJ Frels K Regassa T Baenziger WBM S, Variation for nitrogen use efficiency traits in current and historical great plains hard winter wheat Euphytica 2017 213 87 10.1007/s10681-017-1869-5
Hartung J, Laidig F, Piepho HP (2023) Effects of systematic data reduction on trend estimation from German registration trials. Theor Appl Genet 136:0 10.1007/s00122-023-04266-5
Hasegawa H High-yielding rice cultivars perform best even at reduced nitrogen fertilizer rate Crop Sci 2003 43 921 926 10.2135/cropsci2003.9210
Häußermann U Klement L Breuer L Ullrich A Wechsung G Bach M Nitrogen soil surface budgets for districts in Germany 1995 to 2017 Environ Sci Eur 2020 32 109 10.1186/s12302-020-00382-x
Hawkesford MJ Reducing the reliance on nitrogen fertilizer for wheat production J Cereal Sci 2014 59 276 283 10.1016/j.jcs.2013.12.001 24882935
Hawkesford MJ Riche AB Impacts of G x E x M on nitrogen use efficiency in wheat and future prospects Front Plant Sci 2020 11 1157 10.3389/fpls.2020.01157 32903740
Hitz K Clark AJ Van Sanford DA Identifying nitrogen-use efficient soft red winter wheat lines in high and low nitrogen environments Field Crops Res 2017 200 1 9 10.1016/j.fcr.2016.10.001
Ivic M Grljusic S Plavsin I Dvojkovic K Lovric A Rajkovic B Maric evic M, Cerne, Popovic B, Loncaric Z, Bentley AR, Swarbreck SM, Sarcevic H, Novoselovic D, Variation for nitrogen use efficiency traits in wheat under contrasting nitrogen treatments in South-Eastern Europe Front Plant Sci 2021 12 682333 10.3389/fpls.2021.682333 34868096
Johnson NL Kotz S Kemp AW Univariate discrete distributions 1993 2 New York Wiley
Johnston AE Poulton PR Coleman K Soil organic matter: its importance in sustainable agriculture and carbon dioxide fluxes Adv Agron 2009 101 1 57 10.1016/S0065-2113(08)00801-8
Kelstrup L, Rowarth JS, Williams PH, Ronson C (1996) Nitrogen fixation in peas (Pisum sativum L.), lupins (Lupinus angustifolius L.) and lentils (Lens culinaris Medik.) Proc Agron Soc NZ 26:71–74
Kessel B Schierholt A Becker HC Nitrogen use efficiency in a genetically diverse set of winter oilseed rape (Brassica napus L.) Crop Sci 2012 52 2546 2554 10.2135/cropsci2012.02.0134
Ladha JK Tirol-Padre A Reddy CK Cassman KG Verma S Powlson DS van Kessel C de Richter D Richter DB Chakraborty D Pathak H Global nitrogen budgets in cereals: a 50-year assessment for maize, rice, and wheat production systems Sci Rep 2016 6 19355 10.1038/srep19355 26778035
Laidig F Piepho HP Drobek T Meyer U Genetic and non-genetic long-term trends in 12 different crops in German official variety performance trials and on-farm yield trends Theor Appl Genet 2014 127 2599 2617 10.1007/s00122-014-2402-z 25307935
Laidig F Piepho HP Rentel D Drobek T Meyer U Huesken A Breeding progress, environmental variation and correlation of winter wheat yield and quality traits in German official variety trials and on-farm during 1983–2014 Theor Appl Genet 2017 130 223 245 10.1007/s00122-016-2810-3 27796431
Laidig F Piepho HP Rentel D Drobek T Meyer U Breeding progress, genotypic and environmental variation and correlation of quality traits in malting barley in German official variety trials between 1983 and 2015 Theor Appl Genet 2017 130 2411 2429 10.1007/s00122-017-2967-4 28821914
Laidig F Feike T Klocke B Macholdt J Miedaner T Rentel D Piepho HPP Long-term breeding progress of yield, yield-related, and disease resistance traits in five cereal crops of German variety trials Theor Appl Genet 2021 134 3805 3827 10.1007/s00122-021-03929-5 34652455
Laidig F Feike T Hadasch S Rentel D Klocke B Miedaner T Piepho HP Yield reduction due to diseases and lodging and impact of input intensity on yield in variety trials in five cereal Crops Euphytica 2022 218 150 10.1007/s10681-022-03094-w
Laperche A Devienne-Barret F Maury O Le Gouis J Ney B A simplified conceptual model of carbon/nitrogen functioning for QTL analysis of wheat adaptation to nitrogen deficiency Theor Appl Genet 2006 113 1131 1146 10.1007/s00122-006-0373-4 16909280
Lassaletta L Billen G Grizetti B Anglade J Garnier J 50 year trends in nitrogen use efficiency of world cropping systems: The relationship between yield and nitrogen input to cropland Environ Res Lett 2014 9 105011 10.1088/1748-9326/9/10/105011
Le Gouis J Beghin B Heumez E Pluchard P Genetic differences for nitrogen uptake and nitrogen utilisation efficiencies in winter wheat Eur J Agron 2000 12 163 173 10.1016/S1161-0301(00)00045-9
Mackay IJ Horwell A Garne RJ White J McKee J Philpott H Reanalysis of the historical series of UK variety trials to quantify the contributions of genetic and environmental factors to trends and variability in yield over time Theor Appl Genet 2011 122 225 238 10.1007/s00122-010-1438-y 20835813
Miersch S Gertz A Breuer F Schierholt A Becker HC Influence of the semi-dwarf growth type on Nitrogen use efficiency in winter oilseed rape Crop Sci 2016 56 2952 2961 10.2135/cropsci2016.01.0044
Moll RH Kamprath EJ Jackson WA Analysis and interpretation of factors which contribute to efficiency of nitrogen utilization Agron J 1982 74 562 564 10.2134/agronj1982.00021962007400030037x
Muurinen S Slafer GA Peltonen-Sainio P Breeding effects on nitrogen use efficiency of spring cereals under northern conditions Crop Sci 2006 46 561 568 10.2135/cropsci2005-05-0046
Neuweiler JE Trini J Maurer HP Würschum T Do lower nitrogen fertilization levels require breeding of different types of cultivars in triticale? Theor Appl Genet 2022 135 993 1009 10.1007/s00122-021-04012-9 34958397
Oberforster M, Werteker M (2005) Stickstoff-Effizienz als Parameter des Sortenwertes von Winterweizen. Bericht 60. ALVA-Jahrestagung, 23.-25. Mai, Linz. Arbeitsgemeinschaft Lebensmittel- Veterinär- und Agrarwesen (ALVA), Wien, pp 108–112
Ortiz-Monasterio I Sayre KD Rajaram S McMahon M Genetic progress in wheat yield and nitrogen use efficiency under four N rates Crop Sci 1997 37 898 904 10.2135/cropsci1997.0011183X003700030033x
Piepho HP Allowing for the structure of a designed experiment when estimating and testing trait correlations J Agric Sci 2018 156 59 70 10.1017/S0021859618000059
Piepho HP Möhring J Computing heritability and selection response from unbalanced plant breeding trials Genetics 2007 177 1881 1888 10.1534/genetics.107.074229 18039886
Piepho H-P Laidig F Drobek T Meyer U Dissecting genetic and non-genetic sources of long-term yield in German official variety trials Theor Appl Genet 2014 127 1009 1018 10.1007/s00122-014-2275-1 24553961
Piepho HP Gabriel D Hartung J Büchse A Grosse M Kurz S Laidig F Michel V Proctor I Sedlmeier JE Toppel K Wittenburg D One, two, three: Portable sample size in agricultural research J Agric Sci 2022 160 6 459 482 10.1017/S0021859622000466
Piepho HP (2019) A coefficient of determination (R 2) for generalized linear mixed models. Biometrical Journal 61:860–872–13. 10.1002/bimj.201800270
Pituello C Polese R Morari F Berti A Outcomes from a long-term study on crop residue effects on plant yield and nitrogen use efficiency in contrasting soils Eur J Agron 2016 77 179 187 10.1016/j.eja.2015.11.027
Raun WR Johnson GV Improving nitrogen use efficiency for cereal production Agron J 1999 91 357 363 10.2134/agronj1999.00021962009100030001x
Riedesel L Laidig F Hadasch S Rentel D Hackauf B Piepho HP Feike T Breeding progress reduces carbon footprints of wheat and rye J Clean Prod 2022 377 1 14 10.1016/j.jclepro.2022.134326
Rudary MR (2009) On predictive linear Gaussian models. PhD thesis, ProQuest LLC, Ann Arbor, MI, USA
Ruisi P Giambalvo D Di Miceli G Frenda AS Saia S Amato G Tillage effects on yield and nitrogen fixation of legumes in Mediterranean conditions Agron J 2012 104 1459 1466 10.2134/agronj2012.0070
Rütting T Delin AH S, Efficient use of nitrogen in agriculture Nutr Cycl Agroecosyst 2018 110 1 5 10.1007/s10705-017-9900-8
Sieling K Kage H Apparent fertilizer N recovery and the relationship between grain yield and grain protein concentration of different winter wheat varieties in a long-term field trial Eur J Agron 2021 124 126246 10.1016/j.eja.2021.126246
Sieling K Kage H Winter barley grown in a long-term field trial with a large variation in N supply: Grain yield, yield components, protein concentration and heir trends Eur J Agron 2022 136 126505 10.1016/j.eja.2022.126505
Sinebo W Gretzmacher R Edelbauer A Genotypic variation for nitrogen use efficiency in Ethiopian barley Field Crops Res 2004 85 43 60 10.1016/S0378-4290(03)00135-7
Stahl A Pfeifer M Frisch M Wittkop B Snowdon RJ Recent genetic gains in nitrogen use efficiency in oilseed rape Front Plant Sci 2017 8 963 10.3389/fpls.2017.00963 28638399
Stahl A Vollrath P Samans B Frisch M Wittkop B Snowdon RJ Effect of breeding on nitrogen use efficiency-associated traits in oilseed rape J Exp Bot 2019 70 1969 1986 10.1093/jxb/erz044 30753580
Stewart WM Dibb DW Johnston AE Smyth TJ The contribution of commercial fertilizer nutrients to food production Agron J 2005 97 1 6 10.2134/agronj2005.0001
Sylvester-Bradley R Kindred DR (2009) Analysing nitrogen responses of cereals to prioritize routes to the improvement of nitrogen use efficiency J Exp Bot 2009 60 1939 1951 10.1093/jxb/erp116 19395389
van Grinsven HJM Ebanyat P Glendining M Gu B Hijbeek R Lam SK Lassaletta L Mueller ND Pacheco FS Quemada M Bruulsema TW Jacobsen BH ten Berge HFM Establishing long-term nitrogen response of global cereals to assess sustainable fertilizer rates Nat Food 2022 3 122 132 10.1038/s43016-021-00447-x 37117954
Voss-Fels KP Stahl A Wittkop B Lichthardt C Nagler S Rose T Breeding improves wheat productivity under contrasting agrochemical input levels Nat Plants 2019 5 706 714 10.1038/s41477-019-0445-5 31209285
Weiser C Fuß R Kage H Flessa H Do farmers in Germany exploit the potential yield and nitrogen benefits from preceding oilseed rape in winter wheat cultivation? Arch Agron Soil Sci 2018 64 25 37 10.1080/03650340.2017.1326031
Yan M Pan G Lavallee JM Conant RT Rethinking sources of nitrogen to cereal crops Glob Change Biol 2019 26 191 199 10.1111/gcb.14908
Yang C Bueckert R Schoenau J Diederichsen A Zakeri H Warkentin TD Evaluation of growth and nitrogen fixation of pea nodulation mutants in western Canada Can J Plant Sci 2017 97 1121 1129 10.1139/cjps-2016-0383
Zetzsche H Friedt W Ordon F Breeding progress for pathogen resistance is a second major driver for yield increase in German winter wheat at contrasting N levels Sci Rep 2020 10 20374 10.1038/s41598-020-77200-0 33230232
