A typical prokaryote population sequencing study can now consist of hundreds or thousands of isolates. Interrogating these datasets can provide detailed insights into the genetic structure of ...prokaryotic genomes. We introduce Roary, a tool that rapidly builds large-scale pan genomes, identifying the core and accessory genes. Roary makes construction of the pan genome of thousands of prokaryote samples possible on a standard desktop without compromising on the accuracy of results. Using a single CPU Roary can produce a pan genome consisting of 1000 isolates in 4.5 hours using 13 GB of RAM, with further speedups possible using multiple processors.
Roary is implemented in Perl and is freely available under an open source GPLv3 license from http://sanger-pathogens.github.io/Roary
roary@sanger.ac.uk
Supplementary data are available at Bioinformatics online.
Pathogenic species in the Yersinia genus have historically been targets for research aimed at understanding how bacteria evolve into mammalian pathogens. The advent of large-scale population genomic ...studies has greatly accelerated the progress in this field, and Yersinia pestis, Yersinia pseudotuberculosis and Yersinia enterocolitica have once again acted as model organisms to help shape our understanding of the evolutionary processes involved in pathogenesis. In this Review, we highlight the gene gain, gene loss and genome rearrangement events that have been identified by genomic studies in pathogenic Yersinia species, and we discuss how these findings are changing our understanding of pathogen evolution. Finally, as these traits are also found in the genomes of other species in the Enterobacteriaceae, we suggest that they provide a blueprint for the evolution of enteropathogenic bacteria.
Vancomycin-resistant Enterococcus faecium (VREfm) is an important cause of healthcare-associated infections worldwide. We undertook whole-genome sequencing (WGS) of 495 E. faecium bloodstream ...isolates from 2001-2011 in the United Kingdom and Ireland (UK&I) and 11 E. faecium isolates from a reference collection. Comparison between WGS and multilocus sequence typing (MLST) identified major discrepancies for 17% of isolates, with multiple instances of the same sequence type (ST) being located in genetically distant positions in the WGS tree. This confirms that WGS is superior to MLST for evolutionary analyses and is more accurate than current typing methods used during outbreak investigations. E. faecium has been categorized as belonging to three clades (Clades A1, hospital-associated; A2, animal-associated; and B, community-associated). Phylogenetic analysis of our isolates replicated the distinction between Clade A (97% of isolates) and Clade B but did not support the subdivision of Clade A into Clade A1 and A2. Phylogeographic analyses revealed that Clade A had been introduced multiple times into each hospital referral network or country, indicating frequent movement of E. faecium between regions that rarely share hospital patients. Numerous genetic clusters contained highly related vanA-positive and -negative E. faecium, which implies that control of vancomycin-resistant enterococci (VRE) in hospitals also requires consideration of vancomycin-susceptible E. faecium Our findings reveal the evolution and dissemination of hospital-associated E. faecium in the UK&I and provide evidence for WGS as an instrument for infection control.
Background. Vancomycin-resistant Enterococcus faecium (VREfm) is a leading cause of nosocomial infection. Here, we describe the utility of whole-genome sequencing in defining nosocomial VREfm ...transmission. Methods. A retrospective study at a single hospital in the United Kingdom identified 342 patients with E. faecium bloodstream infection over 7 years. Of these, 293 patients had a stored isolate and formed the basis for the study. The first stored isolate from each case was sequenced (200 VREfm 197 vanA, 2 vanB, and 1 isolate containing both vanA and vanB, 93 vancomycin-susceptible E. faecium) and epidemiological data were collected. Genomes were also available for E. faecium associated with bloodstream infections in 15 patients in neighboring hospitals, and 456 patients across the United Kingdom and Ireland. Results. The majority of infections in the 293 patients were hospital-acquired (n = 249) or healthcare-associated (n = 42). Phylogenetic analysis showed that 291 of 293 isolates resided in a hospital-associated clade that contained numerous discrete clusters of closely related isolates, indicative of multiple introductions into the hospital followed by clonal expansion associated with transmission. Fine-scale analysis of 6 exemplar phylogenetic clusters containing isolates from 93 patients (32%) identified complex transmission routes that spanned numerous wards and years, extending beyond the detection of conventional infection control. These contained both vancomycin-resistant and -susceptible isolates. We also identified closely related isolates from patients at Cambridge University Hospitals NHS Foundation Trust and regional and national hospitals, suggesting interhospital transmission. Conclusions. These findings provide important insights for infection control practice and signpost areas for interventions. We conclude that sequencing represents a powerful tool for the enhanced surveillance and control of nosocomial E. faecium transmission and infection.
To give an update on the molecular epidemiology and global distribution of carbapenemase encoding genes, we subjected 313 carbapenem-resistant
isolated from 114 study centers in 47 countries in five ...world regions, Africa, Asia, Europe, Latin America, and North America, to whole genome sequencing. Numbers of isolates investigated were proportional to the population size of the contributing countries. Molecular epidemiology was investigated using seven-loci and core genome multilocus sequence typing, whole-genome single nucleotide polymorphism phylogenies, and the intrinsic
variant. Carbapenemase encoding genes were identified by multiplex PCR and ResFinder. Among the total of 313 isolates, 289 (92.3%) were assigned to
international clones (IC) IC1-IC8. IC2 predominated with 196 isolates (62.6%) and was spread worldwide, followed by IC5 with 44 isolates (14.1%) mainly confined to Latin America. Six isolates (1.9%) originating from Belgium, Egypt, Italy, and Pakistan represent the novel IC9. Acquired OXA-type carbapenemase genes were found in 300 (96%) isolates with
and
predominating, which constitutes a significant increase compared to our findings from 2010. Metallo-beta-lactamases were rare with seven isolates (2.2%). The distribution of ICs and carbapenemase determinants can vary widely among different geographical regions. IMPORTANCE Carbapenem-resistant
are of increasing public health importance, as they are resistant to last-line antibiotics. International clones with well-characterized resistance genes dominate globally; however, locally, other lineages with different properties may be of importance to consider. This study investigated isolates from a broad geographic origin from 114 hospitals in 47 countries and from five world regions ensuring the greatest possible diversity in an organism known for its propensity for clonal epidemic spread and reflecting the current global epidemiology of carbapenem-resistant
. In Latin America, a lineage different from other geographic regions circulates, with a different resistance gene profile. This knowledge is important to adjust local infection prevention measures. In a global world with migration and increasing use of antimicrobials, multidrug-resistant bacteria will continue to adapt and challenge our healthcare systems worldwide.
Genome sequencing has provided snapshots of the transmission of methicillin-resistant
(MRSA) during suspected outbreaks in isolated hospital wards. Scale-up to populations is now required to ...establish the full potential of this technology for surveillance. We prospectively identified all individuals over a 12-month period who had at least one MRSA-positive sample processed by a routine diagnostic microbiology laboratory in the East of England, which received samples from three hospitals and 75 general practitioner (GP) practices. We sequenced at least 1 MRSA isolate from 1465 individuals (2282 MRSA isolates) and recorded epidemiological data. An integrated epidemiological and phylogenetic analysis revealed 173 transmission clusters containing between 2 and 44 cases and involving 598 people (40.8%). Of these, 118 clusters (371 people) involved hospital contacts alone, 27 clusters (72 people) involved community contacts alone, and 28 clusters (157 people) had both types of contact. Community- and hospital-associated MRSA lineages were equally capable of transmission in the community, with instances of spread in households, long-term care facilities, and GP practices. Our study provides a comprehensive picture of MRSA transmission in a sampled population of 1465 people and suggests the need to review existing infection control policy and practice.
Abstract Carbapenem-resistant Klebsiella pneumoniae (CRKP) are of particular concern due to the spread of antibiotic resistance genes associated with mobile genetic elements. In this study, we ...collected 687 carbapenem-resistant strains recovered among clinical samples from 41 hospitals in nine Southern European countries (2016-2018). We identified 11 major clonal lineages, with most isolates belonging to the high-risk clones ST258/512, ST101, ST11, and ST307. bla KPC-like was the most prevalent carbapenemase-encoding gene (46%), with bla OXA-48 present in 39% of isolates. Through the combination and comparison of this EURECA collection with the previous EuSCAPE collection (2013-2014), we investigated the spread of high-risk clones circulating in Europe exhibiting regional differences. We particularly found bla KPC-like ST258/512 in Greece, Italy, and Spain, bla OXA-48 ST101 in Serbia and Romania, bla NDM ST11 in Greece, and bla OXA-48-like ST14 in Türkiye. Genomic surveillance across Europe thus provides crucial insights for local risk mapping and informs necessary adaptions for implementation of control strategies.
•This article reports the clonal expansion of an emerging ST1 colistin-resistant Acinetobacter baumannii.•Identical colistin-susceptible and colistin-resistant isolates were obtained from the same ...hospital.•This is the first report on colistin-resistant clonal Acinetobacter baumannii in South Africa.
To describe colistin-resistant Acinetobacter baumannii isolates in Cape Town, South Africa.
A. baumannii isolates identified on Vitek 2 Advanced Expert System were collected from Tygerberg Hospital referral laboratory between 2016 and 2017. Colistin resistance was confirmed using broth microdilution and SensiTest. mcr-1–5 were detected using PCR and strain typing was performed by rep-PCR. Whole genome sequencing (WGS) was performed on a subset of isolates to identify chromosomal colistin resistance mechanisms and strain diversity using multilocus sequence typing (MLST) and pairwise single nucleotide polymorphism analyses.
Twenty-six colistin-resistant and six colistin-susceptible A. baumannii were collected separately based on Vitek susceptibility; 20/26 (77%) were confirmed colistin-resistant by broth microdilution. Four colistin-resistant isolates were isolated in 2016 and 16 in 2017, from five healthcare facilities. Thirteen colistin-resistant isolates and eight colistin-susceptible isolates were identical by rep-PCR and MLST (ST1), all from patients admitted to a tertiary hospital during 2017. The remaining colistin-resistant isolates were unrelated.
An increase in colistin-resistant A. baumannii isolates from a tertiary hospital in 2017 appears to be clonal expansion of an emerging colistin-resistant strain. This strain was not detected in 2016 or from other hospitals. Identical colistin-susceptible isolates were also isolated, suggesting relatively recent acquisition of colistin resistance.
Tackling multidrug-resistant Escherichia coli requires evidence from One Health studies that capture numerous potential reservoirs in circumscribed geographic areas.
We conducted a survey of extended ...β-lactamase (ESBL)-producing E. coli isolated from patients, canals and livestock wastewater in eastern Thailand between 2014 and 2015, and analyzed isolates using whole genome sequencing.
The bacterial collection of 149 isolates consisted of 84 isolates from a single hospital and 65 from the hospital sewer, canals and farm wastewater within a 20 km radius. E. coli ST131 predominated the clinical collection (28.6%), but was uncommon in the environment. Genome-based comparison of E. coli from infected patients and their immediate environment indicated low genetic similarity overall between the two, although three clinical-environmental isolate pairs differed by ≤ 5 single nucleotide polymorphisms. Thai E. coli isolates were dispersed throughout a phylogenetic tree containing a global E. coli collection. All Thai ESBL-positive E. coli isolates were multidrug resistant, including high rates of resistance to tobramycin (77.2%), gentamicin (77.2%), ciprofloxacin (67.8%) and trimethoprim (68.5%). ESBL was encoded by six different CTX-M elements and SHV-12. Three isolates from clinical samples (n = 2) or a hospital sewer (n = 1) were resistant to the carbapenem drugs (encoded by NDM-1, NDM-5 or GES-5), and three isolates (clinical (n = 1) and canal water (n = 2)) were resistant to colistin (encoded by mcr-1); no isolates were resistant to both carbapenems and colistin.
Tackling ESBL-producing E. coli in this setting will be challenging based on widespread distribution, but the low prevalence of resistance to carbapenems and colistin suggests that efforts are now required to prevent these from becoming ubiquitous.
The genus Yersinia has been used as a model system to study pathogen evolution. Using whole-genome sequencing of all Yersinia species, we delineate the gene complement of the whole genus and define ...patterns of virulence evolution. Multiple distinct ecological specializations appear to have split pathogenic strains from environmental, nonpathogenic lineages. This split demonstrates that contrary to hypotheses that all pathogenic Yersinia species share a recent common pathogenic ancestor, they have evolved independently but followed parallel evolutionary paths in acquiring the same virulence determinants as well as becoming progressively more limited metabolically. Shared virulence determinants are limited to the virulence plasmid pYV and the attachment invasion locus ail . These acquisitions, together with genomic variations in metabolic pathways, have resulted in the parallel emergence of related pathogens displaying an increasingly specialized lifestyle with a spectrum of virulence potential, an emerging theme in the evolution of other important human pathogens.