In microbial communities, the ecological interactions between species of different populations are responsible for the spatial distributions observed in aggregates (granules, biofilms or flocs). To explore the underlying mechanisms that control these processes, we have developed a mathematical modelling framework able to describe, label and quantify defined spatial structures that arise from microbial and environmental interactions in communities. An artificial system of three populations collaborating or competing in an aggregate is simulated using individual-based modelling under different environmental conditions. In this study, neutralism, competition, commensalism and concurrence of commensalism and competition have been considered. We were able to identify interspecific segregation of communities that appears in competitive environments (columned stratification), and a layered distribution of populations that emerges in commensal (layered stratification). When different ecological interactions were considered in the same aggregate, the resultant spatial distribution was identified as the one controlled by the most limiting substrate. A theoretical modulus was defined, with which we were able to quantify the effect of environmental conditions and ecological interactions to predict the most probable spatial distribution. The specific microbial patterns observed in our results allowed us to identify the optimal spatial organizations for bacteria to thrive when building a microbial community and how this permitted co-existence of populations at different growth rates. Our model reveals that although ecological relationships between different species dictate the distribution of bacteria, the environment controls the final spatial distribution of the community.

Is it possible to find trends between the parameters that define microbial growth to help us explain the vast microbial diversity? Through an extensive database of kinetic parameters of nitrifiers, we analyzed if the dominance of specific populations of nitrifiers could be predicted and explained. We concluded that, in general, higher growth yield (Y_XS) and ammonia affinity (a⁰_NH3) and lower growth rate (µ_max) are observed for ammonia-oxidizing archaea (AOA) than bacteria (AOB), which would explain their considered dominance in oligotrophic environments. However, comammox (CMX), with the maximum energy harvest per mole of ammonia, and some AOB, have higher a⁰_NH3 and lower µ_max than some AOA. Although we were able to correlate the presence of specific terminal oxidases with observed oxygen affinities (a⁰_O2) for nitrite-oxidizing bacteria (NOB), that correlation was not observed for AOB. Moreover, the presumed dominance of AOB over NOB in O₂-limiting environments is discussed. Additionally, lower statistical variance of a⁰_O2 values than for ammonia and nitrite affinities was observed, suggesting nitrogen limitation as a stronger selective pressure. Overall, specific growth strategies within nitrifying groups were not identified through the reported kinetic parameters, which might suggest that mostly, fundamental differences in biochemistry are responsible for underlying kinetic parameters.