Whole genome sequencing options for bacterial strain typing and epidemiologic analysis based on single nucleotide polymorphism versus gene-by-gene–based approaches

The use of next generation sequencing for improving food safety: Translation into practice

Next Generation Sequencing (NGS) combined with powerful bioinformatic approaches are revolutionising food microbiology. Whole genome sequencing (WGS) of single isolates allows the most detailed comparison possible hitherto of individual strains. The two principle approaches for strain discrimination, single nucleotide polymorphism (SNP) analysis and genomic multi-locus sequence typing (MLST) are showing concordant results for phylogenetic clustering and are complementary to each other. Metabarcoding and metagenomics, applied to total DNA isolated from either food materials or the production environment, allows the identification of complete microbial populations. Metagenomics identifies the entire gene content and when coupled to transcriptomics or proteomics, allows the identification of functional capacity and biochemical activity of microbial populations. The focus of this review is on the recent use and future potential of NGS in food microbiology and on current challenges. Guidance is provided for new users, such as public health departments and the food industry, on the implementation of NGS and how to critically interpret results and place them in a broader context. The review aims to promote the broader application of NGS technologies within the food industry as well as highlight knowledge gaps and novel applications of NGS with the aim of driving future research and increasing food safety outputs from its wider use.

A core genome approach that enables prospective and dynamic monitoring of infectious outbreaks

Whole-genome sequencing is increasingly adopted in clinical settings to identify pathogen transmissions, though largely as a retrospective tool. Prospective monitoring, in which samples are continuously added and compared to previous samples, can generate more actionable information. To enable prospective pathogen comparison, genomic relatedness metrics based on single-nucleotide differences must be consistent across time, efficient to compute and reliable for a large variety of samples. The choice of genomic regions to compare, i.e., the core genome, is critical to obtain a good metric. We propose a novel core genome method that selects conserved sequences in the reference genome by comparing its k-mer content to that of publicly available genome assemblies. The conserved-sequence genome is sample set-independent, which enables prospective pathogen monitoring. Based on clinical data sets of 3436 S. aureus, 1362 K. pneumoniae and 348 E. faecium samples, ROC curves demonstrate that the conserved-sequence genome disambiguates same-patient samples better than a core genome consisting of conserved genes. The conserved-sequence genome confirms outbreak samples with high sensitivity: in a set of 2335 S. aureus samples, it correctly identifies 44 out of 44 known outbreak samples, whereas the conserved-gene method confirms 38 known outbreak samples.

Integrating Whole-Genome Sequencing Data Into Quantitative Risk Assessment of Foodborne Antimicrobial Resistance: A Review of Opportunities and Challenges

Whole-genome sequencing (WGS) will soon replace traditional phenotypic methods for routine testing of foodborne antimicrobial resistance (AMR). WGS is expected to improve AMR surveillance by providing a greater understanding of the transmission of resistant bacteria and AMR genes throughout the food chain, and therefore support risk assessment activities. At this stage, it is unclear how WGS data can be integrated into quantitative microbial risk assessment (QMRA) models and whether their integration will impact final risk estimates or the assessment of risk mitigation measures. This review explores opportunities and challenges of integrating WGS data into QMRA models that follow the Codex Alimentarius Guidelines for Risk Analysis of Foodborne AMR. We describe how WGS offers an opportunity to enhance the next-generation of foodborne AMR QMRA modeling. Instead of considering all hazard strains as equally likely to cause disease, WGS data can improve hazard identification by focusing on those strains of highest public health relevance. WGS results can be used to stratify hazards into strains with similar genetic profiles that are expected to behave similarly, e.g., in terms of growth, survival, virulence or response to antimicrobial treatment. The QMRA input distributions can be tailored to each strain accordingly, making it possible to capture the variability in the strains of interest while decreasing the uncertainty in the model. WGS also allows for a more meaningful approach to explore genetic similarity among bacterial populations found at successive stages of the food chain, improving the estimation of the probability and magnitude of exposure to AMR hazards at point of consumption. WGS therefore has the potential to substantially improve the utility of foodborne AMR QMRA models. However, some degree of uncertainty remains in relation to the thresholds of genetic similarity to be used, as well as the degree of correlation between genotypic and phenotypic profiles. The latter could be improved using a functional approach based on prediction of microbial behavior from a combination of 'omics' techniques (e.g., transcriptomics, proteomics and metabolomics). We strongly recommend that methodologies to incorporate WGS data in risk assessment be included in any future revision of the Codex Alimentarius Guidelines for Risk Analysis of Foodborne AMR.

