Genomic comparison between Mycobacterium bovis and Mycobacterium microti and in silico analysis of peptide-based biomarkers for serodiagnosis.

In recent years, there has been an increase in the number of reported cases of Mycobacterium microti infection in various animals, which can interfere with the ante-mortem diagnosis of animal tuberculosis caused by Mycobacterium bovis. In this study, whole genome sequencing (WGS) was used to search for protein-coding genes to distinguish M. microti from M. bovis. In addition, the population structure of the available M. microti genomic WGS datasets is described, including three novel Belgian isolates from infections in alpacas. Candidate genes were identified by examining the presence of the regions of difference and by a pan-genome analysis of the available WGS data. A total of 80 genes showed presence-absence variation between the two species, including genes encoding Proline-Glutamate (PE), Proline-Proline-Glutamate (PPE), and Polymorphic GC-Rich Sequence (PE-PGRS) proteins involved in virulence and host interaction. Filtering based on predicted subcellular localization, sequence homology and predicted antigenicity resulted in 28 proteins out of 80 that were predicted to be potential antigens. As synthetic peptides are less costly and variable than recombinant proteins, an in silico approach was performed to identify linear and discontinuous B-cell epitopes in the selected proteins. From the 28 proteins, 157 B-cell epitope-based peptides were identified that discriminated between M. bovis and M. microti species. Although confirmation by in vitro testing is still required, these candidate synthetic peptides containing B-cell epitopes could potentially be used in serological tests to differentiate cases of M. bovis from M. microti infection, thus reducing misdiagnosis in animal tuberculosis surveillance.

Sheeppox virus genome sequences from the European outbreaks in Spain, Bulgaria, and Greece in 2022-2024.

In 2022-2024, three outbreaks of sheeppox (SPP) were reported in the European Union. These occurred in Spain, Bulgaria, and Greece and had serious economic consequences due to animal losses and trade restrictions. Five sheeppox virus (SPPV) whole-genome sequences (WGSs) were determined from samples collected during these outbreaks and analyzed in the context of all other published WGSs. Sheeppox virus strains can be divided in two, or possibly three, main groups. The isolates from the recent outbreaks belong to clade A2, which includes strains historically circulating in the Middle East and Northern Africa. Sequence divergence was low among the isolates that caused the recent European outbreaks. These results highlight the need for more regular and dense surveillance in under-sampled areas and the use of WGS to increase the chance of pinpointing the origin of an introduction, identifying potential introduction routes, and providing insights into SPPV evolution.

Benchmarking bacterial taxonomic classification using nanopore metagenomics data of several mock communities.

Taxonomic classification is crucial in identifying organisms within diverse microbial communities when using metagenomics shotgun sequencing. While second-generation Illumina sequencing still dominates, third-generation nanopore sequencing promises improved classification through longer reads. However, extensive benchmarking studies on nanopore data are lacking. We systematically evaluated performance of bacterial taxonomic classification for metagenomics nanopore sequencing data for several commonly used classifiers, using standardized reference sequence databases, on the largest collection of publicly available data for defined mock communities thus far (nine samples), representing different research domains and application scopes. Our results categorize classifiers into three categories: low precision/high recall; medium precision/medium recall, and high precision/medium recall. Most fall into the first group, although precision can be improved without excessively penalizing recall with suitable abundance filtering. No definitive 'best' classifier emerges, and classifier selection depends on application scope and practical requirements. Although few classifiers designed for long reads exist, they generally exhibit better performance. Our comprehensive benchmarking provides concrete recommendations, supported by publicly available code for reassessment and fine-tuning by other scientists.

Development and validation of a ddPCR assay to detect and quantify tobacco DNA in smoke and smokeless tobacco and tobacco-free products.

