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.

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.

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.

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.

Listeria monocytogenes isolates from Cornu aspersum snails: Whole genome-based characterization and host-pathogen interactions in a snail infection model.

Even though Listeria monocytogenes is an extensive-studied foodborne pathogen, genome analysis of isolates from snails that may represent a reservoir of L. monocytogenes are still scarce. Here, we use whole-genome sequencing (WGS) to assess the genomic diversity of hypervirulent, virulent and non-virulent phenotypes of 15 L. monocytogenes isolated from snails to unveil their survival, virulence, and host-pathogen mechanisms of interactions in a snail infection model. Most of isolates (66.7%) were characterized as multidrug resistant (MDR) and belonged to clonal complexes (CCs) which are strongly associated with cases of human infection. All isolates contained intact genes associated with invasion and infection while hypervirulent isolates are adapted to host environment, possessing genes which are involved in teichoic acid biosynthesis, peptidoglycan modification and biofilm formation, correlating with their tolerance to haemolymph plasma phenotype and biofilm formation ability. A snail infection model showed that hypervirulent isolates triggered programmed host cell death pathway by increasing up to 30% the circulating apoptotic hemocytes in combination with induced nitrate production and reactive oxygen species (ROS) generation in snails' haemolymph. In contrast, the administration of the non-virulent strain which possesses a truncated mogR gene that regulates flagellar motility gene expression led only to an increase of necrotic non-apoptotic cells. Overall, this study provides significant insights into the genetic diversity of L. monocytogenes from snails, the genomic features of them linked to their hypervirulent/non-virulent phenotype, and the mechanisms of host-pathogen interactions.

Urinary CC16, a potential indicator of lung integrity and inflammation, increases in children after short-term exposure to PM2.5/PM10 and is driven by the CC16 38GG genotype.

Particular matter (PM) exposure is a big hazard for public health, especially for children. Serum CC16 is a well-known biomarker of respiratory health. Urinary CC16 (U-CC16) can be a noninvasive alternative, albeit requiring adequate adjustment for renal handling. Moreover, the SNP CC16 G38A influences CC16 levels. This study aimed to monitor the effect of short-term PM exposure on CC16 levels, measured noninvasively in schoolchildren, using an integrative approach. We used a selection of urine and buccal DNA samples from 86 children stored in an existing biobank. Using a multiple reaction monitoring method, we measured U-CC16, as well as RBP4 (retinol binding protein 4) and ß2M (beta-2-microglobulin), required for adjustment. Buccal DNA samples were used for CC16 G38A genotyping. Linear mixed-effects models were used to find relevant associations between U-CC16 and previously obtained data from recent daily PM ≤ 2.5 or 10 µm exposure (PM2.5, PM10) modeled at the child's residence. Our study showed that exposure to low PM at the child's residence (median levels 18.9 µg/m³ (PM2.5) and 23.6 µg/m³ (PM10)) one day before sampling had an effect on the covariates-adjusted U-CC16 levels. This effect was dependent on the CC16 G38A genotype, due to its strong interaction with the association between PM levels and covariates-adjusted U-CC16 (P = 0.024 (PM2.5); P = 0.061 (PM10)). Only children carrying the 38GG genotype showed an increase of covariates-adjusted U-CC16, measured 24h after exposure, with increasing PM2.5 and PM10 (ß = 0.332; 95% CI: 0.110 to 0.554 and ß = 0.372; 95% CI: 0.101 to 0.643, respectively). To the best of our knowledge, this is the first study using an integrative approach to investigate short-term PM exposure of children, using urine to detect early signs of pulmonary damage, and taking into account important determinants such as the genetic background and adequate adjustment of the measured biomarker in urine.

A systematic review and meta-analysis of host genetic factors associated with influenza severity.

BACKGROUND: The severity of influenza disease can range from mild symptoms to severe respiratory failure and can partly be explained by host genetic factors that predisposes the host to severe influenza. Here, we aimed to summarize the current state of evidence that host genetic variants play a role in the susceptibility to severe influenza infection by conducting a systematic review and performing a meta-analysis for all markers with at least three or more data entries. RESULTS: A total of 34 primary human genetic association studies were identified that investigated a total of 20 different genes. The only significant pooled ORs were retrieved for the rs12252 polymorphism: an overall OR of 1.52 (95% CI [1.06-2.17]) for the rs12252-C allele compared to the rs12252-T allele. A stratified analysis by ethnicity revealed opposite effects in different populations. CONCLUSION: With exception for the rs12252 polymorphism, we could not identify specific genetic polymorphisms to be associated with severe influenza infection in a pooled meta-analysis. This advocates for the use of large, hypothesis-free, genome-wide association studies that account for the polygenic nature and the interactions with other host, pathogen and environmental factors.

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.

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.