Genomic analysis of Staphylococcus aureus along a pork production chain and in the community, Shandong Province, China

Livestock-associated meticillin-resistant Staphylococcus aureus (LA-MRSA) is an increasingly important public health concern worldwide; however, data on LA-MRSA from Asian countries is scarce. As such, a comprehensive molecular epidemiological survey of S. aureus along a pork production chain and in the community was undertaken in Shandong Province, China. spa typing and whole-genome sequencing were used to survey the occurrence and potential transmission of S. aureus in various sectors, including 899 porcine samples (snout or skin swabs, carcass swabs and pork portions), 845 human nasal samples and 239 environmental samples from commercial farms, a slaughterhouse, a pork wholesale market and the surrounding community. MRSA was detected in higher frequencies in samples from two commercial pig farms (pigs, 49%; farm workers, 64%; environmental samples, 16%) than in samples from the slaughterhouse (fatteners, 8.2%; carcasses, 1.1%; operation workers, 0%; environmental samples, 3.8%), the pork wholesale market (pork, 14%; sellers, 0%) and individuals in the community (6.8%). There were significant differences in population structures, antimicrobial susceptibility profiles, and the presence of resistance and virulence genes between human- and pig-associated isolates. The phylogenetic analysis confirmed the dissemination of LA-MRSA between various segments along the pork production chain. However, MRSA of the same sequence type was not found to be disseminated between the commercial farms and the surrounding communities. Furthermore, one MRSA ST398 was observed, and a novel CC9 variant ST3597 was detected within the chain. The high MRSA carriage rates and the emergence of a new MRSA CC9 variant identified in this study highlight the need for MRSA surveillance.

Current strategy for local- to global-level molecular epidemiological characterisation of global antimicrobial resistance surveillance system pathogens

The prime goal of molecular epidemiology is to identify the origin and evolution of pathogens, which can potentially influence the public health worldwide. Traditional methods provide limited information which is not sufficient for outbreak investigation and studying transmission dynamics. The recent advancement of next-generation sequencing had a major impact on molecular epidemiological studies. Currently, whole-genome sequencing (WGS) has become the gold standard typing method, especially for clinically significant pathogens. Here, we aimed to describe the application of appropriate molecular typing methods for global antimicrobial resistance surveillance system pathogens based on the level of discrimination and epidemiological settings. This shows that sequence-based methods such as multi-locus sequence typing (MLST) are widely used due to cost-effectiveness and database accessibility. However, WGS is the only method of choice for studying Escherichia coli and Shigella spp. WGS is shown to have higher discrimination than other methods in typing Klebsiella pneumoniae, Acinetobacter baumannii and Salmonella spp. due to its changing accessory genome content. For Gram positives such as Streptococcus pneumoniae, WGS would be preferable to understand the evolution of the strains. Similarly, for Staphylococcus aureus, combination of MLST, staphylococcal protein A or SCCmec typing along with WGS could be the choice for epidemiological typing of hospital- and community-acquired strains. This review highlights that combinations of different typing methods should be used to get complete information since no one standalone method is sufficient to study the varying genome diversity.

Genus-wide Leptospira core genome multilocus sequence typing for strain taxonomy and global surveillance

Leptospira is a highly heterogeneous bacterial genus that can be divided into three evolutionary lineages and >300 serovars. The causative agents of leptospirosis are responsible of an emerging zoonotic disease worldwide. To advance our understanding of the biodiversity of Leptospira strains at the global level, we evaluated the performance of whole-genome sequencing (WGS) as a genus-wide strain classification and genotyping tool. Herein we propose a set of 545 highly conserved loci as a core genome MLST (cgMLST) genotyping scheme applicable to the entire Leptospira genus, including non-pathogenic species. Evaluation of cgMLST genotyping was undertaken with 509 genomes, including 327 newly sequenced genomes, from diverse species, sources and geographical locations. Phylogenetic analysis showed that cgMLST defines species, clades, subclades, clonal groups and cgMLST sequence types (cgST), with high precision and robustness to missing data. Novel Leptospira species, including a novel subclade named S2 (saprophytes 2), were identified. We defined clonal groups (CG) optimally using a single-linkage clustering threshold of 40 allelic mismatches. While some CGs such as L. interrogans CG6 (serogroup Icterohaemorrhagiae) are globally distributed, others are geographically restricted. cgMLST was congruent with classical MLST schemes, but had greatly improved resolution and broader applicability. Single nucleotide polymorphisms within single cgST groups was limited to <30 SNPs, underlining a potential role for cgMLST in epidemiological surveillance. Finally, cgMLST allowed identification of serogroups and closely related serovars. In conclusion, the proposed cgMLST strategy allows high-resolution genotyping of Leptospira isolates across the phylogenetic breadth of the genus. The unified genomic taxonomy of Leptospira strains, available publicly at http://bigsdb.pasteur.fr/leptospira, will facilitate global harmonization of Leptospira genotyping, strain emergence follow-up and novel collaborative studies of the epidemiology and evolution of this emerging pathogen.

