Skip to main content
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Current Issue
    • Accepted Manuscripts
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About JCM
    • Editor in Chief
    • Editorial Board
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • Subscribe
    • Members
    • Institutions
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
Journal of Clinical Microbiology
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Current Issue
    • Accepted Manuscripts
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About JCM
    • Editor in Chief
    • Editorial Board
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • Subscribe
    • Members
    • Institutions
Bacteriology

Defining and Evaluating a Core Genome Multilocus Sequence Typing Scheme for Whole-Genome Sequence-Based Typing of Listeria monocytogenes

Werner Ruppitsch, Ariane Pietzka, Karola Prior, Stefan Bletz, Haizpea Lasa Fernandez, Franz Allerberger, Dag Harmsen, Alexander Mellmann
D. J. Diekema, Editor
Werner Ruppitsch
aGerman-Austrian Binational Consiliary Laboratory for Listeria, Austrian Agency for Health and Food Safety (AGES), Vienna, Austria
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ariane Pietzka
aGerman-Austrian Binational Consiliary Laboratory for Listeria, Austrian Agency for Health and Food Safety (AGES), Vienna, Austria
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Karola Prior
bDepartment for Periodontology, University Hospital Muenster, Muenster, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Stefan Bletz
cInstitute of Hygiene, University Hospital Muenster, Muenster, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Haizpea Lasa Fernandez
aGerman-Austrian Binational Consiliary Laboratory for Listeria, Austrian Agency for Health and Food Safety (AGES), Vienna, Austria
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Franz Allerberger
aGerman-Austrian Binational Consiliary Laboratory for Listeria, Austrian Agency for Health and Food Safety (AGES), Vienna, Austria
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Dag Harmsen
bDepartment for Periodontology, University Hospital Muenster, Muenster, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alexander Mellmann
cInstitute of Hygiene, University Hospital Muenster, Muenster, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
D. J. Diekema
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/JCM.01193-15
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Whole-genome sequencing (WGS) has emerged today as an ultimate typing tool to characterize Listeria monocytogenes outbreaks. However, data analysis and interlaboratory comparability of WGS data are still challenging for most public health laboratories. Therefore, we have developed and evaluated a new L. monocytogenes typing scheme based on genome-wide gene-by-gene comparisons (core genome multilocus the sequence typing [cgMLST]) to allow for a unique typing nomenclature. Initially, we determined the breadth of the L. monocytogenes population based on MLST data with a Bayesian approach. Based on the genome sequence data of representative isolates for the whole population, cgMLST target genes were defined and reappraised with 67 L. monocytogenes isolates from two outbreaks and serotype reference strains. The Bayesian population analysis generated five L. monocytogenes groups. Using all available NCBI RefSeq genomes (n = 36) and six additionally sequenced strains, all genetic groups were covered. Pairwise comparisons of these 42 genome sequences resulted in 1,701 cgMLST targets present in all 42 genomes with 100% overlap and ≥90% sequence similarity. Overall, ≥99.1% of the cgMLST targets were present in 67 outbreak and serotype reference strains, underlining the representativeness of the cgMLST scheme. Moreover, cgMLST enabled clustering of outbreak isolates with ≤10 alleles difference and unambiguous separation from unrelated outgroup isolates. In conclusion, the novel cgMLST scheme not only improves outbreak investigations but also enables, due to the availability of the automatically curated cgMLST nomenclature, interlaboratory exchange of data that are crucial, especially for rapid responses during transsectorial outbreaks.

INTRODUCTION

Listeria monocytogenes is a facultative anaerobe, a Gram-positive, psychrophilic and salt-tolerant, facultative intracellular pathogen of humans and animals, causing clinical manifestations like gastroenteritis, encephalitis, meningitis, and septicemia. A high hospitalization rate of >90% and a case-fatality rate up to 30% make L. monocytogenes an important human pathogen (1). The characteristic traits (growth at low temperatures, survival of freezing and high-salt and nitrite preservation methods, and biofilm formation) of L. monocytogenes represent a major issue for industrialized food production and facilitate food contamination at several stages of food production (2). Nearly all cases of listeriosis are caused by consumption or use of contaminated food or feed.

The traditional L. monocytogenes serotyping scheme allows the differentiation of 12 serotypes of which 4b, 1/2a, and 1/2b isolates cause about 96% of all reported human listeriosis cases (3). Low discriminatory power, insufficient reproducibility, and antigen sharing between serotypes impede the value of serotyping in outbreak investigations and necessitate more accurate and more discriminatory typing solutions (4).

Pulsed-field gel electrophoresis (PFGE) has been established as the current “gold standard” for L. monocytogenes typing by PulseNet (5, 6) and has been essential for outbreak investigation worldwide (7). However, PFGE is time-consuming, expensive, and difficult to standardize (8, 9). Methods based on DNA sequence analysis appear more promising for fast, accurate, and reproducible strain typing (10). Whereas multilocus sequence typing (MLST) (11, 12) and multi-virulence-locus sequence typing (MVLST) (13) schemes for L. monocytogenes share the characteristics of sequence-based methods, they both lack the discriminative power needed for outbreak investigation of this clonal pathogen (7, 14).

Nowadays, the recent and ongoing evolution of sequencing technologies from Sanger sequencing to next-generation sequencing enables sequence analysis on a whole-genome level. Several studies on a variety of bacterial species have already shown that whole-genome sequence (WGS)-based typing, based either on single nucleotide variants (SNVs) (15, 16) or on gene-by-gene allelic profiling of core genome genes, frequently named core genome MLST (cgMLST) or MLST+ (17, 18), currently represents the ultimate diagnostic tool for strain typing. Recently, we successfully applied a cgMLST typing approach to L. monocytogenes (19). Nevertheless, the broad use of WGS-based approaches is still hampered by the lack of standardized nomenclature that would facilitate global exchange of data, as has already been the reality for classical MLST data (20) for more than a decade.

To achieve a stable cgMLST scheme for L. monocytogenes that can form the basis of a standardized nomenclature for WGS-based L. monocytogenes typing, first we defined an L. monocytogenes core genome gene set representing the genetic diversity within the L. monocytogenes population based on well-characterized reference strains, and second we challenged this scheme for suitability in outbreak investigations using isolates from two outbreaks and sporadic cases.

MATERIALS AND METHODS

Microorganisms and DNA extraction.All strains and genome sequences used for the development of the novel cgMLST L. monocytogenes scheme are listed in Table 1. For subsequent evaluation of the scheme, a total of 67 L. monocytogenes isolates from sporadic cases (n = 8 isolates, that served also as outgroups for the outbreaks with matching serotypes and highly similar or even identical PFGE pattern) and two outbreaks (n = 42) (21–23) with reference strains for all serotypes (n = 17) were used (Table 2). All strains were cultured overnight at 37°C on RAPID'L.Mono agar (Bio-Rad, Vienna, Austria) for species confirmation and subcultivated on Columbia blood agar plates (bioMérieux, Marcy I'Etoile, France) prior to DNA extraction using the GenElute Bacterial Genomic DNA kit (Sigma, St. Louis, MO, USA) according to the manufacturer's instructions.

View this table:
  • View inline
  • View popup
TABLE 1

List of L. monocytogenes strains and genomes used for SeqSphere cgMLST L. monocytogenes target definition

View this table:
  • View inline
  • View popup
TABLE 2

List of L. monocytogenes isolates used for evaluation of the SeqSphere cgMLST L. monocytogenes schemea

Whole-genome sequencing and assembly.Sequencing libraries were prepared using Nextera XT chemistry (Illumina Inc., San Diego, CA, USA) for a 250-bp paired-end sequencing run on an Illumina MiSeq sequencer. Samples were sequenced to aim for minimum coverage of 100-fold using Illumina's recommended standard protocols. The resulting FASTQ files were first quality trimmed and then de novo assembled using the Velvet assembler (24) integrated in Ridom SeqSphere+ software (25) (version 2.3; Ridom GmbH, Münster, Germany). Here, reads were trimmed at their 5′ and 3′ ends until an average base quality of 30 was reached in a window of 20 bases, and the assembly was performed with Velvet version 1.1.04 using optimized k-mer size and coverage cutoff values based on the average length of contigs with >1,000 bp.

BAPS.To determine the overall L. monocytogenes species variation, we applied a Bayesian analysis of population structure (BAPS) (26, 27). All multilocus sequence typing (MLST) data available as of 24 July 2014 (673 sequence types [STs]) were downloaded from the MLST website (14), and all allelic gene sequences per locus were multiple aligned using MUSCLE (28) and finally concatenated for each ST. The BAPS was carried out using the clustering of linked molecular data functionality. Ten runs were performed, setting an upper limit of 20 partitions. Admixture analysis was performed using the following parameters: minimum population size considered, 5; iterations, 50; number of reference individuals simulated from each population, 50; and number of iterations for each reference individual, 10.

cgMLST target gene definition.To determine the cgMLST gene set (named MLST+ in the SeqSphere+ software), a genome-wide gene-by-gene comparison was performed using the MLST+ target definer (version 1.1) function of SeqSphere+ with default parameters. These parameters comprise the following filters to exclude certain genes of the EGD-e reference genome (GenBank accession number NC_003210.1, dated 26 March 2015) from the cgMLST scheme: a minimum length filter that discards all genes shorter than 50 bp; a start codon filter that discards all genes that contain no start codon at the beginning of the gene; a stop codon filter that discards all genes that contain no stop codon or more than one stop codon or that do not have the stop codon at the end of the gene; a homologous gene filter that discards all genes with fragments that occur in multiple copies within a genome (with identity of 90% and >100 bp overlap); and a gene overlap filter that discards the shorter gene from the cgMLST scheme if the two genes affected overlap >4 bp. The remaining genes were then used in a pairwise comparison with BLAST version 2.2.12 (parameters used were word size 11, mismatch penalty −1, match reward 1, gap open costs 5, and gap extension costs 2) with the query L. monocytogenes chromosomes. All genes of the reference genome that were common in all query genomes with a sequence identity of ≥90% and 100% overlap and, with the default parameter stop codon percentage filter turned on, formed the final cgMLST scheme; this discards all genes that have internal stop codons in >20% of the query genomes.

Evaluation of the cgMLST target gene set.To evaluate the applicability and representativeness of the L. monocytogenes cgMLST target gene set, a total of 67 isolates (Table 2) were subsequently analyzed to determine the presence of these target genes. It was assumed that a well-defined cgMLST scheme should cover at least 95% of the cgMLST genes present in all isolates.

