Predicting the effect of mutations to investigate recent events of selection across 60,472 Escherichia coli strains.

Microbial genomics studies focusing on the dynamics of selection have often used a small number of distant genomes. As a result, they could only analyze mutations that had become fixed during the divergence between species. However, thousands of genomes of some species are now available in public databases, thanks to high-throughput sequencing. These data provide a more complete picture of the polymorphisms segregating within a species, offering a unique insight into the processes that shape the recent evolution of a species. In this study, we present GLASS (Gene-Level Amino-acid Score Shift), a selection test that is based on the predicted effects of amino acid changes. By comparing the distribution of effects of mutations observed in a gene to the expectation in the absence of selection, GLASS can quantify the intensity of selection. We applied GLASS to a dataset of 60,472 Escherichia coli strains and used this to reexamine the longstanding debate about the role of essentiality versus expression level in the rate of protein evolution. We found that selection has contrasting short-term and long-term dynamics, with essential genes being subject to strong purifying selection in the short term, while expression level determines the rate of gene evolution in the long term. GLASS also found an overrepresentation of inactivating mutations in specific transcription factors, such as efflux pump repressors, which is consistent with selection for antibiotic resistance. These gene-inactivating polymorphisms do not reach fixation, suggesting another contrast between short-term fitness gains and long-term counterselection.

rRNA operon multiplicity as a bacterial genome stability insurance policy.

Quick growth restart after upon encountering favourable environmental conditions is a major fitness contributor in natural environment. It is widely assumed that the time required to restart growth after nutritional upshift is determined by how long it takes for cells to synthesize enough ribosomes to produce the proteins required to reinitiate growth. Here we show that a reduction in the capacity to synthesize ribosomes by reducing number of ribosomal RNA (rRNA) operons (rrn) causes a longer transition from stationary phase to growth of Escherichia coli primarily due to high mortality rates. Cell death results from DNA replication blockage and massive DNA breakage at the sites of the remaining rrn operons that become overloaded with RNA polymerases (RNAPs). Mortality rates and growth restart duration can be reduced by preventing R-loop formation and improving DNA repair capacity. The same molecular mechanisms determine the duration of the recovery phase after ribosome-damaging stresses, such as antibiotics, exposure to bile salts or high temperature. Our study therefore suggests that a major function of rrn operon multiplicity is to ensure that individual rrn operons are not saturated by RNAPs, which can result in catastrophic chromosome replication failure and cell death during adaptation to environmental fluctuations.

The ability to modulate translation capacity, which resides greatly on a number of ribosomes, provides robustness in fluctuating environments. Because translation is energetically the most expensive process in cells, cells must constantly adapt the rate of ribosome production to resource availability. This is primarily achieved by regulating ribosomal RNA (rRNA) synthesis, to which ribosomal proteins synthesis is adjusted. The multiplicity of rRNA encoding operons per bacterial genome exceeds requirements for the maximal growth rates in non-stress conditions. In this study, the authors provide evidence that a major function of rRNA operon multiplicity is to ensure that individual operons are not saturated by RNA polymerases during adaptation to environmental fluctuations, which can result in catastrophic chromosome replication failure and cell death.

Klebsiella pneumoniae Carbapenemase Variants Resistant to Ceftazidime-Avibactam: an Evolutionary Overview.

First variants of the Klebsiella pneumoniae carbapenemase (KPC), KPC-2 and KPC-3, have encountered a worldwide success, particularly in K. pneumoniae isolates. These beta-lactamases conferred resistance to most beta-lactams including carbapenems but remained susceptible to new beta-lactam/beta-lactamase inhibitors, such as ceftazidime-avibactam. After the marketing of ceftazidime-avibactam, numerous variants of KPC resistant to this association have been described among isolates recovered from clinical samples or derived from experimental studies. In KPC variants resistant to ceftazidime-avibactam, point mutations, insertions and/or deletions have been described in various hot spots. Deciphering the impact of these mutations is crucial, not only from a therapeutic point of view, but also to follow the evolution in time and space of KPC variants resistant to ceftazidime-avibactam. In this review, we describe the mutational landscape of the KPC beta-lactamase toward ceftazidime-avibactam resistance based on a multidisciplinary approach including epidemiology, microbiology, enzymology, and thermodynamics. We show that resistance is associated with three hot spots, with a high representation of insertions and deletions compared with other class A beta-lactamases. Moreover, extension of resistance to ceftazidime-avibactam is associated with a trade-off in the resistance to other beta-lactams and a decrease in enzyme stability. Nevertheless, the high natural stability of KPC could underlay the propensity of this enzyme to acquire in vivo mutations conferring resistance to ceftazidime-avibactam (CAZavi), particularly via insertions and deletions.