Identification of bacterial isolates from a public hospital in Australia using complexity-reduced genotyping

Bacterial identification methods used in routine identification of pathogens in medical microbiology include a combination approach of biochemical tests, mass spectrometry or molecular biology techniques. Extensive publicly-available databases of DNA sequence data from pathogenic bacteria have been amassed in recent years; this provides an opportunity for using bacterial genome sequencing for identification purposes. Whole genome sequencing is increasing in popularity, although at present it remains a relatively expensive approach to bacterial identification and typing. Complexity-reduced bacterial genome sequencing provides an alternative. We evaluate genomic complexity-reduction using restriction enzymes and sequencing to identify bacterial isolates. A total of 165 bacterial isolates from hospital patients in the Australian Capital Territory, between 2013 and 2015 were used in this study. They were identified and typed by the Microbiology Department of Canberra Public Hospital, and represented 14 bacterial species. DNA extractions from these samples were processed using a combination of the restriction enzymes PstI with MseI, PstI with HpaII and MseI with HpaII. The resulting sequences (length 30-69 bp) were aligned against publicly available bacterial genome and plasmid sequences. Results of the alignment were processed using a bioinformatics pipeline developed for this project, Currito3.1 DNA Fragment Analysis Software. All 165 samples were correctly identified to genus and species by each of the three combinations of restriction enzymes. A further 35 samples typed to the level of strain identified and compared for consistency with MLST typing data and in silico MLST data derived from the nearest sequenced candidate reference. The high level of agreement between bacterial identification using complexity-reduced genome sequencing and standard hospital identifications indicating that this new approach is a viable alternative for identification of bacterial isolates derived from pathology specimens. The effectiveness of species identification and in particular, strain typing, depends on access to a comprehensive and taxonomically accurate bacterial genome sequence database containing relevant bacterial species and strains.

Challenges in Vibrio parahaemolyticus infections caused by the pandemic clone

Vibrio Parahaemolyticus infections caused by the pandemic clone have become a global public health issue. The pandemic clone includes over ten sequence types and 49 serotypes. Several markers such as toxRS/new, orf8 and genomic islands were considered specific for pandemic strains, but subsequent studies later confirmed a lack of specificity. Thus, identifying stable indicators for the pandemic clone is still an open question. In recent years, several environmental pandemic strains are growing, constituting a new threat to seafood safety and human health. Traditional methods show limited discrimination in studying the microevolution of pandemic strains. For example, multilocus sequence typing divides many pandemic strains into ST3 type, making it difficult to further distinguish the variability within ST3 strains from different contexts. When using a whole genome sequencing-based technique, strains including those with the same sequence type, could be well separated. Whole genome sequencing-based technology also played important roles in dissecting the evolution process and revealing the mechanism underlying rapid serotype conversion within pandemic strains. In addition, the emergence of multiple-antibiotic resistant pandemic strains needs attention. Altogether, we are facing many challenges posed by pandemic V. parahaemolyticus strains, which need to be resolved in future studies.

Direct vertical transmission of ESBL/pAmpC-producing Escherichia coli limited in poultry production pyramid