The last decade, smoke and smokeless products claiming to be tobacco-free, including herbal cigarettes and herbal shisha, became available on the European market and gained popularity. This study proposes a new digital droplet PCR (ddPCR) method, designed based on a previously developed real-time PCR (qPCR) method being currently used by the U.S. Food and Drug Administration (FDA) to specifically detect the presence of tobacco DNA in targeting a sequence from the Nicotiana tabacum nia-1 gene. To ensure a harmonized and reliable control by enforcement laboratories, both of these qPCR and ddPCR methods were then evaluated and validated for their compliance to an international standard. First, the performance of these PCR-based methods was successfully assessed as specific and sensitive, and in line with minimum performance requirements from international standard. Secondly, the transferability to external laboratory was confirmed for these PCR-based methods. Finally, the applicability of these PCR-based methods was demonstrated using 7 ground tobacco reference materials from the Tobacco Research Center (TRC) Toronto University as well as 6 commercial smokeless and tobacco-free smoke and smokeless products. Based on this study, the previously developed qPCR method was confirmed as complying with international standard, ensuring a efficient and harmonize use by enforcement laboratories for tobacco control on the European market. Moreover, this study proposed to enforcement laboratories the possibility to use a ddPCR method, enabling the simultaneous detection and absolute quantification of tobacco DNA as well as a limited impact of PCR inhibitors.

Towards facilitated interpretation of shotgun metagenomics long-read sequencing data analyzed with KMA for the detection of bacterial pathogens and their antimicrobial resistance genes.

Metagenomic sequencing is a promising method that has the potential to revolutionize the world of pathogen detection and antimicrobial resistance (AMR) surveillance in food-producing environments. However, the analysis of the huge amount of data obtained requires performant bioinformatics tools and databases, with intuitive and straightforward interpretation. In this study, based on long-read metagenomics data of chicken fecal samples with a spike-in mock community, we proposed confidence levels for taxonomic identification and AMR gene detection, with interpretation guidelines, to help with the analysis of the output data generated by KMA, a popular k-mer read alignment tool. Additionally, we demonstrated that the completeness and diversity of the genomes present in the reference databases are key parameters for accurate and easy interpretation of the sequencing data. Finally, we explored whether KMA, in a two-step procedure, can be used to link the detected AMR genes to their bacterial host chromosome, both detected within the same long-reads. The confidence levels were successfully tested on 28 metagenomics datasets which were obtained with sequencing of real and spiked samples from fecal (chicken, pig, and buffalo) or food (minced beef and food enzyme products) origin. The methodology proposed in this study will facilitate the analysis of metagenomics sequencing datasets for KMA users. Ultimately, this will contribute to improvements in the rapid diagnosis and surveillance of pathogens and AMR genes in food-producing environments, as prioritized by the EU.

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.

Strain-level characterization of foodborne pathogens without culture enrichment for outbreak investigation using shotgun metagenomics facilitated with nanopore adaptive sampling.

Introduction: Shotgun metagenomics has previously proven effective in the investigation of foodborne outbreaks by providing rapid and comprehensive insights into the microbial contaminant. However, culture enrichment of the sample has remained a prerequisite, despite the potential impact on pathogen detection resulting from the growth competition. To circumvent the need for culture enrichment, we explored the use of adaptive sampling using various databases for a targeted nanopore sequencing, compared to shotgun metagenomics alone. Methods: The adaptive sampling method was first tested on DNA of mashed potatoes mixed with DNA of a Staphylococcus aureus strain previously associated with a foodborne outbreak. The selective sequencing was used to either deplete the potato sequencing reads or enrich for the pathogen sequencing reads, and compared to a shotgun sequencing. Then, living S. aureus were spiked at 105 CFU into 25 g of mashed potatoes. Three DNA extraction kits were tested, in combination with enrichment using adaptive sampling, following whole genome amplification. After data analysis, the possibility to characterize the contaminant with the different sequencing and extraction methods, without culture enrichment, was assessed. Results: Overall, the adaptive sampling outperformed the shotgun sequencing. While the use of a host removal DNA extraction kit and targeted sequencing using a database of foodborne pathogens allowed rapid detection of the pathogen, the most complete characterization was achieved when using solely a database of S. aureus combined with a conventional DNA extraction kit, enabling accurate placement of the strain on a phylogenetic tree alongside outbreak cases. Discussion: This method shows great potential for strain-level analysis of foodborne outbreaks without the need for culture enrichment, thereby enabling faster investigations and facilitating precise pathogen characterization. The integration of adaptive sampling with metagenomics presents a valuable strategy for more efficient and targeted analysis of microbial communities in foodborne outbreaks, contributing to improved food safety and public health.

