-
Soil properties among different treatments are listed in Table 1. For basic physicochemical and enzymatic properties, compared with UL, TN was significantly increased in NPK, M and MNPK (P < 0.05), AP was significantly increased in NPK and MNPK (P < 0.05), TK was significantly decreased in the other four treatments (P < 0.05). The highest pH was in UL. Compared with CK, NPK, M and MNPK increased contents of TN, AP, AK and CAT, but decreased contents of SM, AN and TP. Compared with CK, TN and CAT were significantly higher in M and MNPK (P < 0.05).
Table 1. Soil biophysicochemical properties in different treatments
Index Treatments UL CK NPK M MNPK SM (g·kg−1) 0.17 ± 0.02ab 0.19 ± 0.02a 0.15 ± 0.01b 0.18 ± 0.01ab 0.15 ± 0.01b pH 8.27 ± 0.09a 7.86 ± 0.49a 7.85 ± 0.42a 8.15 ± 0.06a 8.14 ± 0.13a TN (g·kg−1) 0.27 ± 0.04c 0.32 ± 0.06bc 0.40 ± 0.07b 0.56 ± 0.01a 0.59 ± 0.06a AN (mg·kg−1) 304.05 ± 0.00ab 344.59 ± 11.50a 302.85 ± 30.84ab 315.97 ± 145.67ab 174.08 ± 99.02b TP (g·kg−1) 0.38 ± 0.00a 0.69 ± 0.46a 0.39 ± 0.20a 0.35 ± 0.28a 0.39 ± 0.20a AP (mg·kg−1) 5.09 ± 0.00b 10.93 ± 0.26ab 17.68 ± 5.80a 12.14 ± 0.80ab 18.65 ± 7.16a TK (g·kg−1) 5.61 ± 0.00a 3.72 ± 0.74b 2.88 ± 1.40b 3.20 ± 0.03b 4.01 ± 0.70b AK (mg·kg−1) 80.06 ± 0.00a 74.71 ± 8.34a 77.39 ± 8.34a 86.73 ± 10.08a 86.74 ± 15.17a CAT (ml·g−1) 2.15 ± 0.06a 1.52 ± 0.22b 1.65 ± 0.03b 2.22 ± 0.05a 1.97 ± 0.26a SOC (g·kg−1) 1.74 ± 0.00c 1.55 ± 0.08c 2.07 ± 0.17b 3.35 ± 0.24a 3.12 ± 0.04a MBC (mg·kg−1) 9.85 ± 0.45d 17.06 ± 1.50c 19.72 ± 0.58b 37.78 ± 2.29a 37.45 ± 0.82a DOC (g·kg−1) 0.09 ± 0.01a 0.09 ± 0.02a 0.08 ± 0.01a 0.08 ± 0.01a 0.08 ± 0.01a POC (g·kg−1) 0.24 ± 0.23a 0.28 ± 0.15a 0.52 ± 0.12a 0.60 ± 0.18a 0.71 ± 0.57a EOC (g·kg−1) 1.58 ± 0.22a 1.48 ± 0.58a 1.34 ± 0.35a 1.93 ± 0.46a 1.65 ± 0.69a Values indicate mean ± standard deviations (n = 3). Different letters in each row represent a significant difference among treatments (one-way ANOVA, P < 0.05). For soil carbon relevant properties, NPK, M and MNPK increased POC, and significantly increased SOC, MBC compared with UL and CK (P < 0.05). MBC was significantly higher in four other treatments than UL (P < 0.05).
RubisCO activity and cbbL gene abundance in different treatments
-
Soil RubisCO activity are depicted in Fig.1a and significant differences were observed among treatments. Compared with other treatments, M treatment significantly increased RubisCO activity (P < 0.05). RubisCO activity was the lowest in NPK and had no significant difference among UL, CK and MNPK (P > 0.05).
Figure 1.
(a) Soil RubisCO activity and (b) cbbL gene abundances in different treatments. Values are means (n = 3), and error bars represent standard deviation. Different lowercase letters above columns indicate significant differences (one-way ANOVA, P < 0.05) among treatments.
Soil cbbL gene abundance was depicted in Fig. 1b and was significantly the highest in M (2.07 × 107 copies/g dry soil, P < 0.05), which was consistent with RubisCO activity. Moreover, the ratio of cbbL/16S rRNA in M was also significantly the highest (3.37%) (Supplemental Fig. S2). The lowest cbbL gene abundance was in CK (7.92 × 105 copies/g dry soil). Compared with CK, cbbL gene abundance was increased 4.18, 25.18 and 8.74 times in NPK, M and MNPK, respectively. But there was no statistical difference among UL, CK, NPK and MNPK treatments.
Alpha diversity and composition of carbon-fixing microbial community in different treatments
-
A total of 6,140,983 raw sequences were obtained. Each sample contained 34931 high quality sequences after quality filtering and subsampling (normalizing the sequence number according to the minimum sample). All the sequences were further classified into 17 phyla, 43 classes, 84 orders, 160 families, 361 genera and 3719 OTUs. The rarefaction curve (Supplemental Fig. S3) showed that the current sampling depth included most carbon-fixing microbial taxa in the samples and was sufficient for further analyses.
Based on sequencing data, the α diversity indices (Sobs Richness, Shannon-Weaver Diversity, Pielou Evenness indices) of the carbon-fixing microbes in different treatments are depicted in Supplemental Fig. S4. The highest values of Sobs, Shannon-Weaver and Pielou indices were in MNPK, and Sobs index was significantly higher than the other four treatments (P < 0.05). The Shannon-Weaver and Pielou indices were higher in NPK, M and MNPK than CK. They suggested that fertilization increased soil microbial diversity and evenness, especially in the MNPK, and MNPK also significantly improved soil microbial richness (P < 0.05).
Furthermore, carbon-fixing microbial community composition at phylum level and class level are shown in Fig. 2a & b, respectively. The dominant phyla were Proteobacteria (73.26%−86.23%), Cyanobacteria (5.58%−17.69%) and Actinobacteria (2.16%−5.40%), comprising 91.15%−96.00% of the sequences. The dominant classes were Alphaproteobacteria (43.07%−70.18%), Betaproteobacteria (11.68%−21.60%), Gammaproteobacteria (1.77%−8.55%), Cyanophyceae (4.81%−15.14%), Cyanobacteria (0.64%−1.96%) and Actinobacteria (1.93%−5.04%), comprising 90.36%−95.78% of the sequences. It is worth noting that carbon-fixing fungi and archaea were also detected, which were not reflected in the figure due to relatively low abundance. In addition, the top 10 dominant genera were Bradyrhizobium, Rhodopseudomonas, Noviherbaspirillum, Cyanobium, Variovorax, Devosia, Marichromatium, Mesorhizobium, Nitrobacter and Thermoleptolyngbya, comprising 66.84%−88.26% of the sequences.
Figure 2.
Taxonomic composition of carbon-fixing microbial communities in soils at (a) phylum level and (b) class level.
Supplemental Table S1 showed that the relative abundances of the phyla Proteobacteria, Cyanobacteria and Actinobacteria were significantly different among treatments. Compared with UL, the other four treatments increased Proteobacteria, while significantly decreased Cyanobacteria (P < 0.05). NPK, M and MNPK increased Cyanobacteria and Actinobacteria, but decreased Proteobacteria compared with CK.
Supplemental Table S2 showed that the relative abundances of the classes Alphaproteobacteria, Gammaproteobacteria, Cyanophyceae, Cyanobacteria and Actinobacteria were significantly different among treatments. Compared with UL, the other four treatments increased Alphaproteobacteria, but significantly decreased Gammaproteobacteria, Cyanophyceae and Cyanobacteria (P < 0.05). NPK, M and MNPK increased Cyanobacteria and Actinobacteria, but decreased Alphaproteobacteria compared with CK.
LefSe analysis was applied to ascertain the biomarkers, which were significantly enriched carbon-fixing microbial taxa in certain treatments (Fig. 3a). There were 17 biomarkers with LDA score > 3.5, 16 of which belonged to Proteobacteria phylum. Therefore, reclamation mainly influenced the assembly of proteobacteria, which was consistent with the results of carbon-fixing microbial community composition (Fig. 2a). Compared with UL, class Alphaproteobacteria was significantly enriched in four other treatments; family Bradyrhizobiaceae was significantly enriched in NPK and MNPK; family Comamonadaceae was significantly enriched in M and MNPK; family Phyllobacteriaceae, genus Mesorhizobium and genus Nitrobacter were significantly enriched in MNPK; family Devosiaceae and genus Devosia were significantly enriched in CK. Compared with four other treatments, five biomarkers were significantly enriched in UL, namely class Gammaproteobacteria, order Chromatiales, family Chromatiaceae and genus Marichromatium, i.e. these taxa were significantly decreased in all bioreclamation treatments (CK, NPK, M and MNPK).
Figure 3.
Linear discriminant analysis (LDA) effect size analysis determined biomarkers (a) between UL and other treatments (CK, NPK, M and MNPK) and (b) between CK and fertilization treatments (NPK, M and MNPK). The LDA score indicates the effect size and ranking of each differentially abundant taxon (P < 0.05, LDA score > 3.5, a) (P < 0.05, LDA score > 3.0, b). The ordinate is the taxon with significant difference among groups, and the abscissa is a bar chart to visually show the LDA log score of each taxon. Blue, red, gray and orange bars represent the biomarkers in CK, NPK, M and MNPK having significantly greater abundances than in UL, respectively (a). Green bars represent the biomarkers in UL having significantly greater abundances than in all the other four treatments (a). Gray orange and red bars represent the biomarkers in M MNPK and NPK having significantly greater abundances than in CK, respectively (b). Blue bars represent the biomarkers in CK having significantly greater abundances than in all the other three treatments (b).
There were 17 biomarkers with LDA score > 3.0 (Fig. 3b), 10, five and two of them respectively belonged to Proteobacteria, Actinobacteria and Cyanobacteria phylum. Therefore, fertilization mainly influenced the assembly of these three phyla, which was consistent with the results of taxa composition mentioned in Fig. 2a. Compared with CK, genus Nitrobacter was significantly enriched in NPK M and MNPK; family Phyllobacteriaceae and genus Mesorhizobium were significantly enriched in M and MNPK; class Actinobacteria was significantly enriched in NPK and MNPK; Phylum Cyanobacteria, class Cyanophyceae, family Thiobacillaceae, genus Sulfuritortus, family Thermomonosporaceae and genus Thermomonospora were significantly enriched in NPK. Compared with NPK, M and MNPK, family Devosiaceae and genus Devosia were significantly enriched in CK, which means these taxa were significantly decreased in all fertilization treatments (NPK, M and MNPK).
Carbon-fixing microbial community structure in different treatments
-
Afterwards, PCA and biclustering heatmap analysis were applied to explore β-diversity of carbon-fixing microbes. PCA biplot showed that 35.93% of total variance in carbon-fixing microbial community structure was explained by the first two axes (Fig. 4). It also showed that five treatments were clearly differentiated into three clusters, UL and CK as cluster 1, M and MNPK as cluster 2, and NPK as cluster 3 along the PC2 axis, which indicated that the carbon-fixing microbial community structure among the three clusters was different. The score in PC2 axis significantly differed among each cluster (P < 0.05).
Figure 4.
Principal component analysis (PCA) of the carbon-fixing microbial community in different treatments.
Biclustering heatmap were also applied to explore the community structure of carbon-fixing microbes under different treatments (Fig. 5). Vertical clustering of heatmap showed that five treatments were classified into two clusters: M and MNPK were clustered together, while UL, CK and NPK formed a second cluster, and the second cluster was divided into two subclusters, UL and CK as a subcluster, NPK as another subcluster, which was consistent with PCA. Horizontal clustering of heatmap showed that the carbon-fixing microbial community structure differed among treatments. Compared with UL and CK, Actinomadura was increased in NPK, M and MNPK. Compared with the other three treatments Variovorax, Hydrogenophaga, Novosphingobium and Polumorphum harbored a higher relative abundance in M and MNPK. Pseudonocardia relative abundance was higher in NPK and MNPK than in three other treatments. Compared with four other treatments, Mesorhizobium, Nitrobacter, Cupriavidus, Rhodoferax, Rhodoblastus, Mycobacterium, Sphingomonas and Sinorhizobium were more abundant in MNPK. The relative abundance of Noviherbaspirillum, Marichromatium and Thermoleptolyngbya were higher in UL than in four other treatments.
Figure 5.
Biclustering heatmap of the carbon-fixing microbial distributions of the top 30 abundant genera is present in different treatments. The color intensity of the color lumps represents the abundance of the carbon-fixing microbial genera in different treatments, with red representing higher abundance and blue representing lower abundance.
Correlations between the carbon-fixing microbial communities and soil properties
-
RDA were utilized to investigate the correlation of soil properties with soil carbon-fixing microbial community structure. Soil properties with VIF > 5 were removed before RDA to avoid the effect of collinearity. Then, eight variables, including CAT, EOC, AP, DOC, AK, SM, TP and AN were screened out for RDA.
The RDA (Fig. 6a) showed that the first two axes explained 41.93% and 11.45% of the total variance in soil carbon-fixing microbial community, respectively. The score in RDA2 axis significantly differed among treatments (P < 0.05). M and MNPK were clustered together and obviously separated from UL, CK and NPK, which was consistent with PCA and biclustering heatmap results. The soil properties with the highest explanatory proportion of soil carbon-fixing microbial community structure were SM (r = 0.6288, P = 0.001), CAT (r = 0.5508, P = 0.006), TP (r = 3972, P = 0.048), DOC (r = 0.3875, P = 0.041).
Figure 6.
(a) Redundancy analysis (RDA) linking carbon-fixing microbial communities with environmental variables in different treatments. Arrows represent the correlation between the soil properties and carbon-fixing microbial communities. Variables that are angled at more than 90° of each other have the least correlation. The length of the arrow represents the correlation. Variables that have arrows extended in opposite directions correlate negatively to each other. (b) Diagrams explaining variance partitioning (VPA) show the relative contribution of ecological drivers with VIF < 5 to soil carbon-fixing microbial community structure. The abiotic variables include EOC, AP, DOC, AK, SM, TP and AN; the biotic variables include CAT. The numbers are the percentages of the total variables explained by the factors.
To further unveil the relationship between carbon-fixing microbial community and environmental variables, the explanatory variables included in RDA were divided into two groups (abiotic variables and biotic variables) for VPA. The VPA (Fig. 6b) showed that abiotic variables were the more important factors affecting carbon-fixing microbial community assembly. The abiotic variables and biotic variables explained 22.38% and 0.15% of the variation in the carbon-fixing microbial community, respectively. Moreover, 64.98% of variation remained unexplained.
Based on RDA and PCA, the dominant soil properties (CAT, DOC, TP and EOC) and fertilization treatments with different manure application doses that affected the structure of carbon-fixing microbial community, were further included in SEM analysis. SEM was applied to reveal the causal relationship among soil properties, treatments (manure application dose), RubisCO activity, community composition (loading score on the first PCA axis) (Li et al., 2015), cbbL gene abundance, α diversity (Shannon-Weaver diversity index) (Zhang et al., 2021) and soil carbon accumulation (Fig. 7). It also used to further explain and quantify the contribution of key factors influencing both carbon-fixing microbial community and soil carbon fixation.
Figure 7.
Structural equation model (SEM) shows the causal influences of treatments, DOC, EOC, CAT, TP, RubisCO activity, α diversity of cbbL, community composition (carbon-fixing microbial community), cbbL gene abundance, SOC and MBC. Positive and negative effects are respectively showed in red and green, and significant and non-significant effects are showed with solid and dashed arrow lines, respectively. The standardized coefficients are marked above each path (only marked significant effect paths) and indicate the expected impact of a unit standard-deviation change at one node on units of standard-deviation change in connected nodes. R2 values represent the proportion of the variance explained for each endogenous variable.
This model is completely consistent with our causal hypothesis (χ2 = 13.060, df = 23, P = 0.951, GFI = 0.892, RMSEA = 0.000, AIC = 99.060) and it could explain 99%, 96%, 94%, 57% and 43% variance of the SOC, MBC, community composition, RubisCO activity and α diversity of cbbL, respectively. SOC was directly affected by treatments, MBC, TP, CAT, RubisCO activity, α diversity of cbbL and community composition and indirectly affected by DOC (standardized indirect effects = −0.15) and EOC (standardized indirect effects = 0.06). MBC was directly affected by treatments, CAT, DOC, RubisCO activity, community composition and cbbL gene abundance. Community composition was directly affected by treatments, EOC, DOC and α diversity of cbbL.
-
Our study revealed that reclamation even without fertilization in the initial stage could bring significant shifts in MBC and dominant taxa (i.e., Proteobacteria, Cyanobacteria, Devosia and Marichromatium) and showed a trend of ecological recovery. Moreover, reclamation with fertilization significantly increased SOC and MBC (P < 0.05) and significantly altered carbon-fixing microbial community composition. Among these fertilization treatments, the application of manure (M) is more conducive to increasing carbon-fixing microbial abundance, the cbbL/16S rRNA ratio and RubisCO activity in the current short-term reclamation. SM, CAT, TP and DOC were the key factors significant influencing soil carbon-fixing microbial community structure, which provided a theoretical basis for improving the carbon-fixing potential. Reclamation and fertilization could significantly influence carbon-fixing microbial community structure (P < 0.05) and increase soil carbon storage due to altered soil properties and manure application dose. The fertilization treatments with manure were more conducive to improving the soil carbon-fixing potential. These findings contributed to improving soil fertility and accelerating ecological restoration and reconstruction in the mining area.
-
About this article
Cite this article
Shang Y, Wu M, Zhang J, Meng H, Hong J, et al. 2023. Nutrient enhanced reclamation promoted growth, diversity and activities of carbon fixing microbes in a coal-mining subsistence land. Soil Science and Environment 2:2 doi: 10.48130/SSE-2023-0002
Nutrient enhanced reclamation promoted growth, diversity and activities of carbon fixing microbes in a coal-mining subsistence land
- Received: 02 December 2022
- Accepted: 03 March 2023
- Published online: 03 April 2023
Abstract: Carbon-fixing microbes can potentially improve soil fertility. However, the potential and function of carbon-fixing microbes remains largely uninvestigated in reclaimed soil of coal-mining subsidence areas. In this study, treatments included UL (uncultivated land), CK (maize cultivation without fertilization), NPK (maize cultivation with chemical fertilizer), M (maize cultivation with manure), MNPK (maize cultivation with manure and chemical fertilizer) after 1-year reclamation in a typical coal mining subsidence area. Quantitative PCR, enzyme-linked immunosorbent assay (ELISA) and high-throughput sequencing were employed to investigate the topsoil carbon-fixing microbial biomass, RubisCO activity and community composition. The results showed that the dominant taxa (i.e., Proteobacteria, Cyanobacteria, Devosia and Marichromatium) were significantly changed after reclamation (P < 0.05). Carbon-fixing microbial community structure in fertilization treatments (NPK, M and MNPK) obviously differed from non-fertilizer treatments (UL and CK). Soil organic carbon and microbial biomass carbon were significantly higher in fertilization treatments than non-fertilizer treatments (P < 0.05). M significantly increased RubisCO activity and cbbL gene abundance (P < 0.05), MNPK significantly increased carbon-fixing microbial richness (P < 0.05). Carbon-fixing microbial community structure was strongly influenced by soil moisture, catalase, total phosphorus and dissolved organic carbon. Some environmental factors indirectly influenced SOC by affecting carbon-fixing microbial biomass, diversity and community structure. Our study implies that even short-term (1-year) reclamation and fertilization could significantly influence carbon-fixing microbial community structure and promote soil carbon accumulation, and the fertilization treatments with manure (M and MNPK) were more conducive, which indicated that carbon-fixing microbes were greatly conducive to improve soil fertility in reclaimed mining areas and achieve carbon neutrality.
-
Key words:
- coal-mining /
- reclaimed soil /
- fertilization /
- cbbL gene /
- carbon-fixing microbes /
- RubisCO activity