Extended-spectrum beta-lactamase (ESBL) and plasmidic AmpC (pAmpC) producing Escherichia coli are found in the poultry production even without antibiotic use. The spread of these bacteria has been suggested to occur via imported parent birds, enabling transmission to production level broilers vertically via eggs. We studied transmission of ESBL/pAmpC-producing E. coli and E. coli without antibiotic selection by sampling imported parent birds (n = 450), egg surfaces prior to and after the incubation period (n = 300 and n = 428, respectively) and the laying house environment (n = 20). Samples were additionally taken from embryos (n = 422). To study the prevention of transmission, a competitive exclusion (CE) solution was added onto freshly laid eggs prior to incubation period (n = 150). Results showed carriage of ESBL/pAmpC-producing E. coli in parent birds (26.7%), the environment (5%) and egg surfaces before the incubation period (1.3%), but not from egg surfaces or embryos after the incubation period. Whole genome sequencing revealed ESBL/pAmpC-producing E. coli isolates belonging to clonal lineages ST429 and ST2040. However, the finding of E. coli cultured without antibiotic selection in two (2.2%) embryos strengthens the need to study E. coli transmission in poultry production in more depth. Since ESBL/pAmpC-producing E. coli seem not to persist on egg surfaces, there is no need to use CE solution ex ovo as a prevention method. The results indicate that other routes, such as for example transmission through fomites or horizontal gene transfer by other bacterial species, could be more important than vertical transmission in the spread of resistance in broiler production.

Symposium review: Intramammary infections-Major pathogens and strain-associated complexity

Intramammary infection (IMI) is one of the most costly diseases to the dairy industry. It is primarily due to bacterial infection and the major intramammary pathogens include Escherichia coli, Streptococcus uberis, and Staphylococcus aureus. The severity and outcome of IMI is dependent on several host factors including innate host resistance, energy balance, immune status, parity, and stage of lactation. Additionally, the infecting organism can influence the host immune response and progression of disease. It is increasingly recognized that not only the infecting pathogen species, but also the strain, can affect the transmission, severity, and outcome of IMI. For each of 3 major IMI-associated pathogens, S. aureus, Strep. uberis, and E. coli, specific strains have been identified that are adapted to the intramammary environment. Strain-dependent variation in the host immune response to infection has also been reported. The diversity of strains associated with IMI must be considered if vaccines effective against the full repertoire of mammary pathogenic strains are to be developed. Although important advances have been made recently in understanding the molecular mechanism underpinning strain-specific virulence, further research is required to fully elucidate the cellular and molecular pathogenesis of mammary adapted strains and the role of the strain in influencing the pathophysiology of infection. Improved understanding of molecular pathogenesis of strains associated with bovine IMI will contribute to the development of new control strategies, therapies, and vaccines. The development of enabling technologies such as pathogenomics, transcriptomics, and proteomics can facilitate system-level studies of strain-specific molecular pathogenesis and the identification of key mediators of host-pathogen interactions.

Pulsed Field Gel Electrophoresis: Past, present, and future

Pulsed Field Gel Electrophoresis (PFGE) has been considered for many years the 'gold-standard' for characterizing many pathogenic organisms as well as for subtyping bacterial species causing infection outbreaks. This article reviews the basic principles of PFGE and it includes the main advantages and limitations of the different electrode configurations that have been used in PFGE equipment and their influence on the DNA electrophoretic separation. Remarkably, we summarize here the most relevant theoretical and practical aspects that we have learned for more than 20 years developing and using the miniaturized PFGE systems. We also discussed the theoretical aspects related to DNA migration in PFGE agarose gels. It served as the basis for simulating the DNA electrophoretic patterns in CHEF mini gels and mini-chambers during experimental design and optimization. A critical comparison between standard and miniaturized PFGE systems, as well as the enzymatic and non-enzymatic methods for intact immobilized DNA preparation, is provided throughout the review. The PFGE current applications, advantages, limitations and future challenges of the methodology are also discussed.

Evaluation of a gene-by-gene approach for prospective whole-genome sequencing-based surveillance of multidrug resistant Mycobacterium tuberculosis

Whole-genome sequencing (WGS) offers unprecedented resolution for tracking Mycobacterium tuberculosis transmission and antibiotic-resistance spread. Still, the establishment of standardized WGS-based pipelines and the definition of epidemiological clusters based on genetic relatedness are under discussion. We aimed to implement a dynamic gene-by-gene approach, fully relying on freely available software, for prospective WGS-based tuberculosis surveillance, demonstrating its application for detecting transmission chains by retrospectively analysing all M/XDR strains isolated in 2013-2017 in Portugal. We observed a good correlation between genetic relatedness and epidemiological links, with strongly epilinked clusters displaying mean pairwise allele differences (AD) always below 0.3% (ratio of mean AD over the total number of shared loci between same-cluster strains). This data parallels the genetic distances acquired by the core-SNV analysis, while providing higher resolution and epidemiological concordance than MIRU-VNTR genotyping. The dynamic analysis of strain sub-sets (i.e., increasing the number of shared loci within each sub-set) also strengthens the confidence in detecting epilinked clusters. This gene-by-gene strategy also offers several practical benefits (e.g., reliance on freely-available software, scalability and low computational requirements) that further consolidated its suitability for a timely and robust prospective WGS-based laboratory surveillance of M/XDR-TB cases.