To extract the target genes, the default parameters were used in the SeqSphere+ software: (i) for processing options, “Ignore contigs shorter than 200 bases”; (ii) for scanning options, “Matching scanning thresholds for creating targets from assembled genomes” with “required identity to reference sequence of 90%” and “required alignment to reference sequence with 100%”; and (iii) for BLAST options, word size 11, mismatch penalty −1, match reward 1, gap open costs 5, and gap extension costs 2. In addition, the target genes were assessed for quality, i.e., the absence of frame shifts and ambiguous nucleotides. A core genome gene was considered a “good target” only if all of the above criteria were met, in which case the complete sequence was analyzed in comparison to the EGD-e reference. Alleles for each gene were assigned automatically by the SeqSphere+ software to ensure unique nomenclature. The combination of all alleles in each strain formed an allelic profile that was used to generate minimum spanning trees (MST) using the parameter “pairwise ignore missing values” during distance calculation.

In order to maintain backwards compatibility with classical L. monocytogenes MLST, sequences of the seven genes comprising the allelic profile of the MLST scheme were extracted separately from the genome sequences and queried against the L. monocytogenes MLST database in order to assign classical STs in silico.

Nucleotide sequence accession number.All raw reads generated were submitted to the European Nucleotide Archive (http://www.ebi.ac.uk/ena/) under the study accession number PRJEB6551.

RESULTS

BAPS partition and admixture analysis based on 673 STs resulted in seven partitions (see Table S1 in the supplemental material). As BAPS partitions Lm05 and Lm06 comprised exclusively Listeria innocua species isolates of 43 STs with significant admixtures, these two partitions were excluded from further analysis. For the remaining five partitions, three (partitions Lm01, Lm02, and Lm04) were among the available NCBI RefSeq genome sequences of L. monocytogenes. To achieve complete coverage of the L. monocytogenes population, we sequenced six additional strains from Seeliger's Listeria culture collection (SLCC0717, SLCC0759, SLCC1042, SLCC3287, SLCC4771, and SLCC6263), representing the missing BAPS partitions Lm03 and Lm07 (Table 1). In total, 42 genome sequences, including L. monocytogenes strain EDG-e as reference for core genome gene definition were fed into the MLST+ target definer and resulted in 1,701 genes out of 2,867 genes of strain EGD-e (53.2% of the EDG-e strain chromosome nucleotides) (see Table S2 in the supplemental material).The cgMLST scheme was then challenged with two sets of strains: the first contained 17 serotype reference strains representing all serotypes, genetic lineages, and BAPS partitions to determine its ability to cover the whole L. monocytogenes diversity and the second consisted of 48 isolates from two published outbreaks, including eight outgroup isolates (Table 2). All 17 serotype reference strains had ≥99.1% good cgMLST targets (mean, 99.5%), and for all serotype representatives the correct MLST was obtained. Similarly, for the two outbreaks, all isolates had ≥99.3% good cgMLST targets (mean, 99.7%). The results are summarized in Table 2.

The cgMLST scheme was further evaluated for its usability in outbreak investigation, i.e., whether outbreak isolates could be attributed to the same clone, named cluster type (CT) in the context of cgMLST typing, and clearly separated from the outgroup isolates. Therefore, we determined the maximum number of differing genes within each outbreak that reflect putative microevolutionary events. To facilitate cluster investigations in the future, we finally defined the so-called CT threshold that gives the maximum number of differing alleles that are shared by the same CT. In the two retrospectively analyzed outbreaks, a jellied pork outbreak (JPO) in Austria in the year 2008 and two epidemiologically linked clusters forming the acid curd cheese (Quargel) outbreak (ACCO) in Austria, the Czech Republic, and Germany in the years 2009/2010 (Table 2), detailed analysis resulted in a maximum number of 10 differing alleles (see Table S3 in the supplemental material). cgMLST of seven human and two food isolates from the JPO correctly grouped these isolates together with a maximum of four allelic differences (Fig. 1). Outgroup isolates L2708 (ST249) and L7508 (ST4) exhibited more than 1,000 allelic differences, and reference strains F2365 and LL195 (both ST1) exhibited ≥32 allelic differences (Fig. 1). Extraction of classical MLST targets resulted in STs of all outbreak isolates that were identical to those of ST1 and confirmed the previous Sanger sequencing (Table 2).

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Minimum-spanning tree based on cgMLST allelic profiles of 9 L. monocytogenes isolates (all share ST1) from the jellied pork outbreak (21) and two outgroup isolates L2708 (ST249) and L7508 (ST4) in comparison to reference strains F2365 (GenBank accession number NC_002973) and LL195 (NC_019556) (both ST1) exhibiting the same serotype 4b. Each circle represents an allelic profile based on sequence analysis of 1,701 cgMLST target genes. The numbers on the connecting lines illustrate the numbers of target genes with differing alleles. The different groups of strains are distinguished by the colors of the circles. Closely related genotypes (≤10 allele difference) are shaded in gray. NCBI RefSeq strains are marked with an asterisk.

cgMLST of 33 isolates from the ACCO correctly identified the two different clones (ACCO I and ACCO II) that caused this outbreak (Fig. 2). Within the ACCO I clone, nine isolates were ST403 and five were ST777, a single locus bglA variant of ST403. cgMLST revealed the same dichotomy as classical MLST; the right branch of the ACCO I tree comprised all ST777 isolates (L27-09, L31-09, L32-09, L35-09, and L68-09). All outgroup isolates (MRL-13-00230, LD27-12, and Ro-015) were ST403 with at least 16 allelic differences compared to the ACCO I isolates. ACCO I isolates displayed a maximum of 10 allelic differences from each other (Fig. 2). ACCO II isolates had a maximum of two allelic differences from each other. All ACCO II isolates were correctly assigned to ST398. The three epidemiologically unrelated outgroup isolates (L38-11, 12025641, and 12025647) with an identical PFGE band pattern (data not shown) also exhibited ST398 and had ≥23 allelic differences compared to the ACCO II food isolates (Fig. 2). ACCO I and ACCO II isolates differed in >1,000 alleles from each other (Fig. 2).

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Minimum-spanning tree illustrating the phylogenetic relationship based on the cgMLST allelic profiles of 33 L. monocytogenes isolates from the outbreak associated with acid curd cheese (ACCO) (22, 23) consisting of two clones (ACCO I and ACCO II). Three outgroup isolates per outbreak (with identical PFGE profiles and serotypes) are shown in comparison to the reference strain EGD-e (GenBank accession number NC_003210; ST35). ACCO I isolates L27-09, L31-09, L32-09, L35-09, and L68-09 were ST777; the remaining isolates, including the three ACCO I outgroup isolates were ST403. ACCO II isolates, including the three ACCO II outgroup isolates were all ST398. Each circle represents an allelic profile based on sequence analysis of 1,701 genes. The numbers on the connecting lines illustrate the numbers of target genes with differing alleles. The different groups of strains are distinguished by the colors of the circles. Closely related genotypes (≤10 allele difference) are shaded in gray. The NCBI RefSeq strain is marked with an asterisk.

DISCUSSION

In outbreak situations, a rapid, accurate, and standardized classification of bacterial isolates is essential. Since its introduction in 1998, MLST has become a proof-of-principle method for sequence-based typing methods with a unique centrally curated and thereby standardized nomenclature (20). Building on these experiences nowadays, it is possible to analyze thousands of genes using next-generation sequencing, which dramatically increases discriminatory power and thereby now enables outbreak investigations (18, 19, 29–32). In our study, we were able not only to show that our cgMLST typing scheme is representative for the breadth of the L. monocytogenes population with ≥99.1% successfully extracted cgMLST targets but also to differentiate outbreak from nonoutbreak isolates clearly.

The microevolutionary events within each outbreak and the CT threshold of ≤10 differences warrant further comments. Within the first outbreak, the JPO (21), very few allelic changes were detected and the maximum allelic distance within the outbreak was only four alleles. This high similarity reflects the outbreak situation without much time for intraoutbreak microevolution, because all patients belonged to one travel group and became ill after consuming contaminated jellied pork at an Austrian tavern (21). The ACCO cluster of listeriosis occurred from 2009 until 2010 in Austria, Germany, and the Czech Republic and was caused by contaminated acid curd cheese (Quargel) (22, 23). Further epidemiological and molecular outbreak investigations revealed that two different serotype 1/2a clones with distinct PFGE patterns and inlB STs were responsible for this outbreak (33). Interestingly, a recent study focusing on the comparative genomics of the two outbreak clones revealed significant differences in virulence (34). Again, cgMLST analysis corroborated these findings, and the number of differing alleles among the outbreak clones again reflected the outbreak length. Whereas the ACCO I isolates were found over a period of 8 months and up to 10 different alleles were detected; isolates of the ACCO II were found only during a 3-month period, where at maximum two different alleles were recorded. Therefore, we assume that ACCO I isolates are a representative microevolutionary model for the CT threshold determination to facilitate outbreak investigations using cgMLST. Although the software supports outbreak investigation by providing the CT, this does not release the epidemiologist from thorough investigation.

MLST and cgMLST both use alleles and not nucleotide polymorphisms as units of comparison. Irrespective of the number of nucleotide polymorphisms involved, each allelic change is numbered as a single event; i.e., an allelic change is related to at least one point mutation but can also contain several nucleotide changes. This principle covers the conflicting signals of horizontal and vertical transfer of genetic material and considers the higher frequency of recombination than point mutations in bacteria (30, 32). One major advantage of such an allele-based approach is easy storage and curating the nomenclature in a central database, which is obligatory to guarantee universal nomenclature. For classical MLST, this scenario was one of the key factors to success. However, manual curation of the current MLST databases frequently hampers the rapid use of novel allelic sequences as human intervention is necessary to assign new alleles and STs. With the software solution used here, it was already possible to automatically assign novel cgMLST alleles, after dedicated quality control of the read and assembly data. This automation is crucial as the vast amount of sequencing data is not humanly readable anymore in a reasonable time frame that is needed for effective implementation of hygiene measures during outbreaks. The immediate and automated assignment of novel alleles also enables any software user to access identical nomenclature for L. monocytogenes cgMLST typing, a prerequisite for successful interlaboratory exchange of data. In the future, it is desirable to have an open Internet-based nomenclature server that is able to be interrogated by any software or user (35). The SpaServer (http://spaserver.ridom.de), which automatically hosts the nomenclature of the Staphylococcus aureus protein A gene typing (spa) and now contains >300,000 typing entries originating from >100 countries, might serve as a blueprint for such service (36).

Our approach has one limitation. The analysis is reduced to coding regions only because the second-generation sequencing instruments currently in use produce only relatively short reads that do not assemble the frequently highly repetitive intergenic regions well, leading to faulty assemblies. Therefore, when second-generation sequencing machines are used, focusing on coding regions helps to improve the analytical quality. This might change when third-generation sequencing instruments that produce much longer reads from a single molecule are widely available, preferably as benchtop systems. Nevertheless, the current cgMLST approach will be sustainable as it will maintain backward compatibility with expansion of typing schemes to present typing as we see today with the in silico extraction of classical MLST STs from WGS data.

In conclusion, we established a highly representative cgMLST scheme for WGS-based typing of L. monocytogenes and demonstrated both a high discriminatory power and concordance to previous findings in different outbreak scenarios. The remaining challenge is to establish an Internet-based nomenclature server that can be interrogated like the current MLST servers to facilitate universal global nomenclature for any user.

ACKNOWLEDGMENTS

This work was funded by the European Community's Seventh Framework Program (grant FP7/2007-2013 to D.H.) under grant agreement 278864 in the framework of the EU PathoNGenTrace project, by the Medical Faculty of the University of Münster (grant BD9817044 to AM), and by a Leonardo da Vinci Global Training grant to H.L.F.

D.H. is codeveloper of the Ridom SeqSphere+ software mentioned in the paper, a development of Ridom GmbH (Münster, Germany), which is partially owned by him. The other authors declare no conflicts of interest.

FOOTNOTES

    • Received 4 May 2015.
    • Returned for modification 5 June 2015.
    • Accepted 13 June 2015.
    • Accepted manuscript posted online 1 July 2015.
  • Supplemental material for this article may be found at http://dx.doi.org/10.1128/JCM.01193-15.

  • Copyright © 2015 Ruppitsch et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution-Noncommercial-ShareAlike 3.0 Unported license, which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original author and source are credited.

REFERENCES

  1. 1.↵
    1. Allerberger F,
    2. Huhulescu S
    . 2015. Pregnancy related listeriosis: treatment and control. Expert Rev Anti Infect Ther 13:395−403. doi:10.1586/14787210.2015.1003809.
    OpenUrlCrossRef
  2. 2.↵
    1. Allerberger F,
    2. Bagó Z,
    3. Huhulescu S,
    4. Pietzka A
    . 2015. Listeriosis: the dark side of refrigeration and ensiling, p 249–286. In Sing A (ed), Zoonoses—infections affecting humans and animals. focus on public health aspects Springer Verlag, Heidelberg, Germany.
  3. 3.↵
    1. Kasper S,
    2. Huhulescu S,
    3. Auer B,
    4. Heller I,
    5. Karner F,
    6. Würzner R
    . 2009. Epidemiology of listeriosis in Austria. Wien Klin Wochenschr 121:113–119. doi:10.1007/s00508-008-1130-2.
    OpenUrlCrossRefPubMed
  4. 4.↵
    1. Liu D
    . 2006. Identification, subtyping and virulence determination of Listeria monocytogenes, an important foodborne pathogen. J Med Microbiol 55:645–659. doi:10.1099/jmm.0.46495-0.
    OpenUrlCrossRefPubMedWeb of Science
  5. 5.↵
    1. Graves LM,
    2. Swaminathan B
    . 2001. PulseNet standardized protocol for subtyping Listeria monocytogenes by macrorestriction and pulsed-field gel electrophoresis. Int J Food Microbiol 65:55–62. doi:10.1016/S0168-1605(00)00501-8.
    OpenUrlCrossRefPubMedWeb of Science
  6. 6.↵
    1. Murchan S,
    2. Kaufmann ME,
    3. Deplano A,
    4. de Ryck R,
    5. Struelens M,
    6. Zinn CE,
    7. Fussing V,
    8. Salmenlinna S,
    9. Vuopio-Varkila J,
    10. El Solh N,
    11. Cuny C,
    12. Witte W,
    13. Tassios PT,
    14. Legakis N,
    15. van Leeuwen W,
    16. van Belkum A,
    17. Vindel A,
    18. Laconcha I,
    19. Garaizar J,
    20. Haeggman S,
    21. Olsson-Liljequist B,
    22. Ransjo U,
    23. Coombes G,
    24. Cookson B
    . 2003. Harmonization of pulsed-field gel electrophoresis protocols for epidemiological typing of strains of methicillin-resistant Staphylococcus aureus: a single approach developed by consensus in 10 European laboratories and its application for tracing the spread of related strains. J Clin Microbiol 41:1574–1585. doi:10.1128/JCM.41.4.1574-1585.2003.
    OpenUrlAbstract/FREE Full Text
  7. 7.↵
    1. Haase JK,
    2. Didelot X,
    3. Lecuit M,
    4. Korkeala H, L. monocytogenes MLST Study Group,
    5. Achtman M
    . 2014. The ubiquitous nature of Listeria monocytogenes clones: a large-scale multilocus sequence typing study. Environ Microbiol 16:405–416. doi:10.1111/1462-2920.12342.
    OpenUrlCrossRef
  8. 8.↵
    1. van Belkum A,
    2. van Leeuwen W,
    3. Kaufmann ME,
    4. Cookson B,
    5. Forey F,
    6. Etienne J,
    7. Goering R,
    8. Tenover F,
    9. Steward C,
    10. O'Brien F,
    11. Grubb W,
    12. Tassios P,
    13. Legakis N,
    14. Morvan A,
    15. El Solh N,
    16. de Ryck R,
    17. Struelens M,
    18. Salmenlinna S,
    19. Vuopio-Varkila J,
    20. Kooistra M,
    21. Talens A,
    22. Witte W,
    23. Verbrugh H
    . 1998. Assessment of resolution and intercenter reproducibility of results of genotyping Staphylococcus aureus by pulsed-field gel electrophoresis of SmaI macrorestriction fragments: a multicenter study. J Clin Microbiol 36:1653–1659.
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    1. Tenover FC,
    2. Arbeit R,
    3. Archer G,
    4. Biddle J,
    5. Byrne S,
    6. Goering R,
    7. Hancock G,
    8. Hébert GA,
    9. Hill B,
    10. Hollis R,
    11. Jarvis WR,
    12. Kreiswirth B,
    13. Eisner W,
    14. Maslow J,
    15. McDougal LK,
    16. Miller JM,
    17. Mulligan M,
    18. Pfaller MA
    . 1994. Comparison of traditional and molecular methods of typing isolates of Staphylococcus aureus. J Clin Microbiol 32:407–415.
    OpenUrlAbstract/FREE Full Text
  10. 10.↵
    1. Aires-de-Sousa M,
    2. Boye K,
    3. de Lencastre H,
    4. Deplano A,
    5. Enright MC,
    6. Etienne J,
    7. Friedrich A,
    8. Harmsen D,
    9. Holmes A,
    10. Huijsdens XW,
    11. Kearns AM,
    12. Mellmann A,
    13. Meugnier H,
    14. Rasheed JK,
    15. Spalburg E,
    16. Strommenger B,
    17. Struelens MJ,
    18. Tenover FC,
    19. Thomas J,
    20. Vogel U,
    21. Westh H,
    22. Xu J,
    23. Witte W
    . 2006. High interlaboratory reproducibility of DNA sequence-based typing of bacteria in a multicenter study. J Clin Microbiol 44:619–621. doi:10.1128/JCM.44.2.619-621.2006.
    OpenUrlAbstract/FREE Full Text
  11. 11.↵
    1. Jolley KA,
    2. Chan M-S,
    3. Maiden MCJ
    . 2004. mlstdbNet—distributed multi-locus sequence typing (MLST) databases. BMC Bioinformatics 5:86. doi:10.1186/1471-2105-5-86.
    OpenUrlCrossRefPubMed
  12. 12.↵
    1. Salcedo C,
    2. Arreaza L,
    3. Alcalá B,
    4. de la Fuente L,
    5. Vázquez JA
    . 2003. Development of a multilocus sequence typing method for the analysis of Listeria monocytogenes clones. J Clin Microbiol 41:757–762. doi:10.1128/JCM.41.2.757-762.2003.
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    1. Zhang W,
    2. Jayarao BM,
    3. Knabel SJ
    . 2004. Multi-virulence-locus sequence typing of Listeria monocytogenes. Appl Environ Microbiol 70:913–920. doi:10.1128/AEM.70.2.913-920.2004.
    OpenUrlAbstract/FREE Full Text
  14. 14.↵
    1. Ragon M,
    2. Wirth T,
    3. Hollandt F,
    4. Lavenir R,
    5. Lecuit M,
    6. Le Monnier A,
    7. Brisse S
    . 2008. A new perspective on Listeria monocytogenes evolution. PLoS Pathog 4:e1000146. doi:10.1371/journal.ppat.1000146.
    OpenUrlCrossRefPubMed
  15. 15.↵
    1. Turabelidze G,
    2. Lawrence SJ,
    3. Gao H,
    4. Sodergren E,
    5. Weinstock GM,
    6. Abubucker S,
    7. Wylie T,
    8. Mitreva M,
    9. Shaikh N,
    10. Gautom R,
    11. Tarr PI
    . 2013. Precise dissection of an Escherichia coli O157:H7 outbreak by single nucleotide polymorphism analysis. J Clin Microbiol 51:3950–3954. doi:10.1128/JCM.01930-13.
    OpenUrlAbstract/FREE Full Text
  16. 16.↵
    1. Eyre DW,
    2. Golubchik T,
    3. Gordon NC,
    4. Bowden R,
    5. Piazza P,
    6. Batty EM,
    7. Ip CL,
    8. Wilson DJ,
    9. Didelot X,
    10. O'Connor L,
    11. Lay R,
    12. Buck D,
    13. Kearns AM,
    14. Shaw A,
    15. Paul J,
    16. Wilcox MH,
    17. Donnelly PJ,
    18. Peto TE,
    19. Walker AS,
    20. Crook DW
    . 2012. A pilot study of rapid benchtop sequencing of Staphylococcus aureus and Clostridium difficile for outbreak detection and surveillance. BMJ Open 2:pii=e001124. doi:10.1136/bmjopen-2012-001124.
    OpenUrlCrossRef
  17. 17.↵
    1. Mellmann A,
    2. Harmsen D,
    3. Cummings CA,
    4. Zentz EB,
    5. Leopold SR,
    6. Rico A,
    7. Prior K,
    8. Szczepanowski R,
    9. Ji Y,
    10. Zhang W,
    11. McLaughlin SF,
    12. Henkhaus JK,
    13. Leopold B,
    14. Bielaszewska M,
    15. Prager R,
    16. Brzoska PM,
    17. Moore RL,
    18. Guenther S,
    19. Rothberg JM,
    20. Karch H
    . 2011. Prospective genomic characterization of the German enterohemorrhagic Escherichia coli O104:H4 outbreak by rapid next generation sequencing technology. PLoS One 6:e22751. doi:10.1371/journal.pone.0022751.
    OpenUrlCrossRefPubMed
  18. 18.↵
    1. Maiden MC,
    2. van Rensburg MJ,
    3. Bray JE,
    4. Earle SG,
    5. Ford SA,
    6. Jolley KA,
    7. McCarthy ND
    . 2013. MLST revisited: the gene-by-gene approach to bacterial genomics. Nat Rev Microbiol 11:728–736. doi:10.1038/nrmicro3093.
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Schmid D,
    2. Allerberger F,
    3. Huhulescu S,
    4. Pietzka A,
    5. Amar C,
    6. Kleta S,
    7. Prager R,
    8. Preußel K,
    9. Aichinger E,
    10. Mellmann A
    . 2014. Whole genome sequencing as a tool to investigate a cluster of seven cases of listeriosis in Austria and Germany, 2011-2013. Clin Microbiol Infect 20:431–436. doi:10.1111/1469-0691.12638.
    OpenUrlCrossRefPubMed
  20. 20.↵
    1. Maiden MC,
    2. Bygraves JA,
    3. Feil E,
    4. Morelli G,
    5. Russell JE,
    6. Urwin R,
    7. Zhang Q,
    8. Zhou J,
    9. Zurth K,
    10. Caugant DA,
    11. Feavers IM,
    12. Achtman M,
    13. Spratt BG
    . 1998. Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganisms. Proc Natl Acad Sci U S A 95:3140–3145. doi:10.1073/pnas.95.6.3140.
    OpenUrlAbstract/FREE Full Text
  21. 21.↵
    1. Pichler J,
    2. Much P,
    3. Kasper S,
    4. Fretz R,
    5. Auer B,
    6. Kathan J,
    7. Mann M,
    8. Huhulescu S,
    9. Ruppitsch W,
    10. Pietzka A,
    11. Silberbauer K,
    12. Neumann C,
    13. Gschiel E,
    14. de Martin A,
    15. Schuetz A,
    16. Gindl J,
    17. Neugschwandtner E,
    18. Allerberger F
    . 2009. An outbreak of febrile gastroenteritis associated with jellied pork contaminated with Listeria monocytogenes. Wien Klin Wochenschr 121:149–156. doi:10.1007/s00508-009-1137-3.
    OpenUrlCrossRefPubMed
  22. 22.↵
    1. Fretz R,
    2. Sagel U,
    3. Ruppitsch W,
    4. Pietzka A,
    5. Stoger A,
    6. Huhulescu S,
    7. Heuberger S,
    8. Pichler J,
    9. Much P,
    10. Pfaff G,
    11. Stark K,
    12. Prager R,
    13. Flieger A,
    14. Feenstra O,
    15. Allerberger F
    . 2010. Listeriosis outbreak caused by acid curd cheese Quargel, Austria and Germany 2009. Euro Surveill 15(16):pii=19477. http://www.eurosurveillance.org/ViewArticle.aspx?ArticleId=19543.
    OpenUrl
  23. 23.↵
    1. Fretz R,
    2. Pichler J,
    3. Sagel U,
    4. Much P,
    5. Ruppitsch W,
    6. Pietzka AT,
    7. Stöger A,
    8. Huhulescu S,
    9. Heuberger S,
    10. Appl G,
    11. Werber D,
    12. Stark K,
    13. Prager R,
    14. Flieger A,
    15. Karpísková R,
    16. Pfaff G,
    17. Allerberger F
    . 2010. Update: multinational listeriosis outbreak due to ‘Quargel,’ a sour milk curd cheese, caused by two different L. monocytogenes serotype 1/2a strains, 2009-2010. Euro Surveill 15(16):pii=19543. http://www.eurosurveillance.org/ViewArticle.aspx?ArticleId=19543.
    OpenUrl
  24. 24.↵
    1. Zerbino DR,
    2. Birney E
    . 2008. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res 18:821–829. doi:10.1101/gr.074492.107.
    OpenUrlAbstract/FREE Full Text
  25. 25.↵
    1. Jünemann S,
    2. Sedlazeck FJ,
    3. Prior K,
    4. Albersmeier A,
    5. John U,
    6. Kalinowski J,
    7. Mellmann A,
    8. Goesmann A,
    9. von Haeseler A,
    10. Stoye J,
    11. Harmsen D
    . 2013. Updating benchtop sequencing performance comparison. Nat Biotechnol 31:294–296. doi:10.1038/nbt.2522.
    OpenUrlCrossRefPubMed
  26. 26.↵
    1. Corander J,
    2. Marttinen P,
    3. Sirén J,
    4. Tang J
    . 2008. Enhanced Bayesian modelling in BAPS software for learning genetic structures of populations. BMC Bioinformatics 9:539. doi:10.1186/1471-2105-9-539.
    OpenUrlCrossRefPubMed
  27. 27.↵
    1. van Tonder AJ,
    2. Mistry S,
    3. Bray JE,
    4. Hill DM,
    5. Cody AJ,
    6. Farmer CL,
    7. Klugman KP,
    8. von Gottberg A,
    9. Bentley SD,
    10. Parkhill J,
    11. Jolley KA,
    12. Maiden MC,
    13. Brueggemann AB
    . 2014. Defining the estimated core genome of bacterial populations using a Bayesian decision model. PLoS Comput Biol 10:e1003788. doi:10.1371/journal.pcbi.1003788.
    OpenUrlCrossRefPubMed
  28. 28.↵
    1. Edgar RC
    . 2004. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 32:1792–1797. doi:10.1093/nar/gkh340.
    OpenUrlCrossRefPubMedWeb of Science
  29. 29.↵
    1. Jolley KA,
    2. Maiden MCJ
    . 2014. Using multilocus sequence typing to study bacterial variation: prospects in the genomic era. Future Microbiol 9:623–630. doi:10.2217/fmb.14.24.
    OpenUrlCrossRefPubMedWeb of Science
  30. 30.↵
    1. Kohl TA,
    2. Diel R,
    3. Harmsen D,
    4. Rothgänger J,
    5. Walter KM,
    6. Merker M,
    7. Weniger T,
    8. Niemann S
    . 2014. Whole-genome-based Mycobacterium tuberculosis surveillance: a standardized, portable, and expandable approach. J Clin Microbiol 52:2479–2486. doi:10.1128/JCM.00567-14.
    OpenUrlAbstract/FREE Full Text
  31. 31.↵
    1. Leopold SR,
    2. Goering RV,
    3. Witten A,
    4. Harmsen D,
    5. Mellmann A
    . 2014. Bacterial whole genome sequencing revisited: portable, scalable and standardized analysis for typing and detection of virulence and antibiotic resistance genes. J Clin Microbiol 52:2365–2370. doi:10.1128/JCM.00262-14.
    OpenUrlAbstract/FREE Full Text
  32. 32.↵
    1. Underwood AP,
    2. Jones G,
    3. Mentasti M,
    4. Fry NK,
    5. Harrison TG
    . 2013. Comparison of the Legionella pneumophila population structure as determined by sequence-based typing and whole genome sequencing. BMC Microbiol 13:302. doi:10.1186/1471-2180-13-302.
    OpenUrlCrossRefPubMed
  33. 33.↵
    1. Pietzka AT,
    2. Stöger A,
    3. Huhulescu S,
    4. Allerberger F,
    5. Ruppitsch W
    . 2011. Gene scanning of an internalin B gene fragment using high-resolution melting curve analysis as a tool for rapid typing of Listeria monocytogenes. J Mol Diagn 13:57–63. doi:10.1016/j.jmoldx.2010.11.002.
    OpenUrlCrossRefPubMed
  34. 34.↵
    1. Rychli K,
    2. Müller A,
    3. Zaiser A,
    4. Schoder D,
    5. Allerberger F,
    6. Wagner M,
    7. Schmitz-Esser S
    . 2014. Genome sequencing of Listeria monocytogenes “Quargel” listeriosis outbreak strains reveals two different strains with distinct in vitro virulence potential. PLoS One 9:e89964. doi:10.1371/journal.pone.0089964.
    OpenUrlCrossRefPubMed
  35. 35.↵
    1. Aarestrup FM,
    2. Brown EW,
    3. Detter C,
    4. Gerner-Smidt P,
    5. Gilmour MW,
    6. Harmsen D,
    7. Hendriksen RS,
    8. Hewson R,
    9. Heymann DL,
    10. Johansson K,
    11. Ijaz K,
    12. Keim PS,
    13. Koopmans M,
    14. Kroneman A,
    15. Lo Fo Wong D,
    16. Lund O,
    17. Palm D,
    18. Sawanpanyalert P,
    19. Sobel J,
    20. Schlundt J
    . 2012. Integrating genome-based informatics to modernize global disease monitoring, information sharing, and response. Emerg Infect Dis 18:e1. doi:10.3201/eid/1811.120453.
    OpenUrlCrossRefPubMed
  36. 36.↵
    1. Mellmann A,
    2. Friedrich AW,
    3. Rosenkötter N,
    4. Rothgänger J,
    5. Karch H,
    6. Reintjes R,
    7. Harmsen D
    . 2006. Automated DNA sequence-based early warning system for the detection of methicillin-resistant Staphylococcus aureus outbreaks. PLoS Med 3:e33. doi:10.1371/journal.pmed.0030033.
    OpenUrlCrossRefPubMed
PreviousNext
Back to top
Download PDF
Citation Tools
Defining and Evaluating a Core Genome Multilocus Sequence Typing Scheme for Whole-Genome Sequence-Based Typing of Listeria monocytogenes
Werner Ruppitsch, Ariane Pietzka, Karola Prior, Stefan Bletz, Haizpea Lasa Fernandez, Franz Allerberger, Dag Harmsen, Alexander Mellmann
Journal of Clinical Microbiology Aug 2015, 53 (9) 2869-2876; DOI: 10.1128/JCM.01193-15

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print

Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this Journal of Clinical Microbiology article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Defining and Evaluating a Core Genome Multilocus Sequence Typing Scheme for Whole-Genome Sequence-Based Typing of Listeria monocytogenes
(Your Name) has forwarded a page to you from Journal of Clinical Microbiology
(Your Name) thought you would be interested in this article in Journal of Clinical Microbiology.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Defining and Evaluating a Core Genome Multilocus Sequence Typing Scheme for Whole-Genome Sequence-Based Typing of Listeria monocytogenes
Werner Ruppitsch, Ariane Pietzka, Karola Prior, Stefan Bletz, Haizpea Lasa Fernandez, Franz Allerberger, Dag Harmsen, Alexander Mellmann
Journal of Clinical Microbiology Aug 2015, 53 (9) 2869-2876; DOI: 10.1128/JCM.01193-15
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • MATERIALS AND METHODS
    • RESULTS
    • DISCUSSION
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

Related Articles

Cited By...

About

  • About JCM
  • Editor in Chief
  • Board of Editors
  • Editor Conflicts of Interest
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Article Types
  • Resources for Clinical Microbiologists
  • Ethics
  • Contact Us

Follow #JClinMicro

@ASMicrobiology

       

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

 

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Print ISSN: 0095-1137; Online ISSN: 1098-660X