Hostname: page-component-78c5997874-t5tsf Total loading time: 0 Render date: 2024-11-05T04:15:00.228Z Has data issue: false hasContentIssue false

Soil C:N impacts on soil biological health and consequences on weed control in soybean and corn systems

Published online by Cambridge University Press:  25 March 2024

Maria A. Gannett
Affiliation:
Graduate Student, Cornell University, Ithaca, NY, USA
Aleah L. Butler-Jones
Affiliation:
Graduate Student, Cornell University, Ithaca, NY, USA
Antonio DiTommaso
Affiliation:
Professor, School of Integrative Plant Science, Soil and Crop Science Section, Cornell University, Ithaca, NY, USA
Jed P. Sparks
Affiliation:
Professor, Department of Ecology and Evolutionary Biology, Cornell University, Ithaca, NY, USA
Jenny Kao-Kniffin*
Affiliation:
Associate Professor, School of Integrative Plant Science, Horticulture Section, Cornell University, Ithaca, NY, USA
*
Corresponding author: Jenny Kao-Kniffin; Email: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

Nitrogen availability has an important influence on agricultural weed growth, because many weeds in annual cropping systems are more competitive in high-nitrogen soils. A potential method to control nitrogen availability is through soil carbon amendments, which stimulate soil microbial growth and immobilize nitrogen. Additionally, carbon amendments may alter soil microbial community composition, increase soil biological functioning, and improve soil health. In a 2-yr field experiment in corn (Zea mays L.) and soybean [Glycine max (L.) Merr.], we implemented five amendment treatments to test their ability to alter weed and crop growth through soil nitrogen availability and soil biological functioning. The treatments included: an untreated control, an unamended weed-free control, rye hay adding 3,560 kg C ha−1 and 3,350 kg C ha−1 in 2020 and 2021, respectively, sawdust adding 5,030 kg C ha−1 and 4,350 kg C ha−1 in 2020 and 2021, respectively, and a rye hay and sawdust combined treatment adding 8,590 kg C ha−1 and 7,700 kg C ha−1 in 2020 and 2021, respectively. Each treatment was replicated five times in corn and six times in soybean. Each season, we explored correlations between crop and weed biomass and weed community composition and nitrogen immobilization measured through soil respiration and nitrogen availability. We also explored changes to the soil microbial community composition and soil health as a secondary result of the carbon amendment treatments. Nitrogen availability was lowest in plots treated with the highest C:N amendment. Increasing carbon improved soil health metrics, but the microbial community composition was most affected by the rye hay treatment. Amendments with high C:N reduced weed growth in both soybean and corn plots but only selected for specific weed communities in soybean, leading to improved soybean competitiveness against weeds. In corn, crop growth and weed community composition remained consistent across amendment treatments. Targeted nitrogen immobilization may improve leguminous crop competition in some weed communities as part of an integrated weed management program.

Type
Research Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
© The Author(s), 2024. Published by Cambridge University Press on behalf of Weed Science Society of America

Introduction

An important aspect of weed control is accurate fertilization. In annual cropping systems, the most abundant weeds are nitrophilous, meaning that they grow well in high-nitrogen soils (Costea et al. Reference Costea, Weaver and Tardif2004; Little et al. Reference Little, DiTommaso, Westbrook, Ketterings and Mohler2021; Moreau et al. Reference Moreau, Busset, Matejicek and Munier-Jolain2014). Some of these species can outcompete crops for nitrogen and have higher nitrogen use efficiencies, which can lead to lower crop yields in fertilized soil relative to unfertilized soil at high weed densities (Di Tomaso Reference Di Tomaso1995; Little et al. Reference Little, DiTommaso, Westbrook, Ketterings and Mohler2021). Once a weed accumulates more nutrients than a crop, it may more strongly compete for additional resources, such as light (Di Tomaso Reference Di Tomaso1995). Weeds can be especially competitive over crops with inorganic fertilizer use (Saberali and Mohammadi Reference Saberali and Mohammadi2015).

Previous research has focused on varying the rates, forms, timing, and location of fertilizer applications to reduce nutrient uptake by weed species and enhance the competitive ability of crops over weeds. Studies that vary fertilization rates have produced mixed results depending on the nitrogen responsiveness of the crop and weed species (Blackshaw and Brandt Reference Blackshaw and Brandt2008). In some cases, reduced fertilization increased crop yield by reducing competition from weeds (Blackshaw et al. Reference Blackshaw, Brandt, Janzen, Entz, Grant and Derksen2003; Davis and Liebman Reference Davis and Liebman2001; Wortman et al. Reference Wortman, Davis, Schutte and Lindquist2011). In other cases, there was no effect of nutrient availability on crop yield (Barker et al. Reference Barker, Knezevic, Martin, Walters and Lindquist2006; Wortman et al. Reference Wortman, Davis, Schutte and Lindquist2011). And in still others, crop yield increased with increased fertilization rates, despite increased weed growth (Anderson et al. Reference Anderson, Tanaka, Black and Schweizer1998; Juroszek et al. Reference Juroszek, Drews, Neuhoff and Kopke2004). Little et al. (Reference Little, DiTommaso, Westbrook, Ketterings and Mohler2021) developed a framework to help predict the effects of fertilization on weed competition and crop yield. According to this framework, key sources of variation in the effects of fertilization on crops and weeds include relative responsiveness to added nutrient, relative competitiveness, and relative shading ability under high fertility.

Another strategy to manage the growth of nitrophilous weeds is through nitrogen immobilization. Nitrogen immobilization is a microbially mediated process (Schimel and Bennett Reference Schimel and Bennett2004). Carbon amendments added to the soil increase the growth of the soil microbial community (Zak et al. Reference Zak, Tilman, Parmenter, Rice, Fisher, Vose, Milchunas and Martin1994). As the microbes grow, they also take up available soil nitrogen to maintain their stoichiometric ratio of carbon to nitrogen. This process leaves less nitrogen available to plant roots and may reduce weed growth.

Nitrogen immobilization may contribute to the efficacy of several existing agricultural weed management strategies that utilize high-carbon substrates. Many perennial cropping systems, such as strawberries [Fragaria ×ananassa (Weston) Duchesne ex Rozier ssp. ananassa] (Pritts and Handley Reference Pritts and Handley1998), raspberries (Rubus idaeus L.) (Trinka and Pritts Reference Trinka and Pritts1992), blueberries (Vaccinium corymbosum L.), and apples [Malus domestica (Suckow) Borkh.] (Merwin et al. Reference Merwin, Rosenberger, Engle, Rist and Fargione1995), rely on high-carbon mulches to help manage weeds. Terminated cover crops can also immobilize nitrogen, and nitrogen immobilization may contribute to their weed-suppressive ability (Pittman et al. Reference Pittman, Barney and Flessner2020; Wells et al. Reference Wells, Reberg-Horton, Smith and Grossman2013; Williams et al. Reference Williams, Wells, Dickey, Hu, Maul, Raskin, Reberg-Horton and Mirsky2018). Mechanisms other than nitrogen immobilization, including physical impedance (Creamer et al. Reference Creamer, Bennett, Stinner, Cardina and Regnier1996; Teasdale and Mohler Reference Teasdale and Mohler2000), light quantity and quality (Teasdale Reference Teasdale1993), changing microclimatic conditions (Teasdale et al. Reference Teasdale, Brandsaeter, Calegari and Neto2007), and allelopathy (Scavo and Mauromicale Reference Scavo and Mauromicale2021), certainly contribute to weed suppression by these mulches. Nitrogen immobilization, however, likely plays an important role in all these strategies, and it is worth exploring in isolation. Managing plant growth through nitrogen immobilization has been studied as a tool to restore native plant growth in natural ecosystems invaded by nonnative plants but has not yet been assessed in cropping systems. A more in-depth understanding of how nitrogen immobilization affects crop and weed growth may provide future weed management opportunities in cropping systems.

Adding carbon to soil to induce nitrogen immobilization may also alter the composition of the soil microbial community, soil biological functioning, and soil health (Lal Reference Lal2014, Reference Lal2016). Good soil biological functioning is associated with diverse benefits such as improved aggregate stability, leading to increased water-holding capacity, porosity, and water infiltration; increased nutrient cycling; improved CEC; and increased crop yields (Luo et al. Reference Luo, Li, Friman, Guo, Guo, Shen and Ling2018; Thangarajan et al. Reference Thangarajan, Bolan, Tian, Naidu and Kunhikrishnan2013). The specific microbial composition can also play an important role in modulating plant growth (Chaparro et al. Reference Chaparro, Sheflin, Manter and Vivanco2012). Many bacterial and fungal endophytes, saprophytes, hyperparasites, and arbuscular mycorrhizal fungi can induce plant defense mechanisms, helping protect plants from disease (Raaijmakers et al. Reference Raaijmakers, Paulitz, Steinberg, Alabouvette and Moënne-Loccoz2009). Some microbes are microbial antagonists or microbial parasites, which help protect plants from pathogenic microbes. Many Trichoderma, Gliocladium, and nonpathogenic Fusarium (all Ascomycota) are antagonistic fungi (Harman et al. Reference Harman, Petzoldt, Comis and Chen2004; Raaijmakers et al. Reference Raaijmakers, Paulitz, Steinberg, Alabouvette and Moënne-Loccoz2009). Many of the Pseudomonas (Proteobacteria), Burkholderia (Proteobacteria), Bacillus (Firmicutes), and Actinobacteria are antagonistic bacteria (Barka et al. Reference Barka, Vatsa, Sanchez, Gaveau-Vaillant, Jacquard, Klenk, Clément, Ouhdouch and van Wezel2015; Raaijmakers et al. Reference Raaijmakers, Paulitz, Steinberg, Alabouvette and Moënne-Loccoz2009). Bacillus can also improve macronutrient cycling by solubilizing phosphorous-containing potassium compounds (Song et al. Reference Song, Feng, Brooks, Zhou, Zhou, Zhao and Li2019; Wang et al. Reference Wang, Zhang, Sun, Qi, Chen and Zhao2017). It is likely that stimulating nitrogen immobilization with high carbon amendments will alter soil microbial community composition and soil health functioning, but it is unclear how those changes will affect weed management.

There were two main objectives of this research. The first objective was to isolate and assess the efficacy of nitrogen immobilization as a tool for managing weeds in different field cropping systems. We hypothesized that amendments with higher carbon:nitrogen ratios (C:N) added to the soil would increase nitrogen immobilization and reduce weed growth to a greater extent than lower C:N. This suppression would be greater for nitrophilous weeds, thus altering weed community composition. We predicted that soybean [Glycine max (L.) Merr.] would be more competitive in nitrogen-immobilized soils, as soybeans can fix atmospheric nitrogen and are less reliant on soil nitrogen availability. We predicted that corn (Zea mays L.) growth would decrease as nitrogen immobilization increased, because modern corn varieties are highly responsive to nitrogen (Liu et al. Reference Liu, Hao, Sha, Huang, Guo, Ke, Chen, Yuan and Mi2022) and have higher recommended nitrogen rates than many other crops (Mylavarapu et al. Reference Mylavarapu, Wright and Kidder2021). The second objective of this study was to explore how high rates of carbon addition would change soil microbial community composition and soil health functioning. We hypothesized that higher carbon inputs would stimulate greater microbial growth and activity and alter microbial community composition. We predicted that these compositional changes would be correlated with improved soil health and function and may be a multifunctional benefit associated with nitrogen immobilization.

Material and Methods

Field Set-Up