The Population Genomics of Increased Virulence and Antibiotic Resistance in Human Commensal Escherichia coli over 30 Years in France.

Escherichia coli is a commensal species of the lower intestine but is also a major pathogen causing intestinal and extraintestinal infections that is increasingly prevalent and resistant to antibiotics. Most studies on genomic evolution of E. coli used isolates from infections. Here, instead, we whole-genome sequenced a collection of 403 commensal E. coli isolates from fecal samples of healthy adult volunteers in France (1980 to 2010). These isolates were distributed mainly in phylogroups A and B2 (30% each) and belonged to 152 sequence types (STs), the five most frequent being ST10 (phylogroup A; 16.3%), ST73 and ST95 (phylogroup B2; 6.3 and 5.0%, respectively), ST69 (phylogroup D; 4.2%), and ST59 (phylogroup F; 3.9%), and 224 O:H serotypes. ST and serotype diversity increased over time. The O1, O2, O6, and O25 groups used in bioconjugate O-antigen vaccine against extraintestinal infections were found in 23% of the strains of our collection. The increase in frequency of virulence-associated genes and antibiotic resistance was driven by two evolutionary mechanisms. Evolution of virulence gene frequency was driven by both clonal expansion of STs with more virulence genes ("ST-driven") and increases in gene frequency within STs independent of changes in ST frequencies ("gene-driven"). In contrast, the evolution of resistance was dominated by increases in frequency within STs ("gene-driven"). This study provides a unique picture of the phylogenomic evolution of E. coli in its human commensal habitat over 30 years and will have implications for the development of preventive strategies. IMPORTANCE Escherichia coli is an opportunistic pathogen with the greatest burden of antibiotic resistance, one of the main causes of bacterial infections and an increasing concern in an aging population. Deciphering the evolutionary dynamics of virulence and antibiotic resistance in commensal E. coli is important to understand adaptation and anticipate future changes. The gut of vertebrates is the primary habitat of E. coli and probably where selection for virulence and resistance takes place. Unfortunately, most whole-genome-sequenced strains are isolated from pathogenic conditions. Here, we whole-genome sequenced 403 E. coli commensals isolated from healthy French subjects over a 30-year period. Virulence genes increased in frequency by both clonal expansion of clones carrying them and increases in frequency within clones, whereas resistance genes increased by within-clone increased frequency. Prospective studies of E. coli commensals should be performed worldwide to have a broader picture of evolution and adaptation of this species.

Success of prophylactic antiviral therapy for SARS-CoV-2: Predicted critical efficacies and impact of different drug-specific mechanisms of action.

Repurposed drugs that are safe and immediately available constitute a first line of defense against new viral infections. Despite limited antiviral activity against SARS-CoV-2, several drugs are being tested as medication or as prophylaxis to prevent infection. Using a stochastic model of early phase infection, we evaluate the success of prophylactic treatment with different drug types to prevent viral infection. We find that there exists a critical efficacy that a treatment must reach in order to block viral establishment. Treatment by a combination of drugs reduces the critical efficacy, most effectively by the combination of a drug blocking viral entry into cells and a drug increasing viral clearance. Below the critical efficacy, the risk of infection can nonetheless be reduced. Drugs blocking viral entry into cells or enhancing viral clearance reduce the risk of infection more than drugs that reduce viral production in infected cells. The larger the initial inoculum of infectious virus, the less likely is prevention of an infection. In our model, we find that as long as the viral inoculum is smaller than 10 infectious virus particles, viral infection can be prevented almost certainly with drugs of 90% efficacy (or more). Even when a viral infection cannot be prevented, antivirals delay the time to detectable viral loads. The largest delay of viral infection is achieved by drugs reducing viral production in infected cells. A delay of virus infection flattens the within-host viral dynamic curve, possibly reducing transmission and symptom severity. Thus, antiviral prophylaxis, even with reduced efficacy, could be efficiently used to prevent or alleviate infection in people at high risk.

The interplay between anticancer challenges and the microbial communities from the gut.

