Comparison of Whole-Genome Sequencing Methods for Analysis of Three Methicillin-Resistant Staphylococcus aureus Outbreaks

Comparison of the effectiveness of core genome multilocus sequence typing and polymerase chain reaction-based open reading frame typing in tracing nosocomial methicillin-resistant Staphylococcus aureus transmission

INTRODUCTION

A clonal shift from staphylococcal cassette chromosome mec (SCCmec) type II/ST5 methicillin-resistant Staphylococcus aureus (MRSA) to SCCmec type IV/clonal complex (CC)1 MRSA has occurred rapidly in Japan. Our previous research in a geriatric hospital found SCCmec type IV/CC1 MRSA prevalence in long-term care wards. Due to intensive personal care requirements, frequent contact with healthcare providers can potentially cause unintentional nosocomial MRSA transmission. We performed polymerase chain reaction-based open reading frame typing (POT) and core genome multilocus sequence typing (cgMLST) to investigate the occurrence of nosocomial transmission and to compare the results of these methods.

METHODS

POT and whole genome sequencing were performed in 83 MRSA isolates. Commercial automated software (Ridom SeqSphere+) was used to perform cgMLST. MRSA isolates with 0-8 allelic differences were considered related, and medical records were consulted in these cases.

RESULTS

SCCmec type IV/CC1 MRSA was the most frequently detected clone (n = 56, 67.5 %), which was divided into 14 POT types, followed by SCCmec type I/ST8 (n = 9) and SCCmec type IV/ST8 (n = 8). Identical POT types were found across 7 of 11 wards. However, cgMLST analysis identified only three cases (six strains) of high genetic similarity, indicating nosocomial transmission; only one involved SCCmec type IV/CC1 (two strains). The mean allelic difference in the core genomes between strains with identical POT types in the same ward was 55.3 ± 22.0.

CONCLUSIONS

The cgMLST method proved more effective for identifying nosocomial transmissions compared to POT, highlighting its utility in tracking MRSA spread in healthcare settings.

Evaluation of nationwide analysis surveillance for methicillin-resistant Staphylococcus aureus within Genomic Medicine Sweden

Background. National epidemiological investigations of microbial infections greatly benefit from the increased information gained by whole-genome sequencing (WGS) in combination with standardized approaches for data sharing and analysis.Aim. To evaluate the quality and accuracy of WGS data generated by different laboratories but analysed by joint pipelines to reach a national surveillance approach.Methods. A national methicillin-resistant Staphylococcus aureus (MRSA) collection of 20 strains was distributed to nine participating laboratories that performed in-house procedures for WGS. Raw data were shared and analysed by three pipelines: 1928 Diagnostics, JASEN (GMS pipeline) and CLC-Genomics Workbench. The outcomes were compared according to quality, correct strain identification and genetic distances.Results. One isolate contained intraspecies contamination and was excluded from further analysis. The mean sequencing depth varied between sites and technologies. However, all analysis methods identified 12 strains that belonged to one of five outbreak clusters. The cut-off definition was set to <10 allele differences for core genome multilocus sequence typing (cgMLST) and <20 genetic differences for SNP analysis in a pairwise comparison.Conclusions. MRSA isolates, which are whole genome sequenced by different laboratories and analysed using the same bioinformatic pipelines, yielded comparable results for outbreak clustering for both cgMLST and SNP, using the 1928 analysis pipeline. In this study, JASEN was best suited to analyse Illumina data and CLC to analyse within respective technology. In the future, real-time sharing of data and harmonized analysis within the Genomic Medicine Sweden consortium will further facilitate investigations of outbreaks and transmission routes.

Comparison of non-aureus staphylococcal and mammaliicoccal species found in both composite milk and bulk-tank milk samples of dairy cows collected in tandem

Non-aureus staphylococci and the closely related mammaliicoccal species (NASM) are the most common causes of bovine subclinical mastitis on modern dairy farms and are highly prevalent in bulk-tank milk. The purpose of this study was to determine the distribution of NASM in both composite cow milk (CCM) and bulk-tank milk (BTM) samples collected in tandem in commercial Flemish dairy herds and to estimate the origin of the different (subgroups of) NASM species present in BTM by applying strain typing (random amplification of polymorphic DNA or random amplified DNA [RAPD]). A single cross-sectional sampling was performed over 5 herds that volunteered to participate in the study. Composite cow milk samples (n = 356) were collected from all lactating cows (except those with clinical mastitis) during a milking in tandem with 6 BTM samples per herd sequentially collected immediately post that milking (n = 30). In total, 421 and 80 NASM isolates were recovered and identified by MALDI-TOF mass spectrometry from the CCM and BTM samples, respectively and a total of 21 and 12 different NASM species were identified from CCM and BTM samples, respectively. Staphylococcus cohnii was the most prevalent NASM species found in BTM followed by Staphylococcus haemolyticus, Staphylococcus epidermidis, Mammaliicoccus lentus, and Staphylococcus equorum, whereas from CCM samples the most common species were S. hemolyticus, S. cohnii, S. equorum, S. epidermidis, and Staphylococcus chromogenes. The prevalent NASM species in both CCM and BTM samples was distinct for each herd, corroborating other studies observing a herd-specific NASM microbiota. Random amplified DNA analysis was performed on 9 NASM species (S. chromogenes, S. epidermidis, S. haemolyticus, S. equorum, Mammaliicoccus sciuri, Staphylococcus xylosus, S. cohnii, Staphylococcus debuckii, and M. lentus) because these species were isolated from both sample types in a herd. The same RAPD types were found in both sample types for all NASM species selected for strain typing in varying degrees. When assessing the distribution of NASM species, differences within NASM species should be examined meaning a closer look should be taken at the strain level rather than at the species level only.

Genomic Diversity of Methicillin-Resistant Staphylococcus aureus CC398 Isolates Collected from Diseased Swine in the German National Resistance Monitoring Program GERM-Vet from 2007 to 2019

Livestock-associated methicillin-resistant Staphylococcus aureus (LA-MRSA) clonal complex 398 (CC398) isolates (n = 178) collected in the national resistance monitoring program GERM-Vet from diseased swine in Germany from 2007 to 2019 were investigated for their genomic diversity with a focus on virulence and antimicrobial resistance (AMR) traits. Whole-genome sequencing was followed by molecular typing and sequence analysis. A minimum spanning tree based on core-genome multilocus sequence typing was constructed, and antimicrobial susceptibility testing was performed. Most isolates were assigned to nine clusters. They displayed close phylogenetic relationships but a wide molecular variety, including 13 spa types and 19 known and four novel dru types. Several toxin-encoding genes, including eta, seb, sek, sep, and seq, were detected. The isolates harbored a wide range of AMR properties mirroring the proportions of the classes of antimicrobial agents applied in veterinary medicine in Germany. Multiple novel or rare AMR genes were identified, including the phenicol-lincosamide-oxazolidinone-pleuromutilin-streptogramin A resistance gene cfr, the lincosamide-pleuromutilin-streptogramin A resistance gene vga(C), and the novel macrolide-lincosamide-streptogramin B resistance gene erm(54). Many AMR genes were part of small transposons or plasmids. Clonal and geographical correlations of molecular characteristics and resistance and virulence genes were more frequently observed than temporal relations. In conclusion, this study provides insight into population dynamics of the main epidemic porcine LA-MRSA lineage in Germany over a 13-year-period. The observed comprehensive AMR and virulence properties, most likely resulting from the exchange of genetic material between bacteria, highlighted the importance of LA-MRSA surveillance to prevent further dissemination among swine husbandry facilities and entry into the human community. IMPORTANCE The LA-MRSA-CC398 lineage is known for its low host specificity and frequent multiresistance to antimicrobial agents. Colonized swine and their related surroundings represent a considerable risk of LA-MRSA-CC398 colonization or infection for occupationally exposed people through which such isolates might be further disseminated within the human community. This study provides insight into the diversity of the porcine LA-MRSA-CC398 lineage in Germany. Clonal and geographical correlations of molecular characteristics and resistance and virulence traits were detected and may be associated with the spread of specific isolates through livestock trade, human occupational exposure, or dust emission. The demonstrated genetic variability underlines the lineage's ability to horizontally acquire foreign genetic material. Thus, LA-MRSA-CC398 isolates have the potential to become even more dangerous for various host species, including humans, due to increased virulence and/or limited therapeutic options for infection control. Full-scale LA-MRSA monitoring at the farm, community, and hospital level is therefore essential.

