2.1 Soil samples
In this study, the two rice genotypes Wuyujing 3 (WYJ3) and Wuyunjing 7 (WYJ7), which had been reported to have a weak and strong BNI exudation capacity, respectively (Sun et al. 2016) were selected. Four paddy soil samples with different pH from three provinces in China, including two acidic soils in Jiangxi (JX) and Fujian (FJ) and two alkaline soils in Sichuan (SC1 and SC2) were collected, after harvesting in October 2021. Soil samples were taken from the cultivated layer at a depth of 0–20 cm at three randomly selected points at each site. After being brought back to the laboratory, impurities such as large gravel particles and dead branches were first removed, then passed through a 2 mm sieve and mixed thoroughly. These samples were stored at 4°C and used for pot experiments as well as soil property analysis. Soil properties are detailed in Table 1.
Table 1
Soil properties of the four soils used in this study (average ± standard deviation)
|
pH
|
C/N
|
SOC
|
TN
|
DOC
|
NH4+
|
NO3−
|
|
g/kg
|
mg/kg
|
mg N/kg
|
JX
|
5.09 ± 0.09c
|
10.01 ± 0.47bc
|
14.10 ± 0.82c
|
1.38 ± 0.02b
|
74.01 ± 0.47b
|
10.03 ± 0.38a
|
4.48 ± 0.10c
|
FJ
|
6.00 ± 0.12b
|
20.35 ± 0.24a
|
23.80 ± 0.41a
|
1.18 ± 0.01c
|
81.63 ± 3.68ab
|
1.18 ± 0.40b
|
5.33 ± 0.43c
|
SC1
|
7.96 ± 0.02a
|
10.57 ± 0.23b
|
12.93 ± 0.39c
|
1.21 ± 0.01c
|
69.09 ± 3.65bc
|
1.08 ± 0.16b
|
9.52 ± 0.83b
|
SC2
|
7.94 ± 0.04a
|
9.67 ± 0.06c
|
17.91 ± 0.10b
|
1.84 ± 0.01a
|
87.39 ± 3.64a
|
1.02 ± 0.25b
|
13.95 ± 0.92a
|
JX is a paddy soil collected from Jiangxi, China; FJ is a paddy soil collected from Fujian, China; and SC1 and SC2 are paddy soils collected from Sichuan, China. C/N, SOC: TN ratio; SOC, soil organic carbon; TN, total nitrogen; DOC, dissolved organic carbon |
2.2 Plant materials and growth condition
Rice seeds (WYJ3 and WYJ7) were soaked in 10% H2O2 for 30 min, washed with deionized water, and subsequently placed in an incubator with two foam floating nets at 28°C for 4 days in the dark and after germination seedlings were planted in soil. Normal daylight conditions (light 12 h/dark 12 h, 28°C during the day and 22°C at night), prevailed in the growth chamber throughout the experimental period, with the light field open from 6 am to 6 pm. After 30 days, the seedlings were transplanted into a series of pots with one part used for the 15N urea labeling experiment and the other part for the 15N tracing experiment after the seedlings had grown to the tillering stage.
2.3 15N urea labeling experiment
In this experiment, the pots were divided into eight groups, i.e. WYJ3 and WYJ7 were transplanted into four soils (JX, FJ, SC1, SC2), respectively. After one month of growing rice in seedling pots, seedlings of uniform growth were selected for transplanting, and one seedling was transplanted per pot. Before transplanting, all soils equivalent to 250 g dry weight were weighed into pots and 5 mL of urea solution (5.28% atom of 15N urea) was added to each pot as a base fertilizer at a rate of 40 mg N kg− 1 dry soil. Each pot should be watered and cleared of aquatic plants every two days. Destructive sampling of each pot was carried out on day 35 after transplanting, the rice plants were collected and rinsed, then the soils in the pots were uniformly mixed and extracted, N content and 15N abundance were measured after soil samples had been dried at 60°C. N content, 15N abundance and biomass were determined after plants had been dried at 80°C.
2.4 15N tracing experiment
Each fresh soil sample equivalent to 250 g of dry soil weight was added to each pot before transplanting, followed by the application of (NH4)2SO4 as a base fertilizer (40 mg N kg− 1 soil) and about 2 cm of flooded water. Other managements were the same as in 2.3.
Appropriate modifications were made to the 15N tracing experiment designed by (He et al. 2020). Briefly, on day 30 after transplanting, two sets of 15N isotope markers were set, one for 15NH414NO3 (10.20% atom) and the other for 14NH415NO3 (10.12% atom). 40 mg N kg− 1 soil of 15NH414NO3 or 14NH415NO3 is uniformly added to the soil core by four-hole injection, destructive sampling was then carried out at 0.5, 24, 48, 72 hours to determine the soil NH4+ and NO3− concentration and 15N abundance of each treatment. Dissolved organic carbon (DOC), N content, 15N abundance and biomass of each dry plant were also measured.
2.5 Analytical methods
The DMP-2mV/pH detector (Quark Ltd, Nanjing, China) was used to determine soil pH at a ratio of 1:2.5 (w/v). Soil organic carbon (SOC) content was measured using high temperature exothermic digestion by H2SO4-K2Cr2O7. NH4+ and NO3− were extracted with 1 M KCl at a ratio of 1:5 (w/v) in a shaker for 60 min at 250 rpm under 25°C. DOC was extracted by ultrapure water at a ratio of 1:5 (w/v) and measured by Analyzer Multi N/C (Analytik Jena, Jena, Germany). Concentrations of NH4+ and NO3− were determined using a continuous flow analyser (SA1000, Skalar, Breda, Netherlands). NH4+ and NO3− were extracted using micro-diffusion and absorbed by filter paper with 1 M oxalic acid (Zhang et al. 2017), the 15N abundance of NH4+ and NO3− were then determined by stable isotope ratio mass spectrometry (IRMS 20–22, SerCon, Grewe, UK) after a thorough drying. N content and 15N abundance in plants and soils from pot experiment were determined using a Delta V advantage isotope ratio mass spectrometer (Thermo Scientific, Bremen, Germany).
DNA was extracted from 0.5 g fresh soil after destructive sampling which was according to the manufacturer's instructions (MoBio Laboratories, Inc., Carlsbad, CA, USA). Real-time fluorescent quantitative PCR was performed on biological trios by CFX-96TM Real-Time System (Bio-Rad Laboratories Inc., Hercules, CA, USA) and the abundance of AOA (Arch-amoAF/Arch-amoAR) and AOB (amoA1F/amoA2R) determined by primer sets and thermal conditions are listed in Table S1. Reaction mixtures for AOB and AOA amoA were set and standard curves were established according to (Huang et al. 2015), with amplification efficiencies of 96.8% and 99.4% for AOA and AOB, respectively. Amplification specificity was calculated using melting curve analysis, while copy number was calculated by oven drying basis.
Gross N transformation rates in paddy soils were determined by 15N tracing experiment in combination with the Ntraceplant tool (He et al. 2020) to quantify gross N transformation and plant N uptake rates. NH4+ and NO3− concentrations (mean ± standard deviation) and 15N atomic percent excesses (mean ± standard deviation) measured by the 15N tracing experiment are used as input variables to the tool with N transformation rates calculated by zero-order kinetics, first-order kinetics or Michaelis-Menten kinetics.
2.6 Calculation and statistical analysis
The following equations were used to calculate fertilizer NH4+ retention time (NH4+ Retention time), net NH4+ release rate (Net RNH4), NH4+ immobilization rate (INH4), N use efficiency (NUE) and N loss.
$${NH}_{4}^{+} Retention time=\frac{Fertilizer {NH}_{4}^{+}}{{O}_{NH4}+{I}_{NH4}}$$
$$Net {R}_{NH4}={R}_{NH4}-{A}_{NH4}$$
$${I}_{NH4}={I}_{NH4-Nrec}+{I}_{NH4-Nlab}$$
where NH4+ Retention time refers to fertilizer NH4+ retention time (d), Fertilizer NH4+ is the amount of NH4+ applied in the 15N tracing experiment (40 mg N kg− 1), ONH4 means the rate of oxidation from NH4+ to NO3− (mg N kg− 1 d− 1), INH4 is the rate of NH4+ immobilization (mg N kg− 1 d− 1), RNH4 is the release rate of adsorbed NH4+, ANH4 is the adsorption rate of free NH4+, INH4−Nrec is the immobilization rate of NH4+ to recalcitrant organic N, INH4−Nlab is the immobilization rate of NH4+ to labile organic N.
$$NUE \left(\%\right) =\frac{{C}_{N Plant}\times {A}_{N plant}\times {W}_{plant}}{{C}_{Nf}\times {A}_{Nf}}\times 100$$
$$Soil N Retention \left(\%\right)=\frac{{C}_{N soil}\times {A}_{N soil}\times {W}_{soil}}{{C}_{Nf}\times {A}_{Nf}}\times 100$$
$$N loss \left(\%\right)=100\%-NUE\left(\%\right)-Soil N Retention\left(\%\right)$$
where CN plant is the N content of the plant (g kg− 1), AN plant is the 15N atomic percent excess of N in the plant (% atom), Wplant is the dry weight of the plant (kg), CNf is the amount of the applied N fertilizer (g pot− 1), and ANf is the 15N atomic percent excess of applied N fertilizer (% atom). CN soil is soil N content (g kg− 1), AN soil is 15N atomic percent excess of soil N (% atom), Wsoil is soil dry weight (kg), CNf refers to the amount of the applied N fertilizer (g pot− 1), and ANf is the 15N atomic percent excess of applied N fertilizer (% atom). N loss is calculated by the 15N mass balance method (Liu et al. 2020).
Statistical analyses on all gross N transformation rates were performed by SigmaStat 3.4 (Systat Software, San Jose, CA, USA). One-way ANOVA combined with LSD test was used to assess differences in the same genotype grown in four soils, and t-test was used to determine the level of variation in the same soil after planting different genotypes. Briefly, the mean and standard deviation of each gross transformation rate, as well as actual experimental replicates, were used to test for differences in gross soil N transformation rates and plant N uptake rates. The rest of the soil property analyses were carried out in SPSS 25.0 (IBM crop., USA). Pearson correlation and nonlinear analysis were used to test the correlation between N fertilizer fates, soil gross N transformation rates and soil properties.