Cancer being an increasing burden on human health, the use of anticancer drugs has risen over the last decades. The physiological effects of these drugs are not only perceived by the host's cells but also by the microbial cells it harbors as commensals, notably the gut microbiota. Since the early '50 s, the cytotoxicity of anticancer chemotherapy was evaluated on bacteria revealing some antimicrobial activities that result in an established perturbation of the gut microbiota. This perturbation can affect the host's health through dysbiosis, which can lead to multiple complications, but has also been shown to have a direct effect on the treatment efficiency.We, therefore, conducted a review of literature focusing on this triangular relationship involving the microbial communities from the gut, the host's disease, and the anticancer treatment. We focused specifically on the antimicrobial effects of anticancer chemotherapy, their impact on mutagenesis in bacteria, and the perspectives of using bacteria-based tools to help in the diagnostic and treatment of cancer.

Tempo and mode of genome evolution in a 50,000-generation experiment.

Adaptation by natural selection depends on the rates, effects and interactions of many mutations, making it difficult to determine what proportion of mutations in an evolving lineage are beneficial. Here we analysed 264 complete genomes from 12 Escherichia coli populations to characterize their dynamics over 50,000 generations. The populations that retained the ancestral mutation rate support a model in which most fixed mutations are beneficial, the fraction of beneficial mutations declines as fitness rises, and neutral mutations accumulate at a constant rate. We also compared these populations to mutation-accumulation lines evolved under a bottlenecking regime that minimizes selection. Nonsynonymous mutations, intergenic mutations, insertions and deletions are overrepresented in the long-term populations, further supporting the inference that most mutations that reached high frequency were favoured by selection. These results illuminate the shifting balance of forces that govern genome evolution in populations adapting to a new environment.

CRISPR/Cas9 recombineering-mediated deep mutational scanning of essential genes in Escherichia coli.

Deep mutational scanning can provide significant insights into the function of essential genes in bacteria. Here, we developed a high-throughput method for mutating essential genes of Escherichia coli in their native genetic context. We used Cas9-mediated recombineering to introduce a library of mutations, created by error-prone PCR, within a gene fragment on the genome using a single gRNA pre-validated for high efficiency. Tracking mutation frequency through deep sequencing revealed biases in the position and the number of the introduced mutations. We overcame these biases by increasing the homology arm length and blocking mismatch repair to achieve a mutation efficiency of 85% for non-essential genes and 55% for essential genes. These experiments also improved our understanding of poorly characterized recombineering process using dsDNA donors with single nucleotide changes. Finally, we applied our technology to target rpoB, the beta subunit of RNA polymerase, to study resistance against rifampicin. In a single experiment, we validate multiple biochemical and clinical observations made in the previous decades and provide insights into resistance compensation with the study of double mutants.

Elucidating the molecular architecture of adaptation via evolve and resequence experiments.

Evolve and resequence (E&R) experiments use experimental evolution to adapt populations to a novel environment, then next-generation sequencing to analyse genetic changes. They enable molecular evolution to be monitored in real time on a genome-wide scale. Here, we review the field of E&R experiments across diverse systems, ranging from simple non-living RNA to bacteria, yeast and the complex multicellular organism Drosophila melanogaster. We explore how different evolutionary outcomes in these systems are largely consistent with common population genetics principles. Differences in outcomes across systems are largely explained by different starting population sizes, levels of pre-existing genetic variation, recombination rates and adaptive landscapes. We highlight emerging themes and inconsistencies that future experiments must address.

KPC Beta-Lactamases Are Permissive to Insertions and Deletions Conferring Substrate Spectrum Modifications and Resistance to Ceftazidime-Avibactam.

To explore the mutational possibilities of insertions and deletions (indels) in the Klebsiella pneumoniae carbapenemase (KPC) beta-lactamase, we selected for ceftazidime-avibactam-resistant mutants. Of 96 screened mutants, we obtained 19 indels (2 to 15 amino acids), all located in the loops surrounding the active site. Three antibiotic susceptibility phenotypes emerged: an extended-spectrum-beta-lactamase-like phenotype, an activity restricted to ceftazidime, and a carbapenem-susceptible KPC-like phenotype. Tolerance for indels reflects the evolvability of KPC beta-lactamase, which could challenge the therapeutic management of patients.

Genomic characterization of 16S rRNA methyltransferase-producing Escherichia coli isolates from the Parisian area, France.

