Application of Whole-Genome Sequencing for Bacterial Strain Typing in Molecular Epidemiology

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⁻⁸ 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.     

Whole-Genome Sequencing for High-Resolution Investigation of Methicillin-Resistant Staphylococcus aureus Epidemiology and Genome Plasticity

Methicillin-resistant Staphylococcus aureus (MRSA) infections pose a major challenge in health care, yet the limited heterogeneity within this group hinders molecular investigations of related outbreaks. Pulsed-field gel electrophoresis (PFGE) has been the gold standard approach but is impractical for many clinical laboratories and is often replaced with PCR-based methods. Regardless, both approaches can prove problematic for identifying subclonal outbreaks. Here, we explore the use of whole-genome sequencing for clinical laboratory investigations of MRSA molecular epidemiology. We examine the relationships of 44 MRSA isolates collected over a period of 3 years by using whole-genome sequencing and two PCR-based methods, multilocus variable-number tandem-repeat analysis (MLVA) and spa typing. We find that MLVA offers higher resolution than spa typing, as it resolved 17 versus 12 discrete isolate groups, respectively. In contrast, whole-genome sequencing reproducibly cataloged genomic variants (131,424 different single nucleotide polymorphisms and indels across the strain collection) that uniquely identified each MRSA clone, recapitulating those groups but enabling higher-resolution phylogenetic inferences of the epidemiological relationships. Importantly, whole-genome sequencing detected a significant number of variants, thereby distinguishing between groups that were considered identical by both spa typing (minimum, 1,124 polymorphisms) and MLVA (minimum, 193 polymorphisms); this suggests that these more conventional approaches can lead to false-positive identification of outbreaks due to inappropriate grouping of genetically distinct strains. An analysis of the distribution of variants across the MRSA genome reveals 47 mutational hot spots (comprising ∼2.5% of the genome) that account for 23.5% of the observed polymorphisms, and the use of this selected data set successfully recapitulates most epidemiological relationships in this pathogen group.     

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.     

Characterization of Foodborne Outbreaks of Salmonella enterica Serovar Enteritidis with Whole-Genome Sequencing Single Nucleotide Polymorphism-Based Analysis for Surveillance and Outbreak Detection

Salmonella enterica serovar Enteritidis is a significant cause of gastrointestinal illness in the United States; however, current molecular subtyping methods lack resolution for this highly clonal serovar. Advances in next-generation sequencing technologies have made it possible to examine whole-genome sequencing (WGS) as a potential molecular subtyping tool for outbreak detection and source trace back. Here, we conducted a retrospective analysis of S. Enteritidis isolates from seven epidemiologically confirmed foodborne outbreaks and sporadic isolates (not epidemiologically linked) to determine the utility of WGS to identify outbreaks. A collection of 55 epidemiologically characterized clinical and environmental S. Enteritidis isolates were sequenced. Single nucleotide polymorphism (SNP)-based cluster analysis of the S. Enteritidis genomes revealed well supported clades, with less than four-SNP pairwise diversity, that were concordant with epidemiologically defined outbreaks. Sporadic isolates were an average of 42.5 SNPs distant from the outbreak clusters. Isolates collected from the same patient over several weeks differed by only two SNPs. Our findings show that WGS provided greater resolution between outbreak, sporadic, and suspect isolates than the current gold standard subtyping method, pulsed-field gel electrophoresis (PFGE). Furthermore, results could be obtained in a time frame suitable for surveillance activities, supporting the use of WGS as an outbreak detection and characterization method for S. Enteritidis.     

Tracking Nosocomial Klebsiella pneumoniae Infections and Outbreaks by Whole-Genome Analysis: Small-Scale Italian Scenario within a Single Hospital

Multidrug-resistant (MDR) Klebsiella pneumoniae is one of the most important causes of nosocomial infections worldwide. After the spread of strains resistant to beta-lactams at the end of the previous century, the diffusion of isolates resistant to carbapenems and colistin is now reducing treatment options and the containment of infections. Carbapenem-resistant K. pneumoniae strains have spread rapidly among Italian hospitals, with four subclades of pandemic clonal group 258 (CG258). Here we show that a single Italian hospital has been invaded by three of these subclades within 27 months, thus replicating on a small scale the “Italian scenario.” We identified a single clone responsible for an epidemic outbreak involving seven patients, and we reconstructed its star-like pattern of diffusion within the intensive care unit. This epidemiological picture was obtained through phylogenomic analysis of 16 carbapenem-resistant K. pneumoniae isolates collected in the hospital during a 27-month period, which were added to a database of 319 genomes representing the available global diversity of K. pneumoniae strains. Phenotypic and molecular assays did not reveal virulence or resistance determinants specific for the outbreak isolates. Other factors, rather than selective advantages, might have caused the outbreak. Finally, analyses allowed us to identify a major subclade of CG258 composed of strains bearing the yersiniabactin virulence factor. Our work demonstrates how the use of combined phenotypic, molecular, and whole-genome sequencing techniques can help to identify quickly and to characterize accurately the spread of MDR pathogens.     

