Delineating community outbreaks of Salmonella enterica serovar Typhimurium by use of whole-genome sequencing: insights into genomic variability within an outbreak

Epidemiology and whole genome sequencing of an ongoing point-source Salmonella Agona outbreak associated with sushi consumption in western Sydney, Australia 2015

During May 2015, an increase in Salmonella Agona cases was reported from western Sydney, Australia. We examine the public health actions used to investigate and control this increase. A descriptive case-series investigation was conducted. Six outbreak cases were identified; all had consumed cooked tuna sushi rolls purchased within a western Sydney shopping complex. Onset of illness for outbreak cases occurred between 7 April and 24 May 2015. Salmonella was isolated from food samples collected from the implicated premise and a prohibition order issued. No further cases were identified following this action. Whole genome sequence (WGS) analysis was performed on isolates recovered during this investigation, with additional S. Agona isolates from sporadic-clinical cases and routine food sampling in New South Wales, January to July 2015. Clinical isolates of outbreak cases were indistinguishable from food isolates collected from the implicated sushi outlet. Five additional clinical isolates not originally considered to be linked to the outbreak were genomically similar to outbreak isolates, indicating the point-source contamination may have started before routine surveillance identified an increase. This investigation demonstrated the value of genomics-guided public health action, where near real-time WGS enhanced the resolution of the epidemiological investigation.

Genomic heterogeneity of Salmonella enterica serovar Typhimurium bacteriuria from chronic infection

We sequenced the genomes of 14 sequential Salmonella enterica serovar Typhimurium isolates obtained over a five year period from a patient with persistent Salmonella bacteriuria. The isolates formed five distinct lineages; two of which co-existed over four years. We inferred that the observed within-patient variation resulted from mutation events.

Impact of whole genome sequencing in Public Health reference laboratories

Public Health Microbiology reference laboratories fulfil a critical role in providing overarching testing and surveillance for notifiable, emerging and important pathogens. These duties require the laboratory to possess an extensive repertoire of validated assays and the ability to rapidly respond to novel threats and outbreaks. For these, among other reasons, the ‘one stop shop' approach of whole genome sequencing (WGS) has been embraced by microbiology reference laboratories. The ability to replace multiple labour-intensive assays with a single technique of superior typeability and discrimination at an often competitive price, although not without its challenges, has already begun to change the workflow of Public Health reference laboratories.

Evolution of Variable Number Tandem Repeats and Its Relationship with Genomic Diversity in Salmonella Typhimurium

Salmonella enterica serovar Typhimurium is the most common Salmonella serovar causing human infections in Australia and many other countries. A total of 12,112 S. Typhimurium isolates from New South Wales were analyzed by multi-locus variable number of tandem repeat (VNTR) analysis (MLVA) using five VNTRs from 2007 to 2014. We found that mid ranges of repeat units of 8–14 in VNTR locus STTR5, 6–13 in STTR6, and 9–12 in STTR10 were always predominant in the population (>50%). In vitro passaging experiments using MLVA type carrying extreme length alleles found that the majority of long length alleles mutated to short ones and short length alleles mutated to longer ones. Both data suggest directional mutability of VNTRs toward mid-range repeats. Sequencing of 28 isolates from a newly emerged MLVA type and its five single locus variants revealed that single nucleotide variation between isolates with up to two MLVA differences ranged from 0 to 12 single nucleotide polymorphisms (SNPs). However, there was no relationship between SNP and VNTR differences. A population genetic model of the joint distribution of VNTRs and SNPs variations was used to estimate the mutation rates of the two markers, yielding a ratio of 1 VNTR change to 6.9 SNP changes. When only one VNTR repeat difference was considered, the majority of pairwise SNP difference between isolates were 4 SNPs or fewer. Based on this observation and our previous findings of SNP differences of outbreak isolates, we suggest that investigation of S. Typhimurium community outbreaks should include cases of 1 repeat difference to increase sensitivity. This study offers new insights into the short-term VNTR evolution of S. Typhimurium and its application for epidemiological typing.

