Next Article in Journal
New Antifungal Microbial Pigment Applied to Improve Safety and Quality of Processed Meat-Products
Next Article in Special Issue
Preservation of Underground Microbial Diversity in Ancient Subsurface Deposits (>6 Ma) of the Rio Tinto Basement
Previous Article in Journal
The Effect of Antimicrobial Treatment upon the Gill Bacteriome of Atlantic Salmon (Salmo salar L.) and Progression of Amoebic Gill Disease (AGD) In Vivo
Previous Article in Special Issue
Influence of Temperature and Sulfate Concentration on the Sulfate/Sulfite Reduction Prokaryotic Communities in the Tibetan Hot Springs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Straw Incorporation with Nitrogen Amendment Shapes Bacterial Community Structure in an Iron-Rich Paddy Soil by Altering Nitrogen Reserves

1
College of Environmental Science and Engineering, Yangzhou University, Yangzhou 225127, China
2
Zhenjiang Agricultural Technology Extension Station, Zhenjiang 212009, China
*
Author to whom correspondence should be addressed.
Submission received: 8 March 2021 / Revised: 21 April 2021 / Accepted: 26 April 2021 / Published: 3 May 2021

Abstract

:
Incorporation of crop straw into the soil along with inorganic fertilization is a widespread agricultural practice and is essential in nutrient-scarce soils, such as iron-rich (ferruginous) paddy soils. The responses of soil bacterial communities to straw incorporation under different nitrogen inputs in iron-rich soils remain unclear. Therefore, 6000 kg ha−1 dry wheat (Triticum aestivum L. cv. Zhengmai 12) straw was applied to a rice paddy with and without nitrogen amendment (0, 80, 300, and 450 kg ha−1 N as urea), to investigate its effects on soil fertility and bacterial community structure. Organic matter, total nitrogen, and water contents tended to decrease in straw-incorporated soils with different nitrogen inputs. Proteobacteria was the dominant bacterial phylum across all treatments (26.3–32.5% of total sequences), followed by Chloroflexi, Acidobacteria, and Nitrospirae. Up to 18.0% of all the taxa in the bacterial communities were associated with iron cycling. Straw incorporation with nitrogen amendment increased the relative abundance of iron oxidizers, Gallionellaceae, while decreasing the relative abundance of iron reducers, Geobacteraceae. Bacterial community composition shifted in different treatments, with total nitrogen, water, and Fe(III) contents being the key drivers. Straw incorporation supplemented by 300 kg ha−1 N increased bacterial richness and enhanced all the predicted bacterial functions, so that it is recommended as the optimal nitrogen dosage in practice.

1. Introduction

Straw application is a widespread agronomic practice in developed countries [1,2]. However, in China, burning and discarding still account for a large fraction of the disposed straw under an annual straw productivity of 900 million tons [3]. Such improper agronomic practices have not only led to nutrient losses [4] but also to environmental challenges. Open burning of crop straw, for instance, has raised major concerns with regard to air pollution [5]. In recent years, straw application has become increasingly popular in China [6], and numerous studies have investigated the effects of straw application on soil fertility [7,8,9]. Straw application could improve soil structure and water content, as well as soil nitrogen (N) content and N use efficiency [10]. Previous studies have also reported that long-term straw application enhances carbon (C) sequestration considerably, by 50–100%, in China’s subtropical paddy soils [11,12]. According to Potthoff et al. [13], applying straw could increase soil microbial C, and, in turn, soil C stocks, which are already relatively high in Chinese croplands [14,15].
The effects of straw application on soil nutrients are rather complex. Generally, soil organic matter (SOM) is influenced by soil type and biomass, among other factors. Different types of straw have variable C/N ratios, which show equally varying effects on SOM [16]. In addition, SOM content affects soil physicochemical and biological properties, which could control straw decomposition and nutrient cycling [17]. Due to the high C/N ratios of cereal crop straw (up to 80% in wheat straw), the N in straw is mostly retained by microbes during decomposition [18,19]. Studies have shown that straw application without additional N fertilization decreases rice yield, potentially due to competition between straw decomposition and rice growth, which impairs N uptake by plants [20]. Therefore, straw application to the field is usually supplemented by N fertilizer, and N input rate could influence crop productivity. However, few studies have explored the effects of straw application with different N fertilization rates on soil nutrient conditions.
Soil microorganisms balance soil C and other nutrients, and have been demonstrated to facilitate the maintenance of soil health and crop growth [21,22,23]. Soil microbial communities exhibit diverse responses to straw application. According to Tardy et al. [24], soil bacteria exhibit strong responses to wheat straw incorporation, with both species richness and evenness decreasing transiently, regardless of soil history. Changes in soil nutrient status following straw application could be directly associated with microbial processes and could regulate microbial community composition and function [25,26,27]. However, inconsistent results have been reported with regard to the effects of crop straw application on soil bacterial community structure. Specifically, the effects of straw application could vary with different soil types, fertilizer types, straw types, climatic conditions, and agricultural practices [25,28,29]. Consequently, further studies on the effects of straw application on soil bacterial communities are required.
Iron (Fe) is a redox-sensitive metal that is critical for the cycling of soil organic C and other elements. Soil Fe oxides strongly interact with SOM and facilitate SOM stabilization [30]. Periodic wet and dry conditions in paddy soils lead to continual redox reactions of Fe [31], which could be accelerated by diverse microorganisms. For example, Fe(III) reduction is responsible for a large proportion of C mineralization in soils (up to 80%) [32], which is correlated with the presence of Fe(III)-reducing microorganisms. Furthermore, researchers have demonstrated that Fe(III) reducers might compete for electrons with methanogens, leading to decreased methane (CH4) emissions in rice paddies [33,34]. In addition, Hu et al. [35] showed that Fe(II) addition decreased CH4 emissions induced by rice straw application in flooded paddy soils.
Soil Fe oxides adsorb a wide range of organic anions and cations [36], and the formation of organo-Fe complexes could reduce the microbial degradation of labile organic molecules, such as glucose and citrate [23]. Conversely, organic compounds regulate Fe transformation in soil environments. Adding plant residues could stimulate microbial reduction of native soil Fe oxides and, in turn, influence the associated microbial community structures [37]. Moreover, Fe(II) oxidation and Fe(III) reduction are coupled to N biogeochemical processes in flooded paddy soils [38]. Paddy soils are key habitats for microorganisms involved in Feammox, a major N loss pathway, with potential major bacterial drivers such as Geobacter, Pseudomonas, and Thiobacillus [39,40,41]. Reduced forms of Fe(II) would be oxidized by microaerophilic or nitrate-reducing Fe oxidizers that are ubiquitous in paddy soils, which accelerates the turnover of N [42]. However, it remains unclear how soil bacterial communities respond to straw application supplemented by N fertilizer in paddies with high Fe concentrations.
The Yangtze River Delta in China is a major rice production area. Numerous studies have previously been carried out to explore the effects of straw application on soil fertility and crop yield in this region [43]. However, our understanding of the response of soil bacterial communities to straw application and the optimal N fertilization rate when supplemented with straw is still poor [44], specifically in Fe-rich (ferruginous) paddy soils. Therefore, the objectives of the present study were to investigate the effects of straw incorporation on soil bacterial community structure and to determine the optimal N application rate in combination with straw incorporation for improving soil function in an Fe-rich rice paddy.

2. Materials and Methods

2.1. Study Site

The field experimental site (31°53′36.70″ N, 119°24′2.81″ E) was located in Zhenjiang, Jiangsu Province, China. This site belongs to the northern subtropical climate zone, with an average annual temperature of 15.6 °C. The average annual precipitation is 1088.3 mm and occurs mostly from June to August. Paddy soil is the main soil type at the experimental site. Due to the influence of irrigation and seasonal precipitation, there are frequent shifts in soil oxidation status, which leads to the high accumulation of Fe in the tillage layer. Fe plaque is also common in rice paddies at the experimental site.

2.2. Experimental Design and Soil Sampling

