ORIGINAL RESEARCH article

Front. Environ. Sci., 04 July 2022

Sec. Toxicology, Pollution and the Environment

Volume 10 - 2022 | https://doi.org/10.3389/fenvs.2022.937688

Response of Soil Fertility and Bacterial Community Composition to Vegetation Species in a Coal Mining Subsidence Area: A Survey After 20-Year Reclamation

  • 1. School of Earth and Environment, Anhui University of Science and Technology, Huainan, China

  • 2. Engineering Laboratory of Anhui Province for Comprehensive Utilization of Water and Soil Resources and Construction of Ecological Protection in Mining Area with High Groundwater Level, Anhui University of Science and Technology, Huainan, China

  • 3. State Key Laboratory of Crop Stress Adaptation and Improvement, School of Life Sciences, Henan University, Kaifeng, China

  • 4. Key Laboratory of Tobacco Biology and Processing, Ministry of Agriculture and Rural Affairs, Tobacco Research Institute, Chinese Academy of Agricultural Sciences, Qingdao, China

Article metrics

View details

7

Citations

1,7k

Views

788

Downloads

Abstract

Revegetation is an important restoration approach after reclamation in coal mining subsidence area. However, few studies have paid attention to the impacts of different vegetation species on soil fertility and bacterial community composition in the reclamation area filled with fly ash for a long time. In this study, soil fertility and bacterial community composition were investigated in a non-subsidence area (FCK) and a coal mining subsidence reclamation area restored vegetation with woody plants (pagoda trees, FS; peach trees, FP) and herbaceous plants (wheat-maize rotation, FW) for 20 years. Results showed that topsoil and subsoil nutrients including available K and alkali-hydrolyzable N, soil organic matter, and total carbon in the non-subsidence site were significantly higher than those in reclamation sites. Topsoil fertility indices (SFI) in reclamation sites were lower than that in non-subsidence site, and soil fertility index in FW was higher than that in FS and FP. Moreover, the SFI in topsoil (from −0.24 to 2.08) was significantly higher than that in subsoil (from −1.12 to −0.39). 16S rRNA gene sequencing analysis showed Proteobacteria, Acidobacteria and Actinobacteria were the dominant bacterial phyla in all sites, but the subsoil bacterial alpha indices (Shannon and Simpson) in the non-subsidence site were higher than those in reclamation sites. Principal coordinates analysis and non-metric multidimensional scaling analysis showed that the bacterial community composition in reclamation sites was significantly different from that in the non-subsidence site, and they were similar in the soil with wheat-maize rotation model between non-subsidence and reclamation sites. Co-occurrence network analysis noted that the network of dominant bacterial operational taxonomic units in the subsoil was more complex than that in the topsoil. Redundancy analysis suggested soil pH in topsoil and subsoil was an important driving factor for soil bacteria community composition. Overall, the reclamation site with the wheat-maize rotation model had higher soil fertility and bacterial community composition was similar to that of the non-subsidence site, indicating revegetation with the wheat-maize rotation model is a good reclamation approach to improve soil fertility in coal mining subsidence area.

Introduction

With the development of industry, demands for coal resources are increasing in most countries over the world, especially in developing countries (Feng et al., 2019). For instance, coal consumption in China accounted for 56.8% of the total energy consumption in 2020 (Hu et al., 2022), and coal consumption will continue to grow by 3.4% in the next 20 years (Yu et al., 2022). The extensive exploitation of coal resources has caused serious environmental and social problems including land subsidence, land damage, and poor harvest (Wang et al., 2020; Ma et al., 2021). In high groundwater level areas, coal mining always resulted in land subsidence and changed terrestrial to aquatic ecosystems (Ma et al., 2021). Such dramatic environmental changes have destroyed agricultural lands and threatened the safe production of crops and the survival of local residents. Therefore, it is of great significance to repair the coal mining subsidence area.

Previous literature reviews have noted that fly ash, coal gangue, and other solid wastes can be used as filling materials in the coal mining subsidence areas and recover with surface soil above the filling materials (Ram and Masto, 2014; Hemalatha and Ramaswamy, 2017; Hu, 2019). Fly ash as a by-product of the coal-fired power station was widely used to fill the subsidence area (Hu, 2019). It can also provide lots of essential nutrients (K, Si, Fe, etc.) for plants (Ram and Masto, 2014), and its filling was beneficial for soil carbon sequestration and biodiversity conservation (Pandey, 2013). However, the filling substrates were significantly different from the subsoil and cannot effectively retain soil moisture and fertility, which may result in poor fertility and low soil productivity (Ma et al., 2019). For example, Duo and Hu (2018) reported that the soil quality significantly decreased after reclamation, and the crop yield was only 80.5% of that in the non-subsidence area after 5 years of reclamation (Duo and Hu, 2018). Consequently, the long-term changes in soil quality should be paid more attention in reclamation area filled with solid wastes/amendments.

