Regularized S-Map Reveals Varying Bacterial Interactions

There is a growing awareness that bacterial interactions follow a highly nonlinear pattern in reality. However, it is challenging to track the varying bacterial interactions using pairwise correlation analysis, which fails to explore their potential effects on the behavior of microbes. Here, we utilized a regularized sequential locally weighted global linear map (S-map) to capture the varying interspecific interactions from the time series data of a bacterial community under exposure to nitrite. Our results show that bacterial interactions are highly variable and that asymmetric interactions dominate the interaction pattern in a community. Furthermore, we propose a Jacobian coefficient-based statistical method to predict the harmony level of a bacterial community at each successive ecosystem state. The results show that the bacterial community exhibits a higher harmony level in nitrite-treated samples than in the control group. We show that the community harmony level is positively associated with the specific endogenous respiration rates and biofilm formation of the culture. In addition, the community tends to process lower diversity and structural stability under zero- and high-nitrite stresses. We demonstrate that the harmony level, rather than structural stability, is a useful index for unveiling the underlying mechanism of bacterial performance. Overall, the regularized S-map can help us to understand bacterial interactions in ecosystems more accurately than previous approaches.IMPORTANCE It has long been acknowledged that bacterial interactions play important roles in community structure and function. Revealing the interaction variability can allow an understanding of how bacteria respond to perturbation and why bacterial community performance changes. Such information should improve our skills in engineering bacterial communities (e.g., in a wastewater treatment plant) and achieve better removal performance and lower energy consumption.