Closing the gap: Oxford Nanopore Technologies R10 sequencing allows comparable results to Illumina sequencing for SNP-based outbreak investigation of bacterial pathogens.

Whole-genome sequencing has become the method of choice for bacterial outbreak investigation, with most clinical and public health laboratories currently routinely using short-read Illumina sequencing. Recently, long-read Oxford Nanopore Technologies (ONT) sequencing has gained prominence and may offer advantages over short-read sequencing, particularly with the recent introduction of the R10 chemistry, which promises much lower error rates than the R9 chemistry. However, limited information is available on its performance for bacterial single-nucleotide polymorphism (SNP)-based outbreak investigation. We present an open-source workflow, Prokaryotic Awesome variant Calling Utility (PACU) (https://github.com/BioinformaticsPlatformWIV-ISP/PACU), for constructing SNP phylogenies using Illumina and/or ONT R9/R10 sequencing data. The workflow was evaluated using outbreak data sets of Shiga toxin-producing Escherichia coli and Listeria monocytogenes by comparing ONT R9 and R10 with Illumina data. The performance of each sequencing technology was evaluated not only separately but also by integrating samples sequenced by different technologies/chemistries into the same phylogenomic analysis. Additionally, the minimum sequencing time required to obtain accurate phylogenetic results using nanopore sequencing was evaluated. PACU allowed accurate identification of outbreak clusters for both species using all technologies/chemistries, but ONT R9 results deviated slightly more from the Illumina results. ONT R10 results showed trends very similar to Illumina, and we found that integrating data sets sequenced by either Illumina or ONT R10 for different isolates into the same analysis produced stable and highly accurate phylogenomic results. The resulting phylogenies for these two outbreaks stabilized after ~20 hours of sequencing for ONT R9 and ~8 hours for ONT R10. This study provides a proof of concept for using ONT R10, either in isolation or in combination with Illumina, for rapid and accurate bacterial SNP-based outbreak investigation.

Whole-Genome Sequencing-Based Screening of MRSA in Patients and Healthcare Workers in Public Hospitals in Benin.

Methicillin-resistant Staphylococcus aureus (MRSA) constitutes a serious public health concern, with a considerable impact on patients' health, and substantial healthcare costs. In this study, patients and healthcare workers (HCWs) from six public hospitals in Benin were screened for MRSA. Strains were identified as MRSA using conventional microbiological methods in Benin, and confirmed using matrix-assisted laser desorption/ionization time-of-flight mass spectrometry in Belgium. Whole-genome sequencing (WGS) was used on the confirmed MRSA isolates, to characterize their genomic content and study their relatedness. Amongst the 305 isolates (304 wound swabs and 61 nasal swabs) that were collected from patients and HCWs, we detected 32 and 15 cases of MRSA, respectively. From this collection, 27 high-quality WGS datasets were obtained, which carried numerous genes and mutations associated with antimicrobial resistance. The mecA gene was detected in all the sequenced isolates. These isolates were assigned to five sequence types (STs), with ST8 (55.56%, n = 15/27), ST152 (18.52%, n = 5/27), and ST121 (18.52%, n = 5/27) being the most common. These 27 isolates carried multiple virulence genes, including the genes encoding the Panton-Valentine leukocidin toxin (48.15%, n = 13/27), and the tst gene (29.63%, n = 8/27), associated with toxic shock syndrome. This study highlights the need to implement a multimodal strategy for reducing the risk of the cross-transmission of MRSA in hospitals.

Transforming Shiga toxin-producing Escherichia coli surveillance through whole genome sequencing in food safety practices.

Introduction: Shiga toxin-producing Escherichia coli (STEC) is a gastrointestinal pathogen causing foodborne outbreaks. Whole Genome Sequencing (WGS) in STEC surveillance holds promise in outbreak prevention and confinement, in broadening STEC epidemiology and in contributing to risk assessment and source attribution. However, despite international recommendations, WGS is often restricted to assist outbreak investigation and is not yet fully implemented in food safety surveillance across all European countries, in contrast to for example in the United States. Methods: In this study, WGS was retrospectively applied to isolates collected within the context of Belgian food safety surveillance and combined with data from clinical isolates to evaluate its benefits. A cross-sector WGS-based collection of 754 strains from 1998 to 2020 was analyzed. Results: We confirmed that WGS in food safety surveillance allows accurate detection of genomic relationships between human cases and strains isolated from food samples, including those dispersed over time and geographical locations. Identifying these links can reveal new insights into outbreaks and direct epidemiological investigations to facilitate outbreak management. Complete WGS-based isolate characterization enabled expanding epidemiological insights related to circulating serotypes, virulence genes and antimicrobial resistance across different reservoirs. Moreover, associations between virulence genes and severe disease were determined by incorporating human metadata into the data analysis. Gaps in the surveillance system were identified and suggestions for optimization related to sample centralization, harmonizing isolation methods, and expanding sampling strategies were formulated. Discussion: This study contributes to developing a representative WGS-based collection of circulating STEC strains and by illustrating its benefits, it aims to incite policymakers to support WGS uptake in food safety surveillance.

Retrospective surveillance of viable Bacillus cereus group contaminations in commercial food and feed vitamin B2 products sold on the Belgian market using whole-genome sequencing.

Bacillus cereus is a spore-forming bacterium that occurs as a contaminant in food and feed, occasionally resulting in food poisoning through the production of various toxins. In this study, we retrospectively characterized viable B. cereus sensu lato (s.l.) isolates originating from commercial vitamin B2 feed and food additives collected between 2016 and 2022 by the Belgian Federal Agency for the Safety of the Food Chain from products sold on the Belgian market. In total, 75 collected product samples were cultured on a general medium and, in case of bacterial growth, two isolates per product sample were collected and characterized using whole-genome sequencing (WGS) and subsequently characterized in terms of sequence type (ST), virulence gene profile, antimicrobial resistance (AMR) gene profile, plasmid content, and phylogenomic relationships. Viable B. cereus was identified in 18 of the 75 (24%) tested products, resulting in 36 WGS datasets, which were classified into eleven different STs, with ST165 (n = 10) and ST32 (n = 8) being the most common. All isolates carried multiple genes encoding virulence factors, including cytotoxin K-2 (52.78%) and cereulide (22.22%). Most isolates were predicted to be resistant to beta-lactam antibiotics (100%) and fosfomycin (88.89%), and a subset was predicted to be resistant to streptothricin (30.56%). Phylogenomic analysis revealed that some isolates obtained from different products were closely related or even identical indicating a likely common origin, whereas for some products the two isolates obtained did not show any close relationship to each other or other isolates found in other products. This study reveals that potentially pathogenic and drug-resistant B. cereus s.l. can be present in food and feed vitamin B2 additives that are commercially available, and that more research is warranted to assess whether their presence in these types of products poses a threat to consumers.

Metagenomic Characterization of Multiple Genetically Modified Bacillus Contaminations in Commercial Microbial Fermentation Products.

Genetically modified microorganisms (GMM) are frequently employed for manufacturing microbial fermentation products such as food enzymes or vitamins. Although the fermentation product is required to be pure, GMM contaminations have repeatedly been reported in numerous commercial microbial fermentation produce types, leading to several rapid alerts at the European level. The aim of this study was to investigate the added value of shotgun metagenomic high-throughput sequencing to confirm and extend the results of classical analysis methods for the genomic characterization of unauthorized GMM. By combining short- and long-read metagenomic sequencing, two transgenic constructs were characterized, with insertions of alpha-amylase genes originating from B. amyloliquefaciens and B. licheniformis, respectively, and a transgenic construct with a protease gene insertion originating from B. velezensis, which were all present in all four investigated samples. Additionally, the samples were contaminated with up to three unculturable Bacillus strains, carrying genetic modifications that may hamper their ability to sporulate. Moreover, several samples contained viable Bacillus strains. Altogether these contaminations constitute a considerable load of antimicrobial resistance genes, that may represent a potential public health risk. In conclusion, our study showcases the added value of metagenomics to investigate the quality and safety of complex commercial microbial fermentation products.

Population Analysis of O26 Shiga Toxin-Producing Escherichia coli Causing Hemolytic Uremic Syndrome in Italy, 1989-2020, Through Whole Genome Sequencing.

Shiga toxin-producing Escherichia coli (STEC) belonging to the O26 serogroup represent an important cause of Hemolitic Uremic Syndrome (HUS) in children worldwide. The localization of STEC virulence genes on mobile genetic elements allowed the emergence of clones showing different assets of this accessory genomic fraction. A novel O26 STEC clone belonging to Sequence Type (ST) 29 and harboring stx2a, ehxA and etpD plasmid-borne genes has emerged and spread in Europe since the mid-1990s, while another ST29 clone positive for stx2d and lacking plasmid-borne virulence genes was recently described as emerging in France. In Italy, O26 has been the most frequently detected STEC serogroup from HUS cases since the late 1990s. In this study we describe the genomic characterization and population structure of 144 O26 STEC strains isolated from human sources in Italy in the period 1989-2020. A total of 89 strains belonged to ST21, 52 to ST29, two to ST396 and one to ST4944. ST29 strains started to be isolated from 1999. 24 strains were shown to harbour stx1a, alone (n=20) or in combination with stx2a (n=4). The majority of the strains (n=118) harbored stx2a genes only and the two ST396 strains harbored stx2d. A Hierarchical Clustering on Principal Components (HCPC) analysis, based on the detection of accessory virulence genes, antimicrobial resistance (AMR) genes and plasmid replicons, classified the strains in seven clusters identified with numbers from 1 to 7, containing two, 13, 39, 63, 16, 10 and one strain, respectively. The majority of the genetic features defining the clusters corresponded to plasmid-borne virulence genes, AMR genes and plasmid replicons, highlighting specific assets of plasmid-borne features associated with different clusters. Core genome Multi Locus Sequence Typing grouped ST21 and ST29 strains in three clades each, with each ST29 clade exactly corresponding to one HCPC cluster. Our results showed high conservation of either the core or the accessory genomic fraction in populations of ST29 O26 STEC, differently from what observed in ST21 strains, suggesting that a different selective pressure could drive the evolution of different populations of these pathogens possibly involving different ecological niches.

Retrospective evaluation of routine whole genome sequencing of Mycobacterium tuberculosis at the Belgian National Reference Center, 2019.

OBJECTIVES: Since January 2019, the Belgian National Reference Center for Mycobacteria (NRC) has switched from conventional typing to prospective whole-genome sequencing (WGS) of all submitted Mycobacterium tuberculosis complex (MTB) isolates. The ISO17025 validated procedure starts with semi-automated extraction and purification of gDNA directly from the submitted MGIT tubes, without preceding subculturing. All samples are then sequenced on an Illumina MiSeq sequencer and analyzed using an in-house developed and validated bioinformatics workflow to determine the species and antimicrobial resistance. In this study, we retrospectively compare results obtained via WGS to conventional phenotypic and genotypic testing, for all Belgian MTB strains analyzed in 2019 (n = 306). RESULTS: In all cases, the WGS-based procedure was able to identify correctly the MTB species. Compared to MGIT drug susceptibility testing (DST), the sensitivity and specificity of genetic prediction of resistance to first-line antibiotics were respectively 100 and 99% (rifampicin, RIF), 90.5 and 100% (isoniazid, INH), 100 and 98% (ethambutol, EMB) and 61.1 and 100% (pyrazinamide, PZA). The negative predictive value was above 95% for these four first-line drugs. A positive predictive value of 100% was calculated for INH and PZA, 80% for RIF and 45% for EMB. CONCLUSIONS: Our study confirms the effectiveness of WGS for the rapid detection of M. tuberculosis complex and its drug resistance profiles for first-line drugs even when working directly on MGIT tubes, and supports the introduction of this test into the routine workflow of laboratories performing tuberculosis diagnosis.

Whole-Genome Sequencing-Based Antimicrobial Resistance Characterization and Phylogenomic Investigation of 19 Multidrug-Resistant and Extended-Spectrum Beta-Lactamase-Positive Escherichia coli Strains Collected From Hospital Patients in Benin in 2019.

The increasing worldwide prevalence of extended-spectrum beta-lactamase (ESBL) producing Escherichia coli constitutes a serious threat to global public health. Surgical site infections are associated with high morbidity and mortality rates in developing countries, fueled by the limited availability of effective antibiotics. We used whole-genome sequencing (WGS) to evaluate antimicrobial resistance and the phylogenomic relationships of 19 ESBL-positive E. coli isolates collected from surgical site infections in patients across public hospitals in Benin in 2019. Isolates were identified by MALDI-TOF mass spectrometry and phenotypically tested for susceptibility to 16 antibiotics. Core-genome multi-locus sequence typing and single-nucleotide polymorphism-based phylogenomic methods were used to investigate the relatedness between samples. The broader phylogenetic context was characterized through the inclusion of publicly available genome data. Among the 19 isolates, 13 different sequence types (STs) were observed, including ST131 (n = 2), ST38 (n = 2), ST410 (n = 2), ST405 (n = 2), ST617 (n = 2), and ST1193 (n = 2). The bla CTX-M-15 gene encoding ESBL resistance was found in 15 isolates (78.9%), as well as other genes associated with ESBL, such as bla OXA-1 (n = 14) and bla TEM-1 (n = 9). Additionally, we frequently observed genes encoding resistance against aminoglycosides [aac-(6')-Ib-cr, n = 14], quinolones (qnrS1 , n = 4), tetracyclines [tet(B), n = 14], sulfonamides (sul2, n = 14), and trimethoprim (dfrA17, n = 13). Nonsynonymous chromosomal mutations in the housekeeping genes parC and gyrA associated with resistance to fluoroquinolones were also detected in multiple isolates. Although the phylogenomic investigation did not reveal evidence of hospital-acquired transmissions, we observed two very similar strains collected from patients in different hospitals. By characterizing a set of multidrug-resistant isolates collected from a largely unexplored environment, this study highlights the added value for WGS as an effective early warning system for emerging pathogens and antimicrobial resistance.

Characterization of Genetically Modified Microorganisms Using Short- and Long-Read Whole-Genome Sequencing Reveals Contaminations of Related Origin in Multiple Commercial Food Enzyme Products.

Despite their presence being unauthorized on the European market, contaminations with genetically modified (GM) microorganisms have repeatedly been reported in diverse commercial microbial fermentation produce types. Several of these contaminations are related to a GM Bacillus velezensis used to synthesize a food enzyme protease, for which genomic characterization remains currently incomplete, and it is unknown whether these contaminations have a common origin. In this study, GM B. velezensis isolates from multiple food enzyme products were characterized by short- and long-read whole-genome sequencing (WGS), demonstrating that they harbor a free recombinant pUB110-derived plasmid carrying antimicrobial resistance genes. Additionally, single-nucleotide polymorphism (SNP) and whole-genome based comparative analyses showed that the isolates likely originate from the same parental GM strain. This study highlights the added value of a hybrid WGS approach for accurate genomic characterization of GMM (e.g., genomic location of the transgenic construct), and of SNP-based phylogenomic analysis for source-tracking of GMM.

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.

Large diversity of linezolid-resistant isolates discovered in food-producing animals through linezolid selective monitoring in Belgium in 2019.

BACKGROUND: Linezolid is a critically important antibiotic used to treat human infections caused by MRSA and VRE. While linezolid is not licensed for food-producing animals, linezolid-resistant (LR) isolates have been reported in European countries, including Belgium. OBJECTIVES: To: (i) assess LR occurrence in staphylococci and enterococci isolated from different Belgian food-producing animals in 2019 through selective monitoring; and (ii) investigate the genomes and relatedness of these isolates. METHODS: Faecal samples (n = 1325) and nasal swab samples (n = 148) were analysed with a protocol designed to select LR bacteria, including a 44-48 h incubation period. The presence of LR chromosomal mutations, transferable LR genes and their genetic organizations and other resistance genes, as well as LR isolate relatedness (from this study and the NCBI database) were assessed through WGS. RESULTS: The LR rate differed widely between animal host species, with the highest rates occurring in nasal samples from pigs and sows (25.7% and 20.5%, respectively) and faecal samples from veal calves (16.4%). WGS results showed that LR determinants are present in a large diversity of isolates circulating in the agricultural sector, with some isolates closely related to human isolates, posing a human health risk. CONCLUSIONS: LR dedicated monitoring with WGS analysis could help to better understand the spread of LR. Cross-selection of LR transferable genes through other antibiotic use should be considered in future action plans aimed at combatting antimicrobial resistance and in future objectives for the rational use of antibiotics in a One Health perspective.

Whole-genome-based phylogenomic analysis of the Belgian 2016-2017 influenza A(H3N2) outbreak season allows improved surveillance.

Seasonal influenza epidemics are associated with high mortality and morbidity in the human population. Influenza surveillance is critical for providing information to national influenza programmes and for making vaccine composition predictions. Vaccination prevents viral infections, but rapid influenza evolution results in emerging mutants that differ antigenically from vaccine strains. Current influenza surveillance relies on Sanger sequencing of the haemagglutinin (HA) gene. Its classification according to World Health Organization (WHO) and European Centre for Disease Prevention and Control (ECDC) guidelines is based on combining certain genotypic amino acid mutations and phylogenetic analysis. Next-generation sequencing technologies enable a shift to whole-genome sequencing (WGS) for influenza surveillance, but this requires laboratory workflow adaptations and advanced bioinformatics workflows. In this study, 253 influenza A(H3N2) positive clinical specimens from the 2016-2017 Belgian season underwent WGS using the Illumina MiSeq system. HA-based classification according to WHO/ECDC guidelines did not allow classification of all samples. A new approach, considering the whole genome, was investigated based on using powerful phylogenomic tools including beast and Nextstrain, which substantially improved phylogenetic classification. Moreover, Bayesian inference via beast facilitated reassortment detection by both manual inspection and computational methods, detecting intra-subtype reassortants at an estimated rate of 15 %. Real-time analysis (i.e. as an outbreak is ongoing) via Nextstrain allowed positioning of the Belgian isolates into the globally circulating context. Finally, integration of patient data with phylogenetic groups and reassortment status allowed detection of several associations that would have been missed when solely considering HA, such as hospitalized patients being more likely to be infected with A(H3N2) reassortants, and the possibility to link several phylogenetic groups to disease severity indicators could be relevant for epidemiological monitoring. Our study demonstrates that WGS offers multiple advantages for influenza monitoring in (inter)national influenza surveillance, and proposes an improved methodology. This allows leveraging all information contained in influenza genomes, and allows for more accurate genetic characterization and reassortment detection.

A Bioinformatics Whole-Genome Sequencing Workflow for Clinical Mycobacterium tuberculosis Complex Isolate Analysis, Validated Using a Reference Collection Extensively Characterized with Conventional Methods and In Silico Approaches.

The use of whole-genome sequencing (WGS) for routine typing of bacterial isolates has increased substantially in recent years. For Mycobacterium tuberculosis (MTB), in particular, WGS has the benefit of drastically reducing the time required to generate results compared to most conventional phenotypic methods. Consequently, a multitude of solutions for analyzing WGS MTB data have been developed, but their successful integration in clinical and national reference laboratories is hindered by the requirement for their validation, for which a consensus framework is still largely absent. We developed a bioinformatics workflow for (Illumina) WGS-based routine typing of MTB complex (MTBC) member isolates allowing complete characterization, including (sub)species confirmation and identification (16S, csb/RD, hsp65), single nucleotide polymorphism (SNP)-based antimicrobial resistance (AMR) prediction, and pathogen typing (spoligotyping, SNP barcoding, and core genome multilocus sequence typing). Workflow performance was validated on a per-assay basis using a collection of 238 in-house-sequenced MTBC isolates, extensively characterized with conventional molecular biology-based approaches supplemented with public data. For SNP-based AMR prediction, results from molecular genotyping methods were supplemented with in silico modified data sets, allowing us to greatly increase the set of evaluated mutations. The workflow demonstrated very high performance with performance metrics of >99% for all assays, except for spoligotyping, where sensitivity dropped to ∼90%. The validation framework for our WGS-based bioinformatics workflow can aid in the standardization of bioinformatics tools by the MTB community and other SNP-based applications regardless of the targeted pathogen(s). The bioinformatics workflow is available for academic and nonprofit use through the Galaxy instance of our institute at https://galaxy.sciensano.be.

Phylogenomic Investigation of Increasing Fluoroquinolone Resistance among Belgian Cases of Shigellosis between 2013 and 2018 Indicates Both Travel-Related Imports and Domestic Circulation.

Shigellosis is an acute enteric infection caused mainly by the species Shigella flexneri and Shigella sonnei. Since surveillance of these pathogens indicated an increase in ciprofloxacin-resistant samples collected in Belgium between 2013 and 2018, a subset of 148 samples was analyzed with whole genome sequencing (WGS) to investigate their dispersion and underlying genomic features associated with ciprofloxacin resistance. A comparison between observed phenotypes and WGS-based resistance prediction to ciprofloxacin revealed perfect correspondence for all samples. Core genome multi-locus sequence typing and single nucleotide polymorphism-typing were used for phylogenomic investigation to characterize the spread of these infections within Belgium, supplemented with data from international reference collections to place the Belgian isolates within their global context. For S. flexneri, substantial diversity was observed with ciprofloxacin-resistant isolates assigned to several phylogenetic groups. Besides travel-related imports, several clusters of highly similar Belgian isolates could not be linked directly to international travel suggesting the presence of domestically circulating strains. For S. sonnei, Belgian isolates were all limited to lineage III, and could often be traced back to travel to countries in Asia and Africa, sometimes followed by domestic circulation. For both species, several clusters of isolates obtained exclusively from male patients were observed. Additionally, we illustrated the limitations of conventional serotyping of S. flexneri, which was impacted by serotype switching. This study contributes to a better understanding of the spread of shigellosis within Belgium and internationally, and highlights the added value of WGS for the surveillance of this pathogen.

Validation strategy of a bioinformatics whole genome sequencing workflow for Shiga toxin-producing Escherichia coli using a reference collection extensively characterized with conventional methods.

Whole genome sequencing (WGS) enables complete characterization of bacterial pathogenic isolates at single nucleotide resolution, making it the ultimate tool for routine surveillance and outbreak investigation. The lack of standardization, and the variation regarding bioinformatics workflows and parameters, however, complicates interoperability among (inter)national laboratories. We present a validation strategy applied to a bioinformatics workflow for Illumina data that performs complete characterization of Shiga toxin-producing Escherichia coli (STEC) isolates including antimicrobial resistance prediction, virulence gene detection, serotype prediction, plasmid replicon detection and sequence typing. The workflow supports three commonly used bioinformatics approaches for the detection of genes and alleles: alignment with blast+, kmer-based read mapping with KMA, and direct read mapping with SRST2. A collection of 131 STEC isolates collected from food and human sources, extensively characterized with conventional molecular methods, was used as a validation dataset. Using a validation strategy specifically adopted to WGS, we demonstrated high performance with repeatability, reproducibility, accuracy, precision, sensitivity and specificity above 95 % for the majority of all assays. The WGS workflow is publicly available as a 'push-button' pipeline at https://galaxy.sciensano.be. Our validation strategy and accompanying reference dataset consisting of both conventional and WGS data can be used for characterizing the performance of various bioinformatics workflows and assays, facilitating interoperability between laboratories with different WGS and bioinformatics set-ups.

First detection of a plasmid-encoded New-Delhi metallo-beta-lactamase-1 (NDM-1) producing Acinetobacter baumannii using whole genome sequencing, isolated in a clinical setting in Benin.

BACKGROUND: Carbapenem-resistant Acinetobacter baumannii is considered a top priority pathogen by the World Health Organization for combatting increasing antibiotic resistance and development of new drugs. Since it was originally reported in Klebsiella pneumoniae in 2009, the quick spread of the blaNDM-1 gene encoding a New-Delhi metallo-beta-lactamase-1 (NDM-1) is increasingly recognized as a serious threat. This gene is usually carried by large plasmids and has already been documented in diverse bacterial species, including A. baumannii. Here, we report the first detection of a NDM-1-producing A. baumannii strain isolated in Benin. CASE PRESENTATION: A 31-year-old woman was admitted to a surgical unit with a diagnosis of post-cesarean hematoma. An extensively-drug resistant A. baumannii strain solely susceptible to amikacin, colistin and ciprofloxacin, and resistant to several other antibiotics including ceftazidime, imipenem, meropenem, gentamicin, tobramycin, ceftazidime/avibactam, and sulfamethoxazole-trimethoprim, was isolated from the wound. Production of NDM-1 was demonstrated by immunochromatographic testing. Whole genome sequencing of the isolate confirmed the presence of blaNDM-1, but also antibiotic resistance genes against multiple beta-lactamases and other classes of antibiotics, in addition to several virulence genes. Moreover, the blaNDM-1 gene was found to be present in a Tn125 transposon integrated on a plasmid. CONCLUSIONS: The discovery of this extensively-drug resistant A. baumannii strain carrying blaNDM-1 in Benin is worrying, especially because of its high potential risk of horizontal gene transfer due to being integrated into a transposon located on a plasmid. Strict control and prevention measures should be taken, once NDM-1 positive A. baumannii has been identified to prevent transfer of this resistance gene to other Enterobacterales. Capacity building is required by governmental agencies to provide suitable antibiotic treatment options and strategies, in combination with strengthening laboratory services for detection and surveillance of this pathogen.

Impact of DNA extraction on whole genome sequencing analysis for characterization and relatedness of Shiga toxin-producing Escherichia coli isolates.

Whole genome sequencing (WGS) has proven to be the ultimate tool for bacterial isolate characterization and relatedness determination. However, standardized and harmonized workflows, e.g. for DNA extraction, are required to ensure robust and exchangeable WGS data. Data sharing between (inter)national laboratories is essential to support foodborne pathogen control, including outbreak investigation. This study evaluated eight commercial DNA preparation kits for their potential influence on: (i) DNA quality for Nextera XT library preparation; (ii) MiSeq sequencing (data quality, read mapping against plasmid and chromosome references); and (iii) WGS data analysis, i.e. isolate characterization (serotyping, virulence and antimicrobial resistance genotyping) and phylogenetic relatedness (core genome multilocus sequence typing and single nucleotide polymorphism analysis). Shiga toxin-producing Escherichia coli (STEC) was selected as a case study. Overall, data quality and inferred phylogenetic relationships between isolates were not affected by the DNA extraction kit choice, irrespective of the presence of confounding factors such as EDTA in DNA solution buffers. Nevertheless, completeness of STEC characterization was, although not substantially, influenced by the plasmid extraction performance of the kits, especially when using Nextera XT library preparation. This study contributes to addressing the WGS challenges of standardizing protocols to support data portability and to enable full exploitation of its potential.

The Benefits of Whole Genome Sequencing for Foodborne Outbreak Investigation from the Perspective of a National Reference Laboratory in a Smaller Country.

Gradually, conventional methods for foodborne pathogen typing are replaced by whole genome sequencing (WGS). Despite studies describing the overall benefits, National Reference Laboratories of smaller countries often show slower uptake of WGS, mainly because of significant investments required to generate and analyze data of a limited amount of samples. To facilitate this process and incite policy makers to support its implementation, a Shiga toxin-producing Escherichia coli (STEC) O157:H7 (stx1+, stx2+, eae+) outbreak (2012) and a STEC O157:H7 (stx2+, eae+) outbreak (2013) were retrospectively analyzed using WGS and compared with their conventional investigations. The corresponding results were obtained, with WGS delivering even more information, e.g., on virulence and antimicrobial resistance genotypes. Besides a universal, all-in-one workflow with less hands-on-time (five versus seven actual working days for WGS versus conventional), WGS-based cgMLST-typing demonstrated increased resolution. This enabled an accurate cluster definition, which remained unsolved for the 2013 outbreak, partly due to scarce epidemiological linking with the suspect source. Moreover, it allowed detecting two and one earlier circulating STEC O157:H7 (stx1+, stx2+, eae+) and STEC O157:H7 (stx2+, eae+) strains as closely related to the 2012 and 2013 outbreaks, respectively, which might have further directed epidemiological investigation initially. Although some bottlenecks concerning centralized data-sharing, sampling strategies, and perceived costs should be considered, we delivered a proof-of-concept that even in smaller countries, WGS offers benefits for outbreak investigation, if a sufficient budget is available to ensure its implementation in surveillance. Indeed, applying a database with background isolates is critical in interpreting isolate relationships to outbreaks, and leveraging the true benefit of WGS in outbreak investigation and/or prevention.

Identification of an unauthorized genetically modified bacteria in food enzyme through whole-genome sequencing.

Recently, the unexpected presence of a viable unauthorized genetically modified bacterium in a commercialized food enzyme (protease) product originating from a microbial fermentation process has been notified at the European level (RASFF 2019.3332). This finding was made possible thanks to the use of the next-generation sequencing technology, as reported in this study. Whole-genome sequencing was used to characterize the genetic modification comprising a sequence from the pUB110 shuttle vector (GenBank: M19465.1), harbouring antimicrobial resistance genes conferring a resistance to kanamycine, neomycin and bleomycin, flanked on each side by a sequence coding for a protease (GenBank: WP_032874795.1). In addition, based on these data, two real-time PCR methods, that can be used by enforcement laboratories, specific to this unauthorized genetically modified bacterium were developed and validated. The present study emphasizes the key role that whole-genome sequencing can take for detection of unknown and unauthorized genetically modified microorganisms in commercialized microbial fermentation products intended for the food and feed chain. Moreover, current issues encountered by the Competent Authorities and enforcement laboratories with such unexpected contaminations and the importance of performing official controls were highlighted.