Pilot market surveillance of GMM contaminations in alpha-amylase food enzyme products: A detection strategy strengthened by a newly developed qPCR method targeting a GM Bacillus licheniformis producing alpha-amylase.

Using high-throughput metagenomics on commercial microbial fermentation products, DNA from a new unauthorized genetically modified microorganism (GMM), namely the GM B. licheniformis strain producing alpha-amylase (GMM alpha-amylase2), was recently discovered and characterized. On this basis, a new qPCR method targeting an unnatural association of sequences specific to the GMM alpha-amylase2 strain was designed and developed in this study, allowing to strengthen the current GMM detection strategy. The performance of the newly developed qPCR method was assessed for its specificity and sensitivity to comply with the minimum performance requirements established by the European Network of GMO Laboratories for GMO analysis. Moreover, the transferability of the in house validated qPCR method was demonstrated. Finally, its applicability was confirmed by a pilot market surveillance of GMM contaminations conducted for the first time on 40 alpha-amylase food enzyme products labelled as containing alpha-amylase. This pilot market surveillance allowed also to highlight numerous contaminations with GMM alpha-amylase2, including frequent cross-contaminations with other GMM strains previously characterized. In addition, the presence of full-length AMR genes, raising health concerns, was also reported.

Development of a portable on-site applicable metagenomic data generation workflow for enhanced pathogen and antimicrobial resistance surveillance.

Rapid, accurate and comprehensive diagnostics are essential for outbreak prevention and pathogen surveillance. Real-time, on-site metagenomics on miniaturized devices, such as Oxford Nanopore Technologies MinION sequencing, could provide a promising approach. However, current sample preparation protocols often require substantial equipment and dedicated laboratories, limiting their use. In this study, we developed a rapid on-site applicable DNA extraction and library preparation approach for nanopore sequencing, using portable devices. The optimized method consists of a portable mechanical lysis approach followed by magnetic bead-based DNA purification and automated sequencing library preparation, and resulted in a throughput comparable to a current optimal, laboratory-based protocol using enzymatic digestion to lyse cells. By using spike-in reference communities, we compared the on-site method with other workflows, and demonstrated reliable taxonomic profiling, despite method-specific biases. We also demonstrated the added value of long-read sequencing by recovering reads containing full-length antimicrobial resistance genes, and attributing them to a host species based on the additional genomic information they contain. Our method may provide a rapid, widely-applicable approach for microbial detection and surveillance in a variety of on-site settings.

Comparison of 6 DNA extraction methods for isolation of high yield of high molecular weight DNA suitable for shotgun metagenomics Nanopore sequencing to detect bacteria.

BACKGROUND: Oxford Nanopore Technologies (ONT) offers an accessible platform for long-read sequencing, which improves the reconstruction of genomes and helps to resolve complex genomic contexts, especially in the case of metagenome analysis. To take the best advantage of long-read sequencing, DNA extraction methods must be able to isolate pure high molecular weight (HMW) DNA from complex metagenomics samples, without introducing any bias. New methods released on the market, and protocols developed at the research level, were specifically designed for this application and need to be assessed. RESULTS: In this study, with different bacterial cocktail mixes, analyzed as pure or spiked in a synthetic fecal matrix, we evaluated the performances of 6 DNA extraction methods using various cells lysis and purification techniques, from quick and easy, to more time-consuming and gentle protocols, including a portable method for on-site application. In addition to the comparison of the quality, quantity and purity of the extracted DNA, the performance obtained when doing Nanopore sequencing on a MinION flow cell was also tested. From the obtained results, the Quick-DNA HMW MagBead Kit (Zymo Research) was selected as producing the best yield of pure HMW DNA. Furthermore, this kit allowed an accurate detection, by Nanopore sequencing, of almost all the bacterial species present in a complex mock community. CONCLUSION: Amongst the 6 tested methods, the Quick-DNA HMW MagBead Kit (Zymo Research) was considered as the most suitable for Nanopore sequencing and would be recommended for bacterial metagenomics studies using this technology.

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.

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.