BACKGROUND: The resistance to all aminoglycosides (AGs) conferred by 16S rRNA methyltransferase enzymes (16S-RMTases) is a major public health concern. OBJECTIVES: To characterize the resistance genotype, its genetic environment and plasmid support, and the phylogenetic relatedness of 16S-RMTase-producing Escherichia coli from France. METHODS: We screened 137 E. coli isolates resistant to all clinically relevant AGs from nine Parisian hospitals for 16S-RMTases. WGS was performed on clinical isolates with high-level AG resistance (MIC ≥256 mg/L) and their transformants. RESULTS: Thirty of the 137 AG-resistant E. coli produced 16S-RMTases: 11 ArmA, 18 RmtB and 1 RmtC. The 16S-RMTase producers were also resistant to third-generation cephalosporins (90% due to a blaCTX-M gene), co-trimoxazole, fluoroquinolones and carbapenems (blaNDM and blaVIM genes) in 97%, 83%, 70% and 10% of cases, respectively. Phylogenomic diversity was high in ArmA producers, with 10 different STs, but a similar genetic environment, with the Tn1548 transposon carried by a plasmid closely related to pCTX-M-3 in 6/11 isolates. Conversely, RmtB producers belonged to 12 STs, the most frequent being ST405 and ST complex (STc) 10 (four and four isolates, respectively). The rmtB gene was carried by IncF plasmids in 10 isolates and was found in different genetic environments. The rmtC gene was carried by the pNDM-US plasmid. CONCLUSIONS: ArmA and RmtB are the predominant 16S-RMTases in France, but their spread follows two different patterns: (i) dissemination of a conserved genetic support carrying armA in E. coli with high levels of genomic diversity; and (ii) various genetic environments surrounding rmtB in clonally related E. coli.

Mutator genomes decay, despite sustained fitness gains, in a long-term experiment with bacteria.

Understanding the extreme variation among bacterial genomes remains an unsolved challenge in evolutionary biology, despite long-standing debate about the relative importance of natural selection, mutation, and random drift. A potentially important confounding factor is the variation in mutation rates between lineages and over evolutionary history, which has been documented in several species. Mutation accumulation experiments have shown that hypermutability can erode genomes over short timescales. These results, however, were obtained under conditions of extremely weak selection, casting doubt on their general relevance. Here, we circumvent this limitation by analyzing genomes from mutator populations that arose during a long-term experiment with Escherichia coli, in which populations have been adaptively evolving for >50,000 generations. We develop an analytical framework to quantify the relative contributions of mutation and selection in shaping genomic characteristics, and we validate it using genomes evolved under regimes of high mutation rates with weak selection (mutation accumulation experiments) and low mutation rates with strong selection (natural isolates). Our results show that, despite sustained adaptive evolution in the long-term experiment, the signature of selection is much weaker than that of mutational biases in mutator genomes. This finding suggests that relatively brief periods of hypermutability can play an outsized role in shaping extant bacterial genomes. Overall, these results highlight the importance of genomic draft, in which strong linkage limits the ability of selection to purge deleterious mutations. These insights are also relevant to other biological systems evolving under strong linkage and high mutation rates, including viruses and cancer cells.

Advantage of the F2:A1:B- IncF Pandemic Plasmid over IncC Plasmids in In Vitro Acquisition and Evolution of blaCTX-M Gene-Bearing Plasmids in Escherichia coli.

Despite a fitness cost imposed on bacterial hosts, large conjugative plasmids play a key role in the diffusion of resistance determinants, such as CTX-M extended-spectrum ß-lactamases. Among the large conjugative plasmids, IncF plasmids are the most predominant group, and an F2:A1:B- IncF-type plasmid encoding a CTX-M-15 variant was recently described as being strongly associated with the emerging worldwide Escherichia coli sequence type 131 (ST131)-O25b:H4 H30Rx/C2 sublineage. In this context, we investigated the fitness cost of narrow-range F-type plasmids, including the F2:A1:B- IncF-type CTX-M-15 plasmid, and of broad-range C-type plasmids in the K-12-like J53-2 E. coli strain. Although all plasmids imposed a significant fitness cost to the bacterial host immediately after conjugation, we show, using an experimental-evolution approach, that a negative impact on the fitness of the host strain was maintained throughout 1,120 generations with the IncC-IncR plasmid, regardless of the presence or absence of cefotaxime, in contrast to the F2:A1:B- IncF plasmid, whose cost was alleviated. Many chromosomal and plasmid rearrangements were detected after conjugation in transconjugants carrying the IncC plasmids but not in transconjugants carrying the F2:A1:B- IncF plasmid, except for insertion sequence (IS) mobilization from the fliM gene leading to the restoration of motility of the recipient strains. Only a few mutations occurred on the chromosome of each transconjugant throughout the experimental-evolution assay. Our findings indicate that the F2:A1:B- IncF CTX-M-15 plasmid is well adapted to the E. coli strain studied, contrary to the IncC-IncR CTX-M-15 plasmid, and that such plasmid-host adaptation could participate in the evolutionary success of the CTX-M-15-producing pandemic E. coli ST131-O25b:H4 lineage.

