First detection of a plasmid located carbapenem resistant bla(VIM-1) gene in E. coli isolated from meat products at retail in Belgium in 2015

Carbapenemase-producing Enterobacteriaceae (CPE) confer resistance to antibiotics that are of critical importance to human medicine. There have only been a few reported cases of CPEs in the European food chain. We report the first detection of a carbapenemase-producing Escherichia coli (ST 5869) in the Belgian food chain. Our aim was to characterize the origin of the carbapenem resistance in the E. coli isolate. The isolate was detected during the screening of 178 minced pork samples and was shown to contain the carbapenemase gene bla(VIM-1) by PCR and Sanger sequencing. Whole genome short and long read sequencing (MiSeq and MinION) was performed to characterize the isolate. With a hybrid assembly we reconstructed a 190,205 bp IncA/C2 plasmid containing bla(VIM-1) (S15FP06257_p), in addition to other critically important resistance genes. This plasmid showed only low similarity to plasmids containing bla(VIM-1) previously reported in Germany. Moreover, no sequences existed in the NCBI nucleotide database that completely covered S15FP06257_p. Analysis of the bla(VIM-1) gene cassette demonstrated that it likely originated from an integron of a Klebsiella plasmid reported previously in a clinical isolate in Europe, suggesting that the meat could have been contaminated by human handling in one of the steps of the food chain. This study shows the relevance of fully reconstructing plasmids to characterize their genetic content and to allow source attribution. This is especially important in view of the potential risk of antimicrobial resistance gene transmission through mobile elements as was reported here for the of public health concern bla(VIM-1)

A practical method to implement strain-level metagenomics-based foodborne outbreak investigation and source tracking in routine

The management of a foodborne outbreak depends on the rapid and accurate identification of the responsible food source. Conventional methods based on isolation of the pathogen from the food matrix and target-specific real-time polymerase chain reactions (qPCRs) are used in routine. In recent years, the use of whole genome sequencing (WGS) of bacterial isolates has proven its value to collect relevant information for strain characterization as well as tracing the origin of the contamination by linking the food isolate with the patient’s isolate with high resolution. However, the isolation of a bacterial pathogen from food matrices is often time-consuming and not always successful. Therefore, we aimed to improve outbreak investigation by developing a method that can be implemented in reference laboratories to characterize the pathogen in the food vehicle without its prior isolation and link it back to human cases. We tested and validated a shotgun metagenomics approach by spiking food pathogens in specific food matrices using the Shiga toxin-producing Escherichia coli (STEC) as a case study. Different DNA extraction kits and enrichment procedures were investigated to obtain the most practical workflow. We demonstrated the feasibility of shotgun metagenomics to obtain the same information as in ISO/TS 13136:2012 and WGS of the isolate in parallel by inferring the genome of the contaminant and characterizing it in a shorter timeframe. This was achieved in food samples containing different E. coli strains, including a combination of different STEC strains. For the first time, we also managed to link individual strains from a food product to isolates from human cases, demonstrating the power of shotgun metagenomics for rapid outbreak investigation and source tracking

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

Application of a strain- level shotgun metagenomics approach on food samples : resolution of the source of a Salmonella food-borne outbreak