Field trials were initiated in the spring of 2020 at the Homer C. Thompson Vegetable Research Farm in Freeville, NY, USA (42.516965°N, 76.334245°W). The soil is classified as Howard gravelly loam (loamy-skeletal, mixed, active, mesic Glossic Hapludalfs) (https://websoilsurvey.nrcs.usda.gov/app/). Baseline soil samples were submitted to the Cornell Nutrient Analysis Laboratory at the beginning of the study, and nutrient levels were within recommended ranges for corn and soybean (Supplementary Table S1). The field used for both corn and soybean plots was previously planted to spinach (Spinacia oleracea L.) in 2019 and was seeded with winter rye (Secale cereale L.) in autumn of 2019 and 2020. In early May 2020 and 2021, rye was chopped at the transition phase (about 15 to 20 cm height) and removed from the field, and then the field was plowed to prepare for treatments and planting. This rye hay, in addition to rye planted and cut at the same time in a neighboring field, was used for the rye hay amendment. Sawdust was collected from a local sawmill (42.2658°N, 76.6461°W) and consisted of about 90% white pine (Pinus strobus L.) and hemlock [Tsuga canadensis (L.) Carrière]. The amendments were spread and incorporated into the soil within 3 wk of rye harvest, and therefore the rye hay was still green and wet when used as an amendment (77% and 57% moisture in 2020 and 2021, respectively). Soybean and corn were planted separately. For each crop, plots were arranged in a complete randomized block design. Within each block, there were five treatments: sawdust amended, rye hay amended, rye hay and sawdust combined, nontreated control, and unamended weed-free control. Soil amendment rates and carbon and nitrogen characteristics are listed in Table 1. Amendments were spread starting mid-May, rototilled to about 15-cm depth, and then soybean and corn were planted on May 27 in 2020 and May 24 in 2021 in rows oriented north to south. Rainfall was below average in 2020, so supplemental irrigation was provided as needed. In 2021, rainfall was above average, so no supplemental irrigation was used (Supplementary Figure S1).

Table 1. Carbon (C) and nitrogen (N) content of amendments added to the soil in 2020 and 2021 (Homer C. Thompson Vegetable Research Farm in Freeville, NY, USA)

a To convert to g m−2, multiply by 0.1.

Soybeans were planted at 76-cm between-row and 4-cm within-row spacing. The soybean planting was divided into six blocks. Each block contained all five amendment treatments. Soybean plots were 6-m wide and 4-m long, contained eight crop rows, and included a split-plot treatment in which half of each plot (four crop rows) was inoculated with rhizobia and half remained uninoculated. The soybean variety was 01072245 Roundup Ready 2 Xtend®, Roundup Ready 2 Yield® Technology (Asgrow® AG20X9, Bayer, Whippany, NJ), which was not preinoculated. Soybeans planted in inoculated plots were inoculated with rhizobia at a rate of 9.4 g kg−1 seed (N-Dure™ Soybean, Verdesian, Cary, NC) about 10 min before planting. However, at harvest in 2020, we dug three soybean plants from each plot, washed the roots, and counted root nodulation. There was no difference in nodulation rate by inoculation treatment (df = 1, F-value = 0.004, P-value = 0.95). This finding indicates that there was an endemic population of Bradyrhizobium japonicum in the soil that successfully colonized the non-inoculated soybeans, reducing the effect of the inoculation treatment. Therefore, the split-plot inoculation treatment was sometimes removed from analysis when it was not significant. Removing the split-plot treatment effectively increased the number of soybean crop replicates to 12 for several of the analyses.

There were five blocks of silage corn, and each block contained all five amendment treatments for a total of 25 plots. Each plot was 3-m wide and 4-m long to accommodate four crop rows. The corn was a hybrid blend of 95% ‘FS 46R64VT2P’ and 5% ‘VNS RR’ (VTDouble PRO Rib Complete, Monsanto, St Louis, MO). Corn was planted at 76-cm between-row and 20-cm within-row spacing. At planting, 224.2 kg ha−1 of nitrogen from 10-20-20 banded fertilizer (plant blended, Phelps Supply, Phelps, NY) was applied at a depth of 5 cm, which is on the high end of regional application rates. Although the goal was to immobilize nitrogen, corn is known to require substantial amounts of nitrogen fertilization. We hoped that band fertilizing could provide adequate nutrients to corn plants, while still allowing for nitrogen immobilization to suppress weeds between rows.

Weed Communities

Seeds from four weed species were scattered in each plot in early June 2020 (not in 2021): Powell amaranth (Amaranthus powellii S. Watson), velvetleaf (Abutilon theophrasti Medik.), white clover (Trifolium repens L.), and ivyleaf morningglory (Ipomoea hederacea Jacq.). Seeds had all been collected from local farms and were cold stratified for at least 1 yr at 4 C. This is somewhat late for seedbank spiking; however, we wanted to wait until after plots had been rototilled so as not to bury the seeds too deep within the soil profile. Preweighed seeds from all four species (Supplementary Table S2) were combined into a screw-top bottle (HDPE Plastic Wide Mouth, Thermo Scientific Nalgene, Rochester, NY), mixed with 120 ml of paver sand (Multi-purpose Patio/Paver Sand, Pavestone, Atlanta, GA) to consistently distribute the different-sized seeds within the bottle, and then distributed evenly across each plot. All plots were raked to improve seed–soil contact. Amaranthus powellii and A. theophrasti were expected to have reduced growth in carbon-amended soils, as these species are known to be responsive to nitrogen (Mohler et al. Reference Mohler, Teasdale and DiTommaso2021). Trifolium repens was not expected to show a strong response, as it is a legume and may be less responsive to nitrogen. We were unsure of the response of I. hederacea, which has a high Ellenberg N index, but which showed a reduced response to low nitrogen soils in previous research in our lab. This seedbank spike, however, was largely unsuccessful. All species made up 2.6% or less of the total weed biomass each year, except A. powellii, which had resident populations (Supplementary Table S2).

Weed surveys were conducted in early August of each year of the experiment. A 0.25-m2 quadrat was placed 1 m into the third crop row of each plot from the southeast corner. The quadrats were 0.33-m wide and 0.75-m long and straddled the crop row, with 0.375 m extending into the interrow space. In each quadrat, all weeds with at least one true leaf were identified to species and clipped at the soil surface, dried at 60 C to a constant weight, and then weighed. Only weeds rooted within the quadrat were counted.

Field Measurements

Soil nitrogen availability was measured using a method adapted from the Kellogg Biological Station Long-Term Ecological Research program (http://lter.kbs.msu.edu/protocols/105). Anion and cation resin exchange membranes (Membranes International, Ringwood, NJ) were cut into strips 10-cm long and 2.5-cm wide. To charge the exchange strips, they were soaked in a 0.5 M HCl bath for 1 h, with stirring every 10 min. Then, strips were rinsed with ultrapure water and soaked in a 0.5 M NaHCO3 bath for 5 h, which was changed every hour. Strips were rinsed with ultrapure water a final time and stored at 4 C in sealed plastic bags to prevent drying. Anion and cation strips were charged and stored separately.

Strips were deployed in the field for periods of 3 wk. In 2020, two measurements were taken in late June and early August. In 2021, four measurements were taken in late June, mid-July, early August, and late August. To best measure root available nitrogen, strips were buried vertically 9.5-cm deep and approximately 15 cm on either side of the crops in the between-row space. Care was taken to reduce soil disturbance and to ensure good strip–soil contact. After collection from the field, strips were rinsed with ultrapure water to remove visible soil particles and stored at 4 C until extraction. Any strips that had been partially eaten by rodents in the field were measured to determine the area of strip missing.

Soil respiration was measured using a soil gas flux survey chamber and CO2 analyzer (Smart Chamber and LI-870 CO2/H2O Analyzer, Li-Cor, Lincoln, NE). Soil collars were made from SDR35 PVC pipe, 20 cm in diameter and 10-cm tall. Collars were installed in the field in late July 2020 and early June 2021, leaving about 1 to 3 cm of the collar above the soil surface. The average offset height of each collar was integrated into the measurement of respiration per soil volume over 120 s. Respiration was measured twice in 2020 (late July and August) and four times in 2021 (mid-June, July, August, and early September). Measurements were taken on consecutive mornings for corn and soybean, at least 5 d after collar installation and at least 24 h after removing any weeds from the collar.

Crop Analyses

Corn and soybean samples were harvested in late September to early October 2020 and late September 2021. Samples were collected from the third crop row of each plot relative to the southeast corner, the row next to which weeds had been collected. Within a 2-m section in the middle of the row, all aboveground crop biomass was harvested, counted, and then dried at 60 C to a constant weight. Soybean pods were shelled by hand in 2020 and with a thresher in 2021 (LD 350, Wintersteiger, Innkreis, Austria) to measure seed yield. The entire corn stalk was weighed together. In 2020, the roots of three soybean plants within each harvested row were removed from the soil, carefully washed, and photographed. Root nodules were counted from the photographs.

Lab Analyses

Before harvest, two soil samples were collected and homogenized from each corn and soybean plot. Samples were collected from the between-row space of the two middle rows, using a 2-cm-wide and 13-cm-long soil core sampler (AMS Inc., Americana Falls, ID). One sample from each plot was collected and discarded before keeping samples to “rinse” the core. Samples were stored on ice in the field and then transferred and stored at 3 C (Setpoint temperature control model SP-322, Goldline, North Kingstown, RI) until processed.

Soil was sieved through a 2-mm sieve and further divided for future analyses. One aliquot was frozen at −20 C to be used for DNA extraction, while another aliquot was set out to air dry. Sieves were rinsed with water, wiped with ethanol, and dried between samples. Gravimetric soil moisture was calculated by weighing a 10-g subsample of fresh soil into a tin cap and drying at 105 C in an oven (Isotemp Oven, Fisher Scientific, Suwanee, GA) for 24 h before reweighing.

To measure pH, 3.0 ± 0.1 g of air-dried soil was mixed with 6 ml of deionized water in a 15 ml falcon tube and shaken on a reciprocal shaker (E6000 Medium-Duty, Eberbach, Ann Arbor, MI) at 180 rpm for 60 min. Tubes were removed from the shaker and allowed to sit for at least 10 min before pH measurements were taken with a pH meter (Orion Star A215, Thermo Scientific, Waltham, MA).

Soil carbon and nitrogen were measured using air-dried and finely ground soil. The grinding apparatus was rinsed with ethanol between samples to prevent cross-contamination. About 30 mg of powdered soil was weighed into tin capsules and submitted to the Cornell Stable Isotope Lab for percent carbon (C), nitrogen (N), 13C, and 15N analysis on an isotope ratio mass spectrometer (Delta V, Thermo Scientific) connected to an elemental analyzer (NC2500, Thermo Scientific).

Both anion and cation resin strips were mixed with 70 ml of 1 M KCl in an acid-washed sample cup. Cups were shaken for 24 and 22 h in 2020 and 2021, respectively, at 180 rpm on a reciprocal shaker (MaxQ 4000, Thermo Scientific) then decanted into 15-ml falcon tubes. Collected extract from all samples was stored frozen at −20 C until analysis. Nitrate and ammonium were quantified using the colorimetric method described by Hood-Nowotny et al. (Reference Hood-Nowotny, Umana, Inselbacher, Oswald- Lachouani and Wanek2010). Samples were analyzed on a microplate reader (Synergy HT, Bio-Tek Instruments, Winooski, VT) and compared with a standard curve to calculate concentration.

In 2021, subsets of soil collected from the non-inoculated soybean and corn plots were submitted to the Cornell Soil Health Lab for a soil health assessment. Soil health indicator tests included: aggregate stability, organic matter, autoclaved-citrate extractable (ACE) soil protein index, active carbon, and macro- and micronutrients. Depth to compaction was also measured by recording the depth at which 21.1 kg cm−2 of penetration resistance was reached using a penetrometer (AgraTronix Soil Compaction Tester, Streetsboro, OH) (Duiker Reference Duiker2002).

Microbial Amplification

Bacterial and fungal community composition were assessed using high-throughput amplicon sequencing. DNA was extracted from approximately 200 mg of soil using the Qiagen DNeasy soil DNA extraction kit (Beverly, MA) following the manufacturer’s instructions. The V3-V4 region of the 16S rRNA gene (16S) was targeted using 341F and 805R primers for bacterial identification (Herlemann et al. Reference Herlemann, Labrenz, Jürgens, Bertilsson, Waniek and Andersson2011). The ITS2 region of the fungal internal transcribed spacer gene (ITS) region was targeted using the ITS1F and 58A2R primers for fungal identification (Gardes and Bruns Reference Gardes and Bruns1993; Martin and Rygiewicz Reference Martin and Rygiewicz2005). Amplified samples were cleaned, indexed, and pooled following the procedure detailed by Garcia et al. (Reference Garcia, Gannett, Wei, Cheng, Hu, Sparks, Giovannoni and Kao-Kniffin2022). Pooled samples were sequenced on the Illumina MiSeq at the Cornell Genomics Facility (Ithaca, NY) using a 500-cycle MiSeq Reagent Kit v. 2 for the ITS pool and a 600-cycle MiSeq Reagent Kit v. 3 for the 16S pool. The raw sequencing data for both 16S and ITS regions of the soil microbial community were deposited into the National Center for Biotechnology Information Sequence Read Archive under BioProject ID: PRJNA1005696.

Sequence merging and filtering were performed in Qiime2. Reads were demultiplexed, and paired ends were merged and trimmed of primers using q2-demux. Sequences were then denoised into amplicon sequence variants (ASVs) using DADA2 (Callahan et al. Reference Callahan, McMurdie, Rosen, Han, Johnson and Holmes2016). The ASVs were clustered using de novo q2-vsearch, with 97% shared identity as the cutoff, to create a table of assigned operational taxonomic units (OTUs). The OTU table was used to create a feature table using q2-feature-table. Bacterial and fungal classifiers were trained with a Naïve Bayes classifier trained on the 99% Silva 138 database for bacteria (Bokulich et al. Reference Bokulich, Kaehler, Rideout, Dillon, Bolyen, Knight, Huttley and Caporaso2018; Robeson et al. Reference Robeson, O’Rourke, Kaehler, Ziemski, Dillon, Foster and Bokulich2021) and the UNITE 13.8 database for fungi (Abarenkov et al. Reference Abarenkov, Henrik Nilsson, Larsson, Alexander, Eberhardt, Erland, Høiland, Kjøller, Larsson, Pennanen, Sen, Taylor, Tedersoo, Ursing and Vrålstad2010), and then taxonomy was assigned to sequences in the taxonomy table using classify-sklearn.

Statistical Analyses

All statistical analyses were performed in R v. 4.1.2 (R Core Team 2022). Corn and soybean yields were analyzed separately using linear mixed-effects models (function lmer, in package lme4 v. 1.1-31) with soil amendment treatment, year, weed biomass, and the interaction of weed biomass and amendment as fixed effects. The model for soybean yield also included an inoculation treatment and the interaction of inoculation and amendment. Block was always included as a random effect. The correlation between soybean seed yield and soybean aboveground biomass was calculated using the Spearman method. The same fixed and random effects, except weed biomass, were included in models of total weed biomass across all plots. The crop biomass per meter of crop row and the weed biomass per square meter were summed, and then percent crop biomass of total biomass was calculated. These weed and crop biomass data were collected at different times in the growing season, so the percent crop biomass measure is intended as a coarse way to compare the relative growth of crop and weeds but is not a measure of the carrying capacity of the soil. The crop biomass, weed biomass, and percent crop biomass of total plant growth were analyzed separately for each year, using a linear mixed-effects model that included all the aforementioned covariates except year.

We analyzed the weed community composition, using an indicator species analysis (function multipatt, in package indicspecies v. 1.7.12) based on weed biomass per species per plot, excluding weed-free control plots. To assess weed species traits, representative Ellenberg N index values, seed masses, and specific leaf areas were obtained for each weed species in our data set. Data were preferentially taken from Bàrberi et al. (Reference Bàrberi, Bocci, Carlesi, Armengot, Blanco-Moreno and Sans2018), then from the TRY Plant Trait Database (Kattge et al. Reference Kattge, Bönisch, Díaz, Lavorel, Prentice, Leadley, Tautenhahn, Werner, Aakala, Abedi, Acosta, Adamidis, Adamson, Aiba and Albert2020), then from Mohler et al. (Reference Mohler, Teasdale and DiTommaso2021). In Mohler et al. (Reference Mohler, Teasdale and DiTommaso2021), the Ellenberg N index was not reported but responsiveness to nutrients was categorized as low, moderate, or high. Ellenberg N indexes of 2, 5, and 8 were assigned to those categories, respectively. Specific leaf area of large crabgrass [Digitaria sanguinalis (L.) Scop.] was taken from Garnier et al. (Reference Garnier, Cordonnier, Guillerm and Sonié1997). Representative trait values for each species were multiplied by the number of individuals found in each plot. The resulting values were summed across all weed species, then divided by the total number of weeds found in the plot to calculate the weighted mean trait value for the weed community. Weighted mean trait values were used as dependent variables for linear mixed-effects models, as described earlier.

Soil respiration, nitrate, and ammonium were analyzed for each crop separately. For corn treatments, the amendment was the only fixed effect. For soybean treatments, the fixed effects were amendment, inoculation, and the interaction between amendment and inoculation. Post hoc analyses were performed using Tukey’s HSD (function cld, in package lsmeans v. 2.27-62). Individual soil health indicator test values from the Cornell Soil Health Test were rated on a scale of 0 to 100 based on their desirability (Moebius-Clune et al. Reference Moebius-Clune, Moebius-Clune, Gugino, Idowu, Schindelbeck, Ristow, van Es, Thies, Shayler, McBride, Wolfe and Abawi2016). These ratings were averaged across all indicator tests to create the overall soil health indicator score. This overall soil health score and the raw values (not rated on the 0 to 100 scale) were analyzed using linear mixed-effects models, as described earlier. Raw values were also centered and scaled, and a principal component analysis (PCA) was run (function prcomp, in package stats v. 3.6.2). The first two eigenvectors were plotted to show the dominant relationship between the plots. The rotation data were plotted to show the influence of each soil health indicator test on the overall soil health score.

To prepare data for microbial community analyses, samples with fewer than 1,000 reads were removed. Data were randomly subsampled to give an equal number of reads per sample (function rrarefy, in package vegan v. 2.6-4), and feature counts were converted to percentages. Percent abundance was converted to a Bray-Curtis distance matrix (function vegdist, in package vegan v. 2.5-7). A principal coordinates analysis (PCoA) was performed on the distance matrix (function cmdscale, in package stats v. 4.1.2). A permutational multivariate analysis of variance (PERMANOVA) using 999 permutations was run to determine significant differences between treatments with pseudo P-values (function adonis2, in package vegan v. 2.5-7).

Differential abundances of OTUs based on rye or no-rye amendment were calculated in two ways. The first workflow used a Dirichlet-multinomial model with transformed (centered log-ratio transformation) nonrarefied abundance data (function aldex, in package ALDEx2 v. 1.4.0) (Fernandes et al. Reference Fernandes, Reid, Macklaim, McMurrough, Edgell and Gloor2014). The second workflow detected multivariable associations using linear models of normalized (total-sum scaling) and log-transformed nonrarefied abundance data with a Benjamini-Hochberg correction method to reduce the false discovery rate (function Maaslin2, in package Maaslin2 v. 1.8.0) (Mallick et al. Reference Mallick, Rahnavard, McIver, Ma, Zhang, Nguyen, Tickle, Weingart, Ren, Schwager, Chatterjee, Thompson, Wilkinson, Subramanian and Lu2021). Additionally, a random forest machine learning algorithm was used (function randomForest, in package randomForest v. 4.6-14) to create a predictive model of rye-amended versus no-rye soils. Random forest was run on centered log-ratio values of OTUs with greater than 2% abundance across all the plots.

Results and Discussion

Amendments with High C:N Shift the Competitive Ability of Soybeans over the Weed Community

Soybean Biomass Changed with Soil Amendments

Soybean aboveground biomass was affected by soil amendment treatments (df = 4, F-value = 6.97, P < 0.0001) and year (df = 1, F-value = 6.85, P = 0.01) (Table 2) and was highly correlated with seed yield (R = 0.98). Inoculation treatment, total weed biomass, and their interactions did not affect soybean biomass. All amendment treatments reduced aboveground biomass compared with the weed-free control, except sawdust and rye hay combined in 2020 (df = 4, F-value = 10.06, P < 0.0001) (Figure 1). In 2020 and 2021, the nontreated control and the rye hay–amended soils had the lowest soybean biomass, respectively.

Table 2. P-values, degrees of freedom (df), and F-values of crop and weed aboveground biomass from a type II ANOVA of a linear model for each independent variable included in the model

Figure 1. Soybean crop yield and total weed biomass in soils with different amendments. Bars labeled with different letters indicate significantly different means (P < 0.05) by treatment.

Total Weed Biomass and Weed Species Community Changed with Soil Amendments

Total weed aboveground biomass was also affected by soil amendment treatments (df = 4, F-value = 146.6, P < 0.0001) (Table 2). Rye hay–amended plots consistently had the highest weed biomass. Sawdust-amended plots and sawdust and rye hay combined plots in 2020 did not have higher weed biomass than the weed-free control plots. Rye hay is notable for containing benzoxazinone, which is an allelopathic compound (Du Fall and Solomon Reference Du Fall and Solomon2011; Schulz et al. Reference Schulz, Marocco, Tabaglio, Macias and Molinillo2013). This study did not explore allelopathy, and so it is outside the realm of this discussion. However, we do not believe that allelopathy played a strong role in weed suppression in this experiment, as rye hay–amended plots consistently had the highest weed growth.

Soil Amendments with Lower C:N Ratios Selected for Nitrogen-Responsive Weed Communities

In addition to total weed biomass, the weed community composition changed based on the soil amendment treatments in soybean plots. An indicator species analysis identified several species indicative of soil amendment treatments in 2020 and 2021 (Table 3). These species were identified as having a high specificity, which is a measure of how predictive a species is of a treatment, and/or a high sensitivity of a species to a treatment, which is the estimate of the probability of finding the species in sites belonging to the site group (de Cáceres and Legendre Reference De Cáceres and Legendre2009; Dufrêne and Legendre Reference Dufrêne and Legendre1997). More weed species were identified as indicator species in 2021 than in 2020, suggesting greater species differentiation over time. In 2020, common purslane (Portulaca oleracea L.) had a specificity value of 1.0, indicating it was only found in nontreated control, rye hay, and sawdust and rye hay plots (not sawdust-amended plots). Only Amaranthus spp. were found in all rye hay–amended plots, but all identified species were indicative of rye hay–amended plots. In 2021, this trend held; all identified species were indicative of rye hay–amended plots, but Amaranthus spp. and common lambsquarters (Chenopodium album L.) were indicative of only rye hay–amended plots (but were not found in every rye hay–amended plot). Common chickweed [Stellaria media (L.) Vill.] was only found in rye hay and sawdust and rye hay combined plots. Both witchgrass (Panicum capillare L.) and hairy nightshade (Solanum physalifolium Rusby) were only found in nontreated control, rye hay, and sawdust and rye hay combined plots (not sawdust-amended plots).

Table 3. Significant indicator species identified in plots with different amendment treatments

a High specificity value indicates that species occurs in plots belonging to this treatment only.

b High sensitivity value indicates that species appears in all plots belonging to this treatment.

Along with a change in weed species composition, there was a corresponding change in weed species traits. In 2021, weighted-average weed species traits were affected by soil amendment treatment (Figure 2). The Ellenberg N index (df = 3, F-value = 7.17, P = 0.0007), seed weight (df = 3, F-value = 3.26, P = 0.03), and specific leaf area (df = 3, F-value = 6.14, P = 0.002) all showed a similar trend with the highest value in rye hay–amended plots and the lowest value in sawdust-amended plots, although no treatments were significantly different according to post hoc analysis of seed weight.

Figure 2. Average weed species traits in plots treated with different amendments for corn and soybean plots. Trait values were retrieved from the literature for each species and weighted by the number of that species found in a plot. Box plots labeled with different letters indicate significantly different means (P < 0.05) by treatment. Average weed species traits not significantly affected by soil amendments are not labeled with letters.

Taken together, these data suggest that rye hay, which had a low C:N, quickly selected for specific weed communities. Many of these species, such as Amaranthus spp. and C. album, are highly responsive to nitrogen (Blackshaw et al. Reference Blackshaw, Brandt, Janzen, Entz, Grant and Derksen2003; Costea et al. Reference Costea, Weaver and Tardif2004; Lindsey et al. Reference Lindsey, Warncke, Steinke and Everman2013), and the whole weed community was selected for functional traits associated with r-selected growth: aboveground resource capture and nitrogen responsiveness (Grime Reference Grime1977; Grime and Hunt Reference Grime and Hunt1975; Moreau et al. Reference Moreau, Busset, Matejicek and Munier-Jolain2014). The trend toward higher seed weight in rye hay–amended soils was unexpected, because r-selected species are often associated with smaller seed weights (Booth et al. Reference Booth, Murphy and Swanton2003; Tilman Reference Tilman1985). Low weighted-average Ellenberg N index values and low specific leaf area suggested some corresponding selection of species associated with K-selected growth in sawdust-amended soils, which had the highest C:N. However, selection of K-selected species may be slower, as no indicator species were identified for sawdust-amended plots.

In High Carbon-amended Soils, Percent Soybean Aboveground Biomass Increased

These changes in soybean aboveground biomass and weed community composition and biomass resulted in significant changes in percent soybean biomass of total plant growth measured within a plot based on soil amendment treatments. This percentage does not represent the carrying capacity of the soil, as weed and crop biomass was collected at different times, but may provide an indication of the relative growth of crops and weeds. As the C:N of amendments increased, there was a greater proportion of soybean biomass in both 2020 (df = 4, F-value = 25.96, P < 0.0001) and 2021 (df = 4, F-value = 47.89, P < 0.0001) (Figure 3). This effect was driven both by higher soybean biomass in soils amended with higher C:N amendments, as well as lower weed biomass. Conversely, soils amended with lower C:N amendments always had the highest percentage of weed growth, which was correlated with the shift in weed community traits toward greater response to nitrogen and resource capture.

Figure 3. Percent soybean and corn crop yield relative to the total aboveground plant biomass within a square meter in soils with different amendments. Bars labeled with different letters indicate significantly different means (P < 0.05) by treatment. Years where percent plant biomass was not significantly affected by soil amendments are not labeled with letters.

These data suggest that soybeans had an increased competitive advantage over weeds in soils with high C:N amendments, which is consistent with studies assessing soybean yield in cover crop systems within soils with low available nitrogen (Pittman et al. Reference Pittman, Barney and Flessner2020; Wells et al. Reference Wells, Reberg-Horton, Smith and Grossman2013; Williams et al. Reference Williams, Wells, Dickey, Hu, Maul, Raskin, Reberg-Horton and Mirsky2018). Soybeans can fix up to 300 kg N ha−1 or 95% of their own nitrogen from the atmosphere by forming a symbiotic relationship with rhizobia in the soil (Keyser and Li Reference Keyser and Li1992; Mourtzinis et al. Reference Mourtzinis, Kaur, Orlowski, Shapiro, Lee, Wortmann, Holshouser, Nafziger, Kandel, Niekamp, Ross, Lofton, Vonk, Roozeboom and Thelen2018). This trait can confer a competitive advantage to soybeans over weeds, especially in low-resource soils (van Heemst Reference van Heemst1985). In our study, soil amendment treatment (df = 4, F-value = 2.87, P-value = 0.03), but not rhizobial inoculation (df = 1, F-value = 0.004, P-value = 0.95), affected nodulation rates. The sawdust and rye hay combined treatment had the highest nodulation rates, and weed-free control plots had the lowest nodulation rates. This pattern was correlated more with soil respiration and percent soil carbon than with nitrogen availability. Soil respiration and percent soil carbon were also highest for the sawdust and rye hay combined plots and lowest in the weed-free control plots. Many factors affect biological nitrogen fixation, but it is well established that nitrogen fixation increases as soil available nitrogen rates decrease (Keyser and Li Reference Keyser and Li1992; Siczek and Lipiec Reference Siczek and Lipiec2011). From our data, it is likely that both percent soil carbon and soil nitrogen availability could play a role in modulating soybean nodulation rates.

Amendments Did Not Shift the Competitive Ability of Corn

Corn Biomass Did Not Change with Soil Amendments

Total corn biomass was not significantly affected by any independent variable (Table 2) except when analyzed separately by year. Sawdust and rye hay combined plots reduced corn biomass relative to the weed-free control in 2020 (df = 4, F-value = 5.28, p = 0.007) (Figure 4). However, total weed biomass in corn plots followed the same trend as weeds in soybean plots, consistently varying based on amendment treatment (df = 4, F-value = 68.3, P < 0.0001). In 2020, all amendment treatments had increased weed biomass compared with the weed-free control, but in 2021, only the rye hay and sawdust and rye hay combined plots had higher weed biomass than the weed-free control.

Figure 4. Corn crop yield and total weed biomass in soils with different amendments. Bars labeled with different letters indicate significantly different means (P < 0.05) by treatment. Years where total plant biomass was not significantly affected by soil amendments are not labeled with letters.

Weed Species Community Composition Was Minimally Changed with Soil Amendments

Despite changes in total weed biomass, there was very little evidence for changes in weed community composition based on the soil amendment treatments in corn plots. An indicator species analysis did identify several species indicative of soil amendment treatments in 2020 but none in 2021 (Table 3). In 2020, D. sanguinalis was found in all sawdust and sawdust and rye hay combined plots. Both shepherd’s purse [Capsella bursa-pastoris (L.) Medik.] and P. oleracea followed the same trends in corn plots as they did in soybean plots. Capsella bursa-pastoris was indicative of nontreated control plots and rye hay–amended plots. Portulaca oleracea was indicative of nontreated control plots, rye hay plots, and sawdust and rye hay combined plots (not sawdust-amended plots). Mirroring this minor change in weed species composition, there was also no change in weed species traits (Figure 2).

Reductions in Total Weed Biomass Drove Changes in Percent Corn Aboveground Biomass

Percent corn biomass of total plant growth was not affected by soil amendments in 2020 (df = 4, F-value = 2.88, P = 0.06), but it was in 2021 (df = 4, F-value = 8.36, P = 0.0008) (Figure 3). Weed-free control plots all had the highest percentage of crop growth, followed by sawdust-amended plots. Rye hay–amended plots consistently had the lowest percentage of crop growth. These data suggest that higher C:N amendments did not increase the competitiveness of corn relative to weeds. As in soybean, the percent corn biomass of the total plant community increased with higher C:N amendments. However, this trend was driven by the changing biomass of the weed community, not changes in weed community composition or in corn growth. Similar to soybean plots, total weed growth decreased in soils with higher C:N amendments, but unlike soybean, corn biomass did not vary based on the amendment.

Only one cornfield and one soybean field were included in our study, so we cannot statically compare the two crops. However, we were able to find evidence that soybeans became more competitive with weeds as soil was subjected to higher C:N amendments, but this trend was not found for corn. Crop type can be one of the strongest filters on weed community composition (Smith and Gross Reference Smith and Gross2007). This has been observed for different species included in intercropped systems (Gomez and Gurevitch Reference Gomez and Gurevitch1998) and even for different cultivars of the same crop (de la Fuente et al. Reference De la Fuente, Suarez and Bhersa2006). However, this filter has often been attributed to differences in cultural management practices such as crop rotations (Sosnoskie et al. Reference Sosnoskie, Herms and Cardina2006) or cover cropping (Smith and Gross Reference Smith and Gross2007). Within our study, most management practices were the same for both corn and soybean plots—field preparation, planting, harvesting, irrigation, and data collection were all done at approximately the same time. However, we did band-fertilize corn at planting, which may have reduced differences in crop growth and/or provided space for nitrogen-responsive weeds to grow in plots where nitrogen was otherwise immobilized. Unfortunately, we did not measure between-row and within-row differences in nitrogen availability or weed community, so we cannot confirm whether this was the case. There was also a difference in planting spacing. Rates were chosen based on local optimal spacing, but the denser seeding rate for soybeans may have resulted in greater changes in the weed community composition. Findings from other studies have shown stronger selection pressure on the weed community in soybean fields compared with cornfields (Culpepper Reference Culpepper2006; Rauber et al. Reference Rauber, Demaría, Jobbágy, Arroyo and Poggio2018).

These findings highlight the potential to utilize reverse fertilization in agriculture. In a highly fertilized field, nitrogen-responsive weeds will eventually dominate the weed community. Reverse fertilization temporarily limits nitrogen availability, reducing the competitive advantage of these weeds. Before the weed community has time to shift toward species with greater belowground investment, leguminous crops will have a competitive advantage. This competitive advantage may not occur in nonleguminous crops. However, additional research is needed to determine what effect, if any, band-fertilizing corn in our study had on results.

Indicators of Nitrogen Immobilization in Soils Amended with High C:N Amendments

Soil CO2 was responsive to soil amendment treatments at every sample date except June (soybean soils) and September (corn soils) in 2021 (Table 4). Total respiration pooled across sample dates within a year was also different based on amendment treatment for both corn (df = 5, F-value = 10.4, P = 0.0002 in 2020; df = 5, F-value = 7.6, P = 0.0007 in 2021) and soybean (df = 6, F-value = 13.9, P > 0.0001 in 2020; df = 6, F-value = 8.4, P > 0.0001 in 2021). Respiration was consistently highest for the soil amended with sawdust and rye hay combined (Figure 5), although it was not significantly different from rye hay–amended soil for corn in 2021, sawdust-amended soil for soybean in 2020, or rye hay– or sawdust-amended soil for soybean in 2021. Soil CO2 is a measure of microbial respiration, which is an indicator of soil microbial activity. These consistently high respiration rates in sawdust and rye hay combined plots, which had the greatest total amount of carbon added to the soil, implies that carbon additions are the primary driver of microbial community activity.

Table 4. Mean soil measurements (respiration, nitrate, and ammonium) in soils with different amendment treatments a

a Different letters indicate significantly different means (P < 0.05) by treatment for that sample date.

Figure 5. Mean total soil respiration as affected by amendment treatment in 2020 and 2021. Measurements were taken twice in 2020 (July and August) and four times in 2021 (June, July, August, and September). Bars labeled with different letters indicate significantly different means (P < 0.05) by treatment within year.

Increased microbial respiration did not always correspond with decreased nitrogen availability. Total nitrate measured over each year was significantly different based on amendment treatment for both corn (df = 5, F-value = 3.7, P = 0.03 in 2020; df = 5, F-value = 4.9, P = 0.009 in 2021) and soybean (df = 6, F-value = 3.5, P > 0.01 in 2020; df = 6, F-value = 6.2, P = 0.0004 in 2021). Consistently, total nitrate was lowest in sawdust-amended soil and highest in rye hay–amended soil or unamended weed-free control soil (Figure 6). Total ammonium availability was more variable. Total ammonium was significantly affected by amendment treatment in corn in 2021 (df = 5, F-value = 5.5, P = 0.005) and soybean in 2020 (df = 6, F-value = 4.6, P = 0.003). It was numerically lowest in sawdust-amended soil in 2020, but it was numerically lowest in the nontreated control for corn and in the unamended weed-free control for soybean in 2021. Rye hay–amended soil consistently had the highest ammonium availability, although it was only significantly higher for corn in 2021 and soybean in 2020 (Figure 7).

Figure 6. Mean total nitrate as affected by amendment treatment in 2020 and 2021. Measurements were taken twice in 2020 (July and August) and four times in 2021 (June, July, August, and September). Bars labeled with different letters indicate significantly different means (P < 0.05) by treatment within year.

Figure 7. Mean total ammonium as affected by amendment treatment in 2020 and 2021. Measurements were taken twice in 2020 (July and August) and four times in 2021 (June, July, August, and September). Bars labeled with different letters indicate significantly different means (P < 0.05) by treatment within year.

Although less responsive than soil respiration, both nitrate and ammonium availability were lowest in sawdust-amended plots, indicating that the greatest nitrogen immobilization occurred with the highest C:N amendment. Sawdust and rye hay combined plots, which had the greatest total amount of carbon added to the soil, had greater microbial activity, as indicated by the highest soil respiration rates, but the highest respiration did not correspond to the lowest nitrogen availability. These results suggest that to create a nutrient-limiting environment, it is more important to add proportionally more carbon than nitrogen to the soil. Otherwise, the soil microbes can utilize nitrogen from the amendment source, thus leaving more nitrogen available in the soil. These results are consistent with the literature. Nitrogen immobilization is a microbially mediated process, occurring when nitrogen is replaced by carbon as the limiting nutrient for microbial community growth during substrate decomposition (Barrett and Burke Reference Barrett and Burke2000; Manzoni et al. Reference Manzoni, Jackson, Trofymow and Porporato2008). It is well established that a substrate’s C:N is a critical determinant of the rate and duration of nitrogen immobilization or mineralization processes (Hodge et al. Reference Hodge, Robinson and Fitter2000; Manzoni et al. Reference Manzoni, Trofymow, Jackson and Porporato2010), although there are exceptions to this trend (Hättenschwiler et al. Reference Hättenschwiler, Coq, Barantal and Handa2011).

Future research should focus on how to refine nitrogen immobilization into a precise reverse fertilization tool. In the current study, large quantities of amendments were added to the soil to markedly increase soil carbon, with little regard for decomposition rates. Reverse fertilization, however, may be more effective if little to no nitrogen is available for plants early in the growing season and the soil quickly transitions to nitrogen mineralization once crops have become established. The period of immobilization should correspond with the duration of time that nutrient resources are available in the crop seed, allowing the crop to establish dominance over the weed community before soil nitrogen becomes more readily available. This time frame would be crop specific and may be similar to the critical period of weed control. Understanding immobilization and mineralization rates as well as amendment traits will be crucial to selecting the ideal amendment for each cropping system. Many studies have suggested that organic matter characteristics, in addition to C:N, affect the rate and duration of nitrogen immobilization or mineralization. Lignin content (Aber and Melillo Reference Aber and Melillo1982) and solid-state 13C NMR spectroscopy (Bonanomi et al. Reference Bonanomi, Incerti, Giannino, Mingo, Lanzotti and Mazzoleni2013) have been found to be useful in predicting substrate decomposition rates. Bonanomi et al. (Reference Bonanomi, Sarker, Zotti, Cesarano, Allevato and Mazzoleni2019) suggest a framework dividing amendments into four quadrants based on two scales: carbon complexity/availability and nitrogen content. High C:N amendments with low-complexity carbon may be more appropriate for reverse fertilization.

Microbial Community Composition Was More Closely Linked to Soil Amendments Than to Soil Health Measurements

Carbon Amendment Drives Improvements in Soil Health

One of the main tenets of the soil health paradigm is to recognize the living component of soil and the important functions of soil biota in nutrient retention and cycling (Chaparro et al. Reference Chaparro, Sheflin, Manter and Vivanco2012; Haney et al. Reference Haney, Haney, Smith, Harmel and White2018; Stika Reference Stika2013). Therefore, stimulating microbial growth and nitrogen immobilization with high-carbon amendments may also have multifunctional benefits for improved soil health. Average soil health indicator scores in our study were affected by soil amendment treatment (df = 3, F-value = 5.37, P = 0.003). Sawdust and rye hay combined plots had the highest average soil health scores, followed by sawdust-amended plots, both of which had a high amount of carbon added to the soil (Figure 8A). Unamended control plots had the lowest average soil health scores. Most individual soil health indicator tests contributed to this trend. All individual soil health indicators were significantly affected by soil amendment treatment, except depth to hardpan, pH, aluminum, iron, and sulfur (Table 5). PCA of centered and scaled values of soil health indicator tests demonstrates that the improved health in sawdust and rye hay combined plots cluster together (Figure 8B) and are associated with several biological indicators: organic matter, active carbon, ACE soil protein, percent carbon, and soil respiration. Organic matter, active carbon, and percent carbon are primarily measures of soil carbon. Many biological soil health indicators are measures of soil carbon and carbon cycling (Bünemann et al. Reference Bünemann, Bongiorno, Bai, Creamer, De Deyn, de Goede, Fleskens, Geissen, Kuyper, Mäder, Pulleman, Sukkel, van Groenigen and Brussaard2018; Liptzin et al. Reference Liptzin, Norris, Cappellazzi, Bean, Cope, Greub, Rieke, Tracy, Aberle, Ashworth, Bañuelos Tavarez, Bary, Baumhardt, Borbón Gracia and Brainard2022), because soil carbon is the primary driver of the soil food system and is therefore critically important in increasing soil health (Lal Reference Lal2014, Reference Lal2016). Active carbon is a measure of the potassium permanganate oxidizable carbon, which is a proxy for the food source of microbial communities (Weil et al. Reference Weil, Islam, Stine, Gruver and Samson-Liebig2003). The ACE soil protein index estimates the amount of organically bound, or bioavailable, nitrogen by measuring a range of proteins extracted from soil (Hurisso et al. Reference Hurisso, Moebius-Clune, Culman, Moebius-Clune, Thies and van Es2018).

Figure 8. (A) Average soil health ratings for each plot treated with different amendments. Box plots labeled with different letters indicate significantly different means (P < 0.05) by treatment. (B) Principal component analysis (PCA) of centered and scaled soil health indicator values. Shape indicates crop, size indicates soil health rating, and color indicates amendment treatment. All indicators were significantly affected by treatment (P < 0.05), except indicators outlined in red. ACE, autoclaved-citrate extractable soil protein index.

Table 5. The P-value and mean soil health indicator value by soil amendment treatment a

a Values labeled with different letters indicate significantly different means from other values within the row based on Tukey’s post hoc analysis.

b Autoclaved-citrate extractable (ACE) soil protein index.

Improved soil health ratings in our study were also associated with chemical soil health indicators: phosphorus, potassium, calcium, copper, magnesium, manganese, and zinc. Increases in macro- and micronutrient availability may have come from the added sawdust or rye hay (Sadeghpour et al. Reference Sadeghpour, Adeyemi, Hunter, Luo and Armstrong2021) or increased microbial activity may have improved access to soil nutrients (Song et al. Reference Song, Feng, Brooks, Zhou, Zhou, Zhao and Li2019). Increased soil moisture was also associated with improved average soil health ratings. Soil moisture may have increased due to the increased water-holding capacity in soils with added organic matter (Thangarajan et al. Reference Thangarajan, Bolan, Tian, Naidu and Kunhikrishnan2013). These improved soil biological, chemical, and physical health indicators were well correlated with the amount of total carbon added to the soil.

Microbial Community Composition Was Most Closely Linked to Rye Hay Amendments

In both corn and soybean plots, microbial beta diversity clustered by treatment in a PCoA (Figure 9). The PERMANOVA results indicate significant associations between the microbial community and year (P-value = 0.001 for bacteria; P-value = 0.001 for fungi), soil amendment treatment (P-value = 0.001 for bacteria; P-value = 0.001 for fungi), and crop (P-value = 0.005 for bacteria; P-value = 0.001 for fungi). When beta diversity is analyzed separately for each year, or for only the soils where soil health was measured, soil amendment treatment and crop remain significant (Table 6).

Figure 9. Principal coordinates analysis (PCoA) of the bacterial (A) and fungal (B) communities based on Illumina MiSeq amplicon sequencing. Shape indicates crop, color indicates amendment treatment, and ellipse line type indicates year. Permutational multivariate analysis of variance (PERMANOVA) revealed significant effects of soil amendment (pseudo P-value = 0.001 for bacteria and fungi), crop (pseudo P-value = 0.005 for bacteria; pseudo P-value = 0.001 for fungi), and year (pseudo P-value = 0.001 for bacteria and fungi) on microbial community composition.

Table 6. The P-values and degrees of freedom (df) for bacterial (16S rRNA region) and fungal (ITS2 region) beta diversity, measured with a permutational multivariate analysis of variance (PERMANOVA) on a principal coordinates analysis (PCoA) of a Bray-Curtis distance matrix of operational taxonomic units (OTUs)

Community composition clustered most strongly around the rye hay amendment. Comparisons between rye hay–amended plots (including both rye hay plots and sawdust and rye hay combined plots) and no-rye plots (including both unamended controls and sawdust plots) showed clearer patterns than comparisons among all five amendment treatments. The random forest algorithm accurately divides plots into rye hay–amended or no–rye hay plots using both bacterial and fungal communities. The out of box (OOB) error rates for bacterial communities were 5.41% in 2020 and 3.33% in 2021. The OOB error rates for fungal communities were 8.24% in 2020 and 12% in 2021. These error rates are much smaller than the error rates using all five amendment treatments as categories (55.41% in 2020 and 51.67% in 2021 for bacteria; 24.71% in 2020 and 24% in 2021 for fungi).

The results from Maaslin2 support this trend: OTUs with significantly different abundance between rye-amended and no-rye treatments were identified for bacteria and fungi in both years. Percent abundances in rye and no-rye plots are shown in Figure 10 for the 10 OTUs with the lowest q-values, and taxonomic information for these 10 OTUs is presented in Table 7. Despite the low random forest model error rates and the identification of OTUs with differential abundance between rye-amended and no-rye plots according to Maaslin2, the ALDEx2 approach did not identify any bacterial or fungal OTUs as differentially abundant. ALDEx2 is more conservative in measuring differential abundance than Maaslin2; however, Maaslin2 is more conservative than other methods (Nearing et al. Reference Nearing, Douglas, Hayes, MacDonald, Desai, Allward, Jones, Wright, Dhanani, Comeau and Langille2022). Combining methods (Maaslin2 and random forest) can increase confidence in the conclusion that there is a strong difference in microbial community composition between rye and no-rye plots.

Figure 10. The 10 bacterial (A) and fungal (B) operational taxonomic units (OTUs) whose abundances differed most significantly between rye-amended and no-rye plots. These OTUs were identified with multivariable associations with linear models (Maaslin2). Maaslin2 was run on nonrarefied abundance data, but percentage relative abundance data are shown. Identified taxa are ordered from most (top) to least (bottom) significantly different. Full taxonomic classification for each identified taxon can be found in Table 7.

Table 7. Top 10 bacterial and fungal operational taxonomic units (OTUs) with significantly different abundance between plots amended with rye hay and no-rye plots, as identified by the Microbiome Multivariable Association with Linear Models 2.0 algorithm

a The six-digit OTU identifiers correspond to the same six-digit OTU identifiers in Figure 9.

Exploring the differentially abundant OTUs based on rye amendment in our experiment, many of the microbes that were more abundant in rye-amended soils appear to be beneficial. Two strains of Stenotrophomonas geniculata were significantly more abundant in rye-amended soils, compared with no-rye soils, in 2020. Stenotrophomonas geniculata has been shown to have plant growth-promoting properties in chick pea (Cicer arietinum L.) nodules and corn (Gopalakrishnan et al. Reference Gopalakrishnan, Srinivas, Prakash, Sathya and Vijayabharathi2015; Singh et al. Reference Singh, Singh, Trivedi, Sahu, Paul, Paul and Saxena2020). Stenotrophomonas geniculata can also decompose complex compounds such as the herbicide paraquat (Wu et al. Reference Wu, Wu, Chen and Wu2020), nicotine (Liu et al. Reference Liu, Wang, Huang, Wang, Nie, Jiang, Li, Liu, Xu and Tang2014), and polybutylene succinate/polylactic acid (Srimalanon et al. Reference Srimalanon, Prapagdee and Sombatsompop2020). Pseudomonas vancouverensis was also more abundant in rye-amended soils and has been shown to improve the stress tolerance of red peppers (Capsicum annuum L.) (Samaddar et al. Reference Samaddar, Chatterjee, Roy Choudhury, Ahmed and Sa2019) and tomatoes (Solanum lycopersicum L.) (Subramanian et al. Reference Subramanian, Mageswari, Kim, Lee and Sa2015). An uncultured Staphylococcus, which was more abundant in rye hay plots than no-rye plots, was the only OTU shared between the 2020 and 2021 lists of the top 10 bacterial OTUs with significantly different abundance. In 2021, the three most significantly different bacteria were Bacillus species, which contribute to macronutrient cycling by solubilizing phosphorous-containing potassium compounds (Song et al. Reference Song, Feng, Brooks, Zhou, Zhou, Zhao and Li2019). These three bacteria were more abundant in rye hay–amended soils, although the sixth most significantly different OTU was a Bacillus species more abundant in no-rye soils.

Two fungal OTUs were identified by Maaslin2 as significantly different between rye hay and no-rye plots in both years. Both OTUs were more abundant in rye-amended soils. Papiliotrema laurentii is found in a variety of habitats and can metabolize a variety of carbon substrates. It can control phytopathogenic fungi and improve mycorrhizal colonization, nitrogen retention, and plant growth (de Almeida et al. Reference De Almeida, Ventorim, de Moura Ferreira and da Silveira2022). The other fungus was in the family Nectriaceae, but was not further classified. Trichosporon insectorum, Ascobolus foliicola, and Trichoderma evansii were all more abundant in rye-amended soils in 2021. Trichosporon insectorum is a yeast, often associated with insects, with fungicidal properties (Fuentefria et al. Reference Fuentefria, Suh, Landell, Faganello, Schrank, Vainstein, Blackwell and Valente2008). Ascobolus foliicola is often found on bare ground or decomposing plant matter (Uzun et al. Reference Uzun, Yakar, Karacan and Halil2019), and Trichoderma evansii was isolated from the sapwood of Lophira alata and cola [Cola verticillata (Thonn.) Stapf ex A. Chev.] (Samuels and Ismaiel Reference Samuels and Ismaiel2009).

Rye hay–amended plots did seem to support the growth of beneficial microbial communities, but that change was not directly captured in the soil health indicator tests, which were higher as the total carbon amount added to the soil increased. This finding is consistent with the literature, which shows that it is often difficult to directly connect the whole soil microbiome to soil health functioning (Bünemann et al. Reference Bünemann, Bongiorno, Bai, Creamer, De Deyn, de Goede, Fleskens, Geissen, Kuyper, Mäder, Pulleman, Sukkel, van Groenigen and Brussaard2018). The measures of soil biological health often included in soil health studies are primarily measures of carbon and nitrogen cycling, which are indirect measures of the soil microbial community. It is still expensive to measure the microbial community composition directly, and more research needs to be done to link specific communities to improved soil health. A recent study by Wilhelm et al. (Reference Wilhelm, van Es and Buckley2022) used microbial community composition to predict soil health ratings with almost 70% accuracy. In our study, the total microbial activity was more important for improving soil health ratings and altering plant growth dynamics than the specific microbial community composition.

In summary, this study demonstrates that reverse fertilization can improve the competitive advantage of leguminous crops with multifunctional benefits for soil health. High C:N organic amendments stimulated nitrogen immobilization through microbial growth, and greater amounts of added carbon improved soil health ratings. The soil microbial community composition consistently changed when different amendments were added to the soil; however, these compositional shifts were most closely associated with the presence of rye hay. Although these results are promising, reverse fertilization requires further refinement before it can be adopted by growers. Potential ways to refine this tool include immobilizing nitrogen only in specific locations such as between the crop rows or at the soil surface to inhibit the germination of small-seeded, r-selected weeds. Alternatively, nitrogen may be immobilized only during specific times, such as within the critical period of weed control or for a few years before competitive species begin to dominate. Nitrogen immobilization is an important regulator of plant growth that should be considered when designing weed management strategies in agroecosystems.

Supplementary material

To view supplementary material for this article, please visit https://doi.org/10.1017/wsc.2024.17

Acknowledgments

The authors would like to thank Lynn Johnson for her statistical advice; Steve McKay, Rick Randolph, Micaela Moravek, Adam Frosolone, and Sofia Kashtelyan for their field help; and Zach Butler-Jones for his lab help. We would also like to thank Sophie Westbrook and Connor Lane for edits on previous versions of this article.

Funding

This research was partially funded by the Towards Sustainability Foundation grant number 1458761. This research was supported in part by the intramural research program of the U.S. Department of Agriculture, National Institute of Food and Agriculture, Hatch program 2019-20-166 and a graduate research fellowship from the Foundation for Food and Agriculture Research (FFAR).

Competing interests

No conflicts of interest have been declared.

Footnotes

Associate Editor: William Vencill, University of Georgia

References

Abarenkov, K, Henrik Nilsson, R, Larsson, K-H, Alexander, IJ, Eberhardt, U, Erland, S, Høiland, K, Kjøller, R, Larsson, E, Pennanen, T, Sen, R, Taylor, AFS, Tedersoo, L, Ursing, BM, Vrålstad, T, et al. (2010) The UNITE database for molecular identification of fungi—recent updates and future perspectives. New Phytol 186:281285 10.1111/j.1469-8137.2009.03160.xCrossRefGoogle ScholarPubMed
Aber, JD, Melillo, JM (1982) Nitrogen immobilization in decaying hardwood leaf litter as a function of initial nitrogen and lignin content. Can J Bot 60:22632269 10.1139/b82-277CrossRefGoogle Scholar
Anderson, RL, Tanaka, DL, Black, AL, Schweizer, EE (1998) Weed community and species response to crop rotation, tillage, and nitrogen fertility. Weed Technol 12:531536 10.1017/S0890037X00044262CrossRefGoogle Scholar
Bàrberi, P, Bocci, G, Carlesi, S, Armengot, L, Blanco-Moreno, JM, Sans, FX (2018) Linking species traits to agroecosystem services: a functional analysis of weed communities. Weed Res 58:7688 10.1111/wre.12283CrossRefGoogle Scholar
Barka, EA, Vatsa, P, Sanchez, L, Gaveau-Vaillant, N, Jacquard, C, Klenk, H-P, Clément, C, Ouhdouch, Y, van Wezel, GP (2015) Taxonomy, physiology, and natural products of Actinobacteria. Microbiol Molec Biol Rev 80:143 10.1128/MMBR.00019-15CrossRefGoogle ScholarPubMed
Barker, DC, Knezevic, SZ, Martin, AR, Walters, DT, Lindquist, JL (2006) Effect of nitrogen addition on the comparative productivity of corn and velvetleaf (Abutilon theophrasti). Weed Sci 54:354363 10.1614/WS-05-127R.1CrossRefGoogle Scholar
Barrett, JE, Burke, IC (2000) Potential nitrogen immobilization in grassland soils across a soil organic matter gradient. Soil Biol Biochem 32:17071716 10.1016/S0038-0717(00)00089-4CrossRefGoogle Scholar
Blackshaw, RE, Brandt, RN (2008) Nitrogen fertilizer rate effects on weed competitiveness is species dependent. Weed Sci 56:743747 10.1614/WS-08-065.1CrossRefGoogle Scholar
Blackshaw, RE, Brandt, RN, Janzen, HH, Entz, T, Grant, CA, Derksen, DA (2003) Differential response of weed species to added nitrogen. Weed Sci 51:532539 10.1614/0043-1745(2003)051[0532:DROWST]2.0.CO;2CrossRefGoogle Scholar
Bokulich, NA, Kaehler, BD, Rideout, JR, Dillon, M, Bolyen, E, Knight, R, Huttley, GA, Caporaso, JG (2018) Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2’s Q2-feature-classifier plugin. Microbiome 6:90 10.1186/s40168-018-0470-zCrossRefGoogle ScholarPubMed
Bonanomi, G, Incerti, G, Giannino, F, Mingo, A, Lanzotti, V, Mazzoleni, S (2013) Litter quality assessed by solid state 13C NMR spectroscopy predicts decay rate better than C/N and Lignin/N ratios. Soil Biol Biochem 56:4048 10.1016/j.soilbio.2012.03.003CrossRefGoogle Scholar
Bonanomi, G, Sarker, TC, Zotti, M, Cesarano, G, Allevato, E, Mazzoleni, S (2019) Predicting nitrogen mineralization from organic amendments: beyond C/N ratio by 13C-CPMAS NMR approach. Plant Soil 441:129146 10.1007/s11104-019-04099-6CrossRefGoogle Scholar
Booth, BD, Murphy, SD, Swanton, CJ (2003) Interactions between populations I: competition and allelopathy. Pages 119–137 in Booth BD, Murphy SD, Swanton CJ, eds. Weed Ecology in Natural and Agricultural Systems. Wallingford, UK: CABI.10.1079/9780851995281.0119CrossRefGoogle Scholar
Bünemann, EK, Bongiorno, G, Bai, Z, Creamer, RE, De Deyn, G, de Goede, R, Fleskens, L, Geissen, V, Kuyper, TW, Mäder, P, Pulleman, M, Sukkel, W, van Groenigen, JW, Brussaard, L (2018) Soil quality—a critical review. Soil Biol Biochem 120:105125 10.1016/j.soilbio.2018.01.030CrossRefGoogle Scholar
Callahan, BJ, McMurdie, PJ, Rosen, MJ, Han, AW, Johnson, AJA, Holmes, SP (2016) DADA2: high resolution sample inference from Illumina amplicon data. Nat Methods 13:581583 10.1038/nmeth.3869CrossRefGoogle ScholarPubMed
Chaparro, JM, Sheflin, AM, Manter, DK, Vivanco, JM (2012) Manipulating the soil microbiome to increase soil health and plant fertility. Biol Fertil Soils 48:489499 10.1007/s00374-012-0691-4CrossRefGoogle Scholar
Costea, M, Weaver, SE, Tardif, FJ (2004) The biology of Canadian weeds. 130. Amaranthus retroflexus L., A. powellii S. Watson and A. hybridus L. Can J Plant Sci 84:631–66810.4141/P02-183CrossRefGoogle Scholar
Creamer, NG, Bennett, MA, Stinner, BR, Cardina, J, Regnier, EE (1996) Mechanisms of weed suppression in cover crop-based production systems. HortScience 31:410413 10.21273/HORTSCI.31.3.410CrossRefGoogle Scholar
Culpepper, AS (2006) Glyphosate-induced weed shifts. Weed Technol 20:277281 10.1614/WT-04-155R.1CrossRefGoogle Scholar
Davis, AS, Liebman, M (2001) Nitrogen source influences wild mustard growth and competitive effect on sweet corn. Weed Sci 49:558566 10.1614/0043-1745(2001)049[0558:NSIWMG]2.0.CO;2CrossRefGoogle Scholar
De Almeida, ELM, Ventorim, RZ, de Moura Ferreira, MA, da Silveira, WB (2022) Papiliotrema laurentii: general features and biotechnological applications. Appl Microbiol Biotechnol 106:69636976 10.1007/s00253-022-12208-2CrossRefGoogle ScholarPubMed
De Cáceres, M, Legendre, P (2009) Associations between species and groups of sites: indices and statistical inference. Ecology 90:35663574 10.1890/08-1823.1CrossRefGoogle ScholarPubMed
De la Fuente, EB, Suarez, SA, Bhersa, CM (2006) Soybean weed community composition and richness between 1995 and 2003 in the Rolling Pampas (Argentina). Agric Ecosyst Environ 115:229236 10.1016/j.agee.2006.01.009CrossRefGoogle Scholar
Di Tomaso, JM (1995) Approaches for improving crop competitiveness through the manipulation of fertilization strategies. Weed Sci 43:491497 10.1017/S0043174500081522CrossRefGoogle Scholar
Du Fall, LA, Solomon, PS (2011) Role of cereal secondary metabolites involved in mediating the outcome of plant-pathogen interactions. Metabolites 1:6478 10.3390/metabo1010064CrossRefGoogle ScholarPubMed
Dufrêne, M, Legendre, P (1997) Species assemblages and indicator species: the need for a flexible asymmetrical approach. Ecol Monogr 67:345366 Google Scholar
Duiker, SW (2002) Diagnosing Soil Compaction Using a Penetrometer. University Park, PA: Penn State Extension Agronomy Fact 63. https://extension.psu.edu/diagnosing-soil-compaction-using-a-penetrometer-soil-compaction-tester Accessed: April 1, 2024Google Scholar
Fernandes, AD, Reid, JN, Macklaim, JM, McMurrough, TA, Edgell, DR, Gloor, GB (2014) Unifying the analysis of high-throughput sequencing datasets: characterizing RNA-seq, 16S rRNA gene sequencing and selective growth experiments by compositional data analysis. Microbiome 2:15 10.1186/2049-2618-2-15CrossRefGoogle ScholarPubMed
Fuentefria, AM, Suh, SO, Landell, MF, Faganello, J, Schrank, A, Vainstein, MH, Blackwell, M, Valente, P (2008) Trichosporon insectorum sp. nov., a new anamorphic basidiomycetous killer yeast. Mycol Res 112:9399 10.1016/j.mycres.2007.05.001CrossRefGoogle Scholar
Garcia, J, Gannett, M, Wei, L, Cheng, L, Hu, S, Sparks, J, Giovannoni, J, Kao-Kniffin, J (2022) Selection pressure on the rhizosphere microbiome can alter nitrogen use efficiency and seed yield in Brassica rapa. Commun Biol 5:959 10.1038/s42003-022-03860-5CrossRefGoogle ScholarPubMed
Gardes, M, Bruns, TD (1993) ITS primers with enhanced specificity for basidiomycetes—application to the identification of mycorrhizae and rusts. Mol Ecol 2:113118 10.1111/j.1365-294X.1993.tb00005.xCrossRefGoogle Scholar
Garnier, E, Cordonnier, P, Guillerm, J-L, Sonié, L (1997) Specific leaf area and leaf nitrogen concentration in annual and perennial grass species growing in Mediterranean old-fields. Oecologia 111:490498 10.1007/s004420050262CrossRefGoogle ScholarPubMed
Gomez, P, Gurevitch, J (1998) Weed community responses in a corn-soybean intercrop. Appl Veg Sci 1:281288 10.2307/1478958CrossRefGoogle Scholar
Gopalakrishnan, S, Srinivas, V, Prakash, B, Sathya, A, Vijayabharathi, R (2015) Plant growth-promoting traits of Pseudomonas geniculata isolated from chickpea nodules. 3 Biotech 5:653661 10.1007/s13205-014-0263-4CrossRefGoogle ScholarPubMed
Grime, JP (1977) Evidence for the existence of three primary strategies in plants and its relevance to ecological and evolutionary theory. Am Nat 111:11691194 10.1086/283244CrossRefGoogle Scholar
Grime, JP, Hunt, R (1975) Relative growth-rate: its range and adaptive significance in a local flora. J Ecol 63:393422 10.2307/2258728CrossRefGoogle Scholar
Haney, RL, Haney, EB, Smith, DR, Harmel, RD, White, MJ (2018) The soil health tool—theory and initial broad-scale application. Appl Soil Ecol 125:162168 10.1016/j.apsoil.2017.07.035CrossRefGoogle Scholar
Harman, GE, Petzoldt, R, Comis, A, Chen, J (2004) Interactions between Trichoderma harzianum Strain T22 and maize inbred line Mo17 and effects of these interactions on diseases caused by Pythium ultimum and Colletotrichum graminicola . Phytopathology 94:147153 10.1094/PHYTO.2004.94.2.147CrossRefGoogle ScholarPubMed
Hättenschwiler, S, Coq, S, Barantal, S, Handa, IT (2011) Leaf traits and decomposition in tropical rainforests: revisiting some commonly held views and towards a new hypothesis. New Phytol 189:950965 10.1111/j.1469-8137.2010.03483.xCrossRefGoogle ScholarPubMed
Herlemann, DP, Labrenz, M, Jürgens, K, Bertilsson, S, Waniek, JJ, Andersson, AF (2011) Transitions in bacterial communities along the 2000 km salinity gradient of the Baltic Sea. ISME J 5:15711579 10.1038/ismej.2011.41CrossRefGoogle ScholarPubMed
Hodge, A, Robinson, D, Fitter, A (2000) Are microorganisms more effective than plants at competing for nitrogen? Trends Plant Sci 5:304308 10.1016/S1360-1385(00)01656-3CrossRefGoogle ScholarPubMed
Hood-Nowotny, R, Umana, NH-N, Inselbacher, E, Oswald- Lachouani, P, Wanek, W (2010) Alternative methods for measuring inorganic, organic, and total dissolved nitrogen in soil. Soil Sci Soc Am J 74:10181027 10.2136/sssaj2009.0389CrossRefGoogle Scholar
Hurisso, TT, Moebius-Clune, DJ, Culman, SW, Moebius-Clune, BN, Thies, JE, van Es, HM (2018) Soil protein as a rapid soil health indicator of potentially available organic nitrogen. Agric Environ Lett 3:180006 10.2134/ael2018.02.0006CrossRefGoogle Scholar
Juroszek, P, Drews, S, Neuhoff, D, Kopke, U (2004) Effects of organic fertilizers on the development of weeds and winter wheat. J Plant Dis Prot 19:611618 Google Scholar
Kattge, J, Bönisch, G, Díaz, S, Lavorel, S, Prentice, IC, Leadley, P, Tautenhahn, S, Werner, GDA, Aakala, T, Abedi, M, Acosta, ATR, Adamidis, GC, Adamson, K, Aiba, M, Albert, CH, et al. (2020) TRY plant trait database—enhanced coverage and open access. Glob Change Biol 26:119188 10.1111/gcb.14904CrossRefGoogle ScholarPubMed
Keyser, HH, Li, F (1992) Potential for increasing biological nitrogen fixation in soybean. Pages 119–135 in Ladha JK, George T, Bohlool BB, eds. Biological Nitrogen Fixation for Sustainable Agriculture: Extended Versions of Papers Presented in the Symposium, Role of Biological Nitrogen Fixation in Sustainable Agriculture at the 13th Congress of Soil Science, Kyoto, Japan, 1990, Developments in Plant and Soil Sciences. Dordrecht, Netherlands: SpringerGoogle Scholar
Lal, R (2014) Societal value of soil carbon. J Soil and Water Conserv 69:186A192A 10.2489/jswc.69.6.186ACrossRefGoogle Scholar
Lal, R (2016) Soil health and carbon management. Food Energy Secur 5:212222 10.1002/fes3.96CrossRefGoogle Scholar
Lindsey, LE, Warncke, DD, Steinke, K, Everman, WJ (2013) Fertilizer and population affects nitrogen assimilation of common lambsquarters (Chenopodium album) and redroot pigweed (Amaranthus retroflexus). Weed Sci 61:131–3510.1614/WS-D-12-00094.1CrossRefGoogle Scholar
Liptzin, D, Norris, CE, Cappellazzi, SB, Bean, GM, Cope, M, Greub, KLH, Rieke, EL, Tracy, PW, Aberle, E, Ashworth, A, Bañuelos Tavarez, O, Bary, AI, Baumhardt, RL, Borbón Gracia, A, Brainard, DC, et al. (2022) An evaluation of carbon indicators of soil health in long-term agricultural experiments. Soil Biol Biochem 172:108708 10.1016/j.soilbio.2022.108708CrossRefGoogle Scholar
Little, NG, DiTommaso, A, Westbrook, AS, Ketterings, QM, Mohler, CL (2021) Effects of fertility amendments on weed growth and weed–crop competition: a review. Weed Sci 69:132146 10.1017/wsc.2021.1CrossRefGoogle Scholar
Liu, Y, Wang, L, Huang, K, Wang, W, Nie, X, Jiang, Y, Li, P, Liu, S, Xu, P, Tang, H (2014) Physiological and biochemical characterization of a novel nicotine-degrading bacterium Pseudomonas geniculata N1. PLoS ONE 9:e84399 10.1371/journal.pone.0084399CrossRefGoogle ScholarPubMed
Liu, Z, Hao, Z, Sha, Y, Huang, Y, Guo, W, Ke, L, Chen, F, Yuan, L, Mi, G (2022) High responsiveness of maize grain yield to nitrogen supply is explained by high ear growth rate and efficient ear nitrogen allocation. Field Crops Res 286:108610 10.1016/j.fcr.2022.108610CrossRefGoogle Scholar
Luo, G, Li, L, Friman, V-P, Guo, J, Guo, S, Shen, Q, Ling, N (2018) Organic amendments increase crop yields by improving microbe-mediated soil functioning of agroecosystems: a meta-analysis. Soil Biol Biochem 124:105115 10.1016/j.soilbio.2018.06.002CrossRefGoogle Scholar
Mallick, H, Rahnavard, A, McIver, LJ, Ma, S, Zhang, Y, Nguyen, LH, Tickle, TL, Weingart, G, Ren, B, Schwager, EH, Chatterjee, S, Thompson, KN, Wilkinson, JE, Subramanian, A, Lu, Y, et al. (2021) Multivariable association discovery in population-scale meta-omics studies. PLoS Comput Biol 17:e1009442 10.1371/journal.pcbi.1009442CrossRefGoogle ScholarPubMed
Manzoni, S, Jackson, RB, Trofymow, JA, Porporato, A (2008) The global stoichiometry of litter nitrogen mineralization. Science 321:684686 10.1126/science.1159792CrossRefGoogle ScholarPubMed
Manzoni, S, Trofymow, JA, Jackson, RB, Porporato, A (2010) Stoichiometric controls on carbon, nitrogen, and phosphorus dynamics in decomposing litter. Ecol Monogr 80:89106 10.1890/09-0179.1CrossRefGoogle Scholar
Martin, KJ, Rygiewicz, PT (2005) Fungal-specific PCR primers developed for analysis of the ITS region of environmental DNA extracts. BMC Microbiol 5:28 10.1186/1471-2180-5-28CrossRefGoogle ScholarPubMed
Mylavarapu, RS, Wright, DL, Kidder, G (2021) SL129/SS163: UF/IFAS Standardized Fertilization Recommendations for Agronomic Crops. https://edis.ifas.ufl.edu/publication/SS163. Accessed: April 13, 2023Google Scholar
Merwin, IA, Rosenberger, DA, Engle, CA, Rist, DL, Fargione, M (1995) Comparing mulches, herbicides, and cultivation as orchard groundcover management systems. HortTechnol 5:151158 10.21273/HORTTECH.5.2.151CrossRefGoogle Scholar
Moebius-Clune, BN, Moebius-Clune, DJ, Gugino, BK, Idowu, OJ, Schindelbeck, RR, Ristow, AJ, van Es, HM, Thies, JE, Shayler, HA, McBride, MB, Wolfe, DW, Abawi, GS (2016) Comprehensive Assessment of Soil Health. 3rd ed. Geneva, NY: Cornell University. P 33 Google Scholar
Mohler, CL, Teasdale, JR, DiTommaso, A (2021) Manage Weeds on Your Farm: A Guide to Ecological Strategies. SARE Handbook. 418 p. https://www.sare.org/wp-content/uploads/Manage-Weeds-on-Your-Farm.pdf. Accessed: April 1, 2024Google Scholar
Moreau, D, Busset, H, Matejicek, A, Munier-Jolain, N (2014) The ecophysiological determinants of nitrophily in annual weed species. Weed Res 54:335346 10.1111/wre.12090CrossRefGoogle Scholar
Mourtzinis, S, Kaur, G, Orlowski, JM, Shapiro, CA, Lee, CD, Wortmann, C, Holshouser, D, Nafziger, ED, Kandel, H, Niekamp, J, Ross, WJ, Lofton, J, Vonk, J, Roozeboom, KL, Thelen, KD, et al. (2018) Soybean response to nitrogen application across the United States: a synthesis-analysis. Field Crops Res 215:7482 10.1016/j.fcr.2017.09.035CrossRefGoogle Scholar
Nearing, JT, Douglas, GM, Hayes, MG, MacDonald, J, Desai, DK, Allward, N, Jones, CMA, Wright, RJ, Dhanani, AS, Comeau, AM, Langille, MGI (2022) Microbiome differential abundance methods produce different results across 38 datasets. Nat Commun 13:342 10.1038/s41467-022-28034-zCrossRefGoogle ScholarPubMed
Pittman, KB, Barney, JN, Flessner, ML (2020) Cover crop residue components and their effect on summer annual weed suppression in corn and soybean. Weed Sci 68:301310 10.1017/wsc.2020.16CrossRefGoogle Scholar
Pritts, M, Handley, D, eds (1998) Strawberry Production Guide for the Northeast, Midwest, and Eastern Canada. Ithaca, NY: Natural Resource, Agriculture, and Engineering Service NRAES-88. 162 pGoogle Scholar
Raaijmakers, JM, Paulitz, TC, Steinberg, C, Alabouvette, C, Moënne-Loccoz, Y (2009) The rhizosphere: a playground and battlefield for soilborne pathogens and beneficial microorganisms. Plant Soil 321:341361 10.1007/s11104-008-9568-6CrossRefGoogle Scholar
Rauber, RB, Demaría, MR, Jobbágy, EG, Arroyo, DN, Poggio, SL (2018) Weed communities in semiarid rainfed croplands of central Argentina: Comparison between corn (Zea mays) and soybean (Glycine max) Crops. Weed Sci 66:368–7810.1017/wsc.2017.76CrossRefGoogle Scholar
R Core Team (2022) R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing. https://www.R-project.org Google Scholar
Robeson, MSI, O’Rourke, DR, Kaehler, BD,Ziemski, M, Dillon, MR, Foster, JT, Bokulich, NA (2021) RESCRIPt: Reproducible sequence taxonomy reference database management. PLOS Comput Biol 17:e1009581 10.1371/journal.pcbi.1009581CrossRefGoogle ScholarPubMed
Saberali, SF, Mohammadi, K (2015) Organic amendments application downweight the negative effects of weed competition on the soybean yield. Ecol Eng 82:451458 10.1016/j.ecoleng.2015.05.038CrossRefGoogle Scholar
Sadeghpour, A, Adeyemi, O, Hunter, D, Luo, Y, Armstrong, S (2021) Precision planting impacts on winter cereal rye growth, nutrient uptake, spring soil temperature and adoption cost. Renewable Agric Food Syst 36:328333 10.1017/S1742170520000411CrossRefGoogle Scholar
Samaddar, S, Chatterjee, P, Roy Choudhury, A, Ahmed, S, Sa, T (2019) Interactions between Pseudomonas spp. and their role in improving the red pepper plant growth under salinity stress. Microbiol Res 219:6673 10.1016/j.micres.2018.11.005CrossRefGoogle ScholarPubMed
Samuels, GJ, Ismaiel, A (2009) Trichoderma evansii and T. lieckfeldtiae: two new T. hamatum-like species. Mycologia 101:142156 10.3852/08-161CrossRefGoogle Scholar
Scavo, A, Mauromicale, G (2021) Crop allelopathy for sustainable weed management in agroecosystems: knowing the present with a view to the future. Agronomy 11:2104 10.3390/agronomy11112104CrossRefGoogle Scholar
Schimel, JP, Bennett, J (2004) Nitrogen mineralization: challenges of a changing paradigm. Ecology 85:591602 10.1890/03-8002CrossRefGoogle Scholar
Schulz, M, Marocco, A, Tabaglio, V, Macias, FA, Molinillo, JMG (2013) Benzoxazinoids in rye allelopathy—from discovery to application in sustainable weed control and organic farming. J Chem Ecol 39:154174 10.1007/s10886-013-0235-xCrossRefGoogle Scholar
Siczek, A, Lipiec, J (2011) Soybean nodulation and nitrogen fixation in response to soil compaction and surface straw mulching. Soil Tillage Res 114:5056 10.1016/j.still.2011.04.001CrossRefGoogle Scholar
Singh, S, Singh, UB, Trivedi, M, Sahu, PK, Paul, S, Paul, D, Saxena, AK (2020) Seed biopriming with salt tolerant endophytic Pseudomonas geniculata-modulated biochemical responses provide ecological fitness in maize (Zea mays L.) grown in saline sodic soil. Int J Environ Res Public Health 17:253 10.3390/ijerph17010253CrossRefGoogle Scholar
Smith, RG, Gross, KL (2007) Assembly of weed communities along a crop diversity gradient. J Appl Ecol 44:10461056 10.1111/j.1365-2664.2007.01335.xCrossRefGoogle Scholar
Song, P, Feng, G, Brooks, J, Zhou, B, Zhou, H, Zhao, Z, Li, Y (2019) Environmental risk of chlorine-controlled clogging in drip irrigation system using reclaimed water: the perspective of soil health. J Cleaner Prod 232:14521464 10.1016/j.jclepro.2019.06.050CrossRefGoogle Scholar
Sosnoskie, LM, Herms, CP, Cardina, J (2006) Weed seedbank community composition in a 35-yr-old tillage and rotation experiment Weed Sci 54:263–7310.1614/WS-05-001R2.1CrossRefGoogle Scholar
Srimalanon, P, Prapagdee, B, Sombatsompop, N (2020) Soil inoculation with Pseudomonas geniculata WS3 for accelerating the biodegradation process of in situ compatibilized PBS/PLA blends doped with HPQM. J Polym Environ 28:11381149 10.1007/s10924-020-01670-6CrossRefGoogle Scholar
Stika, JA (2013) Addressing some hurdles to implementing soil health management systems. J Soil Water Conserv 68:99A 10.2489/jswc.68.4.99ACrossRefGoogle Scholar
Subramanian, P, Mageswari, A, Kim, K, Lee, Y, Sa, T (2015) Psychrotolerant endophytic Pseudomonas sp. strains OB155 and OS261 induced chilling resistance in tomato plants (Solanum lycopersicum Mill.) by activation of their antioxidant capacity. Mol Plant-Microbe Interact 28:10731081 10.1094/MPMI-01-15-0021-RCrossRefGoogle ScholarPubMed
Teasdale, JR (1993) Interaction of light, soil moisture, and temperature with weed suppression by hairy vetch residue. Weed Sci 41:4651 10.1017/S0043174500057568CrossRefGoogle Scholar
Teasdale, JR, Brandsaeter, LO, Calegari, A, Neto, FS (2007) Cover crops and weed management. Pages 49–64 in Upadhyaya MK, Blackshaw RE, eds. Non-chemical Weed Management: Principles, Concepts and Technology. Wallingford, UK: CABI10.1079/9781845932909.0049CrossRefGoogle Scholar
Teasdale, JR, Mohler, CL (2000) The quantitative relationship between weed emergence and the physical properties of mulches. Weed Sci 48:385392 10.1614/0043-1745(2000)048[0385:TQRBWE]2.0.CO;2CrossRefGoogle Scholar
Thangarajan, R, Bolan, NS, Tian, G, Naidu, R, Kunhikrishnan, A (2013) Role of organic amendment application on greenhouse gas emission from soil. Sci Total Environ 465:7296 10.1016/j.scitotenv.2013.01.031CrossRefGoogle ScholarPubMed
Tilman, D (1985) The resource-ratio hypothesis of plant succession. Am Nat 125:827–5210.1086/284382CrossRefGoogle Scholar
Trinka, DL, Pritts, MP (1992) Micropropagated raspberry plant establishment responds to weed control practice, row cover use, and fertilizer placement. J Am Soc Hort Sci 117:874880 10.21273/JASHS.117.6.874CrossRefGoogle Scholar
Uzun, Y, Yakar, S, Karacan, I, Halil, Kaya A (2019) New additions to the Turkish Pezizales. Turk J Bot 42:335345 Google Scholar
van Heemst, HDJ (1985) The influence of weed competition on crop yield. Agric Syst 18:8193 10.1016/0308-521X(85)90047-2CrossRefGoogle Scholar
Wang, R, Zhang, H, Sun, L, Qi, G, Chen, S, Zhao, X (2017) Microbial community composition is related to soil biological and chemical properties and bacterial wilt outbreak. Sci Rep 7:343 10.1038/s41598-017-00472-6CrossRefGoogle ScholarPubMed
Weil, RR, Islam, KR, Stine, MA, Gruver, JB, Samson-Liebig, SE (2003) Estimating active carbon for soil quality assessment: a simplified method for laboratory and field use. Am J Altern Agric 18:317 Google Scholar
Wells, MS, Reberg-Horton, SC, Smith, AN, Grossman, JM (2013) The reduction of plant-available nitrogen by cover crop mulches and subsequent effects on soybean performance and weed interference. Agron J 105:539545 10.2134/agronj2012.0396CrossRefGoogle Scholar
Wilhelm, RC, van Es, HM, Buckley, DH (2022) Predicting measures of soil health using the microbiome and supervised machine learning. Soil Biol Biochem 164:108472 10.1016/j.soilbio.2021.108472CrossRefGoogle Scholar
Williams, A, Wells, MS, Dickey, DA, Hu, S, Maul, J, Raskin, DT, Reberg-Horton, SC, Mirsky, SB (2018) Establishing the relationship of soil nitrogen immobilization to cereal rye residues in a mulched system. Plant Soil 426:95107 10.1007/s11104-018-3566-0CrossRefGoogle Scholar
Wortman, SE, Davis, AS, Schutte, BJ, Lindquist, JL (2011) Integrating management of soil nitrogen and weeds. Weed Sci 59:162170 10.1614/WS-D-10-00089.1CrossRefGoogle Scholar
Wu, C, Wu, X, Chen, S, Wu, D (2020) A newly discovered humic-reducing bacterium, Pseudomonas geniculata PQ01, isolated from paddy soil promotes paraquat anaerobic transformation. Front Microbiol 11, https://doi.org/10.3389/fmicb.2020.02003 CrossRefGoogle Scholar
Zak, DR, Tilman, D, Parmenter, RR, Rice, CW, Fisher, FM, Vose, J, Milchunas, D, Martin, CW (1994) Plant production and soil microorganisms in late-successional ecosystems: a continental-scale study. Ecology 75:23332347 10.2307/1940888CrossRefGoogle Scholar
Figure 0

Table 1. Carbon (C) and nitrogen (N) content of amendments added to the soil in 2020 and 2021 (Homer C. Thompson Vegetable Research Farm in Freeville, NY, USA)

Figure 1

Table 2. P-values, degrees of freedom (df), and F-values of crop and weed aboveground biomass from a type II ANOVA of a linear model for each independent variable included in the model

Figure 2

Figure 1. Soybean crop yield and total weed biomass in soils with different amendments. Bars labeled with different letters indicate significantly different means (P < 0.05) by treatment.

Figure 3

Table 3. Significant indicator species identified in plots with different amendment treatments

Figure 4

Figure 2. Average weed species traits in plots treated with different amendments for corn and soybean plots. Trait values were retrieved from the literature for each species and weighted by the number of that species found in a plot. Box plots labeled with different letters indicate significantly different means (P < 0.05) by treatment. Average weed species traits not significantly affected by soil amendments are not labeled with letters.

Figure 5

Figure 3. Percent soybean and corn crop yield relative to the total aboveground plant biomass within a square meter in soils with different amendments. Bars labeled with different letters indicate significantly different means (P < 0.05) by treatment. Years where percent plant biomass was not significantly affected by soil amendments are not labeled with letters.

Figure 6

Figure 4. Corn crop yield and total weed biomass in soils with different amendments. Bars labeled with different letters indicate significantly different means (P < 0.05) by treatment. Years where total plant biomass was not significantly affected by soil amendments are not labeled with letters.

Figure 7

Table 4. Mean soil measurements (respiration, nitrate, and ammonium) in soils with different amendment treatmentsa

Figure 8

Figure 5. Mean total soil respiration as affected by amendment treatment in 2020 and 2021. Measurements were taken twice in 2020 (July and August) and four times in 2021 (June, July, August, and September). Bars labeled with different letters indicate significantly different means (P < 0.05) by treatment within year.

Figure 9

Figure 6. Mean total nitrate as affected by amendment treatment in 2020 and 2021. Measurements were taken twice in 2020 (July and August) and four times in 2021 (June, July, August, and September). Bars labeled with different letters indicate significantly different means (P < 0.05) by treatment within year.

Figure 10

Figure 7. Mean total ammonium as affected by amendment treatment in 2020 and 2021. Measurements were taken twice in 2020 (July and August) and four times in 2021 (June, July, August, and September). Bars labeled with different letters indicate significantly different means (P < 0.05) by treatment within year.

Figure 11

Figure 8. (A) Average soil health ratings for each plot treated with different amendments. Box plots labeled with different letters indicate significantly different means (P < 0.05) by treatment. (B) Principal component analysis (PCA) of centered and scaled soil health indicator values. Shape indicates crop, size indicates soil health rating, and color indicates amendment treatment. All indicators were significantly affected by treatment (P < 0.05), except indicators outlined in red. ACE, autoclaved-citrate extractable soil protein index.

Figure 12

Table 5. The P-value and mean soil health indicator value by soil amendment treatmenta

Figure 13

Figure 9. Principal coordinates analysis (PCoA) of the bacterial (A) and fungal (B) communities based on Illumina MiSeq amplicon sequencing. Shape indicates crop, color indicates amendment treatment, and ellipse line type indicates year. Permutational multivariate analysis of variance (PERMANOVA) revealed significant effects of soil amendment (pseudo P-value = 0.001 for bacteria and fungi), crop (pseudo P-value = 0.005 for bacteria; pseudo P-value = 0.001 for fungi), and year (pseudo P-value = 0.001 for bacteria and fungi) on microbial community composition.

Figure 14

Table 6. The P-values and degrees of freedom (df) for bacterial (16S rRNA region) and fungal (ITS2 region) beta diversity, measured with a permutational multivariate analysis of variance (PERMANOVA) on a principal coordinates analysis (PCoA) of a Bray-Curtis distance matrix of operational taxonomic units (OTUs)

Figure 15

Figure 10. The 10 bacterial (A) and fungal (B) operational taxonomic units (OTUs) whose abundances differed most significantly between rye-amended and no-rye plots. These OTUs were identified with multivariable associations with linear models (Maaslin2). Maaslin2 was run on nonrarefied abundance data, but percentage relative abundance data are shown. Identified taxa are ordered from most (top) to least (bottom) significantly different. Full taxonomic classification for each identified taxon can be found in Table 7.

Figure 16

Table 7. Top 10 bacterial and fungal operational taxonomic units (OTUs) with significantly different abundance between plots amended with rye hay and no-rye plots, as identified by the Microbiome Multivariable Association with Linear Models 2.0 algorithm

Supplementary material: File

Gannett et al. supplementary material

Gannett et al. supplementary material
Download Gannett et al. supplementary material(File)
File 129.8 KB