Soil pH and carbon quality index regulate the biogeochemical cycle couplings of carbon, nitrogen and phosphorus in the profiles of Isohumosols.

Soil biogeochemical cycles are essential for regulating ecosystem functions and services. However, little knowledge has been revealed on microbe-driven biogeochemical processes and their coupling mechanisms in soil profiles. This study investigated the vertical distribution of soil functional composition and their contribution to carbon (C), nitrogen (N) and phosphorus (P) cycling in the humus horizons (A-horizons) and parent material horizons (C-horizons) in Udic and Ustic Isohumosols using shotgun sequencing. Results showed that the diversity and relative abundance of microbial functional genes was influenced by soil horizons and soil types. In A-horizons, the relative abundances of N mineralization and liable C decomposition genes were significantly greater, but the P cycle-related genes, recalcitrant C decomposition and denitrification genes were lower compared to C-horizons. While, Ustic Isohumosols had lower relative abundances of C decomposition genes but higher relative abundances of N mineralization and P cycling-related pathways compared to Udic Isohumosols. The network analysis revealed that C-horizons had more interactions and stronger stability of functional gene networks than in A-horizons. Importantly, our results provide new insights into the potential mechanisms for the coupling processes of soil biogeochemical cycles among C, N and P, which is mediated by specific microbial taxa. Soil pH and carbon quality index (CQI) were two sensitive indicators for regulating the relative abundances and the relationships of functional genes in biogeochemical cycles. This study contributes to a deeper understanding of the ecological functions of soil microorganisms, thus providing a theoretical basis for the exploration and utilization of soil microbial resources and the development of soil ecological control strategies.