Ceftriaxone and Cefotaxime Have Similar Effects on the Intestinal Microbiota in Human Volunteers Treated by Standard-Dose Regimens.

Ceftriaxone has a higher biliary elimination than cefotaxime (40% versus 10%), which may result in a more pronounced impact on the intestinal microbiota. We performed a monocenter, randomized open-label clinical trial in 22 healthy volunteers treated by intravenous ceftriaxone (1 g/24 h) or cefotaxime (1 g/8 h) for 3 days. We collected fecal samples for phenotypic analyses, 16S rRNA gene profiling, and measurement of the antibiotic concentration and compared the groups for the evolution of microbial counts and indices of bacterial diversity over time. Plasma samples were drawn at day 3 for pharmacokinetic analysis. The emergence of 3rd-generation-cephalosporin-resistant Gram-negative enteric bacilli (Enterobacterales), Enterococcus spp., or noncommensal microorganisms was not significantly different between the groups. Both antibiotics reduced the counts of total Gram-negative enteric bacilli and decreased the bacterial diversity, but the differences between the groups were not significant. All but one volunteer from each group exhibited undetectable levels of antibiotic in feces. Plasma pharmacokinetic endpoints were not correlated to alteration of the bacterial diversity of the gut. Both antibiotics markedly impacted the intestinal microbiota, but no significant differences were detected when standard clinical doses were administered for 3 days. This might be related to the similar daily amounts of antibiotics excreted through the bile using a clinical regimen. (This study has been registered at ClinicalTrials.gov under identifier NCT02659033.).

Long-term evolution of the natural isolate of Escherichia coli 536 in the mouse gut colonized after maternal transmission reveals convergence in the constitutive expression of the lactose operon.

In vitro experimental evolution has taught us many lessons on the molecular bases of adaptation. To move towards more natural settings, evolution in the mice gut has been successfully performed. Yet, these experiments suffered from the use of laboratory strains as well as the use of axenic or streptomycin-treated mice to maintain the inoculated strains. To circumvent these limitations, we conducted a one-year experimental evolution in vivo using a natural isolate of E. coli, strain 536, in conditions mimicking as much as possible natural environment with mother-to-offspring microbiota transmission. Mice were then distributed in 24 independent cages and separated into two different diets: a regular one (chow diet, CD) and high-fat and high-sugar one (Western Diet, WD). Genome sequences revealed an early and rapid selection during the breastfeeding period that selected the constitutive expression of the well-characterized lactose operon. E. coli was lost significantly more in CD than WD; however, we could not detect any genomic signature of selection, nor any diet specificities during the later part of the experiments. The apparently neutral evolution presumably due to low population size maintained nevertheless at high frequency the early selected mutations affecting lactose regulation. The rapid loss of lactose operon regulation challenges the idea that plastic gene expression is both optimal and stable in the wild.

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.

First-Step Mutations during Adaptation Restore the Expression of Hundreds of Genes.

The temporal change of phenotypes during the adaptive process remains largely unexplored, as do the genetic changes that affect these phenotypic changes. Here we focused on three mutations that rose to high frequency in the early stages of adaptation within 12 Escherichia coli populations subjected to thermal stress (42 °C). All the mutations were in the rpoB gene, which encodes the RNA polymerase beta subunit. For each mutation, we measured the growth curves and gene expression (mRNAseq) of clones at 42 °C. We also compared growth and gene expression with their ancestor under unstressed (37 °C) and stressed conditions (42 °C). Each of the three mutations changed the expression of hundreds of genes and conferred large fitness advantages, apparently through the restoration of global gene expression from the stressed toward the prestressed state. These three mutations had a similar effect on gene expression as another single mutation in a distinct domain of the rpoB protein. Finally, we compared the phenotypic characteristics of one mutant, I572L, with two high-temperature adapted clones that have this mutation plus additional background mutations. The background mutations increased fitness, but they did not substantially change gene expression. We conclude that early mutations in a global transcriptional regulator cause extensive changes in gene expression, many of which are likely under positive selection for their effect in restoring the prestress physiology.

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.

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.