Dataset Information

An Escherichia coli ST131 pangenome atlas reveals population structure and evolution across 4,071 isolates.

ABSTRACT: Escherichia coli ST131 is a major cause of infection with extensive antimicrobial resistance (AMR) facilitated by widespread beta-lactam antibiotic use. This drug pressure has driven extended-spectrum beta-lactamase (ESBL) gene acquisition and evolution in pathogens, so a clearer resolution of ST131's origin, adaptation and spread is essential. E. coli ST131's ESBL genes are typically embedded in mobile genetic elements (MGEs) that aid transfer to new plasmid or chromosomal locations, which are mobilised further by plasmid conjugation and recombination, resulting in a flexible ESBL, MGE and plasmid composition with a conserved core genome. We used population genomics to trace the evolution of AMR in ST131 more precisely by extracting all available high-quality Illumina HiSeq read libraries to investigate 4,071 globally-sourced genomes, the largest ST131 collection examined so far. We applied rigorous quality-control, genome de novo assembly and ESBL gene screening to resolve ST131's population structure across three genetically distinct Clades (A, B, C) and abundant subclades from the dominant Clade C. We reconstructed their evolutionary relationships across the core and accessory genomes using published reference genomes, long read assemblies and k-mer-based methods to contextualise pangenome diversity. The three main C subclades have co-circulated globally at relatively stable frequencies over time, suggesting attaining an equilibrium after their origin and initial rapid spread. This contrasted with their ESBL genes, which had stronger patterns across time, geography and subclade, and were located at distinct locations across the chromosomes and plasmids between isolates. Within the three C subclades, the core and accessory genome diversity levels were not correlated due to plasmid and MGE activity, unlike patterns between the three main clades, A, B and C. This population genomic study highlights the dynamic nature of the accessory genomes in ST131, suggesting that surveillance should anticipate genetically variable outbreaks with broader antibiotic resistance levels. Our findings emphasise the potential of evolutionary pangenomics to improve our understanding of AMR gene transfer, adaptation and transmission to discover accessory genome changes linked to novel subtypes.

SUBMITTER: Decano AG

PROVIDER: S-EPMC6874702 | biostudies-literature | 2019 Nov

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

An Escherichia coli ST131 pangenome atlas reveals population structure and evolution across 4,071 isolates.

Decano Arun Gonzales AG Downing Tim T

Scientific reports 20191122 1

Escherichia coli ST131 is a major cause of infection with extensive antimicrobial resistance (AMR) facilitated by widespread beta-lactam antibiotic use. This drug pressure has driven extended-spectrum beta-lactamase (ESBL) gene acquisition and evolution in pathogens, so a clearer resolution of ST131's origin, adaptation and spread is essential. E. coli ST131's ESBL genes are typically embedded in mobile genetic elements (MGEs) that aid transfer to new plasmid or chromosomal locations, which are ...[more]

PMID: 31758048

Similar Datasets

Project description:Background: Based on 32 Escherichia coli and Shigella genome sequences, we have developed an E. coli pan-genome microarray. Publicly available genomes were annotated in a consistent manor to define all currently known genes potentially present in the species. The chip design was evaluated by hybridization of DNA from two sequenced E. coli strains, K-12 MG1655 (a commensal) and O157:H7 EDL933 (an enterotoxigenic E. coli). A dual channel and single channel analysis approach was compared for the comparative genomic hybridization experiments. Moreover, the microarray was used to characterize four unsequenced probiotic E. coli strains, currently marketed for beneficial effects on the human gut flora. Results: Based on the genomes included in this study, we were able to group together 2,041 genes that were present in all 32 genomes. Furthermore, we predict that the size of the E. coli core genome will approach ~1,560 essential genes, considerably less than previous estimates. Although any individual E. coli genome contains between 4,000 and 5,000 genes, we identified more than twice as many (11,872) distinct gene groups in the total gene pool (“pan-genome”) examined for microarray design. Benchmarking of the design based on sequenced control strain samples demonstrated a high sensitivity and relatively low false positive rate. Moreover, the array was highly sufficient to investigate the gene content of apathogenic isolates, despite the strong bias towards pathogenic E. coli strains that have been sequenced so far. Our analysis of four probiotic E. coli strains demonstrate that they share a gene pool very similar to the E. coli K-12 strains but also show significant similarity with enteropathogenic strains. Nonetheless, virulence genes were largely absent. Strain-specific genes found in probiotic E. coli but absent in E. coli K12 were most frequently phage-related genes, transposases and other genes related to mobile DNA, and metabolic enzymes or factors that may offer colonization fitness, which together with their asymptomatic nature may explain their nature. Conclusion: This high-density microarray provides an excellent tool for characterizing either DNA content or gene expression from unknown E. coli strains. Keywords: Comparative genomic hybridizations