Time series modelling for wastewater-based epidemiology of COVID-19: A nationwide study in 40 wastewater treatment plants of Belgium, February 2021 to June 2022.

BACKGROUND: Wastewater-based epidemiology (WBE) has been implemented to monitor surges of COVID-19. Yet, multiple factors impede the usefulness of WBE and quantitative adjustment may be required. AIM: We aimed to model the relationship between WBE data and incident COVID-19 cases, while adjusting for confounders and autocorrelation. METHODS: This nationwide WBE study includes data from 40 wastewater treatment plants (WWTPs) in Belgium (02/2021-06/2022). We applied ARIMA-based modelling to assess the effect of daily flow rate, pepper mild mottle virus (PMMoV) concentration, a measure of human faeces in wastewater, and variants (alpha, delta, and omicron strains) on SARS-CoV-2 RNA levels in wastewater. Secondly, adjusted WBE metrics at different lag times were used to predict incident COVID-19 cases. Model selection was based on AICc minimization. RESULTS: In 33/40 WWTPs, RNA levels were best explained by incident cases, flow rate, and PMMoV. Flow rate and PMMoV were associated with -13.0 % (95 % prediction interval: -26.1 to +0.2 %) and +13.0 % (95 % prediction interval: +5.1 to +21.0 %) change in RNA levels per SD increase, respectively. In 38/40 WWTPs, variants did not explain variability in RNA levels independent of cases. Furthermore, our study shows that RNA levels can lead incident cases by at least one week in 15/40 WWTPs. The median population size of leading WWTPs was 85.1 % larger than that of non­leading WWTPs. In 17/40 WWTPs, however, RNA levels did not lead or explain incident cases in addition to autocorrelation. CONCLUSION: This study provides quantitative insights into key determinants of WBE, including the effects of wastewater flow rate, PMMoV, and variants. Substantial inter-WWTP variability was observed in terms of explaining incident cases. These findings are of practical importance to WBE practitioners and show that the early-warning potential of WBE is WWTP-specific and needs validation.

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.

Development of a Droplet Digital PCR to Monitor SARS-CoV-2 Omicron Variant BA.2 in Wastewater Samples.

Wastewater-based surveillance can be used as a complementary method to other SARS-CoV-2 surveillance systems. It allows the emergence and spread of infections and SARS-CoV-2 variants to be monitored in time and place. This study presents an RT-ddPCR method that targets the T19I amino acid mutation in the spike protein of the SARS-CoV-2 genomes, which is specific to the BA.2 variant (omicron). The T19I assay was evaluated both in silico and in vitro for its inclusivity, sensitivity, and specificity. Moreover, wastewater samples were used as a proof of concept to monitor and quantify the emergence of the BA.2 variant from January until May 2022 in the Brussels-Capital Region which covers a population of more than 1.2 million inhabitants. The in silico analysis showed that more than 99% of the BA.2 genomes could be characterized using the T19I assay. Subsequently, the sensitivity and specificity of the T19I assay were successfully experimentally evaluated. Thanks to our specific method design, the positive signal from the mutant probe and wild-type probe of the T19I assay was measured and the proportion of genomes with the T19I mutation, characteristic of the BA.2 mutant, compared to the entire SARS-CoV-2 population was calculated. The applicability of the proposed RT-ddPCR method was evaluated to monitor and quantify the emergence of the BA.2 variant over time. To validate this assay as a proof of concept, the measurement of the proportion of a specific circulating variant with genomes containing the T19I mutation in comparison to the total viral population was carried out in wastewater samples from wastewater treatment plants in the Brussels-Capital Region in the winter and spring of 2022. This emergence and proportional increase in BA.2 genomes correspond to what was observed in the surveillance using respiratory samples; however, the emergence was observed slightly earlier, which suggests that wastewater sampling could be an early warning system and could be an interesting alternative to extensive human testing.

