Generalist bacterial predators are likely to strongly shape many important ecological and evolutionary features of microbial communities, for example by altering the character and pace of molecular evolution, but investigations of such effects are scarce. Here we report how predator-prey interactions alter the evolution of fitness, genomes and phenotypic diversity in coevolving bacterial communities composed of Myxococcus xanthus as predator and Escherichia coli as prey, relative to single-species controls. We show evidence of reciprocal adaptation and demonstrate accelerated genomic evolution specific to coevolving communities, including the rapid appearance of mutator genotypes. Strong parallel evolution unique to the predator-prey communities occurs in both parties, with predators driving adaptation at two prey traits associated with virulence in bacterial pathogens—mucoidy and the outer-membrane protease OmpT. Our results suggest that generalist predatory bacteria are important determinants of how complex microbial communities and their interaction networks evolve in natural habitats.

Microbial genotypes with similarly high proficiency at a cooperative behaviour in genetically pure groups often exhibit fitness inequalities caused by social interaction in mixed groups. Winning competitors in this scenario have been referred to as ‘cheaters’ in some studies. Such interaction-specific fitness inequalities, as well as social exploitation (in which interaction between genotypes increases absolute fitness), might evolve due to selection for competitiveness at the focal behaviour or might arise non-adaptively due to pleiotropy, hitchhiking or genetic drift. The bacterium Myxococcus xanthus sporulates during cooperative development of multicellular fruiting bodies. Using M. xanthus lineages that underwent experimental evolution in allopatry without selection on sporulation, we demonstrate that interaction-specific fitness inequalities and facultative social exploitation during development readily evolved indirectly among descendant lineages. Fitness inequalities between evolved genotypes were not caused by divergence in developmental speed, as faster-developing strains were not over-represented among competition winners. In competitions between ancestors and several evolved strains, all evolved genotypes produced more spores than the ancestors, including losers of evolved-versus-evolved competitions, indicating that adaptation in non-developmental contexts pleiotropically increased competitiveness for spore production. Overall, our results suggest that fitness inequalities caused by social interaction during cooperative processes may often evolve non-adaptively in natural populations.

Predator impacts on prey diversity are often studied among higher organisms over short periods, but microbial predator-prey systems allow examination of prey-diversity dynamics over evolutionary timescales. We previously showed that Escherichia coli commonly evolved minority mucoid phenotypes in response to predation by the bacterial predator Myxococcus xanthus by one time point of a coevolution experiment now named MyxoEE-6. Here we examine mucoid frequencies across several MyxoEE-6 timepoints to discriminate between the hypotheses that mucoids were increasing to fixation, stabilizing around equilibrium frequencies, or heading to loss toward the end of MyxoEE-6. In four focal coevolved prey populations, mucoids rose rapidly early in the experiment and then fluctuated within detectable minority frequency ranges through the end of MyxoEE-6, generating frequency dynamics suggestive of negative frequency-dependent selection. However, a competition experiment between mucoid and non-mucoid clones found a predation-specific advantage of the mucoid clone that was insensitive to frequency over the examined range, leaving the mechanism that maintains minority mucoidy unresolved. The advantage of mucoidy under predation was found to be associated with reduced population size after growth (productivity) in the absence of predators, suggesting a tradeoff between productivity and resistance to predation that we hypothesize may reverse mucoid vs non-mucoid fitness ranks within each MyxoEE-6 cycle. We also found that mucoidy was associated with diverse colony phenotypes and diverse candidate mutations primarily localized in the exopolysaccharide operon yjbEFGH. Collectively, our results show that selection from predatory bacteria can generate apparently stable sympatric phenotypic polymorphisms within coevolving prey populations and also allopatric diversity across populations by selecting for diverse mutations and colony phenotypes associated with mucoidy. More broadly, our results suggest that myxobacterial predation increases long-term diversity within natural microbial communities.