ORIGINAL RESEARCH article
Applicability of a “Multi-Stage Pulse Labeling” 15N Approach to Phenotype N Dynamics in Maize Plant Components during the Growing Season
- 1Department of Agronomy, Kansas State University, Manhattan, KS, United States
- 2Agronomy Department, Purdue University, West Lafayette, IN, United States
- 3Australian Grain Technologies, Adelaide, SA, Australia
• This work utilizes “multi-stage pulse labeling” 15N applications, primarily during reproductive growth stages, as a phenotyping strategy to identify maize hybrids with superior N use efficiency (NUE) under low N conditions.
Research using labeled isotopic N (15N) can precisely quantify fertilizer nitrogen (N) uptake and organ-specific N allocation in field crops such as maize (Zea mays L.). The overall research objective was to study plant N uptake patterns potentially correlated with N use efficiency (NUE) in field-grown maize hybrids using a “multi-stage pulse labeling” 15N phenotyping strategy with an emphasis on the reproductive period. Five hybrids varying in NUE were compared under zero N fertilizer application (0N) plus a moderate rate of 112 kg N ha−1 (112N) in 2013 (2 locations) and 2014 growing seasons. The equivalent of 3.2 (2013) to 2.1 (2014) kg of 15N ha−1, as labeled Ca(15NO3)2, was injected into soil on both sides of consecutive plants at multiple stages between V14 and R5. Aboveground plant biomass was primarily collected in short-term intervals (4–6 days after each 15N application) in both years, and following a single long-term interval (at R6 after 15N injection at R1) in 2014. Averaged across hybrids and site-years, the moderate N rate (112N) increased absolute 15N uptake at all stages; however, plants in the 0N treatment allocated proportionally more 15N to reproductive organs. Before flowering, short-term recovery of 15N (15Nrec) totaled ~0.30 or 0.40 kg kg−1 of the 15N applied, and ~50% of that accumulated 15Nu was found in leaves and 40% in stems. After flowering, plant 15Nrec totaled ~0.30 kg kg−1 of 15N applied, and an average 30% of accumulated 15Nu was present in leaves, 17% in stems, and the remainder—usually the majority—in ears. At the R5 stage, despite a declining overall rate of 15N uptake per GDD thermal unit, plant 15Nrec represented ~0.25 kg kg−1 of 15N applied, of which ~65% was allocated to kernels. Overall long-term 15Nrec during grain filling was ~0.45 and 0.70 kg kg−1 of total 15N applied at R1 with 0 and 112N, respectively, and most (~77%) 15N uptake was found in kernels. The “multi-stage pulse labeling” technique proved to be a robust phenotyping strategy to differentiate reproductive-stage N uptake/allocation patterns to plant organs and maize efficiencies with newly available fertilizer N.
Maize is one of the most important crops cultivated for food and feed production worldwide and it is also important for biofuel production in some countries. Maize requires significant amounts of N to maximize yield (Coelho et al., 1991; Fageria and Baligar, 2005). Application of N fertilizer increased over the years from 1960 to 1980 (Cassman et al., 2002; Tilman et al., 2002; Van Cleemput et al., 2008; Scharf, 2015), but since then overall N rates utilized for maize production in the United States have stabilized (USDA-ERS, 2013). However, maize may still only take up about 50% of total N fertilizer applied (Dobermann and Cassman, 2004; Silva et al., 2006; Van Cleemput et al., 2008; Ciampitti and Vyn, 2014). Increasing the efficiency by which plants use available N is crucial to optimize crop yield potential, reduce the costs of N inputs, and subsequently avoid N losses to the environment (Raun and Johnson, 1999; Stevens et al., 2005).
The quantitative nature of genes controlling N use efficiency (NUE) integrated with the complex N cycles in soil-plant systems and in alternate crop management and environmental scenarios make it challenging to improve NUE (Amado et al., 2002; Mosier et al., 2004; Fageria and Baligar, 2005; Ladha et al., 2005; Mi et al., 2005; Coque et al., 2008). Previous improvements in maize NUE in field research programs have largely been based on grain yield and N uptake evaluations of diverse genetic materials advanced from both traditional and transgenic breeding approaches (Guo et al., 2014) with or without integrated agronomic practices (often involving multiple N rates) in various environmental conditions (Duvick, 1984; Ciampitti and Vyn, 2011). Because NUE is directly related to yield improvement, NUE gains occur when maize plants take up more N from soil and fertilizer N sources [preferably with a high N recovery efficiency (NRE)] and, thereafter, produce maximum grain yield per unit of plant N uptake (known as N internal efficiency; NIE) (Moll et al., 1982; Cassman et al., 2003; Dobermann, 2005; Mi et al., 2005; Coque and Gallais, 2007; Ciampitti and Vyn, 2012).
Raun and Johnson (1999) stated that low NUE (≈33 kg grain kg−1) in crop production was a consequence of the excessive use of N fertilizer. Thus, to optimize crop yield and improve NUE it is necessary to develop strategies that synchronize optimal timings of N application with plant N requirements (Cassman et al., 2002; Tilman et al., 2002; Dobermann, 2005; Fageria and Baligar, 2005; Mueller et al., 2014). Regardless of whether NUE is determined by the simple or the more robust difference method, these methods are unable to distinguish the in-season soil vs. fertilizer sources of N uptake, allocation, and remobilization by plants. The use of isotopic N (15N) allows closer estimation of the current N uptake by plants. Additionally, multiple 15N applications during the growing season can help to precisely determine the fate of the most recent inorganic N uptake and its initial organ-specific allocation inside plants.
Most commonly, 15N experiments involving crop response are developed under controlled environments in greenhouses (Pan et al., 1986; Schmidt and Scrimgeour, 2001; Rimski-Korsakov et al., 2009) or in the field by the use of chambers or lysimeters (Portela et al., 2006; O'Brien et al., 2012). In the greenhouse, Friedrich and Schrader (1979), evaluated maize N remobilization patterns during the grain filling period. They applied labeled 15N as nutrient solution to maize plants during the entire vegetative period up to the R1 stage and harvested plants from 1 to 7 weeks after R1. More recently, Paponov and Engels (2005), evaluated plant N partitioning also in a greenhouse study. They applied 15N to the soil surface 12 days after silking and harvested maize plants 14 and 47 days after silking.
Research with 15N techniques in large-scale field experiments has been limited due to the labeled product costs and elaborate sampling procedure in mass spectrometry analysis, which requires specialized technical assistance. Because of high costs of 15N products, researchers using isotopic N in open and/or closed systems usually apply 15N jointly with N fertilizer treatments at low enrichment rates (~5% atom 15N) or use 15N-depleted fertilizers (Coelho et al., 1991; Boaretto et al., 2004; Gallais et al., 2006).
In some cases, fertilizers incorporating labeled N were applied early- to mid-season, so that studies could examine N accumulation, remobilization, and recovery in maize plants at flowering or physiological maturity (Cliquet et al., 1990a; Ma et al., 1998; Stevens et al., 2005; Gallais et al., 2006; Silva et al., 2006; Ning et al., 2017). Several experiments conducted in the field have also used 15N with traditional N management practices to analyze the “long-term” N accumulation in the plants within a season (Coelho et al., 1991; Gallais et al., 2006; Duete et al., 2009). Gallais et al. (2006), mainly focused on the evaluation of sources of N accumulated in the grain at physiological maturity. They used two 15N application techniques: 15N was applied at either the beginning of stem elongation or at R1 and plants were harvested multiple times after these two applications (at R1, R1+ 15, 25, and 35 days, and at R6). They found that up to 55% of grain N originates from post-silking N uptake, while the remaining comes from N remobilization from the vegetative organs. Such “long-term” approaches seems to be more appropriate to quantify the remobilization of N accumulated in vegetative organs to the grain at maturity and the total N accumulation during a certain period than to estimate the fate of newest N uptake within the plant during the reproductive period.
There has only been a single field study in which 15N was applied to maize more than once in the growing season with plants harvested shortly after the 15N application. In that single-year and single-location experiment, Ta and Weiland (1992), compared two historic hybrids (crosses with a common B73 inbred) with different leaf senescence characteristics at two pre-plant N rates. The authors applied 50 mg of 15N to the soil beside single maize plants at three development stages (V14, R1, and R1 + 22 days) and these individual plants were harvested 3 days later (and at 35, 45, and ~67 days after the R1 application). Their study estimated that ~45% of the 15N applied was absorbed within 3 days and that more of the labeled N went to the ears with later 15N application times, or a longer remobilization period. Biomass and N uptake measurements from single plants, though indicative of trends, are never as precise as those from larger plant samples.
To date, the Ta and Weiland's (1992), study above is the only known report involving 15N applied to field-based maize a single time after the R1 stage. The knowledge gap concerning the most recent N uptake during the reproductive period is important because modern maize hybrids respond to N much differently than hybrids of earlier decades. More recent hybrids generally have higher N internal efficiency, and much more of their total N uptake occurs post-flowering (Ciampitti and Vyn, 2012, 2013; Chen et al., 2015; Mueller and Vyn, 2016). The gains in total N uptake with modern hybrids is not because of any higher leaf N concentrations at flowering (they may even be lower: Chen et al., 2015), but mostly because of the enhanced N accumulation during reproductive growth (Chen et al., 2015; Mueller and Vyn, 2016). Just how late in the grain filling period new N is taken up from the soil with more recent hybrids is still unidentified. Precise measurements of the efficiency with which modern maize hybrids uptake and allocate new N in its organs at multiple reproductive stages remains unknown.
Integrating cutting-edge cropping systems with high precision management and more resource-efficient genotypes are fundamental to increase crop production. New phenotyping approaches are necessary to identify NUE differences among maize hybrids during the growing season, and to better understand the when and where questions of hybrids that are supposedly superior in fertilizer N efficiencies.
The present study attempted to validate the utilization of the “multi-stage pulse labeling” 15N approach in field-grown maize as an advanced phenotyping technique to identify key mechanisms in N uptake, allocation, and partitioning during the growing season among hybrids with possible NUE variation. In this research, we modify earlier methods of tracer 15N utilization to more precisely estimate the fate of recent N uptake, partitioning, and recovery in maize plant components under different levels of N stress throughout the reproductive period. The research objectives were (i) to utilize “multi-stage pulse labeling” 15N application as a high precision phenotyping technique in maize hybrids grown at low N rates; (ii) to quantify 15N proportional allocation, and 15N fertilizer recovery, in plant components just before and at least four times during the reproductive growth period; and (iii) to examine how DM accumulation in discrete growth intervals impacts total plant 15N uptake activity throughout the growing season.
Materials and Methods
Weather Description, Management Practices, and Experimental Design
A non-irrigated research study was conducted during two growing seasons (2013–2014) in the US Midwestern Corn Belt region. In 2013, field experiments were established at the Purdue University Agronomy Center for Research and Education (ACRE) near West Lafayette—IN (Lat 40.486675° Lon 87.004635°, elevation 216 m) and at Pinney-Purdue Agricultural Center (PPAC) near Wanatah—IN (Lat 41.445113° Lon 86.943464°, elevation 222 m). At ACRE, soil was a Chalmers silty—clay loam (Fine—silty, mixed, superactive, mesic Typic Endoaquolls) and at PPAC soil was Sebewa loam (Fine—loamy over sandy or sandy—skeletal, mixed, superactive, mesic Typic Argiaquolls). In 2014, one experiment was established at ACRE (Lat 40.493593°, Lon 86.493593°, elevation 216 m). Soil type was a Raub—Brenton complex (Fine-silty, mixed, superactive, mesic Aquic Argiudolls).
Maize field trials were established following soybean [Glycine max (L.) Merr.] at all three site-years. In 2013, tillage operations prior to planting maize involved fall chisel plow followed by a spring field cultivator. For the 2014 site, fall plus spring strip-tillage on the no-till soybean stubble preceded maize planting (with precision guidance in all operations).
Experiments were planted with a four—row precision planter (Almaco SeedPro 360) with 76 cm row spacing, to achieve a final plant density of ~79,000 plants ha−1. In 2013 plot lengths were shorter (6.7 m) than in 2014 (13.7 m) due to seed supply limitations. Plant populations were evaluated in 5.3 m row length sections of all four rows at the V5 growth stage.
Weather data was acquired on a daily basis from planting to harvest (1 April—May to 31 October 2013 and 2014) from automated weather stations operated by the Indiana State Climate Office proximately located to the research sites. The field experiment at ACRE in 2013 was planted on 14 May and harvested 10 October. The total growing season precipitation was 385 mm; maximum and minimum daily average air temperatures were 26.9 and 14.6°C, respectively, for the entire growing season (Table 1). The field experiment at PPAC in 2013 was planted 5 June and harvested 30 October. Total growing season precipitation was 626 mm; maximum and minimum daily average temperatures were 24.7 and 12.1°C, respectively, for the entire growing season (Table 1). The field experiment at ACRE in 2014 was planted on 25 April and harvested 25 September. Total precipitation was 592 mm and maximum and minimum daily temperatures averaged 26.2 and 13.6°C throughout the growing season (Table 1). Weed pressure was minimal at all locations.
Table 1. Pertinent phenology and climate information of isotopic N experiments conducted at ACRE and PPAC locations in 2013 and at ACRE in 2014.
The two N rate treatments were 0 (0N), representing soil indigenous N pool and high N stress, and 112 kg N ha−1 (112N) as a moderate level of stress; the latter was sidedressed at the V4 stage between corn rows with a DMI Nutri-Placr at 10–12 cm soil depth as Urea Ammonium Nitrate (UAN) (28-0-0). The 0N plots received the same machine pass to avoid variations in plant growth performance due to soil compaction.
Comparisons were made between four modern hybrids varying in NUE with a similar 114 relative maturity range (Dow AgroSciences, Inc., Indianapolis, IN) and one historic hybrid from the 1970's. In this paper, results will be presented as combined means of genotypes to place emphasis on the “multi-stage pulse labeling” technique; a future paper will focus on the methodology's ability to separate hybrid differences.
At all sites, the plot arrangement was a split-plot experimental design, consisting of six replications and two treatment factors. Nitrogen rates were the main plots and hybrids were the sub-plots. The experiment consisted of six replications (because of cost constraints and small plot size, only three replications were used at a time for 15N application).
Statistical analyses were executed with SAS GLM (SAS Institute, 2014) at each growth stage and each year separately. The whole unit error (N rate) was pooled with the subunit error (Hybrid) for all ANOVAs because it was not significant for the majority of analyses (p > 0.25). A combined location analysis of variance was performed in 2013 for each variable for all growth stages where the error variances were homogeneous for the majority of the stages (p > 0.01). Data presentation in this publication is confined to the combined means of both locations in 2013 (as there were no significant treatment differences between locations in a majority of the observations), and the mean of ACRE in 2014. Location by treatment interactions were pooled with the experimental error when not significant (p > 0.01). Fisher's protected least significant difference (LSD) test was used to compare treatment means for each growth stage where the corresponding ANOVA F-test was significant (p ≤ 0.05).
Differences between growth stages were analyzed with Tukey's HSD test (p ≤ 0.05) using the R program (R Development Core Team, 2005).
Method of 15N Application and Biomass Harvesting
Near the beginning of the growing season, single-row 15N microplots were identified in three replications and were comprised of five (2013) or six (2014) consecutive plants at the appropriate density with sufficient undisturbed border plants (i.e., beyond 1 m from the microplot) so as to maintain a uniform canopy surrounding each plant group.
Microplots were labeled with 3.17 kg of 15N ha−1 in 2013 and 2.12 kg of 15N ha−1 in 2014 as Ca (15NO3)2 containing 98 atom % 15N (SIGMA-ALDRICH Co., St. Louis, MO). The labeled fertilizer provided 0.041 or 0.027 g of 15 N plant−1, respectively, for 2013 and 2014 growing seasons which represents a very small portion relative to the main N rate of 112 kg ha−1 (equivalent to 1.42 g of N plant−1). The 15N rate was lowered in 2014 to reduce costs after we learned that we could detect 15N concentration differences in all plant components during the entire 2013 season. Our 15N rates were similar to the 3.0 kg of 15N ha−1 applied by Ta and Weiland (1992) with a much lower plant density of ~60,000 plants ha−1.
Injection of 15N occurred in multiple plant development stages (from late vegetative stage to near the end of the grain filling period (Table 1). Prior to each time of 15N application, the labeled fertilizer was diluted in water (~300 ml of water per gram of Ca (15NO3), placed in 30 ml plastic syringes, properly sealed, and taken to the field. In the field, screwdrivers were used to make holes ~15 cm deep into the soil at a distance of 15 cm perpendicular to the row on both sides of each plant in the microplot. The 15N was injected into the hole.
Immediately after the 15N application, PVC pipes (0.6 diameter and 30 cm length) were installed on top of each hole and 0.8 l of water was applied (for a total of 1.6 l per plant). The added water was intended to ensure immediate availability of labeled nitrate-N to the plant roots.
Application dates for 15N during the late vegetative stages were determined based on when 50% of plants in the entire experiment had fully expanded leaf collar for each targeted stage. Application timing during reproductive stages occurred when over 50% of plants were in early R1 (ER1) (50% silking), R1 (90% silking), R2 (kernel blister), R4 (kernel dough), R5 (kernel dent) (Abendroth et al., 2011). Development stages of biomass harvesting varied slightly among locations and years (Table 1). In general, development stages refer to the time of plant biomass harvest 4–6 days after the 15N application (except for the V16 stage at ACRE in 2013 that indicates the time of 15N application).
At ACRE in 2013, plants treated with 15N at the V16 stage were in the field for 6 days (from 15N application 17 July to biomass harvesting 23 July) and reproductive development was so rapid that 90% of plants in the entire experiment had completed silking by the time of biomass harvest. Plants treated with 15N at the ER1 stage were harvested 4 days later (23 July). Subsequently, the microplots labeled as R1 received 15N when 90% of plants in the experiment had completed silking. For PPAC in 2013, we achieved better separation of the growth stages close to flowering (in part because of cooler temperatures). Thus, because the V16-stage samples at ACRE also covered the tassel formation period we decided it was still appropriate to combine the V16 results at ACRE with the V16 at PPAC. However, ER1 stage was evaluated only at ACRE and the R1 stage represents the combined means for ACRE and PPAC in 2013.
For both locations in 2013, the six actual biomass harvesting times and further 15N evaluations were generally completed at the same growth stages except that the V14 harvesting only occurred at PPAC and the early R1 harvesting occurred only at ACRE. Therefore, the 2013 results from V16, R1, R2, R4, and R5 stages consistently represent the means of the two locations (Table 2). In 2014, when there was only one location (ACRE), one of the pre-R1 stage treatments was substituted for a duplicate R1 application time that was then not harvested until the R6 stage (Table 1).
Table 2. Effects of overall sidedress N rate (0 and 112 kg N ha−1) on 15N fertilizer recovery (15Nrec, kg kg−115N applied) and proportional allocation of 15N uptake (15Np, %) in plant components at multiple development stages in 2013.
Destructive aboveground biomass of maize was collected from the three (2013) or four (2014) middle plants in the previously selected micro-plots from 4 to 6 days after the 15N application. One additional micro-plot application was implemented in 2014 whereby plants receiving 15N at R1 stage were not harvested until the R6 stage (R1/R6). Plant biomass was removed by cutting at the base of the stems with the use of a lopper. Plants were partitioned into stems + leaf sheaths + tassels (stems), leaf blades (leaves), husks + shank + silks (husks), and either whole ear (prior to R4) or cobs and kernels at R5 in 2013, and at R4, R5, and R6 stages in 2014. Stems were chopped into sections immediately after harvest and leaf removal. All samples were dried for 7 days at 60°C, weighed and coarse ground. Stems, leaves and husks were ground to pass a 2 mm screen and cobs and kernels were ground to a pass 4 mm screen size. Sub-samples were taken to the isotope laboratory to proceed with the mass spectrometry analysis.
Calculations of 15N Abundance in Plant Tissue
In order to estimate 15N uptake (15Nu) in plant components, N concentration (Nc) and delta 15N (δ15N ‰) were determined in the mass spectrometer for each sample. Nitrogen concentration represented the amount of 14N + 15N in the sample. Delta 15N (δ15N ‰) was the proportional amount of 15N in parts per mil in the sample relative to an international isotope standard of Air (0.366 atom % 15N) (Shearer and Kohl, 1986; Fry, 2006; Van Cleemput et al., 2008).
Additionally, δ15N ‰ is also described as the 15N atom% excess and can be calculated in the mass spectrometer through the isotopic ratio (IR =15N/14N) in the compounds. Absolute amounts of isotopic N cannot be simply estimated due to its atomic unit (Mariotti, 1983). Therefore, the linear relationship between 15N contents and δ15N values allows accurate estimation based on the isotope ratio (IR) approach in the mass spectrometer (Fry, 2006).
The isotopic ratio in a standard (IRstrd) represents the amount of 15N in the standard relative to the 15N abundance in the atmospheric dinitrogen (N2) defined as 0.366 atom % 15N (Faust, 1981; Shearer and Kohl, 1986; Dawson et al., 2002; Fry, 2006; Van Cleemput et al., 2008). The isotopic ratio in the sample (IRsmp) represents the quantity of 15N in the sample relative to the amount of 15N in a given standard.
The equations used to calculate 15Nu in plant components are as follows:
1. The isotopic ratio in the samples (IRsmp) was calculated through the proportional 15N abundance (δ15N) in the plant samples relative to the 15N abundance in the atmospheric dinitrogen (N2) defined as 0.0036764 ‰ (IRstrd) (Fry, 2006; Van Cleemput et al., 2008) (Equation 1):
2. Converting the IR to the absolute number of 15N atoms in 100 atoms of the total Nc in the sample. The absolute amount of atom % 15N (at15N) was calculated by dividing the proportional amount of 15N in the sample (IRsmp) by the total N in the sample (IRsmp + 1). (Equation 2):
3. Total N uptake (Nu) per plant component was estimated per unit area (kg ha−1) by multiplying the aboveground plant biomass (DM) in kg ha−1 by the N concentration (Nc) for respective plant components (Equation 3).
4. Absolute amounts of 15N uptake (15Nu) in plant components per unit area (kg ha−1) were estimated by multiplying Nu by the 15N found in respective plant component (Equation 4):
5. Proportional 15N uptake (15Np) for each plant component was estimated by dividing the 15Nu for respective plant component by the total plant 15Nu (Equation 5).
6. 15N recovery per plant component (15Nrec) was estimated per unit area (kg−1 kg of 15N applied ha−1) by dividing the total plant 15Nu (kg ha−1) by the amount of labeled fertilizer applied at 3.17 or 2.12 kg of 15N Ca (15NO3)2 per ha, respectively, for 2013 and 2014 (Equation 6).
Proportional Allocation of the Total Aboveground 15N Uptake in Plant Components, and Plant 15N Recovery Efficiency, Per Growth Stage
Results on the fate of the recently applied 15N fertilizer in the above-ground plant components are presented using three parameters: total 15N uptake per unit area (15Nu, kg ha−1), the proportional allocation of the total 15N uptake (15Np, %) in each separate plant component, and plant 15N recovery efficiency (15Nrec, kg kg−115N applied) of the total labeled fertilizer applied (3.2 or 2.1 kg of 15N ha−1) at each application time for 2013 and 2014, respectively. Hybrid differences in 15Nu and 15Np of individual plant components, as well as 15Nrec, for each stage will not be discussed in detail in this paper but the specific results for each hybrid, and the full main, sub and interaction statistics for each stage are available in the Supplementary Tables S3 and S4. Results were averaged across hybrid treatment for each N rate (0 and 112N). Supplementary Figure S1 demonstrates hybrid impacts on 15N allocation to plant components at R1 and R5 stages in 2013.
In 2013, at the V14 stage, mean 15Nrec (averaged across five hybrids) represented 0.25 and 0.30 kg kg−1 of the total 15N applied, respectively for 0 and 112N (Table 2 and Figure 1) and ~50% of that 15Nu was accumulated in the leaves and 50% in the stems (Figure 2). The maximum 15Nrec within the 6-day period following 15N application was observed near the critical period bracketing silking; at the V16 stage plants recovered almost half of the 15N applied (~0.45 kg kg−1) for both N rates, and ~50% of that 15Nu accumulated in the leaves followed by 35% in the stems, 15% in the husks (Table 2 and Figure 2). Only a negligible amount was found in the ear shoots that were just beginning to enlarge (data not shown). Nitrogen rate had no impact on 15Np in leaves; however, the 112N treatment increased 15Np to stems, while decreasing 15Np to husks. At the ER1 stage (when ~50% of plants in the entire experiment had visible extruded silks), maize plants recovered almost 0.40 kg kg−1 of the 15N fertilizer applied, and of this 15N uptake about 48% was allocated to the leaves, 36% to the stems, and 15% to the husks (Table 2 and Figure 2). At this ER1 stage, plants that were stressed due to the lack of N applied at 0N partitioned significantly more 15Np to the leaf and husk components than those at 112N (Table 2 and Figure 2). Similar patterns of 15Np between V16 and ER1 stages was probably due to overlapping days of 15N application and biomass harvesting at the ACRE site. Reducing the time from 15N application to biomass harvesting lowered total 15Nu at ER1 as expected (Figure 2). By dividing total 15Nu by number of days that plants remained in the field (4–6) after 15N application, it was observed that there was a similar daily 15Nu of ~0.30 kg of 15N ha−1 day−1 between those two stages at ACRE location (data not shown).
Figure 1. Total aboveground plant 15N recovery (15Nrec, kg kg−115N applied) per growth stage at 0 and 112N in 2013. Values are the means of five hybrids and two locations in Indiana (ACRE and PPAC) in 2013. Total 15N applied per growth stage was equal to 3.2 kg 15N ha−1. Error bars represent the standard error of the means. Letters represent significant differences between growth stages at the 0.05 level (Tukey HSD).
Figure 2. Partitioning of 15N uptake (15Nu, kg ha−1) per plant component (stem, leaf, husk, and ear) and proportional allocation of 15N uptake (15Np, %) at 0 and 112N over time in 2013. Values are the means of five hybrids and two locations in Indiana (ACRE and PPAC) in 2013. Except for V14 and ER1 representing each location only PPAC or ACRE, respectively. Total 15N fertilizer applied per growth stage was equal to 3.2 kg of 15N ha−1. *Represents significant difference between plant components across N rates at the 0.05 level T-test (LSD). Letters represent significant differences across growth stages at the 0.05 level (Tukey HSD).
In 2013, at the end of R1 stage (i.e., when >90% of plants completed silking), the 15Nrec declined, relative to previous stages, to about 0.32 and 0.34 kg kg−1 of the 15N applied, respectively, for 0 and 112N (Table 2 and Figure 1). At the R2 stage, plants recovered similar amounts of 15N as the R1 stage, but more 15Nu was partitioned to the ears (35 and 27% at 0 and 112N respectively) at the R2 stage than at R1. Additionally, at the R2 stage leaves and ears received equal 15Np in the 0N treatment (~35%). In the R4 stage it was observed that, out of 0.25–0.33 kg kg−1 of the 15N recovered by plants (at 0 and 112N, respectively), about 58 and 51% of the accumulated 15Nu was allocated into the ears, followed by 26–33% to leaves (Table 2 and Figure 2). Although the R5 stage presented the lowest total 15Nu of the growing season, plant 15Nrec was about 0.18 or 0.28 kg kg−1 of the 15N applied, respectively, for 0 and 112N rates. Approximately 68% of this late-stage 15N uptake was allocated to the kernels, followed by ~17% to the leaves, 11% to stems, and 3% to the husks.
In 2014, 15Nrec was nearly constant across all growth stages from V15 to R5 stage (Figure 3). There were larger differences across N rates this year than in 2013. At the R1 stage mean 15Nrec was about 0.21 and 0.42 kg kg−1, at 0 and 112N respectively (Table 3 and Figure 3), and, when averaged over both N rates, ~40% of this 15Nrec accumulated in the leaves, 25% in the stems, 11% in the husks, and a significant amount of 23% was allocated to the ears (Figure 4). At the R2 stage, plants recovered about 0.24 and 0.44 kg kg−1 of the 15N applied, at 0 and 112N respectively, allocating 31% to the leaves, 17% to the stems, and 44% to the ears on average of both N rates (Table 3 and Figure 4). Even with a smaller 15N fertilizer rate applied in 2014, plants recovered ~5% more of 15Nu at R2 in 2014 than in 2013, and a smaller proportion of 15Nu was partitioned to the vegetative components while higher amounts of 15Np were found in ears in 2014 than in 2013. In contrast to 2013, at the R4 stage in 2014 15Nu allocation to the ears (60%) was not statistically different for both N treatments. At R5, lower 15Np was found in the ears (62%) at both N rates in 2014 than in 2013.
Figure 3. Total aboveground plant 15N recovery (15Nrec, kg kg−115N applied) per growth stage at 0 and 112N in 2014. Values are the means of five hybrids and one location in Indiana (ACRE). Total 15N applied per growth stage was equal to 2.12 kg 15N ha−1. Error bars represent the standard error of the means. Letters represent significant differences between growth stages at the 0.05 level (Tukey HSD).
Table 3. Effects of overall sidedress N rate (0 and 112 kg N ha−1) on 15N fertilizer recovery (15Nrec, kg kg−115N applied) and proportional allocation of 15N uptake (15Np, %) in plant components at multiple development stages in 2014.
Figure 4. Partitioning of 15N uptake (15Nu, kg ha−1) per plant component (stem, leaf, husk, and ear) and proportional allocation of 15N uptake (15Np, %) at 0 and 112N over time in 2014. Values are the means of five hybrids and one location in Indiana (ACRE). Total 15N fertilizer applied per growth stage was equal to 2.12 kg of 15N ha−1. *Represents significant difference between plant components across N rate at the 0.05 level T-test (LSD). Letters represent significant differences across growth stages at the 0.05 level (Tukey HSD).
For the N rate factor, plant 15Nrec was almost always N rate dependent (Tables 2, 3). Furthermore, results suggest a strong influence on both 15N accumulation and the 15N allocations to the sink and source organs. A higher N rate generally helped increase 15Nu by maize plants at all development stages. Even though leaves accumulated higher amounts of 15N at the moderate N rate of 112N, before the flowering period in 2013 it was observed that 15Np to leaf components was higher with 0N treatment than at 112N. Conversely, both 15Np and 15Nu in stem components were higher at 112N than with the 0N. In the same year, plants with more N deficiency stress allocated more 15N to their reproductive ear organs at 0N than with the 112N rate. Except for the R5 stage, 15Np to the ears (cobs + kernels) following 15N application during reproductive stages was higher at 112N. In 2014, 15Nu accumulation as well as the proportional allocation of 15Nu to leaves was higher at 112N than at 0N for the entire season. In contrast to the previous year, the allocation of the 15Nu to the stems, husks and cobs was significantly larger at the 0N than with the 112N during the reproductive period.
Total 15N Uptake Rate Per GDD and Total Aboveground Dry Matter Accumulation Per GDD
Through the “multi-stage pulse labeling” 15N method it was also possible to estimate the extent to which plants uptake N during the growing season (Figures 5, 6). The amount of 15N taken up per daily cumulative GDD was calculated by dividing total plant 15Nu (kg ha−1) (Tables S1, S2) by the GDD °C (base 10°C) units accumulated from 15N application to plant biomass harvest (Table 1).
Figure 5. Effects of growth stages on total 15N uptake rate per GDD (15Nu, kg ha−1 GDD °C−1) at two N rates 0 and 112 kg N ha−1. Values are the means of two locations (ACRE13 and PPAC) in 2013. Total 15N applied per growth stage was equal to 3.2 kg 15N ha−1. Error bars represent the standard error of the means. Letters represent significant differences between growth stages at the 0.05 level (Tukey HSD).
Figure 6. Effects of growth stages on total 15N uptake rate per GDD (15Nu, kg ha−1 GDD °C−1) at two N rates 0 and 112 kg N ha−1. Values are the means of one location (ACRE) in 2014. Total 15N applied per growth stage was equal to 2.12 kg 15N ha−1. Error bars represent the standard error of the means. Letters represent significant differences between growth stages at the 0.05 level (Tukey HSD).
In 2013, the 112N treatment increased amounts of 15Nu per GDD after the flowering period, from R2 to R5 stages, relative to the 0N treatment. However, in 2014, 15Nu per GDD was higher at the 112N rate in comparison to the 0N rate for all stages.
For both years, the period of greatest 15Nu per GDD occurred from V14 to R4 stage, and the lowest 15Nu per GDD was observed at the R5 stage for both N rates (Figures 5, 6). In 2013, although plants showed higher 15Nrec at the V16 relative to the reproductive stages (R2, R4, and R5) (Figure 1), the actual rate of 15Nu accumulation per GDD was similar across stages, except that 15Nu per GDD was much lower at the R5 stage (Figure 5).
For both years, the 112N rate significantly increased dry matter accumulation (DM) per GDD in comparison to the 0N treatment (Figures 7, 8). There was no interaction between growth stages and N rate (data not shown). In 2013, the highest rate of DM accumulation was observed from V6 to the R4 stage with about 20 kg ha−1 per GDD on average of both N rates (0 and 112N) (Figure 7). In 2014, rates of DM accumulation per GDD were relatively consistent from V6 to R5 at ~20 kg ha−1 per GDD °C (Figure 8), or for a longer period relative to 2013. Additionally, the apparent lower DM accumulation from R2 to R4 stage could be explained by the lower GDD accumulation during that period.
Figure 7. Effects of growth stage intervals on total DM accumulation per GDD (DM, kg ha−1 GDD °C−1) average of two N rates (0 and 112 kg N ha−1). Numbers in parentheses are the daily cumulative GDD per growth interval (GDD °C day−1). Values are the combined means of two locations (ACRE and PPAC) in 2013. Error bars represent the standard error of the means. Letters represent significant differences between growth stage intervals at the 0.05 level (Tukey HSD).
Figure 8. Effects of growth stage intervals on total DM accumulation per GDD (DM, kg ha−1 GDD °C−1) average of two N rates (0 and 112 kg N ha−1). Numbers in parentheses are the daily cumulative GDD per growth interval (GDD °C day−1). Values are the means of one location (ACRE) in 2014. Error bars represent the standard error of the means. Letters represent significant differences between growth stage intervals at the 0.05 level (Tukey HSD).
Our results suggest that plant N uptake is strongly driven by both the sink strength and source availability. In 2013, plants took up similar rates of the 15N applied up to the R4 stage and DM accumulation began to decrease from R4 to R5 stages. In 2014, with possibly more acute soil N deficiency, 15Nu per GDD decreased even with a continual DM accumulation up to the R5 stage.
Fate of Labeled N Taken Up at Silking and Its Final Allocation at Physiological Maturity
In the “long term” labeling method implemented in 2014, the amount of 15Nu accumulated during the grain filling period of the 15N applied at silking (R1) was measured at physiological maturity in the duplicate R1R6 microplots (Table 3 and Figures 9, 10). At maturity 15Nrec was dramatically affected by the sidedressed N treatments; plants recovered ~0.45 and 0.72 kg kg−1 of the 15N fertilizer applied at R1, respectively, for 0 and 112N (Figure 9). Although the proportional allocation of the 15Nu (15Np) in the whole ear (cob + kernel) was not affected by the N rate, more 15Np to the kernels was observed with 112N (Table 3 and Figure 10). Overall, averaged over both N rates, ~77% of that accumulated 15N uptake at flowering (from the 15N application at R1) was partitioned to the kernels at maturity, followed by ~7% to the leaves, stems, and cobs and 4% to the husks. In the pulse labeled method, total plant 15Nrec averaged ~0.21 and 0.42 kg kg−1 of the 15N applied in the ~5 day period at R1; thus the remaining 15Nrec (0.24 and 0.30 kg kg−1) that was found at R1R6 was the 15Nu that was recovered during the rest of the grain filling period.
Figure 9. Total aboveground plant 15N recovery (15Nrec) at the R1R6 growth stage at 0 and 112N in 2014. Values are the means of five hybrids and one location in Indiana (ACRE). Total 15N applied per growth stage was equal to 2.12 kg 15N ha−1. Error bars represent the standard error of the means. *Represents a significant difference across N rates at the 0.05 level T-test (LSD).
Figure 10. Partitioning of 15N uptake (15Nu, kg ha−1) per plant component (stem, leaf, husk, and ear) and proportional allocation of 15N uptake (15Np, %) at 0 and 112N at the R1R6 stage in 2014. Values are the means of five hybrids and one location in Indiana (ACRE). Total 15N fertilizer applied per growth stage was equal to 2.12 kg of 15N ha−1. *Represents a significant difference of ear 15Np across N rate at the 0.05 level T-test (LSD).
We also used the net change approach to calculate the post-flowering 15N uptake per unit area (post15N, kg ha−1) as (Total plant 15Nu at R6 minus Total plant 15Nu at R1), to estimate remobilized 15Nu per unit area (15N rem, kg ha−1) as (Total plant 15Nu at R1 minus Stover 15Nu at R6), and to estimate the 15N remobilization efficiency (15NremEF, kg kg−1) as (Total plant 15Nu at R1 minus Stover 15Nu at R6 divided by Total plant 15Nu at R1) (Table 4). Stover was the sum of all plant components except kernels. Because we estimated the 15Nu at the silking period (multi-stage pulse labeling method) (Table 3), it was also possible to differentiate the proportion of the 15Nu in the grain that originated from the 15Nu accumulated during the reproductive period [i.e., post-flowering 15N uptake (post15N)], to the amount of 15Nu remobilized (15Nrem) from silking to maturity (Table 4).
Table 4. Consequences of overall N rate on 15N dynamics at grain maturity (R6 stage) following 15N application at the silking period (R1 stage) in 2014.
Pan et al. (1986) calculated Post15N uptake (Post15N) in maize using the net change approach; however, they sampled plants at the R1 stage on the same day that 15N was applied for the “long term” evaluation. Therefore, it might not be possible to estimate post15N and 15Nrem by the net change approach in our study since we applied 15N at the same time for both pulse labeled and long term approaches. Thus, from this method of calculation, the post15N represented about 54–42% of the total 15Nu accumulated from R1 to R6, and the 15Nrem accounted for 22–36% of the total 15Nu at R6, respectively, for 0 and 112N treatments. For both sources of 15Nu to the grain it was found that post15N contributed more at 0N than with 112N. Post 15N averaged ~70 and 53% of the grain 15Nu at 0 and 112N, respectively, while 15N rem contributed ~30–47% to grain 15Nu at R6, respectively, for 0 and 112N. The higher post15NuGr at 0N that was calculated by the net change from R1 to R6 may be explained by the larger effect of the N rate on the actual quantities of 15Nu.
Use of “Multi-Stage Pulse Labeling” 15N Technique to Determine the Fate of the Recent 15N Uptake into Plant Components
In contrast to most studies using labeled N, this research describes the potential implementation of the “multi-stage pulse labeling” 15N application in field-grown maize; this approach involves tracking 15N uptake and identifying changes in the proportional allocation of the recent 15N taken up into plant components at multiple maize development stages when field maize is grown under limited N supply. It is well-known that plants with superior ability to remobilize N to sink components yield more (Cliquet et al., 1990a), and that the higher ability of modern hybrids to continue accumulating N after the flowering period with efficient N allocation to its sink organs has increased both NUE and yield (Ciampitti et al., 2013). For these reasons, we regard the use of the labeled N technique to accurately estimate the fate of the recent N uptake in specific plant components over time, and especially during reproductive growth stages, to be a potentially powerful phenotyping tool.
For this 2-year study using the “multi-stage pulse labeling” technique, and after combining means of five hybrids, at two N rates (0 and 112N), it was observed that prior to flowering the distribution of short-term 15N uptake to leaves and stems was not affected by the N rate treatment as much as the actual amount of 15N uptake (except for V16 stage in 2013). Paponov and Engels (2005) also observed no effect of N rates (873 and 2,626 mg of N per plant) on 15N allocation to plant components during the vegetative and reproductive period. In our study, at late vegetative stages, maize plants primarily allocated 15Nu to their main sinks such as leaves and stems followed later by the husks and immature earshoots. Approximately 50% of the total 15N uptake was allocated to leaves for both N rates (0 and 112N). In contrast to our results, Cliquet et al. (1990b) analyzed N allocation in the below and aboveground maize plant components (8 days after 15N application) and found that stems were a stronger sink component than all plant organs at late vegetative stage; out of the total plant new N uptake, stems accumulated ~45% of the recent N taken up, followed by the leaf blades (22%) and roots (18%). For both N rates in our study, stems and husks appeared to act as transitional destinations of the 15Nu, and ears became an increasingly powerful sink component. Stems appeared to largely function as a “pipeline” where the 15N passed through accumulating primarily into the leaves from V14 to ~R2 stage at first and then to the ears at ~R2 or R4 and R5 stages. Ma et al. (1998) also confirmed the movement of the 15N to major sink component around the silking period. In their study, 3 days after 15N application (by stem infusion) at silking, plants accumulated ~40% of the new N in reproductive organs vs. 30% to the leaves and ~18% to the stems. However, in contrast to our plant partitioning procedure, they included husks as a reproductive organ (which may have resulted in their higher allocation values to reproductive plant parts at R1).
During the grain filling period, the moderate N rate (112N) in our study helped to maintain N contents in maize plants longer in the season, especially in the leaf components, while the 0N treatment showed the dynamics of 15Np in plant components deficient in N. Thus, our moderate N supply (112N) increased absolute amounts of 15Nu in the plants but it did not necessarily increase the allocation of 15Nu to the sink organs. Maize plants increased the allocation of the most recent 15Nu to the developing organs under the 0N rate compared to the 112N rate. Likewise, Paponov and Engels (2005) acknowledged the sink power of the ears over the vegetative plant parts at low N rates.
In agreement with our study, Rajcan and Tollenaar (1999) examined the effects of source: sink ratio on the N uptake during the reproductive period and demonstrated that after flowering the accumulation of N gradually declined in the vegetative plant parts with more severe declines occurring at lower N rates, while N accumulation increased into the ears. Our results show that the primary destination of the most recent 15Nu at both N rates was similar during the vegetative period but higher to the leaf component at R1 (~42%) and R2 stages (~38%), then increased to the ears at R4 (~55%), and kernels at R5 stages (~70%). The latter overall pattern of N partitioning in the plant was also confirmed by other authors using different labeling approaches. Crawford et al. (1982) provided 15N to the plants by nutrient solution just once at the time of pollen shed followed by multiple harvests. Yoneyama et al. (2003) used a leaf feeding labeling method to assess N assimilation in several crops including maize. According to the latter research, in general, N transportation within plants is driven by the two elements of growing plant components and reserve organs, and once N is absorbed in the reserve organs and mature leaves it will probably be reallocated to a developing organ.
The main N treatment influenced the amounts of plant 15N uptake and allocation of 15N at all growth stages; however, it did not alter the overall vegetative distribution of the 15Nu within the plants, and leaves always accumulated more 15N than stems. The multi-stage pulse labeling N approach measured the distribution of the recent N uptake by the plants, demonstrating that even at late grain fill plants were actively taking up N and mainly allocating the most recent N to the ears and/or kernels. Therefore, comparison between N treatments in this study was crucial to observe plant response to the N supply until late in the season and also to examine how plants change allocation patterns when they are source limited or stressed due to lack of N. For both N treatments, the largest amount of 15Nu and 15Nrec within 4–6 days of application was observed at the V16 stage and the lowest at the R5 stage. However, rates of 15N uptake per GDD°C were fairly consistent across stages, and only declined at R5. Correspondingly, 15Nrec was as high as 0.46 (kg kg−1) at 112N to as low as 18 (kg kg−1) at 0N. The latter suggests that further investigations in late vegetative stage or late reproductive stages may help to identify plant processes related to its N uptake efficiency.
Our results indicate that plant organs with lower Nc showed higher 15Np to the same component, and therefore 15Np in the plants operated as a signal of plant organs with the most deficiency in Nc. Paponov and Engels (2005) also found the same tradeoff with higher allocation of the 15Np to grain with reduced grain Nc under low N supply. Detailed evaluations of the proportional allocation of the 15Nu in plant components in our study suggested the late vegetative stage (V15) is a good stage to indicate relative plant total N status and variation in genotypic responses (data not shown). Indeed, 15Np to the leaves at V15 was strongly correlated with the responsiveness of 15Np to N rate treatments for the same plant component later in the season (data not shown).
In summary, for both years, leaves were the main sink component up to the kernel set period (R2) and, from R4 stage onward, more than 50% of the 15N absorbed accumulated in the ears. However, lower soil N availability and, consequently, higher N stress in 2014 prompted higher partitioning of 15Np to ears already at the R2 stage. Plant components that demonstrated higher 15Np had more apparent N deficiency.
Implementation and Considerations of 15N Use in Maize Field Experiments
There are many important factors that should be taken into account for a successful field study using a labeled N approach. Adequate plant sample size in field experiments using labeled N is crucial to identify genotypic responses for 15Nu partitioning to plant components. In a greenhouse study, Pan et al. (1986) addressed the large plant biomass variation in total N concentration or N content determination in plant components of five plants harvested at silking. Our studies were based in the field, and microplots were stringently selected for equal-density representations from the middle 3 or 4 plants. In 2014, after increasing sample size to four plants per growth stage and decreasing the 15N fertilizer rate by a third, we were still able demonstrate efficient assessment of 15Nu dynamics in the plants. Although harvesting more plants is even more preferable, resource costs of labeled N application and analyses are significant barriers.
The experiment layout and distance of plant zones that are used to analyze labeled or unlabeled treatments is critical (Van Cleemput et al., 2008). Several authors have suggested the use of plastic films, chambers, or tarpaulin materials to avoid problems with 15N percolation and with cross-contamination among treatments (Anhar, 2005). In the present study, we tried to represent the natural soil conditions in open environments using microplots within the main plots with no physical barrier. Microplots were installed at least 1 m from each other; the distance and the very low rates of 15N we utilized helped to circumvent contamination of plant zones by labeled N from elsewhere, and also helped prevent growth rate distortions arising from the extra sunlight in remaining plants harvested for biomass in other zones at a later time. The labeled N was carefully applied into the soil close to the root zones (~15 cm) to avoid direct contact of the 15N applied with the plant canopy and that, plus the water application, ensured immediate access of the labeled nitrate to growing roots.
The source of the labeled fertilizer used in the present study [Ca (15NO3)2 98 atom%15N] was intended to make nitrate-N (-N) immediately available to the plants and to minimize possible interactions of ammonium N sources (-N) to the negatively charged soil colloids. Roots can assimilate N by different mechanisms depending on the N source (Yoneyama et al., 2003), and several studies have observed that maize does not have a preferential uptake between mineral N sources (/-N) (Reddy and Reddy, 1993; Crozier et al., 1998). However, Pan et al. (1986) found higher accumulation of the labeled (NH4)2SO4 surface applied at planting in the grain when compared to the labeled Ca(15NO3)2 source which was applied 26 and 52 days after planting. Our study assumed that all labeled N taken up by plants was distributed to a uniform depth on both sides of the maize plants with the expectation that N flux into plant components over time should be consistent for the 3 or 4 consecutive plants that were removed from the 5- to 6-plant micro-plot zones.
The amount of labeled fertilizer applied per plant [0.16 g (2014) or 0.24 g (2013)] was enough to enable instrument detection of the 15N into all plant components, and yet not so much that it changed the dominant soil plus fertilizer N supply pool that met the majority of the maize plant N requirements. The “spike” of the highly enriched labeled fertilizer injected into the soil allowed a precise 15N determination in the mass spectrometer. The additional water after each time of 15N application was crucial to increase the N availability to plants via mass flow while reducing the risk of soil-mineral interactions from the fertilizer 15N applied with soil microorganisms. Labeled N (15NO−3) losses via leaching were expected to be minimal with the 1.6 l of water applied unless there was an intensive rain event between the pulse application and biomass harvest.
In the isotope laboratory, highly enriched plant materials required the use of additional “blank” and standard samples to monitor and avoid 15N carryover between labeled materials during the mass spectrometer analysis. Thus, lower quantities of labeled fertilizer were applied in 2014 because we realized after the first year that extremely enriched plant samples increased the risk of subsequent sample contamination during testing. High costs (i.e., >$35,000 with the labeled fertilizer alone in the 2-year period) made it even more imperative to utilize a lower practical rate of 15N in the second year.
The natural abundance of isotopic N in the atmospheric dinitrogen (N2) is considered equal to 0.366 atom %15N in a ratio (14N:15N) equal to 273:1 (Boaretto et al., 2004; Schepers and Raun, 2008). Similar values were found in the present study by sampling control plants in all hybrid and N rate treatments at the R4–R5 growth stages. These “control” values were used to calculate atom %15N excess (At%15NExcs) in plant components for each labeled sample as the amount of 15N which exceeded that represented in the non-labeled plant tissues originating from the natural abundance in the soil. The At%15NExcs was estimated by subtracting atom%15N (At15N) of each labeled plant component by the N natural abundance found in the present study (0.3683 or 0.3677 Atom %15N) (Equation a). In this work, amounts of pre-existing 15N in the soil based on natural abundance were neglected; however, the latter approach confirmed that soils were not labeled before the current field experiment. Additionally, discrimination from the uptake of natural abundant 15N to the 15N fertilizer was not important to estimate in our evaluations of short-term 15N accumulation in the plants because the large 98% enrichment of the fertilizer 15N applied provided confident identification of the isotope N taken up by plants.
Perspectives on Labeled 15N Approaches in Field Maize Research
The use of isotopic N has been acknowledged for its accuracy as N tracer for several decades (Van Cleemput et al., 2008). The multifaceted use of 15N in open systems is, when properly conducted, highly effective for N balance investigations in real, as distinct from simulated, soil-plant systems under different G × E × M interaction treatments. In this work, through mass spectrometry analysis, we estimated the δ15N in aboveground plant components to calculate uptake, partitioning, allocation, and recovery of N throughout the season. Observations of high correlations between 15N and total N allowed for detailed evaluations of maize plant response to labeled N additions in the root zone at multiple stages of development and at different N rates. After decades evaluating variations of plants δ15N in both laboratory and field experiments using 15N feeding methods, Yoneyama et al. (2003) also showed a strong association of δ15N with N dynamics in several crops.
Several techniques have been implemented using 15N in either open or controlled environments (field, greenhouses, and chambers) during short- or long-term methods of 15N application (Duete et al., 2009; Ning et al., 2017). However, few studies have invested in field studies using highly enriched 15N to address the allocation of N in plant components spanning the reproductive period or other intervals of the growing season in maize.
The “long term” labeled N application has occasionally been implemented as a one-time alternative to differentiate and more precisely estimate post-flowering N uptake and N remobilization from late vegetative stage or silking periods to physiological maturity. This method has helped to quantify N flux within the plants over the period of interest. Mae and Ohira (1981) and Ta and Weiland (1992) calculated the 15N remobilization by subtracting the 15N uptake in the stover at R6 by the 15N in the stover at R1. However, 15N remobilization can also be estimated as the difference between the 15N accumulated in the stover at maturity minus the whole plant 15N uptake at silking (Pan et al., 1986). The latter calculation was used in the present study.
When we compared 15Nu accumulated in the R1R6 zones by the more traditional “long term” evaluation method of 15N application vs. the results of the Pan et al. (1986) study, we observed a much higher level (~77%) of 15Nu accumulated in the kernels at R6 vs. the lower allocation (~51%) of 15N uptake to the grain in the Pan et al. study (Table 3). This may be partially explained by their lower total plant 15N recovery (~75%) when compared to ours (~83%), by their use of older and different genotypes that preferentially accumulated 15N to the stalk and root components, and by their shorter post-flowering period of 33 days in comparison to our 66 days. Furthermore, Pan et al. (1986) found that post 15N accounted for 42% of the total 15Nu accumulated in the plant at physiological maturity, while we found that post 15N contributed 54 and 42% to the total 15Nu accumulated at R6, for 0 and 112N treatments (Table 4). Similarly to our results, the review study of Ciampitti and Vyn (2013) demonstrated that post-flowering N uptake can contribute up to about 56% of the total plant N uptake at R6 stage.
Results from Ma et al. (1998) generally agreed with our current findings for the “long term” 15N accumulation. While we found that 77% of 15Nu in the kernels was taken up from R1 to R6 stage, their results from a one-time 15N application show 65–75% of infused 15N was accumulated in the kernels. However, they estimated 59 and 82% as the amount of remobilized N going to the grain from vegetative components, and a much lower contribution of the post-silking 15N uptake to the total 15Nu at R6 (ranging from 18 to 41%, which was similar to our 42% in the 112N treatment). Bertin and Gallais (2000) indicated a range of 35–65% of the N accumulated in the grain originated from N remobilization. We observed that remobilized 15N accounted for 22 and 36% of the total 15Nu accumulated at R6, for 0 and 112N rates (Table 4).
Furthermore, in attempting to integrate the isotope technique and breeding approaches, Gallais et al. (2007) estimated the contribution of the N remobilization and post N uptake in several maize inbred lines using a one-time 15N application. In their 3-year experiment, comparing different methods of labeled N application, they found that ~83% of that 15Nu was accumulated into the kernels at R6, slightly higher than our 15N proportional allocation ~77%. They also estimated ~25% of the total 15Nu at maturity originated from the post-flowering N uptake and about 62% was coming from the 15N remobilization. Our study found that post 15N accounted for ~50% of the total 15Nu accumulated at R6, and about 30% of that 15N accumulated at R6 originated from the remobilized 15N (Table 4). The difference in results of these two detailed studies may be explained by our 15N application at silking while their 15N application was at the beginning of stem elongation period (which enabled a better estimate of the accumulation in vegetative components), and by their use of 15N fertilizer with much lower enrichment to test an even larger number of genotypes.
This “long term” labeled approach has been broadly used to assess 15N uptake dynamics at multiple plant development stages (Friedrich and Schrader, 1979; Cliquet et al., 1990a,b; Deléens et al., 1994; Ma et al., 1998; Ning et al., 2017). When the influence of the 15N application at planting and at anthesis on improvements of N harvest index (NHI) in wheat (Triticum aestivum L.) was studied, Wuest and Cassman (1992) showed that 15N applied at flowering improved NHI (0.89) by almost 20%. The latter suggests a possibility that the 15N injection at the R1 stage overestimates the real amounts partitioned to the grain during the reproductive period.
Apart from injection method into the soil, Hertenberger and Wanek (2004) suggested that 15N stem infiltration and leaf feeding were the most effective labeling techniques to analyze N in thick-stem plant and grasses species, respectively. More recently, Putz et al. (2011) has advanced the termed “in-situ” labeling technique in plant species cultivated in greenhouse. Through the use of paint brush, they applied 2 mg of 15N Urea (98 atom% 15N) onto the leaves once a day during 5 days and plants were collected at short interval of hours after labeling.
Hence, in contrast to the techniques previously described by authors as “long term” and “short term” labeling approaches, our modified “multi-stage pulse labeling” technique 15N was used as a tracer to evaluate the fate of recently applied 15N into maize plants. The “multi-stage pulse labeling” technique used very small absolute amounts of highly enriched 15N fertilizer (98 atom% 15N) and different plants were routinely exposed to the labeled N for short intervals (just 4–6 days) mainly during the reproductive period. Our approach was able to accurately estimate plant N use in maize and to detect the dynamics of the in-season 15N fertilizer applications at development-stage-specific levels with a special focus on the reproductive period. This technique has high potential as an additional tool for phenotyping evaluation of genotypes with superior N use efficiency.
Although our emphasis was on applying 15N to maize plants during reproductive stages, this research observed a positive association of the allocation of 15Nu at late vegetative stage with the allocation of 15Nu during grain filling. Therefore, we propose that three N plant development stages (V15, R2, and early R5) are promising growth stages to evaluate potential hybrids with superior NUE. Future research studies using the multi-stage pulse labeling approach for phenotyping purposes with limited budgets should prioritize sampling a higher amount of plants (at least four) rather than increasing the rates of labeled N fertilizer beyond those used in our studies.
This research describes and evaluates the effectiveness of intentional multi-stage 15N applications that we term as a “multi-stage pulse labeling” method of 15N application to late vegetative and reproductive-stage maize in field experiments. It is a potentially superior phenotyping technique to identify hybrids superior in NUE because of improved post-silking DM and N gains compared to hybrids of earlier decades. One of the major benefits of the 15N technique is that it precisely differentiates N allocation in plant organs of diverse genotypes under various N management situations, and then our modified approach demonstrates the extent to which plants are still actively taking up N during reproductive growth stages in response to new soil N supply. Even at the R5 stage when minimal new DM is being accumulated, maize plants were still taking up N and mainly allocating it to the kernels (~70%).
Our results suggest that plant N uptake is strongly driven by the sink strength and source availability. Leaves were the main sink component up to the kernel set period (R2) and, from R4 stage onward, more than 50% of the 15N absorbed accumulated in the ears. The moderate N rate (112N) increased absolute 15N uptake at all growth stages relative to the 0N treatment. Plants more stressed due to N deficiency (0N) partitioned more 15Nu to the reproductive organs. The primary allocation of the new 15Nu at both N rates was similar to stems and leaves during the late vegetative period, then higher to the leaf component at R1 (~42%) and R2 stages (~38%), and then increased to the ears at R4 (~55%) and kernels at R5 stages (~70%). However, lower soil N availability and, consequently, higher N stress in 2014 prompted higher partitioning of 15Np to ears already at the R2 stage.
Hence, the multi-stage pulse labeling technique proved to be a reliable and precise phenotyping approach that seed companies and universities can implement to investigate genotypes with particular traits targeted at improved maize N use efficiency.
AOV conducted field experiments, laboratory measurements, and statistics and initial writing. TV conceived of and directed the project, the data interpretation, and its writing and editing, JC assisted with the writing, TC provided institutional support for the project via Dow AgroSciences (where he was a maize genetics research leader in the Trait Product Development Department at the time), and TF provided the Purdue Isotope Lab resources and technical assistance with the data output from the mass spectrometer analyses of isotopic plant tissue samples.
The research was financially supported by the Trait Product Development group of Dow AgroSciences, Inc., Indianapolis, IN (2013–2015) via grant numbers 16-45542727 and 16-45577925.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We express our thanks to the field technicians Terry West and Alicia West, and the staff at the Purdue University research farms for their help with the establishment and management of the field experiments. Thanks to David Gamblin for helping with the mass spectrometer analysis in the Isotope Lab in the Department of Earth and Atmosphere Sciences at Purdue University.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2017.01360/full#supplementary-material
Supplementary 15N Methodology Procedures
More details on the modified “multi-stage pulse labeling” approach can be viewed in Supplementary Section A. Supplementary data of individual hybrid responses can be viewed in Supplementary Section B.
15N, Isotopic Nitrogen; 15Np, Proportional allocation of 15N uptake in plant components; 15Nrec, 15N fertilizer recovery; 15Nu, 15N uptake per unit area; DM, Aboveground plant dry matter; GY, grain yield; Nc, Total Nitrogen concentration; 15Nrem, Remobilized per unit area; 15NremEF, 15N Remobilization efficiency; Nu, Total Nitrogen uptake per unit area; NUE, Nitrogen use efficiency; Post15N, Post silking 15N uptake per unit area.
Amado, T. J. C., Mielniczuk, J., and Aita, C. (2002). Recomendação de adubação nitrogenada para o milho no RS e SC adaptada ao uso de culturas de cobertura do solo, sob plantio direto. (In Portuguese, with English abstract). Rev. Brasil. Ciênc. Solo 26, 241–248. doi: 10.1590/S0100-06832002000100025
Boaretto, A. E., Trivelin, P. C. O., and Muraoka, T. (2004). “Uso de isotopos como tracadores em fertilidade do solo e nutricao de plantas, (In Portuguese),” in Proceedings of the Soil Fertility and Plant Nutrition Conference (FERTBIO), (Lages: Centro de Energia Nuclear na Agricultura, CENA).
Cassman, K. G., Dobermann, A., Walters, D. T., and Yang, H. (2003). Meeting cereal demand while protecting natural resources and improving environmental quality. Annu. Rev. Environ. Resour. 28, 315–358. doi: 10.1146/annurev.energy.28.040202.122858
Chen, K., Kumudini, S. V., Tollenaar, M., and Vyn, T. J. (2015). Plant biomass and nitrogen partitioning changes between silking and maturity in newer versus older maize hybrids. Field Crops Res. 183, 315–328. doi: 10.1016/j.fcr.2015.08.013
Ciampitti, I. A., and Vyn, T. J. (2011). A comprehensive study of plant density consequences on nitrogen uptake dynamics of maize plants from vegetative to reproductive stages. Field Crops Res. 121, 2–18. doi: 10.1016/j.fcr.2010.10.009
Ciampitti, I. A., and Vyn, T. J. (2012). Physiological perspectives of changes over time in maize yield dependency on nitrogen uptake and associated nitrogen efficiencies: a review. Field Crops Res. 133, 48–67. doi: 10.1016/j.fcr.2012.03.008
Ciampitti, I. A., Murrell, S. T., Tuinstra, M., Camberato, J. J., Xia, Y., Friedemann, P., et al. (2013). Physiological dynamics of maize nitrogen uptake and partitioning in response to plant density and N stress factors: I. Vegetative phase. Crop Sci. 53, 2588–2602. doi: 10.2135/cropsci2013.01.0041
Cliquet, J. B., Deléens, E., and Mariotti, A. (1990a). C and N mobilization from stalk and leaves during kernel filling by C and N tracing in Zea mays L. Plant Physiol. 94, 1547–1553. doi: 10.1104/pp.94.4.1547
Cliquet, J. B., Deléens, E., Bousser, A., Martin, M., Lescure, J. C., Prioul, J. L., et al. (1990b). Estimation of carbon and nitrogen allocation during stalk elongation by 13C and 15N tracing in Zea mays L. Plant Physiol. 92, 79–87. doi: 10.1104/pp.92.1.79
Coelho, A. M., Franca, G. C., Bahia, A. F. C., and Guedes, G. A. A. (1991). Balanço de nitrogênio 15N em latossolo vermelho-Escuro, sob vegetação de cerrado, cultivado com milho. (In Portuguese, with English abstract). Revi. Brasil. Ciênc. Solo 15, 187–193.
Coque, M., and Gallais, A. (2007). Genetic variation for nitrogen remobilization and postsilking nitrogen uptake in maize recombinant inbred lines: heritabilities and correlations among traits. Crop Sci. 47, 1787–1796. doi: 10.2135/cropsci2007.02.0096
Coque, M., Martin, A., Veyrieras, J. B., Hirel, B., and Gallais, A. (2008). Genetic variation for N-remobilization and postsilking N-uptake in a set of maize recombinant inbred lines. 3. QTL detection and coincidences. Theor. Appl. Genet. 117, 729–747. doi: 10.1007/s00122-008-0815-2
Crawford, T. W., Rendig, V. V., and Broadbent, F. E. (1982). Sources, fluxes, and sinks of nitrogen during early reproductive growth of maize (Zea mays L.). Plant Physiol. 70, 1654–1660. doi: 10.1104/pp.70.6.1654
Crozier, C. R., King, L. D., and Volk, R. J. (1998). Tracing nitrogen movement in corn production systems in the North Carolina Piedmont: a nitrogen-15 study. Agron. J. 90, 171–177. doi: 10.2134/agronj1998.00021962009000020009x
Deléens, E., Cliquet, J. B., and Prioul, J. L. (1994). Use of 13C and 15N plant label near natural abundance for monitoring carbon and nitrogen partitioning. Aust. J. Plant Physiol. 21, 133–146. doi: 10.1071/PP9940133
Dobermann, A., and Cassman, K. G. (2004). “Environmental dimensions of fertilizer nitrogen: what can be done to increase nitrogen use efficiency and ensure global food security?” in Agriculture and the Nitrogen Cycle: Assessing the Impacts of Fertilizer Use on Food Production and the Environment, eds A. R. Mosier, J. K. Syers, and J. R. Freney (Paris: SCOPE 65), 261–278.
Duete, R. R. C., Muraoka, T., Da Silva, E. C., Ambrosano, E. J., and Trivelin, P. C. O. (2009). Acúmulo de nitrogênio (15N) pelos grãos de milho em função da fonte nitrogenada em latossolo vermelho. (In Portuguese, with English abstract). Bragantia 68, 463–472. doi: 10.1590/S0006-87052009000200021
Duvick, D. N. (1984). “Genetic contributions to yield gains of U.S. hybrid maize, 1930 to 1980,” in Genetic Contributions to Yield Gains of Five Major Crop Plants. CSSA Special Publication No. 7, ed W. R. Fehr (Madison, WI: Crop Science Society of America and American Society of Agronomy), 15–47.
Friedrich, J. W., and Schrader, L. E. (1979). N deprivation in maize during grain-filling. II. Remobilization of 15N and 35S and the relationship between N and S accumulation1. Agron. J. 71, 466–472. doi: 10.2134/agronj1979.00021962007100030021x
Gallais, A., Coque, M., Le Gouis, J., Prioul, J. L., Hirel, B., and Quilléré, I. (2007). Estimating the proportion of nitrogen remobilization and of postsilking nitrogen uptake allocated to maize kernels by nitrogen-15 labeling. Crop Sci. 47, 685–693. doi: 10.2135/cropsci2006.08.0523
Gallais, A., Coque, M., Quilléré, I., Prioul, J. L., and Hirel, B. (2006). Modelling postsilking nitrogen fluxes in maize (Zea mays) using 15N-labelling field experiments. New Phytol. 172, 696–707. doi: 10.1111/j.1469-8137.2006.01890.x
Guo, M., Rupe, M. A., Wei, J., Winkler, C., Goncalves-Butruille, M., Weers, B. P., et al. (2014). Maize ARGOS1 (ZAR1) transgenic alleles increase hybrid maize yield. J. Exp. Bot. 65, 249–260. doi: 10.1093/jxb/ert370
Hertenberger, G., and Wanek, W. (2004). Evaluation of methods to measure differential 15N labeling of soil and root N pools for studies of root exudation. Rapid Commun. Mass Spectrom. 18, 2415–2425. doi: 10.1002/rcm.1615
Ladha, J. K., Pathak, H., Krupnik, T. J., Six, J., and Kessel, C. V. (2005). Efficiency of fertilizer nitrogen in cereal production: retrospects and Prospects. Adv. Agron. 87, 85–176. doi: 10.1016/S0065-2113(05)87003-8
Ma, B. L., Dwyer, L. M., Tollenaar, M., and Smith, D. L. (1998). Stem infusion of nitrogen-15 to quantify nitrogen remobilization in maize. Commun. Soil Sci. Plant Anal. 29, 305–307. doi: 10.1080/00103629809369947
Moll, R. H., Kamprath, E. J., and Jackson, W. A. (1982). Analysis and interpretation of factors which contribute to efficiency of nitrogen utilization. Agron. J. 74, 562–564. doi: 10.2134/agronj1982.00021962007400030037x
Mosier, A. R., Syers, J. K., and Freney, J. R. (eds.) (2004). Agriculture and the Nitrogen Cycle: Assessing the Impacts of Fertilizer Use on Food Production and the Environment. Washington, DC: Island Press.
Mueller, N. D., West, P. C., Gerber, J. S., MacDonald, G. K., Polasky, S., and Foley, J. A. (2014). A tradeoff frontier for global nitrogen use and cereal production. Environ. Res. Lett. 9:054002. doi: 10.1088/1748-9326/9/5/054002
Ning, P., Fritschi, F. B., and Li, C. (2017). Temporal dynamics of post-silking nitrogen fluxes and their effects on grain yield in maize under low to high nitrogen inputs. Field Crops Res. 204, 249–259. doi: 10.1016/j.fcr.2017.01.022
O'Brien, J. M., Hamilton, S. K., Podzikowski, L., and Ostrom, N. (2012). The fate of assimilated nitrogen in streams: an in situ benthic chamber study. Freshw. Biol. 57, 1113–1125. doi: 10.1111/j.1365-2427.2012.02770.x
Pan, W. L., Camberato, J. J., Jackson, W. A., and Moll, R. H. (1986). Utilization of previously accumulated and concurrently absorbed nitrogen during reproductive growth in maize: influence of prolificacy and nitrogen source. Plant Physiol. 82, 247–253. doi: 10.1104/pp.82.1.247
Portela, S. I., Andriulo, A. E., Sasal, M. C., Mary, B., and Jobbágy, E. G. (2006). Fertilizer vs. organic matter contributions to nitrogen leaching in cropping systems of the Pampas: 15N application in field lysimeters. Plant Soil 289, 265–277. doi: 10.1007/s11104-006-9134-z
Putz, B., Drapela, T., Wanek, W., Schmidt, O., Frank, T., and Zaller, J. G. (2011). A simple method for in situ-labelling with 15N and 13C of grassland plant species by foliar brushing. Methods Ecol. Evol. 2, 326–332. doi: 10.1111/j.2041-210X.2010.00072.x
R Development Core Team (2005). R: A Language and Environment for Statistical computing. Vienna: R Foundation for Statistical Computing. Available online at: http://www.R-project.org
Schepers, J. S., and Raun, W. (2008). Nitrogen in Agricultural Systems. Madison, WI: American Society of Agronomy; Crop Science Society of America; and Soil Science Society of America; Agronomy Monographs.
Silva, E. C., Muraoka, T., Buzetti, S., and Trivelin, P. C. O. (2006). Nitrogen management in corn under no-tillage with different cover crops in a Rhodic Hapludox soil. (In Portuguese, with English abstract). Pesquisa Agropecuaria Brasileira (EMBRAPA) Brasilia 41, 477–486. doi: 10.1590/S0100-204X2006000300015
USDA-ERS (2013). United States Department of Agriculture. Economic Research Service. Available online at: https://www.ers.usda.gov/data-products/fertilizer-use-and-price/aspx. (Accessed Feb 20,2017).
Wuest, S. B., and Cassman, K. G. (1992). Fertilizer-nitrogen use efficiency of irrigated wheat: II. partitioning efficiency of preplant versus late-season application. Agron. J. 84, 689–694. doi: 10.2134/agronj1992.00021962008400040029x
Keywords: isotopic N, maize hybrids, plant N uptake, plant component N distribution, N recovery efficiency
Citation: de Oliveira Silva A, Camberato JJ, Coram T, Filley T and Vyn TJ (2017) Applicability of a “Multi-Stage Pulse Labeling” 15N Approach to Phenotype N Dynamics in Maize Plant Components during the Growing Season. Front. Plant Sci. 8:1360. doi: 10.3389/fpls.2017.01360
Received: 16 January 2017; Accepted: 20 July 2017;
Published: 04 August 2017.
Edited by:Marta Wilton Vasconcelos, Universidade Católica Portuguesa, Portugal
Reviewed by:Lixing Yuan, China Agricultural University, China
Jan Kofod Schjoerring, University of Copenhagen, Denmark
Copyright © 2017 de Oliveira Silva, Camberato, Coram, Filley and Vyn. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Tony J. Vyn, email@example.com