Wide-spread limitation of soil organic nitrogen transformations by substrate availability and not by extracellular enzyme content

Proteins constitute the single largest soil organic nitrogen (SON) reservoir and its decomposition drives terrestrial N availability. Protein cleavage by extracellular enzymes is the rate limiting step in the soil organic N cycle and can be controlled by extracellular enzyme production or protein availability/stabilization in soil. Both controls can be affected by geology and land use, as well as be vulnerable to changes in soil temperature and moisture/O2. To explore major controls of soil gross protein depolymerization we sampled six soils from two soil parent materials (calcareous and silicate), where each soil type included three land uses (cropland, pasture and forest). Soil samples were subjected to three temperature treatments (5, 15, 25 °C at 60% water-holding capacity (WHC) and aerobic conditions) or three soil moisture/O2 treatments (30 and 60% WHC at 21% O2, 90% WHC at 1% O2, at 20 °C) in short-term experiments. Samples were incubated for one day in the temperature experiment and for one week in the moisture/O2 experiment. Gross protein depolymerization rates were measured by a novel 15N isotope pool dilution approach. The low temperature sensitivity of gross protein depolymerization, the lack of relationship with protease activity and strong effects of soil texture and pH demonstrate that this process is constrained by organo-mineral associations and not by soil enzyme content. This also became apparent from the inverse effects in calcareous and silicate soils caused by water saturation/O2 limitation. We highlight that the specific soil mineralogy influenced the response of gross depolymerization rates to water saturation/O2 limitation, causing (I) increasing gross depolymerization rates due to release of adsorbed proteins by reductive dissolution of Fe- and Mn-oxyhydroxides in calcareous soils and (II) decreasing gross depolymerization rates due to mobilization of coagulating and toxic Al3+ compounds in silicate soils.