Whole-genome sequencing of the rare Salmonella enterica serovar Anfo isolated from food handlers

Field investigations were conducted after a small cluster of food poisoning involving six cases was reported. While no stool samples were available from the cases for microbiological testing, Salmonella species was found to be present in the stools of food handlers with gastroenteritis symptoms. Four Salmonella isolates recovered from the food handlers were retrospectively investigated at the genome level using whole-genome sequencing (WGS). WGS showed that S. Anfo (antigenic formulae 39:y:1,2), a rarely isolated serovar, caused infections in the food handlers. S. Anfo analysed in this study contained virulence factors required for causing disease. They did not contain any antibiotic resistance genes or plasmid. The epidemiologically related isolates differed to each other by a maximum of one single nucleotide polymorphism. WGS was useful in identifying rare Salmonella serovars and it is potentially more cost-effective than traditional serotyping methods. It can also confidently group epidemiologically related isolates belonging to S. Anfo.

A new DNA sensor system for specific and quantitative detection of mycobacteria

In the current study, we describe a novel DNA sensor system for specific and quantitative detection of mycobacteria, which is the causative agent of tuberculosis. Detection is achieved by using the enzymatic activity of the mycobacterial encoded enzyme topoisomerase IA (TOP1A) as a biomarker. The presented work is the first to describe how the catalytic activities of a member of the type IA family of topoisomerases can be exploited for specific detection of bacteria. The principle for detection relies on a solid support anchored DNA substrate with dual functions namely: (1) the ability to isolate mycobacterial TOP1A from crude samples and (2) the ability to be converted into a closed DNA circle upon reaction with the isolated enzyme. The DNA circle can act as a template for rolling circle amplification generating a tandem repeat product that can be visualized at the single molecule level by fluorescent labelling. This reaction scheme ensures specific, sensitive, and quantitative detection of the mycobacteria TOP1A biomarker as demonstrated by the use of purified mycobacterial TOP1A and extracts from an array of non-mycobacteria and mycobacteria species. When combined with mycobacteriophage induced lysis as a novel way of effective yet gentle extraction of the cellular content from the model Mycobacterium smegmatis, the DNA sensor system allowed detection of mycobacteria in small volumes of cell suspensions. Moreover, it was possible to detect M. smegmatis added to human saliva. Depending on the composition of the sample, we were able to detect 0.6 or 0.9 million colony forming units (CFU) per mL of mycobacteria, which is within the range of clinically relevant infection numbers. We, therefore, believe that the presented assay, which relies on techniques that can be adapted to limited resource settings, may be the first step towards the development of a new point-of-care diagnostic test for tuberculosis.

Whole genome sequencing of Salmonella Chester reveals geographically distinct clusters, Norway, 2000 to 2016

IntroductionDuring summer 2016, Norway observed an increase in Salmonella enterica subsp. enterica serovar Chester cases among travellers to Greece.AimOur aim was to investigate genetic relatedness of S. Chester for surveillance and outbreak detection by core genome multilocus sequence typing (cgMLST) and compare the results to genome mapping.MethodsWe included S. Chester isolates from 51 cases of salmonellosis between 2000 and 2016. Paired-end sequencing (2 × 250 bp) was performed on Illumina MiSeq. Genetic relatedness by cgMLST for Salmonella enterica subsp. enterica, including 3,002 genes and seven housekeeping genes, was compared by reference genome mapping with CSI Phylogeny version 1.4 and conventional MLST.ResultsConfirmed travel history was available for 80% of included cases, to Europe (n = 13), Asia (n = 12) and Africa (n = 16). Isolates were distributed into four phylogenetic clusters corresponding to geographical regions. Sequence type (ST) ST411 and a single-locus variant ST5260 (n = 17) were primarily acquired in southern Europe, ST1954 (n = 15) in Africa, ST343 (n = 11) and ST2063 (n = 8) primarily in Asia. Part of the European cluster was further divided into a Greek (n = 10) and a Cypriot (n = 4) cluster. All isolates in the African cluster displayed resistance to ≥ 1 class of antimicrobials, while resistance was rare in the other clusters.ConclusionWhole genome sequencing of S. Chester in Norway showed four geographically distinct clusters, with a possible outbreak occurring during summer 2016 related to Greece. We recommend public health institutes to implement cgMLST-based real-time Salmonella enterica surveillance for early and accurate detection of future outbreaks and further development of cluster cut-offs.