Delineating Community Outbreaks of Salmonella enterica Serovar Typhimurium by Use of Whole-Genome Sequencing: Insights into Genomic Variability within an Outbreak

Whole-genome next-generation sequencing (NGS) was used to retrospectively examine 57 isolates from five epidemiologically confirmed community outbreaks (numbered 1 to 5) caused by Salmonella enterica serovar Typhimurium phage type DT170. Most of the human and environmental isolates confirmed epidemiologically to be involved in the outbreaks were either genomically identical or differed by one or two single nucleotide polymorphisms (SNPs), with the exception of those in outbreak 1. The isolates from outbreak 1 differed by up to 12 SNPs, which suggests that the food source of the outbreak was contaminated with more than one strain while each of the other four outbreaks was caused by a single strain. In addition, NGS analysis ruled in isolates that were initially not considered to be linked with the outbreak, which increased the total outbreak size by 107%. The mutation process was modeled by using known mutation rates to derive a cutoff value for the number of SNP difference to determine whether or not a case was part of an outbreak. For an outbreak with less than 1 month of ex vivo/in vivo evolution time, the maximum number of SNP differences between isolates is two or four using the lowest or highest mutation rate, respectively. NGS of S. Typhimurium significantly increases the resolution of investigations of community outbreaks. It can also inform a more targeted public health response by providing important supplementary evidence that cases of disease are or are not associated with food-borne outbreaks of S. Typhimurium.     

Comparing Whole-Genome Sequencing with Sanger Sequencing for spa Typing of Methicillin-Resistant Staphylococcus aureus

spa typing of methicillin-resistant Staphylococcus aureus (MRSA) has traditionally been done by PCR amplification and Sanger sequencing of the spa repeat region. At Hvidovre Hospital, Denmark, whole-genome sequencing (WGS) of all MRSA isolates has been performed routinely since January 2013, and an in-house analysis pipeline determines the spa types. Due to national surveillance, all MRSA isolates are sent to Statens Serum Institut, where the spa type is determined by PCR and Sanger sequencing. The purpose of this study was to evaluate the reliability of the spa types obtained by 150-bp paired-end Illumina WGS. MRSA isolates from new MRSA patients in 2013 (n = 699) in the capital region of Denmark were included. We found a 97% agreement between spa types obtained by the two methods. All isolates achieved a spa type by both methods. Nineteen isolates differed in spa types by the two methods, in most cases due to the lack of 24-bp repeats in the whole-genome-sequenced isolates. These related but incorrect spa types should have no consequence in outbreak investigations, since all epidemiologically linked isolates, regardless of spa type, will be included in the single nucleotide polymorphism (SNP) analysis. This will reveal the close relatedness of the spa types. In conclusion, our data show that WGS is a reliable method to determine the spa type of MRSA.     

Prospective Whole-Genome Sequencing Enhances National Surveillance of Listeria monocytogenes

Whole-genome sequencing (WGS) has emerged as a powerful tool for comparing bacterial isolates in outbreak detection and investigation. Here we demonstrate that WGS performed prospectively for national epidemiologic surveillance of Listeria monocytogenes has the capacity to be superior to our current approaches using pulsed-field gel electrophoresis (PFGE), multilocus sequence typing (MLST), multilocus variable-number tandem-repeat analysis (MLVA), binary typing, and serotyping. Initially 423 L. monocytogenes isolates underwent WGS, and comparisons uncovered a diverse genetic population structure derived from three distinct lineages. MLST, binary typing, and serotyping results inferred in silico from the WGS data were highly concordant (>99%) with laboratory typing performed in parallel. However, WGS was able to identify distinct nested clusters within groups of isolates that were otherwise indistinguishable using our current typing methods. Routine WGS was then used for prospective epidemiologic surveillance on a further 97 L. monocytogenes isolates over a 12-month period, which provided a greater level of discrimination than that of conventional typing for inferring linkage to point source outbreaks. A risk-based alert system based on WGS similarity was used to inform epidemiologists required to act on the data. Our experience shows that WGS can be adopted for prospective L. monocytogenes surveillance and investigated for other pathogens relevant to public health.     

Rapid Whole-Genome Sequencing of Mycobacterium tuberculosis Isolates Directly from Clinical Samples