Navigating Microbiological Food Safety in the Era of Whole-Genome Sequencing

The epidemiological investigation of a foodborne outbreak, including identification of related cases, source attribution, and development of intervention strategies, relies heavily on the ability to subtype the etiological agent at a high enough resolution to differentiate related from nonrelated cases. Historically, several different molecular subtyping methods have been used for this purpose; however, emerging techniques, such as single nucleotide polymorphism (SNP)-based techniques, that use whole-genome sequencing (WGS) offer a resolution that was previously not possible. With WGS, unlike traditional subtyping methods that lack complete information, data can be used to elucidate phylogenetic relationships and disease-causing lineages can be tracked and monitored over time. The subtyping resolution and evolutionary context provided by WGS data allow investigators to connect related illnesses that would be missed by traditional techniques. The added advantage of data generated by WGS is that these data can also be used for secondary analyses, such as virulence gene detection, antibiotic resistance gene profiling, synteny comparisons, mobile genetic element identification, and geographic attribution. In addition, several software packages are now available to generate in silico results for traditional molecular subtyping methods from the whole-genome sequence, allowing for efficient comparison with historical databases. Metagenomic approaches using next-generation sequencing have also been successful in the detection of nonculturable foodborne pathogens. This review addresses state-of-the-art techniques in microbial WGS and analysis and then discusses how this technology can be used to help support food safety investigations. Retrospective outbreak investigations using WGS are presented to provide organism-specific examples of the benefits, and challenges, associated with WGS in comparison to traditional molecular subtyping techniques.

Whole genome sequencing of Salmonella Typhimurium illuminates distinct outbreaks caused by an endemic multi-locus variable number tandem repeat analysis type in Australia, 2014

Background

Salmonella Typhimurium (STM) is an important cause of foodborne outbreaks worldwide. Subtyping of STM remains critical to outbreak investigation, yet current techniques (e.g. multilocus variable number tandem repeat analysis, MLVA) may provide insufficient discrimination. Whole genome sequencing (WGS) offers potentially greater discriminatory power to support infectious disease surveillance.

Methods

We performed WGS on 62 STM isolates of a single, endemic MLVA type associated with two epidemiologically independent, food-borne outbreaks along with sporadic cases in New South Wales, Australia, during 2014. Genomes of case and environmental isolates were sequenced using HiSeq (Illumina) and the genetic distance between them was assessed by single nucleotide polymorphism (SNP) analysis. SNP analysis was compared to the epidemiological context.

Results

The WGS analysis supported epidemiological evidence and genomes of within-outbreak isolates were nearly identical. Sporadic cases differed from outbreak cases by a small number of SNPs, although their close relationship to outbreak cases may represent an unidentified common food source that may warrant further public health follow up. Previously unrecognised mini-clusters were detected.

Conclusions

WGS of STM can discriminate foodborne community outbreaks within a single endemic MLVA clone. Our findings support the translation of WGS into public health laboratory surveillance of salmonellosis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-016-0831-3) contains supplementary material, which is available to authorized users.

Genomic Epidemiology: Whole-Genome-Sequencing-Powered Surveillance and Outbreak Investigation of Foodborne Bacterial Pathogens

As we are approaching the twentieth anniversary of PulseNet, a network of public health and regulatory laboratories that has changed the landscape of foodborne illness surveillance through molecular subtyping, public health microbiology is undergoing another transformation brought about by so-called next-generation sequencing (NGS) technologies that have made whole-genome sequencing (WGS) of foodborne bacterial pathogens a realistic and superior alternative to traditional subtyping methods. Routine, real-time, and widespread application of WGS in food safety and public health is on the horizon. Technological, operational, and policy challenges are still present and being addressed by an international and multidisciplinary community of researchers, public health practitioners, and other stakeholders.

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.

