Fangfang Bai
Wei Guo1,2†
Ping Li- 1Institute of Farmland Irrigation, Chinese Academy of Agricultural Sciences, Xinxiang, China
- 2Water Environment Factor Risk Assessment Laboratory of Agricultural Products Quality and Safety, Ministry of Agriculture and Rural Affairs, Xinxiang, China
- 3Agricultural Water Soil Environmental Field Research Station of Xinxiang, Chinese Academy of Agricultural Sciences, Xinxiang, China
Introduction: It is well known that reduced nitrogen application and groundwater depth can change soil microbial communities, but the associated difference in the response of abundant and rare bacterial composition to these local environmental changes remains unclear.
Methods: In this study a lysimeter experiment was carried out to examine the impact of reduced nitrogen and groundwater depth on the composition of abundant and rare bacteria.
Results and discussion: Our results demonstrated that the summer maize field soil species composition of rare bacterial sub-communities was significantly regulated by reduced nitrogen application, groundwater depth change and their interactions. However, only reduced nitrogen application had a significant influence on the species composition of abundant bacterial sub-communities. The structural equation model (SEM) indicated that reduced nitrogen application and groundwater depth change also could indirectly regulate the species composition of abundant and rare bacteria by altering soil attributes. The changes in soil pH and TSN had the most significant effects on the community composition of abundant and rare bacteria, respectively. More importantly, rare bacterial sub-communities were more sensitive to the changes in nitrogen input, groundwater depth and soil factors. Collectively, our study first demonstrated that abundant and rare microbial sub-communities responded differently to reduced nitrogen application and groundwater depth change. This study highlights that summer maize farmland production management should take nitrogen input and groundwater depth into consideration to maintain the compositional stability of soil rare microbial sub-communities.
1. Introduction
Soil microorganisms play a major role in regulating key ecosystem processes and functions (Koranda et al., 2013; Mooshammer et al., 2014). Nitrogen is one of the most important limiting factors for the growth of soil microbes (Sun et al., 2015; Hou et al., 2018; Zhang et al., 2022). Nitrogen fertilizers can not only directly affect soil microbial communities but also indirectly affect soil microbial communities by regulating soil fertility, soil salt content and groundwater salt content (She et al., 2018; Bai et al., 2020; Wu et al., 2020; Bai et al., 2023). However, over-application of nitrogen fertilizers may result in a series of severe environmental problems, including soil salinization, nitrate contamination and groundwater pollution (Zhang et al., 2013; Lasagna et al., 2016). In addition, over-application of nitrogen fertilizers may greatly negatively affect the function of ecosystems by altering soil bacterial diversity and composition (Luo et al., 2016; Cocco et al., 2018; Vergine et al., 2023). Therefore, it is necessary to explore the appropriate nitrogen input amount to maintain sustainable agricultural development (Ju et al., 2009; Liu B.-Y. et al., 2021).
Soil bacterial communities are primarily dominated by a few abundant species, whereas a great number of other species (“rare biosphere”) have extremely low abundance (Pedrós-Alió, 2012; Ji et al., 2020; Wang et al., 2021b). It is widely believed that abundant subcommunities mostly determine ecosystem functioning (Egidi et al., 2019). Most rare communities can be widely present in the microbial communities of various habitats, and they play an important role in regulating ecosystem stability and function (Liu et al., 2015; Jiao et al., 2017). Rare subcommunities account for only the entire microbial community, However, their contribution to the temporal dynamic change of community structure can reach 97% (Shade et al., 2014). The spatial distributions and driving factors of abundant and rare bacteria have been well examined in different ecosystem types, such as agricultural land (Wan et al., 2020; Jiao and Lu, 2020b), forest (Wang et al., 2021a), alpine grassland (Ji et al., 2020), wetland (Wan et al., 2021), and desert (Wang et al., 2021b). Abundant and rare bacterial species may be subjected to different determinants, and thus exhibit distinct ecological patterns (Wan et al., 2020; Jiao and Lu, 2020b; Wang et al., 2021b). Hence, comparing the response of abundant and rare bacterial sub-communities to reduced nitrogen application may help us better predict how soil biodiversity and function respond to global change. However, how differently reduced nitrogen application affects the species composition of abundant and rare bacteria in agricultural soils remains largely unclear.
As one of the most important crop-producing regions of China, the North China Plain has a total area of 300,000 km2. In the past decades, its groundwater consumption increased sharply because of population growth, urbanization and development of industry and agriculture (Changming et al., 2001; Chen et al., 2020; Jing et al., 2023; Xiao et al., 2023), causing various environmental geologic problems, including groundwater table fall, water quality worsening and surface subsidence (Changming et al., 2001; Shi et al., 2020; Gan et al., 2022). Notably, groundwater table can regulate soil microbial species richness and composition (Bai et al., 2020; Wu et al., 2020; Zhang et al., 2020). Besides, groundwater depth change can indirectly influence soil microbial communities by altering soil attributes such as salinity, pH and nutrients (Bai et al., 2020; Zhang et al., 2020). The interactions between reduced nitrogen application and groundwater depth have a significant influence on soil microbial communities (Bai et al., 2020). Furthermore, over-nitrogen form release may lead to severe groundwater contamination (Matzeu et al., 2017; Paladino et al., 2020; Abbasi and Sepaskhah, 2023). However, to date, the effect of reduced nitrogen application, groundwater depth and their interactions on the species composition of abundant and rare bacteria has been rarely studied.
Soil microbial communities may be determined by a wide range of soil attributes such as soil pH, salinity and nutrients (Bai et al., 2020; Wan et al., 2020; Jiao and Lu, 2020b; Wang et al., 2021b). In addition, the effects of different soil factors vary significantly between abundant and rare sub-communities (Jiao and Lu, 2020b; Wang et al., 2021b). For example, the assemblies of abundant and rare sub-communities are regulated by soil pH and soil moisture content, respectively (Wang et al., 2021a). Moreover, nitrogen application amount and groundwater depth and their interactions can cause substantial changes in soil attributes (Bai et al., 2020). We therefore posited that abundant and rare bacteria had different responses to the alteration in soil attributes caused by reduced nitrogen application and groundwater depth. However, few studies have focused on the direct and indirect effect of reduced nitrogen application, groundwater depth and soil attributes on abundant and rare bacteria in agricultural soils.
Therefore, the major objective of our study was to examine the associated difference in the response of abundant and rare bacteria to reduced nitrogen application and groundwater depth in agricultural soils. In this study, a large lysimeter experiment was conducted. Soil abundant and rare bacterial sub-communities were assessed based on 16S rRNA gene sequence data. Specifically, this study was designed to address the following questions: (1) How does the species composition of abundant and rare bacteria respond to reduced nitrogen application, groundwater depth and their interactions? (2) What are the direct and indirect influences of reduced nitrogen application, groundwater depth change and soil factors on soil abundant and rare bacterial composition?
2. Materials and methods
2.1. Experimental site and design
The experiment was carried out in the long-term monitoring lysimeter at the Xinxiang Agricultural Water and Soil Environment Field Scientific Observation and Experiment Station, Chinese Academy of Agricultural Sciences (35°19′N, 113°53′E), Xinxiang City, Henan Province, People’s Republic of China (Figure 1). The planting mode was winter wheat-summer maize two-cropping rotation a year. Its mean annual air temperature and annual precipitation is 14.1°C, 588.8 mm, respectively. Almost all precipitation is concentrated in July, August and September. The mean annual evapotranspiration is approximately 2,000 mm. The experiments were conducted on a large lysimeter platform in 2019, and each lysimeter was evenly filled layer by layer with silt loam soil collected from the neighboring cropland. The groundwater level in the micro-lysimeters was controlled at the constant depths of 2, 3, and 4 m from the soil surface by Mariotte flash.
Figure 1. Above ground part of the lysimeter.
The field site was sown with maize (Zea mays L.), the maize variety tested was Huaiyu 208, which was selected and bred by Henan Huaichuan Seed Industry Co., Ltd., and its approval number was Yu Shenyu 2,013,008, where we examined six treatments which were produced by the interactions between two levels of N input [local conventional N of 300 kg N ha−1 y−1 (N300), reduced N of 240 kg N ha−1 y−1 (N240)] and three groundwater depth levels (GW2, 2 m; GW3, 3 m and GW4, 4 m) (Table 1). Local conventional nitrogen level (300 kg N ha−1 y−1) was used as a control with 100% of normal nitrogen input. Maize was sown on June 10, 2020 and harvested on September 23, with a growth period of 105 days. Forty five millimeter water were irrigated immediately after sown. Excessively shallow groundwater depth will limit the growth of crops and root aerobic respiration, and lead to soil waterlogging (Askri et al., 2014; Soylu et al., 2014). A previous study has observed that 2.5 m is a crucial groundwater depth to control soil salinization (Askri et al., 2014). When the groundwater table is shallow, there are obvious two-ways of water exchange at the interface under the root layer, the upper soil moisture can infiltrate into deep soil and ground water, and crop evapotranspiration can partially utilize the groundwater. When the groundwater table becomes deeper, the soil moisture at the interface mainly migrates downward, and the critical groundwater depth for crop growth in the North China Plain is generally 4 m. Therefore, we used 2 m, 3 m and 4 m as the treatments of groundwater depth in our tests. The experiment followed a random block design with four replicates. Furthermore, 120 kg K2O ha−1 and 150 kg P2O5 ha−1 were applied in six treatments in the form of potassium sulfate and monocalcium phosphate, respectively. Nitrogen fertilizer was urea with 46% nitrogen content. Base fertilizer were applied on June 29 of 2020 (14, BBCH-scale), and topdressing were applied on July 29 of 2020 (51, BBCH-scale). The ratio of base fertilizer and topdressing of nitrogen was 4:6, phosphate fertilizer and potassium fertilizer were applied as base fertilizer at one time. Fertilization immediately followed by irrigation. And ground irrigation was conducted when the soil moisture was lower than 55% of the field water holding rate. A RS-XAJ-100 probe was buried in the 20 cm soil layer to monitor soil moisture and guide irrigation. The whole growth period of maize was irrigated 45 mm on June 10, June 29, July 15, July 29, August 10 and September 4, respectively.
Table 1. The experimental design.
2.2. Sampling and soil chemical properties determination
At the September 24 of 2020, five soil samples (0–20 cm depth) were randomly collected from each lysimeter using an auger and then mixed into a composite sample. A total of 24 composite samples were collected in this study. After being sieved through a 2 mm mesh to remove roots and stones, each composite sample was divided into two following parts: one portion was stored in thermally insulated boxes to measure soil attributes, and the other was used to extract DNA. Afterwards, soil pH (pH), organic carbon content (TOC), total nitrogen content (TSN) and total phosphor (TSP) were determined following the previous methodology (Wang et al., 2020). Soil available phosphorus (AP) was extracted from 0.5 M NaHCO3 and determined by a colorimetric method. Meanwhile, soil electrical conductivity (EC) was measured by soil electrical conductivity. Finally, soil NH4+-N(NHN) and NO3−-N(NON) were determined based on the previous methodology (Li et al., 2019).
2.3. Molecular and bioinformatics analysis
The bacterial DNA was extracted from 0.5 g of well-mixed fresh soil samples using E.Z.N.A. soil DNA kits (OMEGA, United States) following the manufacturer’s instructions. The V3-V4 hypervariable regions of the bacterial 16S rRNA gene were amplified using the universal primers 338F (5′-ACTCCTACGGGAGGCAGCAG-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′). These primers contained a set of 8-nucleotide barcodes sequence unique to each sample. The PCR program was as follows: 94°C for 5 min, 25 cycles at 94°C for 30 s, 55°C for 30 s and 72°C for 30 s with a final extension of 72°C for 10 min. Moreover, PCR amplicons were extracted from 2% agarose gels and purified using the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, United States) based on the manufacturer’s instructions, and quantified using QuantiFluor™ -ST (Promega, United States). PCR made with distilled water without a soil sample DNA template were set as a negative control NTC, if the negative control has amplified bands, this batch of PCR will be redone. Purified amplicons were pooled in equimolar and paired-end sequenced (2 × 300) on an Illumina Miseq PE300 sequencing platform at Beijing Allwegene Tech, Ltd. (Beijing, China).
Bacterial sequences >200 bp and with an average quality score > 20 and no ambiguous base calls were processed within the QIIME package. These high-quality sequences were clustered into operational taxonomic units (OTUs) based on a 97% similarity threshold using uparse. The taxonomy of each 16S gene sequence was analyzed by comparing it against sequences within the SILVA reference database (v.138). Subsequently, we removed the that contained less than 20 reads to avoid random effects on the identification of rare taxa (Jiao and Lu, 2020b; Wan et al., 2021). Based on minimum number of sample sequences, sequences were rarefied at 31,124 sequences from each sample to eliminate the effect of sequencing depths on the analyses. The rare and abundant OTUs were defined following previous study (Jiao and Lu, 2020b; Wan et al., 2021). OTUs with a relative abundance of >0.1% of the total abundance were defined as abundant bacteria, and rare bacteria were those with a relative abundance of <0.01% (Mo et al., 2018; Wang et al., 2021b; Ren et al., 2022).
2.4. Statistical analyses
Bray-Curtis community difference matrices of abundant and rare bacteria (the abundance data was Hellinger-transformed) and standardized environmental Euclidean distance matrix for each variable were used for evaluation. The difference distance of Bray-Curtis community was used to reflect the difference in species composition (β-diversity) among soil bacterial communities. Permutational analysis of variance (PERMANOVA) was used to examine the differences in species composition (β-diversity) of abundant and rare bacteria across nitrogen-treatment and groundwater-depth groups. Moreover, the difference in community composition was visualized through Nonmetric multidimensional scaling (NMDS). Levins’ niche breadth (B) index was used to compare the niche breadth of abundant and rare bacterial sub-communities (Levins, 1968). Particularly, the B-value of each bacterial OTU was estimated using the previous methods (Jiao and Lu, 2020a). Abundance-weighted mean B-values across all bacterial OTUs within the community were employed to quantify community B-values (Bcom). The higher the B-value of community, the stronger the metabolic flexibility (Pandit et al., 2009). Notably, the Levins’ niche breadth (B) index was calculated within the “spaa” R package.
Mantel test (10,000 permutations) was used to explore the influence of soil properties on soil bacterial richness. The structural equation model (SEM) was used to analyze the indirect and direct influences of groundwater depth, nitrogen application amount and soil factors on abundant and rare bacterial composition. In order to reduce the strong collinearity between variables, soil pH, TSP and SAP were removed because of the high correlation between these variables (i.e., Pearson’s r > 0.6). Prior to the analysis, we built a prior structural equation model (Supplementary Figure S1) according to relevant theories. The comparative fit index (CFI), goodness of fit index (GFI), low root square mean error of approximation (RMSEA) and χ2-test were used to fit SEM. Additionally, the “lavaan” program was utilized to perform SEM.
3. Results
3.1. Overall distribution patterns of abundant and rare sub-communities
After chimeric sequences were filtered and removed, a total of 746,976 high-quality sequences were clustered into 4,619 bacterial OTUs (Supplementary Table S2), among which only 203 (4.39%) OTUs with 36,2,101 sequences (50.52%) were identified as abundant bacteria, while a total of 3,096 (67.03%) OTUs with 86,975 sequences (12.13%) were identified as rare bacteria. Abundant sub-communities were primarily dominated by Pseudomonadota (29.17%), Actinobacteriota (25.46%), Acidobacteriota (19.55%), and Chloroflexota (11.18%), while rare sub-communities were primarily dominated by Chloroflexota (14.66%), Pseudomonadota (14.21%), and Actinobacteriota (10.48%) (Supplementary Table S3).
Our results showed that the mean Bcom value for rare sub-communities was significantly lower than that for abundant sub-communities (Figure 2A). Moreover, abundant species had a greater presence than rare species across six treatments. Particularly, 99.05% of the abundant OTUs (201 OTUs) occurred in all samples of six treatments, whereas only 0.32% of the rare OTUs (10 OTUs) occurred in all samples of six treatments (Figure 2B). In addition, abundance-occupancy relationships indicated that rare bacteria had a greater positive association than abundant bacteria.
Figure 2. Difference in mean habitat niche breadths (Bcom, which was estimated to reveal the contributions of species sorting and dispersal limitation to microbial community assembly) between abundant and rare bacterial subcommunities (A), and abundance-occupancy relationships for abundant and rare sub-communities (B); ***p < 0.0001; Wilcoxon rank-sum test.
3.2. Difference in the response of abundant and rare bacterial sub-communities to reduced nitrogen application and groundwater depth change
Results of PERMANOVA and NMDS showed that reduced nitrogen application significantly altered the species composition of rare sub-communities (F = 2.405, p < 0.0001), groundwater depth (F = 1.445, p < 0.01) and their interactions (F = 1.468, p < 0.01; Table 2 and Figure 3B). However, the species composition of abundant sub-communities only had a significant response to reduced nitrogen application (F = 3.724, p < 0.0001), but not to groundwater depth (F = 1.533, p > 0.05) or their interactions (F = 1.712, p > 0.05; Table 2 and Figure 3A). Abundant and rare sub-communities exhibited different shifts under N input and groundwater depth levels.
Table 2. Effects of groundwater depth, reduced nitrogen application and their interactions on the species composition of abundant and rare sub-communities.
Figure 3. Nonmetric multidimensional scaling (NMDS) ordination of abundant (A) and rare (B) sub-community compositional differences (Bray-Curtis distance) under different groundwater depths and reduced nitrogen application. GW2N240, reduced N treatment × groundwater depth at 2 m; GW3N240, reduced N treatment × groundwater depth at 3 m;GW4N240, reduced N treatment × groundwater depth at 4 m; GW2N300, control with 100% of normal nitrogen input × groundwater depth at 2 m; GW3N300, control with 100% of normal nitrogen input × groundwater depth at 3 m; GW4N300, control with 100% of normal nitrogen input × groundwater depth at 4 m.
Wilcoxon rank-sum test demonstrated that rare bacterial sub-communities had significantly greater compositional dissimilarity than abundant bacterial sub-communities (Figure 4A, P < 0.001). Mantel correlation analysis revealed that the species composition of abundant bacterial sub-communities had a significant response to the changes in soil pH, EC, SOC, and TSN (p < 0.05), whereas that of rare sub-communities exhibited a significant response to soil pH, EC, SOC, TSN, and NON (p < 0.05) (Figures 4B, 5). Moreover, we found that the strength of responses (i.e., Mantel R) in abundant and rare sub-community dissimilarity to the change in each soil variable was different (Figure 4B). Among all soil variables, soil pH and TSN played a major role in altering abundant and rare sub-community dissimilarity, respectively (Figures 4B, 5A,B).
Figure 4. Difference in the compositional dissimilarity between abundant and rare sub-communities (A) and the response of the compositional dissimilarity of abundant and rare sub-communities to soil factors (B). ***p < 0.001, **p < 0.01, *p < 0.05, NSp > 0.05.
Figure 5. Relationships of soil abundant and rare bacterial compositional dissimilarity to the variations in soil pH (A) and TSN (B), respectively.
3.3. Influence of reduced nitrogen application, groundwater depth change and soil attributes on soil abundant and rare bacterial composition
The fitted SEM showed that reduced nitrogen application, changes in groundwater depth and soil attributes collectively explained 15.80 and 45.50% of the variation in the community composition of abundant and rare bacteria, respectively (Figures 5A,B). Notably, groundwater depth could directly and indirectly affect the species composition of rare bacteria, but only had an indirect impact on that of abundant bacteria. Contrarily, reduced nitrogen application could directly and indirectly affect the species composition of both abundant and rare bacteria. Moreover, standardized total effects derived from the SEM revealed that the species composition of abundant bacteria was primarily determined by reduced nitrogen application, followed by soil pH, TSN and groundwater depth, while that of rare bacteria was primarily regulated by reduced nitrogen application, followed by TSN, groundwater depth, pH, NON, and SOC (Figures 6, 7).
Figure 6. The species composition of abundant (A) and rare (B) subcommunities is described by SEM in terms of the direct and indirect effects of groundwater depth, reduced nitrogen application, and soil properties (Brey-Curtis community differences). Solid red and dashed gray arrows indicate the significant direct paths (p < 0.05) and indirect paths (p < 0.05), respectively.
Figure 7. The standardized total effects (direct plus indirect effects) of groundwater depth, reduced nitrogen application and soil properties on abundant (A) and rare sub-communities (B) derived from the SEM model.
4. Discussion
4.1. Different responses of abundant and rare bacterial composition to reduced nitrogen application and groundwater depth change
Global environmental change such as nitrogen overload has been recognized as one of the major threats to biodiversity and ecosystems. Nitrogen fertilization can significantly alter soil microbial community composition (She et al., 2018; Bai et al., 2020; Liu W. et al., 2021). Moreover, groundwater plays a key role in shaping soil microbial community composition (Zhang et al., 2020). Previous studies have also found that soil rare and abundant microbial taxa adopt different strategies to cope with environmental changes (Hou et al., 2020; Liang et al., 2020; Jiao and Lu, 2020b). However, how differently rare and abundant bacterial sub-communities respond to reduced nitrogen application, groundwater depth change and their interactions remains largely unclear. In this study, we observed that reduced nitrogen application led to a substantial change in both abundant and rare bacterial composition, which is similar to the findings of our study on denitrifier communities (Bai et al., 2020). However, groundwater depth and its interactions with reduced nitrogen application only significantly influenced the community composition of rare bacteria. SEM further demonstrated that reduced nitrogen application could directly and indirectly affect both abundant and rare sub-communities, while groundwater depth had no direct influence on abundant sub-communities. We also observed that reduced nitrogen application and groundwater depth change resulted in a greater shift in rare species composition than in abundant species composition, which is inconsistent with the findings of a previous report regarding soil microbial responses to climatic changes (Liang et al., 2020). This result implies that abundant bacterial composition is more stable to local environmental changes than rare bacterial composition in the agricultural ecosystem. From the above results, we suggest that the response of soil microbes to reduced nitrogen application and groundwater depth change varies between abundant and rare taxa.
Abundant and rare microbial sub-communities respond differently to the environment (Jiao and Lu, 2020b; Wang et al., 2021b). The difference in functional traits between abundant and rare species may mediate their response to environmental conditions by influencing their fitness and performance (McGill et al., 2006; Webb et al., 2010). Abundant species may have broader tolerances (Kurm et al., 2017; Delgado-Baquerizo et al., 2018), which enables them to perform well under different environmental conditions (Jiao and Lu, 2020b). On the contrary, rare species may have worse performance as they are more specialized in rare resources (Umaña et al., 2015). Hence, our study observed that abundant bacteria were distributed more ubiquitously and exhibited a wider habitat niche breadth than rare bacteria. Consequently, abundant species can occupy diverse niches and have greater tolerance and adaptability to environmental changes than rare species (Jiao and Lu, 2020b). Therefore, reduced nitrogen application and groundwater depth change may select species with higher resistance and resource competitiveness within the community, and filter out stress-avoidant rare species, which may account for why rare sub-communities were more sensitive to reduced nitrogen application and groundwater depth change than abundant sub-communities.
Taxa with small populations may be more susceptible to extinction under rapid global environmental changes (Storch et al., 2018). Rare species play key roles in regulating the health and function of ecosystems (Zhang Y. et al., 2019; Xiong et al., 2021). Given that rare sub-communities were more sensitive to reduced nitrogen application and groundwater depth change than abundant sub-communities, we hypothesized that over-nitrogen application, groundwater depth change and their interactions may greatly negatively influence the health and function of the agricultural ecosystem by altering the community composition of rare microbes. Collectively, we highlight that future cropland management should explore the appropriate nitrogen input and groundwater table levels to maintain the compositional stability of soil microbes, especially rare taxa.
4.2. Different influences of the variation in soil attributes on abundant and rare bacterial composition
Both reduced nitrogen application and groundwater depth change can remarkably alter soil attributes, thereby resulting in a shift in soil microbial composition (Bai et al., 2020). Thus, in this study the composition of abundant and rare sub-communities was altered by the variations in soil attributes caused by reduced nitrogen application and groundwater depth change (Supplementary Table S1), which is consistent with the previous studies on denitrifier communities (Bai et al., 2020). In addition, the response direction and strength to soil attributes may vary significantly between abundant and rare bacteria (Jiao and Lu, 2020b; Wang et al., 2021a). Similarly, both abundant and rare bacterial sub-community composition showed the same direction of responses to the changes in soil attributes, while rare bacteria had stronger response strength (i.e., Mantel R) than abundant bacteria. This result suggests that rare bacterial sub-communities may be more sensitive to the changes in soil attributes caused by reduced nitrogen application and groundwater depth change.
Soil microbial communities were controlled by a wide range of soil variables such as soil pH (Tripathi et al., 2018; Jiao and Lu, 2020b), soil salinity (Zhang K. et al., 2019), and soil fertility (Jiao and Lu, 2020a). Among soil attributes, the community composition of rare bacteria was primarily determined by soil TSN, followed by pH, NON, and SOC, while that of abundant bacteria was only shaped by soil pH and TSN, which may suggest that soil factors affect abundant and rare bacterial sub-communities differently. More importantly, both soil pH and TSN had the most important influence on the composition of abundant and rare bacteria, which indicates that reduced nitrogen application and groundwater depth change alter soil bacterial composition primarily by regulating soil pH and nutrients. Soil pH and nutrient availability can directly affect bacterial growth and interact with abiotic and biotic factors to regulate soil bacterial community structure (Delgado-Baquerizo et al., 2018). In addition, soil pH and nutrient availability determine soil microbial community structure in different ecosystem types (Fierer and Jackson, 2006). For example, soil pH and nutrient availability have the most powerful influence on abundant and rare sub-communities in large-scale agricultural ecosystems (Jiao and Lu, 2020a,b). However, our previous study has observed that soil nutrients and salinity play more important roles in shaping the composition of denitrifying communities (Bai et al., 2020). Here are several interpretations of the results. Firstly, the strong covariation among different variables made it difficult to precisely quantify their influence. Secondly, such discrepancies may be also partly due to the differences in the environmental regime between 1- and 2-year experiments and the different sensitivity of microbial taxa to soil attributes (Chase et al., 2017; Ji et al., 2020; Liang et al., 2020). Taken together, we emphasize the different influences of soil-factor variations on abundant and rare bacterial composition.
4.3. Research limitation and future perspectives
This study provided first-hand evidence that abundant and rare bacterial sub-communities had different responses to reduced nitrogen application, groundwater depth change and their interactions. Rare sub-communities were more sensitive to the changes in soil attributes caused by reduced nitrogen application and groundwater depth change. Notably, our study conducted a short-term experiment (2 years) and only selected two levels of N input, which may lead to a probable biased result. Hence, a longer observation period (>5 years) and more nitrogen input levels and microbial taxa should be considered in further research to better elucidate the intrinsic influence of nitrogen fertilization on the agricultural ecosystem, and the appropriate nitrogen input amount should be accurately determined to maintain the health and function of agricultural ecosystems.
5. Conclusion
This study conducted a comprehensive comparison of the response of abundant and rare bacterial sub-communities to reduced nitrogen application, groundwater depth change, their interactions and soil factors. The community composition of rare bacteria was significantly altered by reduced nitrogen application, groundwater depth change and their interactions in summer maize field, while that of abundant bacteria was regulated by reduced nitrogen application alone. Moreover, SEM demonstrated that reduced nitrogen application and groundwater depth change could also alter the species composition of abundant and rare bacteria by causing significant variations in summer maize soil attributes. The variations in soil pH and TSN caused by reduced nitrogen application and groundwater depth played the most important role in altering the community composition of abundant and rare bacteria, respectively. More importantly, rare bacterial sub-communities were more sensitive to these local environmental changes. Collectively, our findings provide first-hand evidence that abundant and rare microbial sub-communities of summer maize soil respond differently to reduced nitrogen application and groundwater depth change.
Data availability statement
The data presented in the study are deposited in the NCBI repository, accession number PRJNA1021577.
Author contributions
FB, XQ, and PL conceptualized the study. FB carried out field research and gathered information. FB, XQ, and ZD created the techniques. FB and WG wrote the manuscript. DQ was in charge of language editing. All authors read and approved the final manuscript.
Funding
This work was supported by the National Key R&D Program of China (Grant No. 2021YFD1700900), the Agricultural Science and Technology Innovation Program of the Chinese Academy of Agricultural Sciences (Grant No. CAAS-ASTIP), and the Science and Technology Project of Henan Province (Grant No. 212102110233).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1220731/full#supplementary-material
References
Abbasi, M. R., and Sepaskhah, A. R. (2023). Nitrogen leaching and groundwater N contamination risk in saffron/wheat intercropping under different irrigation and soil fertilizers regimes. Sci. Rep. 13:6587. doi: 10.1038/s41598-023-33817-5
Askri, B., Ahmed, A. T., Abichou, T., and Bouhlila, R. (2014). Effects of shallow water table, salinity and frequency of irrigation water on the date palm water use. J. Hydrol. 513, 81–90. doi: 10.1016/j.jhydrol.2014.03.030
Bai, F., Qi, X., Li, P., Du, Z., and Guo, W. (2023). Groundwater depth and nitrogen application amount jointly regulate the water and residual soil nitrate accumulation in agricultural soil profile. Agronomy 13:1163. doi: 10.3390/agronomy13041163
Bai, F., Qi, X., Li, P., Qiao, D., Wang, J., Du, Z., et al. (2020). Depression of groundwater table and reduced nitrogen application jointly regulate the bacterial composition of nir S-type and nir K-type genes in agricultural soil. Water 12:3459. doi: 10.3390/w12123459
Changming, L., Jingjie, Y., and Kendy, E. (2001). Groundwater exploitation and its impact on the environment in the North China plain. Water Int. 26, 265–272. doi: 10.1080/02508060108686913
Chase, A. B., Karaoz, U., Brodie, E. L., Gomez-Lunar, Z., Martiny, A. C., and Martiny, J. B. (2017). Microdiversity of an abundant terrestrial bacterium encompasses extensive variation in ecologically relevant traits. MBio 8, e01809–e01817. doi: 10.1128/mBio.01809-17
Chen, X., Wang, P., Muhammad, T., Xu, Z., and Li, Y. (2020). Subsystem-level groundwater footprint assessment in North China plain – the world’s largest groundwater depression cone. Ecol. Indic. 117:106662. doi: 10.1016/j.ecolind.2020.106662
Cocco, E., Bertora, C., Squartini, A., Delle Vedove, G., Berti, A., Grignani, C., et al. (2018). How shallow water table conditions affect N2O emissions and associated microbial abundances under different nitrogen fertilisations. Agric. Ecosyst. Environ. 261, 1–11. doi: 10.1016/j.agee.2018.03.018
Delgado-Baquerizo, M., Oliverio, A. M., Brewer, T. E., Benavent-González, A., Eldridge, D. J., Bardgett, R. D., et al. (2018). A global atlas of the dominant bacteria found in soil. Science 359, 320–325. doi: 10.1126/science.aap9516
Egidi, E., Delgado-Baquerizo, M., Plett, J. M., Wang, J., Eldridge, D. J., Bardgett, R. D., et al. (2019). A few Ascomycota taxa dominate soil fungal communities worldwide. Nat. Commun. 10:2369. doi: 10.1038/s41467-019-10373-z
Fierer, N., and Jackson, R. B. (2006). The diversity and biogeography of soil bacterial communities. Proc. Natl. Acad. Sci. U. S. A. 103, 626–631. doi: 10.1073/pnas.0507535103
Gan, L., Huang, G., Pei, L., Gan, Y., Liu, C., Yang, M., et al. (2022). Distributions, origins, and health-risk assessment of nitrate in groundwater in typical alluvial-pluvial fans, North China plain. Environ. Sci. Pollut. Res. Int. 29, 17031–17048. doi: 10.1007/s11356-021-17067-4
Hou, S., Ai, C., Zhou, W., Liang, G., and He, P. (2018). Structure and assembly cues for rhizospheric nirK-and nirS-type denitrifier communities in long-term fertilized soils. Soil Biol. Biochem. 119, 32–40. doi: 10.1016/j.soilbio.2018.01.007
Hou, J., Wu, L., Liu, W., Ge, Y., Mu, T., Zhou, T., et al. (2020). Biogeography and diversity patterns of abundant and rare bacterial communities in rice paddy soils across China. Sci. Total Environ. 730:139116. doi: 10.1016/j.scitotenv.2020.139116
Ji, M., Kong, W., Stegen, J., Yue, L., Wang, F., Dong, X., et al. (2020). Distinct assembly mechanisms underlie similar biogeographical patterns of rare and abundant bacteria in Tibetan plateau grassland soils. Environ. Microbiol. 22, 2261–2272. doi: 10.1111/1462-2920.14993
Jiao, S., and Lu, Y. (2020a). Abundant fungi adapt to broader environmental gradients than rare fungi in agricultural fields. Glob. Chang. Biol. 26, 4506–4520. doi: 10.1111/gcb.15130
Jiao, S., and Lu, Y. (2020b). Soil pH and temperature regulate assembly processes of abundant and rare bacterial communities in agricultural ecosystems. Environ. Microbiol. 22, 1052–1065. doi: 10.1111/1462-2920.14815
Jiao, S., Luo, Y., Lu, M., Xiao, X., Lin, Y., Chen, W., et al. (2017). Distinct succession patterns of abundant and rare bacteria in temporal microcosms with pollutants. Environ. Pollut. 225, 497–505. doi: 10.1016/j.envpol.2017.03.015
Jing, H., He, X., Tian, Y., Lancia, M., Cao, G., Crivellari, A., et al. (2023). Comparison and interpretation of data-driven models for simulating site-specific human-impacted groundwater dynamics in the North China plain. J. Hydrol. 616:128751. doi: 10.1016/j.jhydrol.2022.128751
Ju, X.-T., Xing, G.-X., Chen, X.-P., Zhang, S.-L., Zhang, L.-J., Liu, X.-J., et al. (2009). Reducing environmental risk by improving N management in intensive Chinese agricultural systems. Proc. Natl. Acad. Sci. U. S. A. 106, 3041–3046. doi: 10.1073/pnas.0813417106
Koranda, M., Kaiser, C., Fuchslueger, L., Kitzler, B., Sessitsch, A., Zechmeister-Boltenstern, S., et al. (2013). Seasonal variation in functional properties of microbial communities in beech forest soil. Soil Biol. Biochem. 60, 95–104. doi: 10.1016/j.soilbio.2013.01.025
Kurm, V., Van Der Putten, W. H., De Boer, W., Naus-Wiezer, S., and Hol, W. G. (2017). Low abundant soil bacteria can be metabolically versatile and fast growing. Ecology 98, 555–564. doi: 10.1002/ecy.1670
Lasagna, M., De Luca, D. A., and Franchino, E. (2016). Nitrate contamination of groundwater in the western Po plain (Italy): the effects of groundwater and surface water interactions. Environ. Earth Sci. 75, 1–16. doi: 10.1007/s12665-015-5039-6
Levins, R. (1968). Evolution in changing environments: some theoretical explorations. Princeton, NJ: Princeton University Press.
Li, W., Niu, S., Liu, X., and Wang, J. (2019). Short-term response of the soil bacterial community to differing wildfire severity in Pinus tabulaeformis stands. Sci. Rep. 9:1148. doi: 10.1038/s41598-019-38541-7
Liang, Y., Xiao, X., Nuccio, E. E., Yuan, M., Zhang, N., Xue, K., et al. (2020). Differentiation strategies of soil rare and abundant microbial taxa in response to changing climatic regimes. Environ. Microbiol. 22, 1327–1340. doi: 10.1111/1462-2920.14945
Liu, B.-Y., Lin, B.-J., Li, X.-X., Virk, A. L., N'dri Yves, B., Zhao, X., et al. (2021). Appropriate farming practices of summer maize in the North China plain: reducing nitrogen use to promote sustainable agricultural development. Resour. Conserv. Recycl. 175:105889. doi: 10.1016/j.resconrec.2021.105889
Liu, W., Liu, L., Yang, X., Deng, M., Wang, Z., Wang, P., et al. (2021). Long-term nitrogen input alters plant and soil bacterial, but not fungal beta diversity in a semiarid grassland. Glob. Chang. Biol. 27, 3939–3950. doi: 10.1111/gcb.15681
Liu, L., Yang, J., Yu, Z., and Wilkinson, D. M. (2015). The biogeography of abundant and rare bacterioplankton in the lakes and reservoirs of China. ISME J. 9, 2068–2077. doi: 10.1038/ismej.2015.29
Luo, S., Yu, L., Liu, Y., Zhang, Y., Yang, W., Li, Z., et al. (2016). Effects of reduced nitrogen input on productivity and N2O emissions in a sugarcane/soybean intercropping system. Eur. J. Agron. 81, 78–85. doi: 10.1016/j.eja.2016.09.002
Matzeu, A., Secci, R., and Uras, G. (2017). Methodological approach to assessment of groundwater contamination risk in an agricultural area. Agric. Water Manag. 184, 46–58. doi: 10.1016/j.agwat.2017.01.003
Mcgill, B. J., Enquist, B. J., Weiher, E., and Westoby, M. (2006). Rebuilding community ecology from functional traits. Trends Ecol. Evol. 21, 178–185. doi: 10.1016/j.tree.2006.02.002
Mo, Y., Zhang, W., Yang, J., Lin, Y., Yu, Z., and Lin, S. (2018). Biogeographic patterns of abundant and rare bacterioplankton in three subtropical bays resulting from selective and neutral processes. ISME J. 12, 2198–2210. doi: 10.1038/s41396-018-0153-6
Mooshammer, M., Wanek, W., Hämmerle, I., Fuchslueger, L., Hofhansl, F., Knoltsch, A., et al. (2014). Adjustment of microbial nitrogen use efficiency to carbon: nitrogen imbalances regulates soil nitrogen cycling. Nat. Commun. 5, 1–7. doi: 10.1038/ncomms4694
Paladino, O., Massabò, M., and Gandoglia, E. (2020). Assessment of nitrate hazards in Umbria region (Italy) using field datasets: good agriculture practices and farms sustainability. Sustainability 12:9497. doi: 10.3390/su12229497
Pandit, S. N., Kolasa, J., and Cottenie, K. (2009). Contrasts between habitat generalists and specialists: an empirical extension to the basic metacommunity framework. Ecology 90, 2253–2262. doi: 10.1890/08-0851.1
Pedrós-Alió, C. (2012). The rare bacterial biosphere. Annu. Rev. Mar. Sci. 4, 449–466. doi: 10.1146/annurev-marine-120710-100948
Ren, Z., Luo, W., and Zhang, C. (2022). Rare bacterial biosphere is more environmental controlled and deterministically governed than abundant one in sediment of thermokarst lakes across the Qinghai-Tibet plateau. Front. Microbiol. 13:944646. doi: 10.3389/fmicb.2022.944646
Shade, A., Jones, S. E., Caporaso, J. G., Handelsman, J., Knight, R., Fierer, N., et al. (2014). Conditionally rare taxa disproportionately contribute to temporal changes in microbial diversity. MBio 5, e01371–e01314. doi: 10.1128/mBio.01371-14
She, W., Bai, Y., Zhang, Y., Qin, S., Feng, W., Sun, Y., et al. (2018). Resource availability drives responses of soil microbial communities to short-term precipitation and nitrogen addition in a desert Shrubland. Front. Microbiol. 9:186. doi: 10.3389/fmicb.2018.00186
Shi, M., Gong, H., Gao, M., Chen, B., Zhang, S., and Zhou, C. (2020). Recent ground subsidence in the North China plain, China, revealed by sentinel-1A datasets. Remote Sens. 12:3579. doi: 10.3390/rs12213579
Soylu, M. E., Kucharik, C. J., and Loheide Ii, S. P. (2014). Influence of groundwater on plant water use and productivity: development of an integrated ecosystem–variably saturated soil water flow model. Agric. For. Meteorol. 189-190, 198–210. doi: 10.1016/j.agrformet.2014.01.019
Storch, D., Bohdalková, E., and Okie, J. (2018). The more-individuals hypothesis revisited: the role of community abundance in species richness regulation and the productivity–diversity relationship. Ecol. Lett. 21, 920–937. doi: 10.1111/ele.12941
Sun, R., Guo, X., Wang, D., and Chu, H. (2015). Effects of long-term application of chemical and organic fertilizers on the abundance of microbial communities involved in the nitrogen cycle. Appl. Soil Ecol. 95, 171–178. doi: 10.1016/j.apsoil.2015.06.010
Tripathi, B. M., Stegen, J. C., Kim, M., Dong, K., Adams, J. M., and Lee, Y. K. (2018). Soil pH mediates the balance between stochastic and deterministic assembly of bacteria. ISME J. 12, 1072–1083. doi: 10.1038/s41396-018-0082-4
Umaña, M. N., Zhang, C., Cao, M., Lin, L., and Swenson, N. G. (2015). Commonness, rarity, and intraspecific variation in traits and performance in tropical tree seedlings. Ecol. Lett. 18, 1329–1337. doi: 10.1111/ele.12527
Vergine, M., Yuan, X., Zhang, J., Chang, F., Wang, X., Zhang, X., et al. (2023). Effects of nitrogen reduction combined with bio-organic fertilizer on soil bacterial community diversity of red raspberry orchard. PLoS One 18:e0283718. doi: 10.1371/journal.pone.0283718
Wan, W., Gadd, G. M., Yang, Y., Yuan, W., Gu, J., Ye, L., et al. (2021). Environmental adaptation is stronger for abundant rather than rare microorganisms in wetland soils from the Qinghai-Tibet plateau. Mol. Ecol. 30, 2390–2403. doi: 10.1111/mec.15882
Wan, X., Gao, Q., Zhao, J., Feng, J., Van Nostrand, J. D., Yang, Y., et al. (2020). Biogeographic patterns of microbial association networks in paddy soil within eastern China. Soil Biol. Biochem. 142:107696. doi: 10.1016/j.soilbio.2019.107696
Wang, J., Li, M., and Li, J. (2021a). Soil pH and moisture govern the assembly processes of abundant and rare bacterial communities in a dryland montane forest. Environ. Microbiol. Rep. 13, 862–870. doi: 10.1111/1758-2229.13002
Wang, J., Wang, Y., He, N., Ye, Z., Chen, C., Zang, R., et al. (2020). Plant functional traits regulate soil bacterial diversity across temperate deserts. Sci. Total Environ. 715:136976. doi: 10.1016/j.scitotenv.2020.136976
Wang, J., Wang, Y., Li, M., Xu, L., He, N., Yan, P., et al. (2021b). Differential response of abundant and rare bacterial subcommunities to abiotic and biotic gradients across temperate deserts. Sci. Total Environ. 763:142942. doi: 10.1016/j.scitotenv.2020.142942
Webb, C. T., Hoeting, J. A., Ames, G. M., Pyne, M. I., and Leroy Poff, N. (2010). A structured and dynamic framework to advance traits-based theory and prediction in ecology. Ecol. Lett. 13, 267–283. doi: 10.1111/j.1461-0248.2010.01444.x
Wu, T., Gray, A., Liu, G., Kaminski, H., Osi Efa, B., and Aubrey, D. P. (2020). Groundwater depth overrides tree-species effects on the structure of soil microbial communities involved in nitrogen cycling in plantation forests. Forests 11:275. doi: 10.3390/f11030275
Xiao, Y., Liu, H., Senapathi, V., Wang, L., and Li, C. (2023). Editorial: interactions between groundwater and human communities: perspectives on the resources, environments, threats and sustainable development. Front. Environ. Sci. 11:1221837. doi: 10.3389/fenvs.2023.1221837
Xiong, C., He, J. Z., Singh, B. K., Zhu, Y. G., Wang, J. T., Li, P. P., et al. (2021). Rare taxa maintain the stability of crop mycobiomes and ecosystem functions. Environ. Microbiol. 23, 1907–1924. doi: 10.1111/1462-2920.15262
Zhang, D., Cui, R., Fu, B., Yang, Y., Wang, P., Mao, Y., et al. (2020). Shallow groundwater table fluctuations affect bacterial communities and nitrogen functional genes along the soil profile in a vegetable field. Appl. Soil Ecol. 146:103368. doi: 10.1016/j.apsoil.2019.103368
Zhang, Y., Dong, S., Gao, Q., Ganjurjav, H., Wang, X., and Geng, W. (2019). “Rare biosphere” plays important roles in regulating soil available nitrogen and plant biomass in alpine grassland ecosystems under climate changes. Agric. Ecosyst. Environ. 279, 187–193. doi: 10.1016/j.agee.2018.11.025
Zhang, W.-F., Dou, Z.-X., He, P., Ju, X.-T., Powlson, D., Chadwick, D., et al. (2013). New technologies reduce greenhouse gas emissions from nitrogenous fertilizer in China. Proc. Natl. Acad. Sci. 110, 8375–8380. doi: 10.1073/pnas.1210447110
Zhang, K., Shi, Y., Cui, X., Yue, P., Li, K., Liu, X., et al. (2019). Salinity is a key determinant for soil microbial communities in a desert ecosystem. mSystems 4:e00225–18. doi: 10.1128/mSystems.00225-18
Keywords: bacterial composition, abundant and rare bacteria, groundwater depth, reduced nitrogen application, summer maize field
Citation: Bai F, Guo W, Li P, Qiao D, Du Z and Qi X (2023) Different responses of abundant and rare bacterial composition to groundwater depth and reduced nitrogen application in summer maize field. Front. Microbiol. 14:1220731. doi: 10.3389/fmicb.2023.1220731
Edited by:
Raza Waseem, Nanjing Agricultural University, ChinaReviewed by:
Jie Wang, China Agricultural University, ChinaWu Xiong, Nanjing Agricultural University, China
Marika Pellegrini, University of L'Aquila, Italy
Copyright © 2023 Bai, Guo, Li, Qiao, Du and Qi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xuebin Qi, qixuebin@caas.cn
†These authors have contributed equally to this work and share first authorship