- Research
- Open access
- Published:
Genomic epidemiology of a protracted hospital outbreak caused by multidrug-resistant Acinetobacter baumanniiin Birmingham, England
Genome Medicine volume 6, Article number: 70 (2014)
Abstract
Background
Multidrug-resistant Acinetobacter baumannii commonly causes hospital outbreaks. However, within an outbreak, it can be difficult to identify the routes of cross-infection rapidly and accurately enough to inform infection control. Here, we describe a protracted hospital outbreak of multidrug-resistant A. baumannii, in which whole-genome sequencing (WGS) was used to obtain a high-resolution view of the relationships between isolates.
Methods
To delineate and investigate the outbreak, we attempted to genome-sequence 114 isolates that had been assigned to the A. baumannii complex by the Vitek2 system and obtained informative draft genome sequences from 102 of them. Genomes were mapped against an outbreak reference sequence to identify single nucleotide variants (SNVs).
Results
We found that the pulsotype 27 outbreak strain was distinct from all other genome-sequenced strains. Seventy-four isolates from 49 patients could be assigned to the pulsotype 27 outbreak on the basis of genomic similarity, while WGS allowed 18 isolates to be ruled out of the outbreak. Among the pulsotype 27 outbreak isolates, we identified 31 SNVs and seven major genotypic clusters. In two patients, we documented within-host diversity, including mixtures of unrelated strains and within-strain clouds of SNV diversity. By combining WGS and epidemiological data, we reconstructed potential transmission events that linked all but 10 of the patients and confirmed links between clinical and environmental isolates. Identification of a contaminated bed and a burns theatre as sources of transmission led to enhanced environmental decontamination procedures.
Conclusions
WGS is now poised to make an impact on hospital infection prevention and control, delivering cost-effective identification of routes of infection within a clinically relevant timeframe and allowing infection control teams to track, and even prevent, the spread of drug-resistant hospital pathogens.
Background
Acinetobacter baumannii is an important cause of nosocomial infection, particularly ventilator-associated pneumonia and bloodstream infections in critically ill patients, and has a tendency to cause hospital outbreaks [1],[2]. Multidrug-resistant (MDR) and even pan-drug-resistant strains have been reported worldwide [3]. It has also emerged as a threat to casualties of the conflicts in Iraq and Afghanistan, with the secondary problem that strains introduced to hospitals by military personnel can cause cross infection of staff and patients [4]-[9]. Although existing molecular typing methods play an important role in identifying outbreaks [10],[11], they lack the resolution necessary to identify chains and modes of transmission within outbreaks and so can provide only limited guidance to infection control teams on how best to control or terminate an outbreak.
Whole-genome sequencing (WGS) of bacterial isolates provides a promising new method for investigating the epidemiology of outbreaks, particularly when coupled to clinical locational and temporal data [12]-[17]. Here, we describe a protracted hospital outbreak which occurred in Birmingham, England between July 2011 and February 2013 and was caused by a strain of Acinetobacter baumannii belonging to pulse-field gel electrophoresis type (pulsotype) 27. During the outbreak, we used genome sequencing to obtain a high-resolution view of the relationships between isolates, allowing us to reconstruct chains of transmission, confirm or refute epidemiological hypotheses and to provide the infection control team with useful insights into the sources and routes of infection during this outbreak.
Methods
Microbiological investigations
Here, we report a routine and clinically indicated infection control investigation into an outbreak, with no experimentation on human subjects. No additional samples other than those that were clinically relevant were taken from patients and the use of genome sequencing falls under the remit of laboratory method development, which does not need ethical approval. Multidrug-resistant Acinetobacter (MDR-Aci) isolates were obtained from routine clinical samples through culture on blood agar, followed by single-colony isolation. Bacterial identification and antibiotic susceptibility testing were performed in the hospital microbiology laboratory on the Vitek 2 system according to the manufacturer's recommendations (bioMérieux, Basingstoke, UK) [18]. Multidrug resistance was defined as resistance to ≥3 classes of antibiotics (quinolones, extended-spectrum cephalosporins, β-lactam/β-lactamase inhibitor combinations, aminoglycosides and carbapenems).
All MDR-Aci isolates from the Queen Elizabeth Hospital Birmingham during the outbreak period (July 2011 to February 2013) were considered for inclusion in the study. During this period, 65 patients tested positive for MDR-Aci in the clinical laboratory. Patients were numbered consecutively, based on the date of first isolation of MDR-Aci. The initial MDR-Aci isolate from each patient was sent to the Laboratory of HealthCare Associated Infection in Colindale, London for speciation and typing by pulsed-field gel electrophoresis (PFGE) and other molecular methods [10]. When the reference laboratory finds that two or more isolates from the UK share a novel PFGE pattern, the isolates are assigned to a new numerical pulsotype, for example, pulsotype 27 or pulsotype 29.
An attempt was made to propagate isolates from all MDR-Aci-positive patients for genomic analysis. However, isolates from three patients (patients 15, 28 and 38) were lost on sub-culture or contaminated, leaving us with 74 genome-sequenced pulsotype 27 isolates from 58 patients. To examine within-host diversity, multiple isolates were obtained from 13 patients from different body sites and/or at different times. In addition, 18 isolates from 15 patients that had been identified as A. baumannii complex by Vitek 2, but turned out not to belong to the outbreak, were subjected to genome analysis, as were 10 environmental isolates and four control strains, which had been subjected to prolonged subculture in the laboratory. We also genome-sequenced the first pulsotype 27 isolate from the UK (kindly supplied by Jane Turton at the Laboratory of HealthCare Associated Infection), which was recovered in 2006 from a patient that had recently undergone surgery in India.
Genomic and epidemiological investigation
Genomic DNA was extracted from 114 putative Acinetobacter isolates, applying Qiagen 100/G Genomic-tips to 5 to 10 mL of overnight culture. A barcoded fragment library was generated for each isolate using the Nextera Sample Preparation and Nextera Index Kits (Illumina), then sequenced on an Illumina MiSeq, using paired-end (2 × 151 or 2 × 251) protocols, to give a minimum depth of coverage of 10×. We implemented a filtering pipeline that trimmed reads at both ends, removing adaptors and bases with sequencing quality <20, and discarded all reads that mapped to PhiX or that contained Ns or where >20% of bases had a sequencing quality of <20.
The genome of an isolate from a patient early in the outbreak (patient 6) was sequenced on two different sequencing platforms (454 FLX +and Illumina MiSeq), then a hybrid assembly was created to provide a reference genome for the outbreak, using Newbler v2.6 [19]. This assembly consisted of 4,031,405 base pairs, with 160 contigs in total and 126 contigs >500 bp, with an N50 for contigs >500 base pairs of 31,936 base pairs. Five contigs (seq23, 67, 75, 100 and 128), comprising 77,648 base pairs/80 CDSs, were assigned to a cryptic plasmid on the basis of read depth, patterns of absence in some isolates and homology searches.
The outbreak reference genome was compared to all the MDR-Aci genome sequences that were publically available in May 2013, using the Average Nucleotide Identity (ANI) approach to identify the closest genome-sequenced strain [20],[21]. Isolates were assigned to a species on the basis of ANI to reference genomes [20],[21]. For genotypic investigations of potential outbreaks, genome sequences were mapped to the relevant reference genome using Bowtie 2 [22], with default parameters, except that the reads were soft-clipped at the ends to improve the alignment score (option --local).
SNV discovery procedure
After mapping each set of read data to the reference genome as explained above, we processed with SAMtools v0.1.18 [23] (mpileup with default parameters, disabling the probabilistic realignment for the computation of base alignment quality, that is, we used option -B) and filtered it using BCFtools v0.1.17-dev (using the vcfutils.pl varFilter script to find variants with minimum root-mean-square mapping quality of 30, maximum read depth of 10,000 and minimum distance to a gap of 150 bp, that is, approximately one read length). Using custom scripts, we screened these SNV locations to exclude some potentially spurious SNVs by retaining only SNVs which are:
-
not from SNV-dense regions - no more than three SNVs in a 1,001 bp window centred on the SNV location
-
most likely not from repeat regions - coverage less than twice the average isolate's coverage and
-
at least 150 bp from scaffold boundaries.
The alignments of the remaining variant loci were then manually inspected to check quality. For all SNV loci with coverage five-fold or less or with consensus <90%, we performed PCR and Sanger sequencing as a SNV verification step.
For isolates from the pulsotype 27 outbreak, we defined a new major SNV genotype, with a numerical designation (1.0, 2.0 and so on) whenever isolates with the same SNV profile were obtained from more than one patient. Genotypes confined to single patients within the outbreak were defined as minor genotypes and were given serial numerical designations (1.1, 2.1, 2.2 and so on) based on the closest relationship to a major genotype. All sequences from this study are available under ENA Accession number PRJEB4735.
To reconstruct the most parsimonious set of transmission events, we used a custom-built algorithm and script. We considered three modes of transmission: direct ward contact, where donor and recipient are on the same ward at the same time; delayed ward contact, where donor and recipient have been on the same ward but not at the same time; and theatre contact, where donor and recipient have received treatment in the same operating theatre. In the latter two modes, we assumed that MDR-Aci strains could survive in the hospital environment for up to 5 weeks [24].
Results
Description of the Acinetobacter baumannii pulsotype 27 outbreak and the outbreak strain
Queen Elizabeth Hospital Birmingham (QEHB) is a public hospital with around 1,200 beds. British military casualties are usually repatriated here. In July 2011, we saw the first case in a new MDR-Aci outbreak: a military patient, who had been admitted to the hospital after suffering a blast injury in Afghanistan. The outbreak strain was identified in the clinical laboratory as an MDR-Aci showing resistance to multiple classes of antimicrobial agent, including aminoglycosides, fluoroquinolones, carbapenems, cephalosporins, carboxypenicillins and tetracycline. It was variably sensitive to colistin. The reference laboratory identified it as A. baumannii and, using PFGE, assigned it to pulsotype 27, a pulsotype not previously seen at QEHB.
Over the following 80 weeks, isolates of the outbreak strain were recovered from an additional 51 patients, including civilians and military casualties, spanning multiple wards (Figure 1). Isolates were assigned to the outbreak provisionally on the basis of antibiogram and definitively on the basis of pulsotype. However, after we began to apply WGS to putative MDR-Aci isolates in week 40 of the outbreak (April 2012), we also used genome sequence information to determine which isolates belonged to the outbreak. Analysis of a reference genome for the outbreak, created from an early outbreak isolate (isolate 6), revealed that the outbreak strain was distinct from all other well-characterised strains, including strains from previous outbreaks in local hospitals and clusters with other strains in a previously defined, widespread clonal lineage, International Clone I. A comparison with the first British pulsotype 27 isolate, which had been identified by the Laboratory of HealthCare Associated Infection in 2006, showed over 1,200 SNVs, ruling out any close relationship between our 2011 to 2013 outbreak and this historical isolate.
To delineate and investigate the outbreak, we attempted to genome-sequence 114 isolates identified as A. baumannii complex by the Vitek 2 system and obtained informative draft genome sequences from 102 of them. Seventy-four clinical isolates, obtained from 52 patients (including the index case), and 10 environmental isolates showed sufficient genomic similarity to the outbreak reference strain (≤8 SNVs different) to be ruled into the outbreak (Table 1).
Genomics reveals clusters of infection
We identified 31 SNVs in outbreak isolates (Table 2). These SNVs, together with presence/absence of a cryptic plasmid (that is, a 77-kb plasmid with no obvious phenotype), define seven major outbreak genotypes, which fall into a phylogenetic relationship consistent with the timeline of the outbreak (Figure 2). Most of the major genotypes in the outbreak are accompanied by a cloud of one or two SNV variants in a wheel-and-spokes configuration, so that in total there are 32 distinct genotypes of the outbreak strain (Table 1, Figure 2). Laboratory subculture controls, including two different colonies picked after four serial subcultures and a culture that had been subjected to seven freeze-thaw cycles, all had the same SNV genotype as a minimally passaged parent culture, suggesting that SNVs are not readily acquired in the laboratory.
Among the outbreak genomes, we found two protein-coding genes that contained more than one SNV. Four non-synonymous SNVs were identified in the same sensor kinase gene, pmrB, which has been implicated in colistin resistance [25]. None of these SNVs match known colistin resistance-associated mutations. However, three of the changes occur in isolates (1b, 34e and 23) that show decreased susceptibility to colistin (MIC 256 mg/L), suggesting that they might represent new resistance-associated mutations. Isolates that share the fourth SNV in pmrB, which delineates patient 22 from Genotype 2.0 and also occurs in Genotypes 3.0-7.0, retain sensitivity to colistin (MIC ≤8 mg/L).
Non-outbreak isolates, within-host diversity and mixed infections
Genome sequencing revealed that 18 Acinetobacter isolates, although obtained from hospitalised patients during the study period, did not belong to the main outbreak. A close genomic relationship between four isolates, twinned with a 3-week overlap in ward occupancy, provided evidence of a second small outbreak of MDR-Aci, which from PFGE results on three of the isolates could be assigned to pulsotype 29 (Table 3). Similarly, genome comparisons established that seven isolates obtained from five patients were members of the related species, Acinetobacter pittii, but very large pair-wise differences (>20,000 SNVs) between A. pittii isolates from different patients ruled out cross-infection.
From one trauma patient (patient 26), who was hospitalised for over 7 months, we genome-sequenced seven isolates of MDR-Aci obtained from different anatomical sites over a 4-month period and found five SNV variants (Figure 2):
-
The initial isolate, 26a, which was obtained from a sputum sample, falls one SNV away from Genotype 4.0.
-
A blood isolate (26b) taken 8 days later falls within Genotype 4.0.
-
Isolates 26c/d/f, obtained from a series of CSF samples taken approximately 3 months later, fall one SNV away from 26a
-
A second sputum isolate (26e) represents a unique one-SNV variant of genotype 4.0.
Retrieval of a cloud of genotypes from a single patient illustrates the potential for within-host evolution in MDR-Aci, mirroring findings with other hospital pathogens such as Staphylococcus aureus[26],[27].
From yet another CSF sample from patient 26, we isolated a strain of MDR-Aci that was shown to be distinct from the outbreak strain by PFGE typing and by genome sequencing, providing evidence of double infection. We also found evidence of double infection with Acinetobacter in another trauma patient, patient 44, where two isolates, each from a separate wound swab taken on the same day, were identified by genome sequencing as A. pittii and the outbreak strain of A. baumannii.
Routes and chains of transmission within the main MDR-Aci outbreak
We reconstructed transmission events, assuming the most parsimonious transmission paths between patients. Using conventional epidemiological information alone, we identified 273 potential transmission events - an average of approximately five per patient - that might link patients within the outbreak. When genome sequence data were included, we were able to reduce this to a set of 57 potential transmission events. This set linked all but 10 of the pulsotype 27 patients and, in most cases, provided a single most-parsimonious transmission event that explained how a patient acquired the outbreak strain (Table 4).
Early in the outbreak, epidemiological and genomic analyses indicated that transmission occurred primarily as a result of cross-infection between patients located on the same ward at the same time. Thus, all isolates from Genotypes 1.0 and 2.0 and most of the isolates from Genotype 4.0 came from patients who had stayed on the Ward 1. In some cases, long-term contamination of the ward environment was thought to account for transmission and this was confirmed by environmental swabbing in side rooms after patients had been discharged and the room cleaned (Table 1). For example, isolate E1 was recovered a day after patient 44 was discharged; genomic analyses revealed it shared the same SNV profile (Genotype 6.0) as four of the five MDR-Aci isolates from that patient. Similarly isolates E2-4 were taken a day after patient 55 was discharged and were found to show a one-SNV difference from a patient 55 isolate. In both cases, the patients suffered severe burns and each stayed in a single room for the entire hospital stay. Confirmation of contamination of the hospital environment led to a tightening of ward decontamination procedures.
Some outbreak strain acquisitions could not be explained simply by within-ward transmission, so we were forced to consider alternative routes of infection. As the outbreak progressed, we noticed that most of the affected patients made numerous visits to operating theatres: only five were never treated in an operating theatre. One particular theatre, specializing in the treatment of burns patients, was implicated in transmission between patient 34 (donor) and patients 40 and 39 (recipients). Consequently, in week 46 the burns theatre was closed and underwent deep cleaning (that is, decluttering of the operating theatre, followed by cleaning of all patient-associated equipment, non-fixed items, horizontal surfaces, walls, ceilings, ventilation shafts and storage areas with a chlorine-based disinfectant). Although there were several ward-based transmission events in the weeks that followed, no new theatre-acquired cases were observed for the subsequent 6 weeks and, for a time, the outbreak appeared to have ended.
Unfortunately, the outbreak resumed when a burns patient, patient 52, presented with an isolate from Genotype 6.0 in week 70. Initial epidemiological investigations failed to find any plausible direct ward- or theatre-based route of transmission that might link patient 52 with earlier outbreak cases. However, our finding of genotypic identity between the patient 52 isolate and previous outbreak isolates forced us to perform a more thorough epidemiological investigation, which uncovered a vehicle for transmission: patient 52 had occupied a specialised burns care bed that had been previously occupied by another Genotype 6.0 patient, patient 50. This prompted the development of a decontamination protocol for this specialised type of bed.
The outbreak spread to over a dozen new patients during the subsequent 9 weeks. Our suspicion once again focused on the burns theatre as the likely source of infection. This was confirmed when we obtained six isolates (E5-10) from environmental swabs of the burns operating theatre. All isolates from this phase of the outbreak, from patients and the environment, belonged to, or were closely related to, Genotypes 6.0 and 7.0. These findings prompted a second closure of the burns theatre, with deep cleaning in week 76. Following this second deep clean of the theatre the outbreak ceased and no further acquisitions of the strain were identified. The outbreak was formally declared closed in May 2013 when no inpatients were colonised or infected with the outbreak strain and there had been no new acquisitions for a period of 12 weeks.
Discussion
Like many other hospitals, QEHB suffers from serial clonal outbreaks of MDR-Aci, which result from the importations of outbreak strains, often by military patients [7]-[9]. We have described a prolonged outbreak of MDR-Aci, in which bacterial WGS provided a powerful adjunct to conventional laboratory and epidemiological investigations. In so doing, we have built on our previous smaller-scale efforts on the genomic epidemiology of A. baumannii[7] and on the work of others on a polyclonal outbreak genome-sequenced using the 454 platform [28]. However, we have now benefited from the improved ease of use, cost-effectiveness, throughput and accuracy of the Illumina MiSeq platform. Unlike a recent epidemiological analysis of MDR-Aci in a US hospital [29]. We have focused on a single clonal outbreak and have used a whole-genome SNV-based analysis twinned with loss or gain of a plasmid to reconstruct strain evolution as the outbreak progressed.
We have shown that this particular pathogen, Acinetobacter baumannii, undergoes sufficiently rapid genome evolution within a hospital outbreak to allow SNV analyses to reveal modes and routes of transmission and distinguish between alternate transmission scenarios. When we combined evidence from over 30 genotypes with conventional epidemiological data, we were able to identify the sources of several unexplained transmission events and target additional investigations and infection-control interventions to reduce the risk of further transmission.
In this study, we were able to link patient-derived isolates directly to environmental isolates. Early in the outbreak, confirmation of contamination of the patients' surroundings on the ward led to a tightening of ward disinfection procedures. Later, WGS showed that environmental isolates from a burns operating theatre were linked to patient isolates, which triggered closure of the theatre, followed by deep cleaning. Previous reports have highlighted that burns patients are particularly at risk of infection with Acinetobacter baumannii and that burns units are susceptible to outbreaks and have also implicated contaminated mattresses in infection [30],[31].
We have demonstrated that there is considerable genotypic diversity within patient-associated and environmental populations of Acinetobacter. We discovered two examples of what were clearly mixed infections, where more than one species or strain of Acinetobacter was present in serial samples from the same site. In addition, when we genome-sequenced multiple isolates from multiple samples from a chronically infected patient, we found that the individual was carrying a cloud of variants of the outbreak strain, rather than a single SNV genotype. It remains unclear whether such genotypic diversity evolves in situ or whether multiple genotypes are acquired serially or en masse. Although some have suggested that such genotypic diversity precludes the reconstruction of transmission events from genomic data [32], the congruence and consistency we and others [12] see between genotypic and epidemiologic data in such analyses confirm the utility of such approaches. However, the presence of such genotypic diversity in pathogen populations challenges the long held practice of picking and propagating a single example of each colonial morphotype from a primary culture and suggests a role for metagenomic approaches that might capture bacterial genetic heterogeneity directly from the clinical sample without culture [33]-[35].
Although, in this study, WGS was initially largely retrospective, in the latter part of the outbreak, we were able to go from colony to SNV genotype in less than 1 week. Local WGS thus allowed us to rule patients and isolates in or out of the outbreak more quickly than through the reference laboratory. Excluding isolates meant that we could focus our efforts on determining the connections between genetically related cases, rather than trying to connect all cases of MDR-Aci.
Conclusions
WGS is now poised to make an impact on hospital infection prevention and control, delivering cost-effective identification of routes of infection within a clinically relevant timeframe and allowing infection control teams to track, and even prevent, the spread of drug-resistant hospital pathogens.
Authors' contributions
MRH performed bioinformatics analyses of genomes and wrote a script for assessing the parsimonious transmission routes; JZ-MC performed whole genome sequencing and bioinformatics analysis and helped write the paper; CIC helped design the study and performed whole genome sequencing; NC, CB and M-SB supplied isolates and patient metadata; BO oversaw clinical investigation of the outbreak, helped design the study and write the manuscript; MJP conceived, designed and supervised the study and wrote the manuscript. All authors read and approved the final manuscript.
Abbreviations
- MDR-Aci:
-
Multi-drug-resistant Acinetobacter
- MIC:
-
Minimum inhibitory concentration
- PFGE:
-
Pulsed field gel electrophoresis
- QEHB:
-
Queen Elizabeth Hospital Birmingham
- SNV:
-
Single nucleotide variant
- WGS:
-
Whole-genome sequencing
References
Munoz-Price LS, Weinstein RA: Acinetobacter infection. N Engl J Med. 2008, 358: 1271-1281. 10.1056/NEJMra070741.
Towner KJ: Acinetobacter: an old friend, but a new enemy. J Hosp Infect. 2009, 73: 355-363. 10.1016/j.jhin.2009.03.032.
Roca I, Espinal P, Vila-Farres X, Vila J: The Acinetobacter baumannii oxymoron: commensal hospital dweller turned pan-drug-resistant menace. Front Microbiol. 2012, 3: 148-10.3389/fmicb.2012.00148.
O’Shea MK: Acinetobacter in modern warfare. Int J Antimicrob Agents. 2012, 39: 363-375. 10.1016/j.ijantimicag.2012.01.018.
Whitman TJ, Qasba SS, Timpone JG, Babel BS, Kasper MR, English JF, Sanders JW, Hujer KM, Hujer AM, Endimiani A, Eshoo MW, Bonomo RA: Occupational transmission of Acinetobacter baumannii from a United States serviceman wounded in Iraq to a health care worker. Clin Infect Dis. 2008, 47: 439-443. 10.1086/589247.
Hujer KM, Hujer AM, Hulten EA, Bajaksouzian S, Adams JM, Donskey CJ, Ecker DJ, Massire C, Eshoo MW, Sampath R, Thomson JM, Rather PN, Craft DW, Fishbain JT, Ewell AJ, Jacobs MR, Paterson DL, Bonomo RA: Analysis of antibiotic resistance genes in multidrug-resistant Acinetobacter sp. isolates from military and civilian patients treated at the Walter Reed Army Medical Center. Antimicrob Agents Chemother. 2006, 50: 4114-4123. 10.1128/AAC.00778-06.
Lewis T, Loman NJ, Bingle L, Jumaa P, Weinstock GM, Mortiboy D, Pallen MJ: High-throughput whole-genome sequencing to dissect the epidemiology of Acinetobacter baumannii isolates from a hospital outbreak. J Hosp Infect. 2010, 75: 37-41. 10.1016/j.jhin.2010.01.012.
Jones A, Morgan D, Walsh A, Turton J, Livermore D, Pitt T, Green A, Gill M, Mortiboy D: Importation of multidrug-resistant Acinetobacter spp infections with casualties from Iraq. Lancet Infect Dis. 2006, 6: 317-318. 10.1016/S1473-3099(06)70471-6.
Turton JF, Kaufmann ME, Gill MJ, Pike R, Scott PT, Fishbain J, Craft D, Deye G, Riddell S, Lindler LE, Pitt TL: Comparison of Acinetobacter baumannii isolates from the United Kingdom and the United States that were associated with repatriated casualties of the Iraq conflict. J Clin Microbiol. 2006, 44: 2630-2634. 10.1128/JCM.00547-06.
Turton JF, Gabriel SN, Valderrey C, Kaufmann ME, Pitt TL: Use of sequence-based typing and multiplex PCR to identify clonal lineages of outbreak strains of Acinetobacter baumannii. Clin Microbiol Infect. 2007, 13: 807-815. 10.1111/j.1469-0691.2007.01759.x.
Turton JF, Baddal B, Perry C: Use of the accessory genome for characterization and typing of Acinetobacter baumannii. J Clin Microbiol. 2011, 49: 1260-1266. 10.1128/JCM.02335-10.
Robinson ER, Walker TM, Pallen MJ: Genomics and outbreak investigation: from sequence to consequence. Genome Med. 2013, 5: 36-
Reuter S, Ellington MJ, Cartwright EJ, Koser CU, Torok ME, Gouliouris T, Harris SR, Brown NM, Holden MT, Quail M, Parkhill J, Smith GP, Bentley SD, Peacock SJ: Rapid bacterial whole-genome sequencing to enhance diagnostic and public health microbiology. JAMA Intern Med. 2013, 173: 1397-1404. 10.1001/jamainternmed.2013.7734.
Gardy JL, Johnston JC, Ho Sui SJ, Cook VJ, Shah L, Brodkin E, Rempel S, Moore R, Zhao Y, Holt R, Varhol R, Birol I, Lem M, Sharma MK, Elwood K, Jones SJ, Brinkman FS, Brunham RC, Tang P: Whole-genome sequencing and social-network analysis of a tuberculosis outbreak. N Engl J Med. 2011, 364: 730-739. 10.1056/NEJMoa1003176.
Snitkin ES, Zelazny AM, Thomas PJ, Stock F, Henderson DK, Palmore TN, Segre JA: Tracking a hospital outbreak of carbapenem-resistant Klebsiella pneumoniae with whole-genome sequencing. Sci Transl Med. 2012, 4: 148ra116-10.1126/scitranslmed.3004129.
Eyre DW, Cule ML, Wilson DJ, Griffiths D, Vaughan A, O'Connor L, Ip CL, Golubchik T, Batty EM, Finney JM, Wyllie DH, Didelot X, Piazza P, Bowden R, Dingle KE, Harding RM, Crook DW, Wilcox MH, Peto TE, Walker AS: Diverse sources of C. difficile infection identified on whole-genome sequencing. N Engl J Med. 2013, 369: 1195-1205. 10.1056/NEJMoa1216064.
Rohde H, Qin J, Cui Y, Li D, Loman NJ, Hentschke M, Chen W, Pu F, Peng Y, Li J, Xi F, Li S, Li Y, Zhang Z, Yang X, Zhao M, Wang P, Guan Y, Cen Z, Zhao X, Christner M, Kobbe R, Loos S, Oh J, Yang L, Danchin A, Gao GF, Song Y, Li Y, Yang H, Wang J, Xu J, Pallen MJ, Wang J, Aepfelbacher M, Yang R: Open-source genomic analysis of Shiga-toxin-producing E. coli O104:H4. N Engl J Med. 2011, 365: 718-724. 10.1056/NEJMoa1107643.
Joyanes P, Del Carmen CM, Martinez-Martinez L, Perea EJ: Evaluation of the VITEK 2 system for the identification and susceptibility testing of three species of nonfermenting gram-negative rods frequently isolated from clinical samples. J Clin Microbiol. 2001, 39: 3247-3253. 10.1128/JCM.39.9.3247-3253.2001.
Miller JR, Koren S, Sutton G: Assembly algorithms for next-generation sequencing data. Genomics. 2010, 95: 315-327. 10.1016/j.ygeno.2010.03.001.
Goris J, Konstantinidis KT, Klappenbach JA, Coenye T, Vandamme P, Tiedje JM: DNA-DNA hybridization values and their relationship to whole-genome sequence similarities. Int J Syst Evol Microbiol. 2007, 57: 81-91. 10.1099/ijs.0.64483-0.
Chan JZ, Halachev MR, Loman NJ, Constantinidou C, Pallen MJ: Defining bacterial species in the genomic era: insights from the genus Acinetobacter. BMC Microbiol. 2012, 12: 302-10.1186/1471-2180-12-302.
Langmead B, Salzberg SL: Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012, 9: 357-359. 10.1038/nmeth.1923.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R: The sequence alignment/map format and SAMtools. Bioinformatics. 2009, 25: 2078-2079. 10.1093/bioinformatics/btp352.
Jawad A, Seifert H, Snelling AM, Heritage J, Hawkey PM: Survival of Acinetobacter baumannii on dry surfaces: comparison of outbreak and sporadic isolates. J Clin Microbiol. 1998, 36: 1938-1941.
Beceiro A, Moreno A, Fernandez N, Vallejo JA, Aranda J, Adler B, Harper M, Boyce JD, Bou G: Biological cost and impact on virulence of different mechanisms of colistin resistance in Acinetobacter baumannii. Antimicrob Agents Chemother. 2014, 58: 518-526. 10.1128/AAC.01597-13.
Koser CU, Holden MT, Ellington MJ, Cartwright EJ, Brown NM, Ogilvy-Stuart AL, Hsu LY, Chewapreecha C, Croucher NJ, Harris SR, Sanders M, Enright MC, Dougan G, Bentley SD, Parkhill J, Fraser LJ, Betley JR, Schulz-Trieglaff OB, Smith GP, Peacock SJ: Rapid whole-genome sequencing for investigation of a neonatal MRSA outbreak. N Engl J Med. 2012, 366: 2267-2275. 10.1056/NEJMoa1109910.
Young BC, Golubchik T, Batty EM, Fung R, Larner-Svensson H, Votintseva AA, Miller RR, Godwin H, Knox K, Everitt RG, Iqbal Z, Rimmer AJ, Cule M, Ip CL, Didelot X, Harding RM, Donnelly P, Peto TE, Crook DW, Bowden R, Wilson DJ: Evolutionary dynamics of Staphylococcus aureus during progression from carriage to disease. Proc Natl Acad Sci U S A. 2012, 109: 4550-4555. 10.1073/pnas.1113219109.
Snitkin ES, Zelazny AM, Montero CI, Stock F, Mijares L, Murray PR, Segre JA: Genome-wide recombination drives diversification of epidemic strains of Acinetobacter baumannii. Proc Natl Acad Sci U S A. 2011, 108: 13758-13763. 10.1073/pnas.1104404108.
Wright MS, Haft DH, Harkins DM, Perez F, Hujer KM, Bajaksouzian S, Benard MF, Jacobs MR, Bonomo RA, Adams MD: New insights into dissemination and variation of the health care-associated pathogen Acinetobacter baumannii from genomic analysis. MBio. 2014, 5: e00963-13. 10.1128/mBio.00963-13.
Sherertz RJ, Sullivan ML: An outbreak of infections with Acinetobacter calcoaceticus in burn patients: contamination of patients' mattresses. J Infect Dis. 1985, 151: 252-258. 10.1093/infdis/151.2.252.
Babik J, Bodnarova L, Sopko K: Acinetobacter- serious danger for burn patients. Acta Chir Plast. 2008, 50: 27-32.
Worby CJ, Lipsitch M, Hanage WP: Within-host bacterial diversity hinders accurate reconstruction of transmission networks from genomic distance data. PLoS Comput Biol. 2014, 10: e1003549-10.1371/journal.pcbi.1003549.
Loman NJ, Constantinidou C, Christner M, Rohde H, Chan JZ, Quick J, Weir JC, Quince C, Smith GP, Betley JR, Aepfelbacher M, Pallen MJ: A culture-independent sequence-based metagenomics approach to the investigation of an outbreak of Shiga-toxigenic Escherichia coli O104:H4. JAMA. 2013, 309: 1502-1510. 10.1001/jama.2013.3231.
Pallen MJ: Diagnostic metagenomics: potential applications to bacterial, viral and parasitic infections.Parasitology 2014, [Epub ahead of print].,
Schloissnig S, Arumugam M, Sunagawa S, Mitreva M, Tap J, Zhu A, Waller A, Mende DR, Kultima JR, Martin J, Kota K, Sunyaev SR, Weinstock GM, Bork P: Genomic variation landscape of the human gut microbiome. Nature. 2013, 493: 45-50. 10.1038/nature11711.
Acknowledgements
We thank Jane Turton and her team at Public Health England for providing strains and access to routinely collected typing information; we thank Nick Loman, Danny Wilson and Mark Achtman for advice in drafting the manuscript. This paper presents independent research funded by the UK Medical Research Council (MRC) and the National Institute for Health Research (NIHR) Surgical Reconstruction and Microbiology Research Centre (partnership between University Hospitals Birmingham NHS Foundation Trust, the University of Birmingham and the Royal Centre for Defence Medicine). The views expressed are those of the author(s) and not necessarily those of the MRC, NHS, the NIHR or the Department of Health.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
About this article
Cite this article
Halachev, M.R., Chan, J.ZM., Constantinidou, C.I. et al. Genomic epidemiology of a protracted hospital outbreak caused by multidrug-resistant Acinetobacter baumanniiin Birmingham, England. Genome Med 6, 70 (2014). https://doi.org/10.1186/s13073-014-0070-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13073-014-0070-x