Food- borne outbreak investigation currently relies on the time- consuming and challenging bacterial isolation from food, to be able to link food- derived strains to more easily obtained isolates from infected people. When no food isolate can be obtained, the source of the outbreak cannot be unambiguously determined. Shotgun metagenomics approaches applied to the food samples could circumvent this need for isolation from the suspected source, but require downstream strain- level data analysis to be able to accurately link to the human isolate. Until now, this approach has not yet been applied outside research settings to analyse real food- borne outbreak samples. In September 2019, a Salmonella outbreak occurred in a hotel school in Bruges, Belgium, affecting over 200 students and teachers. Following standard procedures, the Belgian National Reference Center for human salmonellosis and the National Reference Laboratory for Salmonella in food and feed used conventional analysis based on isolation, serotyping and MLVA (multilocus variable number tandem repeat analysis) comparison, followed by wholegenome sequencing, to confirm the source of the contamination over 2 weeks after receipt of the sample, which was freshly prepared tartar sauce in a meal cooked at the school. Our team used this outbreak as a case study to deliver a proof of concept for a short- read strain- level shotgun metagenomics approach for source tracking. We received two suspect food samples: the full meal and some freshly made tartar sauce served with this meal, requiring the use of raw eggs. After analysis, we could prove, without isolation, that Salmonella was present in both samples, and we obtained an inferred genome of a Salmonella enterica subsp. enterica serovar Enteritidis that could be linked back to the human isolates of the outbreak in a phylogenetic tree. These metagenomics- derived outbreak strains were separated from sporadic cases as well as from another outbreak circulating in Europe at the same time period. This is, to our knowledge, the first Salmonella food- borne outbreak investigation uniquely linking the food source using a metagenomics approach and this in a fast time frame

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

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

IntroductionShiga 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.MethodsIn 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.ResultsWe 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.DiscussionThis 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

Osteopontin expression identifies a subset of recruited macrophages distinct from Kupffer cells in the fatty liver

Metabolic-associated fatty liver disease (MAFLD) represents a spectrum of disease states ranging from simple steatosis to non-alcoholic steatohepatitis (NASH). Hepatic macrophages, specifically Kupffer cells (KCs), are suggested to play important roles in the pathogenesis of MAFLD through their activation, although the exact roles played by these cells remain unclear. Here, we demonstrated that KCs were reduced in MAFLD being replaced by macrophages originating from the bone marrow. Recruited macrophages existed in two subsets with distinct activation states, either closely resembling homeostatic KCs or lipid-associated macrophages (LAMs) from obese adipose tissue. Hepatic LAMs expressed Osteopontin, a biomarker for patients with NASH, linked with the development of fibrosis. Fitting with this, LAMs were found in regions of the liver with reduced numbers of KCs, characterized by increased Desmin expression. Together, our data highlight considerable heterogeneity within the macrophage pool and suggest a need for more specific macrophage targeting strategies in MAFLD

Stellate cells, hepatocytes, and endothelial cells imprint the Kupffer cell identity on monocytes colonizing the liver macrophage niche

Macrophages are strongly adapted to their tissue of residence. Yet, little is known about the cell-cell interactions that imprint the tissue-specific identities of macrophages in their respective niches. Using conditional depletion of liver Kupffer cells, we traced the developmental stages of monocytes differentiating into Kupffer cells and mapped the cellular interactions imprinting the Kupffer cell identity. Kupffer cell loss induced tumor necrosis factor (TNF)- and interleukin-1 (IL-1) receptor-dependent activation of stellate cells and endothelial cells, resulting in the transient production of chemokines and adhesion molecules orchestrating monocyte engraftment. Engrafted circulating monocytes transmigrated into the perisinusoidal space and acquired the liver-associated transcription factors inhibitor of DNA 3 (ID3) and liver X receptor-alpha (LXR-alpha). Coordinated interactions with hepatocytes induced ID3 expression, whereas endothelial cells and stellate cells induced LXR-alpha via a synergistic NOTCH-BMP pathway. This study shows that the Kupffer cell niche is composed of stellate cells, hepatocytes, and endothelial cells that together imprint the liver-specific macrophage identity

Complementarity of single-nucleus and single-cell RNA sequencing to study the transcriptional heterogeneity of neuroblastoma tumors

C

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 data generated and analyzed during the current study are available in the NCBI SRA repository (https:// www.ncbi.nlm.nih.gov/sra) under accession number PRJNA574887 (in-house sequenced data) and its accession numbers are listed in Supplementary Table S11 online.Sciensano, Belgium and the Belgian Federal Public Service of Health, Food Chain Safety and Environment.http://www.nature.com/srepam2021Genetic