MetaPhlAn 4 profiling of unknown species-level genome bins improves the characterization of diet-associated microbiome changes in mice.

Mouse models are key tools for investigating host-microbiome interactions. However, shotgun metagenomics can only profile a limited fraction of the mouse gut microbiome. Here, we employ a metagenomic profiling method, MetaPhlAn 4, which exploits a large catalog of metagenome-assembled genomes (including 22,718 metagenome-assembled genomes from mice) to improve the profiling of the mouse gut microbiome. We combine 622 samples from eight public datasets and an additional cohort of 97 mouse microbiomes, and we assess the potential of MetaPhlAn 4 to better identify diet-related changes in the host microbiome using a meta-analysis approach. We find multiple, strong, and reproducible diet-related microbial biomarkers, largely increasing those identifiable by other available methods relying only on reference information. The strongest drivers of the diet-induced changes are uncharacterized and previously undetected taxa, confirming the importance of adopting metagenomic methods integrating metagenomic assemblies for comprehensive profiling.

The Arginine Deiminase Operon Is Responsible for a Fitness Trade-Off in Extended-Spectrum-ß-Lactamase-Producing Strains of Escherichia coli.

We previously identified an operon involved in an arginine deiminase (ADI) pathway (arc operon) on a CTX-M-producing plasmid from an O102-ST405 strain of Escherichia coli As the ADI pathway was shown to be involved in the virulence of various Gram-positive bacteria, we tested whether the ADI pathway could be involved in the epidemiological success of extended-spectrum-ß-lactamase (ESBL)-producing E. coli strains. We studied two collections of human E. coli isolated in France (n = 493) and England (n = 1,509) and show that the prevalence of the arc operon (i) is higher in ESBL-producing strains (12.1%) than in nonproducers (2.5%), (ii) is higher in CTX-M-producing strains (16%) than in other ESBL producers (3.5%), and (iii) increased over time in ESBL-producing strains from 0% before 2000 to 43.3% in 2011 to 2012. The arc operon, found in strains from various phylogenetic backgrounds, is carried by IncF plasmids (85%) or chromosomes (15%) in regions framed by numerous insertion sequences, indicating multiple arrivals. Competition experiments showed that the arc operon enhances fitness of the strain in vitro in lysogeny broth with arginine. In vivo competition experiments showed that the arc operon is advantageous for the strain in a mouse model of urinary tract infection (UTI), whereas it is a burden in a mouse model of intestinal colonization. In summary, we have identified a trait linked to CTX-M-producing strains that is responsible for a trade-off between two main E. coli lifestyles, UTI and gut commensalism. This trait alone cannot explain the wide spread of ESBLs in E. coli but merits epidemiological surveillance.

Enterobacteriaceae are essential for the modulation of colitis severity by fungi.

BACKGROUND: Host-microbe balance maintains intestinal homeostasis and strongly influences inflammatory conditions such as inflammatory bowel diseases (IBD). Here we focused on bacteria-fungi interactions and their implications on intestinal inflammation, a poorly understood area. METHODS: Dextran sodium sulfate (DSS)-induced colitis was assessed in mice treated with vancomycin (targeting gram-positive bacteria) or colistin (targeting Enterobacteriaceae) and supplemented with either Saccharomyces boulardii CNCM I-745 or Candida albicans. Inflammation severity as well as bacterial and fungal microbiota compositions was monitored. RESULTS: While S. boulardii improved DSS-induced colitis and C. albicans worsened it in untreated settings, antibiotic treatment strongly modified DSS susceptibility and effects of fungi on colitis. Vancomycin-treated mice were fully protected from colitis, while colistin-treated mice retained colitis phenotype but were not affected anymore by administration of fungi. Antibacterial treatments not only influenced bacterial populations but also had indirect effects on fungal microbiota. Correlations between bacterial and fungal relative abundance were dramatically decreased in colistin-treated mice compared to vancomycin-treated and control mice, suggesting that colistin-sensitive bacteria are involved in interactions with fungi. Restoration of the Enterobacteriaceae population by administrating colistin-resistant Escherichia coli reestablished both beneficial effects of S. boulardii and pathogenic effects of C. albicans on colitis severity. This effect was at least partly mediated by an improved gut colonization by fungi. CONCLUSIONS: Fungal colonization of the gut is affected by the Enterobacteriaceae population, indirectly modifying effects of mycobiome on the host. This finding provides new insights into the role of inter-kingdom functional interactions in intestinal physiopathology and potentially in IBD.

Extended-spectrum ß-lactamase-encoding genes are spreading on a wide range of Escherichia coli plasmids existing prior to the use of third-generation cephalosporins.

To understand the evolutionary dynamics of extended-spectrum ß-lactamase (ESBL)-encoding genes in Escherichia coli, we undertook a comparative genomic analysis of 116 whole plasmid sequences of human or animal origin isolated over a period spanning before and after the use of third-generation cephalosporins (3GCs) using a gene-sharing network approach. The plasmids included 82 conjugative, 22 mobilizable and 9 non-transferable plasmids and 3 P-like bacteriophages. ESBL-encoding genes were found on 64 conjugative, 6 mobilizable, 2 non-transferable plasmids and 2 P1-like bacteriophages, indicating that these last three types of mobile elements also play a role, albeit modest, in the diffusion of the ESBLs. The network analysis showed that the plasmids clustered according to their genome backbone type, but not by origin or period of isolation or by antibiotic-resistance type, including type of ESBL-encoding gene. There was no association between the type of plasmid and the phylogenetic history of the parental strains. Finer scale analysis of the more abundant clusters IncF and IncI1 showed that ESBL-encoding plasmids and plasmids isolated before the use of 3GCs had the same diversity and phylogenetic history, and that acquisition of ESBL-encoding genes had occurred during multiple independent events. Moreover, the blaCTX-M-15 gene, unlike other CTX-M genes, was inserted at a hot spot in a blaTEM-1-Tn2 transposon. These findings showed that ESBL-encoding genes have arrived on wide range of pre-existing plasmids and that the successful spread of blaCTX-M-15 seems to be favoured by the presence of well-adapted IncF plasmids that carry a Tn2-blaTEM-1 transposon.

Impaired Aryl Hydrocarbon Receptor Ligand Production by the Gut Microbiota Is a Key Factor in Metabolic Syndrome.

The extent to which microbiota alterations define or influence the outcome of metabolic diseases is still unclear, but the byproducts of microbiota metabolism are known to have an important role in mediating the host-microbiota interaction. Here, we identify that in both pre-clinical and clinical settings, metabolic syndrome is associated with the reduced capacity of the microbiota to metabolize tryptophan into derivatives that are able to activate the aryl hydrocarbon receptor. This alteration is not merely an effect of the disease as supplementation with AhR agonist or a Lactobacillus strain, with a high AhR ligand-production capacity, leads to improvement of both dietary- and genetic-induced metabolic impairments, particularly glucose dysmetabolism and liver steatosis, through improvement of intestinal barrier function and secretion of the incretin hormone GLP-1. These results highlight the role of gut microbiota-derived metabolites as a biomarker and as a basis for novel preventative or therapeutic interventions for metabolic disorders.

Diversity of the auxotrophic requirements in natural isolates of Escherichia coli.

Isolates of Escherichia coli, except Shigella, are generally prototrophic; they do not require any growth factors to grow in mineral medium. However, a nicotinic acid requirement is common among B2 phylogroup STc95 O18 E. coli clone strains. Nicotinic acid is a precursor of nicotinamide adenine dinucleotide (NAD), an essential molecule that plays central role in cellular metabolism. The defect in NAD synthesis of these strains is due to alterations in de novo biosynthesis pathway nadB gene. Here, by studying growth on minimal medium with glycolytic (glucose) or gluconeogenic (pyruvate or succinate) substrates as the carbon supply in a large panel of E. coli natural isolates representative of the species diversity, we identify an absolute nicotinic acid requirement in non-STc95 strains due in one case to a nadA inactivation. The growth on glucose medium of some extraintestinal pathogenic E. coli strains belonging to various non-O18 B2 phylogroup STc95 clones is restored either by aspartate or nicotinate, demonstrating that the nicotinic acid requirement can also be due to an intracellular aspartate depletion. The auxotrophic requirements depend on the carbon source available in the environment. Moreover, some strains prototrophic in glucose medium become auxotrophic in succinate medium, and conversely, some strainsauxotrophic in glucose medium become prototrophic in succinate medium. Finally, a partial depletion of intracellular aspartate can be observed in some prototrophic strains belonging to various phylogroups. The observed more or less significant depletion according to isolates may be due to differences in tricarboxylic acid cycle enzyme activities. These metabolic defects could be involved in the adaptation of E. coli to its various niches.

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.

Interactions between genotype and environment drive the metabolic phenotype within Escherichia coli isolates.

To gain insights into the adaptation of the Escherichia coli species to different environments, we monitored protein abundances using quantitative proteomics and measurements of enzymatic activities of central metabolism in a set of five representative strains grown in four contrasted culture media including human urine. Two hundred and thirty seven proteins representative of the genome-scale metabolic network were identified and classified into pathway categories. We found that nutrient resources shape the general orientation of metabolism through coordinated changes in the average abundances of proteins and in enzymatic activities that all belong to the same pathway category. For example, each culture medium induces a specific oxidative response whatever the strain. On the contrary, differences between strains concern isolated proteins and enzymes within pathway categories in single environments. Our study confirms the predominance of genotype by environment interactions at the proteomic and enzyme activity levels. The buffering of genetic variation when considering life-history traits suggests a multiplicity of evolutionary strategies. For instance, the uropathogenic isolate CFT073 shows a deregulation of iron demand and increased oxidative stress response.

The rpoS gene is predominantly inactivated during laboratory storage and undergoes source-sink evolution in Escherichia coli species.

The rpoS gene codes for an alternative RNA polymerase sigma factor, which acts as a general regulator of the stress response. Inactivating alleles of rpoS in collections of natural Escherichia coli isolates have been observed at very variable frequencies, from less than 1% to more than 70% of strains. rpoS is easily inactivated in nutrient-deprived environments such as stab storage, which makes it difficult to determine the true frequency of rpoS inactivation in nature. We studied the evolutionary history of rpoS and compared it to the phylogenetic history of bacteria in two collections of 82 human commensal and extraintestinal E. coli strains. These strains were representative of the phylogenetic diversity of the species and differed only by their storage conditions. In both collections, the phylogenetic histories of rpoS and of the strains were congruent, indicating that horizontal gene transfer had not occurred at the rpoS locus, and rpoS was under strong purifying selection, with a ratio of the nonsynonymous mutation rate (Ka) to the synonymous substitution rate (Ks) substantially smaller than 1. Stab storage was associated with a high frequency of inactivating alleles, whereas almost no amino acid sequence variation was observed in RpoS in the collection studied directly after isolation of the strains from the host. Furthermore, the accumulation of variations in rpoS was typical of source-sink dynamics. In conclusion, rpoS is rarely inactivated in natural E. coli isolates within their mammalian hosts, probably because such strains rapidly become evolutionary dead ends. Our data should encourage bacteriologists to freeze isolates immediately and to avoid the use of stab storage.

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.

Complexity of Escherichia coli bacteremia pathophysiology evidenced by comparison of isolates from blood and portal of entry within single patients.

The portal of entry of Escherichia coli bacteremia, a frequent and severe disease, is most commonly the urinary tract followed by the digestive tract. Recent reports have evidenced the presence of several distinct E. coli clones within a single patient suffering of extra-intestinal infection. To explore the relationships between the blood and portal of entry strains, we thoroughly studied 98 bacteremic patients from the French prospective COLIBAFI cohort. In these patients, we compared genotypically and phenotypically E. coli strains isolated from the blood and the suspected portal of entry [non-urinary pus (n=52) and urine (n=52)]. We found genetically distinct strains exhibiting distinct antibiotypes in the blood and pus samples (8 patients; 15%) and the blood and urine samples (2 patients; 3.8%) (p=0.09). These data highlight the complexity of pathophysiology of E. coli bacteremia and should be taken into consideration when strain antibiotic susceptibility is tested, especially in bacteremia of pus origin.

Capturing the mutational landscape of the beta-lactamase TEM-1.

Adaptation proceeds through the selection of mutations. The distribution of mutant fitness effect and the forces shaping this distribution are therefore keys to predict the evolutionary fate of organisms and their constituents such as enzymes. Here, by producing and sequencing a comprehensive collection of 10,000 mutants, we explore the mutational landscape of one enzyme involved in the spread of antibiotic resistance, the beta-lactamase TEM-1. We measured mutation impact on the enzyme activity through the estimation of amoxicillin minimum inhibitory concentration on a subset of 990 mutants carrying a unique missense mutation, representing 64% of possible amino acid changes in that protein reachable by point mutation. We established that mutation type, solvent accessibility of residues, and the predicted effect of mutations on protein stability primarily determined alone or in combination changes in minimum inhibitory concentration of mutants. Moreover, we were able to capture the drastic modification of the mutational landscape induced by a single stabilizing point mutation (M182T) by a simple model of protein stability. This work thereby provides an integrated framework to study mutation effects and a tool to understand/define better the epistatic interactions.

Commensal Escherichia coli strains in Guiana reveal a high genetic diversity with host-dependant population structure.

We undertook a large-scale epidemiological survey of commensal Escherichia coli in Trois-Sauts, an isolated village located in the south of French Guiana where human population exchanges are restricted and source of antibiotics controlled. Stools from 162 Wayampi Amerindians and rectal swabs from 33 human associated and 198 wild animals were collected in the close proximity of the village. The prevalence of E. coli was decreasing from humans (100%) to human associated (64%) and wild (45%) animals. A clear genetic structure between these three E. coli populations was observed with human strains belonging very rarely to B2 phylogroup (3.7%), exhibiting few virulence genes and bacteriocins but being antibiotic resistant whereas wild animal strains were characterized by 46.1% of B2 phylogroup belonging, with very unique and infrequent sequence types, numerous extraintestinal genes and bacteriocins but no antibiotic resistance; the human-associated animal strains being intermediate. Furthermore, an unexpected genetic diversity was observed among the strains, as the housekeeping gene nucleotide diversity per site of the Trois-Sauts's strains was higher than the one of reference strains representative of the known species diversity. The existence of such E. coli structured phylogenetic diversity within various hosts of a single localization has never been reported.

Variation in endogenous oxidative stress in Escherichia coli natural isolates during growth in urine.

BACKGROUND: Uropathogenic strains of Escherichia coli cause symptomatic infections whereas asymptomatic bacteriuria (ABU) strains are well adapted for growth in the human urinary tract, where they establish long-term bacteriuria. Human urine is a very complex growth medium that could be perceived by certain bacteria as a stressful environment. To investigate a possible imbalance between endogenous oxidative response and antioxidant mechanisms, lipid oxidative damage estimated as thiobarbituric acid reactive substances (TBARS) content was evaluated in twenty-one E. coli belonging to various pathovars and phylogenetic groups. Antioxidant defense mechanisms were also analysed. RESULTS: During exponential growth in urine, TBARS level differs between strains, without correlation with the ability to grow in urine which was similarly limited for commensal, ABU and uropathogenic strains. In addition, no correlation between TBARS level and the phylogroup or pathogenic group is apparent. The growth of ABU strain 83972 was associated with a high level of TBARS and more active antioxidant defenses that reduce the imbalance. CONCLUSIONS: Our results indicate that growth capacity in urine is not a property of ABU strains. However, E. coli isolates respond very differently to this stressful environment. In strain ABU 83972, on one hand, the increased level of endogenous reactive oxygen species may be responsible for adaptive mutations. On the other hand, a more active antioxidant defense system could increase the capacity to colonize the bladder.

Animal and human pathogenic Escherichia coli strains share common genetic backgrounds.

Escherichia coli is a versatile species encompassing both commensals of the digestive tracts of many vertebrates, including humans, and pathogenic strains causing various intra- and extraintestinal infections. Despite extensive gene flow between strains, the E. coli species has a globally clonal population structure, consisting of distinct phylogenetic groups. Little is known about the relationships between phylogenetic groups and host specificity. We therefore used multilocus sequence typing (MLST) to investigate phylogenetic relationships and evaluated the virulence gene content of 35 E. coli strains representative of the diverse diseases encountered in domestic animals. We compared these strains with a panel of 101 human pathogenic and 98 non-human and human commensal strains representative of the phylogenetic and pathovar diversity of this species. A global factorial analysis of correspondence indicated that extraintestinal infections were caused mostly by phylogenetic group B2 strains, whereas intraintestinal infections were caused mostly by phylogenetic group A/B1/E strains, with strains responsible from extraintestinal or intraintestinal infections having specific virulence factors. It was not possible to distinguish between strains of human and animal origin. A detailed phylogenetic analysis of the MLST data showed that numerous pathogenic animal and human strains are very closely related, and had a number of virulence genes in common. However, a set of specific adhesins was identified in animal non-B2 group strains of all pathotypes. In conclusion, human and animal pathogenic strains share common genetic backgrounds, but non-B2 strains of different origins seem to have different sets of adhesins that could be involved in host specificity.