Anthony Underwood

Complete Nucleotide Sequences of Plasmids pEK204, pEK499, and pEK516, Encoding CTX-M Enzymes in Three Major Escherichia coli Lineages from the United Kingdom, All Belonging to the International O25:H4-ST131 Clone

ABSTRACT We determined the complete nucleotide sequences of three plasmids that encode CTX-M extended-spectrum β-lactamases (ESBLs) in pulsed-field gel electrophoresis-defined United Kingdom variants (strains A, C, and D) of the internationally prevalent Escherichia coli O25:H4-ST131 clone. Plasmid pEK499 (strain A; 117,536 bp) was a fusion of type FII and FIA replicons and harbored the following 10 antibiotic resistance genes conferring resistance to eight antibiotic classes: blaCTX-M-15, blaOXA-1, blaTEM-1,aac6′-Ib-cr, mph(A), catB4, tet(A), and the integron-borne dfrA7, aadA5, and sulI genes. pEK516 (strain D; 64,471 bp) belonged to incompatibility group IncFII and carried seven antibiotic resistance genes: blaCTX-M-15, blaOXA-1, blaTEM-1, aac6′-Ib-cr, catB4, and tet(A), all as in pEK499. It also carried aac3-IIa, conferring gentamicin resistance, and was highly related to pC15-1a, a plasmid encoding the CTX-M-15 enzyme in Canada. By contrast, pEK204 (strain C; 93,732 bp) belonged to incompatibility group IncI1 and carried only two resistance genes, blaCTX-M-3 and blaTEM-1. It probably arose by the transposition of Tn3 and ISEcp1-blaCTX-M-3 elements into a pCOLIb-P9-like plasmid. We conclude that (i) United Kingdom variants of the successful E. coli ST131 clone have acquired different plasmids encoding CTX-M ESBLs on separate occasions, (ii) the blaCTX-M-3 and blaCTX-M-15 genes on pEK204 and pEK499/pEK516 represent separate escape events, and (iii) IncFII plasmids harboring blaCTX-M-15 have played a crucial role in the global spread of CTX-M-15 ESBLs in E. coli.

Whole-genome comparison of two Acinetobacter baumannii isolates from a single patient, where resistance developed during tigecycline therapy

OBJECTIVES The whole genomes of two Acinetobacter baumannii isolates recovered from a single patient were sequenced to gain insight into the nature and extent of genomic plasticity in this important nosocomial pathogen over the course of a short infection. The first, AB210, was recovered before tigecycline therapy and was susceptible to this agent; the second, AB211, was recovered after therapy and was resistant. METHODS DNA from AB210 was sequenced by 454 GS FLX pyrosequencing according to the standard protocol for whole-genome shotgun sequencing, producing ∼250 bp fragment reads. AB211 was shotgun sequenced using the Illumina Genetic Analyzer to produce fragment reads of exactly 36 bp. Single nucleotide polymorphisms (SNPs) and large deletions detected in AB211 in relation to AB210 were confirmed by PCR and DNA sequencing. RESULTS Automated gene prediction detected 3850 putative coding sequences (CDSs). Sequence analysis demonstrated the presence of plasmids pAB0057 and pACICU2 in both isolates. Eighteen putative SNPs were detected between the pre- and post-therapy isolates, AB210 and AB211. Three contigs in AB210 were not covered by reads in AB211, representing three deletions of ∼15, 44 and 17 kb. CONCLUSIONS This study demonstrates that significant differences were detectable between two bacterial isolates recovered 1 week apart from the same patient, and reveals the potential of whole-genome sequencing as a tool for elucidating the processes responsible for changes in antibiotic susceptibility profiles.

Public Health Value of Next-Generation DNA Sequencing of Enterohemorrhagic Escherichia coli Isolates from an Outbreak

ABSTRACT In 2009, an outbreak of enterohemorrhagic Escherichia coli (EHEC) on an open farm infected 93 persons, and approximately 22% of these individuals developed hemolytic-uremic syndrome (HUS). Genome sequencing was used to investigate outbreak-derived animal and human EHEC isolates. Phylogeny based on the whole-genome sequence was used to place outbreak isolates in the context of the overall E. coli species and the O157:H7 sequence type 11 (ST11) subgroup. Four informative single nucleotide polymorphisms (SNPs) were identified and used to design an assay to type 122 other outbreak isolates. The SNP phylogeny demonstrated that the outbreak strain was from a lineage distinct from previously reported O157:H7 ST11 EHEC and was not a member of the hypervirulent clade 8. The strain harbored determinants for two Stx2 verotoxins and other putative virulence factors. When linked to the epidemiological information, the sequence data indicate that gross contamination of a single outbreak strain occurred across the farm prior to the first clinical report of HUS. The most likely explanation for these results is that a single successful strain of EHEC spread from a single introduction through the farm by clonal expansion and that contamination of the environment (including the possible colonization of several animals) led ultimately to human cases.

Evolutionary dynamics of local pandemic H1N1/2009 influenza virus lineages revealed by whole-genome analysis

ABSTRACT Virus gene sequencing and phylogenetics can be used to study the epidemiological dynamics of rapidly evolving viruses. With complete genome data, it becomes possible to identify and trace individual transmission chains of viruses such as influenza virus during the course of an epidemic. Here we sequenced 153 pandemic influenza H1N1/09 virus genomes from United Kingdom isolates from the first (127 isolates) and second (26 isolates) waves of the 2009 pandemic and used their sequences, dates of isolation, and geographical locations to infer the genetic epidemiology of the epidemic in the United Kingdom. We demonstrate that the epidemic in the United Kingdom was composed of many cocirculating lineages, among which at least 13 were exclusively or predominantly United Kingdom clusters. The estimated divergence times of two of the clusters predate the detection of pandemic H1N1/09 virus in the United Kingdom, suggesting that the pandemic H1N1/09 virus was already circulating in the United Kingdom before the first clinical case. Crucially, three clusters contain isolates from the second wave of infections in the United Kingdom, two of which represent chains of transmission that appear to have persisted within the United Kingdom between the first and second waves. This demonstrates that whole-genome analysis can track in fine detail the behavior of individual influenza virus lineages during the course of a single epidemic or pandemic.

Whole-Genome Sequencing for National Surveillance of Shiga Toxin–Producing Escherichia coli O157

BACKGROUND National surveillance of gastrointestinal pathogens, such as Shiga toxin-producing Escherichia coli O157 (STEC O157), is key to rapidly identifying linked cases in the distributed food network to facilitate public health interventions. In this study, we used whole-genome sequencing (WGS) as a tool to inform national surveillance of STEC O157 in terms of identifying linked cases and clusters and guiding epidemiological investigation. METHODS We retrospectively analyzed 334 isolates randomly sampled from 1002 strains of STEC O157 received by the Gastrointestinal Bacteria Reference Unit at Public Health England, Colindale, in 2012. The genetic distance between each isolate, as estimated by WGS, was calculated and phylogenetic methods were used to place strains in an evolutionary context. RESULTS Estimates of linked clusters representing STEC O157 outbreaks in England and Wales increased by 2-fold when WGS was used instead of traditional typing techniques. The previously unidentified clusters were often widely geographically distributed and small in size. Phylogenetic analysis facilitated identification of temporally distinct cases sharing common exposures and delineating those that shared epidemiological and temporal links. Comparison with multi locus variable number tandem repeat analysis (MLVA) showed that although MLVA is as sensitive as WGS, WGS provides a more timely resolution to outbreak clustering. CONCLUSIONS WGS has come of age as a molecular typing tool to inform national surveillance of STEC O157; it can be used in real time to provide the highest strain-level resolution for outbreak investigation. WGS allows linked cases to be identified with unprecedented specificity and sensitivity that will facilitate targeted and appropriate public health investigations.