Evaluation of the added value of viral genomic information for predicting severity of influenza infection.

BACKGROUND: The severity of an influenza infection is influenced by both host and viral characteristics. This study aims to assess the relevance of viral genomic data for the prediction of severe influenza A(H3N2) infections among patients hospitalized for severe acute respiratory infection (SARI), in view of risk assessment and patient management. METHODS: 160 A(H3N2) influenza positive samples from the 2016-2017 season originating from the Belgian SARI surveillance were selected for whole genome sequencing. Predictor variables for severity were selected using a penalized elastic net logistic regression model from a combined host and genomic dataset, including patient information and nucleotide mutations identified in the viral genome. The goodness-of-fit of the model combining host and genomic data was compared using a likelihood-ratio test with the model including host data only. Internal validation of model discrimination was conducted by calculating the optimism-adjusted area under the Receiver Operating Characteristic curve (AUC) for both models. RESULTS: The model including viral mutations in addition to the host characteristics had an improved fit ([Formula: see text]=12.03, df = 3, p = 0.007). The optimism-adjusted AUC increased from 0.671 to 0.732. CONCLUSIONS: Adding genomic data (selected season-specific mutations in the viral genome) to the model containing host characteristics improved the prediction of severe influenza infection among hospitalized SARI patients, thereby offering the potential for translation into a prospective strategy to perform early season risk assessment or to guide individual patient management.

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.

Isolation of Drug-Resistant Gallibacterium anatis from Calves with Unresponsive Bronchopneumonia, Belgium.

Gallibacterium anatis is an opportunistic pathogen, previously associated with deaths in poultry, domestic birds, and occasionally humans. We obtained G. anatis isolates from bronchoalveolar lavage samples of 10 calves with bronchopneumonia unresponsive to antimicrobial therapy. Collected isolates were multidrug-resistant to extensively drug-resistant, exhibiting resistance against 5-7 classes of antimicrobial drugs. Whole-genome sequencing revealed 24 different antimicrobial-resistance determinants, including genes not previously described in the Gallibacterium genus or even the Pasteurellaceae family, such as aadA23, blaCARB-8, tet(Y), and qnrD1. Some resistance genes were closely linked in resistance gene cassettes with either transposases in close proximity or situated on putative mobile elements or predicted plasmids. Single-nucleotide polymorphism genotyping revealed large genetic variation between the G. anatis isolates, including isolates retrieved from the same farm. G. anatis might play a hitherto unrecognized role as a respiratory pathogen and resistance gene reservoir in cattle and has unknown zoonotic potential.

A multiplex oligonucleotide ligation-PCR method for the genoserotyping of common Salmonella using a liquid bead suspension assay.

Salmonella is a major pathogen having a public health and economic impact in both humans and animals. Six serotypes of the Salmonella genus are mentioned in the Belgian and European regulation as to be rapidly excluded from the food chain (EU regulation N°2160/2003, Belgian royal decree 27/04/2017). The reference method for Salmonella serotyping, including slide-agglutination and biochemical tests, is time-consuming, expensive, not always objective, and therefore does not match the fast identification criteria required by the legislation. In this study, a molecular method, using genetic markers detected by Multiplex Oligonucleotide Ligation - PCR and Luminex technology, was developed for the identification of the 6 Salmonella serotypes and their variants subjected to an official control. The resulting method was validated with the analysis of 971 Salmonella isolated from different matrixes (human, animal, food or environment) and 33 non-Salmonella strains. The results were compared with the reference identifications, achieving an accuracy of 99.7%. The cost-effective high-throughput genoserotyping assay is performed in 1 day and generates objective results, thanks to the automatic interpretation of raw data using a barcode system. In conclusion, it is fully adapted to the implementation in first line laboratories and meets the requirements of the regulation.

Strain-Level Metagenomic Data Analysis of Enriched In Vitro and In Silico Spiked Food Samples: Paving the Way towards a Culture-Free Foodborne Outbreak Investigation Using STEC as a Case Study.