Development of Digital Droplet PCR Targeting the Influenza H3N2 Oseltamivir-Resistant E119V Mutation and Its Performance through the Use of Reverse Genetics Mutants.

The monitoring of antiviral-resistant influenza virus strains is important for public health given the availability and use of neuraminidase inhibitors and other antivirals to treat infected patients. Naturally occurring oseltamivir-resistant seasonal H3N2 influenza virus strains often carry a glutamate-to-valine substitution at position 119 in the neuraminidase (E119V-NA). Early detection of resistant influenza viruses is important for patient management and for the rapid containment of antiviral resistance. The neuraminidase inhibition assay allows the phenotypical identification of resistant strains; however, this test often has limited sensitivity with high variability depending on the virus strain, drugs and assays. Once a mutation such as E119V-NA is known, highly sensitive PCR-based genotypic assays can be used to identify the prevalence of such mutant influenza viruses in clinical samples. In this study, based on an existing reverse transcriptase real-time PCR (RT-qPCR) assay, we developed a reverse transcriptase droplet digital PCR assay (RT-ddPCR) to detect and quantify the frequency of the E119V-NA mutation. Furthermore, reverse genetics viruses carrying this mutation were created to test the performance of the RT-ddPCR assay and compare it to the standard phenotypic NA assay. We also discuss the advantage of using an RT-ddPCR instead of qPCR method in the context of viral diagnostics and surveillance.

Using a combination of short- and long-read sequencing to investigate the diversity in plasmid- and chromosomally encoded extended-spectrum beta-lactamases (ESBLs) in clinical Shigella and Salmonella isolates in Belgium.

For antimicrobial resistance (AMR) surveillance, it is important not only to detect AMR genes, but also to determine their plasmidic or chromosomal location, as this will impact their spread differently. Whole-genome sequencing (WGS) is increasingly used for AMR surveillance. However, determining the genetic context of AMR genes using only short-read sequencing is complicated. The combination with long-read sequencing offers a potential solution, as it allows hybrid assemblies. Nevertheless, its use in surveillance has so far been limited. This study aimed to demonstrate its added value for AMR surveillance based on a case study of extended-spectrum beta-lactamases (ESBLs). ESBL genes have been reported to occur also on plasmids. To gain insight into the diversity and genetic context of ESBL genes detected in clinical isolates received by the Belgian National Reference Center between 2013 and 2018, 100 ESBL-producing Shigella and 31 ESBL-producing Salmonella were sequenced with MiSeq and a representative selection of 20 Shigella and six Salmonella isolates additionally with MinION technology, allowing hybrid assembly. The bla CTX-M-15 gene was found to be responsible for a rapid rise in the ESBL Shigella phenotype from 2017. This gene was mostly detected on multi-resistance-carrying IncFII plasmids. Based on clustering, these plasmids were determined to be distinct from the circulating plasmids before 2017. They were spread to different Shigella species and within Shigella sonnei between multiple genotypes. Another similar IncFII plasmid was detected after 2017 containing bla CTX-M-27 for which only clonal expansion occurred. Matches of up to 99 % to plasmids of various bacterial hosts from all over the world were found, but global alignments indicated that direct or recent ESBL-plasmid transfers did not occur. It is most likely that travellers introduced these in Belgium and subsequently spread them domestically. However, a clear link to a specific country could not be made. Moreover, integration of bla CTX-M in the chromosome of two Shigella isolates was determined for the first time, and shown to be related to ISEcp1. In contrast, in Salmonella, ESBL genes were only found on plasmids, of which bla CTX-M-55 and IncHI2 were the most prevalent, respectively. No matching ESBL plasmids or cassettes were detected between clinical Shigella and Salmonella isolates. The hybrid assembly data allowed us to check the accuracy of plasmid prediction tools. MOB-suite showed the highest accuracy. However, these tools cannot replace the accuracy of long-read and hybrid assemblies. This study illustrates the added value of hybrid assemblies for AMR surveillance and shows that a strategy where even just representative isolates of a collection used for hybrid assemblies could improve international AMR surveillance as it allows plasmid tracking.