Usefulness of High-Quality Core Genome Single-Nucleotide Variant Analysis for Subtyping the Highly Clonal and the Most Prevalent Salmonella enterica Serovar Heidelberg Clone in the Context of Outbreak Investigations

Salmonella enterica serovar Heidelberg is the second most frequently occurring serovar in Quebec and the third-most prevalent in Canada. Given that conventional pulsed-field gel electrophoresis (PFGE) subtyping for common Salmonella serovars, such as S. Heidelberg, yields identical subtypes for the majority of isolates recovered, public health laboratories are desperate for new subtyping tools to resolve highly clonal S. Heidelberg strains involved in outbreak events. As PFGE was unable to discriminate isolates from three epidemiologically distinct outbreaks in Quebec, this study was conducted to evaluate whole-genome sequencing (WGS) and phylogenetic analysis as an alternative to conventional subtyping tools. Genomes of 46 isolates from 3 Quebec outbreaks (2012, 2013, and 2014) supported by strong epidemiological evidence were sequenced and analyzed using a high-quality core genome single-nucleotide variant (hqSNV) bioinformatics approach (SNV phylogenomics [SNVphyl] pipeline). Outbreaks were indistinguishable by conventional PFGE subtyping, exhibiting the same PFGE pattern (SHEXAI.0001/SHEBNI.0001). Phylogenetic analysis based on hqSNVs extracted from WGS separated the outbreak isolates into three distinct groups, 100% concordant with the epidemiological data. The minimum and maximum number of hqSNVs between isolates from the same outbreak was 0 and 4, respectively, while >59 hqSNVs were measured between 2 previously indistinguishable outbreaks having the same PFGE and phage type, thus corroborating their distinction as separate unrelated outbreaks. This study demonstrates that despite the previously reported high clonality of this serovar, the WGS-based hqSNV approach is a superior typing method, capable of resolving events that were previously indistinguishable using classic subtyping tools.

Laboratory-Acquired Infection with Salmonella enterica Serovar Typhimurium Exposed by Whole-Genome Sequencing

Despite advances in laboratory design, professional training, and workplace biosafety guidelines, laboratory-acquired infections continue to occur. Effective tools are required to investigate cases and prevent future illness. Here, we demonstrate the value of whole-genome sequencing as a tool for the identification and source attribution of laboratory-acquired salmonellosis.

It Is Not All about Single Nucleotide Polymorphisms: Comparison of Mobile Genetic Elements and Deletions in Listeria monocytogenes Genomes Links Cases of Hospital-Acquired Listeriosis to the Environmental Source

The control of food-borne outbreaks caused by Listeria monocytogenes in humans relies on the timely identification of food or environmental sources and the differentiation of outbreak-related isolates from unrelated ones. This study illustrates the utility of whole-genome sequencing for examining the link between clinical and environmental isolates of L. monocytogenes associated with an outbreak of hospital-acquired listeriosis in Sydney, Australia. Comparative genomic analysis confirmed an epidemiological link between the three clinical and two environmental isolates. Single nucleotide polymorphism (SNP) analysis showed that only two SNPs separated the three human outbreak isolates, which differed by 19 to 20 SNPs from the environmental isolates and 71 to >10,000 SNPs from sporadic L. monocytogenes isolates. The chromosomes of all human outbreak isolates and the two suspected environmental isolates were syntenic. In contrast to the genomes of background sporadic isolates, all epidemiologically linked isolates contained two novel prophages and a previously unreported clustered regularly interspaced short palindromic repeat (CRISPR)-associated (Cas) locus subtype sequence. The mobile genetic element (MGE) profile of these isolates was distinct from that of the other serotype 1/2b reference strains and sporadic isolates. The identification of SNPs and clonally distinctive MGEs strengthened evidence to distinguish outbreak-related isolates of L. monocytogenes from cocirculating endemic strains.

Genomic Variability of Serial Human Isolates of Salmonella enterica Serovar Typhimurium Associated with Prolonged Carriage