Introduction
Between 50 and 90% of soil N occurs in the form of high molecular weight (HMW) peptidic N, comprising (glyco)proteins and peptide side chains in peptidoglycan (Martens and This is an open access article under the CC BY license (https://creativecommons.org/licenses/BY/4.0/). Loeffelmann, 2003), the rest deriving from amino sugar polymers (e.g. chitin and glycan strands in peptidoglycan), aromatic N compounds (e.g. nucleic acids) and to date uncharacterized organic N compounds (Schulten and Schnitzer, 1997). These HMW-organic N compounds represent the largest soil N reservoir but are not directly available for microorganisms and plants due to molecular size constraints on cellular uptake mechanisms (Schimel and Bennett, 2004). Depolymerization of these complex organic N compounds is therefore considered the rate limiting step fueling the terrestrial N cycle and allowing organisms to utilize soil organic N compounds (Schimel and Bennett, 2004;Jan et al., 2009;Jones et al., 2009;Hu et al., 2017). The breakdown of the largest fraction of HMW-organic N, proteinaceous substances, to low molecular weight organic compounds such as oligopeptides and amino acids is catalyzed by extracellular proteases and peptidases, secreted into the soil by micro-organisms and plants (Adamczyk et al., 2008;Fuka et al., 2008;Vranova et al., 2013). The products of this breakdown represent a highly dynamic and easily available pool of organic N in soils with half-life times of minutes to hours (Jones et al., 2009;Wanek et al., 2010). In decomposing litter and in soils, protein depolymerization rates by far exceed those of inorganic N processes such as organic N mineralization and nitrification i.e. by > 8 fold Prommer et al., 2014;Wild et al., 2015), showing that organic N cycling exceeds inorganic N cycling in the terrestrial N cycle. Moreover, isotope tracer studies showed that uptake of intact amino acids and small peptides contributes significantly to the N nutrition of soil microorganisms and plants (Gallet-Budynek et al., 2009;Hill et al., 2011;Farrell et al., 2012), short-cutting the soil inorganic N cycle. However, despite its importance and first insights into the process of gross protein depolymerization in soils, we still lack a thorough mechanistic understanding of major controls on this process.
In principle enzyme-catalyzed depolymerization processes in soils can be controlled either by extracellular enzyme content (enzyme limitation) or by resource availability (substrate limitation, Wallenstein and Weintraub (2008)). Both are strongly affected by environmental conditions, such as by soil temperature and moisture controlling substrate diffusion and catalytic rates, and by soil physicochemical factors that influence protein availability and microbial community structure (Fig. 1).
Soil extracellular enzyme content, measured as potential enzyme activities at substrate saturation and optimal pH, depends on (i) microbial community structure and activity and (ii) the microbial N and C demand, both affecting the production of extracellular enzymes, and (iii) sorption of enzymes/protein substrates in soils changing their activity and decreasing proteolytic decomposition (Quiquampoix, 2000). All these parameters are affected by land use and parent material (Fig. 1). The set of excreted proteases is shaped by the microbial community structure (Fuka et al., 2008). Soil peptidases are produced by archaea, bacteria, and fungi, belonging to ∼150 peptidase subfamilies, though the occurrence and abundance of the genes coding for secreted peptidases varies across the three microbial kingdoms (Nguyen et al., 2017). Land use is known to be an important factor influencing soil microbial community structure and bacterial to fungal biomass ratios (Fierer and Jackson, 2006;Jangid et al., 2008;Lauber et al., 2009). Moreover, land-use related N inputs by synthetic fertilizers or manure may further suppress (inorganic N) or stimulate (organic N) soil protease production (Geisseler and Horwath, 2008). The production and potential activity of extracellular proteases is therefore driven by substrate supply (Brzostek and Finzi, 2011), litter quality (Wickings et al., 2012) and consequently microbial C and N demands (Geisseler and Horwath, 2008). Moreover, soil parent material shapes soil physicochemical properties such as soil pH and texture which in turn affect microbial community structure and the excreted enzyme complement (Fierer and Jackson, 2006;Kallenbach et al., 2016;Rousk et al., 2010), but also determine the number of sorption sites for enzymes and protein substrates.
Soil protein availability depends (i) on the input of labile proteins by plant litter, microbial turnover and manure, and (ii) on stabilization mechanisms (sorption/desorption equilibria) and protein accessibility (aggregate formation), both of which are strongly affected by soil parent material and land-use. Land-use and the respective plant communities affect litter quality and quantity (Fanin and Bertrand, 2016) and thereby the labile protein input into the soil ecosystem. Microbial C limitation due to reduced litter input in cropland soils promotes organic N mineralization as microbes decompose available organic N compounds to meet their energy demands and excrete excess N as NH 4 + (Mooshammer et al., 2014). This is further enhanced by tillage or other disturbances in cropland soils which promote microbial decomposition of soil organic matter (Six and Jastrow, 2002). Microbial turnover in soils also adds proteins from lysing cells (Schimel and Bennett, 2004) and turnover varies with vegetation type, soil depth and fertilizer regime (Spohn et al., 2016;Spohn and Widdig, 2017). Though proteinaceous N makes up the bulk of N in soils and litter, in soils protein availability is strongly constrained by organo-mineral associations (Kögel-Knabner et al., 2008), rendering proteins inaccessible to proteolytic attack (Quiquampoix, 2000).
Fe-and Al-oxides and phyllosilicates are assumed to be the main sorption sites of soil organic matter and to form stable organo-mineral complexes by ligand exchange or electrostatic interactions (Kaiser and Guggenberger, 2003). Both, the net charge of sorbent (mineral phase) and adsorbent (protein) are determined by soil pH, due to protonation of functional groups with decreasing pH. Adsorption of proteins on mineral surfaces is highest at soil pH close to the isoelectric point, where the net charge of a protein is zero (Quiquampoix and Ratcliffe, 1992;Gu et al., 1994). Protein sorption capacity/strength of the soil matrix therefore depends on soil texture, soil mineral composition and soil pH and thereby is strongly affected by soil parent material. Moreover, Mn-oxides, like birnessite catalyze the abiotic hydrolysis of peptide bonds and thus could either increase depolymerization rates by abiotic cleavage of proteins or decrease them by fragmentation and deactivation of proteolytic enzymes (Reardon et al., 2016). We therefore expect strong effects of parent material and land-use on the protein depolymerization process, but rigorous tests thereof are lacking.
Moreover, so far no studies have tested the effects of soil temperature and soil moisture on gross protein depolymerization rates. With increasing soil temperature, catalytic rates of extracellular peptidases increase exponentially with Q 10 values of ∼2, while temperature effects are weaker on diffusion, sorption and desorption processes of soil organic matter and protein substrates (Q 10 ∼1, Conant et al., 2011). Recent studies have shown that extracellular peptidases are less temperature sensitive than C or P cycling enzymes, resulting in a decoupling of protein depolymerization from C and P acquisition at increasing temperatures (Brzostek and Finzi, 2012;Steinweg et al., 2018). Depending on enzyme-or substratelimitation of gross protein depolymerization, we therefore would expect rates to increase with Q 10 reflecting extracellular enzyme responses or to remain unaltered, respectively (Steinweg et al., 2012).
Soil moisture is also a strong driver of in situ decomposition processes. Diffusion of catalytic enzymes and substrates throughout soil is determined by the soil water content (Moldrup et al., 2001) and modified by soil texture. Soil water content also directly affects soil microbial life. At low soil water content, soil microorganisms are exposed to osmotic stress and consequently reduce their internal water potential by accumulation of osmolytes like amino compounds (Csonka, 1989;Witteveen and Visser, 1995;Schimel et al., 2007). In contrast, at high soil water contents, close to field capacity, oxygen becomes a limiting factor causing a shift from aerobic to anaerobic respiration and fermentation. Since, microbial activity is constrained by the lower energetic efficiency of anaerobic respiration and fermentation compared to aerobic respiration (Ponnamperuma, 1972;Picek et al., 2000), this leads to reduced microbial SOM decomposition in water-logged soils in the long term (Kögel-Knabner et al., 2010;Schädel et al., 2016). However, in the short term reductive dissolution of Fe-and Mn-oxyhydroxides causes a release of stabilized SOM and thereby increases decomposition (Kögel-Knabner et al., 2010).
Here we present the first comprehensive study on the above cited potential controls of the protein depolymerization process, based on two incubation experiments, assessing temperature and soil moisture/O 2 effects on gross protein depolymerization and amino acid uptake rates in soils differing in soil parent material (calcareous versus silicate) and land use (forests, pastures and croplands). To assess short-term temperature and moisture effects on gross protein depolymerization (and amino acid uptake) rates, we conducted short-term incubation experiments at three levels of temperature (5, 15 and 25 °C at 60% WHC) and moisture (30 and 60% WHC at 21% O 2 , and 90% WHC at 21% O 2 , all at 20 °C). We hypothesized, that (i) gross protein depolymerization is not enzyme but mainly substrate limited and that protein availability and protein sorption is strongly affected by soil parent material and land use. Second (ii) if protein depolymerization were enzyme limited, higher temperatures will induce higher gross rates due to increased enzyme activity while if substrate-limitation prevails, the constraint by substrate availability will suppress the temperature sensitivity of protein depolymerization. Finally, we hypothesized that (iii) an increase in soil moisture will cause increased protein depolymerization rates due to enhanced diffusion of substrates and enzymes throughout the soil. In contrast, at high soil water contents in association with low oxygen concentrations, we expected that protein depolymerization would be less substrate limited due to dissolution of organo-mineral complexes and microbial amino acid uptake will decrease due to reduced microbial activity under suboxic conditions.

Site description and soil sampling
Soil samples were collected in June 2016 in the Enns valley in the Eastern Alps (Austria) at 700m a.s.l. Mean annual air temperature and precipitation are respectively 7 °C and 1014 mm (Starz et al., 2011). Two different soil parent materials were selected: (1)  At each of the six sites four replicate soil samples were taken with a root corer (diameter 7 cm) to 15 cm soil depth. Prior to soil sampling vegetation and organic horizons were removed. The four replicates were sampled to cover heterogeneities in site topography (hilltop, upper and lower slope, bottom of slope) and were kept and analyzed separately. All soils were first sieved to 4 mm, roots and stones were removed by hand, then homogenized by sieving to 2 mm and stored at 4 °C until further analyses.

Basic soil parameters and microbial biomass
Soil texture, carbonate content, base saturation (BS), effective cation-exchange capacity (CEC eff ) , and exchangeable Na + , K + , Mg 2+ , Ca 2+ , Fe 3+ , Al 3+ and H + were determined by the soil analysis laboratory of the Federal Office for Food Safety (AGES) in Vienna, Austria in air dried soil samples. All analyses were carried out according to Austrian and European standards (ÖNORM). To determine soil water content, sieved soils were dried at 60 °C for 48 h. Soil water-holding capacity (WHC) was determined by repeatedly saturating 10 g field-moist soil with deionized water and draining in between for 2.5 h in a funnel with an ash-free cellulose filter paper. Soil pH was measured in ultra-pure water and in 10 mM CaCl 2 solution (soil:solution ratio = 1:5 (w:v)) using an ISFET pH sensor (Sentron, Leek, The Netherlands). For total soil organic C (SOC) and total soil N (TN), carbonate was removed in soil aliquots with 2 M HCl, subsequently samples were oven-dried, ground in a ball mill (MM 200,Retsch,Germany) and analyzed by an Elemental Analyzer (Carlo Erba 1110, CE Instruments) coupled to a Delta Plus Isotope Ratio Mass Spectrometer (Finnigan MAT, Germany) via a Conflo III interface (Thermo Fisher, Austria). For dissolved C and N pools, 4 g aliquots of field-moist soil were extracted with 1 M KCl (1:5 w:v) for 1 h, filtered through ash-free cellulose filters and stored at −20 °C until analysis. Dissolved organic C (DOC) and total dissolved N (TDN) were measured by a TOC/TN analyzer (TOC-VCPH/ TNM-1, Shimadzu, Austria). Ammonium and nitrate concentrations were measured in the same extracts by colorimetric assays as described by Hood-Nowotny et al. (2010). Dissolved organic N (DON) was calculated as TDN minus ammonium and nitrate. Total primary amines (hereafter, free amino acid concentration (FAA)) were measured fluorimetrically using the OPAME assay (Jones et al.,2002;Darrouzet-Nardi et al., 2013) as modified by Prommer et al. (2014). Interference of NH 4 + was corrected using the NH 4 + concentrations determined via colorimetric assays, and N compounds other than amino acids (mainly free amino sugars) contributed only between 3 and 15% to primary amine-N in the same soils (Hu et al., 2018). Microbial biomass C and N were estimated using the chloroformfumigation extraction method (Brookes et al., 1985). Soils were fumigated with chloroform for 48 h and extracted (1:5 (w:v)) with 1 M KCl. Soil total P (TP) and soil total inorganic P (TIP) were determined in 0.5 M H 2 SO 4 extracts of ignited (450 °C, 4 °C, Kuo (1996)) and untreated soil aliquots by malachite green measurements of reactive phosphate (Lajtha et al., 1999). Total soil organic P (TOP) was calculated as the difference of TP -TIP. Microbial P was estimated using the chloroform-fumigation extraction method (Brookes et al., 1985). Soils were fumigated with chloroform for 48 h, extracted (1:5 (w:v)) with 0.5 M NaHCO 3 for 1 h and filtered through ash-free cellulose filters. Dissolved inorganic P (DIP, Olsen P) was measured in 0.5 M NaHCO 3 extracts (1:7.5 (w:v), pH 8.5) by the malachite green method. Total dissolved P was measured following acid persulfate digestion (Lajtha et al., 1999) and dissolved organic P (DOP) was calculated as the difference of P concentration in non-digested and digested samples. Soil microbial community composition was analyzed by phospholipid fatty acids (PLFAs) according to the procedure by Kaiser et al. (2010) and Hu et al. (2018).

NaOH extractable protein
Proteins in soils are stabilized by sorption to humic substances, clay minerals and Fe-/Alhydroxides (Schulten and Schnitzer, 1997;Lehmann and Kleber, 2015). NaOH extracts free and loosely bound proteins but not proteins stabilized in metal-organic complexes (Wattel-Koekkoek et al., 2001). We assumed that the more labile fraction extractable by NaOH therefore represents the protein fraction available for enzymatic cleavage. Aliquots of 2 g field-moist soil were extracted with 0.5 M NaOH for 2 h in an ultrasonic bath (160 W, Sonorex RK510, Germany) and subsequently for 16 h on a rotary shaker (1:10 (w:v)). Extracts were centrifuged for 15 min at 10.000×g at room temperature and total free amino acid (FAA) and ammonium concentrations measured fluorimetrically and colorimetrically in the supernatant as mentioned above. In a second aliquot of the supernatant, proteins were hydrolyzed with 6 M HCl at 90 °C overnight. The hydrolyzed extracts were dried under a gentle N 2 stream, dissolved in ultra-pure water, and dried again in a vacuum concentrator (SC 110 SpeedVac, Savant, USA) to remove residual HCl. After re-dissolving in ultra-pure water the amino acid and ammonium concentrations were re-measured as above. The difference between the FAA concentrations in the hydrolyzed extracts (total bound amino acids) and in the non-hydrolyzed extracts (total free amino acids) was used as a measure for NaOH-extractable protein.

Incubation experiments
To assess short-term temperature and moisture/O 2 sensitivities of organic N cycling processes and enzyme activities, we conducted two incubation experiments. One week prior to starting the temperature experiment and isotope pool dilution measurements, soils were adjusted to 60% WHC and stored at 15 °C in polyethylene bags. Soils were aerated every second day and water content was adjusted when necessary. For the temperature sensitivity measurements aliquots of 4 g pre-incubated soil were transferred to 50 mL polypropylene centrifuge tubes and soils were allowed to adjust to the respective temperature (5, 15 and 25 °C) overnight. The moisture/O 2 sensitivities of organic N cycling processes were studied in a second incubation experiment, run four weeks after the first experiment. Aliquots of soil (stored at 15 °C and 60% WHC before) were first adjusted to the respective soil moisture contents (30, 60 and 90% WHC), weighed into polypropylene tubes and pre-incubated for 1 week at the respective WHC and O 2 contents at 20 °C. Because water logging in nature comes with suboxic conditions causing microbial oxygen limitation, the 90% WHC treatment aliquots were incubated in a suboxic chamber at 1% O 2 . All other moisture treatments were run under normal air (21% O 2 ). To avoid changes in O 2 concentration during the pre-incubation phase, incubation vials were opened every second day for a few minutes. The isotope pool dilution experiments were started immediately after preincubation.

15 N isotope pool dilution (IPD) experiments
Gross protein depolymerization and amino acid uptake rates were assessed by an isotope pool dilution (IPD) assay as described by Wanek et al. (2010) and Noll et al. (2019) using 15 N labeled amino acids (U-15N-98 at% 15 N amino acid mixture from crude algal protein, Cambridge Isotope Laboratories, Radeberg, Germany). Pool sizes of FAA were determined immediately before starting the IPD experiments. The added 15 N-label amounts were adjusted to account for ∼20% of the FAA pool size, avoiding a strong fertilization effect by tracer addition. Preliminary experiments showed that the added 15 N-FAAs are rapidly immobilized in soils. Therefore, we used highly enriched 15 N labeled amino acids (98 at %) as tracers added to the soils. The tracer was dissolved and diluted with ultra-pure water to the respective concentrations and 400 μL of label solution were added drop wise to 4 g soil with a pipette and vigorous shaking to guarantee good mixing. Soil samples were incubated at the respective temperatures (5, 15 and 25 or 20 °C) for 15 and 45 min. Preliminary experiments showed that the recovery of the amino acid label decreases at constant rates between 15 and 60 min after tracer addition (Noll et al., 2019). The incubation was stopped by adding 20 mL of 1 M KCl solution at room temperature. Hu et al. (2017) showed that enzyme activity is effectively quenched with KCl at 4° and 15 °C in the absence of formaldehyde. Soils were extracted for 1 h on a rotary shaker and then filtered through ashfree cellulose filters. The charged nature of amino acids renders them to bind to positively (anion-exchange mechanisms) or negatively charged sites in soils (cation-exchange mechanisms) and therefore to be unavailable for extraction by water. 1 M KCl has been shown to achieve higher extraction efficiency than 0.01 M CaSO 4 or water (Noll et al., 2019). Moreover, Dippold and Kuzyakov (2013) showed in experiments with life soils and soils inhibited with HgCl 2 and NaN 3 , that only 3-6% of added alanine was sorbed to the soil matrix after 36 h and 83-90% of the sorbed alanine could be extracted with 0.5 M CaCl 2 .
Extracts were stored at −20 °C until further analyses. Isotope analyses were performed based on a novel isotope ratio mass spectrometry (IRMS) protocol adopted from Zhang and Altabet (2008)  Therefore aliquots (2 mL) of soil extracts were transferred to 12 mL exetainer glass vials with PTFE septa. Thereafter, 50 μL 60 mM KBr and 20 μL 10 M NaOH were added and exetainers were closed tightly. Finally, 50 μL NaClO (30 mM) were added with a syringe and samples were incubated for 30 min at 50 °C. To deactivate excess NaClO, 100 μL 0.4 M Na-metaarsenite was added after incubation. Subsequently NO 2 − was reduced by NaN 3 to N 2 O by adding 140 μL NaN 3 (20 mM; Sigma) in 50% acetic acid (Merck) with a gas tight syringe and samples were incubated for 30 min on an orbital shaker at room temperature (Lachouani et al., 2010). Finally, NaN 3 was deactivated by addition of 80 μL 6 M NaOH. Free amino acid concentrations and 15 N enrichments in the produced N 2 O were measured by PT-IRMS, and concentrations (equimolar 20 FAA mix) and isotope enrichments (natural abundance to 20 at% 15 N) calibrated by respective standards (Lachouani et al., 2010). PT-IRMS analyses were performed on a purge-and-trap isotope ratio mass spectrometer (Finnigan Delta V Advantage IRMS, Thermo Fisher, Germany) consisting of a Gasbench II headspace analyzer (Thermo Fisher, Germany) with a cryo-focusing unit.

Potential extracellular peptidase activity
Leucine-amino peptidase activities (EC 3.4.11.1) were measured in 2 mL soil slurries using a fluorometric microtiter plate assay with 50 μL L-leucine-7-amido-4-methyl coumarin (AMC-leucine, 1 mM) as a substrate in Na-acetate buffer (100 mM, pH 5.5, Kaiser et al., 2010). The samples were run in triplicates and were incubated for 2 h at the respective temperatures (5, 15 and 25 °C for the temperature experiment, and 20 °C for the moisture experiment) and measured every 30 min. Fluorescence was measured with a TECAN InfiniteR M200 (Austria) at an excitation wavelength of 365 nm and an emission wavelength of 450 nm, and was corrected for sample blank and quenching prior to calculations of AMC concentration.

Data analysis
Gross protein depolymerization rates (GP) and gross microbial amino acid uptake rates (GI) were calculated following the equations by Kirkham and Bartholomew (1954) and : (2) where N t1 and N t2 are the concentrations of FAA-N at the time points t 1 (15 min) and t 2 (60 min). 15 N enrichments in amino acids at the time points of termination are expressed as at % 15 N t1 and at% 15 N t2 , and at % 15 N b is the background 15 N abundance (0.366 at% 15 N) in non-labeled samples.
Temperature sensitivities of gross protein depolymerization rates, amino acid uptake rates and peptidase activity between 5 and 25 °C were calculated by a log linear Q 10 function (Janssens and Pilegaard, 2003): where R is the measured rate or enzyme activity, Q 10 is the temperature sensitivity, T is the incubation temperature and b a fitted coefficient. Q 10 and b where determined by non-linear least square estimates.

Statistics
Statistics and estimations of Q 10 values were performed with R version 3.1.3 R core (R Development Core Team, 2008). Data were log transformed or rank normalized if not normally distributed or hetero-scedastic. Correlations were done using Pearson's product moment correlation. In addition, we performed generalized linear mixed models (GLMM) of the influence of soil parent material on the relationship of selected soil parameters and gross protein depolymerization rates using the "lmer" function included in the "lme4" package (Bates et al., 2014). Soil parent material was specified as random factor. Regressions between gross protein depolymerization rates and NaOH-extractable protein were done using nonlinear least-squares estimates of the parameters in a nonlinear model. Differences in gross rates and peptidase activities between soil parent material, land use or the respective temperature and moisture/O 2 treatments were addressed by three-way analysis of variance (ANOVA). If effects were significant (p < 0.05) differences between groups were tested with Tukey's HSD. Significant differences in Q 10 values between parent material and land use were addressed by two-way ANOVA followed by Tukey's HSD. We used soil parameters which correlated significantly with gross protein depolymerization rates or free amino acid contents to develop a base multiple regression model. The most parsimonious model was determined by step wise backward/forward multiple regression analysis using Akaike's Information Criterion (AIC). R-squared values were obtained by regressing model predictions against observed values. Community composition based on PLFA analyses was visualized by nonmetric multidimensional scaling (NMDS) ordinations of individual fatty acid methyl esters. Protein depolymerization and soil parameters were fitted onto the ordination using the 'envfit' function in the R package 'vegan' (Oksanen et al., 2016). Significance of correlations was tested with 999 permutations.

Basic soil parameters and organic N pool sizes
The studied soils showed significant differences in physicochemical soil parameters and N pool sizes depending on parent material and land use (Table 1). Clay and silt contents were higher in calcareous than in silicate soils, while silicate soils had significantly higher sand contents than calcareous soils (p < 0.001, ANOVA). Soil pH, cation-exchange capacity (CEC eff ), base saturation (BS) and exchangeable Ca 2+ were higher in calcareous soils than in silicate soils (p < 0.001, ANOVA). Amorphous Al (Al o ) contents were significantly higher in silicate soils than in calcareous soils (p < 0.001, ANOVA), while calcareous soils had significantly higher crystalline Fe (Fe d-o ) contents (p < 0.01, ANOVA). We also found a strong land use effect (p < 0.001, ANOVA) with increasing soil pH in the order forest < pasture < cropland. NaOH-extractable protein-N was significantly higher in calcareous soils than in silicate soils (F = 16.2, p < 0.001), and ranged between 11 and 66% of total soil N ( Table 1). We also observed a significant effect of land use (F = 15.9, p < 0.001) represented by increasing NaOH-extractable protein-N contents in the order cropland < forest < pasture. Total free amino acids (FAA) were higher in silicate soils than in calcareous soils (p < 0.001) and accounted for 27 ± 10% of the dissolved organic N pool in silicate soils and for 13 ± 5% in calcareous soils, but were unaffected by land use. Significant effects of soil parent material and land use were found for total, extractable and microbial C and N pools (Table S1). PLFA contents were highest in calcareous soils, but did not vary with land use, and there was no significant parent material or land use effect on fungi:bacteria ratios (Table  S1).

Effects of parent material and land use on soil processes
Gross protein depolymerization rates ranged between 24 and 280 μg N g d. w. −1 d −1 and showed significant parent material (F = 18.21, p < 0.001) and land use effects (F = 10.45, p < 0.001, Table S2). In detail, gross protein depolymerization rates were significantly higher in silicate soils than in calcareous soils, and in pasture soils compared to cropland and forest soils. Highly similar patterns were observed for gross rates of microbial amino acid uptake.
To examine the effects of soil physicochemical parameters on gross rates of protein depolymerization and microbial amino acid uptake, we performed linear regression analyses (Table 2). Gross protein depolymerization and amino acid uptake rates were significantly negatively correlated with soil properties such as clay and silt content, exchangeable Ca 2+ and CEC eff and strongly positively correlated with sand content, DON and FAA-N, but not with microbial biomass C and N or PLFA contents (Table 2). Considering the strong parent material effects, we used generalized linear mixed models of gross protein depolymerization with soil parent material as random factor. We obtained significant results for the mentioned soil parameters, including the significant positive effect of NaOH-extractable protein and negative ones of CEC eff and clay content (Table S3). The relation of NaOH-extractable protein-N and gross protein depolymerization was fitted to a quadratic polynomial model using a generalized non-linear least squares fit, which explained 70% of the variation in silicate soils and 60% of the variation in calcareous soils (Fig. 2). However, across and within all soil types no significant correlations were found between potential peptidase activity and gross protein depolymerization rates. Following linear regression analyses, multiple regression models were tested for gross protein depolymerization, and the final most parsimonious model included soil N tot (positive effect), NaOH-extractable protein (positive effect), CEC eff (negative effect) and crystalline Fe (negative effect, Table S4). The model explained about 86% of the variance in gross protein depolymerization (r 2 = 0.86, F = 28.11, p < 0.001). Gross microbial amino acid uptake rates across different parent materials and land use types (Table 2) showed similar correlations as described above for gross protein depolymerization rates. Significant effects on amino acid up-take rates were observed for parent material (F = 19.3, p < 0.001) and land use (F = 13.5, p < 0.001) and their interaction (F = 8.6, p = 0.001). Furthermore, amino acid uptake rates were strongly correlated with protein depolymerization rates across all sites (Table 2). Potential leucine-amino peptidase activity was positively correlated to PLFA contents, microbial C and N, NaOH-extractable protein and soil pH (Table 2) while negative correlations were found with DOC, soil C:N ratio, Fe o , Al o , Al d and exchangeable Al 3+ . Soil NaOH-extractable protein was positively correlated to peptidase activity, microbial C and N as well as gram negative PLFA markers and fungal PLFAs while clay content and Al o negatively affected NaOH-extractable protein ( Table 2). Microbial community composition obtained from PLFA analyses correlated with soil pH and leucine-amino peptidase activity, CEC eff , DOP, total N and organic C but not with protein depolymerization rates (NMDS, Fig. S1).

Soil temperature effects
Incubation temperature (5, 15 25 °C) had a weak positive effect on gross protein depolymerization rates compared to strong effects of parent material and land use according to three-way ANOVA (Table 3, Fig. 3.). We also found a weak interaction effect of land use x temperature as cropland soils showed an increase of gross rates with temperature, while forest soils showed a minimum at 15 °C (Fig. 3). Interaction effects of soil parent material x temperature were not significant. Tukey's HSD did not reveal significant differences between single treatments across soils. Within and across all temperature treatments, peptidase activity was not related to depolymerization rates. However, we found a strong positive temperature effect on peptidase activity as well as weaker parent material and land use effects (Table 3, Fig. 3). Peptidase activities increased in the order forest < arable = pasture and were higher in calcareous soils than in silicate soils (Fig. 3). Gross amino acid uptake rates showed a similar pattern as depolymerization rates across the temperature treatments (Fig. 3), with a weak but significant positive temperature effect, and strong effects of parent material and land use (Table 3).
Results from a two-way ANOVA showed significant land use effects on Q 10 values of gross protein depolymerization (F = 6.5, p = 0.007) but no effect of soil parent material. Q 10 values of gross protein depolymerization (and amino acid uptake) showed an increase in the order forest (0.8-1.1) < pasture (1.2-1.4) < crop land (1.6-2.2; Table 4). Q 10 values of amino acid uptake rates (F = 7.6, p = 0.004) and peptidase activity (F = 9.7, p < 0.005) were also affected by land use and significant effects of parent material were observed for the Q 10 values of peptidase activity (F = 7.9, p < 0.05). Mean Q 10 values of peptidase activity ranged between 1.3 and 1.7 and were slightly higher in calcareous soils than in silicate soils (p = 0.01, Table 4). Q 10 values of gross protein depolymerization (1.6-2.2) therefore matched (or even exceeded) those of amino peptidase activity in cropland soils, were slightly lower in pasture soils (1.2-1.4 for protein depolymerization compared to 1.6-1.7 for peptidase activity) and no temperature sensitivity was apparent for forest soils (0.8-1.1) though Q 10 values of peptidase activity in forests ranged between 1.3 and 1.4 (Table 4).

Soil moisture effects
Prior to IPD experiments soils were pre-incubated at the respective soil moisture (30, 60, 90% WHC) and O 2 concentration (21% O 2 at 30 and 60% WHC, 1% O 2 at 90% WHC) for 7 days to allow equilibration of biochemical processes. After pre-incubation free amino acid concentrations (measured 1 day before starting the IPD experiment) were significantly higher in suboxic treatments (90% WHC, 1% O 2 ) than in oxic treatments (30 and 60% WHC, 21% O 2 , Fig. 4d), whereas nitrate was almost entirely depleted in all suboxic treatments (Fig. S2). Three-way ANOVA of gross protein depolymerization with the factors parent material, land use and soil moisture/O 2 showed significant main effects of parent material, land use and soil moisture/O 2 , though interaction effects of parent material x soil moisture/O 2 had the highest explanatory power (Table 5). More specifically, gross depolymerization rates did not differ significantly between 30 and 60% WHC across all parent material and land use types (Fig. 4a). In oxic treatments (30 and 60% WHC) we found similar patterns as described for the temperature treatments including lower gross depolymerization rates in calcareous soils and a parent material dependent land use effect (Fig. 4a). This pattern changed for the suboxic treatment at 90% WHC. In silicate soils at 90% WHC/1% O 2 gross protein depolymerization rates decreased significantly across all land use types (p < 0.05, Tukey HSD), and the average decreases accounted for 54-74% of the rates at 60% WHC. In calcareous soils at 90% WHC/1% O 2 gross protein depolymerization rates increased more than two-fold in cropland soils (p < 0.05, Tukey HSD) compared to the 60% WHC treatment and in pasture soils (p < 0.01, Tukey HSD) compared to the 30% WHC treatment, while no significant changes were observed for the calcareous forest soil.
We found similar effects of parent material and land use on gross amino acid uptake rates as we reported for gross depolymerization rates but no significant moisture effect/O 2 (Table 5). On average gross amino acid uptake rates in the 90% WHC/1%O 2 WHC treatment decreased in silicate soils by 46-68% and increased by 2-to 3-fold in calcareous cropland (p < 0.05, Tukey HSD) and pasture soils (p < 0.01, Tukey HSD, Fig. 4b, Table 5) compared to the 60% WHC treatments. In contrast to gross depolymerization and amino acid uptake rates, parent material had the strongest explanatory effect on peptidase activity, followed by land use and soil moisture/O 2 (Table 5). Differences in peptidase activities between soil moisture/O 2 treatments were significant and showed a slight increase between 60 and 90% WHC (p < 0.001, Tukey HSD). Again peptidase activities were not related to gross protein depolymerization rates within or across soil moisture treatments.

Parent material and land use effects on soil organic N cycling
Our results demonstrate that gross protein depolymerization as catalyzed by extracellular enzymes varied considerably with soil parent material and land use. Soil parent material explained 29% of the variability in protein depolymerization and 23% of the variability in microbial amino acid uptake (Table S2). In detail, gross depolymerization rates responded negatively to clay and silt content, exchangeable Ca 2+ and CEC eff , resulting in lower gross rates in calcareous than in silicate soils (Fig. 3a). The dampened response of gross protein depolymerization rates to increasing NaOH extractable protein contents (representing a large proportion of the total soil protein pool) in calcareous compared to silicate soils (Table 1) suggests that proteins are less available for enzymatic cleavage, likely due to stronger sorption to mineral surfaces in calcareous soils (Rillig et al., 2007). This is supported by the negative effects of CEC eff and crystalline Fe on gross protein depolymerization rates (Table  S4). The studied calcareous soils had higher clay and silt contents than the silicate soils (Table 1 and Table S1), indicating a higher soil surface area and a higher proportion of more strongly sorbed or occluded protein. Moreover, in calcareous soils exchangeable Ca 2+ contents were higher, allowing enhanced Ca 2+ ion bridging between negatively charged proteins and negatively charged mineral surfaces to promote protein sorption (Lützow et al., 2006). These findings and the missing relation between gross depolymerization rates and potential peptidase activity support the hypothesis that gross rates of organic N cycling are substrate limited rather than enzyme limited, triggered by varying mechanisms of protein sorption and occlusion. In line with substrate limitation, Prommer et al. (2014) showed that protein depolymerization decreased after addition of biochar likely through protein sorption in micropores rendering them non-accessible to microbiota and extracellular enzymes. Soil mineralogy affects organic N cycling by stabilization of proteins, but also by abiotic fragmentation of proteins on Mn-oxides. Recent studies showed that proteins can be abiotically degraded in the presence of birnessite, a pedogenic Mn-oxide produced under alternating oxic and suboxic conditions (Russo et al., 2009;Reardon et al., 2016). Soils developed on limestone often contain considerable amounts of Mn-oxides (Dixon and White, 2002), however exchangeable Mn concentrations were not significantly different between the studied calcareous (0.07 ± 0.07 cmol c kg −1 ) and silicate soils (0.06 ± 0.04 cmol c kg −1 ). Though intriguing, it is unclear if abiotic degradation of proteins might enhance organic N availability due to liberation of small peptides and amino acids, or even decreases organic N availability due to deactivation of extracellular proteolytic enzymes (Reardon et al., 2016).
Moreover, the studied soils covered three land use types (forest, pasture, cropland) which differ in organic matter quality, quantity and type of N input, cycling of N and other essential nutrients such as P and K, and in microbial community structure (Tables 1 and S1). Land use explained about 30% of the variability of gross protein depolymerization and microbial amino acid uptake (Table S2), with highest NaOH-extractable protein contents and gross protein depolymerization rates in pasture soils, likely induced by high organic N inputs by cattle or sheep in the form of slurry or feces. Protein availability is also influenced by plantmicrobe-mineral interactions. Root exudates, like organic acids, weaken organo-mineral interactions and thereby release mineral-associated proteins. Moreover, excretion of labile root exudates might enhance microbial N mining due to stimulated microbial activity (Jilling et al., 2018). Since pasture soils are often characterized by high root densities, root exudates might hamper protein stabilization in pasture soils compared to cropland soils. Low protein depolymerization and amino acid uptake rates in cropland soils corresponded to lowest NaOH-extractable protein contents; however in cropland soils Q 10 values of gross depolymerization were similar to those of peptidase activity, indicating enzyme limitation rather than substrate limitation (Table 4). The lower depolymerization rates in forest soils compared to pasture soils might be caused by complexation of proteins with organic substances such as tannins rendering them less accessible for soil proteases (Joanisse et al., 2007). At low soil pH < 4.5, as in the silicate forest soils, formation of organo-metalcomplexes also contributes to the stabilization of SOM and of proteins (Nierop et al., 2002). This might also explain the negative effect of organic C on gross depolymerization rates observed in our model (Table S4). Although microbial community composition was driven by soil pH and had a significant effect on peptidase activity, we found no effect on protein depolymerization rates, supporting the crucial role of soil physicochemical properties and substrate availability on organic N transformation rates in situ.

Temperature effects on soil organic N cycling
Microbial and enzyme activities are generally assumed to be accelerated by higher temperatures, causing faster decomposition of soil organic matter (Kirschbaum, 1995). In our study the short-term experiment temperature changes explained only 7% of the variability of gross protein depolymerization rates and 11% of the variability of microbial amino acid uptake rates (Table 3). In contrast, potential peptidase activity increased significantly with temperature for all soils, with an overall Q 10 value of 1.53 (range 1.3-1.7, Table 4) and temperature explained 72% of the variability (Table 3). Similar Q 10 values of potential peptidase activity around 1.6 have also been reported for alpine organic and mineral forest soils (Koch et al., 2007;Schindlbacher et al., 2015). Low temperature sensitivity of N cycling processes are in accordance with results from soil warming experiments, reporting also no significant response of organic N mineralization rates to instantaneous and long-term increased soil temperatures (Schindlbacher et al., 2015).
Moreover, the temperature sensitivity of protein depolymerization of fungal and bacterial cell wall material (e.g. chitin and peptidoglycan) measured on the same soils was lower than that of exochitinase (Hu et al., 2018). However, in the here studied forest soils increasing potential peptidase activity with temperature did not lead to faster protein depolymerization in the short term, indicating that protein stabilization/occlusion was the main control, likely fostered by protein sorption and by complexation of SOM and of proteins with tannins or metal ions (Nierop et al., 2002;Joanisse et al., 2007;Hu et al., 2018).

Soil moisture effects on soil organic N cycling
Soil moisture had no distinct effect on gross protein depolymerization and microbial amino acid uptake rates ( Fig. 4a and b). Non-significant changes in gross protein depolymerization rates between 30 and 60% WHC indicate that diffusion of proteins and enzymes was not the limiting factor within the studied soil moisture range (30-60% WHC) which most likely also does not translate into a stimulation of depolymerization at higher soil moisture i.e. 90% WHC. Most notable, we observed strong inverse responses of protein depolymerization rates to suboxic conditions in calcareous and silicate soils, i.e. increases in calcareous soils and decreases in silicate soils. To explain these adverse suboxic effects on protein depolymerization the following mechanisms might be relevant: (i) energy limitation induced by lower energetic efficiency of anaerobic compared to aerobic respiration and subsequent reduction of microbial N mining (Kögel-Knabner et al., 2010), (ii) soil pH increases due to redox reactions with alternate electron acceptors such as NO 3 − , Mn 4+ , Fe 3+ , SO 4 2− , and CO 2 (Thompson et al., 2006;Grybos et al., 2009), (iii) release of adsorbed organic matter including proteins due to reductive dissolution of Fe and Mn oxyhydroxides (Herndon et al., 2017), and (iv) release of Al 3+ with negative effects on protein availability and enzyme conformation (Garcidueñas Piña and Cervantes, 1996;Nierop et al., 2002).
(i) Energy limitation: The first possible relevant mechanism relies on the lower energy efficiency of anaerobic metabolism and decomposition processes, eventually causing energy limitation of the microbial communities. This might be followed by decreases in growth and respiration, in extracellular enzyme production for N mining and inhibition of protein depolymerization, while utilization of FAA as C source might increase due to energy limitation. Though intriguing, and important for long-term responses, we found little support for the generality of this mechanism to explain changes in protein depolymerization in the short-term here. The increase of TFAA concentrations in all soils during preincubation at suboxic conditions (90% WHC, Fig. 4d) indicates an imbalance between production rates of amino acids by protein depolymerization and microbial amino acid utilization, and provides direct evidence that amino acids were not increasingly used as a C source. Moreover, in soils from both parent materials, microbial respiration and microbial biomass remained unaltered under suboxic conditions while microbial growth and microbial C use efficiency as reported by Zheng et al. (2019) in associated experiments declined strongly.

(ii)
Increasing pH: The second mechanism is based on the increase of soil pH by 0.3-2 units due to denitrification (Grybos et al., 2009) and reduction of Fe(III)and Mn(IV)-oxyhydroxides (Thompson et al., 2006), both effectively consuming H + . We found an almost complete depletion of the soil NO 3 − pool by dissimilatory nitrate reduction after 7 days at 1% O 2 , indicating increases in soil pH in all soils. Increasing pH negatively affects the sorption capacity of organic matter to Fe-oxides (Gu et al., 1994;Kaiser and Zech, 1999) and this is assumed to be a further control on organic matter desorption under anoxic conditions besides reductive dissolution of Fe oxyhydroxides (Grybos et al., 2009). Optimal pH for sorption of proteins on clay minerals depends on the isoelectric point (IEP) of the adsorbed protein i.e. the pH where they have zero net charge. In laboratory experiments adsorption of albumin showed an optimum close to the IEP, decreasing above that due to repulsion by clays of more negatively charged proteins (Quiquampoix and Ratcliffe, 1992;Staunton and Quiquampoix, 1994).
Since Fe-oxyhydroxides are the main sorption sites of organic matter in soils (Kaiser and Guggenberger, 2003), desorption of proteins due to rising pH might explain the observed increase of protein depolymerization in calcareous soils, but not the patterns found in silicate soils.
(iii) Fe-Mn oxide reduction: In the studied soils the NO 3 − pool was exhausted after pre-incubation at 1% O 2 (Fig. S2), forcing heterotrophic microbes from dissimilatory nitrate reduction (e.g. denitrification) to anoxic Mn(IV) and Fe(III) respiration. As Fe(III) is reduced to Fe(II), Fe-oxyhydroxides are dissolved, and thereby strongly sorbed organic matter is released and becomes bioavailable (Herndon et al., 2017), which can enhance the pH-dependent solubilization process of proteins under suboxic conditions. A combination of increasing pH and Fe(III) reduction therefore might explain the stimulation of protein depolymerization in calcareous soils under O 2 limitation. Considering the different parent materials of the studied soils, the specific mineral composition might influence the response of gross depolymerization to anoxia. The studied silicate soils developed on sericite-quartz and phyllite which have a higher content in easily decomposable Al-silicates such as muscovite and chlorites (Keil and Neubauer, 2011). In contrast, the studied calcareous soils developed on limestone and dolomite deposits originating from the northern calcareous alps, with enhanced contents of kaolinite and Fe-(hydr)oxides such as goethite and hematite (Egli et al., 2008). The differences in pedogenesis resulted in higher dithionite-extractable and crystalline Fe in the studied calcareous soils compared to silicate soils and higher oxalate and dithionite extractable Al in silicate soils ( Table 1). The silicate soils showed higher depolymerization rates at the same level of NaOH-extractable protein compared to calcareous soils (Fig. 2), indicating that proteins are more strongly sorbed in calcareous than in silicate soils. This might have also contributed to greater effects of increasing pH and Fe(III) respiration on calcareous soils where protein depolymerization is more strongly constrained by protein sorption than in silicate soils.
(iv) Al 3+ release: In contrast to Fe-oxyhydroxides, poorly crystalline Al phases, the most reactive sites for DOC complexation (Davis and Gloor, 1981;Kaiser, 2003;Kothawala et al., 2009), are not directly affected by redox changes. However, due to reductive dissolution of Fe-oxyhydroxides Al 3+ and at pH > 4.5, Al(OH) 3 are released from binding sites and might inhibit proteolytic enzymes through conformational changes as well as decreased protein availability by complexation and precipitation of Al-humus complexes (Nierop et al., 2002). Protein depolymerization rates may also decrease due to inhibition of soil microorganisms by elevated Al 3+ under acidic conditions (Garcidueñas Piña and Cervantes, 1996;Scheel et al., 2007), such as found in the silicate soils studied here. In the studied soils this negative Al-effect most likely has overridden the stimulating effect of anoxic conditions on protein availability and depolymerization as demonstrated in calcareous cropland and pasture soils. This is strongly supported by the negative correlations of peptidase activity with exchangeable Al 3− and amorphous Al contents (Table 1).

Conclusions
The presented study aimed at deciphering major controls of gross protein depolymerization and microbial amino acid uptake rates in soils considering temperature, soil moisture/O 2 concentration, parent material and land use. Our results showed that cleavage of soil proteins and microbial amino acid uptake are tightly co-regulated and are mainly controlled by substrate availability, resulting in low sensitivity to short-term temperature and moisture changes compared to e.g. N mineralization. The prominent role of substrate availability is also underlined by the relatively small effects of potential peptidase activities on protein depolymerization rates. The different response of protein depolymerization rates to suboxic conditions in calcareous and silicate soils implies that in calcareous soils Fe-oxyhydroxides are the main sorption sites for soil proteins while on the contrary protein stabilization in silicate soils is rather driven by organo-metal complexes. Our findings imply that soil physicochemical properties, such as pH, mineral composition and clay content caused by different soil parent materials are the key controls of organic N cycling over biotic controls such as microbial activity, microbial community structure or potential enzyme activity.

Supplementary Material
Refer to Web version on PubMed Central for supplementary material.
Jilling A, Keiluweit M, Contosta AR, Frey S, Schimel J, Schnecker J, Smith RG, Tiemann L, Grandy Zhang L, Altabet MA. Amino-group-specific natural abundance nitrogen isotope ratio analysis in amino acids. Schematic of potential controls on gross protein depolymerization in soils including effects of parent material (blue), land use (green), and soil temperature and moisture/O 2 (grey). (For interpretation of the references to color in this figure legend, the reader is referred to the Web version of this article.) Noll   a) Relationship between gross protein depolymerization (at 15 °C) and NaOH-extractable protein contents in calcareous and silicate soils. Solid lines represent 2nd polynomial regression model for silicate (r 2 = 0.721, p < 0.001, n = 12) and calcareous soils (r 2 = 0.598, p < 0.010, n = 12). Regression models were significantly different between silicate and calcareous soils (p < 0.001). b) Relationship between gross protein depolymerization (at 15 °C) and the contribution of NaOH-extractable protein-N to total soil N contents in calcareous and silicate soils. Solid lines represent 2nd polynomial regression models for silicate (r 2 = 0.672, p < 0.01, n = 12) and calcareous soils (r 2 = 0.504, p < 0.05, n = 12). Regression models were significantly different between silicate and calcareous soils (p < 0.001, ANOVA). (a) Gross protein depolymerization, (b) microbial amino acid uptake and (c) peptidase activities at 5, 15 and 25 °C (mean ± s.d., n = 4) in 6 soils differing in parent material (calcareous, silicate) and land use (croplands, pastures, forests). Statistical analyses are provided in Table 3. (a) Gross protein depolymerization, (b) microbial amino acid uptake, (c) peptidase activities and (d) free amino acid contents at different moisture/O 2 levels (30% WHC/21% O 2 , 60% WHC/21% O 2 and 90% WHC/1% O 2 ) at 20 °C in 6 soils differing in parent material (calcareous, silicate) and land use (croplands, pastures, forests; mean ± s.d.), n = 4). Statistical analyses are provided in Table 5.  Table 1 Initial soil physicochemical and biological properties (mean ± s.d., n = 4). Effects of soil parent material (P), land use (L) and their interaction (PxL) were tested by two-way ANOVA. Asterisks represent significance levels.  Three-way ANOVA of gross protein depolymerization rates, gross microbial amino acid uptake rates and peptidase activities for 6 soils differing in parent material and land use, exposed to three soil temperature levels, and their respective interactions (n = 4).  Table 5 3-Way ANOVA of gross protein depolymerization rates, gross microbial amino acid uptake rates and peptidase activities for 6 soils differing in parent material and land use, exposed to three soil moisture/O 2 treatments