A rice paddy with an area of 0.1 ha was selected for straw incorporation (Figure 1). A randomized complete block design was adopted for the experiments. Five fertilization treatments with different N rates were applied under wheat straw incorporation (ST0, ST1, ST2, and ST3) or no straw incorporation (T1; Table 1). The wheat straw applied was from the previous season. It contained 460 g kg−1 total C, 3.17 g kg−1 total N, 0.38 g kg−1 total phosphorus (TP), 14.64 g kg−1 total potassium (TK), 1.4 g kg−1 total sodium, 1.80 g kg−1 total calcium, 1.20 g kg−1 total magnesium, 1.94 g kg−1 total sulfur (S), and 3.20 g kg−1 total chlorine. The straw was chopped into small pieces (<2 mm in length) and mixed thoroughly with the surface soil via ploughing. Additionally, N, P, and K were applied as basal fertilizers together with straw (Table 1). Each treatment had three triplicates, yielding 15 experimental plots (30 m2). The rice variety under cultivation was Oryza sativa subsp. Keng cv. Nanjing 9108. Rice seedlings were transplanted on June 10th and harvested on October 30th. The experiment was repeated each year between 2017 and 2019. The results presented in this paper were observed in the third year (2019).
Soil samples were obtained at the rice ripening stage, right before draining the field. A minimum of 10 soil cores (5 cm in diameter and 15 cm in depth) from each plot were obtained randomly. The topsoil (0–2 cm) was removed to minimize the potential influence of debris. The samples obtained from one plot were pooled, mixed uniformly, and then divided into two portions. One portion was transported on ice to the laboratory and frozen immediately at −80 °C for use in subsequent DNA extraction and sequencing analyses, and the other portion was air-dried for the determination of soil physicochemical properties. The above procedures were repeated for each plot and a total of 15 samples were taken.

2.3. Soil Physicochemical Analysis

Soil samples were air-dried and sieved with 0.85 mm and 0.15 mm sieves for use in the determination of pH and other properties, respectively. Deionized water and soil were mixed at a ratio of 2.5:1 (v/w) and vortexed; soil pH and electrical conductivity (EC) were determined using a digital pH meter (Sanxin S731; Shanghai San-Xin Instruments, Shanghai, China). Gravimetric water content (WC) was measured as the difference in WC between the moist soil and the soil dried at 105 °C for 24 h. SOC, TN, and TP contents were determined according to standard testing procedures [45]. Briefly, SOC was determined using the potassium dichromate volumetric-external heating method; TN was determined using the Kjeldahl method; TP was determined using the HClO4-H2SO4 extraction–Mo-Sb anti-spectrophotometric method. Fe(II) and Fe(III) contents were determined using the 0.5 N HCl extraction and phenathroline spectrophotometric method [46]. Non-crystal Fe oxide content was determined using the dithionite-citrate-bicarbonate (DCB) extraction technique and phenathroline method [47].

2.4. DNA Extraction, PCR Amplification, and Illumina MiSeq Sequencing