Salmonella enterica serovar Typhimurium is an important foodborne human pathogen that often causes self-limiting but severe gastroenteritis. Prolonged excretion of S. Typhimurium after the infection can lead to secondary transmissions. However, little is known about within-host genomic variation in bacteria associated with asymptomatic shedding. Genomes of 35 longitudinal isolates of S. Typhimurium recovered from 11 patients (children and adults) with culture-confirmed gastroenteritis were sequenced. There were three or four isolates obtained from each patient. Single nucleotide polymorphisms (SNPs) were analyzed in these isolates, which were recovered between 1 and 279 days after the initial diagnosis. Limited genomic variation (5 SNPs or fewer) was associated with short- and long-term carriage of S. Typhimurium. None of the isolates was shown to be due to reinfection. SNPs occurred randomly, and the majority of the SNPs were nonsynonymous. Two nonsense mutations were observed. A nonsense mutation in flhC rendered the isolate nonmotile, whereas the significance of a nonsense mutation in yihV is unknown. The estimated mutation rate is 1.49 × 10(-6) substitution per site per year. S. Typhimurium isolates excreted in stools following acute gastroenteritis in children and adults demonstrated limited genomic variability over time, regardless of the duration of carriage. These findings have important implications for the detection of possible transmission events suspected by public health genomic surveillance of S. Typhimurium infections.

Whole-genome sequencing targets drug-resistant bacterial infections

During the past two decades, the technological progress of whole-genome sequencing (WGS) had changed the fields of Environmental Microbiology and Biotechnology, and, currently, is changing the underlying principles, approaches, and fundamentals of Public Health, Epidemiology, Health Economics, and national productivity. Today’s WGS technologies are able to compete with conventional techniques in cost, speed, accuracy, and resolution for day-to-day control of infectious diseases and outbreaks in clinical laboratories and in long-term epidemiological investigations. WGS gives rise to an exciting future direction for personalized Genomic Epidemiology. One of the most vital and growing public health problems is the emerging and re-emerging of multidrug-resistant (MDR) bacterial infections in the communities and healthcare settings, reinforced by a decline in antimicrobial drug discovery. In recent years, retrospective analysis provided by WGS has had a great impact on the identification and tracking of MDR microorganisms in hospitals and communities. The obtained genomic data are also important for developing novel easy-to-use diagnostic assays for clinics, as well as for antibiotic and therapeutic development at both the personal and population levels. At present, this technology has been successfully applied as an addendum to the real-time diagnostic methods currently used in clinical laboratories. However, the significance of WGS for public health may increase if: (a) unified and user-friendly bioinformatics toolsets for easy data interpretation and management are established, and (b) standards for data validation and verification are developed. Herein, we review the current and future impact of this technology on diagnosis, prevention, treatment, and control of MDR infectious bacteria in clinics and on the global scale.

A multiplex oligonucleotide ligation-PCR as a complementary tool for subtyping of Salmonella Typhimurium

Subtyping below the serovar level is essential for surveillance and outbreak detection and investigation of Salmonella enterica subsp. enterica serovar Typhimurium (S. Typhimurium) and its monophasic variant 1,4,[5],12:i:- (S. 1,4,[5],12:i:-), frequent causes of foodborne infections. In an attempt to overcome the intrinsic shortcomings of currently used subtyping techniques, a multiplex oligonucleotide ligation-PCR (MOL-PCR) assay was developed which combines different types of molecular markers in a high-throughput microsphere suspension array. The 52 molecular markers include prophage genes, amplified fragment length polymorphism (AFLP) elements, Salmonella genomic island 1 (SGI1), allantoinase gene allB, MLVA locus STTR10, antibiotic resistance genes, single nucleotide polymorphisms (SNPs) and phase 2 flagellar gene fljB. The in vitro stability of these markers was confirmed in a serial passage experiment. The validation of the MOL-PCR assay for subtyping of S. Typhimurium and S. 1,4,[5],12:i:- on 519 isolates shows that the method is rapid, reproducible, flexible, accessible, easy to use and relatively inexpensive. Additionally, a 100 % typeability and a discriminatory power equivalent to that of phage typing were observed, and epidemiological concordance was assessed on isolates of 2 different outbreaks. Furthermore, a data analysis method is provided so that the MOL-PCR assay allows for objective, computerised data analysis and data interpretation of which the results can be easily exchanged between different laboratories in an international surveillance network.

Use of Alignment-Free Phylogenetics for Rapid Genome Sequence-Based Typing of Helicobacter pylori Virulence Markers and Antibiotic Susceptibility

Whole-genome sequencing is becoming a leading technology in the typing and epidemiology of microbial pathogens, but the increase in genomic information necessitates significant investment in bioinformatic resources and expertise, and currently used methodologies struggle with genetically heterogeneous bacteria such as the human gastric pathogen Helicobacter pylori. Here we demonstrate that the alignment-free analysis method feature frequency profiling (FFP) can be used to rapidly construct phylogenetic trees of draft bacterial genome sequences on a standard desktop computer and that coupling with in silico genotyping methods gives useful information for comparative and clinical genomic and molecular epidemiology applications. FFP-based phylogenetic trees of seven gastric Helicobacter species matched those obtained by analysis of 16S rRNA genes and ribosomal proteins, and FFP- and core genome single nucleotide polymorphism-based analysis of 63 H. pylori genomes again showed comparable phylogenetic clustering, consistent with genomotypes assigned by using multilocus sequence typing (MLST). Analysis of 377 H. pylori genomes highlighted the conservation of genomotypes and linkage with phylogeographic characteristics and predicted the presence of an incomplete or nonfunctional cag pathogenicity island in 18/276 genomes. In silico analysis of antibiotic susceptibility markers suggests that most H. pylori hspAmerind and hspEAsia isolates are predicted to carry the T₂₈₁₂C mutation potentially conferring low-level clarithromycin resistance, while levels of metronidazole resistance were similar in all multilocus sequence types. In conclusion, the use of FFP phylogenetic clustering and in silico genotyping allows determination of genome evolution and phylogeographic clustering and can contribute to clinical microbiology by genomotyping for outbreak management and the prediction of pathogenic potential and antibiotic susceptibility.

Defining the Core Genome of Salmonella enterica Serovar Typhimurium for Genomic Surveillance and Epidemiological Typing

Salmonella enterica serovar Typhimurium is the most common Salmonella serovar causing foodborne infections in Australia and many other countries. Twenty-one S. Typhimurium strains from Salmonella reference collection A (SARA) were analyzed using Illumina high-throughput genome sequencing. Single nucleotide polymorphisms (SNPs) in 21 SARA strains ranged from 46 to 11,916 SNPs, with an average of 1,577 SNPs per strain. Together with 47 strains selected from publicly available S. Typhimurium genomes, the S. Typhimurium core genes (STCG) were determined. The STCG consist of 3,846 genes, a set that is much larger than that of the 2,882 Salmonella core genes (SCG) found previously. The STCG together with 1,576 core intergenic regions (IGRs) were defined as the S. Typhimurium core genome. Using 93 S. Typhimurium genomes from 13 epidemiologically confirmed community outbreaks, we demonstrated that typing based on the S. Typhimurium core genome (STCG plus core IGRs) provides superior resolution and higher discriminatory power than that based on SCG for outbreak investigation and molecular epidemiology of S. Typhimurium. STCG and STCG plus core IGR typing achieved 100% separation of all outbreaks compared to that of SCG typing, which failed to separate isolates from two outbreaks from background isolates. Defining the S. Typhimurium core genome allows standardization of genes/regions to be used for high-resolution epidemiological typing and genomic surveillance of S. Typhimurium.

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, andS. 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.