Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
Sci Transl Med. Author manuscript; available in PMC 2013 August 22.
Published in final edited form as:
PMCID: PMC3521604

Tracking a Hospital Outbreak of Carbapenem-Resistant Klebsiella pneumoniae with Whole-Genome Sequencing


The Gram-negative bacteria Klebsiella pneumoniae is a major cause of nosocomial infections, primarily among immunocompromised patients. The emergence of strains resistant to carbapenems has left few treatment options, making infection containment critical. In 2011, the U.S. National Institutes of Health Clinical Center experienced an outbreak of carbapenem-resistant K. pneumoniae that affected 18 patients, 11 of whom died. Whole-genome sequencing was performed on K. pneumoniae isolates to gain insight into why the outbreak progressed despite early implementation of infection control procedures. Integrated genomic and epidemiological analysis traced the outbreak to three independent transmissions from a single patient who was discharged 3 weeks before the next case became clinically apparent. Additional genomic comparisons provided evidence for unexpected transmission routes, with subsequent mining of epidemiological data pointing to possible explanations for these transmissions. Our analysis demonstrates that integration of genomic and epidemiological data can yield actionable insights and facilitate the control of nosocomial transmission.


The bacterial pathogen Klebsiella pneumoniae is responsible for roughly 15% of Gram-negative infections in hospital intensive care units (ICUs) (1), primarily affecting immunocompromised patients (2). In recent years, the threat posed by K. pneumoniae has markedly increased with the emergence of strains resistant to carbapenem antibiotics (3)and their worldwide dissemination (4, 5). Infections caused by carbapenem-resistant strains have few treatment options (6, 7)and areassociated with mortality rates upwards of 50% (8, 9). Although multiple resistance mechanisms have been identified (10), carbapenem resistance in the United States is primarily caused by the plasmid-encoded K. pneumoniae carbapenemase (KPC) gene (5).

Exacerbating the problems associated with the emergence of these highly resistant strains of K. pneumoniae is their propensity to cause outbreaks in health care institutions (1113). As with most nosocomial pathogens, multiple-drug resistance offers inherent selective advantage (11), which allows such organisms to persist both in the flora of hospitalized patients and in the hospital environment, in which antibiotic usage is widespread. More specific to K. pneumoniae is its capacity to silently colonize patients or hospital personnel (14), that is, by establishing residence in the gastrointestinal tract without causing any signs of infection. Individuals can be silently colonized or asymptomatic carriers for long periods of time, with detection of these carriers often proving difficult (12). These silent carriers act as reservoirs for continued transmission that make spread difficult to control and outbreaks difficult to stop (13). In addition, K. pneumoniae can survive for several hours on the hands of hospital personnel, which likely facilitates nosocomial spread (15).

Effective control of K. pneumoniae outbreaks requires a detailed understanding of how transmission occurs. Molecular typing approaches, such as pulsed-field gel electrophoresis and multi-locus sequence typing, have been used to classify K. pneumoniae andthusunderstandits local and global dissemination (16, 17). However, the high K. pneumoniae clonality (16, 18) creates difficulty in tracking outbreaks because available methods may not provide sufficient resolution to distinguish intrainstitutional spread from introduction of closely related strains from other health care facilities. In the United States, KPC-K. pneumoniae isolates are highly clonal, with 70% belonging to sequence type (ST) 258 (16).

With rapid technological advances, whole-genome sequencing is emerging as the gold standard in bacterial typing (19, 20). Success in tracking worldwide epidemics (2124), regional outbreaks (23, 25), food-borne outbreaks (26, 27), and bioterrorism agents (28) has demonstrated that the fine resolution provided by whole-genome sequencing facilitates our understanding of the spread of infectious agents. The continued improvements in turnaround time and accessibility of DNA sequencing technologies are now approaching a point where genomic data can be generated in a clinically relevant time frame. Genome sequencing has been applied recently to nosocomial strains (2932)but with limited study size or scope in reconstructing transmission links during the course of the outbreak.

Here, we applied whole-genome sequencing to track an outbreak of carbapenem-resistant K. pneumoniae at the U.S. National Institutes of Health (NIH) Clinical Center that colonized 18 patients, with 6 deaths attributable to K. pneumoniae infection. We developed an algorithm to reconstruct the outbreak transmissions based on whole-genome sequencing of isolates and epidemiological data that tracked the location of patients throughout their hospital stays.


Overview of outbreak

On 13 June 2011, patient 1 was transferred to our ICU from a hospital in New York City and was discharged on 15 July. She was known to be colonized with carbapenem-resistant K. pneumoniae and was immediately placed in enhanced contact isolation in a private room in which staff and visitors are required to don gloves and gowns for entry. During her stay, she spent two 24-hour stints in the ICU. Another instance of a KPC-K. pneumoniae isolate was not observed in clinical or surveil-lance cultures until August 5 when it was cultured from the tracheal aspirate of patient 2. Whether the two KPC-K. pneumoniae infections were related was unclear, especially because the two patients were never housed in the same ward at the same time (Fig. 1, A and B). Repetitive polymerase chain reaction (PCR) and pulsed-field gel electrophoresis (fig. S1) identified both patient isolates as belonging to the epidemic ST 258 K. pneumoniae lineage. However, because of the ubiquity of ST 258 in U.S. hospitals (16) and the time lag between patient presentations, both nosocomial transmission and independent introductions were deemed possible.

Fig. 1
Patient location and overlap during the outbreak. (A) Timeline of first positive cultures of the outbreak strain for the 18 affected patients. (B) Patient traces for each of the 18 patients shown in (A). Black lines, first positive culture; blue lines, ...

In the following weeks, an average of 1 new case of colonization or active infection with KPC-K. pneumoniae was detected each week in the NIH Clinical Center, with a total of 17 cases as of 1 January 2012 (table S2). During the outbreak, 10 of the 17 KPC-K. pneumoniae- colonized patients died, with 6 deaths attributable to K. pneumoniae infection and 4 to the underlying disease (Table 1). The outbreak was ultimately contained by implementing strict cohorting of colonized patients to minimize sharing of hospital equipment and of care providers between outbreak patients and the other patients in the hospital (see Materials and Methods for details). To identify asymptomatically colonized patients, three rounds of rectal surveillance cultures were performed on all hospital patients, augmented with more frequent rectal surveillance cultures performed on patients in selected hospital wards (see Materials and Methods). Although implementation of rigorous infection control procedures ultimately halted the outbreak, we sought to understand better how the outbreak progressed to learn how to control future outbreaks of K. pneumonia more effectively.

Table 1
Clinical histories of patients colonized with the outbreak of KPC-producing K. pneumonia strain.

Genome sequencing and comparison of outbreak isolates

Deciphering the transmission events between patients solely on the basis of the epidemiologic data proved extremely difficult because extensive overlap of patients within the hospital wards, particularly the ICU, supported numerous outbreak scenarios (Fig. 1C). To learn more about carbapenem-resistant K. pneumoniae transmission, first we performed whole-genome sequencing of an isolate from the index patient. An additional six isolates were cultured from four body sites over the index patient's 4-week stay at the NIH Clinical Center. Because the index patient had been colonized for several months, we examined potential genetic heterogeneity among these seven isolates to reconstruct putative transmission events more accurately. Sequencing of variants of the index patient's isolates revealed a total of seven single nucleotide variants (SNVs) among the seven isolates (Fig. 2A). The four urine isolates matched the NTUH-K2044 K. pneumoniae reference strain at each of the seven variant sites, supporting this as the “ancestral” genotype. By contrast, the groin and bronchoalveolar lavage (BAL) isolates both exhibited three SNVs compared to the urine isolates, and the throat isolate exhibited three different SNVs compared to the urine, distinct from the groin and BAL isolates. Note that the groin isolate also contained one SNV in addition to the BAL isolate. Identification of this genetic heterogeneity within the isolates recovered from the index patient proved immensely valuable in reconstructing transmission and demonstrated the importance of culturing and sequencing isolates from multiple body sites over time from long-term colonized patients.

Fig. 2
SNVs identified in K. pneumoniae genomes. (A) DNA sequence variation among isolates taken from patient 1 while at the NIH Clinical Center is shown as a heat map, with isolates in chronological order on the y axis and variants on the x axis. Gray, ancestral ...

We next performed whole-genome sequencing of a single KPC-K. pneumoniae isolate from each of the 18 affected patients, in hopes of using the SNVs found in their genomes to infer transmission links. Genomic comparisons revealed a total of 41 SNVs among the 18 outbreak patient isolates, that is, the isolates varied at a total of 41 sites distributed over their 6 million base pair (bp) genomes (table S4).

To address the central question of whether patient 1 initiated the outbreak and, if so, how she was linked to other patients, we compared the 18 isolates' genomes. Grouping KPC–K. pneumoniae outbreak patient isolates on the basis of the patterns of shared variants partitioned them into two large clusters and a third cluster consisting only of patient 8 (Fig. 2B). Inspection of the variants common among these clusters revealed that members of clusters 1 and 3 share the three variants present in patient 1's BAL and groin isolates, whereas the isolates in cluster 2 share the three variants present in patient 1's throat isolate (Fig. 2B). These data therefore suggest that genetically distinct isolates from patient 1 were transmitted to other patients. Moreover, the data indicate that not only is patient 1 linked to the outbreak but also three independent transmission events from this patient led to hospital-wide dissemination of the outbreak strain.

Inference of most likely transmission route

Next, we sought to extract insights from the genetic data more rigorously and to reconstruct a putative transmission map. To distinguish among scenarios that were equally likely based on the genetic data, we used patient overlap in hospital wards and deemed the most probable transmission sequence as the one requiring the minimal length of undetected silent colonization (see Materials and Methods for details of algorithm implementation). The transmission map generated by our approach (Fig. 3) unravels the complexity of Fig. 1C and represents an outbreak progression that is most parsimonious with respect to both genetic and patient overlap data. Comparisons of this map to transmission maps generated on the basis of either genetic (fig. S4) or patient trace data (fig. S5) alone revealed a marked effect of incorporating both types of data.

Fig. 3
Putative map of K. pneumoniae transmission during outbreak. The transmission map was constructed with genetic and patient trace data, as detailed in Materials and Methods. Nodes represent patients, and arrows indicate a transmission event directly or ...

Consistent with our observations above, the integrated map identified three transmissions stemming from the index patient. However, the first transmission was predicted to go through patient 3, and not patient 2, despite patient 2's presentation with infection 10 days earlier than patient 3's. This scenario was supported by the genetic information (patient 2's isolate contains an SNV not found in patient 1'sor 3's) and also by the patient trace, which indicates that patient 3, but not patient 2, overlapped with patient 1 in the ICU (Fig. 1B).

A second transmission from patient 1 was predicted to go through patient 4 and is based entirely on genetic data. Because no direct contact occurred between patients 1 and 4, we looked for evidence of patients acting as asymptomatic intermediates in a transmission between these two patients. Knowledge of patient locations in the hospital allowed us to identify patients whose movements within the hospital make them candidates for being silent transmission vectors. Specifically, we identified all patients who overlapped with the index patient and then with patient 4 before he cultured positive. Among the 1115 patients at the NIH Clinical Center during the outbreak, there were only 5 patients who overlapped with patients 1 and 4 and could have acted as vectors for transmission (Fig. 4). Patients B and D are especially compelling because of their extensive overlap with patients 1 and 4, but neither one cultured positive with surveillance cultures. Although the asymptomatic carrier could have been colonized below detection level or could have been an untested health care provider, this type of mining of epidemiological data has the potential to identify a handful of candidates who merit additional surveillance cultures and/or placement in contact isolation. As genome sequencing becomes even faster, such insights could be obtained in real time to perform targeted, thorough surveillance of patients of interest.

Fig. 4
Patient trace for patients 1 and 4 and all patients who could have acted as transmission intermediates between them. Patient traces for the 1115 patients present in the hospital during the outbreak were mined to identify those patients whose location ...

Four other links in the patient transmission chart cannot be explained by geographic overlap in the same ward (Fig. 1C). These also suggest more complicated transmission modes, perhaps through asymptomatic patients who were never detected or via health care personnel or equipment. Our epidemiologic investigation revealed that no single staff member, group of health care personnel, or procedure was a likely source of transmission to the nonoverlapping patients. However, the possibility of transmission through inanimate objects could not be ruledout and was in fact bolstered by culturing the out-break strain from six sink drains and from a ventilator. The ventilator had been thoroughly cleaned after being used on patient 6, twice with a quaternary ammonium compound and once with bleach. The ventilator isolate had only one SNV distinguishing it from patient 6's original isolate, confirming that the isolate likely survived the cleaning process. All other environmental cultures were negative, with the exception of those taken from the room of patient 8 before decontamination.

Identification of mutations associated with colistin resistance

Finally, we compared genome sequences of outbreak isolates to gain insight into evolution of antibiotic resistance during the course of the outbreak. The index patient's isolate was highly resistant to multiple antibiotics, with minimum inhibitory concentration (MIC) values in the susceptible range for only gentamicin, tigecycline, and colistin (table S3). During the course of the outbreak, isolates with resistant MICs for all of these drugs were observed, leaving no effective therapeutic options for some patients. We focused on colistin resistance because it appeared to evolve independently in all three (only two are described) patient clusters, indicating that the outbreak strain was on the verge of becoming colistin-resistant upon arrival at the NIH Clinical Center. Patient 8's initial isolate was colistin-resistant, and patient 2's isolate acquired colistin resistance during the course of treatment. The resistant isolate from patient 2 had two mutations distinguishing it from its susceptible ancestor, and the isolate from patient 8 had six unique mutations relative to that of patient 1's (Table 2). Patient 2's colistin-resistant isolate had a 2-bp insertion in the coding region of a putative microcin transporter, resulting in a truncated protein. In Salmonella enterica (33) and Escherichia coli (34), mutations disabling the microcin transporter are sufficient to confer resistance to antimicrobial peptides in vitro, although not to colistin. Of the six mutations unique to patient 8's isolate, four result in amino acid changes, three of which are in putative membrane proteins. Thus, in both cases, mutations altering membrane function, possibly reducing drug influx, are associated with colistin resistance.

Table 2
Mutations unique to colistin-resistant isolates.


The increasing speed and decreasing cost of high-throughput DNA sequencing technologies are enabling its application to the practice of medicine (35). Here, we tested whether genome sequencing could help to unravel a nosocomial outbreak and affect hospital infection control decisions. We sequenced patient and environmental isolates within a clinically relevant turnaround time during a hospital KPC–K. pneumoniae outbreak. Among the key insights provided by sequencing were that (i) the outbreak was monoclonal, despite a 3-week interval between the index case and the identification of subsequent cases, (ii) transmission likely occurred from at least two different sites of the index patient, (iii) at least three independent transmission events from the index patient led to two major clusters of colonized patients, (iv) one patient could be linked to a contaminated ventilator, and (v) a small number of putative resistance mutations could be identified in newly colistin-resistant isolates. Moreover, by combining the genetic data with patient trace data and by implementing an algorithm that accounted for K. pneumoniae's capacity for silent colonization, we reconstructed the most likely outbreak transmission route.

Our results point to ways in which carbapenem-resistant K. pneumoniae outbreaks may be controlled more effectively in the future. First, although additional sequencing studies capturing the genetic diversity of clinical strains of K. pneumoniae will be necessary to verify that this outbreak did not stem from multiple introductions, the genetic diversity that accumulated just during the course of the outbreak suggests that a monoclonal outbreak is most likely. The finding that the outbreak likely stemmed from an index patient who left the hospital 3 weeks before the strain was isolated from another patient provides direct evidence that K. pneumoniae's ability to silently colonize patients allows an outbreak to develop stealthily over an extended period of time (14). The potential for undetected transmission provides support for extensive and frequent surveillance during and after discharge of a carbapenem-resistant K. pneumoniae–colonized patient. Second, we think it is likely that groin and throat surveillance cultures performed on patients residing in the ICU concurrent with the index patient failed to detect that one or more of these patients had indeed become colonized. Whether these early cultures were ineffective because of the site of culturing (groin and throat instead of rectum) or because of less sensitive culturing techniques used early in the outbreak (MacConkey agar instead of CHROMagar) is unclear. Nevertheless, our results demonstrate the importance of having ongoing, effective surveillance protocols in place before outbreaks occur. Third, the finding that K. pneumoniae could survive on hospital equipment that had undergone rigorous cleaning suggests that K. pneumoniae has a high degree of environmental stability under some circumstances (13) and indicates that decontamination should be verified after the cleaning process.

In addition to supporting the implementation of these specific infection control procedures, our results suggest several ways in which whole-genome sequencing of outbreak isolates in real time could guide future infection control efforts. First, genetic data can allow for the identification of unexpected modes of transmission. For instance, it was initially assumed that patient 4 was colonized during a 24-hour stay in the ICU, during which he overlapped with patients 2 and 3. However, sequence analysis demonstrated that patient 4's isolate could not have come from patients 2 or 3 and must have derived from an independent transmission chain from patient 1. This finding motivated the search for an intermediate patient to explain the transmission from patient 1 to patient 4, ultimately revealing four highly plausible candidates. If we had such knowledge in real time, these putative intermediate patients could have undergone more rigorous surveillance culturing to identify KPC–K. pneumoniae colonization and/or been placed on contact isolation, potentially terminating a silent transmission chain. Second, genomic sequencing may distinguish between alternate transmission scenarios, which may be critical for surmising the scope of an outbreak within the hospital. For example, patients 15 and 16 cultured positive for KPC-producing K. pneumoniae within 1 day of each other while residing in the same non-ICU ward. There was no obvious epidemiological link connecting these patients to the outbreak, which raised the possibility that patient 16 had brought a new strain to the hospital from a long-term care facility in which he had recently had a prolonged stay. Genomic sequencing revealed that patient 15's and 16's isolates both matched the dominant ICU strain, suggesting that transmission had occurred within our hospital to a non-ICU ward. Finally, genetic data can link patients directly to environmental or infrastructure isolates. Such findings can motivate refinement in cleaning and decontamination procedures by providing insight as to how and when contamination occurred.

Beyond applications to outbreak containment, we foresee that future applications of real-time sequencing have the potential to transform both the control and the treatment of hospital infection. In both epidemic and endemic settings, real-time genomic sequencing can ultimately provide a powerful tool to define the nosocomial epidemiology of important health care–associated pathogens with heretofore unprecedented precision. Sequencing serial isolates in patients undergoing antibiotic therapy can allow unprecedented insight into how bacterial populations respond to treatments and the evolution of resistance (36, 37). Furthermore, by observing how bacteria evolve in different patient populations and in response to different treatments, optimal therapeutic regimes may be tailored, eliminating those associated with the development of resistance. Finally, genetic variation underlying changes in bacterial virulence or tissue specificity may be identified. With the identification of markers of these key phenotypic attributes may come the ability for infection control to anticipate and outpace the evolution of infectious agents.


Description of outbreak

In June 2011, a patient, known to be colonized with KPC–K. pneumoniae (hereafter called the index patient), was admitted to the NIH Clinical Center, a 243-bed clinical research hospital. After this patient's two overnight stays in the ICU, the hospital experienced a largely ICU-based outbreak of carbapenem-resistant K. pneumoniae despite early implementation of extensive infection control measures described below.

Active surveillance cultures

Immediately after index patient's second ICU stay, we collected two sets of throat and groin cultures (BBL CultureSwab with Stuart's Transport Medium, Becton-Dickinson) on all contemporaneous ICU patients. After detection of the second KPC–K. pneumoniae case, surveillance groin and throat cultures were performed on all ICU patients at admission and twice weekly thereafter. We believe that at least in the instance of patient 3, these cultures failed to identify a colonization. Surveillance culture–negative patients transferred to other units remained in enhanced contact isolation (described below) until they had additional negative surveillance cultures. Although rectal cultures are more sensitive for KPC–K. pneumoniae detection, our clinical microbiology laboratory was only able to validate this methodology 4 weeks into the outbreak when appropriate selective medium was acquired. At this time, we implemented rectal surveillance cultures in addition to groin and throat cultures for patients in the ICU. Rectal screening was performed in high-risk wards as cases were detected and three times in the entire hospital patient population between August 2011 and January 2012.

