A population phylogenomic analysis of the origin and spread of Escherichia coli sequence type 131 (ST131)

The incidence of infections caused by extraintestinal Escherichia coli (ExPEC) is rising globally due to their increasing resistance to standard antibiotics. This results in the use of broader-spectrum drugs, prolonged patient ill-health and more nosocomial infections. E. coli sequence type 131 (ST131) is the predominant ExPEC clone worldwide. The antimicrobial resistance (AMR) gene repertoire of ST131 is evolving rapidly due to the widespread use of β-lactam (bla) antibiotics. Here, we performed a genomic investigation of an ST131 outbreak in a long-term care facility (LTCF) to describe transmission, within-host clonal diversity, genetic diversity of antibiotic resistance and the evolution of ST131 in the LTCF over a seven-year period. We analyzed the population structure and inferred the genealogical history of the LTCF isolates in the context of local hospital and global collections of ST131 to elucidate the epidemiology of ST131. We confirmed our initial hypotheses by reconstructing the evolutionary history of a much larger population consisting of >4000 global ST131 genomes This provided a deeper resolution of their evolutionary trajectories and the adaptive mechanisms of AMR driven by their ESBL genes, particularly cefotaximase (blaCTX). We further investigated the intersection of the AMR genes (AMRGs) found in ST131 with that of the human microbiome to understand the extent of their loss, gain and spread across different bacterial species. Across all strains, a large number of ST131’s AMRGs were found in a total of 794 genes in the human microbiome. Various gene families were represented, including transporters, transcription factors, β-lactamases and cell wall biosynthesis enzymes. To establish the main culprit for the dynamic nature of the blaCTX-M genes, we performed long read sequencing using a GridION X5 instrument. Analysis of long read-only assemblies revealed a clear and robust result on the genetic flanking context of blaCTX-M genes in both plasmid and chromosomes. Overall, our findings underpin the tremendous potential power for improving our current treatment of bacterial infections using high-throughput analysis of whole genome sequence data

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

This work was funded by a DCU O’Hare Ph.D. fellowship and a DCU Enhancing Performance grant.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.Publisher PDFPeer reviewe

Complete Assembly of Escherichia coli Sequence Type 131 Genomes Using Long Reads Demonstrates Antibiotic Resistance Gene Variation within Diverse Plasmid and Chromosomal Contexts.

The incidence of infections caused by extraintestinal Escherichia coli (ExPEC) is rising globally, which is a major public health concern. ExPEC strains that are resistant to antimicrobials have been associated with excess mortality, prolonged hospital stays, and higher health care costs. E. coli sequence type 131 (ST131) is a major ExPEC clonal group worldwide, with variable plasmid composition, and has an array of genes enabling antimicrobial resistance (AMR). ST131 isolates frequently encode the AMR genes blaCTX-M-14, blaCTX-M-15, and blaCTX-M-27, which are often rearranged, amplified, and translocated by mobile genetic elements (MGEs). Short DNA reads do not fully resolve the architecture of repetitive elements on plasmids to allow MGE structures encoding blaCTX-M genes to be fully determined. Here, we performed long-read sequencing to decipher the genome structures of six E. coli ST131 isolates from six patients. Most long-read assemblies generated entire chromosomes and plasmids as single contigs, in contrast to more fragmented assemblies created with short reads alone. The long-read assemblies highlighted diverse accessory genomes with blaCTX-M-15, blaCTX-M-14, and blaCTX-M-27 genes identified in three, one, and one isolates, respectively. One sample had no blaCTX-M gene. Two samples had chromosomal blaCTX-M-14 and blaCTX-M-15 genes, and the latter was at three distinct locations, likely transposed by the adjacent MGEs: ISEcp1, IS903B, and Tn2 This study showed that AMR genes exist in multiple different chromosomal and plasmid contexts, even between closely related isolates within a clonal group such as E. coli ST131.IMPORTANCE Drug-resistant bacteria are a major cause of illness worldwide, and a specific subtype called Escherichia coli ST131 causes a significant number of these infections. ST131 bacteria become resistant to treatments by modifying their DNA and by transferring genes among one another via large packages of genes called plasmids, like a game of pass-the-parcel. Tackling infections more effectively requires a better understanding of what plasmids are being exchanged and their exact contents. To achieve this, we applied new high-resolution DNA sequencing technology to six ST131 samples from infected patients and compared the output to that of an existing approach. A combination of methods shows that drug resistance genes on plasmids are highly mobile because they can jump into ST131's chromosomes. We found that the plasmids are very elastic and undergo extensive rearrangements even in closely related samples. This application of DNA sequencing technologies illustrates at a new level the highly dynamic nature of ST131 genomes

Whole genome sequencing and molecular epidemiology of clinical isolates of Staphylococcus aureus from Algeria

