ABSTRACT
Nosocomial infections pose a significant threat to patient health; however, the gold standard laboratory method for determining bacterial relatedness (pulsed-field gel electrophoresis [PFGE]) remains essentially unchanged 20 years after its introduction. Here, we explored bacterial whole-genome sequencing (WGS) as an alternative approach for molecular strain typing. We compared WGS to PFGE for investigating presumptive outbreaks involving three important pathogens: vancomycin-resistant Enterococcus faecium (n = 19), methicillin-resistant Staphylococcus aureus (n = 17), and Acinetobacter baumannii (n = 15). WGS was highly reproducible (average ≤ 0.39 differences between technical replicates), which enabled a functional, quantitative definition for determining clonality. Strain relatedness data determined by PFGE and WGS roughly correlated, but the resolution of WGS was superior (P = 5.6 × 10−8 to 0.016). Several discordant results were noted between the methods. A total of 28.9% of isolates which were indistinguishable by PFGE were nonclonal by WGS. For A. baumannii, a species known to undergo rapid horizontal gene transfer, 16.2% of isolate pairs considered nonidentical by PFGE were clonal by WGS. Sequencing whole bacterial genomes with single-nucleotide resolution demonstrates that PFGE is prone to false-positive and false-negative results and suggests the need for a new gold standard approach for molecular epidemiological strain typing.
INTRODUCTION
Hospital-acquired infections cause significant patient morbidity and mortality around the world (1–3) and incur substantial burdens with respect to health care costs; current estimates exceed $9.7 billion yearly in the United States alone (3). In addition to preventative measures (2), effective infection control is achieved through active surveillance of infection and recognition of burgeoning disease outbreaks that can be interrupted through health care interventions (1). Molecular epidemiological investigations have become increasingly integral to these functions and rely on methods that can rapidly and effectively assess genetic relatedness among strains (1, 4, 5).
Pulsed-field gel electrophoresis (PFGE) has long been considered the gold standard approach for such investigations (6). Developed at a time when phenotypic strain typing approaches were common (7, 8), PFGE was among the first technologies to propel the era of molecular epidemiology (8). Broadly applicable to bacteria of many different genera, PFGE entails the liberation of intact genomic DNA by in situ lysis within an agarose plug, restriction enzyme digestion, and separation of fragments using pulsed-field gel electrophoresis (7). Pairwise comparisons are made among strains' banding patterns, which are interpreted according to the number of dissimilar restriction fragments that are observed: “indistinguishable” (no difference), “closely related,” “possibly related,” and “different” (signaling the greatest amount of strain divergence) (7). Although alternative methods have been developed to enable molecular typing of bacteria (1, 9, 10) or of particular species (11), none has displaced PFGE as a gold standard approach (8), in part because the approach provides a measurement of differences spanning the entire genome (1).
Despite its utility, there are legitimate criticisms of PFGE. PFGE is both time-consuming and labor-intensive and is technically challenging for many multipurpose clinical laboratories to implement (12, 13). Reproducibility among different laboratories can be difficult to achieve (8, 13–15), and disparities between predicted and achieved PFGE banding patterns have been reported (16). Processing isolates using different restriction enzymes can yield distinct PFGE patterns and, subsequently, discordant interpretations of strain relationships (16). Genetic diversity in closely related outbreak strains (1, 16) or clonally dominant endemic strains (17) may be insufficient to permit their correct interpretation, and some isolates with different phenotypes share the same PFGE pattern (18). PFGE results may correlate poorly with the true relatedness of isolates (16, 19, 20), especially among more distantly related strains (21), and at least one study reported that greater discriminatory power could be obtained through alternative molecular approaches (22). Due in part to these concerns, alternative, amplification-based molecular methods for strain typing of bacteria have been developed, including multilocus variable-number tandem-repeat analysis (MLVA) (23), multilocus sequence typing (MLST) (24), and randomly amplified polymorphic DNA (RAPD) typing (25).
Whole-genome sequencing (WGS) represents a relatively new and increasingly accessible means for tracking disease outbreaks that has garnered success in multiple applied contexts (26–34). Using massively parallel (or “next-generation”) DNA sequencing technologies, it is now possible to examine the complete or nearly complete genomes of bacterial isolates. WGS can theoretically distinguish strains which differ at only a single nucleotide and, in the limited number of studies where direct comparisons with PFGE have been performed, has provided greater resolution (28, 35, 36). However, previous work has focused on individual bacterial species, has employed nonuniform validation methodologies, and has utilized heterogeneous library preparation, sequencing, and data analysis methods, making it difficult to integrate and generalize results across studies or to evaluate its suitability for clinical use.
Here we explored the utility of WGS as a strain typing approach for the clinical laboratory using a single, universal protocol (encompassing library preparation, sequencing, and data analysis) for 3 biologically and genomically distinct organisms: methicillin-resistant Staphylococcus aureus (MRSA), vancomycin-resistant Enterococcus (VRE), and multidrug-resistant Acinetobacter baumannii. We examined isolates sent to our laboratory for investigation of suspected clinical outbreaks and systematically evaluated WGS results with matched PFGE data.
MATERIALS AND METHODS
Strains.Three different clinically relevant organisms with phenotypical and genomic properties were selected for this study (Table 1). All strain collections were part of outbreak investigations and were collected from patients with suspected epidemiological relationships. MRSA isolates were derived from sputum (15 isolates), stool (1 isolate), and urine (1 isolate). A. baumannii strains originated from sputum (7 isolates), bronchoalveolar lavage fluid (3 isolates), tissue (2 isolates), a wound (1 isolate), a chest tube drain (1 isolate), and urine (1 isolate). Information on the source of VRE strains was not available. For PFGE analyses, all strains were freshly harvested from primary, single-colony isolates at the time of the outbreak investigation. Strains were cryopreserved until the time of whole-genome sequencing. All strains were maintained on blood agar at 37°C when actively grown. Control strains not spatiotemporally related to the outbreak investigation were selected for each strain collection and were included in all analyses.
Characteristics of organisms studied
PFGE.Agarose plugs were prepared using a rapid method as previously described (37). Genomic DNA was digested with apaI (New England BioLabs, Ipswich, MA) for VRE and A. baumannii or with smaI (New England BioLabs) for MRSA. Restriction fragments were resolved by PFGE using a temperature-controlled Chef-DR III system (Bio-Rad, Hercules, CA). Gels were photographed using a Bio-Rad Gel Doc 1000 system and were interpreted according to standard guidelines (7).
Whole-genome sequencing.DNA was purified using an UltraClean microbial DNA isolation kit (Mo-Bio, Carlsbad, CA). Genomic DNA (30 to 100 ng) was digested for 60 to 90 min at 37°C in a 10-μl volume using 0.3 μl NEBNext double-stranded DNA (dsDNA) Fragmentase (New England BioLabs). DNA was end repaired and A-tailed in a 40-μl reaction mixture containing 1× Rapid ligation buffer (Enzymatics Inc., Beverly, MA), 0.1675 mM (each) deoxynucleoside triphosphate (dNTP) (New England BioLabs), 0.1 μl Escherichia coli DNA polymerase I (New England BioLabs), 0.5 μl T4 polynucleotide kinase (New England BioLabs), and 0.02 μl Taq DNA polymerase (New England BioLabs) and was incubated at 37°C for 30 min and 72°C for 20 min. Annealed Y-adaptors (0.2 μM) (5′-[PO4]GATCGGAAGAGCGGTTCAGCAGGAATGCCGAG-3′ and 5′-ACACTCTTTCCCTACACGACGCTCTTCCGATCT-3′) were added and ligated at 25°C for 20 min using T4 DNA ligase in Rapid ligation buffer (Enzymatics Inc.). Following purification with Agencourt AMPure XP beads (Beckman Coulter, Brea, CA), PCR amplification was performed with KAPA HiFi HotStart ReadyMix (Kapa Biosystems, Wilmington, MA) using primer PRECAP_FWD_AMP_COMMON (5′-AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGC-3′) and sample-specific indexed primers (5′-CAAGCAGAAGACGGCATACGAGATXXXXXXXXCGGTCTCGGCATTCCTGCTGAACCG-3′, where XXXXXXXX indicates the 8-bp index). Cycling conditions were 95°C for 3 min, 10 cycles of 98°C for 20 s, 65°C for 15 s, and 72°C for 1 min, and 1 cycle of 72°C for 5 min. PCR products were purified using AMPure beads, pooled in equimolar amounts, and sequenced on a MiSeq sequencer (Illumina, San Diego, California) using 150-bp paired-end reads with a custom index primer (5′-AGATCGGAAGAGCGGTTCAGCAGGAATGCCGAGACCG-3′). Oligonucleotides were synthesized by IDT. Average read depths obtained are summarized in Table 1.
Data analysis.Adaptors were trimmed and PCR duplicates removed using Fastq-Mcf (http://code.google.com/p/ea-utils/) with skew filtering disabled and other parameters at their defaults. Sequence reads were aligned to completed reference genomes for Enterococcus faecium DO (GenBank identifier [ID] 389867183), Acinetobacter baumannii strain AYE (GenBank ID 169147133), or Staphylococcus aureus USA300 FPR3757 (GenBank ID 87125858) using bwa (v0.6.2) (38) and samtools (v0.1.19) (39). Reads with a mapping quality value of less than 10 were discarded. Calling of single nucleotide variants and small insertions and deletions (indels) was performed using samtools with a haploid genome model and a minimum variant frequency value of 0.5. Variants supported by fewer than 15 reads or a likelihood score of less than 200 were masked as “unknown” data. All-by-all pairwise distance matrices were constructed by comparing sites of variation among isolates, masking sites at which one or both isolates displayed “unknown” data or less than 15× read coverage, and counting only variant sites for which both isolates in the comparison could be confidently genotyped. Pairwise distances were expressed as the absolute number of passing variant sites which distinguished such pairs, with indels weighted the same as single nucleotide variants. Technical replicates were separately considered in order to determine assay variability, and then paired replicates were merged prior to performing intraisolate comparisons. Approximately maximum-likelihood phylogenetic trees were constructed using FastTree 2.1 (40) (see Fig. S1 in the supplemental material).
Analysis of discrepant strains.One VRE isolate evidencing poor mapping statistics despite adequate sequence reads was subjected to additional analysis. A de novo assembly was produced for this isolate using AbySS v1.3.5 (41), with parameters empirically optimized to maximize the N50 statistic (the length for which all contigs of equal or larger size contain half the sum of the entire assembly). To perform taxonomic classification of the isolate, contigs in the assembly were compared to the NCBI nonredundant nucleotide database using a BLAST (42) search. The top match for each contig was recorded, as well as the length of that contig in nucleotides. The organism for which the greatest cumulative numbers of nucleotides were identified as the top matches was examined, as were matches corresponding to 70% or more of the matched nucleotides associated with the top hit.
Nucleotide sequence accession number.Sequence reads generated in this study are available from the NCBI Sequence Read Archive (SRA; http://www.ncbi.nlm.nih.gov/sra) under accession no. SRP042634.
RESULTS
Reproducibility of WGS.Inaccuracies in whole-genome sequence data may result from errors that occur during library preparation (43), sequencing (44), and bioinformatic analysis (45). We limited the numbers of artifacts through laboratory procedures, including the use of high-fidelity DNA polymerase and limiting PCR cycling during library preparation, and performed the sequencing analysis using the highly accurate Illumina platform (46). Bioinformatics artifacts were minimized by discarding reads with low mapping quality, using a model-based variant caller tuned to a haploid genome (39), and by imposing quality filtering on variant calls.
We assessed the reproducibility of WGS using technical replicates; duplicate sequencing libraries prepared from separate aliquots of the same genomic DNA were sequenced and analyzed in parallel (Table 2). Technical replicates were produced for all specimens utilized in this study (51 isolates), allowing robust assessment of the inherent variability of the assay. Of the 51 paired technical replicates, 35 were 100% identical at the genomic sequence level. In general, duplicate genome sequences of bacterial isolates differed by less than a single variant, with an average of 0.39 (standard deviation of 0.63) variants distinguishing among paired controls. We did not observe a statistically significant difference in the burden of technical errors in comparing organisms with the highest and lowest frequencies of discordant replicates (A. baumannii and VRE, respectively; P = 0.21 by a 2-tailed t test), indicating that the levels of reproducibility of the assay were similar across all three species. There was only a weak correlation between genome size and the rate of technical errors (R2 = 0.69). We conclude that our protocol is highly reproducible and generalizable across the three different organisms studied.
Reproducibility of WGS
We consequently set a threshold for interpreting isolates as truly clonal, or genomically indistinguishable, corresponding to detection of 3 or fewer pairwise differences between them (equal to the average technical replicate error frequency + 4.15 standard deviations). Given a normal distribution for the number of sequencing errors encountered per isolate, this threshold is sufficiently sensitive that instances in which isolate differences were attributable to technical artifacts alone would be limited to less than <0.0017%.
Correlation of WGS and PFGE.For each of the strain collections, we considered every possible pairwise comparison of isolates. We identified the PFGE category assignment for each pairwise comparison and tabulated the number of genomic variants identified by WGS which distinguished between the two isolates (Fig. 1 and Table 3).
Numbers of disparate genomic polymorphisms detected among pairwise comparisons of technical replicates and different PFGE classification categories. Pairwise comparisons of isolates are stratified according to their status as whole-genome sequencing technical replicates (T) or classification by PFGE as “indistinguishable” (I), “closely related” (C), “possibly related” (P), or “different” (D). The numbers of genomic variants identified by whole-genome sequencing which distinguished isolate pairs are plotted along the y axis. Results are separately shown for VRE (A), MRSA (B), and A. baumannii (C). Red shading indicates a pairwise comparison with a count of genomic differences within 1 standard deviation of the mean observed in a different PFGE category for the same organism.
Numbers of genomic polymorphisms detected among PFGE categories
The number of genomic differences separating technical replicates was substantially lower than that observed among strains categorized as “indistinguishable” by PFGE for all organisms studied (VRE, P = 5.6 × 10−8; MRSA, P = 0.0004; A. baumannii, P = 0.007 [using 2-tailed t tests]). This finding indicates that WGS is significantly better able to resolve closely related isolates than PFGE.
Additionally, for all organisms, the number of variants identified by WGS correlated roughly with the degree of dissimilarity determined by PFGE. However, the relationship between the two measures was not directly proportional.
In other regards, we observed apparently species-specific differences.
VRE.VRE demonstrated the most concordant results of comparisons between PFGE and WGS (Fig. 1A and Table 3). The average number of genomic variants distinguishing strains increased consistently with the degree of dissimilarity registered by PFGE. PFGE categories were discrete and well-defined by the number of genomic differences they displayed (“indistinguishable” versus “closely related,” P = 7.6 × 10−9;“closely related” versus “possibly related,” P = 1.9 × 10−6; “possibly related” versus “different,” P = 3.9 × 10−14), with no overlap observed among isolates in different PFGE categories (as defined by isolates demonstrating a number of polymorphisms falling within 1 standard deviation of the mean observed for a different PFGE category).
One VRE isolate demonstrated poor (<0.05%) alignment coverage of the reference genome despite adequate numbers of sequencing reads (>667,000) and was excluded from primary analysis. Unlike the other isolates, further investigation revealed that this organism best matched Enterococcus faecalis, rather than E. faecium. Although the PFGE pattern derived from the strain marked it as different from the others, identifying the erroneous taxonomic classification of this organism would not have been achievable by PFGE.
MRSA.The average number of genomic polymorphisms that distinguished pairs of MRSA strains correlated less well with the PFGE typing category than was the case for VRE strains, and the distribution of variants observed in each category was greater than that observed for VRE strains (Fig. 1B and Table 3). The average number of pairwise differences among strains classified as “indistinguishable” was considerably higher than that observed for VRE isolates bearing the same PFGE designation (P = 6.9 × 10−4 [2-tailed t test]). Paired comparisons of “possibly related” isolate pairs had on average significantly fewer genomic differences than paired comparisons among strains classified as “closely related” by PFGE (149.2 variants versus 467.1 variants; P = 1.9 × 10−6 [2-tailed t test]), a category defined as demonstrating less divergence.
The numbers of genomic differences separating isolate pairs within specific PFGE groups showed significant overlap. A total of 10 pairwise comparisons in the “closely related” category (43.5% of pairwise comparisons in that category) had a count of genomic differences (range, 152 to 264 variants) that fell within 1 standard deviation of the mean observed for the “possibly related” category (149.22 ± 119.2 variants). Four comparisons in the “different” category (11.4% of such pairwise comparisons; range, 118 to 160 variants) fell within this range and, surprisingly, were also less than 1 standard deviation from the mean for the “closely related” category (232.87 variants). Similarly, 4 “possibly related” isolate pairs (6.9%; range, 5 to 32 variants) were within 1 standard deviation of the mean observed in the “indistinguishable” category (18.15 ± 17.81 variants).
The distributions of pairwise genomic distances observed for the “indistinguishable,” “closely related,” and “different” categories were bimodal (Fig. 1B), indicating an underlying population structure, not evident by PFGE, that is readily detected by WGS with single-nucleotide resolution (see Fig. S1 in the supplemental material).
A. baumannii.Although summary statistics for the number of genomic differences defining A. baumannii PFGE categories were comparable to those calculated for the other organisms (Table 3), the PFGE results were discordant with the whole-genome data (Fig. 1C).
The numbers of genomic differences that differentiated PFGE categories were least discrete for A. baumannii, and there was substantial overlap in the numbers of polymorphisms demonstrated among separate PFGE categories. The “indistinguishable” and “closely related” PFGE classifications were not significantly different based on the numbers of genomic polymorphisms they exhibited (P = 0.206 [2-tailed t test]). For all 5 (100%) isolate pairs classified as “indistinguishable,” the counts of the distinguishing polymorphisms that they displayed (range, 3 to 8) were within 1 standard deviation of the mean for values exhibited by the “closely related” classification (5.77 ± 3.27 variants). Counts of distinguishing genomic variants for all 5 pairwise comparisons from the “indistinguishable” category and all 44 (100%) paired comparisons deemed in the “closely related” category (range, 0 to 13 variants) fell within 1 standard deviation of the mean value observed for the “possibly related” PFGE classification (296.96 ± 380.45 variants). The “different” category was the category best separated from other PFGE groups; however, two pairs of outliers (793 and 794 variants) fell within the range of genomic differences compatible with the “possibly related” group.
The pairwise comparisons in the “possibly related” category were distributed almost equally across two groups (Fig. 1C). One group was characterized by a higher number of distinguishing mutations and the other by a much lower number, similar to the number of differences observed among the “closely related” PFGE classification.
Evaluation of clonality by PFGE and WGS.Finally, we compared the abilities of WGS and PFGE to resolve clonal isolates (Table 4), a critical determination in outbreak investigations (1).
Numbers of discrepancies between whole-genome sequencing and PFGE for paired-strain comparisons
Given our empirical definition of clonality, a substantial fraction (28.9%) of the total pairwise comparisons interpreted as “indistinguishable” by PFGE were nonclonal by WGS. This fraction ranged from a low of 14.1% for VRE to 75.0% for MRSA (P < 0.0001 [2-sample Z-test]), suggesting organism-specific bias.
None of the pairwise strain comparisons deemed “different” by PFGE were defined as clonal by WGS for any of the three organisms. For MRSA and VRE, WGS did not identify clonal relationships among the pairs characterized as “closely related,” “possibly related,” or “different” by PFGE. For A. baumannii, in contrast, 27.3% of the strains identified as “closely related” by PFGE and 14.8% of the “possibly related” strains were clonal by WGS.
DISCUSSION
Here we compared the performance of WGS to that of PFGE for bacterial strain typing, using a scalable and species-independent protocol for sequencing and analysis. The study revealed a number of important differences between these two methods with respect to their suitability for application in clinical molecular epidemiology investigations.
Clinical surveillance and investigation of nosocomial outbreaks are primarily dependent on determining whether biogeographically related pairs of strains are clonal and thus on whether they could have originated from the same source or are genomically distinct and could therefore reflect independent transmission events (1). To enable derivation of a functional definition of clonality (genomic identity) from whole-genome sequence data, we set a maximum of 3 distinguishing variants in order to rule out spurious differences resulting from technical artifacts. Although this quantitative definition should be applicable across many different bacterial species, in theory, the probability of sequencing errors should be proportional to the genome size and, to some extent, to the genomic architecture and composition (43); thus, it may need to be altered for evolutionarily divergent or unusual species. The error rates observed in this study were sufficiently low (Table 3) that the results from most (82%) of the technical replicates were 100% identical. The error rates were nonlimiting in our study, and yet further improvements may be achieved by incorporating methods such as PCR-free library preparation (43) or limitation of variant calls to particular portions of microbial genomes (29), though the latter strategy could theoretically negatively impact sensitivity.
Regardless, the issue of what truly constitutes a “bacterial clone” (47) or a “clonal outbreak” (48) is surprisingly complex. Strains may undergo diversification over the course of infection within an individual patient (31); thus, isolates which are genetically distinct but very closely related may reflect a common origin. Presently, there is no established threshold defining the number of genome-wide polymorphisms which identify isolates as belonging to an outbreak (48). Various studies have estimated that the genomes of most bacterial strains accumulate single nucleotide polymorphisms at a rate of roughly 2 to 10 per year, depending on the organism (31, 32, 48, 49). However, organism-specific guidelines for interpreting WGS data for molecular epidemiology investigations will require dedicated, large scale, and long-term observational studies of defined outbreaks among patients, involving multiple isolates from the same patient, in order to determine the rates at which particular species accrue variants during carriage in human hosts (28, 31). For the time being, and in keeping with the existing classification framework developed for PFGE, we propose that strain comparisons based on WGS should be interpreted generally with respect to the following three actionable outcomes: strains are to be considered “genomically indistinguishable” if they are separated by 3 or fewer variants, “closely related” if they are separated by up to 12 variants (based on existing studies of inter- and intrapatient strain variability [31]), and “unrelated” if they are distinguishable by 13 or more variants. As with interpreting the results of any strain typing, molecular epidemiological investigations will require integration of these classifications with additional information and context, notably spatiotemporal information about isolate collection, in order to determine the plausibility of potential transmission events (5). Detailed phylogenomic reconstructions of strain relationships (see Fig. S1 in the supplemental material) are achievable with whole-genome sequence data and provide insight into the likely order of transmission events (26), providing additional information that is useful in outbreak investigations when coupled with biogeographic and temporal data.
It is worth noting differences between WGS and PFGE with respect to DNA acquisition events. Horizontal gene transfer among strains, mediated by plasmids and other mobile genetic elements, and recombination may significantly alter the genomic and phenotypic properties of a strain through single, discrete events. In the context of molecular strain typing, such events may result in shifts to PFGE patterns that are not proportional to the time or degree of strain divergence. In contrast, point mutations detectable by WGS accumulate at a far more predictable rate (47, 50) and consequently serve as a more reliable molecular clock for molecular epidemiology reconstructions (51). Although more-sophisticated whole-genome sequence analysis methods could potentially identify and catalog novel DNA elements acquired among strains (52), on principle alone, assessing the continuum of strain relatedness may consequently be best achieved by WGS and comparison of shared genomic regions among strains.
As a case in point, A. baumannii is known to undergo significant amounts of horizontal gene transfer, even over brief timescales that are compatible with disease outbreaks (53). In comparison to the other organisms studied here, we observed an elevated number of discrepancies between WGS and PFGE results for A. baumannii (Table 4). Of particular note, a substantial fraction of the pairwise comparisons deemed “closely related” by PFGE (12 of 44) were clonal by next-generation sequencing (Table 4). Horizontal DNA transfer among isolates or from external sources is a plausible explanation for these inconsistencies: identification of single nucleotide variants and indels by WGS would ignore changes resulting from gene acquisition or loss, effectively circumventing the confounding effects of horizontal gene transfer, whereas, in contrast, PFGE would be more susceptible to them (47), with the unintended consequence of exaggerating the magnitude of divergence among strains.
Our data indicate that WGS has significantly improved resolving power compared to PFGE (Tables 2 and 3). On average, 467 mutations had accumulated between pairs of MRSA strains, and 822 mutations between pairs of VRE isolates, before a difference was registered by PFGE. A minimum of 152 and 767 genomic polymorphisms were observed among pairs of MRSA and VRE isolates resolved as nonidentical by PFGE, respectively, while up to 47 variants distinguished MRSA isolate pairs categorized as “indistinguishable” by PFGE. Restriction patterns reflect variations in only a relatively small fraction of an organism's genome (1), and our findings suggest that relatively large numbers of sequence variants can accumulate in bacterial genomes before the underlying differences are registered by PFGE.
Importantly, in light of the numbers of genomic polymorphisms detectable by WGS that were able to distinguish strains deemed identical by PFGE (Tables 3 and 4), we conclude that PFGE is prone to false-positive results (i.e., to declaring genomically distinct strains to be clonally related). In the analysis of A. baumannii isolates, PFGE also demonstrated results that are consistent with false negatives (strains which are essentially genomically identical but which are classified as different from one another). In the context of infection control, either kind of erroneous result would carry adverse consequences. False-negative results would prevent the identification of true outbreaks, an obvious threat to patient well-being. False-positive results may conversely indicate an outbreak or nosocomial transmission when there is actually none, resulting in unnecessary patient isolation and other interventions (54). Such risks are not merely hypothetical: recent work suggests that most MRSA infections observed in hospitals are not the result of nosocomial transmission but rather originate from independent infections by distinct strains (29, 33), a conclusion which is anecdotally consistent with the investigation of a suspected MRSA outbreak in this study (Table 4). Similar patterns have been reported for Clostridium difficile infections (55).
Given these results, we argue that WGS should replace PFGE as the gold standard method for bacterial strain typing in molecular epidemiology applications. Analytically, WGS appears to offer superior resolution of strain types and is less prone to false-positive and false-negative findings. Further, genomic differences distinguishing strains can be precisely measured and are highly reproducible, allowing high-resolution inference of phylogenomic relationships (see Fig. S1 in the supplemental material) and of the sequence of transmission events. There are additional practical benefits and considerations. The digital nature of sequencing data allows sharing of sequence information among laboratories, potentially enabling a system capable of identifying large-scale outbreaks analogous to existing PFGE databases such as PulseNet (6). Genomic-sequencing libraries are prepared with protocols similar to current PCR-based strain typing protocols (52), requiring less technical skill than PFGE. WGS can enable exploration of isolates' virulence genes, antibiotic resistance mechanisms, and other medically relevant factors (56), concordantly with molecular epidemiology investigation. The per-sample cost of next-generation sequencing drops rapidly when amortized over even a modest number of specimens: here, up to 30 isolates were multiplexed for a total per-sample reagent cost of ∼$35. WGS and analysis can be completed rapidly and are likely to become more rapid (57), potentially enabling real-time epidemiological investigations: the methods implemented here are compatible with a turnaround time of 3 working days, and the computational analysis methods implemented here are scalable to encompass even large numbers of specimens. Although present barriers to the universal adoption of WGS by clinical laboratories include relatively high costs of instrumentation and a lack of bioinformatic expertise (52, 58), these beneficial capabilities of WGS will potentially enable epidemiological investigation of bacterial outbreaks in real time and with an unprecedented ability to define strain relationships, significantly impacting the practice of infection control and patient outcomes and justifying the additional efforts required for its implementation.
ACKNOWLEDGMENT
This work was supported in part by grant UL1TR000423 from the National Center For Advancing Translational Sciences of the National Institutes of Health.
FOOTNOTES
- Received 2 December 2014.
- Returned for modification 13 January 2015.
- Accepted 20 January 2015.
- Accepted manuscript posted online 28 January 2015.
- Address correspondence to Stephen J. Salipante, stevesal{at}uw.edu.
D.J.S. and L.A.C. contributed equally to this article.
Citation Salipante SJ, SenGupta DJ, Cummings LA, Land TA, Hoogestraat DR, Cookson BT. 2015. Application of whole-genome sequencing for bacterial strain typing in molecular epidemiology. J Clin Microbiol 53:1072–1079. doi:10.1128/JCM.03385-14.
Supplemental material for this article may be found at http://dx.doi.org/10.1128/JCM.03385-14.
REFERENCES
- Copyright © 2015, American Society for Microbiology. All Rights Reserved.
