Interplay of Darwinian and frequency-dependent selection in the host-associated microbial populations

In order to analyze the microevolutionary processes in host-associated microorganisms, we simulated the dynamics of rhizobia populations composed of a parental strain and its mutants possessing the altered fitness within “plant-soil” system. The population dynamics was presented as a series of cycles (each one involves “soil→rhizosphere→nodules→soil” succession) described using recurrent equations. For representing the selection and mutation pressures, we used a universal approach based on calculating the shifts in the genetic ratios of competing bacterial genotypes within the particular habitats and across several habitats. Analysis of the model demonstrated that a balanced polymorphism may be established in rhizobia population: mutants with an improved fitness do not supplant completely the parental strain while mutants with a decreased fitness may be maintained stably. This polymorphism is caused by a rescue of low-fitted genotypes via negative frequency-dependent selection (FDS) that is implemented during inoculation of nodules and balances the Darwinian selection that occurs during multiplication or extinction of bacteria at different habitats. The most diverse populations are formed if the rhizobia are equally successful in soil and nodules, while a marked preference for any of these habitats results in the decrease of diversity. Our simulation suggests that FDS can maintain the mutualistic rhizobia–legume interactions under the stress conditions deleterious for surviving the bacterial strains capable for intensive N2 fixation. Genetic consequences of releasing the modified rhizobia strains may be addressed using the presented model.