Targeted High-Throughput Sequencing Enables the Detection of Single Nucleotide Variations in CRISPR/Cas9 Gene-Edited Organisms.

Similar to genetically modified organisms (GMOs) produced by classical genetic engineering, gene-edited (GE) organisms and their derived food/feed products commercialized on the European Union market fall within the scope of European Union Directive 2001/18/EC. Consequently, their control in the food/feed chain by GMO enforcement laboratories is required by the competent authorities to guarantee food/feed safety and traceability (2003/1829/EC; 2003/1830/EC). However, their detection is potentially challenging at both the analytical and interpretation levels since this requires methodological approaches that can target and detect a specific single nucleotide variation (SNV) introduced into a GE organism. In this study, we propose a targeted high-throughput sequencing approach, including (i) a prior PCR-based enrichment step to amplify regions of interest, (ii) a sequencing step, and (iii) a data analysis methodology to identify SNVs of interest. To investigate if the performance of this targeted high-throughput sequencing approach is compatible with the performance criteria used in the GMO detection field, several samples containing different percentages of a GE rice line carrying a single adenosine insertion in OsMADS26 were prepared and analyzed. The SNV of interest in samples containing the GE rice line could successfully be detected, both at high and low percentages. No impact related to food processing or to the presence of other crop species was observed. The present proof-of-concept study has allowed us to deliver the first experimental-based evidence indicating that the proposed targeted high-throughput sequencing approach may constitute, in the future, a specific and sensitive tool to support the safety and traceability of the food/feed chain regarding GE plants carrying SNVs.

Noninvasive integrative approach applied to children in the context of recent air pollution exposure demonstrates association between fractional exhaled nitric oxide (FeNO) and urinary CC16.

Exposure to the air pollutant particulate matter (PM) is associated with increased risks of respiratory diseases and enhancement of airway inflammation in children. In the context of large scale air pollution studies, it can be challenging to measure fractional exhaled nitric oxide (FeNO) as indicator of lung inflammation. Urinary CC16 (U-CC16) is a potential biomarker of increased lung permeability and toxicity, increasing following short-term PM2.5 exposure. The single nucleotide polymorphism (SNP) CC16 G38A (rs3741240) affects CC16 levels and respiratory health. Our study aimed at assessing the use of U-CC16 (incl. CC16 G38A from saliva) as potential alternative for FeNO by investigating their mutual correlation in children exposed to PM. Samples from a small-scale study conducted in 42 children from urban (n = 19) and rural (n = 23) schools examined at two time points, were analysed. When considering recent (lag1) low level exposure to PM2.5 as air pollution measurement, we found that U-CC16 was positively associated with FeNO (ß = 0.23; 95% CI [-0.01; 0.47]; p = 0.06) in an adjusted analysis using a linear mixed effects model. Further, we observed a positive association between PM2.5 and FeNO (ß = 0.56; 95% CI [0.02; 1.09]; p = 0.04) and higher FeNO in urban school children as compared to rural school children (ß = 0.72; 95% CI [0.12; 1.31]; p = 0.02). Although more investigations are needed, our results suggest that inflammatory responses evidenced by increased FeNO are accompanied by potential increased lung epithelium permeability and injury, evidenced by increased U-CC16. In future large scale studies, where FeNO measurement is less feasible, the integrated analysis of U-CC16 and CC16 G38A, using noninvasive samples, might be a suitable alternative to assess the impact of air pollution exposure on the respiratory health of children, which is critical for policy development at population level.

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.