Culture-independent diagnostics, such as metagenomic shotgun sequencing of food samples, could not only reduce the turnaround time of samples in an outbreak investigation, but also allow the detection of multi-species and multi-strain outbreaks. For successful foodborne outbreak investigation using a metagenomic approach, it is, however, necessary to bioinformatically separate the genomes of individual strains, including strains belonging to the same species, present in a microbial community, which has up until now not been demonstrated for this application. The current work shows the feasibility of strain-level metagenomics of enriched food matrix samples making use of data analysis tools that classify reads against a sequence database. It includes a brief comparison of two database-based read classification tools, Sigma and Sparse, using a mock community obtained by in vitro spiking minced meat with a Shiga toxin-producing Escherichia coli (STEC) isolate originating from a described outbreak. The more optimal tool Sigma was further evaluated using in silico simulated metagenomic data to explore the possibilities and limitations of this data analysis approach. The performed analysis allowed us to link the pathogenic strains from food samples to human isolates previously collected during the same outbreak, demonstrating that the metagenomic approach could be applied for the rapid source tracking of foodborne outbreaks. To our knowledge, this is the first study demonstrating a data analysis approach for detailed characterization and phylogenetic placement of multiple bacterial strains of one species from shotgun metagenomic WGS data of an enriched food sample.

Use of Whole Genome Sequencing Data for a First in Silico Specificity Evaluation of the RT-qPCR Assays Used for SARS-CoV-2 Detection.

The current COronaVIrus Disease 2019 (COVID-19) pandemic started in December 2019. COVID-19 cases are confirmed by the detection of SARS-CoV-2 RNA in biological samples by RT-qPCR. However, limited numbers of SARS-CoV-2 genomes were available when the first RT-qPCR methods were developed in January 2020 for initial in silico specificity evaluation and to verify whether the targeted loci are highly conserved. Now that more whole genome data have become available, we used the bioinformatics tool SCREENED and a total of 4755 publicly available SARS-CoV-2 genomes, downloaded at two different time points, to evaluate the specificity of 12 RT-qPCR tests (consisting of a total of 30 primers and probe sets) used for SARS-CoV-2 detection and the impact of the virus' genetic evolution on four of them. The exclusivity of these methods was also assessed using the human reference genome and 2624 closely related other respiratory viral genomes. The specificity of the assays was generally good and stable over time. An exception is the first method developed by the China Center for Disease Control and prevention (CDC), which exhibits three primer mismatches present in 358 SARS-CoV-2 genomes sequenced mainly in Europe from February 2020 onwards. The best results were obtained for the assay of Chan et al. (2020) targeting the gene coding for the spiking protein (S). This demonstrates that our user-friendly strategy can be used for a first in silico specificity evaluation of future RT-qPCR tests, as well as verifying that the former methods are still capable of detecting circulating SARS-CoV-2 variants.

The genetic structure of the Belgian population.

BACKGROUND: National and international efforts like the 1000 Genomes Project are leading to increasing insights in the genetic structure of populations worldwide. Variation between different populations necessitates access to population-based genetic reference datasets. These data, which are important not only in clinical settings but also to potentiate future transitions towards a more personalized public health approach, are currently not available for the Belgian population. RESULTS: To obtain a representative genetic dataset of the Belgian population, participants in the 2013 National Health Interview Survey (NHIS) were invited to donate saliva samples for DNA analysis. DNA was isolated and single nucleotide polymorphisms (SNPs) were determined using a genome-wide SNP array of around 300,000 sites, resulting in a high-quality dataset of 189 samples that was used for further analysis. A principal component analysis demonstrated the typical European genetic constitution of the Belgian population, as compared to other continents. Within Europe, the Belgian population could be clearly distinguished from other European populations. Furthermore, obvious signs from recent migration were found, mainly from Southern Europe and Africa, corresponding with migration trends from the past decades. Within Belgium, a small north-west to south-east gradient in genetic variability was noted, with differences between Flanders and Wallonia. CONCLUSIONS: This is the first study on the genetic structure of the Belgian population and its regional variation. The Belgian genetic structure mirrors its geographic location in Europe with regional differences and clear signs of recent migration.

Application of whole genome data for in silico evaluation of primers and probes routinely employed for the detection of viral species by RT-qPCR using dengue virus as a case study.