The rapid identification of antimicrobial resistance is essential for effective treatment of highly resistant Mycobacterium tuberculosis. Whole-genome sequencing provides comprehensive data on resistance mutations and strain typing for monitoring transmission, but unlike for conventional molecular tests, this has previously been achievable only from cultures of M. tuberculosis. Here we describe a method utilizing biotinylated RNA baits designed specifically for M. tuberculosis DNA to capture full M. tuberculosis genomes directly from infected sputum samples, allowing whole-genome sequencing without the requirement of culture. This was carried out on 24 smear-positive sputum samples, collected from the United Kingdom and Lithuania where a matched culture sample was available, and 2 samples that had failed to grow in culture. M. tuberculosis sequencing data were obtained directly from all 24 smear-positive culture-positive sputa, of which 20 were of high quality (>20× depth and >90% of the genome covered). Results were compared with those of conventional molecular and culture-based methods, and high levels of concordance between phenotypical resistance and predicted resistance based on genotype were observed. High-quality sequence data were obtained from one smear-positive culture-negative case. This study demonstrated for the first time the successful and accurate sequencing of M. tuberculosis genomes directly from uncultured sputa. Identification of known resistance mutations within a week of sample receipt offers the prospect for personalized rather than empirical treatment of drug-resistant tuberculosis, including the use of antimicrobial-sparing regimens, leading to improved outcomes.     

Isolation and Enrichment of Cryptosporidium DNA and Verification of DNA Purity for Whole-Genome Sequencing

Whole-genome sequencing of Cryptosporidium spp. is hampered by difficulties in obtaining sufficient, highly pure genomic DNA from clinical specimens. In this study, we developed procedures for the isolation and enrichment of Cryptosporidium genomic DNA from fecal specimens and verification of DNA purity for whole-genome sequencing. The isolation and enrichment of genomic DNA were achieved by a combination of three oocyst purification steps and whole-genome amplification (WGA) of DNA from purified oocysts. Quantitative PCR (qPCR) analysis of WGA products was used as an initial quality assessment of amplified genomic DNA. The purity of WGA products was assessed by Sanger sequencing of cloned products. Next-generation sequencing tools were used in final evaluations of genome coverage and of the extent of contamination. Altogether, 24 fecal specimens of Cryptosporidium parvum, C. hominis, C. andersoni, C. ubiquitum, C. tyzzeri, and Cryptosporidium chipmunk genotype I were processed with the procedures. As expected, WGA products with low (<16.0) threshold cycle (C_T) values yielded mostly Cryptosporidium sequences in Sanger sequencing. The cloning-sequencing analysis, however, showed significant contamination in 5 WGA products (proportion of positive colonies derived from Cryptosporidium genomic DNA, ≤25%). Following this strategy, 20 WGA products from six Cryptosporidium species or genotypes with low (mostly <14.0) C_T values were submitted to whole-genome sequencing, generating sequence data covering 94.5% to 99.7% of Cryptosporidium genomes, with mostly minor contamination from bacterial, fungal, and host DNA. These results suggest that the described strategy can be used effectively for the isolation and enrichment of Cryptosporidium DNA from fecal specimens for whole-genome sequencing.     

Rapid Whole-Genome Sequencing for Detection and Characterization of Microorganisms Directly from Clinical Samples

Whole-genome sequencing (WGS) is becoming available as a routine tool for clinical microbiology. If applied directly on clinical samples, this could further reduce diagnostic times and thereby improve control and treatment. A major bottleneck is the availability of fast and reliable bioinformatic tools. This study was conducted to evaluate the applicability of WGS directly on clinical samples and to develop easy-to-use bioinformatic tools for the analysis of sequencing data. Thirty-five random urine samples from patients with suspected urinary tract infections were examined using conventional microbiology, WGS of isolated bacteria, and direct sequencing on pellets from the urine samples. A rapid method for analyzing the sequence data was developed. Bacteria were cultivated from 19 samples but in pure cultures from only 17 samples. WGS improved the identification of the cultivated bacteria, and almost complete agreement was observed between phenotypic and predicted antimicrobial susceptibilities. Complete agreement was observed between species identification, multilocus sequence typing, and phylogenetic relationships for Escherichia coli and Enterococcus faecalis isolates when the results of WGS of cultured isolates and urine samples were directly compared. Sequencing directly from the urine enabled bacterial identification in polymicrobial samples. Additional putative pathogenic strains were observed in some culture-negative samples. WGS directly on clinical samples can provide clinically relevant information and drastically reduce diagnostic times. This may prove very useful, but the need for data analysis is still a hurdle to clinical implementation. To overcome this problem, a publicly available bioinformatic tool was developed in this study.     

Clinical Application of Whole-Genome Sequencing To Inform Treatment for Multidrug-Resistant Tuberculosis Cases

The treatment of drug-resistant tuberculosis cases is challenging, as drug options are limited, and the existing diagnostics are inadequate. Whole-genome sequencing (WGS) has been used in a clinical setting to investigate six cases of suspected extensively drug-resistant Mycobacterium tuberculosis (XDR-TB) encountered at a London teaching hospital between 2008 and 2014. Sixteen isolates from six suspected XDR-TB cases were sequenced; five cases were analyzed in a clinically relevant time frame, with one case sequenced retrospectively. WGS identified mutations in the M. tuberculosis genes associated with antibiotic resistance that are likely to be responsible for the phenotypic resistance. Thus, an evidence base was developed to inform the clinical decisions made around antibiotic treatment over prolonged periods. All strains in this study belonged to the East Asian (Beijing) lineage, and the strain relatedness was consistent with the expectations from the case histories, confirming one contact transmission event. We demonstrate that WGS data can be produced in a clinically relevant time scale some weeks before drug sensitivity testing (DST) data are available, and they actively help clinical decision-making through the assessment of whether an isolate (i) has a particular resistance mutation where there are absent or contradictory DST results, (ii) has no further resistance markers and therefore is unlikely to be XDR, or (iii) is identical to an isolate of known resistance (i.e., a likely transmission event). A small number of discrepancies between the genotypic predictions and phenotypic DST results are discussed in the wider context of the interpretation and reporting of WGS results.     

Detection and Whole-Genome Sequencing of Carbapenemase-Producing Aeromonas hydrophila Isolates from Routine Perirectal Surveillance Culture

Perirectal surveillance cultures and a stool culture grew Aeromonas species from three patients over a 6-week period and were without epidemiological links. Detection of the bla_KPC-2 gene in one isolate prompted inclusion of non-Enterobacteriaceae in our surveillance culture workup. Whole-genome sequencing confirmed that the isolates were unrelated and provided data for Aeromonas reference genomes.     

Integration of Whole-Genome Sequencing into Infection Control Practices: the Potential and the Hurdles

Microbial whole-genome sequencing (WGS) is poised to transform many of the currently used approaches in medical microbiology. Recent reports on the application of WGS to understand genetic evolution and reconstruct transmission pathways have provided valuable information that will influence infection control practices. While this technology holds great promise, obstacles to full implementation remain. Two articles in this issue of the Journal of Clinical Microbiology (S. Octavia, Q. Wang, M. M. Tanaka, S. Kaur, V. Sintchenko, and R. Lan, J Clin Microbiol 53:1063–1071, 2015, doi:10.1128/JCM.03235-14, and S. J. Salipante, D. J. SenGupta, L. A. Cummings, T. A. Land, D. R. Hoogestraat, and B. T. Cookson, J Clin Microbiol 53:1072–1079, 2015, doi:10.1128/JCM.03385-14) describe the breadth of application of WGS to the field of clinical epidemiology.     

Comparison of Multilocus Variable-Number Tandem-Repeat Analysis and Whole-Genome Sequencing for Investigation of Clostridium difficile Transmission

No study to date has compared multilocus variable-number tandem-repeat analysis (MLVA) and whole-genome sequencing (WGS) in an investigation of the transmission of Clostridium difficile infection. Isolates from 61 adults with ongoing and/or recurrent C. difficile infections and 17 asymptomatic carriage episodes in children (201 samples), as well as from 61 suspected outbreaks affecting 2 to 41 patients in 31 hospitals in the United Kingdom (300 samples), underwent 7-locus MLVA and WGS in parallel. When the first and last samples from the same individual taken for a median (interquartile range [IQR]) of 63 days (43 to 105 days) apart were compared, the estimated rates of the evolution of single nucleotide variants (SNVs), summed tandem-repeat differences (STRDs), and locus variants (LVs) were 0.79 (95% confidence interval [CI], 0.00 to 1.75), 1.63 (95% CI, 0.00 to 3.59), and 1.21 (95% CI, 0.00 to 2.67)/called genome/year, respectively. Differences of >2 SNVs and >10 STRDs have been used to exclude direct case-to-case transmission. With the first serial sample per individual being used to assess discriminatory power, across all pairs of samples sharing a PCR ribotype, 192/283 (68%) differed by >10 STRDs and 217/283 (77%) by >2 SNVs. Among all pairs of cases from the same suspected outbreak, 1,190/1,488 (80%) pairs had concordant results using >2 SNVs and >10 STRDs to exclude transmission. For the discordant pairs, 229 (15%) had ≥2 SNVs but ≤10 STRDs, and 69 (5%) had ≤2 SNVs but ≥10 STRDs. Discordant pairs had higher numbers of LVs than concordant pairs, supporting the more diverse measure in each type of discordant pair. Conclusions on whether the potential outbreaks were confirmed were concordant in 58/61 (95%) investigations. Overall findings using MLVA and WGS were very similar despite the fact that they analyzed different parts of the bacterial genome. With improvements in WGS technology, it is likely that MLVA locus data will be available from WGS in the near future.