Hospitalized Pets as a Source of Carbapenem-Resistance

The massive and irrational use of antibiotics in livestock productions has fostered the occurrence and spread of resistance to “old class antimicrobials.” To cope with that phenomenon, some regulations have been already enforced in the member states of the European Union. However, a role of livestock animals in the relatively recent alerts on the rapid worldwide increase of resistance to last-choice antimicrobials as carbapenems is very unlikely. Conversely, these antimicrobials are increasingly administered in veterinary hospitals whose role in spreading bacteria or mobile genetic elements has not adequately been addressed so far. A cross-sectional study was carried out on 105 hospitalized and 100 non-hospitalized pets with the aim of measuring the prevalence of carbapenem-resistant Gram-negative bacteria (GNB) colonizing dogs and cats, either hospitalized or not hospitalized and estimating the relative odds. Stool samples were inoculated on MacConkey agar plates containing 1 mg/L imipenem which were then incubated aerobically at 37°C ± 1 for 48 h. Isolated bacteria were identified first by Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry and were confirmed by 16S rRNA sequencing. The genetic basis of resistance was investigated using PCR methods, gene or whole genome sequencing (WGS). The prevalence of pets harboring carbapenem-resistant bacteria was 11.4 and 1.0% in hospitalized and not-hospitalized animals, respectively, with an odds ratio of 12.8 (p < 0.01). One pet carried two diverse isolates. Overall, 14 gram-negative non-fermenting bacteria, specifically, one Acinetobacter radioresistens, five Acinetobacter baumannii, six Pseudomonas aeruginosa and two Stenotrophomonas maltophilia were isolated. The Acinetobacter species carried acquired carbapenemases genes encoded by bla_NDM-1 and bla_OXA-23. In contrast, Pseudomonas phenotypic resistance was associated with the presence of mutations in the oprD gene. Notably, inherent carbapenem-resistant isolates of S. maltophilia were also resistant to the first-line recommended chemotherapeutic trimethoprim/sulfamethoxazole. This study estimates the risk of colonization by carbapenem-resistant non-fermenting GNB in pets hospitalized in veterinary tertiary care centers and highlights their potential role in spreading resistance genes among the animal and human community. Public health authorities should consider extending surveillance systems and putting the release of critical antibiotics under more strict control in order to manage the infection/colonization of pets in veterinary settings.

USA300 methicillin-resistant Staphylococcus aureus in Stockholm, Sweden, from 2008 to 2016

Methicillin-resistant Staphylococcus aureus (MRSA) USA300 isolates have been recognized globally, not only in community but also in healthcare settings. USA300 isolates were initially resistant only to methicillin, but resistance to non-β-lactams has emerged with time. To evaluate the prevalence and antimicrobial susceptibility of USA300 isolates in Stockholm, we conducted a nine-year retrospective study. Of 5359 consecutive MRSA cases in Stockholm, isolates from 285 cases were USA300 strains according to the pulsed-field gel electrophoresis pattern. Of these cases, repeated isolates with altered antibiotic resistance patterns were observed in six individuals. Therefore, antimicrobial susceptibility testing was performed on totally 291 isolates. To study the phylogenetic relatedness of isolates in transmission events and genomic resistance traits, 35 isolates were further studied by whole genome sequencing (WGS). The incidence of MRSA was increased from 17.6 per 100,000 inhabitants in 2008 to 37.3 per 100,000 inhabitants in 2016, while the proportion of USA300 cases declined from 6.6% in 2008 to 2.6% in 2016. Among the USA300 isolates, 73.5% were community-associated, 21.3% healthcare-associated, and 5.2% had unknown acquisition. The highest resistance rate among non-β-lactams was found in erythromycin (86%), followed by fluoroquinolones (68-69%). 57% of the isolates were resistant to both erythromycin and fluoroquinolone. Simultaneous resistance to four non-β-lactam antibiotic classes was found in six isolates. Four isolates were susceptible to all non-β-lactam antibiotics. Ceftaroline, daptomycin, linezolid, mupirocin, rifampicin, teicoplanin, telavancin, trimethoprim-sulfamethoxazole and vancomycin retained full activity in the study. WGS analysis indicated that isolates from an outbreak were phylogenetically closely related. In conclusion, USA300 MRSA isolates in Stockholm have neither been limited to the community setting, nor remained susceptible to non-β-lactam agents. WGS is becoming a useful tool in tracing transmission events. The results herein provide the most up-to-date and comprehensive information regarding status of USA300 strains in this geographic area.