Infection control interventions

The index patient was placed in enhanced contact isolation upon arrival at our hospital and throughout her stay. During the outbreak, all ICU patients, regardless of colonization status, were placed on enhanced contact isolation, including strict enforcement of fastidious hand hygiene before entering and upon leaving patient rooms, universal use of gowns and gloves for all staff and visitors entering patient rooms, restrictions on patient activity outside rooms, restrictions on staff and visitor traffic, dedicating equipment for single-patient use (when feasible), extensive cleaning of shared equipment, and double-cleaning of vacated rooms with bleach. A private company was hired to decontaminate the ICU and all rooms and equipment used by KPC–K. pneumoniae patients with hydrogen peroxide vapor (38).

ICU patients positive for KPC–K. pneumoniae were transferred to a cohorted ICU after we recognized the fourth case in the outbreak. Patient cohorting included geographic separation of KPC–K. pneumoniae–colonized patients from other hospital patients. Staff cohorting included the use of a dedicated cadre of staff, including nursing, respiratory therapy, housekeeping, and physical therapy staff, to care only for KPC–K. pneumoniae–colonized patients. This cohorted staff did not provide care to noncohorted patients during a given shift. Health care personnel were assigned around the clock to monitor and enforce adherence to infection control precautions (39). A second non-ICU cohorted area, as described above, was established for patients who were not critically ill. Physicians were not cohorted and were allowed to care for both KPC–K. pneumoniae–colonized and other patients.