Funding: This research was funded by the Algerian Ministry of Higher Education and Scientific Research (DGRSDT/MESRS), and by grants from the Wellcome Trust (098731/z/11/z to St Andrews University Bioinformatics Unit) and to the Chief Scientists Office (SIRN10 to M.T.G.H.), and the National Institute for Health Research, Medical Research Council and the Department of Health and Social Care (MR/S004785/1 to M.T.G.H.); this award is also part of the EDCTP2 program supported by the European Union.Staphylococcus aureus is an important pathogen responsible for various healthcare- and community-acquired infections. In this study, whole genome sequencing (WGS) was used to genotype S. aureus clinical isolates from two hospitals in Algeria and to characterize their genetic determinants of antimicrobial resistance. Seventeen S. aureus isolates were included in this study. WGS, single-nucleotide polymorphism (SNP)-based phylogenetic analysis, in silico multilocus sequence typing (MLST), spa and staphylococcal cassette chromosome mec (SCCmec) typing and in silico antimicrobial resistance profiling were performed. Phenotypic antibiotic susceptibility testing was performed using the Vitek 2 system and the disk diffusion method. The isolates were separated into sequence types (STs), with ST80 being predominant; five clonal complexes (CCs); four spa types (t044, t127, t368, t386); and two SCCmec types (IVc and IVa). Whole genome analysis revealed the presence of the resistance genes mecA, blaZ, ermC, fusB, fusC, tetK, aph(3′)-IIIa and aad(6) and mutations conferring resistance in the genes parC and fusA. The rate of multidrug resistance (MDR) was 64%. This work provides a high-resolution characterization of methicillin-resistant S. aureus (MRSA) and methicillin-susceptible S. aureus (MSSA) isolates and emphasizes the importance of continuous surveillance to monitor the spread of S. aureus in healthcare settings in the country.Publisher PDFPeer reviewe

Pan-resistome characterization of uropathogenic Escherichia coli and Klebsiella pneumoniae strains circulating in Uganda and Kenya, isolated from 2017-2018

Funding: The Holistic Approach to Unravel Antibacterial Resistance in East Africa is a 3-year Global Context Consortia Award (MR/S004785/1) funded by the National Institute for Health Research, Medical Research Council and the Department of Health and Social Care. The award is also part of the EDCTP2 program supported by the European Union.Urinary tract infection (UTI) develops after a pathogen adheres to the inner lining of the urinary tract. Cases of UTIs are predominantly caused by several Gram-negative bacteria and account for high morbidity in the clinical and community settings. Of greater concern are the strains carrying antimicrobial resistance (AMR)-conferring genes. The gravity of a UTI is also determined by a spectrum of other virulence factors. This study represents a pilot project to investigate the burden of AMR among uropathogens in East Africa. We examined bacterial samples isolated in 2017–2018 from in- and out-patients in Kenya (KY) and Uganda (UG) that presented with clinical symptoms of UTI. We reconstructed the evolutionary history of the strains, investigated their population structure, and performed comparative analysis their pangenome contents. We found 55 Escherichia coli and 19 Klebsiella pneumoniae strains confirmed uropathogenic following screening for the prevalence of UTI virulence genes including fimH, iutA, feoA/B/C, mrkD, and foc. We identified 18 different sequence types in E. coli population while all K. pneumoniae strains belong to ST11. The most prevalent E. coli sequence types were ST131 (26%), ST335/1193 (10%), and ST10 (6%). Diverse plasmid types were observed in both collections such as Incompatibility (IncF/IncH/IncQ1/IncX4) and Col groups. Pangenome analysis of each set revealed a total of 2862 and 3464 genes comprised the core genome of E. coli and K. pneumoniae population, respectively. Among these are acquired AMR determinants including fluoroquinolone resistance-conferring genes aac(3)-Ib-cr and other significant genes: aad, tet, sul1, sul2, and cat, which are associated with aminoglycoside, tetracycline, sulfonamide, and chloramphenicol resistance, respectively. Accessory genomes of both species collections were detected several β-lactamase genes, blaCTX-M, blaTEM and blaOXA, or blaNDM. Overall, 93% are multi-drug resistant in the E. coli collection while 100% of the K. pneumoniae strains contained genes that are associated with resistance to three or more antibiotic classes. Our findings illustrate the abundant acquired resistome and virulome repertoire in uropathogenic E. coli and K. pneumoniae, which are mainly disseminated via clonal and horizontal transfer, circulating in the East African region. We further demonstrate here that routine genomic surveillance is necessary for high-resolution bacterial epidemiology of these important AMR pathogens.Publisher PDFPeer reviewe

Evaluation of sequence hybridisation for respiratory viruses using the Twist Biosciences Respiratory panel and the OneCodex Respiratory Virus sequence analysis workflow