Complete Genome Sequence and Analysis of a ST573 Multidrug-Resistant Methicillin-Resistant Staphylococcus aureus SauR3 Clinical Isolate from Terengganu, Malaysia

Methicillin-resistant Staphylococcus aureus (MRSA) is a World Health Organization-listed priority pathogen. Scarce genomic data are available for MRSA isolates from Malaysia. Here, we present the complete genome sequence of a multidrug-resistant MRSA strain SauR3, isolated from the blood of a 6-year-old patient hospitalized in Terengganu, Malaysia, in 2016. S. aureus SauR3 was resistant to five antimicrobial classes comprising nine antibiotics. The genome was sequenced on the Illumina and Oxford Nanopore platforms and hybrid assembly was performed to obtain its complete genome sequence. The SauR3 genome consists of a circular chromosome of 2,800,017 bp and three plasmids designated pSauR3-1 (42,928 bp), pSauR3-2 (3011 bp), and pSauR3-3 (2473 bp). SauR3 belongs to sequence type 573 (ST573), a rarely reported sequence type of the staphylococcal clonal complex 1 (CC1) lineage, and harbors a variant of the staphylococcal cassette chromosome mec (SCCmec) type V (5C2&5) element which also contains the aac(6')-aph(2″) aminoglycoside-resistance genes. pSauR3-1 harbors several antibiotic resistance genes in a 14,095 bp genomic island (GI), previously reported in the chromosome of other staphylococci. pSauR3-2 is cryptic, whereas pSauR3-3 encodes the ermC gene that mediates inducible resistance to macrolide-lincosamide-streptogramin B (iMLS_B). The SauR3 genome can potentially be used as a reference genome for other ST573 isolates.

Outbreak investigation including molecular characterization of community associated methicillin-resistant Staphylococcus aureus in a primary and secondary school in Eastern Switzerland

At our tertiary children's hospital, infections with newly detected methicillin-resistant Staphylococcus aureus (MRSA) among children attending primary (age 6-12 years) and secondary school (age 13-16 years) nearly doubled in 2018 compared to previous years. This observation initiated an epidemiological outbreak investigation including phenotypic (susceptibility testing) and genotypic (whole genome sequencing) characterization of the isolates. In addition, a cross-sectional study was conducted to determine source of the outbreak, colonization frequency and to identify risk factors for transmission using a questionnaire. As a result, 49 individuals were detected with 57 corresponding isolates. Based on the case definition combined with whole genome sequencing, a core cluster was identified that shared common genetic features and a similar antimicrobial susceptibility pattern (efflux-mediated macrolide resistance, tetracycline susceptibility along with presence of Panton-Valentine leukocidin). Epidemiologic evaluation identified a distinct school as a common risk factor. However, the source of the clustered infections within that school could not be further specified. No further cases could be detected after decolonization of infected and colonized children.

MGTdb: a web service and database for studying the global and local genomic epidemiology of bacterial pathogens

Multilevel genome typing (MGT) enables the genomic characterization of bacterial isolates and the relationships among them. The MGT system describes an isolate using multiple multilocus sequence typing (MLST) schemes, referred to as levels. Thus, for a new isolate, sequence types (STs) assigned at multiple precisely defined levels can be used to type isolates at multiple resolutions. The MGT designation for isolates is stable, and the assignment is faster than the existing approaches. MGT’s utility has been demonstrated in multiple species. This paper presents a publicly accessible web service called MGTdb, which enables the assignment of MGT STs to isolates, along with their storage, retrieval and analysis. The MGTdb web service enables upload of genome data as sequence reads or alleles, which are processed and assigned MGT identifiers. Additionally, any newly sequenced isolates deposited in the National Center for Biotechnology Information’s Sequence Read Archive are also regularly retrieved (currently daily), processed, assigned MGT identifiers and made publicly available in MGTdb. Interactive visualization tools are presented to assist analysis, along with capabilities to download publicly available isolates and assignments for use with external software. MGTdb is currently available for Salmonella enterica serovars Typhimurium and Enteritidis and Vibrio cholerae. We demonstrate the usability of MGTdb through three case studies — to study the long-term national surveillance of S. Typhimurium, the local epidemiology and outbreaks of S. Typhimurium, and the global epidemiology of V. cholerae. Thus, MGTdb enables epidemiological and microbiological investigations at multiple levels of resolution for all publicly available isolates of these pathogens.

Database URL: https://mgtdb.unsw.edu.au

Prevalence and WGS-based characteristics of MRSA isolates in hospitals in Shanghai, China

Methicillin-resistant Staphylococcus aureus (MRSA) isolates remain a serious threat to global health despite a decrease in MRSA infections since 2005. MRSA isolates exhibit great diversity worldwide, and their lineages show geographic variation. In this study, we used whole genome sequencing (WGS) to analyze antibiotic resistance genes and virulence genes, spa, staphylococcal cassette chromosome mec, sequence types (STs), and core genome multilocus sequence typing (cgMLST) of MRSA isolates from patients and environmental surface in hospitals in China to determine their prevalence and molecular traits. The highest number of infections by MRSA isolates was observed in patients aged ≥60 years (69.8%, P < 0.05). We identified a total of 19 STs from 162 MRSA isolates from patients. A significant increase was observed in the incidence of ST764-t002-II MRSA infection, which is replacing ST5-t002-II MRSA as the predominant ST. Similarly, isolates from environmental surface were predominantly ST764-t002-II (47%). Notably, most ST764 isolates (97.7%) carried seb, but not arginine catabolic mobile element (ACME), which differed from ST764 isolates in Japan and Thailand. The potential danger of spreading requires rigorous surveillance of emerging ST764 MRSA isolates. We also found higher resistance to seven antimicrobials [OXA, cefoxitin (FOX), ERY, CFZ, ciprofloxacin (CIP), levofloxacin (LEV), and moxifloxacin (MXF)]. Resistance to gentamicin (38.3%), tetracycline (55.9%), and minocycline (41.5%) were also common. Phenotypic resistance to antimicrobials was associated with resistance genes to its content, and cgMLST clustering suggested a strong link between these strains. Overall, our findings revealed the prevalence and molecular characteristics of MRSA isolates in Shanghai, China, providing a theoretical reference for preventing and controlling clonal transmission of MRSA isolates in hospitals in China.

Metabolic Characteristics of Porcine LA-MRSA CC398 and CC9 Isolates from Germany and China via Biolog Phenotype MicroArrayTM

Livestock-associated methicillin-resistant Staphylococcus aureus (LA-MRSA) is an important zoonotic pathogen, often multi-resistant to antimicrobial agents. Among swine, LA-MRSA of clonal complex (CC) 398 dominates in Europe, Australia and the Americas, while LA-MRSA-CC9 is the main epidemic lineage in Asia. Here, we comparatively investigated the metabolic properties of rare and widespread porcine LA-MRSA isolates from Germany and China using Biolog Phenotype MicroArray technology to evaluate if metabolic variations could have played a role in the development of two different epidemic LA-MRSA clones in swine. Overall, we were able to characterize the isolates' metabolic profiles and show their tolerance to varying environmental conditions. Sparse partial least squares discriminant analysis (sPLS-DA) supported the detection of the most informative substrates and/or conditions that revealed metabolic differences between the LA-MRSA lineages. The Chinese LA-MRSA-CC9 isolates displayed unique characteristics, such as a consistently delayed onset of cellular respiration, and increased, reduced or absent usage of several nutrients. These possibly unfavorable metabolic properties might promote the ongoing gradual replacement of the current epidemic LA-MRSA-CC9 clone in China with the emerging LA-MRSA-CC398 lineage through livestock trade and occupational exposure. Due to the enhanced pathogenicity of the LA-MRSA-CC398 clone, the public health risk posed by LA-MRSA from swine might increase further.

Development and Evaluation of a Core Genome Multilocus Sequencing Typing (cgMLST) Scheme for Serratia marcescens Molecular Surveillance and Outbreak Investigations

Serratia marcescens can cause a range of severe infections and contributes to nosocomial outbreaks. Although whole-genome sequencing (WGS)-based typing is the standard method for molecular surveillance and outbreak investigation, there is no standardized analytic scheme for S. marcescens core genome multilocus sequence typing (cgMLST). Here, the development and evaluation of a S. marcescens cgMLST scheme is reported with the goal of enabling a standardized methodology and typing nomenclature. Four hundred ninety-one high-quality S. marcescens WGS data sets were extracted from public databases and-using the genomic sequence of NCBI reference strain S. marcescens Db11 (NZ_HG326223.1) as a starting point-all Db11 genes present in ≥97% data sets used to create a cgMLST scheme. The novel scheme was evaluated using WGS data from 24 outbreak investigations (n = 175 isolates) distributed over three continents. Analysis of Db11 genes within the 491 data sets identified 2,692 target genes present in ≥97% of genomes (mean, 99.1%; median, 99.9%). These genes formed the novel cgMLST scheme, covering 47.8% of nucleotides in the Db11 genome. Analyzing 175 isolates from 24 outbreaks using the novel scheme gave comparable results to previous typing efforts for both general groupings and allelic distances within clusters. In summary, a novel cgMLST scheme for S. marcescens was developed and evaluated. The scheme and its associated nomenclature will improve standardization of typing efforts for molecular surveillance and outbreak investigation, allowing better understanding of S. marcescens genomic epidemiology and facilitating interlaboratory comparisons.

The power, potential, benefits, and challenges of implementing high-throughput sequencing in food safety systems

The development and application of modern sequencing technologies have led to many new improvements in food safety and public health. With unprecedented resolution and big data, high-throughput sequencing (HTS) has enabled food safety specialists to sequence marker genes, whole genomes, and transcriptomes of microorganisms almost in real-time. These data reveal not only the identity of a pathogen or an organism of interest in the food supply but its virulence potential and functional characteristics. HTS of amplicons, allow better characterization of the microbial communities associated with food and the environment. New and powerful bioinformatics tools, algorithms, and machine learning allow for development of new models to predict and tackle important events such as foodborne disease outbreaks. Despite its potential, the integration of HTS into current food safety systems is far from complete. Government agencies have embraced this new technology, and use it for disease diagnostics, food safety inspections, and outbreak investigations. However, adoption and application of HTS by the food industry have been comparatively slow, sporadic, and fragmented. Incorporation of HTS by food manufacturers in their food safety programs could reinforce the design and verification of effectiveness of control measures by providing greater insight into the characteristics, origin, relatedness, and evolution of microorganisms in our foods and environment. Here, we discuss this new technology, its power, and potential. A brief history of implementation by public health agencies is presented, as are the benefits and challenges for the food industry, and its future in the context of food safety.

Evaluation of within-host evolution of methicillin-resistant Staphylococcus aureus (MRSA) by comparing cgMLST and SNP analysis approaches

Whole genome sequencing (WGS) of methicillin-resistant Staphylococcus aureus (MRSA) provides high-resolution typing, facilitating surveillance and outbreak investigations. The aim of this study was to evaluate the genomic variation rate in MRSA, by comparing commonly used core genome multilocus sequencing (cgMLST) against single nucleotide polymorphism (SNP) analyses. WGS was performed on 95 MRSA isolates, collected from 20 carriers during years 2003–2019. To assess variation and methodological-related differences, two different cgMLST schemes were obtained using Ridom SeqSphere+ and the cloud-based 1928 platform. In addition, two SNP methods, 1928 platform and Northern Arizona SNP Pipeline (NASP) were used. The cgMLST using Ridom SeqSphere+ and 1928 showed a median of 5.0 and 2.0 allele variants/year, respectively. In the SNP analysis, performed with two reference genomes COL and Newman, 1928 showed a median of 13 and 24 SNPs (including presumed recombination) and 3.8 respectively 4.0 SNPs (without recombination) per individual/year. Accordantly, NASP showed a median of 5.5 and 5.8 SNPs per individual/year. In conclusion, an estimated genomic variation rate of 2.0–5.8 genetic events per year (without recombination), is suggested as a general guideline to be used at clinical laboratories for surveillance and outbreak investigations independently of analysis approach used.

Genomic Investigation Reveals a Community Typhoid Outbreak Caused by Contaminated Drinking Water in China, 2016

Typhoid fever is a life-threatening disease caused by Salmonella enterica serovar Typhi (S. Typhi) and remains a significant public health burden in developing countries. In China, typhoid fever is endemic with a limited number of reported outbreaks. Recently, Chinese local Center for Disease Prevention and Control is starting to apply whole genome sequencing for tracking the source of outbreak isolates. In this study, we conducted a retrospective investigation into a community outbreak of typhoid fever in Lanling, China, in 2016. A total of 26 S. Typhi isolates were recovered from the drinking water (n = 1) and patients' blood (n = 24) and stool (n = 1). Phylogenetic analysis indicated the persistence of the outbreak isolates in drinking water for more than 3 months. The genomic comparison demonstrated a high similarity between the isolate from water and isolates from patients in their genomic content, virulence gene profiles, and antimicrobial resistance gene profile, indicating the S. Typhi isolate from drinking water was responsible for the examined outbreak. The result of pulsed-field gel electrophoresis (PFGE) revealed these isolates had identical PFGE pattern, indicating they are clonal variants. Additionally, phylogeographical analysis of global S. Typhi isolates suggested the outbreak isolates are evolutionarily linked to the isolates from the United Kingdom and Vietnam. Taken together, this study highlights the drinking water and international travel as critical control points of mitigating the outbreak, emphasizing the necessity of regular monitoring of this pathogen in China.

Genetic Characterization of Staphylococcus aureus From Subclinical Mastitis Cases in Dairy Cows in Rwanda

Whole-genome sequencing was carried out on 30 Staphylococcus (S.) aureus isolates from dairy cows with subclinical mastitis from all five provinces of Rwanda. Twenty-five of the isolates produced enough sequence to be analyzed using core genome multilocus sequence typing (cg-MLST). The isolates group into three main clusters. The largest cluster contain isolates of sequence type (ST) 152 (n = 6) and the closely related ST1633 (n = 2). These sequence types have previously mainly been encountered in humans. The isolates of the second-largest cluster belong to ST5477 (n = 5),so far exclusively isolated from cows in Rwanda. The third cluster consists of isolates of ST97 (n = 4), which is a well-known bovine-adapted sequence type. These three clusters were all widespread over the country. Isolates of the usually human-adapted sequence types 1 (n = 2) and 5 (n= 1) were found and a single isolate of ST2430, previously found among humans in Africa. Finally, four isolates of novel sequence types were found: ST7108 (n = 2), ST7109 (n = 1), and ST7110 (n = 1). The blaZ penicillin resistance gene was found in 84% of the isolates and was in all cases corroborated by phenotypic resistance determination. Five (20%) of the isolates carried a tetracycline resistance gene, tet(K) or tetM, and three of these five also displayed phenotypic resistance while two isolates carried a tetM-gene but were yet tetracycline susceptible. Seven (28%) isolates carried the dfrG gene conferring resistance to trimethoprim. Four of these isolates indeed were resistant to trimethoprim while three isolates were sensitive. The str gene conferring resistance to aminoglycosides was found in three isolates; however, none of these displayed resistance to gentamycin. Our data revealed a high diversity of the sequence types of S. aureus isolates from cows with subclinical mastitis in Rwanda. Two major clusters of ST97 and ST5477 are likely to be bovine adapted and cause mastitis while the third cluster of ST152 usually have been found in humans and may signify a recent transmission of these types from human to cows, for example from hand milking. The high prevalence of this sequence type among dairy cows may pose zoonotic threat. The sequence types were widely distributed without any geographic correlation. Penicillin resistance, the most common type of resistance with a prevalence over 80%, but also tetracycline and trimethoprim resistance were displayed by several isolates.

Whole Genome Sequencing Provides an Added Value to the Investigation of Staphylococcal Food Poisoning Outbreaks

Through staphylococcal enterotoxin (SE) production, Staphylococcus aureus is a common cause of food poisoning. Detection of staphylococcal food poisoning (SFP) is mostly performed using immunoassays, which, however, only detect five of 27 SEs described to date. Polymerase chain reactions are, therefore, frequently used in complement to identify a bigger arsenal of SE at the gene level (se) but are labor-intensive. Complete se profiling of isolates from different sources, i.e., food and human cases, is, however, important to provide an indication of their potential link within foodborne outbreak investigation. In addition to complete se gene profiling, relatedness between isolates is determined with more certainty using pulsed-field gel electrophoresis, Staphylococcus protein A gene typing and other methods, but these are shown to lack resolution. We evaluated how whole genome sequencing (WGS) can offer a solution to these shortcomings. By WGS analysis of a selection of S. aureus isolates, including some belonging to a confirmed foodborne outbreak, its added value as the ultimate multiplexing method was demonstrated. In contrast to PCR-based se gene detection for which primers are sometimes shown to be non-specific, WGS enabled complete se gene profiling with high performance, provided that a database containing reference sequences for all se genes was constructed and employed. The custom compiled database and applied parameters were made publicly available in an online user-friendly interface. As an all-in-one approach with high resolution, WGS additionally allowed inferring correct isolate relationships. The different DNA extraction kits that were tested affected neither se gene profiling nor relatedness determination, which is interesting for data sharing during SFP outbreak investigation. Although confirming the production of enterotoxins remains important for SFP investigation, we delivered a proof-of-concept that WGS is a valid alternative and/or complementary tool for outbreak investigation.

Comparison of Whole Genome Sequencing and Repetitive Element PCR for Multidrug- Resistant Pseudomonas aeruginosa Strain Typing

Hospital-acquired infections pose significant costly global challenges to patient care. Rapid and sensitive methods to identify potential outbreaks are integral to infection control measures. Whole genome sequencing (WGS)-based bacterial strain typing provides higher discriminatory power over standard nucleotide banding pattern-based methods like repetitive sequence-based PCR (rep-PCR). However, integration of WGS into clinical epidemiology is limited by the lack of consensus in methodology and data analysis/interpretation. In this study, WGS was performed on 22 multidrug-resistant Pseudomonas aeruginosa (MDR-PA) genomic DNA using the Illumina MiSeq platform. Resulting high quality reads were analyzed for phylogenetic relatedness using whole genome multi-locus sequence typing (wgMLST)-based BIOMÉRIEUX EPISEQ CS (EpiSeq; Durham, NC) and Single Nucleotide Variant PHYLogenomics (SNVPhyl). WGS-based results were compared to conventional MLST and archived rep-PCR results. Rep-PCR identified three independent clonal clusters of MDR-PA. Only one clonal cluster identified by rep-PCR, an endemic strain within the pediatric cystic fibrosis population at Texas Children's Hospital, was concordantly identified using wgMLST (EpiSeq) and SNVPhyl. Results were highly consistent between the three sequence-based analyses (conventional MLST, wgMLST, and SNVPhyl), and these results remained consistent with the addition of 74 MDR-PA genomes. These WGS-based methods provided greater resolution for strain discrimination than rep-PCR or standard MLST classification, and the ease of use of EpiSeq renders it clinically viable for analysis, interpretation, and reporting of WGS-based strain typing.

A Combined Phenotypic-Genotypic Predictive Algorithm for In Vitro Detection of Bicarbonate: β-Lactam Sensitization among Methicillin-Resistant Staphylococcus aureus (MRSA)

Antimicrobial susceptibility testing (AST) is routinely used to establish predictive antibiotic resistance metrics to guide the treatment of bacterial pathogens. Recently, a novel phenotype termed "bicarbonate (NaHCO₃)-responsiveness" was identified in a relatively high frequency of clinical MRSA strains, wherein isolates demonstrate in vitro "susceptibility" to standard β-lactams (oxacillin [OXA]; cefazolin [CFZ]) in the presence of NaHCO₃, and in vivo susceptibility to these β-lactams in experimental endocarditis models. We investigated whether a targeted phenotypic-genotypic screening of MRSA could rule in or rule out NaHCO₃ susceptibility upfront. We studied 30 well-characterized clinical MRSA bloodstream isolates, including 15 MIC-susceptible to CFZ and OXA in NaHCO₃-supplemented Mueller-Hinton Broth (MHB); and 15 MIC-resistant to both β-lactams in this media. Using a two-tiered strategy, isolates were first screened by standard disk diffusion for susceptibility to a combination of amoxicillin-clavulanate [AMC]. Isolates then underwent genomic sequence typing: MLST (clonal complex [CC]); agr; SCCmec; and mecA promoter and coding region. The combination of AMC disk susceptibility testing plus mecA and spa genotyping was able to predict MRSA strains that were more or less likely to be NaHCO₃-responsive in vitro, with a high degree of sensitivity and specificity. Validation of this screening algorithm was performed in six strains from the overall cohort using an ex vivo model of endocarditis. This ex vivo model recapitulated the in vitro predictions of NaHCO₃-responsiveness vs. nonresponsiveness above in five of the six strains.

Comparative Transcriptomic Analysis of Staphylococcus aureus Associated with Periprosthetic Joint Infection under in Vivo and in Vitro Conditions

Transcriptomic analysis can provide insight as to how Staphylococcus aureus adapts to the environmental niche of periprosthetic joint infection (PJI), a challenging clinical infection. Here, in vivo RNA expression of eight S. aureus PJIs was compared with expression of the corresponding isolates in planktonic culture using a total RNA-sequencing approach. Expression varied among isolates, with a common trend showing increased expression of several ica-independent biofilm formation genes, including sdr, fnb, ebpS, and aaa; genes encoding enzymes and toxins, including coa, nuc, hlb, and hlgA/B/C; and genes facilitating acquisition of iron via the iron-binding molecule siderophore B (snb) and heme consumption protein (isd) pathways in PJI. Several antimicrobial resistance determinants were detected; although their presence correlated with phenotypic susceptibility of the associated isolates, no difference in expression between in vivo and in vitro conditions was identified.

Hospital transmission of borderline oxacillin-resistant Staphylococcus aureus evaluated by whole-genome sequencing

Introduction

Staphylococcus aureus is a major cause of hospital infections worldwide. Awareness towards methicillin-resistant S. aureus (MRSA) infections is high but attention towards borderline oxacillin-resistant S. aureus (BORSA) is limited, possibly due to an underestimated clinical relevance, presumption of low incidence and diagnostic limitations.

Gap statement

BORSA surveillance has not been routinely implemented, and thus consensus with regard to a definition and infection control measures is lacking.

Aim

Our goals were to investigate the occurrence, molecular characteristics and clinical manifestations of BORSA infections in the hospital setting.

Methodology

Following an increased incidence in 2016, BORSA cases in 2014/2016 (in our institution) were more specifically evaluated. Medical records were reviewed to investigate epidemiological links, clinical characteristics and outcomes. Resistance and virulence markers were assessed by whole genome sequencing (WGS). Conventional methods: amplified fragment length polymorphism (AFLP) ; multilocus sequence typing (MLST) and multiple locus variable-number tandem repeat analysis (MLVA) were compared with core genome MLST (cgMLST) and whole-genome single nucleotide polymorphism (wgSNP) analysis to confirm genetic clusters.

Results

From 2009 to 2013, BORSA comprised 0.1 % of all clinical S. aureus strains. In 2016, the incidence was six-fold higher in comparison to the baseline. Whole-genome SNP and cgMLST confirmed two BORSA clusters among patients with dermatological conditions. Patients with BORSA presented with skin infections, and one case developed a severe invasive infection with a fatal outcome. Infection control measures successfully prevented further transmission in both clusters. WGS findings showed that BORSA strains carried multiple resistance and virulence genes with increased pathogenic potential.

Conclusion

WGS and cgMLST effectively characterized and confirmed BORSA clusters among at-risk patients with clinical manifestations ranging from mild skin infections to life-threatening bacteraemia. Clinical awareness and active monitoring are therefore warranted for the timely implementation of infection control measures to prevent BORSA transmission in high-risk patients.

Household Transmission of Community-Associated Methicillin-Resistant Staphylococcus Aureus

Currently, the mechanism of community-associated methicillin-resistant Staphylococcus aureus (CA-MRSA) transmission mechanism is unclear; however, it must be considered in conjunction with asymptomatic S. aureus strains colonization dynamics. This epidemiological study aimed to determine the role of the household in CA-MRSA transmission in China. Five patients with culture-confirmed CA-MRSA infection and five control patients were recruited from the Sir Run Run Shaw Hospital in Zhejiang, China, between December 2019 and January 2020. The household members of the patients, their pets, and environmental surfaces were sampled and screened for MRSA colonization. Mass spectrometry identification and antimicrobial susceptibility testing were performed on the MRSA isolates. Whole-genome sequencing and core genome multilocus sequence typing (cgMLST) were performed to determine the origin and transmission of the MRSA isolates in the households. Overall, 14 S. aureus-positive specimens (14.1%, 14/99) were obtained from the five households of patients with CA-MRSA infections, of which 12 (85.7%) were MRSA. The overall positivity of MRSA was 12.1% (12/99) among the samples from the CA-MRSA households, while no MRSA isolates were detected in the five control households. Most MRSA isolates belonged to epidemic CA-MRSA clones, such as ST59 (15/35, 42.9%) and ST508 (15/35, 42.9%). The cgMLST results confirmed that MRSA was transmitted among patients, contacts, and pets in the households and was present on environmental surfaces in the CA-MRSA patients' households. In conclusion, the study revealed that the home environment was an important MRSA reservoir. Therefore, focusing on MRSA decolonization in patients alone is not sufficient for infection control of CA-MRSA.

Genetic changes are introduced by repeated exposure of Salmonella spiked in low water activity and high fat matrix to heat

WGS is used to define if isolates are "in" or "out" of an outbreak and/or microbial root cause investigation. No threshold of genetic differences is fixed and the conclusions on similarity between isolates are mainly based on the knowledge generated from previous outbreak investigations and reported mutation rates. Mutation rates in Salmonella when exposed to food processing conditions are lacking. Thus, in this study, the ability of heat and dry stress to cause genetic changes in two Salmonella serotypes frequently isolated from low moisture foods was investigated. S. enterica serovars S. Agona ATCC 51,957 and S. Mbandaka NCTC 7892 (ATCC 51,958) were repeatedly exposed to heat (90 °C for 5 min) in a low water activity and high fat matrix. No increased fitness of the strains was observed after 10 repeated heat treatments. However, genetic changes were introduced and the number of genetic differences increased with every heat treatment cycle. The genetic changes appeared randomly in the genome and were responsible for a population of diverse isolates with 0 to 28 allelic differences (0 to 38 SNPs) between them. This knowledge is key to interpret WGS results for source tracking investigations as part of a root cause analysis in a contamination event as isolates are exposed to stress conditions.

Rapid Open-Source SNP-Based Clustering Offers an Alternative to Core Genome MLST for Outbreak Tracing in a Hospital Setting

Traditional genotyping methods for infection control of antimicrobial-resistant bacteria in healthcare settings have been supplemented by whole-genome sequencing (WGS), often relying on a gene-based approach, e.g., core genome multilocus sequence typing (cgMLST), to cluster-related samples. In this study, we compared clusters of methicillin-resistant Staphylococcus aureus (MRSA) and Enterococcus faecium analyzed with the commercial cgMLST software Ridom SeqSphere+ and with an open-source single-nucleotide polymorphism (SNP)-based phylogenetic analysis pipeline (PAPABAC). A total of 5,655 MRSA and 2,572 E. faecium patient isolates, collected between 2013 and 2018, were processed. Clusters of 1,844 MRSA and 1,355 E. faecium isolates were compared to cgMLST results, and epidemiological data were included when available. The phylogenies inferred by the two different technologies were highly concordant, and the MRSA SNP tree re-captured known hospital-related outbreaks and epidemiologically linked samples. PAPABAC has the advantage over Ridom SeqSphere+ to generate stable, referable clusters without the need for sequence assembly, and it is a free-of-charge, open-source alternative to the commercial software.

Rapid nanopore-based DNA sequencing protocol of antibiotic-resistant bacteria for use in surveillance and outbreak investigation

Outbreak investigations are essential to control and prevent the dissemination of pathogens. This study developed and validated a complete analysis protocol for faster and more accurate surveillance and outbreak investigations of antibiotic-resistant microbes based on Oxford Nanopore Technologies (ONT) DNA whole-genome sequencing. The protocol was developed using 42 methicillin-resistant Staphylococcus aureus (MRSA) isolates identified from former well-characterized outbreaks. The validation of the protocol was performed using Illumina technology (MiSeq, Illumina). Additionally, a real-time outbreak investigation of six clinical S. aureus isolates was conducted to test the ONT-based protocol. The suggested protocol includes: (1) a 20 h sequencing run; (2) identification of the sequence type (ST); (3) de novo genome assembly; (4) polishing of the draft genomes; and (5) phylogenetic analysis based on SNPs. After the sequencing run, it was possible to identify the ST in 2 h (20 min per isolate). Assemblies were achieved after 4 h (40 min per isolate) while the polishing was carried out in 7 min per isolate (42 min in total). The phylogenetic analysis took 0.6 h to confirm an outbreak. Overall, the developed protocol was able to at least discard an outbreak in 27 h (mean) after the bacterial identification and less than 33 h to confirm it. All these estimated times were calculated considering the average time for six MRSA isolates per sequencing run. During the real-time S. aureus outbreak investigation, the protocol was able to identify two outbreaks in less than 31 h. The suggested protocol enables identification of outbreaks in early stages using a portable and low-cost device along with a streamlined downstream analysis, therefore having the potential to be incorporated in routine surveillance analysis workflows. In addition, further analysis may include identification of virulence and antibiotic resistance genes for improved pathogen characterization.

A Bioinformatic Pipeline for Improved Genome Analysis and Clustering of Isolates during Outbreaks of Legionnaires' Disease

Legionnaires' disease, a severe lung infection caused by the bacterium Legionella pneumophila, occurs as single cases or in outbreaks that are actively tracked by public health departments. To determine the point source of an outbreak, clinical isolates need to be compared to environmental samples to find matching isolates. One confounding factor is the genome plasticity of L. pneumophila, making an exact sequence comparison by whole-genome sequencing (WGS) challenging. Here, we present a WGS analysis pipeline, LegioCluster, that is designed to circumvent this problem by automatically selecting the best matching reference genome prior to mapping and variant calling. This approach reduces the number of false-positive variant calls, maximizes the fraction of all genomes that are being compared, and naturally clusters the isolates according to their reference strain. Isolates that are too distant from any genome in the database are added to the list of candidate references, thereby creating a new cluster. Short insertions or deletions are considered in addition to single-nucleotide polymorphisms for increased discriminatory power. This manuscript describes the use of this automated and "locked down" bioinformatic pipeline deployed at the New York State Department of Health's Wadsworth Center for investigating relatedness between clinical and environmental isolates. A similar pipeline has not been widely available for use to support these critically important public health investigations.

Staphylococcus aureus whole genome sequence-based susceptibility and resistance prediction using a clinically amenable workflow

We used graphical user interface-based automated analytical tools from Next Gen Diagnostics (Mountain View, CA) and 1928 Diagnostics (Gothenburg, Sweden) to analyze whole genome sequence (WGS) data from 102 unique blood culture isolates of Staphylococcus aureus to predict antimicrobial susceptibly, with results compared to those of phenotypic susceptibility testing. Of 916 isolate/antibiotic combinations analyzed using the Next Gen Diagnostics tool, there were 9 discrepancies between WGS predictions and phenotypic susceptibility/resistance, including 8 for clindamycin and 1 for minocycline. Of 612 isolate/antibiotic combinations analyzed using the 1928 Diagnostics tool, there were 13 discrepancies between WGS predictions and phenotypic susceptibility/resistance, including 9 for clindamycin, 3 for trimethoprim-sulfamethoxazole, and 1 for rifampin. Trimethoprim-sulfamethoxazole was not assessed by Next Gen Diagnostics, and minocycline was not assessed by 1928 Diagnostics. There was complete concordance between phenotypic susceptibility/resistance and genotypic prediction of susceptibility/resistance using both analytical platforms for oxacillin, vancomycin, and mupirocin, as well as by the Next Gen Diagnostics analytical tool for levofloxacin (the 1928 Diagnostics tool did not assess levofloxacin). These results suggest that, from a performance standpoint, with some caveats, automatic bioinformatics tools may be acceptable to predict susceptibility and resistance to a panel of antibiotics for S. aureus.

Comparison of core-genome MLST, coreSNP and PFGE methods for Klebsiella pneumoniae cluster analysis

In this work we compared the most frequently used Klebsiella pneumoniae typing methods: PFGE, cgMLST and coreSNP. We evaluated the discriminatory power of the three methods to confirm or exclude nosocomial transmission on K. pneumoniae strains isolated from January to December 2017, in the framework of the routine surveillance for multidrug-resistant organisms at the San Raffaele Hospital, in Milan. We compared the results of the different methods to the results of epidemiological investigation. Our results showed that cgMLST and coreSNP are more discriminant than PFGE, and that both approaches are suitable for transmission analyses. cgMLST appeared to be inferior to coreSNP in the K. pneumoniae CG258 phylogenetic reconstruction. Indeed, we found that the phylogenetic reconstruction based on cgMLST genes wrongly clustered ST258 clade1 and clade2 strains, conversely properly assigned by coreSNP approach. In conclusion, this study provides evidences supporting the reliability of both cgMLST and coreSNP for hospital surveillance programs and highlights the limits of cgMLST scheme genes for phylogenetic reconstructions.

Fecal non-aureus Staphylococci are a potential cause of bovine intramammary infection

The presence of non-aureus staphylococci (NAS) in bovine rectal feces has recently been described. Similar to other mastitis causing pathogens, shedding of NAS in the environment could result in intramammary infection. The objective of this study was to investigate whether NAS strains present in feces can cause intramammary infection, likely via teat apex colonization. During a cross-sectional study in 5 dairy herds, samples were collected from the habitats quarter milk, teat apices, and rectal feces from 25%, 10%, and 25% of the lactating cows, respectively, with a cow serving as the source of one type of sample only. Samples from clinical mastitis cases were continuously collected during the 1-year study period as well. The 6 most prevalent NAS species, Staphylococcus (S.) chromogenes, S. cohnii, S. devriesei, S. equorum, S. haemolyticus, and S. hominis, were further subtyped by random amplification of polymorphic deoxyribonucleic acid polymerase chain reaction (RAPD-PCR), when the same NAS species was present in the same herd in the three habitats. For S. chromogenes, S. cohnii, S. devriesei, and S. haemolyticus, the same RAPD type was found in rectal feces, teat apices, and quarter milk, indicating that fecal NAS can infect the mammary gland. For S. hominis and S. equorum, we were unable to confirm the presence of the same RAPD types in the three habitats.

Concordance of SNP- and allele-based typing workflows in the context of a large-scale international Salmonella Enteritidis outbreak investigation

A large European multi-country Salmonella enterica serovar Enteritidis outbreak associated with Polish eggs was characterized by whole-genome sequencing (WGS)-based analysis, with various European institutes using different analysis workflows to identify isolates potentially related to the outbreak. The objective of our study was to compare the output of six of these different typing workflows (distance matrices of either SNP-based or allele-based workflows) in terms of cluster detection and concordance. To this end, we analysed a set of 180 isolates coming from confirmed and probable outbreak cases, which were representative of the genetic variation within the outbreak, supplemented with 22 unrelated contemporaneous S. enterica serovar Enteritidis isolates. Since the definition of a cluster cut-off based on genetic distance requires prior knowledge on the evolutionary processes that govern the bacterial populations in question, we used a variety of hierarchical clustering methods (single, average and complete) and selected the optimal number of clusters based on the consensus of the silhouette, Dunn2, and McClain-Rao internal validation indices. External validation was done by calculating the concordance with the WGS-based case definition (SNP-address) for this outbreak using the Fowlkes-Mallows index. Our analysis indicates that with complete-linkage hierarchical clustering combined with the optimal number of clusters, as defined by three internal validity indices, the six different allele- and SNP-based typing workflows generate clusters with similar compositions. Furthermore, we show that even in the absence of coordinated typing procedures, but by using an unsupervised machine learning methodology for cluster delineation, the various workflows that are currently in use by six European public-health authorities can identify concordant clusters of genetically related S. enterica serovar Enteritidis isolates; thus, providing public-health researchers with comparable tools for detection of infectious-disease outbreaks.

Core genome MLST and resistome analysis of Klebsiella pneumoniae using a clinically amenable workflow

Whole genome sequencing (WGS) is replacing traditional microbiological typing methods for investigation of outbreaks in clinical settings. Here, we used a clinical microbiology laboratory core genome multilocus sequence typing (cgMLST) workflow to analyze 40 isolates of K. pneumoniae which are part of the Antimicrobial Resistance Leadership Group (ARLG) isolate collection, alongside 10 Mayo Clinic K. pneumoniae isolates, comparing results to those of pulsed-field gel electrophoresis (PFGE). Additionally, we used the WGS data to predict phenotypic antimicrobial susceptibility (AST). Thirty-one of 40 ARLG K. pneumoniae isolates belonged to the same PFGE type, all of which, alongside 3 isolates of different PFGE types, formed a large cluster by cgMLST. PFGE and cgMLST were completely concordant for the 10 Mayo Clinic K. pneumoniae isolates. For AST prediction, the overall agreement between phenotypic AST and genotypic prediction was 95.6%.

Pulsed-Field Gel Electrophoresis Typing of Staphylococcus aureus Strains

Pulsed-field gel electrophoresis (PFGE) is one of the most effective genotyping methods for the study of genetic relatedness of microorganisms as well as for epidemiologic investigation of outbreaks caused by pathogenic bacteria, including Staphylococcus aureus. The technique relies on the size patterns of large DNA fragments generated by the cleavage of intact bacterial chromosomes with a rare cutting restriction enzyme, and subsequently resolved by pulsed-field electrophoresis with periodic changes of the orientation of the electrical field across the gel. The high discriminatory power, improved reproducibility by standardization of experimental protocols and data analysis guidelines, and establishment of a national PFGE database of S. aureus profiles have made it a valuable means for global tracking of S. aureus infection sources and determination of genetic relatedness of the strains.

Tracking the environmental dissemination of carbapenem-resistant Klebsiella pneumoniae using whole genome sequencing

The emergence and dissemination of infections caused by carbapenem-resistant Klebsiella pneumoniae (CRKP) are of great concern worldwide, as there are limited options for their treatment. Thus, in this study, whole-genome sequencing (WGS) was applied to assess CRKP distribution and dissemination from hospital settings to the aquatic environment in order to identify the extent of the problem. Samples were collected from hospital wastewaters and receiving water bodies. Susceptible K. pneumoniae and CRKP were enumerated and isolated using standard methods. Seventeen CRKP were DNA-sequenced using an Illumina HiSeq X™ platform. De novo assembly and annotation were performed using SPAdes and RAST, respectively. The study analysed antibiotic resistance traits (antibiotic resistant genes, mobile genetic elements, and virulence genes) in CRKP isolates. Although influent of wastewater harboured the highest CRKP, wastewater treatment plants were efficient in reducing the threat. In terms of resistance per matrix, benthic sediment proved to harbour more CRKP (22.88%) versus susceptible K. pneumoniae, as revealed by their resistant quotient analysis, while effluent of wastewaters (4.21%) and water bodies (4.64%) had the lowest CRKP loads. The disseminating CRKP consisted of six sequence types (ST) - ST307 (n = 7), a novel ST3559 (n = 5), ST15 (n = 2), and one isolate of each of ST39, 152 and 298. All CRKP isolates harboured β-lactams (bla_CTX-M-15 and bla_OXA-1), quinolone (oqxA and oqxB) and fosfomycin (fosA) resistance genes as well as virulence genes. This study highlights the dissemination of 'high' importance and novel ST CRKP from hospital wastewater to waterbodies. This is concerning, particularly in the African context where a sizable number of people still rely on direct water resources for household use, including drinking. Further research is needed to systematically track the occurrence and distribution of these bacteria so as to mitigate their threat.

Using Core-genome Multilocus Sequence Typing to Monitor the Changing Epidemiology of Methicillin-resistant Staphylococcus aureus in a Teaching Hospital

Background

This study was performed to elucidate the changing epidemiology of methicillin-resistant Staphylococcus aureus (MRSA) in a Chinese teaching hospital.

Methods

We performed whole-genome sequencing and phylogenetic analysis on MRSA isolates recovered at Sir Run Run Shaw Hospital (SRRSH) from 2013 to 2015. MRSA isolates underwent susceptibility testing, resistance determinant screening, virulence assays and clinical data analysis.

Results

A total of 292 MRSA infections were identified during the study period. Of these, most belonged to sequence type (ST) 5 (51.4%; 150 of 292) and ST59 (23.3%; 68 of 292). The proportion of ST5 MRSA decreased from 68.3% in 2013 to 32.1% in 2015, and the proportion of ST59 MRSA increased from 8.9% to 41.0%. Core-genome phylogenetic analysis revealed that ST59 MRSA isolates exhibited wider genotypic diversity than ST5 MRSA isolates in both healthcare-onset and community-onset groups. Minimum spanning trees showed that a cluster of ST5 MRSA was circulating at SRRSH, whereas the ST59 MRSA infections were mostly sporadic. ST59 MRSA caused 45.2% of the 93 skin and soft-tissue infection cases and displayed lower levofloxacin (11.8%) and ciprofloxacin (19.1%) resistance rates than the ST239 and ST5 MRSA isolates. ST59 healthcare-onset MRSA displayed enhanced virulence in the skin infection model and hemolysis assays. Notably, these isolates had virulence levels similar to those of classic community-associated MRSA strains.

Conclusions

ST59 MRSA strains with high virulence potential have been replacing ST5 MRSA in predominance in SRRSH and causing nosocomial infections. Whole-genome sequencing is a powerful tool to monitor changes in the epidemiology of MRSA in hospitals.

Molecular Characterization of Equine Staphylococcus aureus Isolates Exhibiting Reduced Oxacillin Susceptibility

The detection of borderline oxacillin-resistant Staphylococcus aureus (BORSA) represents a challenge to both, veterinary and human laboratories. Between 2015 and 2017, 19 equine S. aureus with elevated minimal inhibitory concentrations for oxacillin were detected in routine diagnostics. The aim of this study was to characterize these isolates to identify factors possibly associated with the BORSA phenotype. All S. aureus were subjected to antimicrobial susceptibility testing and whole genome sequencing (WGS). A quantifiable β-lactamase activity assay was performed for a representative subset of 13 isolates. The WGS data analysis of the 19 BORSA isolates identified two different genomic lineages, sequence type (ST) 1 and ST1660. The core genome multilocus sequence typing (cgMLST) revealed a close relatedness of all isolates belonging to either ST1 or ST1660. The WGS analysis identified the resistance genes aadD, dfrG, tet(L), and/or blaZ and aacA-aphD. Phenotypic resistance to penicillins, aminoglycosides, tetracyclines, fluoroquinolones and sulfamethoxazole/trimethoprim was observed in the respective isolates. For the penicillin-binding proteins 1-4, amino acid substitutions were predicted using WGS data. Since neither transglycosylase nor transpeptidase domains were affected, these alterations might not explain the BORSA phenotype. Moreover, β-lactamase activity was found to be associated with an inducible blaZ gene. The lineage-specific differences regarding the expression profiles were noted.

Staphylococcus aureus Infections in Malaysia: A Review of Antimicrobial Resistance and Characteristics of the Clinical Isolates, 1990-2017

Staphylococcus aureus is an important nosocomial pathogen and its multidrug resistant strains, particularly methicillin-resistant S. aureus (MRSA), poses a serious threat to public health due to its limited therapeutic options. The increasing MRSA resistance towards vancomycin, which is the current drug of last resort, gives a great challenge to the treatment and management of MRSA infections. While vancomycin resistance among Malaysian MRSA isolates has yet to be documented, a case of vancomycin resistant S. aureus has been reported in our neighboring country, Indonesia. In this review, we present the antimicrobial resistance profiles of S. aureus clinical isolates in Malaysia with data obtained from the Malaysian National Surveillance on Antimicrobial Resistance (NSAR) reports as well as various peer-reviewed published records spanning a period of nearly three decades (1990-2017). We also review the clonal types and characteristics of Malaysian S. aureus isolates, where hospital-associated (HA) MRSA isolates tend to carry staphylococcal cassette chromosome mec (SCCmec) type III and were of sequence type (ST)239, whereas community-associated (CA) isolates are mostly SCCmec type IV/V and ST30. More comprehensive surveillance data that include molecular epidemiological data would enable further in-depth understanding of Malaysian S. aureus isolates.

Surgical Site Infections Caused by Highly Virulent Methicillin-Resistant Staphylococcus aureus Sequence Type 398, China

We identified 2 methicillin-resistant Staphylococcus aureus strains of sequence type 398 from surgical site infections in China. Genetic analysis and clinical data from these strains suggested that they were human-related but sporadic. Hemolysis analysis and mouse-skin infection models indicated a high virulence potential for these strains.

Use of Whole-Genome Sequencing for Food Safety and Public Health in the United States

Whole-genome sequencing (WGS) is increasingly used by food regulatory and public health agencies in the United States to facilitate the detection, investigation, and control of foodborne bacterial outbreaks, and food regulatory and other activities in support of food safety. WGS has added a level of precision to the surveillance leading to faster and more efficient decision making in the preparedness and response to foodborne infections. In this review, we report the history of WGS technology at the Centers for Disease Control & Prevention (CDC), the Food and Drug Administration (FDA), and the United States Department of Agriculture's Food Safety and Inspection Service (USDA/FSIS) as it applies to food safety. The basic principle of the method, the analysis, and interpretation of the data are explained as is its major strengths and limitations. We also describe the benefits and possibilities of the WGS technology to the food industry throughout the farm-to-fork continuum and the prospects of metagenomic sequencing applied directly to the sample specimen with or without pre-enrichment culture.

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.

Yersinia kristensenii subsp. rochesterensis subsp. nov., isolated from human feces

A single bacterial isolate, EPLC-04^T, was isolated from human feces and identified as representing a member of the genus Yersinia on the basis of phenotypic characteristics, matrix assisted laser desorption ionization time-of-flight mass spectrometry and partial 16S rRNA gene sequencing. The isolate's phenotypic profile differed from that described for the most closely related species, Yersinia kristensenii, by exhibiting lipase production and lacking pyrazinamidase activity. Multiple genetic targets, including the complete (1465 bp) 16S rRNA gene sequence and partial sequences of groEL (539 bp), gyrB (935 bp), glnA (525 bp) and recA (535 bp) indicated that the isolate exhibited 98.91, 92.16, 90.81, 92.78 and 89.01 % identity with Yersinia aldovae, 98.98, 91.99, 90.17, 89.77 and 89.55 % identity with Yersinia intermedia, and 99.66, 98.11, 98.50, 98.49 and 98.51 % identity with Y. kristensenii, respectively. Phylogenetic reconstructions based on the combination of the four housekeeping genes indicated that the isolate formed a unique branch, supported by a bootstrap value of 100 %. Digital DNA-DNA homology and 16S rRNA gene sequencing identified EPLC-04^T as representing Y. kristensenii. However, the unique phenotypic traits and results of phylogenetic analysis indicate that it represents a novel subspecies of Y. kristensenii. The name Yersinia kristenseniisubsp. rochesterensis subsp. nov. is proposed for this novel taxon (type strain EPLC-04^T=ATCC BAA-2637^T, DSMZ 28595^T).

The Continuing Threat of Methicillin-Resistant Staphylococcus aureus

Staphylococcus aureus has been an exceptionally successful pathogen, which is still relevant in modern age-medicine due to its adaptability and tenacity. This bacterium may be a causative agent in a plethora of infections, owing to its abundance (in the environment and in the normal flora) and the variety of virulence factors that it possesses. Methicillin-resistant S. aureus (MRSA) strains—first described in 1961—are characterized by an altered penicillin-binding protein (PBP2a/c) and resistance to all penicillins, cephalosporins, and carbapenems, which makes the β-lactam armamentarium clinically ineffective. The acquisition of additional resistance determinants further complicates their eradication; therefore, MRSA can be considered as the first representative of multidrug-resistant bacteria. Based on 230 references, the aim of this review is to recap the history, the emergence, and clinical features of various MRSA infections (hospital-, community-, and livestock-associated), and to summarize the current advances regarding MRSA screening, typing, and therapeutic options (including lipoglycopeptides, oxazolidinones, anti-MRSA cephalosporins, novel pleuromutilin-, tetracycline- and quinolone-derivatives, daptomycin, fusidic acid, in addition to drug candidates in the development phase), both for an audience of clinical microbiologists and infectious disease specialists.

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.

Whole-genome sequencing for methicillin-resistant Staphylococcus aureus (MRSA) outbreak investigation in a neonatal intensive care unit

OBJECTIVE

To evaluate whole-genome sequencing (WGS) as a molecular typing tool for MRSA outbreak investigation.

DESIGN

Investigation of MRSA colonization/infection in a neonatal intensive care unit (NICU) over 3 years (2014-2017).

SETTING

Single-center level IV NICU.PatientsNICU infants and healthcare workers (HCWs).

METHODS

Infants were screened for MRSA using a swab of the anterior nares, axilla, and groin, initially by targeted (ring) screening, and later by universal weekly screening. Clinical cultures were collected as indicated. HCWs were screened once using swabs of the anterior nares. MRSA isolates were typed using WGS with core-genome multilocus sequence typing (cgMLST) analysis and by pulsed-field gel electrophoresis (PFGE). Colonized and infected infants and HCWs were decolonized. Control strategies included reinforcement of hand hygiene, use of contact precautions, cohorting, enhanced environmental cleaning, and remodeling of the NICU.

RESULTS

We identified 64 MRSA-positive infants: 53 (83%) by screening and 11 (17%) by clinical cultures. Of 85 screened HCWs, 5 (6%) were MRSA positive. WGS of MRSA isolates identified 2 large clusters (WGS groups 1 and 2), 1 small cluster (WGS group 3), and 8 unrelated isolates. PFGE failed to distinguish WGS group 2 and 3 isolates. WGS groups 1 and 2 were codistributed over time. HCW MRSA isolates were primarily in WGS group 1. New infant MRSA cases declined after implementation of the control interventions.

CONCLUSION

We identified 2 contemporaneous MRSA outbreaks alongside sporadic cases in a NICU. WGS was used to determine strain relatedness at a higher resolution than PFGE and was useful in guiding efforts to control MRSA transmission.

Routine Whole-Genome Sequencing for Outbreak Investigations of Staphylococcus aureus in a National Reference Center

The French National Reference Center for Staphylococci currently uses DNA arrays and spa typing for the initial epidemiological characterization of Staphylococcus aureus strains. We here describe the use of whole-genome sequencing (WGS) to investigate retrospectively four distinct and virulent S. aureus lineages [clonal complexes (CCs): CC1, CC5, CC8, CC30] involved in hospital and community outbreaks or sporadic infections in France. We used a WGS bioinformatics pipeline based on de novo assembly (reference-free approach), single nucleotide polymorphism analysis, and on the inclusion of epidemiological markers. We examined the phylogeographic diversity of the French dominant hospital-acquired CC8-MRSA (methicillin-resistant S. aureus) Lyon clone through WGS analysis which did not demonstrate evidence of large-scale geographic clustering. We analyzed sporadic cases along with two outbreaks of a CC1-MSSA (methicillin-susceptible S. aureus) clone containing the Panton–Valentine leukocidin (PVL) and results showed that two sporadic cases were closely related. We investigated an outbreak of PVL-positive CC30-MSSA in a school environment and were able to reconstruct the transmission history between eight families. We explored different outbreaks among newborns due to the CC5-MRSA Geraldine clone and we found evidence of an unsuspected link between two otherwise distinct outbreaks. Here, WGS provides the resolving power to disprove transmission events indicated by conventional methods (same sequence type, spa type, toxin profile, and antibiotic resistance profile) and, most importantly, WGS can reveal unsuspected transmission events. Therefore, WGS allows to better describe and understand outbreaks and (inter-)national dissemination of S. aureus lineages. Our findings underscore the importance of adding WGS for (inter-)national surveillance of infections caused by virulent clones of S. aureus but also substantiate the fact that technological optimization at the bioinformatics level is still urgently needed for routine use. However, the greatest limitation of WGS analysis is the completeness and the correctness of the reference database being used and the conversion of floods of data into actionable results. The WGS bioinformatics pipeline (EpiSeq^TM) we used here can easily generate a uniform database and associated metadata for epidemiological applications.

Comparison of SNP-based subtyping workflows for bacterial isolates using WGS data, applied to Salmonella enterica serotype Typhimurium and serotype 1,4,[5],12:i:

Whole genome sequencing represents a promising new technology for subtyping of bacterial pathogens. Besides the technological advances which have pushed the approach forward, the last years have been marked by considerable evolution of the whole genome sequencing data analysis methods. Prior to application of the technology as a routine epidemiological typing tool, however, reliable and efficient data analysis strategies need to be identified among the wide variety of the emerged methodologies. In this work, we have compared three existing SNP-based subtyping workflows using a benchmark dataset of 32 Salmonella enterica subsp. enterica serovar Typhimurium and serovar 1,4,[5],12:i:- isolates including five isolates from a confirmed outbreak and three isolates obtained from the same patient at different time points. The analysis was carried out using the original (high-coverage) and a down-sampled (low-coverage) datasets and two different reference genomes. All three tested workflows, namely CSI Phylogeny-based workflow, CFSAN-based workflow and PHEnix-based workflow, were able to correctly group the confirmed outbreak isolates and isolates from the same patient with all combinations of reference genomes and datasets. However, the workflows differed strongly with respect to the SNP distances between isolates and sensitivity towards sequencing coverage, which could be linked to the specific data analysis strategies used therein. To demonstrate the effect of particular data analysis steps, several modifications of the existing workflows were also tested. This allowed us to propose data analysis schemes most suitable for routine SNP-based subtyping applied to S. Typhimurium and S. 1,4,[5],12:i:-. Results presented in this study illustrate the importance of using correct data analysis strategies and to define benchmark and fine-tune parameters applied within routine data analysis pipelines to obtain optimal results.

Microbiome at the Frontier of Personalized Medicine

The genomic revolution promises to transform our approach to treat patients by individualizing treatments, reducing adverse events, and decreasing health care costs. The early advances using this have been realized primarily by optimizing preventive and therapeutic approaches in cancer using human genome sequencing. The ability to characterize the microbiome, which includes all the microbes that reside within and upon us and all their genetic elements, using next-generation sequencing allows us to now incorporate this important contributor to human disease into developing new preventive and therapeutic strategies. In this review we highlight the importance of the microbiome in all aspects of human disease, including pathogenesis, phenotype, prognosis, and response to treatment, as well as their role as diagnostic and therapeutic biomarkers. We provide a role for next-generation sequencing in both precise microbial identification of infectious diseases and characterization of microbial communities and their function. Taken together, the microbiome is emerging as an integral part of precision medicine approach as it not only contributes to interindividual variability in all aspects of a disease but also represents a potentially modifiable factor that is amenable to targeting by therapeutics.

Molecular epidemiology of Staphylococcus aureus bacteremia in a single large Minnesota medical center in 2015 as assessed using MLST, core genome MLST and spa typing

Staphylococcus aureus is a leading cause of bacteremia in hospitalized patients. Whether or not S. aureus bacteremia (SAB) is associated with clonality, implicating potential nosocomial transmission, has not, however, been investigated. Herein, we examined the epidemiology of SAB using whole genome sequencing (WGS). 152 SAB isolates collected over the course of 2015 at a single large Minnesota medical center were studied. Staphylococcus protein A (spa) typing was performed by PCR/Sanger sequencing; multilocus sequence typing (MLST) and core genome MLST (cgMLST) were determined by WGS. Forty-eight isolates (32%) were methicillin–resistant S. aureus (MRSA). The isolates encompassed 66 spa types, clustered into 11 spa clonal complexes (CCs) and 10 singleton types. 88% of 48 MRSA isolates belonged to spa CC-002 or -008. Methicillin-susceptible S. aureus (MSSA) isolates were more genotypically diverse, with 61% distributed across four spa CCs (CC-002, CC-012, CC-008 and CC-084). By MLST, there was 31 sequence types (STs), including 18 divided into 6 CCs and 13 singleton STs. Amongst MSSA isolates, the common MLST clones were CC5 (23%), CC30 (19%), CC8 (15%) and CC15 (11%). Common MRSA clones were CC5 (67%) and CC8 (25%); there were no MRSA isolates in CC45 or CC30. By cgMLST analysis, there were 9 allelic differences between two isolates, with the remaining 150 isolates differing from each other by over 40 alleles. The two isolates were retroactively epidemiologically linked by medical record review. Overall, cgMLST analysis resulted in higher resolution epidemiological typing than did multilocus sequence or spa typing.