Outbreak investigation

Data were extracted from clinical charts detailing patients' clinical history, contact with staff members or ancillary departments, and movements throughout the hospital. Patient traces were conducted with TheraDoc (Hospira) to establish temporal and geographic overlap among patients. Throughout the outbreak, 164 environmental surveillance cultures were collected with extensive sampling on high-touch surfaces in patient rooms, the ICU, inpatient units with detected transmission, and ancillary departments of the hospital. Sink drain cultures were collected with swabs inserted through the sieve.

Detection of KPC-positive K. pneumoniae

Early screening cultures for K. pneumoniae were performed with MacConkey agar (Remel), and then starting in September 2011 with more sensitive CHROMagar KPC. Species level identification was done with matrix-assisted laser desorption/ionization–time-of flight (MALDITOF) microflex LT mass spectrometer (Bruker Daltonics Inc.), and carbapenemase production was assessed by the modified Hodge test. Antibiotic susceptibility testing was performed with the BD Phoenix (Becton-Dickinson) and Etest (bioMérieux) systems.

Genome sequencing, assembly, and annotation

Genomic libraries were constructed with the Roche 454 Titanium Kit (Roche 454 Life Sciences). Sequencing was performed on the Roche/454 XLR instrument with an average N50 contig size of 154,336 bp and an average of 135 large contigs. Sequencing assembly statistics for each isolate are provided in table S1. Contig assembly was executed with the gsAssembler v.2.3. Genome annotation was performed with the publicly available Prokaryotic Genomes Automatic Annotation Pipeline.

Detection and filtering of nucleotide differences

Large contigs from each genome assembly were ordered and oriented relative to a finished reference K. pneumoniae genome (NTUH-K2044) and then stitched together to form a pseudo-chromosome with the Mauve contig mover (40). An initial list of SNVs among the outbreak isolates was generated from the Mauve alignments with the export single-nucleotide polymorphism (SNP) function (41). SNVs were filtered to retain only high-confidence SNVs as detailed in the Supplementary Materials.

Constructing the putative transmission map

To construct the most probable transmission map, we integrated both genetic and epidemiological information, building on an approach reported in (42). The most parsimonious transmission map was generated by finding the set of links, spanning all patients with the minimum total genetic distance, with patient trace information being used to differentiate between transmission scenarios that were equally probable based on the genetic data. We incorporated epidemiological data quantitatively by taking into consideration how K. pneumoniae nosocomial outbreaks manifest. Specifically, we accounted for these aspects: (i) outbreaks are thought to spread primarily through patient-to-patient transmission via hospital personnel, equipment, and infrastructure (13); (ii) silent colonization of patients may result in detection of colonization in a recipient well after the actual transmission event (14); and (iii) silent colonization of patients may result in instances in which a source patient cultures positive after another patient to whom s/he transmitted. To account for patient-to-patient transmission as the most likely route, we penalized links with patient A transmitting to patient B if they did not overlap in the same ward before patient B grew the outbreak organism. While allowing for silent colonization, we down-weighted links based on the total number of days of silent colonization required (fig. S6). Despite the reported sensitivity of rectal surveillance (43), for the purpose of constructing the transmission map, we chose to not preclude patients from being carriers based on negative rectal surveillance. This decision was made because the true sensitivity of rectal surveillance remains unknown. Although previous studies show that rectal surveillance is the most sensitive assay, there is no way to accurately quantify false negatives. Evidence for sensitivity not being 100% comes from instances in which rectal surveillance cultures within single patients fluctuate between positive and negative (for example, patient 9 in fig. S2).

Despite the possible limitations of rectal surveillance, it is of value to determine to what degree an effective surveillance approach can simplify outbreak reconstruction. To this end, we performed a parallel analysis in which we assumed that a negative rectal surveillance precluded silent colonization. Incorporating this constraint greatly reduced the number of possible transmission links between patients (fig. S3). However, this reconstruction is not a priori more accurate, underscoring the imperative to determine the true sensitivity and specificity of different surveillance approaches. See Supplementary Materials for details of algorithm implementation.