Additionally, revegetation in reclamation area can effectively improve soil fertility and restore soil ecosystem functions (Zhou et al., 2020). It is essential to establish a self-sustaining revegetation ecosystem by selecting suitable plant species in the light of the condition of the subsidence areas (Kohler et al., 2016; Qu et al., 2017; Li and Liber, 2018). Numerous plants including woody plants (Acacia catechu, Dalbergia sissoo, Pongamia pinnata and Terminalia bellirica, etc.) (Srivastava et al., 2014) and herbaceous plants (grass, wheat, maize, Ricinus communis, Vigna radiata, Saccharum munja, etc.) (Pandey et al., 2012; Pandey, 2013; Yang et al., 2020) have been reported to revegetate in coal mining subsidence areas. However, little information is available on the impacts of vegetation species on soil microbial community composition in reclamation areas filled with fly ash. Soil microbial communities are critical to organic matter decomposition, nutrient cycling, long-term ecosystem stability and revegetation (Lehmann et al., 2011; Chen et al., 2017; Sun et al., 2019), and it is considered as an important indicator for the ecological restoration effect in coal mining subsidence area (Xu et al., 2019; Yuan et al., 2021). Natural succession of reclamation areas is a slow biogeochemical process and needs a long term before reaching a stable state (Ma et al., 2021). However, most studies about ecological restoration concentrated on a relatively short time (1–15 years) (Liu et al., 2018; Ma et al., 2019; Wang et al., 2020), and few studies were available on the evaluation of soil fertility and microbiome composition in reclamation areas.

In this study, a coal mining subsidence area filled with fly ash in eastern China was selected as the study area and soil samples were collected in four sites planted with woody plants (pagoda tree and peach tree) and herbaceous plants (wheat-maize rotation model) for 20 years. Here, we aim to 1) investigate soil physicochemical properties and fertility in reclamation areas; 2) evaluate the impacts of reclamation on soil bacterial community composition and biodiversity.

Materials and Methods

Site Description

This study was conducted in a coal mining subsidence area (∼100 ha, 116°85′-116°88′ E, 33°97′-33°98′ N) filled with fly ash in the eastern Huaibei city, Anhui Province, China (Supplementary Figure S1). The climate is characterized with a warm temperate semi-humid monsoon climate. The average annual precipitation and average annual temperature are 850 mm and 15.9°C in this area, respectively. Since high groundwater level existed above the coal bed, resulting in surface subsidence and land damage after long-term coal mining. In 2000, the surface soil (40 cm) was stripped and filled with fly ash in the subsidence area, and then covered the surface soil above the fly ash. During the next 20 years, most reclamation area was adopted a wheat-maize rotation model with two crops per year (FW), and part of the area was planted with pagoda trees (Sophora japonica Linn) (FS) and peach trees (Prunus persica L.) (FP), respectively. In addition, a non-subsidence site with a wheat-maize rotation model was also selected as the background area (FCK), which was 2.50 km to the reclamation area.

Sample Collection and Physicochemical Analysis

In November 2020, topsoil (0–20 cm) and subsoil (20–40 cm) samples were collected in each area. All soil samples were passed through a 2 mm sieve. Half of the soil samples were dried at room temperature for chemical analysis, and the other samples were stored at −40°C for bioanalysis.

Soil pH was analyzed in a soil-to-water ratio of 1:2.50 (w/v) by pH electrode (PHS-3C, Shanghai, China). Soil organic matter (SOM), available phosphorus (AP), total phosphorus (TP), alkali-hydrolyzable nitrogen (AN), available potassium (AK) and total potassium (TK) were measured according to Lu (2000). Detailed information about the analysis methods is provided in Supplementary Text S1. An elemental analyzer (Vario Macro Cube, Elementar Company, Germany) was used to determine soil total nitrogen (TN) and total carbon (TC) concentrations (Wang et al., 2021).

Soil Biological Analysis