Genomic DNA was extracted from 0.5 g freeze-dried soil samples using the MoBio PowerSoil DNA Isolation Kit (QIAGEN Inc., Valencia, CA, USA) according to the manufacturer’s protocols. The quality and quantity of the DNA samples were checked using a spectrophotometer (Nano Drop ND2000; Thermo Scientific, Wilmington, DE, USA).
To amplify the hypervariable V3–V4 regions of the 16S rRNA gene, the bacterial primer set 338F (5′-ACTCCTACGGGAGGCAGCAG-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′) was used [48]. A PCR reaction contained a 25 μL mixture with 10 μL MilliQ water, 5 μL 5 × FastPfu Buffer, 2 μL 2.5 mM dNTPs, 1.0 μL of primer 338F (5 μM), 1.0 μL of primer 806R (5 μM), 0.5 μL FastPfu Polymerase, 10 ng template DNA, and 0.25 μL bovine serum albumin. The amplification conditions were as follows: DNA denaturation at 94 °C for 3 min, followed by 30 cycles of amplification (94 °C for 45 s, 50 °C for 45 s, and 72 °C for 45 s), and final extension for 10 min at 72° C. PCR reactions were performed in triplicate. The PCR products were gel-purified using an AxyPrepDNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, USA) and quantified using QuantiFluor™-ST (Promega, Madison, WI, USA) according to the manufacturers’ instructions. Sequencing of the purified PCR products was conducted on an Illumina MiSeq platform (Illumina, San Diego, CA, USA) by Majorbio Bio-Pharm Technology Co. Ltd. (Shanghai, China) according to standard protocols, with 250-bp paired-end reads generated.
The sequences obtained were analyzed using Mothur v1.30.2 [49]. Raw sequences were assigned to samples according to the barcodes, and those that did not perfectly match the primer/barcode or were <250 bp in length were removed. Chimeras were removed using the UCHIME algorithm with default parameters implemented in Mothur. The unique sequences were merged, chopped to achieve similar lengths, and then aligned against the SILVA132 16S rRNA database (https://www.arb-silva.de; accessed on 2 September 2020) on the Mothur platform to obtain taxonomic information on bacterial communities. Non-bacterial reads were further removed. Operational taxonomic units (OTUs) were clustered based on 97% similarity [50]. The number of sequences in each sample was normalized to the minimum number of sequences per sample. Bacterial alpha-diversity was estimated using the observed species (Sobs), ACE, and Shannon indexes generated based on the OTU counts.

2.5. Data Analysis

Bioinformatics analysis of bacterial sequence data was performed using a free online platform, Majorbio I-Sanger Cloud (http://www.i-sanger.com; accessed on 5 February 2021). One-way analysis of variance (ANOVA) followed by the Tukey–Kramer test was performed to investigate differences in soil physiochemical properties, bacterial alpha-diversity indexes, and dominant taxa abundances among treatments, with p < 0.05 considered to indicate significant difference. Non-metric multidimensional scaling (NMDS) analysis based on the Bray–Curtis distances was performed on bacterial community composition at both the phylum and genus levels, with the significance of differences tested by the analysis of similarities (ANOSIM). The influences of environmental factors on bacterial community structure were estimated using redundancy analysis (RDA). The correlations between abundant bacterial taxa (the top 15 phyla and 30 genera) and soil properties were calculated based on Spearman’s correlation. Functional genes and metabolic pathways were predicted with Phylogenetic Investigation of Communities by Re-construction of Unobserved States (PICRUSt v1.1.0; https://github.com/picrust; accessed on 8 February 2021).

3. Results

3.1. Soil Physicochemical Properties

The soil samples were acidic, with an average pH of 5.05 (Table 2). There was no statistically significant difference in soil pH among the treatments. However, the ST1 treatment with straw and low N had the highest EC, which was significantly higher than that of the T1 treatment with low N alone (p < 0.05). In contrast, the highest SOM (32.29 g kg−1), TP (0.60 g kg−1), TN (1.75 g kg−1), and WC (65.94%) contents were observed in the T1 treatment. Compared with T1, straw incorporation without N (ST0) and with low N (ST1) decreased SOM; surprisingly, straw incorporation with N (ST1 to ST3) decreased TN and WC (p < 0.05).
The rice paddy was rich in Fe and rust-colored Fe oxides were commonly observed in soil cores. The dynamics of soil Fe content were rather complex (Table 1). The bioavailable Fe and non-crystal Fe oxide contents were not significantly influenced by the different treatments. However, Fe(II) content was markedly low in the ST2 treatment, and Fe(III) content was relatively low in the ST0, ST1, and ST2 treatments, whereas the lowest DCB-Fe content was observed in the T1 treatment.

3.2. Soil Bacterial Community Diversity

The calculated alpha-diversity indexes of bacterial communities in soil samples associated with different treatments are listed in Table 3. The Sobs and ACE indexes represent bacterial species richness, while the Shannon index represents bacterial species evenness; higher values of these indexes suggest greater bacterial richness and lower heterogeneity, respectively. Overall, the coverage between the samples was similar, ranging from 95.9% to 97.2%, indicating that the sequences were representative. The ST2 treatment had the highest Sobs and ACE index values. The highest Shannon index values were observed in the T1 treatment, and they were significantly higher than those in ST3 treatments (p < 0.05). Generally, straw incorporation alone or with low N decreased bacterial richness slightly; straw incorporation with moderate N increased bacterial richness but decreased evenness, whereas straw incorporation with high N decreased both bacterial richness and evenness.
NMDS analysis was performed to visualize variations in bacterial community composition among the treatments (Figure S1). A shift in bacterial community composition was observed following straw incorporation and the bacterial community in T1 samples was clearly separated from the remaining treatment groups at both the phylum and genus levels. An increase in N application also induced a shift in bacterial community composition under straw incorporation, with ST3 forming an individual cluster and clearly separating from the communities in ST0, ST1, and ST2 samples, which were not entirely separate. Furthermore, ANOSIM test results revealed that the bacterial community compositions in T1 and ST3 were significantly different (p < 0.05) from those of the remaining treatment groups.

3.3. Taxonomic Composition of Bacterial Communities

Based on the taxonomic classification of 16S rRNA sequences, 54 bacterial phyla were recognized across all samples. Figure 2 presents the taxonomic assignments of the sequences at the phylum or family level. Proteobacteria was the most dominant phylum in all samples and accounted for 28.5% of the total sequences (Figure 2A). The second most abundant phylum was Chloroflexi (25.1%), followed by Acidobacteria (16.0%), Nitrospirae (6.72%), Actinobacteria (3.98%), Bacteroidetes (3.41%), and Gemmatimonadetes (3.32%). These seven dominant phyla accounted for 87.0% of the total bacterial community.
Proteobacteria had the highest relative abundance in the ST3 treatment (32.5%) of the total sequences, and its lowest relative abundance was observed in the ST0 treatment (26.3%). Within the phylum Proteobacteria, the majority (72%) of the sequences belonged to the classes Betaproteobacteria and Detaproteobacteria. Alphaproteobacteria and Gammaproteobacteria accounted for 17.9% and 8.48% of the Proteobacterial sequences, respectively. The distribution patterns of Chloroflexi were opposite to those of Proteobacteria, with the highest relative abundance observed in the ST0 treatment (29.3%) and the lowest in the ST3 treatment (19.8%). The relative abundances of the dominant phyla varied somewhat among the treatments; however, significant differences were observed only in the less dominant phyla, including Verrucomicrobia, Ignavibacteriae, Planctomycetes, AC1, and Peregrinibacteria (p < 0.05; Figure S2A). The highest relative abundances of Verrucomicrobia and Planctomycetes were observed in the ST1 treatment, while the relative abundance of Ignavibacteriae was the highest in the ST3 treatment.
Bacterial families with relative abundances > 3% included Anaerolineaceae, norank_c_SBR2076, norank_c_Nitrospira, norank_c_Acidobacteria, norank_c_KD4-96, Solibacteraceae_Subgroup_3, Gallionellaceae, and Acdiobacteriaceae_Subgroup_1 (Figure 2B). The top two dominant families belonged to the phylum Chloroflexi. Among the dominant families, the relative abundance of Anaerolineaceae was the highest in the T1 treatment (7.7%). Norank_c_Nitrospira and Gallionellaceae had their highest relative abundances in the ST3 treatment (9.02% and 5.56%, respectively). In addition, the relative abundances of unclassified norank_c_SBR2076 and norank_c_KD4-96, belonging to the phylum Chloroflexi, were significantly lower in ST3 than in ST0 and ST1 treatments, respectively (p < 0.05).
More than 18.0% of the total sequences belonged to Fe-cycling-related bacteria. The relative abundance of Gallionellaceae, a typical Fe(II)-oxidizing bacterial group, was significantly higher in the ST3 treatment than in the ST0 treatment (p < 0.05). The relative abundance of Geobacteraceae, a family containing many Fe(III)-reducing bacteria, was significantly higher in the T1 treatment than in the ST1 treatment (p < 0.05; Figure S2B).

3.4. Influence of Environmental Factors on Bacterial Community Structure

The relationships between the distribution of the major bacterial groups in rice paddy soils and the soil physicochemical properties in different treatments were analyzed by RDA (Figure 3). The first two axes explained 78.13% of the total variance in the bacterial community (RDA1 = 50.44% and RDA2 = 28.19%) at the phylum level. The bacterial community structure differed among the treatments, with T1 and ST3 samples clearly separated from the remaining treatment groups, suggesting that there were distinct geochemical processes influencing bacterial community structure in the two treatments. Bacterial community compositions in ST0, ST1, and ST2 samples were relatively closely related. Soil TN content (r2 = 0.53, p = 0.015), water content (r2 = 0.44, p = 0.03), and Fe(III) content (r2 = 0.66, p = 0.004) were the key environmental factors influencing the bacterial community composition. Soil bioavailable Fe(III) and SOC contents shaped bacterial community structure in the ST3 treatment, while soil pH and DCB-Fe content were more associated with the bacterial community composition in the ST2 treatment.
To further reveal the influence of environmental factors on soil bacterial community structure, the relationships between individual bacterial groups and soil properties were analyzed using Spearman correlation coefficients (Table 4). The relative abundances of Proteobacteria, Chloroflexi, Acidobacteria, Actinobacteria, Gemmatimonadetes, Ignavibacteriae, unclassified_k_norank_d_Bacteria, and Spirochaetae were negatively correlated with soil TN content (p < 0.05). The relative abundances of Chloroflexi, Actinobacteria, Gemmatimonadetes, and Spirochaetae were negatively correlated with SOM content (p < 0.05). In addition, the relative abundance of unclassified_k_norank_d_Bacteria was negatively correlated with soil Fe(II) content (p < 0.05), and the relative abundance of Spirochaetae was negatively correlated with soil DCB-Fe content (p < 0.05). The relative abundances of Chloroflexi and Actinobacteria were negatively correlated with soil Fe(III) content, whereas the relative abundance of Verrucomicrobia was positively correlated with this soil property (p < 0.05).

3.5. Predicted Functions of Bacterial Communities

PICRUSt analysis was performed to predict the potential functions of soil bacterial communities in the different treatments (Figure 4). The abundances of major functional pathways, including metabolism, cellular processes, environmental information processing, and genetic information processing, were compared among the treatments. Notably, the ST2 treatment had the highest abundances for all functional pathways, followed by ST1 and ST0 treatments. The T1 treatment had the lowest abundance of all functional pathways (Figure 4A). Further analyses on the energy metabolism pathways revealed that the predicted pathways associated with C fixation, methane production, and oxidative phosphorylation were relatively more abundant (Figure 4B). In addition, the ST2 treatment had the highest abundances of pathways associated with S, N, and CH4 metabolism, C fixation, photosynthesis, and oxidative phosphorylation. Overall, straw incorporation with N amendment was beneficial for cell growth and metabolism, and the positive effect first increased with increasing N input, and then dropped when a high N rate (450 kg ha−1) was applied.

4. Discussion

4.1. Effects of Straw Incorporation with N Amendment on Paddy Soil Fertility

Straw application has a compound effect on soil nutrient status. On one hand, it could increase SOM [20,25] and stimulate soil microbial growth and activities [44], which, in turn, accelerate the release of diverse nutrients. Conversely, excess C input could cause N immobilization through microbial activities, and, in turn, competition for N between plants and microbes [51]. However, according to our results, straw incorporation did not increase SOM content in the rice paddy, which is inconsistent with the findings of many previous studies [11,12,13,14,15]. One of the reasons could be that the ameliorating effects of straw incorporation occurred in the short term and microbial and plant growth took up much more organic matter in the long term. This might even be associated with SOM decomposition stimulated by straw input, which could explain the highest SOM content in the T1 treatment that had no straw incorporation. Similarly, strong regional variation in the effects of straw application on SOC stocks, and a lack of increase in SOC stocks in regions with low SOC densities have been reported [52]. The decrease in soil TN content with straw incorporation could be attributed to straw-derived increase in N use efficiency and an increase in N loss under high N input [53,54].
Compared with the T1, soil WC dropped following straw incorporation, which could be due to a change in soil structure because of straw decomposition and enhanced bacterial activity. In a previous study, straw application accompanied with moderate N (225 kg ha−1) achieved the maximum grain yield of winter wheat with the improvement of soil structure [55]. Here, we observed that N input eventually decreased soil TN content, and the optimal N application rate, which resulted in the highest rice yield (Supplementary Table S2), could be ~300 kg ha−1. Rice yield decreased in the ST3 treatment, which had the highest input, suggesting that the optimal N application rate in the ST3 treatment should not be exceeded.

4.2. Effects of Straw Incorporation with N Amendment on Bacterial Community Structure

Soil bacterial communities regulate the biogeochemical processes of soil elements, while crop straw provides energy and nutrients required for the growth of soil bacteria [23]. Previous research has shown that straw application influences soil bacterial community composition profoundly [25]. Straw application not only provides rich nutrients but also improves soil structure, thereby increasing soil bacterial diversity [56]. In the present study, however, straw incorporation alone or with low N (80 kg ha−1) decreased soil bacterial richness slightly; in contrast, straw incorporation amended with moderate N (300 kg ha−1) increased bacterial richness, which mildly decreased again when N amendment was increased to a high rate (450 kg ha−1). High N could stimulate the growth of less abundant copiotrophic bacterial groups that are K strategists, and thus decrease the relative abundances of other dominant bacterial groups [57,58] despite an increase in C input [59,60]. The lowest bacterial diversity index values were observed in the ST3 treatment, which had the highest N input. This is probably due to the facilitation of major bacterial groups following an increase in N input.
Across all treatments, Proteobacteria was the most dominant bacterial phylum, followed by Chloroflexi, Acidobacteria, and Nitrospirae, which is consistent with the findings of previous studies in paddy soils [61,62]. Proteobacteria is composed of numerous classes that are sensitive to copiotrophic conditions (K strategy), such as Beta- and Delta-proteobacteria [63]. In the present study, Beta- and Delta-proteobacteria were only slightly more abundant in the ST2 and ST3 treatment when compared with the other treatments, indicating that only with moderate to high N supplementation could straw incorporation increase the abundance of such classes of Proteobacteria. Many taxa of this phylum are known to promote plant growth via improving soil functions and regulating plant immune systems [64], and thus contribute to rice production.
Only a small proportion of the soil microbial populations could benefit from straw incorporation as a C source either directly or indirectly as predators of primary degraders [65,66]. Straw breakdown requires the activities of various microbial groups [67] and low soil pH, which would alter the nutrient status for bacterial communities. Nevertheless, soil pH exhibited minimal change following straw incorporation in the present study. High N input decreased the relative abundances of several potential Fe(III)-reducing groups, including norank_c_SBR2076 and norank_c_KD4-96 of the phylum Chloroflexi, in addition to norank_c_Acidobacteria. Meanwhile, high N input increased the relative abundance of an Fe(II)-oxidizing group, Gallionellaceae, and a potentially keystone group, nonrank_c_Nitrospira with functions in Fe uptake and nitrite oxidation [68]. This could be related to the enhanced plant growth that increased the radial oxygen (O2) loss favorable for Fe(II)-oxidizing bacteria.

4.3. Relationship between Bacterial Community Structure and Soil Environment

Microorganisms enhance soil functions and productivity through their participation in nutrient cycling and organic matter turnover [69]. Conversely, soil environmental factors influence the community composition of soil microorganisms. In the present study, we observed that soil bacterial community composition in the samples without straw was distinct from that of samples with straw incorporation. Extra nutrients introduced by straw incorporation could have altered bacterial community composition. Bu et al. [28] also showed that soil bacterial communities were separated into two groups between a straw-incorporated treatment and non-straw-incorporated treatment in a rice–rice–rapeseed rotation system. Among the treatments with straw incorporation, the bacterial community composition in ST3 samples was distinct from that of the others, demonstrating that high N input had strong effects on soil bacterial community composition [57,59].
In this study, the paddy soil had very high Fe contents and rust-colored Fe oxides were common in the field. The rapid Fe cycling associated with straw decomposition in the Fe-rich paddy soil could be driven by bacterial communities. According to our RDA results, Fe(III) content was a key environmental factor, similar to TN content, influencing soil bacterial community structure. A previous study also showed that Fe oxides could influence soil bacterial community structure in arsenic-contaminated soils [70]. Therefore, we specifically explored the bacterial groups associated with Fe cycling.
Not surprisingly, Gallionellaceae, widely recognized Fe(II)-oxidizing bacteria, were abundant in all samples. Gallionellaceae are microaerophilic bacteria that grow chemoautotrophically using Fe(II) as an energy source and carbon dioxide as a C source under low O2 conditions, which are particularly important traits at the oxic–anoxic interfaces of Fe-rich environments [71]. According to the results of Spearman correlation analysis, Gallionella abundance was positively correlated with the content of non-crystal Fe oxides (Table S1). In addition, straw incorporation with N amendment increased the relative abundance of Gallionellaceae while decreasing the relative abundance of Geobacteraceae; these variations were probably caused by the decrease in soil WC after treatment.
Considering bacterial functions, the ST2 treatment had the highest abundances of pathways associated with S, N, and CH4 metabolism, C fixation, photosynthesis, and oxidative phosphorylation. Although there were no clear trends with regard to the effects of straw incorporation on soil bacterial community structure under different N inputs, straw incorporation with N amendment enhanced cell growth and metabolism, with ST2 (straw incorporation + 300 kg ha−1 N) being the optimal treatment. Many studies have explored the effects of straw application on the responses to total soil microbial biomass and community structure, with less attention directed at soil microbial diversity and function [58]. Therefore, future studies should focus not only on the effects of straw application on phylogenetic diversity of bacteria but also functional diversity of bacteria, to enhance our understanding of soil ecosystem function following straw application.

5. Conclusions

The results of the present study echo the observation that straw incorporation has complex effects on soil nutrient status and bacterial community structure in Fe-rich paddies. Straw incorporation, with application, decreased SOM, TN, and WC contents, without affecting pH and Fe bioavailability. The relative abundances of dominant bacterial phyla, including Proteobacteria, Chloroflexi, Acidobacteria, Nitrospirae, and Actinobacteria, were all negatively correlated with soil TN, TP, Fe(II), and WC contents. The bacterial communities comprised a large proportion of Fe cycling-related bacteria and responded distinctively to different treatments. Soil TN, WC, and Fe(III) contents were the key environmental factors driving the variation in bacterial community structure. In general, straw incorporation with N amendment at a moderate rate (300 kg ha−1) most effectively improved rice yield, increased bacterial richness, and enriched the predicted functional pathways associated with cell growth and metabolism. Therefore, 300 kg ha−1 N is identified and recommended as the optimal N application rate for supplementing straw incorporation in rice paddies.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/microorganisms9050988/s1, Figure S1: Non-metric multidimensional scaling analysis (NMDS) based on the relative abundance of microbial phylum (A) and genus showing different microbial community structure, Figure S2: The composition of bacteria community in different treatments at phylum (A), family (B) and genus (C) levels, Table S1: Spearman correlation analysis between soil physiochemical parameters and the relative abundance of genus, Table S2: Rice yield components in different treatments.

Author Contributions

Conceptualization, J.W., and Y.M.; methodology, J.W. and G.W.; software, J.W., Y.M., and X.Q.; validation, J.W.; formal analysis, J.W.; investigation, J.W., Y.M., L.D., X.Q.; writing—original draft preparation, J.W.; visualization, J.W., Y.M., and G.W.; supervision, J.W.; project administration, J.W. and G.W.; funding acquisition, J.W. and X.Q. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the “National Key R&D Program of China” (2017YFD0200107), the “Natural Science Foundation of China” (41701093), and the “China Postdoctoral Science Foundation” (2017M611928).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are openly available in the NCBI Sequence Read Archive, with the accession number PRJNA726334.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cook, H.F.; Valdes, G.S.B.; Lee, H.C. Mulch effects on rainfall interception, soil physical characteristics and temperature under Zea mays L. Soil Tillage Res. 2006, 91, 227–235. [Google Scholar] [CrossRef]
  2. Lal, R. Beyond Copenhagen: Mitigating climate change and achieving food security through soil carbon sequestration. Glob. Food Sec. 2010, 2, 169–177. [Google Scholar] [CrossRef]
  3. Feng, J.; Wang, B.; Zhang, D.; Chu, S.H.; Zhi, Y.; Hayat, K.; Wang, J.C.; Chen, X.F.; Hui, N.; Zhou, P. Streptomyces griseorubens JSD-1 promotes rice straw composting efficiency in industrial-scale fermenter: Evaluation of change in physicochemical properties and microbial community. Bioresour. Technol. 2021, 321, 124465. [Google Scholar] [CrossRef] [PubMed]
  4. Medina, J.; Monreal, C.; Barea, J.M.; Arriagada, C.; Borie, F.; Cornejo, P. Crop residue stabilization and application to agricultural and degraded soils: A review. Waste Manag. 2015, 42, 41–54. [Google Scholar] [CrossRef]
  5. Chen, L.; Zhang, J.B.; Zhao, B.Z.; Yan, P.; Zhou, G.X.; Xin, X.L. Effects of straw amendment and moisture on microbial communities in Chinese fluvo-aquic soil. J. Soils Sediments 2014, 14, 1829–1840. [Google Scholar] [CrossRef]
  6. Li, H.; Dai, M.W.; Dai, S.L.; Dong, X.J. Current status and environment impact of direct straw return in China’s cropland-A review. Ecotoxicol. Environ. Saf. 2018, 159, 293–300. [Google Scholar] [CrossRef] [PubMed]
  7. Ding, L.L.; Cheng, H.; Liu, Z.F.; Ren, W.W. Experimental warming on the rice-wheat rotation agro-ecosystem. Plant Sci. J. 2013, 31, 49–56. [Google Scholar] [CrossRef]
  8. Zhu, L.Q.; Hu, N.J.; Yang, M.F.; Zhan, X.H.; Zhang, Z.W. Effects of different tillage and straw return on soil organic carbon in a rice-wheat rotation system. PLoS ONE 2014, 9, e0088900. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Yamashita, T.; Feiner, H.; Bettina, J.; Helfrich, M.; Ludwig, B. Organic matter in density fractions of water-stable aggregates in silty soils: Effect of land use. Soil Biol. Biochem. 2006, 38, 3222–3234. [Google Scholar] [CrossRef]
  10. Yang, H.S.; Yang, B.; Dai, Y.J.; Xu, M.M.; Koide, R.T.; Wang, X.H.; Liu, J.; Bian, X.M. Soil nitrogen retention is increased by ditch-buried straw return in a rice-wheat rotation system. Eur. J. Agron. 2015, 69, 52–58. [Google Scholar] [CrossRef]
  11. Wu, J. Carbon accumulation in paddy ecosystems in subtropical China: Evidence from landscape studies. Eur. J. Soil Sci. 2011, 62, 29–34. [Google Scholar] [CrossRef]
  12. Majumder, B.; Mandal, B.; Bandyopadhyay, P.; Gangopadhyay, A.; Mani, P.; Kundu, A.; Mazumdar, D. Organic amendments influence soil organic carbon pools and rice-wheat productivity. Soil Sci. Soc. Am. J. 2008, 72, 775–785. [Google Scholar] [CrossRef]
  13. Potthoff, M.; Dyckmans, J.; Flessa, H.; Muhs, A.; Beese, F.; Joergensen, R.G. Dynamics of maize (Zea mays L.) leaf straw mineralization as affected by the presence of soil and the availability of nitrogen. Soil Biol. Biochem. 2005, 37, 1259–1266. [Google Scholar] [CrossRef]
  14. Wang, J.; Wang, X.; Xu, M.; Feng, G.; Zhang, W.; Lu, C. Crop yield and soil organic matter after long-term straw return to soil in China. Nutr. Cycl. Agroecosyst. 2015, 102, 371–381. [Google Scholar] [CrossRef]
  15. Zhao, Y.C.; Wang, M.Y.; Hu, S.J.; Zhang, X.D.; Ouyang, Z.; Zhang, G.L.; Huang, B.; Zhao, S.W.; Wu, J.S.; Xie, D.T.; et al. Economics-and policy-driven organic carbon input enhancement dominates soil organic carbon accumulation in Chinese croplands. Proc. Natl. Acad. Sci. USA 2018, 115, 4045–4050. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Zhao, X.M.; He, L.; Zhang, Z.D.; Wang, H.B.; Zhao, L.P. Simulation of accumulation and mineralization (CO2 release) of organic carbon in chernozem under different straw return ways after corn harvesting. Soil Tillage Res. 2016, 156, 148–154. [Google Scholar] [CrossRef]
  17. Ryan, M.; Terry, J.; Thomas, W. Bacterial metataxonomic profile and putative functional behavior associated with C and N cycle processes remain altered for decades after forest harvest. Soil Biol. Biochem. 2018, 119, 184–193. [Google Scholar]
  18. Olk, D.C.; Samson, M.I.; Gapas, P. Inhibition of nitrogen mineralization in young humic fractions by anaerobic decomposition of rice crop residues. Eur. J. Soil Sci. 2007, 58, 270–281. [Google Scholar] [CrossRef]
  19. Cucu, M.A.; Pullicino, D.; Maurino, V.; Bonifacio, E.; Romani, M.; Celi, L. Influence of redox conditions and rice straw incorporation on nitrogen availability in fertilized paddy soils. Biol. Fertil. Soils 2014, 50, 75–764. [Google Scholar] [CrossRef]
  20. Xia, L.L.; Lam, S.K.; Wolf, B.; Kiese, R.; Chen, D.L.; Butterbach-Bahl, K. Trade-offs between soil carbon sequestration and reactive nitrogen losses under straw return in global agroecosystems. Glob. Chang. Biol. 2018, 24, 5919–5932. [Google Scholar] [CrossRef]
  21. Breulmann, M.; Schulz, E.; Weißhuhn, K.; Buscot, F. Impact of the plant community composition on labile soil organic carbon, soil microbial activity and community structure in semi-natural grassland ecosystems of different productivity. Plant Soil 2012, 352, 253–265. [Google Scholar] [CrossRef]
  22. Xia, L.L.; Lam, S.K.; Yan, X.Y.; Chen, D.L. How Does Recycling of Livestock Manure in Agroecosystems Affect Crop Productivity, Reactive Nitrogen Losses, and Soil Carbon Balance? Environ. Sci. Technol. 2017, 51, 7450–7457. [Google Scholar] [CrossRef] [PubMed]
  23. Bai, Z.; Caspari, T.; Gonzalez, M.R.; Batjes, N.H.; Toth, Z. Effects of agricultural management practices on soil quality. Agric. Ecosyst. Environ. 2018, 265, 1–7. [Google Scholar] [CrossRef]
  24. Tardy, V.; Chabbi, A.; Charrier, X.; Berranger, C.D.; Reignier, T.; Dequiedt, S.; Faivre-Primot, C.; Terrat, S.; Ranjard, L.; Maron, P. Land use history shifts in situ fungal and bacterial successions following wheat straw input into the soil. PLoS ONE 2015, 10, e0130672. [Google Scholar] [CrossRef] [Green Version]
  25. Chen, Y.L.; Xin, L.; Liu, J.T.; Yuan, M.Z.; Liu, S.T.; Jiang, W.; Chen, J.P. Changes in bacterial community of soil induced by long-term straw returning. Sci. Agric. 2017, 74, 349–356. [Google Scholar] [CrossRef]
  26. Zhao, X.; Wang, S.Q.; Xing, G.X. Nitrification, acidification, and nitrogen leaching from subtropical cropland soils as affected by rice straw-based biochar: Laboratory incubation and column leaching studies. J. Soils Sediments 2014, 14, 471–482. [Google Scholar] [CrossRef]
  27. Huang, X.Z.; Wang, C.; Liu, Q.; Zhu, Z.K.; Lynn, T.M.; Shen, J.L.; Whiteley, A.S.; Kumaresan, D.; Ge, T.D.; Wu, J.S. Abundance of microbial CO2-fixing genes during the late rice season in a long-term management paddy field amended with straw and straw-derived biochar. Can. J. Soil Sci. 2018, 98, 306–316. [Google Scholar] [CrossRef]
  28. Bu, R.Y.; Ren, T.; Lei, M.J.; Liu, B.; Li, X.K.; Cong, R.H.; Zhang, Y.Y.; Lu, J.W. Tillage and straw-returning practices effect on soil dissolved organic matter, aggregate fraction and bacteria community under rice-rice-rapeseed rotation system. Agric. Ecosyst. Environ. 2020, 287. [Google Scholar] [CrossRef]
  29. Miao, S.J.; Zhao, H.F.; Pan, S.Q.; Zhang, N.; Qiao, Y.F. Elevated CO2 weakens the shift in bacterial community structure in response to 8-year soybean straw return in the same experiment. Int. J. Phytoremed. 2020, 1–6. [Google Scholar] [CrossRef]
  30. Pronk, G.J.; Heister, K.; Kögel-Knabner, I. Iron oxides as major available interface component in loamy arable topsoils. Soil Sci. Soc. Am. J. 2011, 75, 2158–2168. [Google Scholar] [CrossRef]
  31. Wang, J.; Muyzer, G.; Bodelier, P.L.E.; Laanbroek, H.J. Diversity of iron oxidizers in wetland soils revealed by novel 16S rRNA primers targeting Gallionella-related bacteria. ISME J. 2009, 3, 715–725. [Google Scholar] [CrossRef] [PubMed]
  32. Hartland, A.; Andersen, M.S.; Hamilton, D.P. Phosphorus and arsenic distributions in a seasonally stratified, iron- and manganese-rich lake: Microbiological and geochemical controls. Environ. Chem. 2015, 12, 708–724. [Google Scholar] [CrossRef]
  33. Pan, W.; Kan, J.; Inamdar, S.; Chen, C.; Sparks, D. Dissimilatory microbial iron reduction release doc (dissolved organic carbon) from carbon-ferrihydrite association. Soil Biol. Biochem. 2016, 103, 232–240. [Google Scholar] [CrossRef]
  34. Wang, N.; Yu, J.G.; Zhao, Y.H.; Chang, Z.Z.; Shi, X.X.; Ma, L.Q.; Li, H.B. Straw enhanced CO2 and CH4 but decreased N2O emissions from flooded paddy soils: Changes in microbial community compositions. Atmos. Environ. 2018, 174, 171–179. [Google Scholar] [CrossRef]
  35. Hu, J.; Wu, H.; Sun, Z.; Peng, Q.A.; Hu, R. Ferrous iron (Fe2+) addition decreases methane emissions induced by rice straw in flooded paddy soils. ACS Earth Space Chem. 2020, 4, 843–853. [Google Scholar] [CrossRef]
  36. Violante, A.; Barberis, E.; Pigna, M.; Boero, V. Factors affecting the formation, nature, and properties of iron precipitation products at the soil–root interface. J. Plant Nutr. 2003, 26, 1889–1908. [Google Scholar] [CrossRef]
  37. Xu, X.W.; Wang, P.; Zhang, J.; Chen, C.; Wang, Z.P.; Kopittke, P.M.; Kretzschmar, R.; Zhao, F.J. Microbial sulfate reduction decreases arsenic mobilization in flooded paddy soils with high potential for microbial Fe reduction. Environ. Pollut. 2019, 251, 952–960. [Google Scholar] [CrossRef]
  38. Weber, K.A.; Achenbach, L.A.; Coates, J.D. Microorganisms pumping iron: Anaerobic microbialiron oxidation and reduction. Nat. Rev. Microbiol. 2006, 4, 752. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Li, H.; Su, J.Q.; Yang, X.R.; Zhou, G.W.; Lassen, S.B.; Zhu, Y.G. RNA stable isotope probing of potential Feammox population in paddy soil. EST 2019, 53, 4841–4849. [Google Scholar] [CrossRef]
  40. Ding, B.J.; Zhang, H.; Luo, W.Q.; Sun, S.Y.; Cheng, F.; Li, Z.K. Nitrogen loss through denitrification, anammox and Feammox in a paddy soil. Sci. Total Environ. 2021, 145601. [Google Scholar] [CrossRef] [PubMed]
  41. Tan, X.; Xie, G.J.; Nie, W.B.; Xing, D.F.; Liu, B.F.; Ding, J.; Ren, N.Q. Fe(III)-mediated anaerobic ammonium oxidation: A novel microbial nitrogen cycle pathway and potential applications. Crit. Rev. Environ. Sci. Technol. 2021, 1–33. [Google Scholar] [CrossRef]
  42. Zhu, Y.G.; Duan, G.L.; Chen, B.D.; Peng, X.H.; Chen, Z.L.; Sun, G.X. Mineral weathering and element cycling in soil-microorganism-plant system. Sci. China Earth Sci. 2014, 57, 888–896. [Google Scholar] [CrossRef]
  43. Xia, L.L.; Lam, S.K.; Wang, S.W.; Zhou, W.; Chen, D.L.; Yan, X.Y. Optimizing nitrogen fertilization rate to enhance soil carbon storage and decrease nitrogen pollution in paddy ecosystems with simultaneous straw incorporation. Agric. Ecosyst. Environ. 2020, 298. [Google Scholar] [CrossRef]
  44. Huang, T.; Yang, H.; Huang, C.; Ju, X. Effects of nitrogen management and straw return on soil organic carbon sequestration and aggregate-associated carbon. Eur. J. Soil Sci. 2018, 69, 913–923. [Google Scholar] [CrossRef]
  45. Bao, S.D. Soil Agro-chemistrical Analysis, 3rd ed.; China Agriculture Press: Beijing, China, 2007; pp. 268–270, 389–391. [Google Scholar]
  46. Kostka, J.E.; Iii, G.W.L. Partitioning and speciation of solid phase iron in saltmarsh sediments. Geochim. Cosmochim. Acta 1994, 58, 1701–1710. [Google Scholar] [CrossRef]
  47. Mehra, O.P.; Jackson, M.L. Iron oxide removal from soils and clays by a dithionite-citrate system buffered with sodium bicarbonate. Clays Clay Miner 1960, 7, 317–327. [Google Scholar] [CrossRef]
  48. Dennis, K.L.; Wang, Y.; Blatner, N.R.; Wang, S.; Abdulrahman, S.; Trudeau, E.; Roers, A.; Weaver, C.T.; Lee, J.J.; Gilbert, J.A. Adenomatous polyps are driven by microbe-instigated focal inflammation and are controlled by IL-10-producing T cells. Cancer Res. 2013, 73, 5905–5913. [Google Scholar] [CrossRef] [Green Version]
  49. Schloss, P.D.; Westcott, S.L.; Ryabin, T.; Hall, J.R.; Hartmann, M.; Hollister, E.B.; Lesniewski, R.A.; Oakley, B.B.; Parks, D.H.; Robinson, C.J.; et al. Introducing Mothur: Open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 2009, 75, 7537–7541. [Google Scholar] [CrossRef] [Green Version]
  50. Stackebrandt, E.; Goebel, B. Taxonomic note: A place for DNA-DNA reassociation and 16S rRNA sequence analysis in the present species definition in bacteriology. Int. J. Syst. Bacteriol. 1994, 44, 846–849. [Google Scholar] [CrossRef] [Green Version]
  51. Kumar, K.; Goh, K.M. Nitrogen release from crop residues and organic amendments as affected by biochemical composition. Commun. Soil Sci. Plan. 2003, 34, 2441–2460. [Google Scholar] [CrossRef]
  52. Wang, S.C.; Zhao, Y.W.; Wang, J.Z.; Zhu, P.; Cui, X.; Han, X.Z.; Xu, M.G.; Lu, C.A. The efficiency of long-term straw return to sequester organic carbon in Northeast China’s cropland. J. Integr. Agric. 2018, 17, 436–448. [Google Scholar] [CrossRef]
  53. Li, S.; Li, M.; Wang, G.; Sun, X.; Hu, Z. Compositional and chemical characteristics of dissolved organic matter in various types of cropped and natural Chinese soils. Chem. Biol. Technol. Agric. 2019, 6, 1–11. [Google Scholar] [CrossRef]
  54. Zhang, J.; Hu, K.; Li, K.; Zheng, C.; Li, B. Simulating the effects of long-term discontinuous and continuous fertilization with straw return on crop yields and soil organic carbon dynamics using the DNDC model. Soil Tillage Res. 2017, 165, 302–314. [Google Scholar] [CrossRef]
  55. Xu, X.; Pang, D.W.; Chen, J.; Luo, Y.L.; Zheng, M.J.; Yin, Y.P.; Li, Y.X.; Li, Y.; Wang, Z.L. Straw return accompany with low nitrogen moderately promoted deep root. Field Crops Res. 2018, 221, 71–80. [Google Scholar] [CrossRef]
  56. Zeng, X.N.; Song, Q.L.; Yu, S.; Zhou, M.; Feng, Y.J.; Lai, Y.C.; Wang, Q. Effects of fertilization and straw incorporation on bacterial communities in black soil, Northeastern China. Appl. Ecol. Environ. Res. 2020, 18, 7535–7550. [Google Scholar] [CrossRef]
  57. Fierer, N.; Christian, L.L.; Ramirez, K.S.; Zaneveld, J.; Bradford, M.A.; Knight, R. Comparative metagenomic, phylogenetic and physiological analyses of soil microbial communities across nitrogen gradients. ISME J. 2012, 6, 1007–1017. [Google Scholar] [CrossRef] [Green Version]
  58. Geisseler, D.; Scow, K.M. Long-term effects of mineral fertilizers on soil microorganisms—A review. Soil Biol. Biochem. 2014, 75, 54–63. [Google Scholar] [CrossRef]
  59. Suding, K.N.; Collins, S.L.; Gough, L.; Clark, C.; Cleland, E.E.; Gross, K.L.; Milchunas, D.G.; Pennings, S.; Chapin, F.S. Functional-and abundance-based mechanisms explain diversity loss due to N fertilization. Proc. Natl. Acad. Sci. USA 2005, 102, 4387–4392. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Ramirez, K.S.; Craine, J.M.; Fierer, N. Consistent effects of nitrogen amendments on soil microbial communities and processes across biomes. Glob. Chang. Biol. 2012, 18, 1918–1927. [Google Scholar] [CrossRef]
  61. Zhu, J.; Peng, H.; Ji, X.H.; Li, C.J.; Li, S.N. Effects of reduced inorganic fertilization and rice straw recovery on soil enzyme activities and bacterial community in double-rice paddy soils. Eur. J. Soil Biol. 2019, 94. [Google Scholar] [CrossRef]
  62. Liu, L.; Ding, M.J.; Zhou, L.K.; Chen, Y.; Li, H.P.; Zhang, F.M.; Li, G.; Zhou, Z.F.; Zhang, Y.; Zhou, X.X. Effects of different rice straw on soil microbial community structure. Agron. J. 2020. [Google Scholar] [CrossRef]
  63. Zhou, G.P.; Gao, S.J.; Chang, D.N.; Rees, R.M.; Cao, W.D. Using milk vetch (Astragalus sinicus L.) to promote rice straw decomposition by regulating enzyme activity and bacterial community. Bioresour. Technol. 2021, 319, 124–215. [Google Scholar] [CrossRef] [PubMed]
  64. Xu, Y.; Ge, Y.; Song, J.; Christopher, R. Assembly of root-associated microbial community of typical rice cultivars in different soil types. Biol. Fert. Soils 2020, 56, 1–12. [Google Scholar] [CrossRef]
  65. Bastida, F.; Torres, I.F.; Hernández, T.; García, C. The impacts of organic amendments: Do they confer stability against drought on the soil microbial community? Soil Biol. Biochem. 2017, 113, 173–183. [Google Scholar] [CrossRef]
  66. Liu, J.J.; Yao, Q.; Li, Y.S.; Zhang, W.; Wang, G.H. Continuous cropping of soybean alters the bulk and rhizospheric soil fungal communities in a Mollisol of Northeast PR China. Land Degrad. Dev. 2019, 30, 725–1738. [Google Scholar] [CrossRef]
  67. Paterson, E.; Osler, G.; Dawson, L.A.; Gebbing, T.; Sim, A.; Ord, B. Labile and recalcitrant plant fractions are utilised by distinct microbial communities in soil: Independent of the presence of roots and mycorrhizal fungi. Soil Biol. Biochem. 2008, 40, 1103–1111. [Google Scholar] [CrossRef]
  68. Baker, B.J.; Sheik, C.S.; Taylor, C.A.; Jain, S.; Bhasi, A.; Cavalcoli, J.D.; Dick, G.J. Community transcriptomic assembly reveals microbes that contribute to deep-sea carbon and nitrogen cycling. ISME J. 2013, 7, 1962–1973. [Google Scholar] [CrossRef] [Green Version]
  69. Harris, J. Soil microbial communities and restoration ecology: Facilitators or followers? Science 2009, 325, 573–574. [Google Scholar] [CrossRef]
  70. Jeong, S.; Hong, J.K.; Jho, E.H.; Nam, K. Interaction among soil physicochemical properties, bacterial community structure, and arsenic contamination: Clay-induced change in long-term arsenic contaminated soils. J. Hazard. Mater. 2019, 378, 120729. [Google Scholar] [CrossRef]
  71. Li, X.; Mou, S.; Chen, Y.; Liu, T.; Dong, J.; Li, F. Microaerobic Fe(II) oxidation coupled to carbon assimilation processes driven by microbes from paddy soil. Sci. China-Earth Sci. 2019, 62, 1719–1729. [Google Scholar] [CrossRef]
Figure 1. Geographical location of the experimental site. (A) Map of Jiangsu Province, China; (B) map of Zhenjiang City; (C) photos of paddy field.
Figure 1. Geographical location of the experimental site. (A) Map of Jiangsu Province, China; (B) map of Zhenjiang City; (C) photos of paddy field.
Microorganisms 09 00988 g001
Figure 2. Relative abundances of major bacterial taxa present at the phylum (A) and family (B) levels in paddy soils under different treatments. Only phylogenetic groups represented by more than 1% of the total reads at the family level are shown. T1, nitrogen fertilization alone (80 kg ha−1); ST0, straw incorporation alone (6000 kg ha−1); ST1, straw incorporation with low nitrogen (80 kg ha−1); ST2, straw incorporation with moderate nitrogen (300 kg ha−1); ST3, straw incorporation with high nitrogen (450 kg ha−1).
Figure 2. Relative abundances of major bacterial taxa present at the phylum (A) and family (B) levels in paddy soils under different treatments. Only phylogenetic groups represented by more than 1% of the total reads at the family level are shown. T1, nitrogen fertilization alone (80 kg ha−1); ST0, straw incorporation alone (6000 kg ha−1); ST1, straw incorporation with low nitrogen (80 kg ha−1); ST2, straw incorporation with moderate nitrogen (300 kg ha−1); ST3, straw incorporation with high nitrogen (450 kg ha−1).
Microorganisms 09 00988 g002
Figure 3. Results of redundancy analysis of bacterial community structure (OTU assignment: distance cutoff 0.03) retrieved from Illumina sequencing of the 16S rRNA genes in relation to soil properties and depth level. SOC: soil organic carbon; EC: electrical conductivity; TN: total nitrogen; TP: total phosphorus; WC: water content; Fe(II): ferrous iron; Fe(III): ferric iron; DCB-Fe: dithionite-citrate-bicarbonate-extracted iron. Treatment abbreviations are defined in Figure 2 legend.
Figure 3. Results of redundancy analysis of bacterial community structure (OTU assignment: distance cutoff 0.03) retrieved from Illumina sequencing of the 16S rRNA genes in relation to soil properties and depth level. SOC: soil organic carbon; EC: electrical conductivity; TN: total nitrogen; TP: total phosphorus; WC: water content; Fe(II): ferrous iron; Fe(III): ferric iron; DCB-Fe: dithionite-citrate-bicarbonate-extracted iron. Treatment abbreviations are defined in Figure 2 legend.
Microorganisms 09 00988 g003
Figure 4. Kyoto Encyclopedia Genes and Genomes-predicted metabolic pathways (A) and energy metabolism pathways (B) of soil bacterial communities in different treatments. Treatment abbreviations are defined in Figure 2 legend.
Figure 4. Kyoto Encyclopedia Genes and Genomes-predicted metabolic pathways (A) and energy metabolism pathways (B) of soil bacterial communities in different treatments. Treatment abbreviations are defined in Figure 2 legend.
Microorganisms 09 00988 g004
Table 1. Fertilization treatments applied in the experiments (unit: kg ha−1).
Table 1. Fertilization treatments applied in the experiments (unit: kg ha−1).
Treatments Wheat Straw Nitrogen in the Form of UreaPhosphorus in the Form of Calcium SuperphosphatePotassium in the Form Potassium Chloride
T10808080
ST0600008080
ST16000808080
ST260003008080
ST36000450 *8080
* Local farmer practice.
Table 2. Summary of the major soil physicochemical properties in different treatments.
Table 2. Summary of the major soil physicochemical properties in different treatments.
TreatmentspHECSOM
(g kg−1)
TP
(g kg−1)
TN
(g kg−1)
WC
(%)
Fe(II)
(g kg−1)
Fe(III)
(g kg−1)
DCB-Fe
(g kg−1)
T15.05 ± 0.08 a 154.00 ± 7.00 b32.29 ± 2.61 a 0.60 ± 0.03 a 1.75 ± 0.15 a 65.94 ± 9.31 a3.81 ± 0.17 a 13.12 ± 1.34 a74.38 ± 1.08 a
ST05.06 ± 0.02 a 164.67 ± 23.01 ab27.13 ± 1.79 b 0.58 ± 0.26 a 1.57 ± 0.26 ab 58.16 ± 9.51 ab3.95 ± 1.66 a10.39 ± 1.58 a79.24 ± 7.96 a
ST15.11 ± 0.22 a188.33 ± 16.07 a27.13 ± 0.40 b 0.57 ± 0.04 a 1.46 ± 0.07 b47.71 ± 1.95 b3.53 ± 1.28 a11.84 ± 3.49 a77.18 ± 3.08 a
ST25.05 ± 0.09 a 158.00 ± 14.11 ab30.21 ± 1.71 ab0.57 ± 0.02 a1.40 ± 0.01 b45.95 ± 4.36 b2.57 ± 0.17 a10.53 ± 3.17 a78.47 ± 11.07 a
ST35.00 ± 0.15 a 179.67 ± 19.09ab30.99 ± 2.86 ab0.55 ± 0.00 a1.47 ± 0.08 b50.84 ± 6.84 b4.86 ± 1.89 a14.17 ± 2.34 a77.92 ± 4.44 a
EC: electrical conductivity; SOC: soil organic matter; TP: total phosphorus; TN: total nitrogen; WC: gravimetric water content; DCB-Fe: dithionite-citrate-bicarbonate-extracted iron. Treatment abbreviations are defined in Table 1. Different letters in the same column indicate statistically significant differences among the samples by one-way ANOVA (Tukey–Kramer, p < 0.05).
Table 3. Alpha-diversity and richness estimates of soil microbial communities based on 97% similarity OTU clusters.
Table 3. Alpha-diversity and richness estimates of soil microbial communities based on 97% similarity OTU clusters.
TreatmentsSobsACEShannonCoverage
T12841.3 ± 96.03 b3812.9 ± 47.71 ab6.86 ± 0.02 a0.959
ST02726.0 ± 85.58 ab3759.6 ± 188.9 a6.85 ± 0.07 ab0.966
ST12771.3 ± 101.53 ab3712 ± 99.25 ab6.79 ± 0.04 ab0.968
ST22970.7 ± 100.55 a3843.8 ± 99.4 a6.81 ± 0.14 ab0.972
ST32667.0 ± 74.73 b3620.2 ± 39.83 b6.66 ± 0.09 b0.964
Treatment abbreviations are defined in Table 1. Different letters in the same column indicate statistically significant differences among the samples based on one-way ANOVA (Tukey–Kramer, p < 0.05).
Table 4. Spearman correlation coefficients between the relative abundances of dominant bacterial phyla and soil physiochemical properties.
Table 4. Spearman correlation coefficients between the relative abundances of dominant bacterial phyla and soil physiochemical properties.
PhylumpHECSOCTPTNWCFe(II)Fe(III)DCB-Fe
Proteobacteria0.24 0.21 0.00 −0.63 * −0.78 *** −0.58 * −0.46 −0.23 0.14
Chloroflexi0.24 0.06 −0.54 * −0.02 −0.57 * −0.31 −0.48 −0.71 ** −0.19
Acidobacteria0.12 0.15 −0.33 −0.04 −0.70 ** −0.49 −0.51 −0.32 −0.03
Nitrospirae−0.05 0.30 0.20 −0.33 −0.40 −0.46 −0.41 0.10 0.44
Actinobacteria0.44 0.07 −0.61 * −0.29 −0.61 * −0.35 −0.36 −0.85 *** −0.07
Bacteroidetes−0.31 −0.36 0.51 −0.40 −0.27 0.19 −0.14 0.07 −0.12
Gemmatimonadetes0.42 0.50 −0.62 * −0.45 −0.74 ** −0.80 *** −0.24 −0.41 0.20
Verrucomicrobia−0.53 * 0.25 0.08 0.39 0.09 −0.08 0.00 0.69 ** 0.39
Ignavibacteriae0.14 0.13 −0.10 −0.57 * −0.80 *** −0.36 −0.33 −0.12 0.05
Parcubacteria0.17 0.48 * −0.28 0.11 −0.31 −0.43 −0.04 0.32 0.46
Unclassified_k_norank_d_Bacteria0.15 0.50 * −0.63 0.01 −0.63 * −0.63 * −0.66 ** −0.30 0.28
Armatimonadetes−0.04 0.28 −0.11 0.45 −0.03 −0.26 −0.03 0.43 0.33
Latescibacteria0.32 0.58 * −0.46 0.11 −0.29 −0.58 * −0.50 −0.13 0.14
Spirochaetae0.18 −0.34 0.04 −0.54 * −0.52 * 0.02 −0.27 −0.44 −0.55 *
Planctomycetes0.03 0.54 * −0.54 * 0.13 −0.32 −0.41 0.10 0.27 0.21
EC: electrical conductivity; SOC: soil organic carbon; TP: total phosphorus; TN: total nitrogen; WC: gravimetric water content; Fe(II), ferrous iron; Fe(III), ferric iron; DCB-Fe: dithionite-citrate-bicarbonate-extracted iron. Asterisks indicate significant difference (* 0.01 < p ≤ 0.05, ** 0.001 < p ≤ 0.01, *** p ≤ 0.001).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, J.; Ma, Y.; Di, L.; Qian, X.; Wang, G. Straw Incorporation with Nitrogen Amendment Shapes Bacterial Community Structure in an Iron-Rich Paddy Soil by Altering Nitrogen Reserves. Microorganisms 2021, 9, 988. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms9050988

AMA Style

Wang J, Ma Y, Di L, Qian X, Wang G. Straw Incorporation with Nitrogen Amendment Shapes Bacterial Community Structure in an Iron-Rich Paddy Soil by Altering Nitrogen Reserves. Microorganisms. 2021; 9(5):988. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms9050988

Chicago/Turabian Style

Wang, Juanjuan, Yao Ma, Lin Di, Xiaoqing Qian, and Guiliang Wang. 2021. "Straw Incorporation with Nitrogen Amendment Shapes Bacterial Community Structure in an Iron-Rich Paddy Soil by Altering Nitrogen Reserves" Microorganisms 9, no. 5: 988. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms9050988

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop