Hypermutator strains of Pseudomonas aeruginosa reveal novel pathways of resistance to combinations of cephalosporin antibiotics and beta-lactamase inhibitors.

Hypermutation due to DNA mismatch repair (MMR) deficiencies can accelerate the development of antibiotic resistance in Pseudomonas aeruginosa. Whether hypermutators generate resistance through predominantly similar molecular mechanisms to wild-type (WT) strains is not fully understood. Here, we show that MMR-deficient P. aeruginosa can evolve resistance to important broad-spectrum cephalosporin/beta-lactamase inhibitor combination antibiotics through novel mechanisms not commonly observed in WT lineages. Using whole-genome sequencing (WGS) and transcriptional profiling of isolates that underwent in vitro adaptation to ceftazidime/avibactam (CZA), we characterized the detailed sequence of mutational and transcriptional changes underlying the development of resistance. Surprisingly, MMR-deficient lineages rapidly developed high-level resistance (>256 µg/mL) largely without corresponding fixed mutations or transcriptional changes in well-established resistance genes. Further investigation revealed that these isolates had paradoxically generated an early inactivating mutation in the mexB gene of the MexAB-OprM efflux pump, a primary mediator of CZA resistance in P. aeruginosa, potentially driving an evolutionary search for alternative resistance mechanisms. In addition to alterations in a number of genes not known to be associated with resistance, 2 mutations were observed in the operon encoding the RND efflux pump MexVW. These mutations resulted in a 4- to 6-fold increase in resistance to ceftazidime, CZA, cefepime, and ceftolozane-tazobactam when engineered into a WT strain, demonstrating a potentially important and previously unappreciated mechanism of resistance to these antibiotics in P. aeruginosa. Our results suggest that MMR-deficient isolates may rapidly evolve novel resistance mechanisms, sometimes with complex dynamics that reflect gene inactivation that occurs with hypermutation. The apparent ease with which hypermutators may switch to alternative resistance mechanisms for which antibiotics have not been developed may carry important clinical implications.

Escherichia coli Genomic Diversity within Extraintestinal Acute Infections Argues for Adaptive Evolution at Play.

Adaptive processes in chronic bacterial infections are well described, but much less is known about the processes at play during acute infections. Here, by sequencing seven randomly selected isolates per patient, we analyzed Escherichia coli populations from three acute extraintestinal infections in adults (meningitis, pyelonephritis, and peritonitis), in which a high-mutation-rate isolate or mutator isolate was found. The isolates of single patients displayed between a few dozen and more than 200 independent mutations, with up to half being specific to the mutator isolate. Multiple signs of positive selection were evidenced: a high ratio of nonsynonymous to synonymous mutations (Ka /Ks ratio) and strong mutational convergence within and between patients, some of them at loci well known for their adaptive potential, such as rpoS, rbsR, fimH, and fliC For all patients, the mutator isolate was likely due to a large deletion of a methyl-directed mismatch repair gene, and in two instances, the deletion extended to genes involved in some genetic convergence, suggesting potential coselection. Intrinsic extraintestinal virulence assessed in a mouse model of sepsis showed variable patterns of virulence ranging from non-mouse killer to mouse killer for the isolates from single patients. However, genomic signature and gene inactivation experiments did not establish a link between a single gene and the capacity to kill mice, highlighting the complex and multifactorial nature of the virulence. Altogether, these data indicate that E. coli isolates are adapting under strong selective pressure when colonizing an extraintestinal site.IMPORTANCE Little is known about the dynamics of adaptation in acute bacterial infections. By sequencing multiple isolates from monoclonal extraintestinal Escherichia coli infections in several patients, we were able to uncover traces of selection taking place at short time scales compared to chronic infection. High genomic diversity was observed in the patient isolates, with an excess of nonsynonymous mutations, and the comparison within and between different infections showed patterns of convergence at the gene level, both constituting strong signs of adaptation. The genes targeted were coding mostly for proteins involved in global regulation, metabolism, and adhesion/motility. Moreover, virulence assessed in a mouse model of sepsis was variable among the isolates of single patients, but this difference was left unexplained at the molecular level. This work gives us clues about the E. coli lifestyle transition between commensalism and pathogenicity.

Evolution of a Dominant Natural Isolate of Escherichia coli in the Human Gut over the Course of a Year Suggests a Neutral Evolution with Reduced Effective Population Size.

In vitro and in vivo evolution experiments on Escherichia coli revealed several principles of bacterial adaptation. However, few data are available in the literature describing the behavior of E. coli in its natural environment. We attempted here to study the evolution in the human gut of a commensal dominant E. coli clone, ED1a belonging to the B2 phylogroup, through a longitudinal genomic study. We sequenced 24 isolates sampled at three different time points within a healthy individual over almost a year. We computed a mutation rate of 6.90 × 10-7 mutations per base per year of the chromosome for E. coli ED1a in healthy human gut. We observed very limited genomic diversity and could not detect any evidence of selection, in contrast to what is observed in experimental evolution over a similar length of time. We therefore suggest that ED1a, being well adapted to the healthy human gut, evolves mostly neutrally with a low effective population size (Ne of ≈500 to 1,700).IMPORTANCE In this study, we follow the genomic fate of a dominant clone of Escherichia coli in the human gut of a healthy individual over about a year. We could compute a low annual mutation rate that supports low diversity, and we could not retrieve any clear signature of selection. These observations support a neutral evolution of E. coli in the human gut, compatible with a very limited effective population size that deviates drastically with the observations made previously in experimental evolution.

