Abstract
Host-associated resident microbiota can protect their host from pathogens-a community-level trait called colonization resistance. The effect of the diversity of the resident community in previous studies has shown contradictory results, with higher diversity either strengthening or weakening colonization resistance. To control the confounding factors that may lead to such contradictions, we use mathematical simulations with a focus on species interactions and their impact on colonization resistance. We use a mediator-explicit model that accounts for metabolite-mediated interactions to perform in silico invasion experiments. We show that the relationship between colonization resistance and species richness of the resident community is not monotonic because it depends on two underlying trends as the richness of the resident community increases: a decrease in instances of augmentation (invader species added, without driving out resident species) and an increase in instances of displacement (invader species added, driving out some of the resident species). These trends hold consistently under different parameters, regardless of the number of compounds that mediate interactions between species or the proportion of the facilitative versus inhibitory interactions among species. Our results show a positive correlation between resistance and diversity in low-richness communities and a negative correlation in high-richness communities, offering an explanation for the seemingly contradictory trend in the resistance-diversity relationship in previous reports.