Soil whole DNA was extracted and purified according to the method of Cui H. et al. (2018). The concentration and purity of DNA were checked using NanoDrop One (Thermo Fisher Scientific, MA, United States). The standard bacterial primer set (515F and 806R) combined with barcode sequences (Caporaso et al., 2011) were used to amplify 16S rRNA genes (Cui H. et al., 2018). The Illumina Hiseq2500 platform was used for PE250 sequencing of the constructed amplicon library. Detailed information about the analysis methods is provided in Supplementary Text S2.

To clarify the associations among the soil key operational taxonomic units (OTUs), we selected the top 30 OTUs with the highest relative abundance in all sites (Hill et al., 2016). Spearman’s correlation coefficient was chosen, and the threshold value of the correlation coefficient was 0.4. The network was analyzed by the tool of the Wekemo Bioincloud (https://www.bioincloud.tech). To evaluate the contributions of microbial populations to the overall community structure, the bipartite networks created in Cytoscape (3.8.2) were used to describe the associations between the OTUs and different sites (Supplementary Text S3) (Smoot et al., 2011; Jiao et al., 2017).

Statistical Analysis

The igraph and vegan in R (4.1.0) were used for statistical analysis The factor analysis was used to calculate the soil fertility index (SFI) (Supplementary Text S4) (Li et al., 2014; Guo et al., 2018). The relationships between soil physicochemical properties and microbial alpha diversity indexes were detected by Pearson correlation analysis. The phylum to genus-level abundance information was extracted by R to visualize using Origin 2021. Richness of OTUs, Chao1, ACE, Shannon, Simpson, and goods coverage indices were applied to characterize the bacterial alpha diversity. Principal coordinates analysis (PCoA) and non-metric multidimensional scaling (NMDS) were performed in the R software. The relationship between environmental parameters and bacteria was analyzed by redundancy analysis (RDA) which was performed by Lingbo MicroClass cloud application (http://www.cloud.biomicroclass.com/cloudplatform/softpage/CCA). The linear discriminant analysis (LDA) effect size (LEfSe) method (http://www.cloud.biomicroclass.com/CloudPlatform/SoftPage/LF2) was performed to identify multiple taxonomical levels of potential biomarkers (LDA score threshold >3.0, factor Kruskal–Wallis and Wilcoxon test with an alpha value of 0.05), the analysis strategy of multi-group classification analysis was used one against one (Chen et al., 2020).

Results and Discussion

Soil Physicochemical Properties and Fertility

Soil physicochemical properties affected by different vegetation species were shown in Table 1. Topsoil was alkalinity in this area and its pH ranged from 8.14 to 8.69. Soil pH in FS and FP (wood land) were 0.11–0.53 units higher than those of FCK and FW (wheat field), suggesting that vegetation species influence soil pH, which was consistent with the findings of Deng et al. (2019) and Luo et al. (2022). Generally, the concentrations of WC, AK, TN, AN, SOM, and TC in FCK were the highest, and followed by FW (Table 1). For example, in the topsoil, WC in FCK was 0.98 1.26 and 0.62 times higher than that in FS, FP, and FW, respectively, and 0.64, 0.98 and 0.41 times higher for SOM. These results suggested filling reclamation caused soil water loss and fertility degradation, which also have been reported in previous study (Qu et al., 2018). However, concentrations of plant nutrients (including AP, SOM, and TC) of topsoil in FW were significantly higher than that in FS and FP (p < 0.05, Table 1), which suggested the wheat-maize rotation model was a good approach for the enhancement of soil fertility.

TABLE 1

ParameterspHWC (%)TP (g kg−1)AP (mg kg−1)AK (mg kg−1)TK (g kg−1)TN (g kg−1)AN (mg kg−1)SOM (g kg−1)TC (g kg−1)
0–20 cmFCK8.38 ± 0.04b16.5 ± 0.86a0.71 ± 0.06a9.50 ± 0.78b328 ± 10a22.6 ± 0.58a1.09 ± 0.13a253 ± 9a29.8 ± 1.08a25.6 ± 2.04a
FS8.49 ± 0.07a8.36 ± 2.18bc0.67 ± 0.04a9.07 ± 2.15b292 ± 22b24.3 ± 0.78a0.66 ± 0.10b116 ± 16b18.1 ± 2.91bc18.0 ± 0.87bc
FP8.57 ± 0.03a7.30 ± 1.10c0.71 ± 0.10a9.15 ± 1.14b254 ± 25c17.3 ± 0.36b0.47 ± 0.05c68 ± 5c15.0 ± 2.27c15.7 ± 1.21c
FW8.14 ± 0.08c10.2 ± 1.29b0.76 ± 0.02a14.6 ± 0.28a290 ± 12b18.0 ± 0.59b0.82 ± 0.02b125 ± 5b21.1 ± 2.81b19.7 ± 1.55b
20–40 cmFCK8.69 ± 0.02a19.7 ± 0.6a0.72 ± 0.1a10.4 ± 3.04a281 ± 13a24.2 ± 0.32a0.55 ± 0.06a126 ± 27a17.7 ± 0.21a20.6 ± 0.20a
FS8.62 ± 0.06ab9.74 ± 1.83b0.76 ± 0.06a10.3 ± 0.15a226 ± 24b18.2 ± 0.68b0.48 ± 0.08ab70 ± 1b13.7 ± 2.62b17.6 ± 1.68b
FP8.55 ± 0.06c10.2 ± 1.69b0.76 ± 0.08a8.65 ± 0.29a217 ± 21b18.1 ± 0.27b0.48 ± 0.04ab56 ± 7b10.9 ± 2.10b17.7 ± 0.42b
FW8.63 ± 0.05ab17.7 ± 2.41a0.76 ± 0.05a9.53 ± 2.02a201 ± 7b17.3 ± 0.68b0.42 ± 0.04c67 ± 4b11.8 ± 2.30bc16.4 ± 0.69b

Impacts of vegetation species on the soil physicochemical properties.

Notes: Values are means ± standard deviation (n = 3) and different letters within the same column denote significant differences (p < 0.05) among different sites. WC, soil water content; TN, total nitrogen; TP, total phosphorus; TK, total potassium; AN, alkali-hydrolyzable nitrogen; AP, available phosphorus; AK, available potassium; SOM, soil organic matter; TC, total carbon. Site abbreviations are defined in Figure 1.

Similar to topsoil, most subsoil plant nutrients concentrations (including AK, TK, TN, AN, SOM, and TC) in FCK were higher than those in the other sites (Table 1). Moreover, subsoil pH (8.35–8.69) was higher than that of topsoil (8.05–8.57), and a similar trend was found for soil moisture content. In contrast to soil pH, contents of AK, TN, AN, and SOM in topsoil were higher than those in the subsoil. For example, topsoil SOM in FCK, FS, FP, and FW were 0.69, 0.32, 0.38, and 0.79 times higher than those of subsoil, respectively. In addition, the principal component analysis showed the eigenvalues of PC1 and PC2 were 3.61 and 1.47, respectively (Supplementary Table S1). The variance contribution rates were 51.6 and 21.0%, which can explain 72.6% of the total variance contribution. The score matrix of the PCA analysis showed that the loadings of TN, AK, and AN in PC1 exceeded 0.87 (Supplementary Table S2), indicating that these factors in PC1 play a significant role in soil fertility changes. SFI in the topsoil of FCK (2.08) was significantly greater than those of FS (−0.24) and FP (−0.56) (p < 0.05) (Supplementary Table S3), but there was no noticeable difference between FCK and FW (1.58) (p > 0.05, Supplementary Table S3). These results suggested that the soil fertility level in reclamation areas with woody plants was lower than that in the non-subsidence area with agricultural crops, and they also noted that the wheat-maize rotation model is an effective method to improve the soil fertility in the reclamation area.

Correlation analysis showed AK, AN, SOM, and TC were positively correlated with each other (p < 0.05, Supplementary Figure S2), especially for the positive correlation between TC and TN (p < 0.01). It may be due to the addition of nitrogen fertilizers that can promote plant growth and carbon accumulation (Qi et al., 2020). Moreover, we found the C/N ratios (Supplementary Figure S3) in the topsoil of wood land (33.33–33.78, FS and FT) were higher than those in the agricultural field (23.47–23.98, FCK and FW), which was also reported by Gao et al. (2014). This may be due to the intensive cultivation and agricultural management practices in agricultural fields resulting in the mineralization and loss of SOM to the atmosphere (Lal, 2002; Gao et al., 2014).

Soil Bacterial Community Composition

In all sites, 360,000 high-quality sequences and 7,141 different OTUs were obtained. The abundance of the dominant bacterial phyla/classes between reclamation and non-subsidence areas was shown in Figure 1 and Supplementary Figure S4. Proteobacteria, Acidobacteria, and Actinobacteria were the most dominant bacterial phyla in the topsoil and subsoil, but their percentages were significantly different between topsoil and subsoil, suggesting the bacterial community composition was significantly different. For example, Proteobacteria in the topsoil of FCK accounted for 27.5%, which was 0.09 times lower than that in the subsoil of FW, and Acidobacteria in the subsoil of FS was higher compared to that in the topsoil of FS. Moreover, Proteobacteria in the topsoil of FW accounted for 40.1%, which was 0.46, 0.68 and 0.73 times higher than that in FCK, FS and FP, respectively. Similar to Wang et al. (2017), Gammaproteobacteria and Alphaproteobacteria were the dominant bacterial classes in the topsoil of the reclamation area, and accounted for 13.3 and 9.24% in all sites (Figure 1; Supplementary Table S4). Previous studies noted Gammaproteobacteria and Alphaproteobacteria participated in the process of nitrogen fixation (Madigan et al., 1984; Wolińska et al., 2020), thus their existence with large amounts in the reclamation area favored the enhancement of soil quality. In addition, the highest relative abundance of bacterial phylum and classes of FW was higher than that of FCK, FP and FS, indicating that the wheat-maize rotation model is a good selection to reconstruct the bacterial community in the reclamation area.

FIGURE 1

FIGURE 1

Relative abundances of dominant bacterial phyla/classes in topsoil (A) and subsoil (B) among different sites. Relative abundances are based on the proportional frequencies of 16S rDNA sequences that could be classified at the phyla/classes level. Note: All, the whole soil samples; FCK, non-subsidence site (wheat-maize rotation); FS, reclamation site with pagoda trees; FP, reclamation site with peach trees; FW, reclamation site with wheat-maize rotation.

The bacterial richness in the reclamation sites (FS, FP and FW) for the topsoil were not significantly different from that in the non-subsidence site (FCK) (Table 2). There were no significant variations in the Chao1 (5,038–5,255) and ACE (4,934–5,195) indices of topsoil and subsoil among the four sites yet. These results indicated topsoil bacterial richness indices in the reclamation sites had reached a steady state within 20 years compared to the non-subsidence site. Previous studies have also reported soil bacterial communities reached a steady state within 15 and 20 years in reclamation areas (Li et al., 2018; Zhao et al., 2019). Moreover, Shannon and Simpson indices in the topsoil of non-subsidence site (FCK) were significantly higher than those in reclamation sites (Table 2), suggesting that different bacterial diversity existed between non-subsidence and reclamation sites. For example, Shannon indices in the topsoil of FK (non-subsidence site) were higher by 0.41, 0.14 and 0.29 than those in FS, FP and FW (reclamation sites) (p < 0.05). However, Shannon and Simpson indices for the subsoil had no significant differences between non-subsidence (FK) and reclamation sites (FS, FP and FW).

TABLE 2

ParametersRichness of OTUsAlpha diversity
Chao1ACEShannonSimpsonGoods coverage
0–20 cmFCK4,542 ± 41a5,150 ± 164a5,050 ± 72a10.4 ± 0.08a0.9982 ± 0.0001a0.9938 ± 0.0004a
FS4,523 ± 145a5,102 ± 166a5,011 ± 154a10.0 ± 0.11c0.9963 ± 0.0013b0.9937 ± 0.0002a
FP4,580 ± 81a5,090 ± 57a5,008 ± 81a10.3 ± 0.08ab0.9974 ± 0.0002ab0.9942 ± 0.0003a
FW4,601 ± 91a5,122 ± 137a5,056 ± 117a10.1 ± 0.14bc0.9966 ± 0.0006b0.9941 ± 0.0003a
20–40 cmFCK4,410 ± 298b5,038 ± 285a4,934 ± 271a10.2 ± 0.21a0.9974 ± 0.0006a0.9934 ± 0.0001c
FS4,608 ± 183ab5,100 ± 189a5,036 ± 163a10.1 ± 0.11a0.9965 ± 0.0006ab0.9941 ± 0.0000ab
FP4,580 ± 74ab5,153 ± 6a5,021 ± 44a10.1 ± 0.03a0.9961 ± 0.0003b0.9940 ± 0.0002b
FW4,810 ± 102a5,255 ± 106a5,195 ± 104a10.3 ± 0.11a0.9972 ± 0.0002a0.9944 ± 0.0001a

Impacts of vegetation species on the richness of OTU and alpha diversity indices.

Notes: Values are means ± standard deviation (n = 3) and different letters within the same column denote significant differences (p < 0.05) among different sites. Site abbreviations are defined in Figure 1.

As shown in Figure 2, PCoA and NMDS analysis noted that microbial composition in woody plant revegetated sites (FS, FP) clustered together in the topsoil and subsoil, and a similar result was found by the cluster analysis at the phylum level (Supplementary Figure S4). Because woody plants had a robust root system and considerable amounts of exudates and mucilage which can strongly absorb soil nutrients (Yang et al., 2020), resulting in an substantial impact on the soil microbial community. These findings suggested that vegetation species significantly influenced soil bacterial community composition (Kohler et al., 2016). Moreover, the PCoA and NMDS (Figure 2) showed the soil bacterial community composition in the non-subsidence site was significantly differentiated from reclamation sites, indicating that reclamation changes the soil environment and influences the soil microbial community structure. It may be attributed to soil fertility in FCK being higher than that in the reclamation sites. In agreement with our study, Xu et al. (2020) also reported that soil organic carbon, inorganic nitrogen and other nutrients were the key factors affecting soil microbial diversity. Moreover, the PCoA and NMDS analysis indicated bacterial community composition in FW was the closest to the site of FCK (Figure 2) compared to other reclamation sites, thus the bacterial community composition in the two sites was similar and which may be due to the same wheat-maize rotation cultivation model. In addition, PCoA and NMDS analysis markedly separated topsoil and subsoil into two clusters (Figure 2), which indicated that the bacterial community composition in topsoil was different from that in the subsoil at each site, and it was in agreement with the result in Figure 1.

FIGURE 2

FIGURE 2

Principal coordinates analysis (PCoA) ordination (A) and non-metric multidimensional scaling (NMDS) ordination (B) based on Bray-Curtis distance at OTU level showing the effects of different vegetation species on bacterial communities among different sites. The circles represent topsoil. The triangles represent subsoil. The colors of the nodes represent different sites. Note: Site abbreviations are defined in Figure 1.

Identification of Biomarkers and Co-Occurrence Networks

LEfSe was used to identify the different biomarkers. As shown in Supplementary Figure S5, more bacterial taxa (93 clades, 15 phyla, 25 classes, 39 orders, and 27 families) in the topsoil were identified than that in the subsoil (78 clades, 12 phyla, 24 classes, 28 orders, and 25 families), which indicated that soil bacteria in the topsoil were more sensitive to the environmental disturbance than those in the subsoil. These biomarkers can be regarded as key microbial taxa shaping the metabolic threshold of the microbial community (Cui Y. et al., 2018). Chen et al. (2020) also noted biomarkers vary depending on the different environments and the relative abundance of core microorganisms. Moreover, the results of LDA revealed 37, 27, 13 and 16 bacterial phylotypes in the topsoil for FCK, FS, FP, and FW (Supplementary Figure S6A), and 27, 10, 30 and 11 bacterial phylotypes in the subsoil for FCK, FS, FP, and FW (Supplementary Figure S6B), respectively. The taxa and quantities of biomarkers were different between the non-subsidence and reclamation sites, which may be caused by the various soil quality and vegetation species.

To describe the interactions among microbial communities, we performed a network analysis of the top 30 OTUs in abundance (Figure 3) and matched them in the NCBI database (Supplementary Table S5, S6). In the topsoil, OTU31 had 9 significant positive correlations and 1 significant negative correlation with other OTUs (p < 0.05), and OTU126 and OTU23 had 13 and 11 significant positive correlations, and 2 and 4 significant negative correlations with other OTUs (p < 0.05) in the subsoil (Figure 3). The above three OTUs all significantly correlated with the top 30 OTUs, thus they can be recognized as keystone taxa and play important roles in the soil reclamation process (Ma et al., 2020). Moreover, the network of dominant bacterial OTUs in the subsoil was more complex than that in the topsoil (Figure 3), suggesting the dominant OTU in the subsoil were more closely connected than those in the topsoil.

FIGURE 3

FIGURE 3

The network of dominant bacterial OTUs of topsoil (A) and subsoil (B) in all sites. The nodes in the network represent the top 30 OTUs. The size of each node is proportional to the relative abundance. The red edge represents the positive correlation between OTUs, the blue edge represents the negative correlation. The colors of the nodes represent different soil bacteria phyla. Note: Site abbreviations are defined in Figure 1.

Nowadays, rare microorganisms have attracted more attention for their important ecological role in maintaining microbial diversity and mediating biogeochemical processes (Sauret et al., 2014; Jiao et al., 2017). To analyze the connection between rare microorganisms and different sites. The contribution of microbial taxa to the community structure and the associations between genus-level OTUs and sites were evaluated by the bipartite association network (Supplementary Figure S7). Similar to Cao et al. (2021), we found the number of OTUs classified as genera in bacterial communities was at a low level (less than 15%). We screened a total of 587 nodes and 4315 edges in all sites, and they can be divided into 10 clusters and several keystone taxa according to the relationship between bacteria and sites. In co-occurrence network analysis, intermediate centrality represents the influence potential of an individual on the interaction between other nodes of the network (Greenblum et al., 2012). The core cluster was consisted by 451 nodes, and it accounted for 76.5% of the total nodes. The bipartite network analysis noted that most of OTUs were cross combinations, meaning similar taxa were presented among different sites, but their abundances were significantly different among all sites. In addition, several OTUs were connected to only a few sites, which may have specific environmental significance. For example, 27F-1492R_uncultured delta proteobacterium was a common genus of FCK and FW in the topsoil, which may be a keystone taxon in wheat-maize rotation mode.

Impacts of Soil Physicochemical Properties on Soil Microbial Community

Figure 4 showed the relationship between the top ten bacterial phyla and soil physicochemical properties. The first and second axes of RDA explained 89.9 and 6.44% of the variation in the bacterial community at the OTU level in the topsoil, while 65.8 and 16.2% were explained in the subsoil. (Figure 4). The total variation in bacterial community explained by the two axes exceeded 80% in both soil layers. Bacterial composition in all sites was significantly influenced by pH (p = 0.044) in the topsoil, and was significantly influenced by pH (p = 0.003), WC (p = 0.001), AK (p = 0.010) and TK (p = 0.003) in the subsoil (Supplementary Table S7). Moreover, Acidobacteria in the topsoil was positively correlated with pH, and Actinobacteria and Bacteroidetes were positively correlated with SOM, TN, TP, and AP (Figure 4). For the subsoil, Rokubacteria was positively correlated with pH, while Chloroflexi and Acidobacteria were positively correlated with AN, SOM, AP, TC, AK, and TN. These results suggested soil pH was a key driving factor for bacterial community structure (Ma et al., 2020; Ma et al., 2021), and further confirmed that soil nutrients can significantly influence soil bacterial communities in different reclamation areas (Wan and He, 2020).

FIGURE 4

FIGURE 4

Redundancy analysis (RDA) based on the bacterial community composition and soil characteristics in topsoil (A) and subsoil (B) among different sites. The red arrows represent soil physicochemical properties. Arrow length indicates the strength of correlation with a given ordination axis. The circles represent soil samples from different sites. The diamonds represent the major soil bacteria phylum. Note: Site abbreviations are defined in Figure 1. Soil physicochemical property abbreviations are defined in Table 1.

In the topsoil, WC was positively correlated with the Shannon and Simpson indices (p < 0.05), but WC was negatively correlated with the Richness of OTUs (p < 0.05) in the subsoil (Table 3). For subsoil, soil pH was markedly negatively correlated with and Richness of OTUs (p < 0.05). Moreover, AN and SOM in the topsoil were positively correlated with alpha diversity indices, this may be because carbon and nitrogen sources are essential nutrients for microbial survival, and they are also important driving factors for the evolution of soil bacterial diversity (Liu et al., 2020). Interestingly, we found that the correlations between alpha-diversity indices and physicochemical properties (TP, AP, SOM, and pH) in the topsoil overall differed from those in the subsoil, which indicated that soil microorganisms in different soil layers were sensitive to the changes in soil physicochemical properties.

TABLE 3

ParameterspHWCTPAPAKTKTNANSOMTC
0–20 cmRichness of OTUs0.27−0.06−0.20−0.25−0.090.10−0.090.010.09−0.16
Chao10.230.05−0.20−0.290.070.190.050.170.130.00
ACE0.210.05−0.21−0.230.000.170.020.110.13−0.04
Shannon0.090.66*0.21−0.24−0.09−0.170.330.250.480.40
Simpson0.240.66*0.18−0.45−0.190.020.290.310.340.41
Goods coverage0.33−0.20−0.09−0.21−0.35−0.10−0.28−0.230.04−0.37
20–40 cmRichness of OTUs−0.60*−0.310.450.060.17−0.560.10−0.37−0.180.06
Chao1−0.57−0.240.470.070.24−0.350.19−0.34−0.100.24
ACE−0.57−0.250.460.130.24−0.410.18−0.28−0.060.15
Shannon−0.170.180.21−0.010.07−0.13−0.10−0.24−0.04−0.08
Simpson0.300.62*−0.120.400.130.48−0.130.220.37−0.03
Goods coverage−0.53−0.390.29−0.20−0.13−0.89***−0.19−0.52−0.49−0.32

Pearson’s correlation between alpha diversity indices of bacterial community and soil physicochemical properties.

Notes: Soil physicochemical property abbreviations are defined in Table 1. *, **, and *** indicate significant correlation at p < 0.05, p < 0.01, and p < 0.001.

Conclusion

In this study, we found soil fertility in the non-subsidence site was significantly higher than that in reclamation sites, and topsoil SFI was higher than that in the subsoil. Furthermore, topsoil SFI in FW was higher than that in FS and FP. Proteobacteria, Acidobacteria and Actinobacteria were the dominant bacterial phyla in all sites, and topsoil bacterial Shannon and Simpson indices in FCK were significantly higher than those in reclamation sites. Soil bacterial community in reclamation sites was different from that in non-subsidence site based on the PCoA and NMDS analysis, and the bacterial community composition in FW was more similar to FCK compared with FS and FP. Moreover, co-occurrence network analysis showed topsoil had a more complex network of dominant OTUs compared with that of subsoil. Our study suggested that soil pH was an important driving factor for soil bacterial community composition. Results of this study suggested soil fertility in reclamation sites with a wheat-maize rotation model was higher and it had similar bacterial community composition to that in the non-subsidence site. Therefore, the wheat-maize rotation model is a good choice to improve soil fertility and soil bacterial community composition in coal mining subsidence area filled with fly ash.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions

Conceptualization: YW, YF, YS, and XZ. Methodology: YW, YF, YS, and XZ. Formal analysis: QW, SZ, and YS. Funding acquisition: YF, XZ, and YW. Project administration: YF, XZ, and YW. Resources: YW and YF. Supervision: XZ, YF, and YW. Validation: XZ, YW, and YF. Investigation: YW, YF, QW, SZ, and XZ. Visualization: YF, QW, SZ, YS, and YW. Writing—original draft: YW, YF, SZ, YS, and XZ Writing—review and editing: YW, YF, YS, and XZ.

Funding

This research was supported by the Natural Science Foundation of Universities of Anhui Province (KJ2020ZD35), the Graduate Innovation Fund Project of Anhui University of Science and Technology (2021CX2011), the Technology R&D Project of Huaibei Mining Group (2020-113 and 2022-103), the National Natural Science Foundation of China (31901195), and the Shandong Provincial Natural Science Foundation (ZR2019BD062).

Conflict of interest

The authors declare that this study received funding from Huibei Mining (Group) Co., Ltd. The funder was not involved in the study design, collection, analysis, interpretation of data, the writing of this article or the decision to submit it for publication.

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/fenvs.2022.937688/full#supplementary-material

References

Summary

Keywords

soil bacterial community composition, reclamation, coal mining subsidence area, revegetation, bacterial alpha diversity

Citation

Wang Y, Fan Y, Wang Q, Zhang S, Shi Y and Zheng X (2022) Response of Soil Fertility and Bacterial Community Composition to Vegetation Species in a Coal Mining Subsidence Area: A Survey After 20-Year Reclamation. Front. Environ. Sci. 10:937688. doi: 10.3389/fenvs.2022.937688

Received

06 May 2022

Accepted

16 May 2022

Published

04 July 2022

Volume

10 - 2022

Edited by

Zhu Li, Institute of Soil Science (CAS), China

Reviewed by

Yi Cheng, Nanjing Normal University, China

Chao Ma, Anhui Agricultural University, China

Updates

Copyright

*Correspondence: Yuchao Fan, ; Xuebo Zheng,

This article was submitted to Toxicology, Pollution and the Environment, a section of the journal Frontiers in Environmental Science

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics