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.

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 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.

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.

The decoupling between genetic structure and metabolic phenotypes in Escherichia coli leads to continuous phenotypic diversity

To assess the extent of intra-species diversity and the links between phylogeny, lifestyle (habitat and pathogenicity) and phenotype, we assayed the growth yield on 95 carbon sources of 168 Escherichia strains. We also correlated the growth capacities of 14 E. coli strains with the presence/absence of enzyme-coding genes. Globally, we found that the genetic distance, based on multilocus sequence typing data, was a weak indicator of the metabolic phenotypic distance. Besides, lifestyle and phylogroup had almost no impact on the growth yield of non-Shigella E. coli strains. In these strains, the presence/absence of the metabolic pathways, which was linked to the phylogeny, explained most of the growth capacities. However, few discrepancies blurred the link between metabolic phenotypic distance and metabolic pathway distance. This study shows that a prokaryotic species structured into well-defined genetic and lifestyle groups can yet exhibit continuous phenotypic diversity, possibly caused by gene regulatory effects.

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.