The Impact of Pandemic Influenza H1N1 on Health-Related Quality of Life: A Prospective Population-Based Study

Background While the H1N1v influenza pandemic in 2009 was clinically mild, with a low case-fatality rate, the overall disease burden measured in quality-adjusted life years (QALY) lost has not been estimated. Such a measure would allow comparison with other diseases and assessment of the cost-effectiveness of pandemic control measures. Methods and Findings Cases of H1N1v confirmed by polymerase chain reaction (PCR) and PCR negative cases with similar influenza-like illness (ILI controls) in 7 regions of England were sent two questionnaires, one within a week of symptom onset and one two weeks later, requesting information on duration of illness, work loss and antiviral use together with EQ-5D questionnaires. Results were compared with those for seasonal influenza from a systematic literature review. A total QALY loss for the 2009 pandemic in England was calculated based on the estimated total clinical cases and reported deaths. A total of 655 questionnaires were sent and 296 (45%) returned. Symptoms and average illness duration were similar between confirmed cases and ILI controls (8.8 days and 8.7 days respectively). Days off work were greater for cases than ILI controls (7.3 and 4.9 days respectively, p = 0.003). The quality-adjusted life days lost was 2.92 for confirmed cases and 2.74 for ILI controls, with a reduction in QALY loss after prompt use of antivirals in confirmed cases. The overall QALY loss in the pandemic was estimated at 28,126 QALYs (22,267 discounted) of which 40% was due to deaths (24% with discounting). Conclusion Given the global public health significance of influenza, it is remarkable that no previous prospective study of the QALY loss of influenza using standardised and well validated methods has been performed. Although the QALY loss was minor for individual patients, the estimated total burden of influenza over the pandemic was substantial when compared to other infectious diseases.

Rapid Identification of Major Escherichia coli Sequence Types Causing Urinary Tract and Bloodstream Infections

ABSTRACT Escherichia coli sequence types (STs) 69, 73, 95, and 131 are collectively responsible for a large proportion of E. coli urinary tract and bloodstream infections, and they differ markedly in their antibiotic susceptibilities. Here, we describe a novel PCR method to rapidly detect and distinguish these lineages. Three hundred eighteen published E. coli genomes were compared in order to identify signature sequences unique to each of the four major STs. The specificities of these sequences were assessed in silico by seeking them in an additional 98 genomes. A PCR assay was designed to amplify size-distinguishable fragments unique to the four lineages and was validated using 515 E. coli isolates of known STs. Genome comparisons identified 22 regions ranging in size from 335 bp to 26.5 kb that are unique to one or more of the four predominant E. coli STs, with two to 10 specific regions per ST. These regions predominantly harbor genes encoding hypothetical proteins and are within or adjacent to prophage sequences. Most (13/22) were highly conserved (>96.5% identity) in the genomes of their respective ST. The new assay correctly identified all 142 representatives of the four major STs in the validation set (n = 515), with only two ST12 isolates misidentified as ST95. Compared with MLST, the assay has 100% sensitivity and 99.5% specificity. The rapid identification of major extraintestinal E. coli STs will benefit future epidemiological studies and could be developed to tailor antibiotic therapy to the different susceptibilities of these dominant lineages.

Comparison of the Legionella pneumophila population structure as determined by sequence-based typing and whole genome sequencing