Respiratory viral infections are a major global clinical problem, and rapid, cheap, scalable and agnostic diagnostic tests that capture genome-level information on viral variation are urgently needed. Metagenomic approaches would be ideal, but remain currently limited in that much of the genetic content in respiratory samples is human, and amplifying and sequencing the viral/pathogen component in an unbiased manner is challenging. PCR-based tests, including those which detect multiple pathogens, are already widely used, but do not capture information on strain-level variation; tests with larger viral repertoires are also expensive on a per-test basis. One intermediate approach is the use of large panels of viral probes or ‘baits’, which target or ‘capture’ sequences representing complete genomes amongst several different common viral pathogens; these are then amplified, sequenced and analysed with a sequence analysis workflow. Here we evaluate one such commercial bait capture method (the Twist Bioscience Respiratory Virus Research Panel) and sequence analysis workflow (OneCodex), using control (simulated) and patient samples head-to-head with a validated multiplex PCR clinical diagnostic test (BioFire FilmArray). We highlight the limited sensitivity and specificity of the joint Twist Bioscience/OneCodex approach, which are further reduced by shortening workflow times and increasing sample throughput to reduce per-sample costs. These issues with performance may be driven by aspects of both the laboratory (e.g. capacity to enrich for viruses present in low numbers), bioinformatics methods used (e.g. a limited viral reference database) and thresholds adopted for calling a virus as present or absent. As a result, this workflow would require further optimization prior to any implementation for respiratory virus characterization in a routine diagnostic healthcare setting

Molecular characterizations of the coagulase-negative staphylococci species causing urinary tract infection in Tanzania : a laboratory-based cross-sectional study

Funding: This study is part of the Holistic Approach to Unravel Antibacterial Resistance in East Africa (HATUA) project funded by the National Institute for Health Research, Medical Research Council and the Department of Health and Social Care, Award (MR/S004785/1).Background: There is a growing body of evidence on the potential involvement of coagulase-negative Staphylococci (CoNS) in causing urinary tract infections (UTIs). The aim of this study was to delineate virulence potential, antimicrobial resistance genes, and sequence types of CoNS isolated from patients with UTI symptoms and pyuria in Tanzania. Methods: CoNS from patients with UTI symptoms and more than 125 leucocytes/μL were retrieved, subcultured, and whole-genome sequenced. Results: Out of 65 CoNS isolates, 8 species of CoNS were identified; Staphylococcus haemolyticus, n = 27 (41.5%), and Staphylococcus epidermidis, n = 24 (36.9%), were predominant. The majority of S. haemolyticus were sequence type (ST) 30, with 8 new ST138-145 reported, while the majority of S. epidermidis were typed as ST490 with 7 new ST1184-1190 reported. Sixty isolates (92.3%) had either one or multiple antimicrobial resistance genes. The most frequently detected resistance genes were 53 (21%) dfrG, 32 (12.9%) blaZ, and 26 (10.5%) mecA genes conferring resistance to trimethoprim, penicillin, and methicillin, respectively. Out of 65 isolates, 59 (90.8%) had virulence genes associated with UTI, with a predominance of the icaC 47 (46.5%) and icaA 14 (13.9%) genes. Conclusion: S. haemolyticus and S. epidermidis harboring icaC, dfrG, blaZ, and mecA genes were the predominant CoNS causing UTI in Tanzania. Laboratories should carefully interpret the significant bacteriuria due to CoNS in relation to UTI symptoms and pyuria before labeling them as contaminants. Follow-up studies to document the outcome of the treated patients is needed to add more evidence that CoNS are UTI pathogens.Publisher PDFPeer reviewe

A population phylogenomic analysis of the origin and spread of Escherichia coli sequence type 131 (ST131)

The incidence of infections caused by extraintestinal Escherichia coli (ExPEC) is rising globally due to their increasing resistance to standard antibiotics. This results in the use of broader-spectrum drugs, prolonged patient ill-health and more nosocomial infections. E. coli sequence type 131 (ST131) is the predominant ExPEC clone worldwide. The antimicrobial resistance (AMR) gene repertoire of ST131 is evolving rapidly due to the widespread use of β-lactam (bla) antibiotics. Here, we performed a genomic investigation of an ST131 outbreak in a long-term care facility (LTCF) to describe transmission, within-host clonal diversity, genetic diversity of antibiotic resistance and the evolution of ST131 in the LTCF over a seven-year period. We analyzed the population structure and inferred the genealogical history of the LTCF isolates in the context of local hospital and global collections of ST131 to elucidate the epidemiology of ST131. We confirmed our initial hypotheses by reconstructing the evolutionary history of a much larger population consisting of >4000 global ST131 genomes This provided a deeper resolution of their evolutionary trajectories and the adaptive mechanisms of AMR driven by their ESBL genes, particularly cefotaximase (blaCTX). We further investigated the intersection of the AMR genes (AMRGs) found in ST131 with that of the human microbiome to understand the extent of their loss, gain and spread across different bacterial species. Across all strains, a large number of ST131’s AMRGs were found in a total of 794 genes in the human microbiome. Various gene families were represented, including transporters, transcription factors, β-lactamases and cell wall biosynthesis enzymes. To establish the main culprit for the dynamic nature of the blaCTX-M genes, we performed long read sequencing using a GridION X5 instrument. Analysis of long read-only assemblies revealed a clear and robust result on the genetic flanking context of blaCTX-M genes in both plasmid and chromosomes. Overall, our findings underpin the tremendous potential power for improving our current treatment of bacterial infections using high-throughput analysis of whole genome sequence data

An <i>Escherichia coli </i>ST131 pangenome atlas reveals population structure and evolution across 4,071 isolates

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