- Department of Animal Science and Technology, Northeast Agricultural University, Harbin, China
Changes in land use types can significantly affect soil porperties and microbial community composition in many areas. However, the underlying mechanism of shift in bacterial communities link to soil properties is still unclear. In this study, Illumina high-throughput sequencing was used to analyze the changes of soil bacterial communities in different land use types in a forest-grassland transition zone, North China. There are two different land use types: grassland (G) and cultivated land (CL). Meanwhile, cultivated land includes cultivated of 10 years (CL10) or 20 years (CL20). Compared with G, CL decreased soil pH, SOC and TN, and significantly increased soil EC, P and K, and soil properties varied significantly with different cultivation years. Grassland reclamation increases the diversity of bacterial communities, the relative abundance of Proteobacteria, Gemmatimonadetes and Bacteroidetes increased, while that of Actinobacteria, Acidobacteria, Rokubacteria and Verrucomicrobia decreased. However, the relative abundance of Proteobacteria decreased and the relative abundance of Chloroflexi and Nitrospirae increased with the increase of cultivated land years. Mantel test and RDA analysis showed that TP, AP, SOC and EC were the main factors affecting the diversity of composition of bacterial communities. In conclusion, soil properties and bacterial communities were significantly altered after long-term cultivation. This study provides data support for land use and grassland ecological protection in this region.
Introduction
Soil is one of the most biodiverse habitats on Earth (Fierer and Jackson, 2006; Bender et al., 2016), and soil microorganisms are major components of soil (Roesch et al., 2007), which play an indispensable role in ecosystem function and sustainability (Paul and Clark, 1989; Comer and Perkins, 2020). For example, soil microbes not only drive nutrient and carbon cycles (Kowalchuk and Stephen, 2001; van der Heijden et al., 2008), but also directly and indirectly involved in many other ecosystem services, such as the formation of stable soil aggregates and soil structures to control erosion (Lynch and Bragg, 1985; Rillig and Mummey, 2006). Although there are some functional redundancies in microbial diversity, higher diversity offers a higher possibility for the existence of some alternative species. These alternative species can make up for the loss of function caused by the disappearance of some species under the conditions of temporal and spatial changes, thereby maintaining the function of the ecosystem (Yachi and Loreau, 1999). Meanwhile, a recent study showed that some ecological functions were performed jointly by multiple microorganisms, so functional redundancy may disappear with more functions (Wagg et al., 2019). However, with the decrease of soil microbial diversity, some functions of the ecosystem will be inhibited. For example, reductions in the diversity of methanotrophs (pmoA genes) and denitrifiers (nosZ 369 genes) could have potential negative consequences for climate regulation on Earth by increasing the amount of methane and N2O released to the atmosphere (Trivedi et al., 2019). The diversity of soil microbial community is important in determining soil health (Zhao et al., 2019), which are the foundation of soil function (Emmerling et al., 2002; van der Heijden et al., 2008; Nielsen et al., 2015).
Changes in land use types result from the interaction of human activities and natural factors, which may have important implications for soil ecosystems (Ji et al., 2020). Changes in land use types not only affect soil physical and chemical properties, but also soil biodiversity and soil ecological functions (Lauber et al., 2008; Rasche et al., 2011; Richter et al., 2018). Throughout human history, grasslands have been intensely modified and fragmented by agriculture (Gibson, 2009). According to FAO (2009), global food demand will increase by 70% in the coming decades, which will increase the demand for intensively grown crops (Godfray et al., 2010). Land use intensification from primitive grassland to intensive agricultural systems generally increases soil pH, resulting in decreased soil C content, lower water retention, and poorer soil structure (Finn et al., 2017). For example, Long-term monoculture of cotton or maize in grassland soils significantly decrease soil organic carbon (SOC), macroaggregates, catabolic diversity and microbial biomass (Acosta-Martínez et al., 2010; Kösters et al., 2013). Therefore, it is very important to understand the changes of soil properties and microorganisms after grassland is transformed into cultivated land for monitoring soil quality.
Changes in soil properties (e.g., soil pH, SOC and nutrient content, etc.) can affect the composition and structure of microbial communities (Lauber et al., 2008; Rasche et al., 2011; Richter et al., 2018). Alterations in plant residue and plant root exudation inputs induced by land use changes will affect the soil environmental, and ultimately influence soil microbial community composition and structure (Xu et al., 2021). Agricultural intensification is considered to be one of the biggest threats to global biodiversity (Convention on Biological Diversity, 2010), it has a great impact on soil health and microbial diversity (Degens et al., 2000). For example, Chen et al. (2022) found that conversion of grassland to farmland resulted in a decrease in bacterial biomass. In particular, changes in different tillage regimes can alter plant species and associated soil properties, with lasting effects on soil microbial community composition (Buckley and Schmidt, 2001; Jangid et al., 2008; Yu et al., 2011). For example, biological nitrogen fixation is the most basic physiological process of natural grassland, but the development of intensive agriculture is largely based on industrially produced mineral fertilizers, which leads to simpler soil food webs and fewer functional groups, and reduces overall soil microbial richness (McDaniel et al., 2014; Tsiafouli et al., 2015; Bender et al., 2016). In addition, different kinds of microorganisms have different abilities to deal with various nutrient forms in soil, fertilization will affect their growth competitiveness, thus affecting the diversity, biomass and activity of soil microbial community. Studies have shown that long-term nitrogen application will lead to changes in the overall bacterial community and composition of individual bacterial communities (Rousk et al., 2010), such as Actinobacteria (Jenkins et al., 2009) and Acidobacteria (Zhao et al., 2014). However, the effects of grassland conversion to cultivated land on soil properties and composition and structure of soil microbial community, as well as the main regulators driving the distribution of soil microbial communities across different land use types are still unclear, so understanding the impact of land use change and environmental drivers on soil-microbiome responses is important for soil microbial diversity and soil health.
The forest-grassland transitional zone is a typical ecotone with fertile soil and rich species. Plant community in the transitional zone have a high degree of biodiversity and competitiveness, which are extremely sensitive to environmental changes and have poor ability to resist interference (Anchorena and Cingolani, 2004). Therefore, it is of great significance to understand the impact of human disturbance on the forest and grassland transition zone for the protection and management it. In this study, the mountainous meadow is located in the forest and grassland transition zone between the Great Khingan Mountains and Hulun Buir steppe, which belongs to the semi-humid and semi-arid transition zone with typical transition characteristics of vegetation, soil, land use, climate and other environmental factors, and is a relatively sensitive ecological fragile zone. In the past, many tracts of grassland have been cultivated for years as agricultural land. Therefore, it is necessary to study the impact of long-term cultivation on soil microbial properties. The objectives of this study were to (1) study the effects of long-term cultivation (10-year cultivated land and 20-year cultivated land converted from grassland) on soil properties and the structure, and diversity and composition of soil bacterial communities and (2) clarify the key factors controlling bacterial response to long-term cultivation in Inner Mongolia meadow steppe.
Materials and methods
Study site
This study area is located at the western foot of the Great Hinggan Mountains, a forest and grassland transition zone between the Great Khingan Mountains and Hulun Buir Grassland, belonging to Yakeshi City, Inner Mongolia Autonomous Region (49°19´ N, 119°56′ E; 666-680 m asl), which has a cold temperate continental monsoon climate. The annual average precipitation is 350–400 mm, and the precipitation period is mostly concentrated in June to September. The average annual temperature was −1.5°C, the highest and lowest temperatures were 36.2°C and − 48.5°C, respectively. The annual accumulated temperature ≥ 10°C is 1700–2,300°C, the annual average sunshine is 2,511–2,663 h, the frost-free period is about 76–95 days, and the soil type is mainly chernozum soil. The grassland type of the study area is mountain meadow, with rich plant species, high plant community coverage and high grassland productivity. The main grassland plants are Sanguisorba officinalis, Vicia sepium L., Lathyrus sativus, Carex spp., etc.
Experiment design
Three plots in each land use types were selected from the mountain meadows on the outskirts of Yakeshi, Inner Mongolia. The plots were grassland, 10-year cultivated land, and 20-year cultivated land, denoted as G, CL10, and CL20, respectively. The grassland has been mowed for many years without grazing or human improvement. The cultivated lands are continuous tillage, no abandoned in the middle, planting crops mainly for wheat, potatoes, rape, etc., with long-term fertilization. Soil at different depths of 0–10 cm and 10–20 cm (denoted as D10 and D20) was collected.
Soil sampling and analysis
The samples were collected in July 2021. Three transect lines were set inside G, CL10 and CL20, and three soil sampling points were randomly set on each line. A soil corer (with a diameter of 4 cm) was used to sample the soil with five-point sampling method. Soil samples were taken at different depths of 0–10 cm (D10) and 10–20 cm (D20). Three samples from the same line and the same depth were thoroughly mixed to remove root and soil intrusions, quickly packed into sterile sealed bags, refrigerated in an ice box and quickly transported to the laboratory. The samples were divided into two parts. One part was air-dried, ground and sieved to measure soil physical and chemical properties; other part was stored at −80°C for soil microbial analysis.
Soil nutrient determination
Soil pH and electrical conductivity (EC; an important indicator of salinity) were determined with soil (air-dried)/water (1:5, w/v) suspensions (Bao, 2008), and soil organic matter (SOM) was assessed by K2CrO7 volumetric method (External heating method; Sparks et al., 1996). Total nitrogen (TN) and alkali-hydrolyzed nitrogen (AN) were measured by Kjeldahl method and alkali-hydrolyzed diffusion method (Lu, 1999). Phosphorus (P) was extracted by HClO4-H2SO4 method, then added the molybdenum-antimony anti-reagent, colorimetric was performed at 880 nm or 700 nm wavelength, and calculated the total phosphorus content (Bao, 2008). Available phosphorus (AP) was extracted with 0.5 M NaHCO3 (pH 8.5), and was determined using molybdenum blue method (UV-752 Shanghai, China; Bao, 2008).
Soil DNA extraction, Illumina sequencing and data analysis
According to the manufacturer, the total DNA was extracted from 5 g of soil (−80°C) using the TGuide S96 magnetic bead method soil genomic DNA extraction kit [Tiangen Biochemical Technology (Beijing) Co., Ltd., model: DP812]. The total DNA concentration and quality was detected using a microplate reader (manufacturer: GeneCompang Limited, model synergy HTX).
According to the construction of bacterial 16S rRNA amplification library, primers 338F (5’-ACTCCTACGGGAGGCAGCA-3’) and 806R (5’- GGACTACHVGGGTWTCTAAT-3’) were used to amplifies hypervariable regions V3 and V4. PCR amplification was carried out in 0.3 μl Vn F (10 μm; Suzhou Hongxun Biotechnology Co., Ltd.), 0.3 μl Vn R (10 μm; Suzhou Hongxun Biotechnology Co., Ltd.), 5 μl KOD FX Neo Buffer (Beijing Bailing Gram Biotechnology LLC Co., Ltd.), 2 μl dNTP (2 mM each), 0.2 μl KOD FX Neo and 50 ng of template DNA in a 10 μl mixture. The amplification process was pre-denaturation at 95°C for 5 min, 25 cycles of 95°C for 30 s, 50°C annealing for 30 s, 72°C for 40 s, and 72°C for 7 min of extension. Each sample was amplified in triplicate. and PCR products was further purified by VAHTSTM DNA Clean Beads magnetic beads. After purification, agarose with a concentration of 1.8% was used for electrophoresis detection, and the concentration (Qubit) was quantified, and the qualified samples were mixed. The mixed samples were recovered by using the Monarch DNA gel recovery kit, and the quality was checked by the Qsep-400 method. Finally, the Illumina novaseq6000 (novaseq6000, illumina) platform was used for sequencing and microbial community analysis.
Using FLASH (version 1.2.11) software (Magoc and Salzberg, 2011), according to the minimum overlap length of 10 bp and the maximum allowable mismatch ratio in the overlap region of 0.2, the reads of each sample were spliced. In order to obtain high-quality Tag data, we use Trimmomatic (version 0.33) to filter tags whose average quality value is less than 20 and whose length is less than 75% of the Tag length (Bolger et al., 2014). Then split the query sequence into chunks without overlap, and use UCHIME (version 8.1) to compare the database to detect and delete chimeras (Edgar et al., 2011). Sequences were clustered (OTU) at the 97% similarity level, and OTUs were filtered with a threshold of 0.005% of all sequences sequenced using USEARCH (version 10.0; Edgar, 2013). Then used RDP software (Version 2.2) to compare with Silva database, and the bacteria were classified at a confidence interval of 80% (Wang et al., 2007; Quast et al., 2013).
Statistical analysis
All statistical analyses were performed using R (4.1.2) software. Multivariate analysis of variance was used to analyze the effects of land use type, soil depth and their interactions on soil physical and chemical properties and bacterial communities. In order to calculate the α diversity of each sample, we used the “picante” package and BMKCloud to calculate Chao1, Simpson diversity, and Shannon’s evenness. Two-way ANOVA and TukeyHSD test in the “agricolae” package was used to analyze the significant differences between groups in soil physicochemical properties and soil bacterial α diversity. Spearman rank correlation coefficient was used to calculate the correlation between soil physicochemical properties and soil bacterial α diversity. Based on the “vegdist” package, we calculated the Bray-Curtis distance matrix between OTUs of different samples, analyze the similarity of 999 permutations (ANOSIM), and PCoA visualization was performed on soil bacteria of different sites and depths. The linear discriminant analysis (LDA) effect size (LEfSe) method was used to identify features with significant differences in bacteria between the three sites and different depths by biomarkers, and evaluate the effect size of each feature with a threshold of 4.0 and a significant α of 0.5. To investigate the relationship between soil properties and bacterial communities, we performed RDA analysis of soil properties and bacterial communities using “vegdist” packs and Mantel tests using the “vegan” and “ggcor_master” packages to infer the potential relationship between microbial community composition and measured soil properties. In order to further study, the effect of soil properties on the differences in bacterial community abundance and phylum, order, and genus levels in different sites, Spearman rank correlation coefficient was used.
Results
Effects of land use types on soil properties
PCoA results showed that soil properties of G, CL10, CL20 and different soil layers were significantly separated, indicating that soil properties of different land use types and different soil layers were significantly different (R2 = 0.86299, p = 0.001; Figure 1). Different land use types had significant effects on soil properties (p < 0.001), soil depth had effects on TP, AP, TK, AK and pH, and years of cultivated land had effects on AN, TP, AP, TK, AK, SOC and EC (Table 1). In the D10 and D20 soil layers, the TN, SOC, pH and C:N of G were the highest (p < 0.05). Compared with G and CL10, the contents of AN and TK in each layer of CL20 were the highest (p < 0.05). The AN, TP, AP, AK and EC of CL10 were the highest (p < 0.05). With the increase of cultivation years, the contents of AN and TK in each soil layer increased (Table 2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1001781/fmicb-13-1001781-HTML/image_m/fmicb-13-1001781-g001.jpg)
Figure 1. Principal co-ordinates analysis (PCoA) of soil properties among different land use types and depths.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1001781/fmicb-13-1001781-HTML/image_m/fmicb-13-1001781-t001.jpg)
Table 1. Anova analysis of the effects of land use type (LUT), soil depth (Depth) and cultivated land years (Year) on soil properties and bacterial diversity and relative abundance.
Effects of land use types on soil microbial community diversity and composition
A total of 1,440,225 bacterial sequences were obtained after Illumina high-throughput sequencing of 18 samples from 6 groups of samples, and 1,397,395 effective reads were obtained after filtration and removal of low-quality sequences. The effective sequence length was mainly distributed in the range of 400–450 bp, and the average ratio of effective data to original data was 97.03%, with the maximum value 98.7% and the minimum 94.97%. The Q20 of all samples ranged from 98.84 to 99.03%, with an average of 98.95%. There were 31,094 bacterial OTUs in total, and the OTUs of each sample ranged from 1,519 to 2,442, with an average of 1765, and the homology threshold was 97%. The Good’s Coverage index of all samples was higher than 99.18%, indicating that the sequencing quality was good.
Land use types had significant effects including bacteria Chao1 richness and Shannon’s evenness (p < 0.05), but soil depth and years of cultivated land had no effects including Chao1 richness, Simpson diversity, Shannon’s evenness and Good’s coverage (Table 1). Compared with G, the Chao1 richness and Shannon’s evenness in CL were significantly increased, but there were no significant differences among different years of cultivated land. Overall, grassland reclamation increased bacterial richness and evenness, which were not affected by years of cultivated land (Figure 2; Supplementary Table S1). Spearman rank correlation analysis showed that Chao1 richness was positively correlated with AN, EC and TK, and negatively correlated with pH, TN and C:N. Bacterial Shannon’s evenness was positively correlated with AN and TK, and negatively correlated with pH and TN (Figure 2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1001781/fmicb-13-1001781-HTML/image_m/fmicb-13-1001781-g002.jpg)
Figure 2. Effects of different land use types, depths, cultivated land years and soil properties on soil bacterial community diversity[Values are the means ± SEs (n = 3). Different letters represent significant differences between the means (p < 0.05).]
Land use type had a significant effect on bacterial β diversity (p < 0.001), while soil depth and years of cultivated land had no effect on bacterial β diversity (Table 1). Based on the Bray-Curtis distance matrix, we performed PCoA visualization of soil bacteria at different land use types and depths. PCo1 and PCo2 explained 46.51 and 31.14% of the variation in bacterial community composition, respectively (Figure 3). The results of Permanova analysis showed that there were significant differences in bacterial community composition among different land use types (R2 = 0.55142, p = 0.005). Along the directions of PCo1 and PCo2, the bacterial communities in G were far away from those in CL, and the separation of soil bacterial communities between CL10 and CL20 was not significant, indicating that the soil bacterial communities in G and CL were significantly different, and the composition of soil bacterial communities in CL10 and CL20 was similar.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1001781/fmicb-13-1001781-HTML/image_m/fmicb-13-1001781-g003.jpg)
Figure 3. Principal co-ordinates analysis (PCoA) of bacterial communities under different land use types and depths.
At the phylum level, land use type had a significant effect on the relative abundance of dominant phylum bacteria, while cultivated land year had a significant effect on the relative abundance of Proteobacteria, Chloroflexi, and Nitrospirae, while soil depth had no effect on the relative abundance of dominant phylum bacteria (Table 1). Soil bacterial phyla in CL10 and CL20 had significantly higher relative abundances of Proteobacteria, Gemmatimonadetes and Bacteroidetes, and significantly lower relative abundance of Actinobacteria, Acidobacteria, Verrucomicrobia and Rokubacteria compared with G. However, the relative abundance of Proteobacteria decreased and the relative abundance of Chloroflexi and Nitrospirae increased with the increase of cultivated land years. Acidobacteria was the most abundant dominant phylum in G, and Proteobacteria was the dominant bacterial phylum in CL10 and CL20 (Figure 4; Supplementary Table S2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1001781/fmicb-13-1001781-HTML/image_m/fmicb-13-1001781-g004.jpg)
Figure 4. Relative abundance of dominant species in phylum, order and genus level under different land use types and depths.
At the order level, except for Phycisphaerales and Rhizobiales, different land use types had significant effects on the relative abundance of other dominant bacteria. Except for the years of cultivated land had significant effects on the relative abundance of uncultured_bacterium_c_Subgroup_6, soil depth and years of cultivated land had no significant effect on the relative abundance of other dominant bacteria (Table 1). The relative abundance of the top 10 dominant bacteria in all soil samples was greater than 48%. The dominant bacteria with the highest abundance in G-D10, G-D20, CL10-D10, CL10-D20, CL10-D10, CL10-D10, CL10-D20, CL20-D10 and CL20-D20 samples were uncultured_bacterium_c_Subgroup_6 (12.61%), uncultured_bacterium_c_Subgroup_6 (13.31%), Sphingomonadales (9.81%), Sphingomonadales (10.95%), uncultured_bacterium_c_Subgroup_6 (9.78%) and uncultured_bacterium_c_Subgroup_6 (10.38%), respectively. The relative abundance of dominant orders was significantly different between CL and G. The relative abundance of Rokubacteriales, Chthoniobacterales and uncultured_bacterium_c_Subgroup_6 in G was higher than that in CL, however, the relative abundances of Xanthomonadales, Sphingomonadales, Gemmatimonadales and Betaproteobacteriales were lower (Figure 4; Supplementary Table S3).
At the genus level, except uncultured_Bacterium_c_KD4-96 and Gemmatimonas, different land use types had significant effects on the relative abundance of other dominant bacteria. The years of cultivated land had significant effect on the relative abundance of uncultured_Bacterium_C_KD4-96 and uncultured_Bacterium_c_Subgroup_6, while the soil depth had no significant effect on the relative abundance of dominant bacteria (Table 1). The relative abundances of the top 10 dominant genera were greater than 28% in all soil samples. The dominant bacteria with the highest abundance in G-D10, G-D20, CL10-D10, CL10-D20, CL20-D10 and CL20-D20 samples were uncultured_bacterium_c_Subgroup_6 (12.61%) and Uncultured_Bacterium_c_Subgroup_6 (13.31%), uncultured_bacterium_c_Subgroup_6 (9.34%), Sphingomonadales (7.87%), uncultured_bacterium_c_Subgroup (9.78%) and uncultured_bacterium_c_Subgroup_6 (10.38%). Compared with G, the relative abundance of uncultured_bacterium_o_Rokubacteriales, Candidatus_Udaeobacter and RB41 in CL soil was significantly decreased. However, the relative abundances of uncultured_Bacterium_f_Burkholderiaceae, uncultured_bacterium_f_Gemmatimonadaceae and Sphingomonas increased significantly. With the increase of cultivated land years, the relative abundance of uncultured_Bacterium_c_KD4-96 and uncultured_Bacterium_c_Subgroup_6 increased significantly (Figure 4; Supplementary Table S4).
LEfSe was used to identify species with significant differences in abundance in each soil. The results showed that there were 10, 19, 7, 10 and 2 bacterial populations of G-10,G-20,CL10-10, and CL20-20 samples with significant differences, respectively. The species with significant differences in soil abundance were the most abundant in G, indicating that the structure and diversity of bacteria in grassland soil were significantly different from that in CL. It is worth noting that with the increase of cultivated years, the species with significant differences in CL20 decreased significantly, mainly Chloroflexi and Firmicutes (Figure 5).
Relationships between soil properties and soil bacterial community
Mantel test was used to analyze the relationship between bacterial community composition and soil properties. The results showed that TP, AP, SOC and EC were the main factors affecting the composition of soil bacterial communities (Table 3). This result was further supported by RDA analysis, with the first RDA1 axis and the RDA1 axis explaining 60.08 and 17.58% of the variation in overall bacterial community composition, respectively. RDA results showed that TN was a soil factor that significantly affected bacterial community composition among all soil traits (p = 0.009; Figure 6; Table 4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1001781/fmicb-13-1001781-HTML/image_m/fmicb-13-1001781-t003.jpg)
Table 3. Mantel test of the effects of soil properties on bacteria community composition under different land use types and depths.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1001781/fmicb-13-1001781-HTML/image_m/fmicb-13-1001781-g006.jpg)
Figure 6. Redundancy analysis (RDA) of bacterial community composition and soil properties in different land use types and depths.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1001781/fmicb-13-1001781-HTML/image_m/fmicb-13-1001781-t004.jpg)
Table 4. Redundancy analysis (RDA) of the influence of soil properties on bacterial community and bacterial dominant phylum, order and genus levels in different land use types and depths.
RDA was used to analyze the relationship between the relative abundance of dominant bacteria at phylum, order and genus level and soil properties. At the phylum level, the first two axes of PDA accounted for 92.01 and 4.13% of the relative bacterial abundance, respectively. At the order level, the first two axes of RDA accounted for 80.75 and 8.35% of the relative bacterial abundance, respectively. At the genus level, the first two axes of RDA accounted for 85.25 and 8.91% of the relative abundance of bacteria, respectively (Supplementary Figure S1). TN was the main factor affecting the relative abundance of dominant bacteria at phylum, order and genus levels (p = 0.005, p = 0.022, p = 0.001; Table 4). Spearman correlation analysis showed that the dominant bacterial phylum Proteobacteria were negatively correlated with TN and pH. Actinobacteria is negatively correlated with AK, EC, AP and TP, while positively correlated with pH. The dominant bacterial order Rokubacteriales were negatively correlated with AK, EC, AP and TP, but positively correlated with pH. Betaproteobacteriales were negatively correlated with EC, AP and TP, but positively correlated with pH. The dominant bacterial genus uncultured_bacterium_c_Subgroup_6 was negatively correlated with EC, AP and TP, but positively correlated with TN, pH and SOC. Sphingomonadales was negatively correlated with AP and TP (Figure 7).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1001781/fmicb-13-1001781-HTML/image_m/fmicb-13-1001781-g007.jpg)
Figure 7. Spearman correlation analysis was conducted to analyze the correlation between the relative abundance of dominant bacterial groups and soil properties.
Discussion
Effects of land use types on soil properties
Soil properties have high spatial variability among land use types (Jangid et al., 2008; Lauber et al., 2008; Zhang et al., 2020). In our study, we found that the soil properties of G and CL were significantly different. Compared with CL, G had higher SOC content. The reasons for this are, on the one hand, the removal of biomass at harvest time in farmland (Mendez-Millan et al., 2014), and on the other hand, CL is more affected by human activities than grassland, which makes soil loose and plant residues in it easier to decompose (Trivedi et al., 2016). In addition, the conversion from G to CL, due to the extensive use of fertilizers (e.g., urea) and the different uptake of ions by vegetation, crops have a higher uptake of cations, releasing protons (H+) into the soil solution and reducing the pH of the soil (Barak et al., 1997; Kunhikrishnan et al., 2016; Flores-Rentería et al., 2020). This acidification process accelerates soil SOC dissolution. However, regardless of CL or G, the SOC content in this study area is relatively higher than that in other areas. It is well known that low temperatures can hinder rapid mineralization of SOC due to the temperature sensitivity of microbial communities (Feng et al., 2017). In the Yakeshi area of Inner Mongolia, the latitude is high, the average annual temperature is only −1.5°C, and the soil temperature is low for a long time. It is decomposed by microorganisms and becomes nutrients that can be directly absorbed and utilized by plant roots. Therefore, the content of SOC in the soil gradually accumulates. We found that compared with G, the contents of N, P and K in CL were higher, and there were no significant differences in different cultivated years (p > 0.05). Our results were consistent with previous studies (Ahn et al., 2012; Zhang et al., 2014), which may be attributed to the fact that the heavy application of fertilizers in agricultural management measures (such as long-term use of nitrogen fertilizer) offset the loss of nutrients (Ramirez et al., 2010; Rousk et al., 2010; Smith et al., 2015). In addition, a recent meta-analysis also showed that adding one or more crops in field rotation can effectively increase soil N content (McDaniel et al., 2014). With the increase of cultivated years, soil N, P and K contents increased, but pH decreased. The results were consistent with previous studies (Zhang et al., 2014; Zhou et al., 2015), which may be attributed to the accumulation of nutrients caused by long-term application of chemical fertilizers, which acidified the soil (Guo et al., 2010). In addition, it was found that the increase of soil depth and the decrease of soil pH indicated that the plant roots had a higher uptake of cations and released protons (H+) into the soil solution. Kugbe and Issahaku (2015) and Fetene and Amera (2018) also found a strong positive relationship between pH and basic cations. The decrease of basic cations and the increase of H+ decreased the pH of the deep soil. The contents of N, P and K in different soil depths were significantly different, which may be due to the reduction of surface litter accumulation by mowing of grassland and harvesting of cultivated plants.
Effects of land use types on soil microbial diversity and structure
Changes in land use types can directly alter soil physicochemical properties and structure, thereby affecting the diversity of soil microbial communities (Liu et al., 2021). We initially hypothesized that changes in soil properties after grassland reclamation would affect soil microbial community characteristics. Through PCoA analysis, we found that the characteristics of soil bacterial communities were significantly different among the plots. Previous studies have shown that agriculture may reduce soil microbial diversity (Zhang et al., 2020), in contrast, our study found that the richness and diversity of soil bacterial communities increased after grassland reclamation. It is well known that many factors affect soil microbial diversity, such as soil properties (Fierer and Jackson, 2006), land management systems (Carbonetto et al., 2014), and soil vegetation types (Stephan et al., 2000). Soil nutrient limitations in natural grasslands may affect the growth of microorganisms, but the application of fertilizers in agricultural fields creates the possibility for the growth of new species to some extent. At the same time, the regular application of chemical fertilizers in CL also exacerbates the temporal and spatial changes of soil physical and chemical properties, which may also create a certain niche for the growth of microorganisms (Lee-Cruz et al., 2013; Figuerola et al., 2015; Trivedi et al., 2016). In addition, frequent disturbance of the soil environment by certain agricultural practices (such as tillage) in the field actually helps to promote bacterial community diversity and system resilience (Connell, 1978).
Changes in land use types can not only affect microbial diversity, but also effectively shape the structure of soil biota by altering the overall abundance of taxa (Wang et al., 2020). This study found that the abundances of dominant phyla, orders, and genera were significantly different between G and CL, but there was little difference between different cultivated land years. This may be due to the long-term fertilization of cultivated land dominant genera to adapt to the environment. However, both G and CL dominant phyla are Proteobacteria and Acidobacteria, which is consistent with the findings of Zhang et al. (2019). The abundance of Proteobacteria in G was significantly lower than that in CL. Studies have shown that Proteobacteria have various nutritional methods, not only can they maintain their survival through heterotrophic methods, but many special genera can also carry out autotrophic life, such as using photosynthetic pigments to obtain energy from sunlight for metabolism, and energy acquisition methods are diverse. Fertilization allows them to adapt to complex environmental conditions (Mukhopadhya et al., 2012), so soil matrix changes caused by regular fertilization of cultivated land increase the abundance of Proteobacteria in the soil (Goldfarb et al., 2011; Li et al., 2014). Compared with nutrient-rich CL, the abundance of Acidobacteria in G was higher. Acidobacteria is a kind of oligotrophic bacteria, which has unique characteristics of resistance to impoverishment and low temperature, so that it can maintain high metabolic activity under low temperature conditions, the rapid utilization and absorption of OM in the soil. In cold areas, they can rapidly decompose OM by secreting a large number of enzymes related to carbon and nitrogen metabolism (such as β-glucosidase and fibrinose hydrolase, etc.) to obtain energy for their own growth, which is an important dominant bacteria phyla prevalent in cold soil (Yergeau et al., 2010; Peng et al., 2017). In addition, in our study, we found that uncultured_bacterium_c_Subgroup_6 was dominant at the G order and genus level, which is consistent with previous findings (Nacke et al., 2011; Lin et al., 2019). Naether et al. (2012) performed 16S rRNA sequencing on 27 grasslands and 30 woodlands in Germany, indicating that the grassland soils are dominated by the Subgroup_6 subgroup of Acidobacteria. LefSe analysis showed that there were significant differences in the dominant phyla of different land use types, and the species with significant differences in G were significantly higher than those in CL, and the difference became more significant with the increase of cultivated land years. This indicates that the diversity and structure of soil bacteria have been significantly changed after grassland reclamation.
Effects of soil nutrients on soil microbial community
The richness and diversity of different soil microorganisms are related to specific soil properties (George et al., 2019). Numerous studies have shown that soil pH is related to bacterial community composition (Rousk et al., 2010; Zhao et al., 2014). Our correlation showed that the bacterial richness index was positively correlated with AN and TK, and negatively correlated with pH, TN and C:N. The evenness index was positively correlated with AN and TK, and negatively correlated with pH. Since most bacterial groups exhibit relatively narrow growth tolerance, soil pH may directly influence changes in bacterial communities (Rousk et al., 2010). After grassland reclamation, the contents of N, P and K in CL soil were significantly higher than those in G. Soil properties change is a powerful driver of soil microbial diversity (Fierer, 2017; Tripathi et al., 2018). Therefore, the α-diversity of soil bacteria in CL was significantly higher than that in G. In addition, Mantel test showed that the bacterial community was affected by TN content. Studies showed that soil bacteria were affected by TN content and played a key role in the change of community composition and structure (Thomson et al., 2015; Zhou et al., 2015). Soil TN content can directly affect soil microorganisms, or indirectly affect soil microorganisms by changing soil carbon availability, C:N ratio and soil pH (Treseder, 2008; Serna-Chavez et al., 2013). Therefore, TN content may be second only to pH in influencing soil microbial community composition.
In addition, RDA analysis and Spearman correlation analysis showed that the dominant phylum, order and genus were related to the contents of AK, EC, AP, TP and SOC in soil. Studies have shown that soil P content can affect the relative abundance of different bacteria (Chhabra et al., 2012; Tan et al., 2012). For example, P application reduces soil pH, resulting in a decrease in the abundance of acidic bacteria due to a positive correlation between the abundance of acidic bacteria and soil pH (Jones et al., 2009). At the same time, the increase in soil P content resulted in an increase in the relative abundance of proteobacteria and oligotrophic bacteria in soil (Fierer et al., 2012). In this study, we also found that the relative abundance of some dominant bacteria phyla, orders and genera was negatively correlated with soil salinity (EC), which was consistent with previous studies (van Horn et al., 2014; Xie et al., 2017). With the increase of soil salinity, only a limited number of bacterial groups may be able to withstand the considerable pressure placed on microbial cells in high-salinity soils, leading to changes in the relative abundance of different bacteria and community diversity (Rath et al., 2018). In addition, we also found that the dominant genus uncultured_bacterium_c_Subgroup_6 was positively correlated with SOC, TN, C:N and pH, and negatively correlated with AP and TP. Previous studies have shown that c_Subgroup_6 is relatively abundant in soils with high nutrient levels, and there is a linear positive correlation between soil pH and soil carbon content (Fierer et al., 2007; Jones et al., 2009; Navarrete et al., 2010; Liu et al., 2016). However, Navarrete et al. (2015) and Zhang et al. (2021) showed that c_Subgroup_6 was negatively correlated with soil pH. C_Subgroup_6 correlated positively or negatively with soil pH, suggesting that the same subgroup of Acidobacteria may behave differently in different types of soil (Naether et al., 2012). While SOC is a key resource necessary for most terrestrial microbial communities, SOC levels may affect the abundance of soil bacterial facies and the diversity of bacterial communities (Cai et al., 2016). In summary, our study further demonstrates that different bacteria have different sensitivity to different soil properties, and that soil (microenvironment) heterogeneity is a factor that maintains the diversity of soil microbial communities.
Conclusion
In conclusion, our study found that the reclamation of mountainous meadow in Inner Mongolia changed soil properties, and with the extension of tillage time, soil nutrients were lost and fertility decreased. Changes in land use types further lead to differences in microbial communities, and there are significant differences in bacterial diversity and community structure between different years. Soil depth had a significant effect on soil properties, but had no effect on soil bacterial communities. Grassland reclamation increases the diversity of bacterial communities, the relative abundance of Proteobacteria, Gemmatimonadetes and Bacteroidetes increased, while that of Actinobacteria, Acidobacteria, Rokubacteria and Verrucomicrobia decreased. And long-term cultivation increased the relative abundance of Chloroflexi and Nitrospirae, while decreased that of Proteobacteria. In addition, we found that the bacterial community is mainly affected by EC, SOC, AP and TP, and different bacteria are affected by different soil properties, which are the factors that maintain the diversity of soil microbial community. Our results provide important information for further understanding soil microbial community changes caused by soil microenvironment changes, which is of great significance for conservation of soil diversity.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: NCBI, PRJNA866211.
Author contributions
ZyB, LZ, and MW conceived and designed this study. ZB drafted the original manuscript. ZjB provided very constructive suggestions for revisions. ZB, LZ, MW, and AJ contributed to the sampling and data analysis. All authors contributed to the article and approved the submitted version.
Funding
This project was supported by the Special Fund Project for Returnees Studying Abroad of Heilongjiang Province (2021294).
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.2022.1001781/full#supplementary-material
SUPPLEMENTARY FIGURE S1 | Redundancy analysis (RDA) of bacterial composition and soil properties at dominant phylum, order and genus levels in soils of different land use types and dept.
References
Acosta-Martínez, V., Bell, C. W., Morris, B. E., Zak, J. C., and Allen, V. G. (2010). Long-term soil microbial community and enzyme activity responses to an integrated cropping-livestock system in a semi-arid region. Agric. Ecosyst. Environ. 137, 231–240. doi: 10.1016/j.agee.2010.02.008
Ahn, J., Song, J., Kim, B., Kim, M., Joa, J. H., and Weon, H. (2012). Characterization of the bacterial and archaeal communities in rice field soils subjected to long-term fertilization practices. J. Microbiol. 50, 754–765. doi: 10.1007/s12275-012-2409-6
Anchorena, J., and Cingolani, A. M. (2004). Identifying habitat types in a disturbed area of the forest-steppe ecotone of Patagonia. Plant Ecol. 158, 97–112. doi: 10.1023/A:1014768822737
Barak, P., Jobe, B., Krueger, A. R., Peterson, L. A., and Laird, D. A. (1997). Effects of long-term soil acidification due to nitrogen fertilizer inputs in Wisconsin. Plant Soil 197, 61–69. doi: 10.1023/a:1004297607070
Bender, S. F., Wagg, C., and van der Heijden, M. G. (2016). An underground revolution: biodiversity and soil ecological engineering for agricultural sustainability. Trends Ecol. Evol. 31, 440–452. doi: 10.1016/j.tree.2016.02.016
Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170
Buckley, D. H., and Schmidt, T. M. (2001). The structure of microbial communities in soil and the lasting impact of cultivation. Microb. Ecol. 42, 11–21. doi: 10.1007/s002480000108
Cai, W., Li, Y., Wang, P., Niu, L., Zhang, W., and Wang, C. (2016). Revealing the relationship between microbial community structure in natural biofilms and the pollution level in urban rivers: a case study in the Qinhuai River basin, Yangtze River Delta. Water Sci. Technol. 74, 1163–1176. doi: 10.2166/wst.2016.224
Carbonetto, B., Rascovan, N., Alvarez, R., Mentaberry, A. N., and Vázquez, M. P. (2014). Structure, composition and metagenomic profile of soil microbiomes associated to agricultural land use and tillage Systems in Argentine Pampas. PLoS One 9:e99949. doi: 10.1371/journal.pone.0099949
Chen, Q., Yang, F., and Cheng, X. (2022). Effects of land use change type on soil microbial attributes and their controls: data synthesis. Ecol. Indic. 138:108852. doi: 10.1016/j.ecolind.2022.108852
Chhabra, S., Brazil, D., Morrissey, J. P., Burke, J. I., O'Gara, F., and Dowling, D. N. (2012). Fertilization management affects the alkaline phosphatase bacterial community in barley rhizosphere soil. Biol. Fertil. Soils 49, 31–39. doi: 10.1007/s00374-012-0693-2
Comer, J. A., and Perkins, L. B. (2020). Resistance of the soil microbial community to land-surface disturbances of high-intensity winter grazing and wildfire. J. Environ. Manag. 279:111596. doi: 10.1016/j.jenvman.2020.111596
Connell, J. H. (1978). Diversity in tropical rain forests and coral reefs. Science 199, 1302–1310. doi: 10.1126/science.199.4335.1302
Convention on Biological Diversity (2010). Global biodiversity outlook 3. Canada: UNEP/Earthprint: Montreál.
Degens, B. P., Schipper, L. A., Sparling, G. P., and Vojvodic-Vukovic, M. (2000). Decreases in organic C reserves in soils can reduce the catabolic diversity of soil microbial communities. Soil Biol.Biochem. 32, 189–196. doi: 10.1016/S0038-0717(99)00141-8
Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/nmeth.2604
Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200. doi: 10.1093/bioinformatics/btr381
Emmerling, C., Schloter, M., Hartmann, A., and Kandeler, E. (2002). Functional diversity of soil organisms: a review of recent research activities in Germany. J. Plant Nutr. Soil Sci. 165, 408–420. doi: 10.1002/1522-2624(200208)165:4<408::AID-JPLN408>3.0.CO;2-3
Feng, W., Liang, J., Hale, L., Jung, C. G., Chen, J., Zhou, J., et al. (2017). Enhanced decomposition of stable soil organic carbon and microbial catabolic potentials by long-term field warming. Glob. Chang. Biol. 23, 4765–4776. doi: 10.1111/gcb.13755
Fetene, E. M., and Amera, M. (2018). The effects of land use types and soil depth on soil properties of Agedit watershed, Northwest Ethiopia. Ethiopian J. Sci. Technol. 11, 39–56. doi: 10.4314/ejst.v11i1.4
Fierer, N. (2017). Embracing the unknown: disentangling the complexities of the soil microbiome. Nat. Rev. Microbiol. 15, 579–590. doi: 10.1038/nrmicro.2017.87
Fierer, N., Bradford, M. A., and Jackson, R. B. (2007). Toward an ecological classification of soil bacteria. Ecology 88, 1354–1364. doi: 10.1890/05-1839
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
Fierer, N., Lauber, C. L., Ramirez, K. S., Zaneveld, J. R., Bradford, M. A., and Knight, R. (2012). Comparative metagenomic, phylogenetic and physiological analyses of soil microbial communities across nitrogen gradients. ISME J. 6, 1007–1017. doi: 10.1038/ismej.2011.159
Figuerola, E. L., Guerrero, L. D., Türkowsky, D., Wall, L. G., and Erijman, L. (2015). Crop monoculture rather than agriculture reduces the spatial turnover of soil bacterial communities at a regional scale. Environ. Microbiol. 17, 678–688. doi: 10.1111/1462-2920.12497
Finn, D., Kopittke, P. M., Dennis, P. G., and Dalal, R. C. (2017). Microbial energy and matter transformation in agricultural soils. Soil Biol. Biochem. 111, 176–192. doi: 10.1016/j.soilbio.2017.04.010
Flores-Rentería, D., Sánchez-Gallén, I., Morales-Rojas, D., Larsen, J., and Ãlvarez-Sanchez, J. (2020). Changes in the abundance and composition of a microbial community associated with land use change in a Mexican tropical rain Forest. J. Soil Sci. Plant Nutr. 20, 1144–1155. doi: 10.1007/s42729-020-00200-6
George, P. B., Lallias, D., Creer, S., Seaton, F. M., Kenny, J. G., Eccles, R., et al. (2019). Divergent national-scale trends of microbial and animal biodiversity revealed across diverse temperate soil ecosystems. Nat. Commun. 10:1107. doi: 10.1038/s41467-019-09031-1
Godfray, H. C., Beddington, J. R., Crute, I. R., Haddad, L. J., Lawrence, D., Muir, J. F., et al. (2010). Food security: the challenge of feeding 9 billion people. Science 327, 812–818. doi: 10.1126/science.1185383
Goldfarb, K. C., Karaoz, U., Hanson, C. A., Santee, C. A., Bradford, M. A., Treseder, K. K., et al. (2011). Differential growth responses of soil bacterial taxa to carbon substrates of varying chemical recalcitrance. Front. Microbiol. 2:94. doi: 10.3389/fmicb.2011.00094
Guo, J. H., Liu, X. J., Zhang, Y., Shen, J. L., Han, W. X., Zhang, W. F., et al. (2010). Significant acidification in major Chinese croplands. Science 327, 1008–1010. doi: 10.1126/science.1182570
Jangid, K., Williams, M., Franzluebbers, A. J., Sanderlin, J. S., Reeves, J., Jenkins, M. B., et al. (2008). Relative impacts of land-use, management intensity and fertilization upon soil microbial community structure in agricultural systems. Soil Biol. Biochem. 40, 2843–2853. doi: 10.1016/J.SOILBIO.2008.07.030
Jenkins, S. N., Waite, I. S., Blackburn, A., Husband, R., Rushton, S. P., Manning, D. A., et al. (2009). Actinobacterial community dynamics in long term managed grasslands. Antonie Van Leeuwenhoek 95, 319–334. doi: 10.1007/s10482-009-9317-8
Ji, L., Guo, X., Sun, Q., and Wang, R. (2020). Soil microbial biomass and activity in different land use types in Ningxia irrigation area. Ecol. Env. Sci. 29, 516–524. doi: 10.16258/j.cnki.1674-5906.2020.03.011
Jones, R. T., Robeson, M. S., Lauber, C. L., Hamady, M., Knight, R., and Fierer, N. (2009). A comprehensive survey of soil acidobacterial diversity using pyrosequencing and clone library analyses. ISME J. 3, 442–453. doi: 10.1038/ismej.2008.127
Kösters, R., Preger, A. C., Preez, C. C., and Amelung, W. (2013). Reaggregation dynamics of degraded cropland soils with prolonged secondary pasture management in the south African Highveld. Geoderma 192, 173–181. doi: 10.1016/j.geoderma.2012.07.011
Kowalchuk, G. A., and Stephen, J. (2001). Ammonia-oxidizing bacteria: a model for molecular microbial ecology. Annu. Rev. Microbiol. 55, 485–529. doi: 10.1146/annurev.micro.55.1.485
Kugbe, J., and Issahaku, Z. (2015). Effects of soil conservation technologies in improving soil productivity in northern Ghana. J. Soil Sci. Environ. 6, 158–167. doi: 10.5897/jssem15.0490
Kunhikrishnan, A., Thangarajan, R., Bolan, N. S., Xu, Y., Mandal, S., Gleeson, D. B., et al. (2016). Functional relationships of soil acidification, liming, and greenhouse gas flux. Adv. Agron. 139, 1–71. doi: 10.1016/bs.agron.2016.05.001
Lauber, C. L., Strickland, M. S., Bradford, M. A., and Fierer, N. (2008). The influence of soil properties on the structure of bacterial and fungal communities across land-use types. Soil Biol. Biochem. 40, 2407–2415. doi: 10.1016/j.soilbio.2008.05.021
Lee-Cruz, L., Edwards, D. P., Tripathi, B. M., and Adams, J. M. (2013). Impact of logging and Forest conversion to oil palm plantations on soil bacterial communities in Borneo. Appl. Environ. Microbiol. 79, 7290–7297. doi: 10.1128/aem.02541-13
Li, X., Rui, J., Mao, Y., Yannarell, A. C., and Mackie, R. I. (2014). Dynamics of the bacterial community structure in the rhizosphere of a maize cultivar. Soil Biol. Biochem. 68, 392–401. doi: 10.1016/j.soilbio.2013.10.017
Lin, Y. T., Lin, Y. F., Tsai, I. J., Chang, E. H., Jien, S. H., Lin, Y. J., et al. (2019). Structure and diversity of soil bacterial communities in Ofshore Islands. Sci.Rep. 9:4689. doi: 10.1038/s41598-019-41170-9
Liu, J. J., Sui, Y. Y., Yu, Z. H., Yao, Q., Shi, Y., Chu, H. Y., et al. (2016). Diversity and distribution patterns of acidobacterial communities in the black soil zone of Northeast China. Soil Biol. Biochem. 95, 212–222. doi: 10.1016/j.soilbio.2015.12.021
Liu, K., Xue, Y. Q., Zhu, L., Xu, F., Zhu, Z., Zhang, T., et al. (2021). Effect of different land use types on the diversity of soil bacterial community in the coastal zone of Jialing River. Environ. Sci. 2021, 1–12. doi: 10.13227/j.hjkx.202106174
Lu, R. (1999). Soil agricultural chemical analysis method. (Beijing: Chinese Agriculture and Sciences Press).
Lynch, J. M., and Bragg, E. (1985). “Microorganisms and soil aggregate stability,” in Advances in soil science. Vol. 17. ed. B. A. Stewart (Germany: Springer), 133–171.
Magoc, T., and Salzberg, S. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963. doi: 10.1093/bioinformatics/btr507
McDaniel, M. D., Tiemann, L. K., and Grandy, A. S. (2014). Does agricultural crop diversity enhance soil microbial biomass and organic matter dynamics? A meta-analysis. Ecol. Appl. 24, 560–570. doi: 10.1890/13-0616.1
Mendez-Millan, M., Nguyen Tu, T. T., Balesdent, J., Derenne, S., Derrien, D., Egasse, C., et al. (2014). Compound-specific 13C and 14C measurements improve the understanding of soil organic matter dynamics. Biogeochemistry 118, 205–223. doi: 10.1007/s10533-013-9920-7
Mukhopadhya, I., Hansen, R., El-Omar, E. M., and Hold, G. L. (2012). IBD—what role do Proteobacteria play? Nat. Rev. Gastroenterol. Hepatol. 9, 219–230. doi: 10.1038/nrgastro.2012.14
Nacke, H., Thürmer, A., Wollherr, A., Will, C., Hodac, L., Herold, N., et al. (2011). Pyrosequencing-based assessment of bacterial community structure along different management types in German forest and grassland soils. PLoS One 6:e17000. doi: 10.1371/journal.pone.0017000
Naether, A., Foesel, B. U., Naegele, V., Wust, P. K., Weinert, J., Bonkowski, M., et al. (2012). Environmental factors affect acidobacterial communities below the subgroup level in grassland and forest soils. Appl. Environ. Microbiol. 78, 7398–7406. doi: 10.1128/aem.01325-12
Navarrete, A. A., Cannavan, F. S., Taketani, R. G., and Tsai, S. M. (2010). A molecular survey of the diversity of microbial communities in different Amazonian agricultural model systems. Diversity 2, 787–809. doi: 10.3390/d2050787
Navarrete, A. A., Venturini, A. M., Meyer, K. M., Klein, A. M., Tiedje, J. M., Bohannan, B. J. M., et al. (2015). Differential response of Acidobacteria subgroups to Forest-to-pasture conversion and their biogeographic patterns in the Western Brazilian Amazon. Front. Microbiol. 6:443. doi: 10.3389/fmicb.2015.01443
Nielsen, U. N., Wall, D. H., and Six, J. (2015). Soil biodiversity and the environment. Annu. Rev. Environ. Resour. 40, 63–90. doi: 10.1146/annurev-environ-102014-021257
Paul, E. A., and Clark, F. E. (1989). Soil microbiology and biochemistry. Academic Press, San Diego, USA.
Peng, M., Jia, H., and Wang, Q. (2017). The effect of land use on bacterial communities in saline-alkali soil. Curr. Microbiol. 74, 325–333. doi: 10.1007/s00284-017-1195-0
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219
Ramirez, K. S., Lauber, C. L., Knight, R., Bradford, M. A., and Fierer, N. (2010). Consistent effects of nitrogen fertilization on soil bacterial communities in contrasting systems. Ecology 91, 3463–3470. doi: 10.1890/10-0426.1
Rasche, F., Knapp, D. J., Kaiser, C., Koranda, M., Kitzler, B., Zechmeister-Boltenstern, S., et al. (2011). Seasonality and resource availability control bacterial and archaeal communities in soils of a temperate beech forest. ISME J. 5, 389–402. doi: 10.1038/ismej.2010.138
Rath, K. M., Fierer, N., Murphy, D. V., and Rousk, J. (2018). Linking bacterial community composition to soil salinity along environmental gradients. ISME J. 13, 836–846. doi: 10.1038/s41396-018-0313-8
Richter, A. S., Schöning, I., Kahl, T., Bauhus, J., and Ruess, L. (2018). Regional environmental conditions shape microbial community structure stronger than local forest management intensity. For. Ecol. Manag. 409, 250–259. doi: 10.1016/J.FORECO.2017.11.027
Rillig, M. C., and Mummey, D. L. (2006). Mycorrhizas and soil structure. New Phytol. 171, 41–53. doi: 10.1111/J.1469-8137.2006.01750.X
Roesch, L. F., Fulthorpe, R., Riva, A., Casella, G., Hadwin, A. M., Kent, A. D., et al. (2007). Pyrosequencing enumerates and contrasts soil microbial diversity. ISME J. 1, 283–290. doi: 10.1038/ismej.2007.53
Rousk, J., Bååth, E., Brookes, P. C., Lauber, C. L., Lozupone, C. A., Caporaso, J. G., et al. (2010). Soil bacterial and fungal communities across a pH gradient in an arable soil. ISME J. 4, 1340–1351. doi: 10.1038/ismej.2010.58
Serna-Chavez, H. M., Fierer, N., and Bodegom, P. M. (2013). Global drivers and patterns of microbial abundance in soil. Glob. Ecol. Biogeogr. 22, 1162–1172. doi: 10.1111/geb.12070
Smith, P., House, J. I., Bustamante, M. M., Sobocka, J., Harper, R., Pan, G., et al. (2015). Global change pressures on soils from land use and management. Glob. Chang. Biol. 22, 1008–1028. doi: 10.1111/gcb.13068
Sparks, D. L., Page, A., Helmke, P., Loeppert, R., Soltanpour, P., Tabatabai, M., et al. (1996). Methods of Soil Analysis. Part 3-Chemical Methods. Soil Science Society of America Inc.
Stephan, A., Meyer, A. H., and Schmid, B. (2000). Plant diversity affects culturable soil bacteria in experimental grassland communities. J. Ecol. 88, 988–998. doi: 10.1046/j.1365-2745.2000.00510.x
Tan, H., Barret, M., Mooij, M. J., Rice, O., Morrissey, J. P., Dobson, A. D., et al. (2012). Long-term phosphorus fertilisation increased the diversity of the total bacterial community and the phoD phosphorus mineraliser group in pasture soils. Biol. Fertil. 49, 661–672. doi: 10.1007/s00374-012-0755-5
Thomson, B. C., Tisserant, É., Plassart, P., Uroz, S., Griffiths, R. I., Hannula, S. E., et al. (2015). Soil conditions and land use intensification effects on soil microbial communities across a range of European field sites. Soil Biol. Biochem. 88, 403–413. doi: 10.1016/j.soilbio.2015.06.012
Treseder, K. K. (2008). Nitrogen additions and microbial biomass: a meta-analysis of ecosystem studies. Ecol. Lett. 11, 1111–1120. doi: 10.1111/j.1461-0248.2008.01230.x
Tripathi, B. M., Stegen, J. C., Kim, M., Dong, K., Adams, J. M., and Lee, Y. (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
Trivedi, P., Delgado-Baquerizo, M., Anderson, I. C., and Singh, B. K. (2016). Response of soil properties and microbial communities to agriculture: implications for primary productivity and soil health indicators. Front. Plant Sci. 7:990. doi: 10.3389/fpls.2016.00990
Trivedi, C., Delgado-Baquerizo, M., Hamonts, K., Lai, K., Reich, P. B., and Singh, B. K. (2019). Losses in microbial functional diversity reduce the rate of key soil processes. Soil Biol. Biochem. 135, 267–274. doi: 10.1016/j.soilbio.2019.05.008
Tsiafouli, M. A., Thébault, É., Sgardelis, S. P., Ruiter, P. C., Putten, W. H., Birkhofer, K., et al. (2015). Intensive agriculture reduces soil biodiversity across Europe. Glob. Chang. Biol. 21, 973–985. doi: 10.1111/gcb.12752
van der Heijden, M. G., Bardgett, R. D., and van Straalen, N. M. (2008). The unseen majority: soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems. Ecol. Lett. 11, 296–310. doi: 10.1111/J.1461-0248.2007.01139.X
van Horn, D., Okie, J. G., Buelow, H. N., Gooseff, M. N., Barrett, J. E., and Takacs-Vesbach, C. D. (2014). Soil microbial responses to increased moisture and organic resources along a salinity gradient in a Polar Desert. Appl. Environ. Microbiol. 80, 3034–3043. doi: 10.1128/AEM.03414-13
Wagg, C., Schlaeppi, K., Banerjee, S., Kuramae, E. E., and van der Heijden, M. G. (2019). Fungal-bacterial diversity and microbiome complexity predict ecosystem functioning. Nat. Commun. 10:4841. doi: 10.1038/s41467-019-12798-y
Wang, Q., Garrity, G. M., Tiedje, J., and Cole, J. R. (2007). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/AEM.00062-07
Wang, X., Li, F., Tang, K., Wang, Y., Suri, G., Bai, Z., et al. (2020). Land use alters relationships of grassland productivity with plant and arthropod diversity in inner Mongolian grassland. Ecol. Appl. 30:e02052. doi: 10.1002/eap.2052
Xie, K., Deng, Y., Zhang, S., Zhang, W., Liu, J., Xie, Y., et al. (2017). Prokaryotic community distribution along an ecological gradient of salinity in surface and subsurface saline soils. Sci. Rep. 7:13332. doi: 10.1038/s41598-017-13608-5
Xu, A., Liu, J., Guo, Z., Wang, C., Pan, K., Zhang, F., et al. (2021). Soil microbial community composition but not diversity is affected by land-use types in the agro-pastoral ecotone undergoing frequent conversions between cropland and grassland. Geoderma 401:115165. doi: 10.1016/J.GEODERMA.2021.115165
Yachi, S., and Loreau, M. (1999). Biodiversity and ecosystem productivity in a fluctuating environment: the insurance hypothesis. Proc. Natl. Acad. Sci. U. S. A. 96, 1463–1468. doi: 10.1073/PNAS.96.4.1463
Yergeau, É., Hogues, H., Whyte, L. G., and Greer, C. W. (2010). The functional potential of high Arctic permafrost revealed by metagenomic sequencing, qPCR and microarray analyses. ISME J. 4, 1206–1214. doi: 10.1038/ismej.2010.41
Yu, Z., Wang, G., Jin, J., Liu, J., and Liu, X. (2011). Soil microbial communities are affected more by land use than seasonal variation in restored grassland and cultivated Mollisols in Northeast China. Eur. J. Soil Biol. 47, 357–363. doi: 10.1016/J.EJSOBI.2011.09.001
Zhang, H., Ding, W., He, X., Yu, H., Fan, J., and Liu, D. (2014). Influence of 20–year organic and inorganic fertilization on organic carbon accumulation and microbial community structure of aggregates in an intensively cultivated Sandy loam soil. PLoS One 9:e92733. doi: 10.1371/journal.pone.0092733
Zhang, Y. T., Ding, K., Yrjälä, K., Liu, H. Y., Tong, Z. K., and Zhang, J. H. (2021). Introduction of broadleaf species into monospecific Cunninghamia lanceolata plantations changed the soil Acidobacteria subgroups composition and nitrogen-cycling gene abundances. Plant Soil 467, 29–46. doi: 10.1007/s11104-021-05014-8
Zhang, X., Xu, M., and Shi, F. (2020). Impact of typical agricultural land use on the characteristics of soil microbial communities in the Nyingchi region of southeastern Tibet. J. Agro-Env. Sci. 39, 331–342.
Zhang, M., Zhang, W., Bai, S. H., Niu, Y., Hu, D., Ji, H., et al. (2019). Minor increases in Phyllostachys edulis (Moso bamboo) biomass despite evident alterations of soil bacterial community structure after phosphorus fertilization alone: based on field studies at different altitudes. For. Ecol. Manag. 451:117561. doi: 10.1016/j.foreco.2019.117561
Zhao, X., Huang, J., Lu, J., and Sun, Y. (2019). Study on the influence of soil microbial community on the long-term heavy metal pollution of different land use types and depth layers in mine. Ecotoxicol. Environ. Saf. 170, 218–226. doi: 10.1016/j.ecoenv.2018.11.136
Zhao, J., Ni, T., Li, Y., Xiong, W., Ran, W., Shen, B., et al. (2014). Responses of bacterial communities in arable soils in a Rice-wheat cropping system to different fertilizer regimes and sampling times. PLoS One 9:e85301. doi: 10.1371/journal.pone.0085301
Keywords: forest-grass transition zone, grassland, land use, soil properties, soil bacteria
Citation: Bai Z, Zheng L, Bai Z, Jia A and Wang M (2022) Long-term cultivation alter soil bacterial community in a forest-grassland transition zone. Front. Microbiol. 13:1001781. doi: 10.3389/fmicb.2022.1001781
Edited by:
Bruno Tilocca, Magna Græcia University, ItalyReviewed by:
Gaosen Zhang, Chinese Academy of Sciences, ChinaCunqiang Ma, Huazhong Agricultural University, China
Copyright © 2022 Bai, Zheng, Bai, Jia and Wang. 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: Mingjun Wang, wangmingjun@neau.edu.cn