Abstract
Understanding interspecies interactions is essential to predict the response of microbial communities to exogenous perturbation. Herein, rhizospheric and bulk soils were collected from five developmental stages of soybean, which grew in soils receiving 16-year nitrogen inputs. Bacterial communities and functional profiles were examined using high-throughput sequencing and quantitative PCR, respectively. The objective of this study was to identify the key bacterial interactions that influenced community dynamics and functions. We found that the stages of soybean development outcompeted nitrogen fertilization management in shaping bacterial community structure, while fertilization treatments significantly shaped the abundance distribution of nitrogen functional genes. Temporal variations in bacterial abundances increased in bulk soils, especially at the stage of soybean branching, which helps to infer underlying negative interspecies interactions. Members of Cyanobacteria and Actinobacteria actively engaged in inter-phylum negative interactions in bulk soils and soybean rhizosphere, respectively. Furthermore, the negative interactions between nitrogen-fixing functional groups and the reduction of nifH gene abundance were coupled during soybean development, which may help to explain the linkages between population dynamics and functions. Overall, these findings highlight the importance of inter-phylum negative interactions in shaping the correlation patterns of bacterial communities and in determining soil functional potential.