Project description:BackgroundEscherichia coli (E. coli) has been one of the most studied model organisms in the history of life sciences. Initially thought just to be commensal bacteria, E. coli has shown wide phenotypic diversity including pathogenic isolates with great relevance to public health. Though pangenome analysis has been attempted several times, there is no systematic functional characterization of the E. coli subgroups according to the gene profile.ResultsSystematically scanning for optimal parametrization, we have built the E. coli pangenome from 1324 complete genomes. The pangenome size is estimated to be ~25,000 gene families (GFs). Whereas the core genome diminishes as more genomes are added, the softcore genome (≥95% of strains) is stable with ~3000 GFs regardless of the total number of genomes. Apparently, the softcore genome (with a 92% or 95% generation threshold) can define the genome of a bacterial species listing the critically relevant, evolutionarily most conserved or important classes of GFs. Unsupervised clustering of common E. coli sequence types using the presence/absence GF matrix reveals distinct characteristics of E. coli phylogroups B1, B2, and E. We highlight the bi-lineage nature of B1, the variation of the secretion and of the iron acquisition systems in ST11 (E), and the incorporation of a highly conserved prophage into the genome of ST131 (B2). The tail structure of the prophage is evolutionarily related to R2-pyocin (a tailocin) from Pseudomonas aeruginosa PAO1. We hypothesize that this molecular machinery is highly likely to play an important role in protecting its own colonies; thus, contributing towards the rapid rise of pandemic E. coli ST131.ConclusionsThis study has explored the optimized pangenome development in E. coli. We provide complete GF lists and the pangenome matrix as supplementary data for further studies. We identified biological characteristics of different E. coli subtypes, specifically for phylogroups B1, B2, and E. We found an operon-like genome region coding for a tailocin specific for ST131 strains. The latter is a potential killer weapon providing pandemic E. coli ST131 with an advantage in inter-bacterial competition and, suggestively, explains their dominance as human pathogen among E. coli strains.

Project description:Background: Based on 32 Escherichia coli and Shigella genome sequences, we have developed an E. coli pan-genome microarray. Publicly available genomes were annotated in a consistent manor to define all currently known genes potentially present in the species. The chip design was evaluated by hybridization of DNA from two sequenced E. coli strains, K-12 MG1655 (a commensal) and O157:H7 EDL933 (an enterotoxigenic E. coli). A dual channel and single channel analysis approach was compared for the comparative genomic hybridization experiments. Moreover, the microarray was used to characterize four unsequenced probiotic E. coli strains, currently marketed for beneficial effects on the human gut flora. Results: Based on the genomes included in this study, we were able to group together 2,041 genes that were present in all 32 genomes. Furthermore, we predict that the size of the E. coli core genome will approach ~1,560 essential genes, considerably less than previous estimates. Although any individual E. coli genome contains between 4,000 and 5,000 genes, we identified more than twice as many (11,872) distinct gene groups in the total gene pool (âpan-genomeâ) examined for microarray design. Benchmarking of the design based on sequenced control strain samples demonstrated a high sensitivity and relatively low false positive rate. Moreover, the array was highly sufficient to investigate the gene content of apathogenic isolates, despite the strong bias towards pathogenic E. coli strains that have been sequenced so far. Our analysis of four probiotic E. coli strains demonstrate that they share a gene pool very similar to the E. coli K-12 strains but also show significant similarity with enteropathogenic strains. Nonetheless, virulence genes were largely absent. Strain-specific genes found in probiotic E. coli but absent in E. coli K12 were most frequently phage-related genes, transposases and other genes related to mobile DNA, and metabolic enzymes or factors that may offer colonization fitness, which together with their asymptomatic nature may explain their nature. Conclusion: This high-density microarray provides an excellent tool for characterizing either DNA content or gene expression from unknown E. coli strains. Factorial design: Each of four test samples (G 1/2, G3/10, G 4/9, G5) are co-hybridized with two control strain samples (K-12 MG1655 and O157:H7 EDL933). Additional replicate co-hybridizations are included of the two control strain samples (O157:H7 EDL933 vs. K-12 MG1655).