Supplementary Material

Sup Materials


We thank S. Holland, P. Murray, F. Candotti, and H. Kong for their thoughtful discussions and advice, and T. Connor and J. Parkhill for sharing unpublished sequence data.

Funding: Supported by the National Human Genome Research Institute and NIH Clinical Center Intramural Research Programs and by an NIH Director's Challenge Award for genome sequencing. E.S.S. is supported by a Pharmacology Research Associate Training Fellowship, National Institute of General Medical Sciences.



Author contributions: E.S.S., A.M.Z., D.K.H., T.N.P., and J.A.S. conceived the study. NISC performed genome sequencing. P.J.T. performed bioinformatics annotation and coordinated sequencing studies. F.S. performed repetitive PCR and pulsed-field gel electrophoresis. E.S.S. performed all data analysis. E.S.S., D.K.H., T.N.P., and J.A.S. wrote the manuscript.

Competing interests: The authors declare that they have no competing interests.

Data and materials availability: The results of the Whole Genome Shotgun project have been deposited at DNA Database of Japan/European Molecular Biology Laboratory/GenBank under accession numbers AJZU00000000, AJZV00000000, AJZW00000000, AJZX00000000, AJZY00000000, AJZZ00000000, AKAA00000000, AKAB00000000, AKAC00000000, AKAD00000000, AKAE00000000, AKAF00000000, AKAG00000000, AKAH00000000, AKAI00000000, AKAJ00000000, AKAK00000000, AKAL00000000, AKAM00000000, and AKAN00000000.


