Site A and site B did not have significantly different NO3−, soil moisture, or pH (Wilcoxon test: P = 0.898, P = 0.530, P = 0.654, respectively, Table 1). NH4+ was significantly and on average 122% higher in site A than site B (Wilcoxon test: P = 0.029; Table 1). Concentrations of NH4+ and NO3− did not differ across vegetation types or depths (Kruskal-Wallis test for vegetation: Chi-squared = 1.754, P = 0.416 for NH4+, Chi-squared = 2.799, P = 0.247 for NO3−; Wilcoxon test for depths: P = 0.457 for NH4+, P = 0.380 for NO3−; Table 1 and Figure 3A, B). Soil moisture was significantly different between depths (Wilcoxon test: P < 0.0001) but not vegetation types (Kruskal-Wallis test: Chi-squared = 0.692, P = 0.708, Table 1). Soil pH was significantly different across both vegetation types and depth (Kruskal-Wallis test for vegetation: chi-squared = 10.672, P = 0.005; Wilcoxon test for depth: P < 0.001). The post-hoc Dunn’s test showed that BG and CG had significantly higher pH than SB (P = 0.041 and 0.005, respectively).
Table 1
Comparison of soil chemical properties between sites and vegetation.
Site | Coordinates | Vegetation | Moisture content (%) | pH | NO3− (µg/ g dry soil) | NH4+ (µg/ g dry soil) |
| N 40.180747 W 112.452441 | Bunchgrass | 8.17 ± 1.20 | 9.02 ± 0.21 | 0.84 ± 0.63 | 1.08 ± 1.16 |
Site A | Cheatgrass | 8.35 ± 0.71 | 9.60 ± 0.28 | 0.00 ± 0.00 | 1.08 ± 1.16 |
| Sagebrush | 8.29 ± 1.01 | 8.83 ± 0.19 | 0.08 ± 0.20 | 1.06 ± 0.98 |
| N 40.180445 W 112.451753 | Bunchgrass | 7.79 ± 1.01 | 9.16 ± 0.25 | 0.60 ± 0.43 | 0.52 ± 0.37 |
Site B | Cheatgrass | 8.59 ± 1.05 | 9.43 ± 0.37 | 0.10 ± 0.24 | 0.23 ± 0.14 |
| Sagebrush | 7.85 ± 0.70 | 8.90 ± 0.21 | 0.22 ± 0.37 | 0.37 ± 0.21 |
The total richness (i.e., number of unique ASVs) was 9,755 and 2,045 for 16S and ITS, respectively. The average richness per soil sample was 1,540 (sd = 306) for 16S and 213 (sd = 28) for ITS, after rarefaction (Figure 2A, B). At the phylum level, the soil prokaryotic communities (Supplementary figure 3) were dominated by Actinobacteria (47.4%), Proteobacteria (19.4%), Chloroflexi (7.2%), Thaumarchaeota (6.1%), and Acidobacteria (5.2%). The fungal communities (Supplementary figure 4) were dominated by Ascomycota (72.1%), Basidiomycota (11.7%), Mortierellomycota (9.8%), Olpidiomycota (3.5%), and Glomeromycota (2.0%).
Since microbial diversity metrics did not significantly differ between both sites (t-test: P = 0.549 for 16S richness, P = 0.840 for 16S Shannon, P = 0.409 for ITS richness, P = 0.384 for ITS Shannon), samples from two sites were not analyzed separately. Prokaryotic richness showed no significant differences among vegetation or depths (ANOVA: F = 0.9442, P = 0.401 for vegetation, F = 0.047, P = 0.829 for depth; Figure 2A), and the same pattern were observed for fungi (ANOVA: F = 1.635, P = 0.212 for vegetation, F = 0.010, P = 0.928 for depth; Figure 2C). Shannon diversity showed a similar nonsignificant pattern (ANOVA for 16S: F = 0.159, P = 0.854 for vegetation, F = 0.220, P = 0.642 for depth; ANOVA for ITS: F = 1.279, P = 0.294 for vegetation, F = 1.341, P = 0.255 for depth; Supplementary Figure 1).
Prokaryotic community composition was significantly different between vegetation types and soil depths (PERMANOVA: R2 = 0.143, P < 0.001 for vegetation, R2 = 0.144, P < 0.001 for depth; Figure 2B; Supplementary Figures 3 and 4). Fungal communities were only significantly different across vegetation types but not depths (PERMANOVA: R2 = 0.145, P < 0.001 for vegetation, R2 = 0.029, P = 0.296 for depth; Figure 2D). Soil moisture and pH were the most important predictors of the composition of prokaryotic communities (Figure 1B). Fungal community composition was mainly determined by pH (Figure 1D).
Among the eight nitrogen cycle-related bacterial groups, the relative abundance of four functional groups (nitrogen fixers, denitrifiers, nitrogen respiratory bacteria, and ureolytic bacteria) were significantly different among vegetation types (Kruskal-Wallis test; Figure 3). Nitrogen fixers (mainly Herbaspirillum spp., Bradyrhizobium spp., and Nostoc spp.) and ureolytic bacteria (mainly Roseomonas spp., Mesorhizobium spp., Singulisphaera spp., Massilia spp., and Methylobacterium spp.) had higher relative abundance in CG samples (Dunn’s test for nitrogen fixers: P = 0.037 between CG and BG; Dunn’s test for ureolytic bacteria: P = 0.005 between CG and SB). Nitrogen respiratory bacteria (Rhodoplanes spp., and Nitrobacter spp) and denitrifiers (Pseudomonas fluorescens, and Rhodoplanes spp) had lower relative abundance in CG (Dunn’s test for nitrogen respiratory bacteria: P = 0.012 between CG and SB; Dunn’s test for denitrifiers: P = 0.021 between CG and SB). After FDR correction for multiple testing, only ureolytic bacteria was significantly different across vegetation. For fungi, the relative abundances of dung saprotrophs and leaf saprotrophs were significantly different among vegetation types after FDR correction (Kruskal-Wallis test; Supplementary Figure 2).