Background: Transforming coniferous plantation into broadleaved or mixed broadleaved-coniferous plantations is the tendency of forest management strategies in subtropical China. However, the effects of this conversion on soil phosphorus (P) cycling microbial functional genes are still unknown.Methods: Soil samples were collected from 0–20, 20–40, and 40–60 cm (topsoil, middle layer, and subsoil, respectively) under coniferous Pinus massoniana (PM), broadleaved Erythrophleum fordii (EF), and their mixed (PM/EF) plantation in subtropical China. Used metagenomic sequencing to examine the alterations of relative abundances and molecular ecological network structure of soil P-cycling functional genes after the conversion of plantations.Results: The composition of P-cycling genes in the topsoil of PM stand was significantly different from that of PM/EF and EF stands (p < 0.05), and total phosphorus (TP) was the main factor causing this difference. After transforming PM plantation into EF plantation, the relative abundances of P solubilization and mineralization genes significantly increased in the topsoil and middle layer with the decrease of soil TP content. The abundances of P-starvation response regulation genes also significantly increased in the subsoil (p < 0.05), which may have been influenced by soil organic carbon (SOC). The dominant genes in all soil layers under three plantations were phoR, glpP, gcd, ppk, and ppx. Transforming PM into EF plantation apparently increased gcd abundance in the topsoil (p < 0.05), with TP and NO₃⁻-N being the main influencing factors. After transforming PM into PM/EF plantations, the molecular ecological network structure of P-cycling genes was more complex; moreover, the key genes in the network were modified with the transformation of PM plantation.Conclusion: Transforming PM into EF plantation mainly improved the phosphate solubilizing potential of microorganisms at topsoil, while transforming PM into PM/EF plantation may have enhanced structural stability of microbial P-cycling genes react to environmental changes.

The soil profile is a strong and complex physicochemical gradient that greatly affects bacterial community structure and function between soil layers. However, little is known about molecular ecological network structure and bacterial community function under differing soil profiles in planted forests. Four typical native tree species (Pinus massoniana Lamb., Castanopsis hystrix Miq., Mytilaria laosensis Lec., and Michelia macclurei Dandy) plantations were selected from subtropical China as the research object. We evaluated molecular ecological network structure as well as potential function of the soil bacterial community at different soil depths (0–20, 20–40, and 40–60 cm) within native tree plantations. Our results showed that (1) compared to the topsoil (0–20 cm), the bacterial molecular ecological network scale increased within the middle layer (20–40 cm) and the subsoil (40–60 cm), and the interaction between species was stronger; (2) module hubs and connectors were the key bacterial groups in each soil layer and increased with increasing soil depth; (3) the dominant functional groups of the bacterial communities in each soil layer were chemoheterotrophy, aerobic chemoheterotrophy, cellulolysis, ureolysis, nitrogen fixation, and nitrate reduction, and they were related to soil carbon and nitrogen cycling; and (4) the different molecular ecological network structures along with relative bacterial functional group abundances among diverse soil layers were mainly affected by soil organic carbon (SOC), NO₃⁻-N, NH₄⁺-N, available phosphorus (AP), and total phosphorus (TP). Our study provides a theoretical foundation for bacterial community structure together with function within soil profiles of native tree plantations in subtropical regions.