Soil samples
Two greenhouse soil samples were collected from Suzhou, Jiangsu Province, China. On one soil, tomato and strawberry were grown continuously for more than 10 years, with a high accumulation of NO3− (on average 82.83 mg N kg−1) (HS). The other was previously planted with tomato and cucumber and lay fallow for several years, resulting in a low NO3− concentration (on average 8.97 mg N kg−1) (LS). Four plots (1 m ⋅ 1 m) were randomly selected at each sampling site. The surface soil (0-20 cm) was collected and pooled together in equal amounts. Then soil samples were immediately passed through a 2 mm sieve and roots were removed. Each sample was divided into three sub-samples. One was air-dried to determine soil chemical properties (Table 1), one was stored in a refrigerator at -80 oC for the measurement of soil microbial properties, and another sub-sample was stored at 4 oC to perform 15N tracing studies.
Table 1
soil properties of initial HS and LS (means ± SD, n=3)
| pH | EC µS cm−1 | C/N | SOC g kg−1 | TN g kg−1 | DOC mg kg−1 | NH4+ mg kg−1 | NO3− mg kg−1 |
LS | 7.76±0.03a | 55.42±1.39b | 8.81±0.12a | 6.73±0.07b | 0.76±0.02b | 143.50±8.49a | 1.05±0.13a | 8.97±0.06b |
HS | 5.81±0.03b | 430.43±5.76a | 7.56±0.20b | 12.19±0.12a | 1.61±0.06a | 154.70±7.23a | 1.02±0.03a | 82.83±1.01a |
| TP g kg−1 | TK g kg−1 | AP mg kg−1 | AK mg kg−1 | Bacteria | Fungi | AOA | AOB |
log10 gene copies g−1 dry soil |
LS | 0.52±0.04b | 16.88±0.79b | 18.97±1.15b | 110.59±3.50b | 9.97±0.02a | 8.84±0.01b | 7.79±0.02b | 6.63±0.06b |
HS | 1.40±0.06a | 18.54±0.00a | 191.33±3.14a | 260.08±3.50a | 9.88±0.04a | 9.03±0.01a | 8.07±0.02a | 6.96±0.04a |
LS, low NO3− accumulation soil; HS, high NO3− accumulation soil. EC, electrical conductivity; SOC, soil organic carbon; TN, total nitrogen; C/N, SOC: TN ratio; DOC, dissolved organic carbon; NH4+, ammonium; NO3−, nitrate; TP, total phosphorus; TK, total kalium; AP, available phosphorus; AK, available kalium; AOA, ammonia-oxidizing archaea; AOB, ammonia-oxidizing bacteria. The different lowercase letters indicate the differences between LS and HS in the same property (p < 0.05). |
The 15N tracing study
Two vegetables (tomato and cucumber) and two greenhouse soils (HS and LS) were selected for 15N tracing pot experiments. The experiment was carried out with four treatments, i.e. (1) HST, tomato planted in HS; (2) LST, tomato planted in LS; (3) HSC, cucumber planted in HS; (4) LSC, cucumber planted in LS. Three sub-treatments were set up for each treatment, i.e. without plant (CK), planting for 15 days (D15) and 26 days (D26). The seeds of tomato and cucumber were surface-sterilized with 5% alcohol and then rinsed with deionized water. The sterilized seeds were germinated on petri dishes placing in an incubator with 30 oC. The germinated seeds were transplanted to the seeding plug tray containing nutritional material (peat: vermiculite: perlite, 6:3:1). After tomato and cucumber seedlings had fully developed their first true leaf, the nutritional materials of seeding roots were removed and seedlings with uniform growth rate were selected and transplanted to pots with 200 g of soil each (oven-dry basis). In addition, a control treatment (without plants) was set up for each soil. Hoagland’s nutrient solution was added to each pot. Topdressing (only N fertilizer) with a rate of 10 mg N kg−1 was applied every 6-days. All pots were cultured at daytime temperatures of 30-35 oC, night temperatures of 20-25 oC and under artificial light intensity 5000 Lux (12 h light and 12 h dark cycle). Furthermore, water was added each day to maintain the soil water content at 60-70% water holding capacity and the air was humidified with a humidifier to create an environmental condition of the greenhouse, e.g. high temperature and humidity.
15 days (D15) and 26 days (D26) after transplanting the seedlings, respectively, 15N tracing pot experiments were conducted. Briefly, two 15N labeled treatments (three repetitions for each treatment), i.e. 15N labeled NH4+ (15NH4NO3, 10.20 atom%) and 15N labeled NO3− (NH415NO3, 10.12 atom%) were set up. A 2 mL 15NH4NO3 or NH415NO3 solution was injected uniformly into each pot using a four-needle injection technique, to give a final concentration of 30 mg NH4+-N kg−1 soil and 30 mg NO3−-N kg−1 soil. At 0.5, 24, 48, 72 h after addition of 15N labeled solutions, soil and whole plant were collected. NH4+ and NO3− in soil were extracted using 1 M potassium chloride solution (KCl) at a solution:soil ratio of 5:1 to determine their concentrations and 15N abundances. All plants were soaked in 1 M KCl solution for five minutes to remove residual inorganic N, and thoroughly rinsed with deionized water. Plant samples were then de-enzymed at 105 oC for 0.5 h, dried at 80 oC to a constant weight, grinded, and passed through a 0.15 mm sieve to measure the N concentrations and 15N abundances. In addition, soil pH, dissolved organic C (DOC) concentration, the abundances of bacteria, fungi, ammonia-oxidizing archaea (AOA) and ammonia-oxidizing bacteria (AOB) were also determined.
Analytical methods
Soil pH was measured in a 2.5:1(v/w) water/soil ratio using a DMP-2mV/pH detector (Quark Ltd., Nanjing, China). Soil electrical conductivity (EC) was determined by a conductivity detector (Kang Yi Crops., Nanjing, China) at a 5:1 (v/w) water/soil ratio. Soil organic carbon (SOC) was analyzed by wet digestion with H2SO4-K2Cr2O7 and total N (TN) was determined by semi-micro Kjeldahl digestion using Se, CuSO4, and K2SO4 as catalysts. Soil dissolved organic C (DOC) was extracted with 1 M KCl solution in a 5:1 KCl solution:soil ratio, and then measured using an Analyzer Multi N/C (Analytic Jena, Jena, Germany). Total soil phosphorus (TP) was determined with NaOH fusion and Mo-Sb colorimetry, available phosphorus (AP) with NaHCO3 digestion and Mo-Sb colorimetry, total potassium (TK) with NaOH fusion and flame photometer, as well as available kalium (AK) with CH3COONH4 digestion and flame photometer (for details see Bao 2000). The soil was extracted with 1 M KCl solution, and concentrations of NH4+ and NO3− in the extracted solution were analyzed by an automated continuous flow wet chemistry analyzer (SA1000, Skalar, Netherlands). A micro-diffusion method with magnesium oxide and Devarda’s alloy combining with oxalic acid absorption was used to separate NH4+ and NO3− in the extracts. Subsequently, the 15N abundances of NH4+ and NO3− were determined with an Isotope-Ratio Mass Spectrometry (IRMS) system (Europa Scientific Integra, Crewe, UK). Plant 15N abundance was analyzed using the IRMS system with a Finnigan Flash EA 1112 elemental analyzer (Thermo Scientific, Bremen, Germany).
DNA was extracted from 0.25 g fresh soil using the Power Soil DNA Isolation Kit (Bioer Technology Co., Ltd) according to the manufacturer’s instructions. The qPCR (quantitative real-time polymerase chain reaction) amplifications were performed in eight-well tubes with a CFX96TM Real-Time System (Bio-Rad Laboratories Inc., Hercules, CA, USA). The abundances of bacteria (Eub338F/Eub518R), fungi (ITS1F/ ITS2R), AOA (Arch-amoAF/Arch-amoAR) and AOB (amoA1F/amoA2R) were measured using the primer sets and thermal conditions listed in Table S1 (Supporting information). The reaction mixture for bacteria, fungi, AOA and AOB was set according to Zhao et al. (2017). The standard curves were established according to Huang et al. (2015) and Zhao et al. (2017), and the amplification efficiencies for bacteria, fungi, AOA and AOB were 115%, 98%, 100%, 94% respectively. Amplification specificity was evaluated by melting-curve analysis and copy numbers were then converted to oven-dry basis.
The 15N tracing tool, Ntraceplant, containing six N pools and twelve processes (Fig. 1), was used to simultaneously quantify soil gross transformation and plant N uptake rates (He et al. 2020). The Markov Chain Monte Carlo (MCMC) algorithm, able to calculate true parameter uncertainties, was used for data analysis. The misfit function in the model between the modeled and observed data f(m) considers the variance of the individual observations (see eqt. 3 in (Müller et al. 2007)). The concentrations and 15N abundances of soil NH4+, NO3−, and plant N (mean ± standard deviation) from two 15N labeled treatments were supplied to Ntraceplant. The average gross N transformation rates over the experimental period were quantified by simultaneously optimizing the kinetic parameters, set to zero-order, first-order, or Michaelis-Menten kinetics, by minimizing the misfit between the modeled and observed values, and expressed in units of mg N kg−1 soil d−1. The most suitable model was selected according to Akaikes Information Criterion (AIC) (Cox et al. 2006) within the MATLAB computing environment (Version 7, The MathWorks Inc.) (for details see He et al. 2020).
Due to the faster growth rate of cucumber than tomato, the N uptake rates of cucumber were higher than that of tomato, and thus, NH4+ and NO3− in soil planted with cucumber could not be detected at 72 h after addition of 15N labeled solution. Therefore, the data applied to the model included concentrations and 15N abundances of soil NH4+, NO3−, and plant N (mean ± standard deviation) of three sets of data (i.e. 0.5, 24, 48 h) in soil planted cucumber, while four sets of data (0.5, 24, 48, 72 h) in soils planted with tomato and without plants.
Calculations and statistical analyses
Total gross N mineralization rates (M) = MNlab + MNrec
Total gross NH4+ immobilization rates (INH4) = INH4Nlab + INH4Nrec
Total gross nitrification rates (N) = ONH4 + ONrec
Total gross microbial N immobilization rates (ITN)= INH4 + INO3
Total plant N uptake rates (UTN) = UNH4 + UNO3
where MNlab = mineralization rate of labile organic N pool to NH4+, MNrec = mineralization rate of recalcitrant organic N pool to NH4+, INH4Nlab = immobilization rate of NH4+ to labile organic N pool, INH4Nrec = immobilization rate of NH4+ to recalcitrant organic N pool, INO3 = immobilization rate of NO3− to organic N pool, ONH4 = oxidation of NH4+ to NO3−, ONrec = oxidation of organic N pool to NO3−, UNH4 = plant NH4+ uptake rate, and UNO3 = plant NO3− uptake rate.
Significant differences of each N transformation rate between means of CK, D15 and D26 in each treatment were estimated by one-way ANOVA (analysis of variance). Based on the actual experimental repetitions (three repetitions), SigmaPlot (Version 3.5) and Excel 2019 were used to calculate the LSD (least significant differences) at 0.05 significance level for each N transformation rate, which was the most conservative way to calculate LSDs (Müller et al. 2011). The T-test analysis was used to examine the differences in soil properties between studied two soils. The relationships among N transformation rates were examined using Spearman correlation analyses with two-sided tests.