Typing and Species Identification of Clinical Klebsiella Isolates by Fourier Transform Infrared Spectroscopy and Matrix-Assisted Laser Desorption Ionization-Time of Flight Mass Spectrometry

Klebsiella pneumoniae and related species are frequent causes of nosocomial infections and outbreaks. Therefore, quick and reliable strain typing is crucial for the detection of transmission routes in the hospital. The aim of this study was to evaluate Fourier transform infrared spectroscopy (FTIR) and matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) as rapid methods for typing clinical Klebsiella isolates in comparison to whole-genome sequencing (WGS), which was considered the gold standard for typing and identification. Here, 68 clinical Klebsiella strains were analyzed by WGS, FTIR, and MALDI-TOF MS. FTIR showed high discriminatory power in comparison to the WGS reference, whereas MALDI-TOF MS exhibited a low ability to type the isolates. MALDI-TOF mass spectra were further analyzed for peaks that showed high specificity for different Klebsiella species. Phylogenetic analysis revealed that the Klebsiella isolates comprised three different species: K. pneumoniae, K. variicola, and K. quasipneumoniae Genome analysis showed that MALDI-TOF MS can be used to distinguish K. pneumoniae from K. variicola due to shifts of certain mass peaks. The peaks were tentatively identified as three ribosomal proteins (S15p, L28p, L31p) and one stress response protein (YjbJ), which exhibit amino acid differences between the two species. Overall, FTIR has high discriminatory power to recognize the clonal relationship of isolates, thus representing a valuable tool for rapid outbreak analysis and for the detection of transmission events due to fast turnaround times and low costs per sample. Furthermore, specific amino acid substitutions allow the discrimination of K. pneumoniae and K. variicola by MALDI-TOF MS.

Genomic Epidemiology of Campylobacter jejuni Transmission in Israel

Objectives: Campylobacter jejuni is responsible for 80% of Campylobacter infections in Israel, a country with a high incidence reaching 91/100,000 population. We studied the phylogeny, diversity and prevalence of virulence factors using whole genome sequencing (WGS) of a national sample of C. jejuni clinical, food, and animal isolates collected over a 10-year period (2003-2012). Methods: C. jejuni isolates (n = 263) were subject to WGS using Illumina sequencing (PE 250bpx2). Raw reads and de novo assemblies were analyzed with the BioNumerics whole genome MLST (wgMLST) pipeline. Reads were screened for 71 virulence genes by the SRST2 script. Allelic profiles were analyzed to create minimum spanning trees and allelic core distances were investigated to determine a reliable cutoff for strain determination. Results: wgMLST analysis of 263 C. jejuni isolates indicated significant diversity among the prevalent clonal complexes (CCs) with CC-21 and CC-353 being the most diverse, and CC-574 the most clonal. Within CC-21, sequence type (ST)-1359 created a separate clade. Human, poultry and bovine isolates clustered together across the different STs. Forty four percent of studied isolates were assigned to 29 genetic clusters. Temporal and geographical relatedness were found among the minority of clusters, while most phylogenetically associated cases appeared diffuse and unassociated epidemiologically. The majority of virulence factors were highly prevalent across the dataset and not associated with genotype, source of isolation or invasiveness. Conversely, all 13 genes associated with type VI secretion system (T6SS) were lineage-related and identified in only 18% of the isolates. T6SS was detected in 95.2% of ST-1359, a common type in Israel. Conclusions: wgMLST supported the assessment that poultry and cattle are likely food sources of infection in Israel. Substantial genetic clustering among C. jejuni isolates suggested multiple point source and diffuse outbreaks that were previously unreported in Israel. The high prevalence of T6SS among ST-1359 isolates is unique to Israel, and requires further investigation. This study exemplifies the importance of studying foodborne pathogens using advanced genomic approaches across the entire spectrum of One Health.