1. Lockhart SR, Abramson MA, Beekmann SE, Gallagher G, Riedel S, Diekema DJ, Quinn JP, Doern GV. Antimicrobial resistance among gram-negative bacilli causing infections in intensive care unit patients in the United States between 1993 and 2004. J. Clin. Microbiol. 2007;45:3352–3359. [PMC free article] [PubMed]
2. Montgomerie JZ. Epidemiology of Klebsiella and hospital-associated infections. Rev. Infect. Dis. 1979;1:736–753. [PubMed]
3. Yigit H, Queenan AM, Anderson GJ, Domenech-Sanchez A, Biddle JW, Steward CD, Alberti S, Bush K, Tenover FC. Novel carbapenem-hydrolyzing β-lactamase, KPC-1, from a carbapenem-resistant strain of Klebsiella pneumoniae. Antimicrob. Agents Chemother. 2001;45:1151–1161. [PMC free article] [PubMed]
4. Nordmann P, Cuzon G, Naas T. The real threat of Klebsiella pneumoniae carbapenemase-producing bacteria. Lancet Infect. Dis. 2009;9:228–236. [PubMed]
5. Gupta N, Limbago BM, Patel JB, Kallen AJ. Carbapenem-resistant Enterobacteriaceae: epidemiology and prevention. Clin. Infect. Dis. 2011;53:60–67. [PubMed]
6. Hirsch EB, Tam VH. Detection and treatment options for Klebsiella pneumoniae carbapenemases (KPCs): An emerging cause of multidrug-resistant infection. J. Antimicrob. Chemother. 2010;65:1119–1125. [PubMed]
7. Livermore DM, Warner M, Mushtaq S, Doumith M, Zhang J, Woodford N. What remains against carbapenem-resistant Enterobacteriaceae? Evaluation of chloramphenicol, ciprofloxacin, colistin, fosfomycin, minocycline, nitrofurantoin, temocillin and tigecycline. Int. J. Antimicrob. Agents. 2011;37:415–419. [PubMed]
8. Patel G, Huprikar S, Factor SH, Jenkins SG, Calfee DP. Outcomes of carbapenem-resistant Klebsiella pneumoniae infection and the impact of antimicrobial and adjunctive therapies. Infect. Control Hosp. Epidemiol. 2008;29:1099–1106. [PubMed]
9. Ben-David D, Kordevani R, Keller N, Tal I, Marzel A, Gal-Mor O, Maor Y, Rahav G. Outcome of carbapenem resistant Klebsiella pneumoniae bloodstream infections. Clin. Microbiol. Infect. 2012;18:54–60. [PubMed]
10. Bush K. Alarming β-lactamase-mediated resistance in multidrug-resistant Enterobacteriaceae. Curr. Opin. Microbiol. 2010;13:558–564. [PubMed]
11. Asensio A, Oliver A, González-Diego P, Baquero F, Pérez-Díaz JC, Ros P, Cobo J, Palacios M, Lasheras D, Cantón R. Outbreak of a multiresistant Klebsiella pneumoniae strain in an intensive care unit: Antibiotic use as risk factor for colonization and infection. Clin. Infect. Dis. 2000;30:55–60. [PubMed]
12. Peña C, Pujol M, Ardanuy C, Ricart A, Pallares R, Liñares J, Ariza J, Gudiol F. Epidemiology and successful control of a large outbreak due to Klebsiella pneumoniae producing extended-spectrum β-lactamases. Antimicrob. Agents Chemother. 1998;42:53–58. [PMC free article] [PubMed]
13. Jarvis WR, Munn VP, Highsmith AK, Culver DH, Hughes JM. The epidemiology of nosocomial infections caused by Klebsiella pneumoniae. Infect. Control. 1985;6:68–74. [PubMed]
14. Selden R, Lee S, Wang WLL, Bennet JV, Eickhoff TC. Nosocomial Klebsiella infections: Intestinal colonization as a reservoir. Ann. Intern. Med. 1971;74:657–664. [PubMed]
15. Casewell M, Phillips I. Hands as route of transmission for Klebsiella species. Br. Med. J. 1977;2:1315–1317. [PMC free article] [PubMed]
16. Kitchel B, Rasheed JK, Patel JB, Srinivasan A, Navon-Venezia S, Carmeli Y, Brolund A, Giske CG. Molecular epidemiology of KPC-producing Klebsiella pneumoniae isolates in the United States: Clonal expansion of multilocus sequence type 258. Antimicrob. Agents Chemother. 2009;53:3365–3370. [PMC free article] [PubMed]
17. Diancourt L, Passet V, Verhoef J, Grimont PAD, Brisse S. Multilocus sequence typing of Klebsiella pneumoniae nosocomial isolates. J. Clin. Microbiol. 2005;43:4178–4182. [PMC free article] [PubMed]
18. Won SY, Munoz-Price LS, Lolans K, Hota B, Weinstein RA, Hayden MK. Centers for Disease Control and Prevention Epicenter Program, Emergence and rapid regional spread of Klebsiella pneumoniae carbapenemase–producing Enterobacteriaceae. Clin. Infect. Dis. 2011;53:532–540. [PubMed]
19. Parkhill J, Wren BW. Bacterial epidemiology and biology—Lessons from genome sequencing. Genome Biol. 2011;12:230. [PMC free article] [PubMed]
20. Pallen MJ, Loman NJ. Are diagnostic and public health bacteriology ready to become branches of genomic medicine? Genome Med. 2011;3:53. [PMC free article] [PubMed]
21. Harris SR, Feil EJ, Holden MTG, Quail MA, Nickerson EK, Chantratita N, Gardete S, Tavares A, Day N, Lindsay JA, Edgeworth JD, de Lencastre H, Parkhill J, Peacock SJ, Bentley SD. Evolution of MRSA during hospital transmission and intercontinental spread. Science. 2010;327:469–474. [PMC free article] [PubMed]
22. Mutreja A, Kim DW, Thomson NR, Connor TR, Lee JH, Kariuki S, Croucher NJ, Choi SY, Harris SR, Lebens M, Niyogi SK, Kim EJ, Ramamurthy T, Chun J, Wood JLN, Clemens JD, Czerkinsky C, Nair GB, Holmgren J, Parkhill J, Dougan G. Evidence for several waves of global transmission in the seventh cholera pandemic. Nature. 2011;477:462–465. [PMC free article] [PubMed]
23. Beres SB, Carroll RK, Shea PR, Sitkiewicz I, Martinez-Gutierrez JC, Low DE, McGeer A, Willey BM, Green K, Tyrrell GJ, Goldman TD, Feldgarden M, Birren BW, Fofanov Y, Boos J, Wheaton WD, Honisch C, Musser JM. Molecular complexity of successive bacterial epidemics deconvoluted by comparative pathogenomics. Proc. Natl. Acad. Sci. U.S.A. 2010;107:4371–4376. [PubMed]
24. Chin CS, Sorenson J, Harris JB, Robins WP, Charles RC, Jean-Charles RR, Bullard J, Webster DR, Kasarskis A, Peluso P, Paxinos EE, Yamaichi Y, Calderwood SB, Mekalanos JJ, Schadt EE, Waldor MK. The origin of the Haitian cholera outbreak strain. N. Engl. J. Med. 2011;364:33–42. [PMC free article] [PubMed]
25. 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 SJM, Brinkman FSL, Brunham RC, Tang P. Whole-genome sequencing and social-network analysis of a tuberculosis outbreak. N. Engl. J. Med. 2011;364:730–739. [PubMed]
26. Lienau EK, Strain E, Wang C, Zheng J, Ottesen AR, Keys CE, Hammack TS, Musser SM, Brown EW, Allard MW, Cao G, Meng J, Stones R. Identification of a salmonellosis outbreak by means of molecular sequencing. N. Engl. J. Med. 2011;364:981–982. [PubMed]
27. Rasko DA, Webster DR, Sahl JW, Bashir A, Boisen N, Scheutz F, Paxinos EE, Sebra R, Chin CS, Iliopoulos D, Klammer A, Peluso P, Lee L, Kislyuk AO, Bullard J, Kasarskis A, Wang S, Eid J, Rank D, Redman JC, Steyert SR, Frimodt-Møller J, Struve C, Petersen AM, Krogfelt KA, Nataro JP, Schadt EE, Waldor MK. Origins of the E. coli strain causing an outbreak of hemolytic–uremic syndrome in Germany. N. Engl. J. Med. 2011;365:709–717. [PMC free article] [PubMed]
28. Rasko DA, Worsham PL, Abshire TG, Stanley ST, Bannan JD, Wilson MR, Langham RJ, Decker RS, Jiang L, Read TD, Phillippy AM, Salzberg SL, Pop M, VanErt MN, Kenefic LJ, Keim PS, Fraser-Liggett CM, Ravel J. Bacillus anthracis comparative genome analysis in support of the Amerithrax investigation. Proc. Natl. Acad. Sci. U.S.A. 2011;108:5027–5032. [PubMed]
29. Eyre DW, Golubchik T, Gordon NC, Bowden R, Piazza P, Batty EM, Ip CLC, Wilson DJ, Didelot X, O'Connor L, Lay R, Buck D, Kearns AM, Shaw A, Paul J, Wilcox MH, Donnelly PJ, Peto TEA, Walker AS, Crook DW. A pilot study of rapid benchtop sequencing of Staphylococcus aureus and Clostridium difficile for outbreak detection and surveillance. BMJ Open. 2012;2:e001124. [PMC free article] [PubMed]
30. Lieberman TD, Michel JB, Aingaran M, Potter-Bynoe G, Roux D, Davis MR, Jr., Skurnik D, Leiby N, LiPuma JJ, Goldberg JB, McAdam AJ, Priebe GP, Kishony R. Parallel bacterial evolution within multiple patients identifies candidate pathogenicity genes. Nat. Genet. 2011;43:1275–1280. [PMC free article] [PubMed]
31. Köser CU, Holden MTG, Ellington MJ, Cartwright EJP, 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. [PMC free article] [PubMed]
32. 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. [PubMed]
33. Pränting M, Negrea A, Rhen M, Andersson DI. Mechanism and fitness costs of PR-39 resistance in Salmonella enterica serovar typhimurium LT2. Antimicrob. Agents Chemother. 2008;52:2734–2741. [PMC free article] [PubMed]
34. Mattiuzzo M, Bandiera A, Gennaro R, Benincasa M, Pacor S, Antcheva N, Scocchi M. Role of the Escherichia coli SbmA in the antimicrobial activity of proline-rich peptides. Mol. Microbiol. 2007;66:151–163. [PubMed]
35. Feero WG, Guttmacher AE, Collins FS. Genomic medicine—An updated primer. N. Engl. J. Med. 2010;362:2001–2011. [PubMed]
36. Comas I, Borrell S, Roetzer A, Rose G, Malla B, Kato-Maeda M, Galagan J, Niemann S, Gagneux S. Whole-genome sequencing of rifampicin-resistant Mycobacterium tuberculosis strains identifies compensatory mutations in RNA polymerase genes. Nat. Genet. 2011;44:106–110. [PMC free article] [PubMed]
37. Arias CA, Panesso D, McGrath DM, Qin X, Mojica MF, Miller C, Diaz L, Tran TT, Rincon S, Barbu EM, Reyes J, Roh JH, Lobos E, Sodergren E, Pasqualini R, Arap W, Quinn JP, Shamoo Y, Murray BE, Weinstock GM. Genetic basis for in vivo daptomycin resistance in enterococci. N. Engl. J. Med. 2011;365:892–900. [PMC free article] [PubMed]
38. Otter JA, Yezli S, Schouten MA, van Zanten ARH, Houmes-Zielman G, Nohlmans-Paulssen MKE. Hydrogen peroxide vapor decontamination of an intensive care unit to remove environmental reservoirs of multidrug-resistant gram-negative rods during an outbreak. Am. J. Infect. Control. 2010;38:754–756. [PubMed]
39. Palmore TN, Michelin AV, Bordner M, Odom RT, Stock F, Sinaii N, Fedorko DP, Murray PR, Henderson DK. Use of adherence monitors as part of a team approach to control clonal spread of multidrug-resistant Acinetobacter baumannii in a research hospital. Infect. Control Hosp. Epidemiol. 2011;32:1166–1172. [PubMed]
40. Rissman AI, Mau B, Biehl BS, Darling AE, Glasner JD, Perna NT. Reordering contigs of draft genomes using the Mauve aligner. Bioinformatics. 2009;25:2071–2073. [PMC free article] [PubMed]
41. Darling AC, Mau B, Blattner FR, Perna NT. Mauve: Multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 2004;14:1394–1403. [PubMed]
42. Jombart T, Eggo RM, Dodd PJ, Balloux F. Reconstructing disease outbreaks from genetic data: A graph approach. Heredity. 2011;106:383–390. [PMC free article] [PubMed]
43. Wiener-Well Y, Rudensky B, Yinnon AM, Kopuit P, Schlesinger Y, Broide E, Lachish T, Raveh D. Carriage rate of carbapenem-resistant Klebsiella pneumoniae in hospitalised patients during a national outbreak. J. Hosp. Infect. 2010;74:344–349. [PubMed]
44. Fouts DE. Phage_Finder: Automated identification and classification of prophage regions in complete bacterial genome sequences. Nucleic Acids Res. 2006;34:5839–5851. [PubMed]
45. Delcher AL, Phillippy A, Carlton J, Salzberg SL. Fast algorithms for large-scale genome alignment and comparison. Nucleic Acids Res. 2002;30:2478–2483. [PMC free article] [PubMed]
46. Edmonds J. Optimum branchings. J. Res. Natl. Bur. Standards. 1967;9:233–240.