Project description:UnlabelledEscherichia coli ST131 is the most frequently isolated fluoroquinolone-resistant (FQR) E. coli clone worldwide and a major cause of urinary tract and bloodstream infections. Although originally identified through its association with the CTX-M-15 extended-spectrum β-lactamase resistance gene, global genomic epidemiology studies have failed to resolve the geographical and temporal origin of the ST131 ancestor. Here, we developed a framework for the reanalysis of publically available genomes from different countries and used this data set to reconstruct the evolutionary steps that led to the emergence of FQR ST131. Using Bayesian estimation, we show that point mutations in chromosomal genes that confer FQR coincide with the first clinical use of fluoroquinolone in 1986 and illustrate the impact of this pivotal event on the rapid population expansion of ST131 worldwide from an apparent origin in North America. Furthermore, we identify virulence factor acquisition events that predate the development of FQR, suggesting that the gain of virulence-associated genes followed by the tandem development of antibiotic resistance primed the successful global dissemination of ST131.ImportanceEscherichia coli sequence type 131 (ST131) is a recently emerged and globally disseminated multidrug-resistant clone frequently associated with human urinary tract and bloodstream infections. In this study, we have used two large publically available genomic data sets to define a number of critical steps in the evolution of this important pathogen. We show that resistance to fluoroquinolones, a class of broad-spectrum antibiotic used extensively in human medicine and veterinary practice, developed in ST131 soon after the introduction of these antibiotics in the United States, most likely in North America. We also mapped the acquisition of several fitness and virulence determinants by ST131 and demonstrate these events occurred prior to the development of fluoroquinolone resistance. Thus, ST131 has emerged by stealth, first acquiring genes associated with an increased capacity to cause human infection, and then gaining a resistance armory that has driven its massive population expansion across the globe.

Project description:Escherichia coli sequence type 131 (ST131) has emerged rapidly to become the most prevalent extraintestinal pathogenic E. coli clones in circulation today. Previous investigations appeared to exonerate retail meat as a source of human exposure to ST131; however, these studies focused mainly on extensively multidrug-resistant ST131 strains, which typically carry allele 30 of the fimH type 1 fimbrial adhesin gene (ST131-H30). To estimate the frequency of extraintestinal human infections arising from foodborne ST131 strains without bias toward particular sublineages or phenotypes, we conducted a 1-year prospective study of E. coli from meat products and clinical cultures in Flagstaff, Arizona. We characterized all isolates by multilocus sequence typing, fimH typing, and core genome phylogenetic analyses, and we screened isolates for avian-associated ColV plasmids as an indication of poultry adaptation. E. coli was isolated from 79.8% of the 2,452 meat samples and 72.4% of the 1,735 culture-positive clinical samples. Twenty-seven meat isolates were ST131 and belonged almost exclusively (n = 25) to the ST131-H22 lineage. All but 1 of the 25 H22 meat isolates were from poultry products, and all but 2 carried poultry-associated ColV plasmids. Of the 1,188 contemporaneous human clinical E. coli isolates, 24 were ST131-H22, one-quarter of which occurred in the same high-resolution phylogenetic clades as the ST131-H22 meat isolates and carried ColV plasmids. Molecular clock analysis of an international ST131-H22 genome collection suggested that ColV plasmids have been acquired at least six times since the 1940s and that poultry-to-human transmission is not limited to the United States.IMPORTANCEE. coli ST131 is an important extraintestinal pathogen that can colonize the gastrointestinal tracts of humans and food animals. Here, we combined detection of accessory traits associated with avian adaptation (ColV plasmids) with high-resolution phylogenetics to quantify the portion of human infections caused by ST131 strains of food animal origin. Our results suggest that one ST131 sublineage-ST131-H22-has become established in poultry populations around the world and that meat may serve as a vehicle for human exposure and infection. ST131-H22 is just one of many E. coli lineages that may be transmitted from food animals to humans. Additional studies that combine detection of host-associated accessory elements with phylogenetics may allow us to quantify the total fraction of human extraintestinal infections attributable to food animal E. coli strains.

Dataset Information

An Escherichia coli ST131 pangenome atlas reveals population structure and evolution across 4,071 isolates.

Publications

An Escherichia coli ST131 pangenome atlas reveals population structure and evolution across 4,071 isolates.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets