Bacterial Communities, Network Complexity, and Multifunctionality Affected by Soil Types in Northeastern China

The Northeast China Plain (NCP) is the country’s most important grain-producing area. Unraveling how bacterial communities in this region assemble and distribute according to soil type is essential for sustainable agricultural development and optimizing the precise management of soil resources. In this study, 106 soil samples were collected from three typical zonal soil types (black calcium soil (BCS), black soil (BS), and dark brown soil (DBS)) spanning from west to east in the NCP. By combining soil field surveys and high-throughput microbial sequencing analysis, we found that bacterial diversity and community structure differed significantly by soil type. Proteobacteria, Gemmatimonadetes, and Acidobacteria were enriched in BCS, BS, and DBS, respectively. Compared to BSC and DBS, BS had the highest nutrient concentration and most neutral pH values, which may recruit more diverse bacterial communities and construct a more connected ecological network. Network analysis further identified Burkholderiales, Sphingomonadales, and SC_I_84 as key hubs in BS, BCS, and BCS, respectively. The majority of classified hubs consistent with the results of the linear discriminant analysis effect size belonged to the predominant biomarkers. Redundancy and Mantel test analyses revealed that the bacterial composition in various soil types showed distinctive responses to heterogeneity in soil physicochemical properties. Soil pH and TP were the primary factors shaping the soil bacterial community structure in these three soil types on the NCP. Moreover, bacterial composition and diversity were strongly related to changes in soil multifunctionality in BCS, and the relative abundances of three classes (TM1, Opitutae, and Deinococci) were the most important biotic variables for predicting BCS ecosystem multifunctionality. In summary, our results suggest that soil type variation has a strong influence in terms of shaping bacterial community structure and affecting soil multifunctionality. Correspondingly, diverse co-occurrence patterns were observed in different soil types.