BACKGROUND: Viral infection by dengue virus is a major public health problem in tropical countries. Early diagnosis and detection are increasingly based on quantitative reverse transcriptase real-time polymerase chain reaction (RT-qPCR) directed against genomic regions conserved between different isolates. Genetic variation can however result in mismatches of primers and probes with their targeted nucleic acid regions. Whole genome sequencing allows to characterize and track such changes, which in turn enables to evaluate, optimize, and (re-)design novel and existing RT-qPCR methods. The immense amount of available sequence data renders this however a labour-intensive and complex task. RESULTS: We present a bioinformatics approach that enables in silico evaluation of primers and probes intended for routinely employed RT-qPCR methods. This approach is based on analysing large amounts of publically available whole genome data, by first employing BLASTN to mine the genomic regions targeted by the RT-qPCR method(s), and afterwards using BLASTN-SHORT to evaluate whether primers and probes will anneal based on a set of simple in silico criteria. Using dengue virus as a case study, we evaluated 18 published RT-qPCR methods using more than 3000 publically available genomes in the NCBI Virus Variation Resource, and provide a systematic overview of method performance based on in silico sensitivity and specificity. CONCLUSIONS: We provide a comprehensive overview of dengue virus RT-qPCR method performance that will aid appropriate method selection allowing to take specific measures that aim to contain and prevent viral spread in afflicted regions. Notably, we find that primer-template mismatches at their 3' end may represent a general issue for dengue virus RT-qPCR detection methods that merits more attention in their development process. Our approach is also available as a public tool, and demonstrates how utilizing genomic data can provide meaningful insights in an applied public health setting such as the detection of viral species in human diagnostics.

Isolation of Burkholderia pseudomallei from a Pet Green Iguana, Belgium.

We isolated Burkholderia pseudomallei, the causative agent of melioidosis, from liver granulomas of a pet green iguana (Iguana iguana) in Belgium. This case highlights a risk for imported green iguanas acting as a reservoir for introduction of this high-threat, zoonotic pathogen into nonendemic regions.

Development and validation of an integrated DNA walking strategy to detect GMO expressing cry genes.

BACKGROUND: Recently, an integrated DNA walking strategy has been proposed to prove the presence of GMO via the characterisation of sequences of interest, including their transgene flanking regions and the unnatural associations of elements in their transgenic cassettes. To this end, the p35S, tNOS and t35S pCAMBIA elements have been selected as key targets, allowing the coverage of most of GMO, EU authorized or not. In the present study, a bidirectional DNA walking method anchored on the CryAb/c genes is proposed with the aim to cover additional GMO and additional sequences of interest. RESULTS: The performance of the proposed bidirectional DNA walking method anchored on the CryAb/c genes has been evaluated in a first time for its feasibility using several GM events possessing these CryAb/c genes. Afterwards, its sensitivity has been investigated through low concentrations of targets (as low as 20 HGE). In addition, to illustrate its applicability, the entire workflow has been tested on a sample mimicking food/feed matrices analysed in GMO routine analysis. CONCLUSION: Given the successful assessment of its performance, the present bidirectional DNA walking method anchored on the CryAb/c genes can easily be implemented in GMO routine analysis by the enforcement laboratories and allows completing the entire DNA walking strategy in targeting an additional transgenic element frequently found in GMO.

Detection and discrimination of five E. coli pathotypes using a combinatory SYBR® Green qPCR screening system.

A detection and discrimination system for five Escherichia coli pathotypes, based on a combination of 13 SYBR® Green qPCR, has been developed, i.e., combinatory SYBR® Green qPCR screening system for pathogenic E. coli (CoSYPS Path E. coli). It allows the discrimination on isolates and the screening of potential presence in food of the following pathotypes of E. coli: shigatoxigenic (STEC) (including enterohemorrhagic (EHEC)), enteropathogenic (EPEC), enteroaggregative (EAggEC), enteroaggregative shigatoxigenic (EAggSTEC), and enteroinvasive (EIEC) E. coli. The SYBR® Green qPCR assays target the uidA, ipaH, eae, aggR, aaiC, stx1, and stx2 genes. uidA controls for E. coli presence and all the other genes are specific targets of E. coli pathotypes. For each gene, two primer pairs have been designed to guarantee a sufficient detection even in case of deletion or polymorphisms in the target gene. Moreover, all the qPCR have been designed to be run together in a single analytical PCR plate. This study includes the primer pairs' design, in silico and in situ selectivity, sensitivity, repeatability, and reproducibility evaluation of the 13 SYBR® Green qPCR assays. Each target displayed a selectivity of 100%. The limit of detection of the 13 assays is between 1 and 10 genomic copies. Their repeatability and reproducibility comply with the European requirements. As a preliminary feasibility study on food, the CoSYPS Path E. coli system was subsequently evaluated on four food matrices artificially contaminated with pathogenic E. coli. It allowed the detection of an initial contamination level as low as 2 to 7 cfu of STEC/25 g of food matrix after 24 h of enrichment.