Core Genome Multilocus Sequence Typing and Single Nucleotide Polymorphism Analysis in the Epidemiology of Brucella melitensis Infections

The use of whole-genome sequencing (WGS) using next-generation sequencing (NGS) technology has become a widely accepted method for microbiology laboratories in the application of molecular typing for outbreak tracing and genomic epidemiology. Several studies demonstrated the usefulness of WGS data analysis through single-nucleotide polymorphism (SNP) calling from a reference sequence analysis for Brucella melitensis, whereas gene-by-gene comparison through core-genome multilocus sequence typing (cgMLST) has not been explored so far. The current study developed an allele-based cgMLST method and compared its performance to that of the genome-wide SNP approach and the traditional multilocus variable-number tandem repeat analysis (MLVA) on a defined sample collection. The data set was comprised of 37 epidemiologically linked animal cases of brucellosis as well as 71 isolates with unknown epidemiological status, composed of human and animal samples collected in Italy. The cgMLST scheme generated in this study contained 2,704 targets of the B. melitensis 16M reference genome. We established the potential criteria necessary for inclusion of an isolate into a brucellosis outbreak cluster to be ≤6 loci in the cgMLST and ≤7 in WGS SNP analysis. Higher phylogenetic distance resolution was achieved with cgMLST and SNP analysis than with MLVA, particularly for strains belonging to the same lineage, thereby allowing diverse and unrelated genotypes to be identified with greater confidence. The application of a cgMLST scheme to the characterization of B. melitensis strains provided insights into the epidemiology of this pathogen, and it is a candidate to be a benchmark tool for outbreak investigations in human and animal brucellosis.

Intra-Hospital, Inter-Hospital and Intercontinental Spread of ST78 MRSA From Two Neonatal Intensive Care Unit Outbreaks Established Using Whole-Genome Sequencing

From 2009 to 2011 [transmission period (TP) 1] and 2014 to 2017 (TP2), two outbreaks involving community-associated clonal complex (CC) 88-MRSA spa types t186 and t786, respectively, occurred in the Neonatal Intensive Care Unit (NICU) of an Irish hospital (H1). This study investigated the relatedness of these isolates, their relationship to other CC88 MRSA from Ireland and their likely geographic origin, using whole-genome sequencing (WGS). All 28 CC88-MRSA isolates identified at the Irish National MRSA Reference Laboratory between 2009 and 2017 were investigated including 20 H1 patient isolates, two H1 isolates recovered from a single healthcare worker (HCW) 2 years apart, three patient isolates from a second hospital (H2) and one patient isolate from each of three different hospitals (H3, H4, and H5). All isolates underwent DNA microarray profiling. Thirteen international isolates with similar microarray profiles to at least one Irish isolate were selected from an extensive global database. All isolates underwent Illumina MiSeq WGS. The majority of Irish isolates (25/28; all H1 isolates, two H2 isolates and the H3 isolate) were identified as ST78-MRSA-IVa and formed a large cluster, exhibiting 1–71 pairwise allelic differences, in a whole-genome MLST-based minimum spanning tree (MST) involving all Irish isolates. A H1/H2, H1/H3, and H1 HCW/patient isolate pair each exhibited one allelic difference. The TP2 isolates were characterised by a different spa type and the loss of hsdS. The three remaining Irish isolates (from H2, H4, and H5) were identified as ST88-MRSA-IVa and dispersed at the opposite end of the MST, exhibiting 81–211 pairwise allelic differences. Core-genome MLST and sequence-based plasmid analysis revealed the recent shared ancestry of Irish and Australian ST78-MRSA-IVa, and of Irish and French/Egyptian ST88-MRSA-IVa. This study revealed the homogeneity of isolates recovered during two NICU outbreaks (despite spa type and hsdS carriage variances), HCW involvement in the outbreak transmission chain and the strain's spread to two other Irish hospitals. The outbreak strain, CC88/ST78-MRSA-IVa, was likely imported from Australia, where it is prevalent. CC88/ST88-MRSA-IVa was also identified in Irish hospitals and was likely imported from Africa, where it is predominant, and/or a country with a large population of African descent.