BackgroundLegionella pneumophila is an opportunistic pathogen of humans where the source of infection is usually from contaminated man-made water systems. When an outbreak of Legionnaires’ disease caused by L. pneumophila occurs, it is necessary to discover the source of infection. A seven allele sequence-based typing scheme (SBT) has been very successful in providing the means to attribute outbreaks of L. pneumophila to a particular source or sources. Particular sequence types described by this scheme are known to exhibit specific phenotypes. For instance some types are seen often in clinical cases but are rarely isolated from the environment and vice versa. Of those causing human disease some types are thought to be more likely to cause more severe disease. It is possible that the genetic basis for these differences are vertically inherited and associated with particular genetic lineages within the population. In order to provide a framework within which to test this hypothesis and others relating to the population biology of L. pneumophila, a set of genomes covering the known diversity of the organism is required.ResultsFirstly, this study describes a means to group L. pneumophila strains into pragmatic clusters, using a methodology that takes into consideration the genetic forces operating on the population. These clusters can be used as a standardised nomenclature, so those wishing to describe a group of strains can do so. Secondly, the clusters generated from the first part of the study were used to select strains rationally for whole genome sequencing (WGS). The data generated was used to compare phylogenies derived from SBT and WGS. In general the SBT sequence type (ST) accurately reflects the whole genome-based genotype. Where there are exceptions and recombination has resulted in the ST no longer reflecting the genetic lineage described by the whole genome sequence, the clustering technique employed detects these sequence types as being admixed, indicating their mixed inheritance.ConclusionsWe conclude that SBT is usually a good proxy for the genetic lineage described by the whole genome, and therefore utility of SBT is still suitable until the technology and economics of high throughput sequencing reach the point where routine WGS of L. pneumophila isolates for outbreak investigation is feasible.

Clusters of genetically similar isolates of Pseudomonas aeruginosa from multiple hospitals in the UK

Variable number tandem repeat (VNTR) analysis at nine loci of isolates of Pseudomonas aeruginosa submitted to the national reference laboratory from UK hospitals, from over 2000 patients, between June 2010 and June 2012 revealed four widely found types that collectively were received from approximately a fifth of patients, including from those with cystic fibrosis. These types were also prevalent among related submissions from the clinical environment and were received from up to 54 (out of 143) hospitals. Multi-locus sequence typing and blaOXA-50-like sequencing confirmed the clonal relationship within each cluster, and representatives from multiple centres clustered within about 70 % by pulsed-field gel electrophoresis. Illumina sequencing of 12 isolates of cluster A of VNTR profile 8, 3, 4, 5, 2, 3, 5, 2, x (where the repeat number at the last, most discriminatory locus is variable) revealed a large number of variably present targets in the accessory genome and seven of these were sought by PCR among a larger set of isolates. Representatives from patients within a single centre mostly had distinct accessory gene profiles, suggesting that these patients acquired the strain independently, while those with clear epidemiological links shared the same profile. Profiles also varied between representatives from different centres. Epidemiological investigations of widely found types such as these require the use of finer-typing methods, which increasingly will be informed by next generation sequencing.

Evolutionary pathways of the pandemic influenza A (H1N1) 2009 in the UK.

The emergence of the influenza (H1N1) 2009 virus provided a unique opportunity to study the evolution of a pandemic virus following its introduction into the human population. Virological and clinical surveillance in the UK were comprehensive during the first and second waves of the pandemic in 2009, with extensive laboratory confirmation of infection allowing a detailed sampling of representative circulating viruses. We sequenced the complete coding region of the haemagglutinin (HA) segment of 685 H1N1 pandemic viruses selected without bias during two waves of pandemic in the UK (April-December 2009). Phylogenetic analysis showed that although temporal accumulation of amino acid changes was observed in the HA sequences, the overall diversity was less than that typically seen for seasonal influenza A H1N1 or H3N2. There was co-circulation of multiple variants as characterised by signature amino acid changes in the HA. A specific substitution (S203T) became predominant both in UK and global isolates. No antigenic drift occurred during 2009 as viruses with greater than four-fold reduction in their haemagglutination inhibition (HI) titre (“low reactors”) were detected in a low proportion (3%) and occurred sporadically. Although some limited antigenic divergence in viruses with four-fold reduction in HI titre might be related to the presence of 203T, additional studies are needed to test this hypothesis.