Using long-term experimental evolution to uncover the patterns and determinants of molecular evolution of an Escherichia coli natural isolate in the streptomycin-treated mouse gut.

Although microbial ecology of the gut is now a major focus of interest, little is known about the molecular determinants of microbial adaptation in the gut. Experimental evolution coupled with whole-genome sequencing can provide insights of the adaptive process. In vitro experiments have revealed some conserved patterns: intermediate convergence, and epistatic interactions between beneficial mutations and mutations in global regulators. To test the relevance of these patterns and to identify the selective pressures acting in vivo, we have performed a long-term adaptation of an E. coli natural isolate, the streptomycin-resistant strain 536, in the digestive tract of streptomycin-treated mice. After a year of evolution, a clone from 15 replicates was sequenced. Consistently with in vitro observations, the identified mutations revealed a strong pattern of convergence at the mutation, gene, operon and functional levels. Yet, the rate of molecular evolution was lower than in in vitro, and no mutations in global regulators were recovered. More specific targets were observed: the dgo operon, involved in the galactonate pathway that improved growth on D-galactonate, and rluD and gidB, implicated in the maturation of the ribosomes, which mutations improved growth only in the presence of streptomycin. As in vitro, the nonrandom associations of mutations within the same pathways suggested a role of epistasis in shaping the adaptive landscape. Overall, we show that 'evolve and sequence' approach coupled with an analysis of convergence, when applied to a natural isolate, can be used to study adaptation in vivo and uncover the specific selective pressures of that environment.

Links between Transcription, Environmental Adaptation and Gene Variability in Escherichia coli: Correlations between Gene Expression and Gene Variability Reflect Growth Efficiencies.

Gene expression is known to be the principle factor explaining how fast genes evolve. Highly transcribed genes evolve slowly because any negative impact caused by a particular mutation is magnified by protein abundance. However, gene expression is a phenotype that depends both on the environment and on the strains or species. We studied this phenotypic plasticity by analyzing the transcriptome profiles of four Escherichia coli strains grown in three different culture media, and explored how expression variability was linked to gene allelic diversity. Genes whose expression changed according to the media and not to the strains were less polymorphic than other genes. Genes for which transcription depended predominantly on the strain were more polymorphic than other genes and were involved in sensing and responding to environmental changes, with an overrepresentation of two-component system genes. Surprisingly, we found that the correlation between transcription and gene diversity was highly variable among growth conditions and could be used to quantify growth efficiency of a strain in a medium. Genetic variability was found to increase with gene expression in poor growth conditions. As such conditions are also characterized by down-regulation of all DNA repair systems, including transcription-coupled repair, we suggest that gene expression under stressful conditions may be mutagenic and thus leads to a variability in mutation rate among genes in the genome which contributes to the pattern of protein evolution.

Emergence of Antimicrobial-Resistant Escherichia coli of Animal Origin Spreading in Humans.

In the context of the great concern about the impact of human activities on the environment, we studied 403 commensal Escherichia coli/Escherichia clade strains isolated from several animal and human populations that have variable contacts to one another. Multilocus sequence typing (MLST) showed a decrease of diversity 1) in strains isolated from animals that had an increasing contact with humans and 2) in all strains that had increased antimicrobial resistance. A specific B1 phylogroup clonal complex (CC87, Institut Pasteur schema nomenclature) of animal origin was identified and characterized as being responsible for the increased antimicrobial resistance prevalence observed in strains from the environments with a high human-mediated antimicrobial pressure. CC87 strains have a high capacity of acquiring and disseminating resistance genes with specific metabolic and genetic determinants as demonstrated by high-throughput sequencing and phenotyping. They are good mouse gut colonizers but are not virulent. Our data confirm the predominant role of human activities in the emergence of antimicrobial resistance in the environmental bacterial strains and unveil a particular E. coli clonal complex of animal origin capable of spreading antimicrobial resistance to other members of microbial communities.

The conserved nhaAR operon is drastically divergent between B2 and non-B2 Escherichia coli and is involved in extra-intestinal virulence.

The Escherichia coli species is divided in phylogenetic groups that differ in their virulence and commensal distribution. Strains belonging to the B2 group are involved in extra-intestinal pathologies but also appear to be more prevalent as commensals among human occidental populations. To investigate the genetic specificities of B2 sub-group, we used 128 sequenced genomes and identified genes of the core genome that showed marked difference between B2 and non-B2 genomes. We focused on the gene and its surrounding region with the strongest divergence between B2 and non-B2, the antiporter gene nhaA. This gene is part of the nhaAR operon, which is in the core genome but flanked by mobile regions, and is involved in growth at high pH and high sodium concentrations. Consistently, we found that a panel of non-B2 strains grew faster than B2 at high pH and high sodium concentrations. However, we could not identify differences in expression of the nhaAR operon using fluorescence reporter plasmids. Furthermore, the operon deletion had no differential impact between B2 and non-B2 strains, and did not result in a fitness modification in a murine model of gut colonization. Nevertheless, sequence analysis and experiments in a murine model of septicemia revealed that recombination in nhaA among B2 strains was observed in strains with low virulence. Finally, nhaA and nhaAR operon deletions drastically decreased virulence in one B2 strain. This effect of nhaAR deletion appeared to be stronger than deletion of all pathogenicity islands. Thus, a population genetic approach allowed us to identify an operon in the core genome without strong effect in commensalism but with an important role in extra-intestinal virulence, a landmark of the B2 strains.