MCR-1 is a lipid A modifying enzyme that confers resistance to the antibiotic colistin. Here, we analyse the impact of MCR-1 expression on E. coli morphology, fitness, competitiveness, immune stimulation and virulence. Increased expression of mcr-1 results in decreased growth rate, cell viability, competitive ability and significant degradation in cell membrane and cytoplasmic structures, compared to expression of catalytically inactive MCR-1 (E246A) or MCR-1 soluble component. Lipopolysaccharide (LPS) extracted from mcr-1 strains induces lower production of IL-6 and TNF, when compared to control LPS. Compared to their parent strains, high-level colistin resistance mutants (HLCRMs) show reduced fitness (relative fitness is 0.41–0.78) and highly attenuated virulence in a Galleria mellonella infection model. Furthermore, HLCRMs are more susceptible to most antibiotics than their respective parent strains. Our results show that the bacterium is challenged to find a delicate equilibrium between expression of MCR-1-mediated colistin resistance and minimalizing toxicity and thus ensuring cell survival.

SignificancePathogens are omnipresent and by definition detrimental to their hosts. Pathogens thus exert high selection on their hosts, which, if adapting, can exert similar levels of selection on the pathogen, resulting in ongoing cycles of reciprocal adaptation between the antagonists. Such coevolutionary interactions have a central influence on the evolution of organisms. Surprisingly, we still know little about the exact selection dynamics and the genome regions involved. Our study uses a controlled experimental approach with an animal host to dissect coevolutionary selection. We find that distinct selective processes underlie rapid coadaptation in the two antagonists, including antagonistic frequency-dependent selection on toxin gene copy number in the pathogen, while the host response is likely influenced by changes in multiple genome regions.

BackgroundHost-parasite coevolution is generally believed to follow Red Queen dynamics consisting of ongoing oscillations in the frequencies of interacting host and parasite alleles. This belief is founded on previous theoretical work, which assumes infinite or constant population size. To what extent are such sustained oscillations realistic?ResultsHere, we use a related mathematical modeling approach to demonstrate that ongoing Red Queen dynamics is unlikely. In fact, they collapse rapidly when two critical pieces of realism are acknowledged: (i) population size fluctuations, caused by the antagonism of the interaction in concordance with the Lotka-Volterra relationship; and (ii) stochasticity, acting in any finite population. Together, these two factors cause fast allele fixation. Fixation is not restricted to common alleles, as expected from drift, but also seen for originally rare alleles under a wide parameter space, potentially facilitating spread of novel variants.ConclusionOur results call for a paradigm shift in our understanding of host-parasite coevolution, strongly suggesting that these are driven by recurrent selective sweeps rather than continuous allele oscillations.

Reciprocal coevolution between host and pathogen is widely seen as a major driver of evolution and biological innovation. Yet, to date, the underlying genetic mechanisms and associated trait functions that are unique to rapid coevolutionary change are generally unknown. We here combined experimental evolution of the bacterial biocontrol agent Bacillus thuringiensis and its nematode host Caenorhabditis elegans with large-scale phenotyping, whole genome analysis, and functional genetics to demonstrate the selective benefit of pathogen virulence and the underlying toxin genes during the adaptation process. We show that: (i) high virulence was specifically favoured during pathogen–host coevolution rather than pathogen one-sided adaptation to a nonchanging host or to an environment without host; (ii) the pathogen genotype BT-679 with known nematocidal toxin genes and high virulence specifically swept to fixation in all of the independent replicate populations under coevolution but only some under one-sided adaptation; (iii) high virulence in the BT-679-dominated populations correlated with elevated copy numbers of the plasmid containing the nematocidal toxin genes; (iv) loss of virulence in a toxin-plasmid lacking BT-679 isolate was reconstituted by genetic reintroduction or external addition of the toxins. We conclude that sustained coevolution is distinct from unidirectional selection in shaping the pathogen's genome and life history characteristics. To our knowledge, this study is the first to characterize the pathogen genes involved in coevolutionary adaptation in an animal host–pathogen interaction system.