Cooperation between autotrophic and heterotrophic denitrifiers under low C/N ratios revealed by individual-based modelling

Denitrifying biofilms, in which autotrophic denitrifiers (AD) and heterotrophic denitrifiers (HD) coexist, play a crucial role in removing nitrate from water or wastewater. However, it is difficult to elucidate the interactions between HD and AD through sequencing-based experimental methods. Here, we developed an individual-based model to describe the interspecies dynamics and priority effects between sulfur-based AD (Thiobacillus denitrificans) and HD (Thauera phenylcarboxya) under different C/N ratios. In test I (coexistence simulation), AD and HD were initially inoculated at a ratio of 1:1. The simulation results showed excellent denitrification performance and a coaggregation pattern of denitrifiers, indicating that cooperation was the predominant interaction at a C/N ratio of 0.25 to 1.5. In test II (invasion simulation), in which only one type of denitrifier was initially inoculated and the other was added at the invasion time, denitrifiers exhibited a stratification pattern in biofilms. When HD invaded AD, the final HD abundance decreased with increasing invasion time, indicating an enhanced priority effect. When AD invaded HD, insufficient organic carbon sources weakened the priority effect by limiting the growth of HD populations. This study reveals the interaction between autotrophic and heterotrophic denitrifiers, providing guidance for optimizing wastewater treatment process.