Datasets for benchmarking antimicrobial resistance genes in bacterial metagenomic and whole genome sequencing.

Whole genome sequencing (WGS) is a key tool in identifying and characterising disease-associated bacteria across clinical, agricultural, and environmental contexts. One increasingly common use of genomic and metagenomic sequencing is in identifying the type and range of antimicrobial resistance (AMR) genes present in bacterial isolates in order to make predictions regarding their AMR phenotype. However, there are a large number of alternative bioinformatics software and pipelines available, which can lead to dissimilar results. It is, therefore, vital that researchers carefully evaluate their genomic and metagenomic AMR analysis methods using a common dataset. To this end, as part of the Microbial Bioinformatics Hackathon and Workshop 2021, a 'gold standard' reference genomic and simulated metagenomic dataset was generated containing raw sequence reads mapped against their corresponding reference genome from a range of 174 potentially pathogenic bacteria. These datasets and their accompanying metadata are freely available for use in benchmarking studies of bacteria and their antimicrobial resistance genes and will help improve tool development for the identification of AMR genes in complex samples.

Harmonization of whole-genome sequencing for outbreak surveillance of Enterobacteriaceae and Enterococci.

Whole-genome sequencing (WGS) is becoming the de facto standard for bacterial typing and outbreak surveillance of resistant bacterial pathogens. However, interoperability for WGS of bacterial outbreaks is poorly understood. We hypothesized that harmonization of WGS for outbreak surveillance is achievable through the use of identical protocols for both data generation and data analysis. A set of 30 bacterial isolates, comprising of various species belonging to the Enterobacteriaceae family and Enterococcus genera, were selected and sequenced using the same protocol on the Illumina MiSeq platform in each individual centre. All generated sequencing data were analysed by one centre using BioNumerics (6.7.3) for (i) genotyping origin of replications and antimicrobial resistance genes, (ii) core-genome multi-locus sequence typing (cgMLST) for Escherichia coli and Klebsiella pneumoniae and whole-genome multi-locus sequencing typing (wgMLST) for all species. Additionally, a split k-mer analysis was performed to determine the number of SNPs between samples. A precision of 99.0% and an accuracy of 99.2% was achieved for genotyping. Based on cgMLST, a discrepant allele was called only in 2/27 and 3/15 comparisons between two genomes, for E. coli and K. pneumoniae, respectively. Based on wgMLST, the number of discrepant alleles ranged from 0 to 7 (average 1.6). For SNPs, this ranged from 0 to 11 SNPs (average 3.4). Furthermore, we demonstrate that using different de novo assemblers to analyse the same dataset introduces up to 150 SNPs, which surpasses most thresholds for bacterial outbreaks. This shows the importance of harmonization of data-processing surveillance of bacterial outbreaks. In summary, multi-centre WGS for bacterial surveillance is achievable, but only if protocols are harmonized.

Higher Prevalence of Bacteroides fragilis in Crohn's Disease Exacerbations and Strain-Dependent Increase of Epithelial Resistance.

Bacteroides fragilis has previously been linked to Crohn's disease (CD) exacerbations, but results are inconsistent and underlying mechanisms unknown. This study investigates the epidemiology of B. fragilis and its virulence factors bft (enterotoxin) and ubiquitin among 181 CD patients and the impact on the intestinal epithelial barrier in vitro. The prevalence of B. fragilis was significantly higher in active (n = 69/88, 78.4%) as compared to remissive (n = 58/93, 62.4%, p = 0.018) CD patients. Moreover, B. fragilis was associated with intestinal strictures. Interestingly, the intestinal barrier function, as examined by transepithelial electrical resistance (TEER) measurements of Caco-2 monolayers, increased when exposed to secretomes of bft-positive (bft-1 and bft-2 isotype; increased TEER ∼160%, p < 0.001) but not when exposed to bft-negative strains. Whole metagenome sequencing and metabolomics, respectively, identified nine coding sequences and two metabolites that discriminated TEER-increasing from non-TEER-increasing strains. This study revealed a higher B. fragilis prevalence during exacerbation. Surprisingly, bft-positive secretomes increased epithelial resistance, but we excluded Bft as the likely causative factor.

Genetic analysis of plasmid-encoded mcr-1 resistance in Enterobacteriaceae derived from poultry meat in the Netherlands.

BACKGROUND: Colistin is classified as the highest priority and critically important antimicrobial for human medicine by WHO as it is the last resort agent for treatment of carbapenem-resistant Enterobacteriaceae in humans. Additional research is necessary to elucidate the genetic structure of mcr-1 resistance genes, commonly found on plasmids, using WGS. OBJECTIVES: To map and compare the genetic characteristics of 35 mcr-1-mediated colistin-resistant Enterobacteriaceae isolated from chicken meat to highlight the genetic variation of the mcr-1-containing plasmids. METHODS: Sequencing was performed using Illumina HiSeq2500, Novaseq6000 and ONT's GridION. GridION data was locally basecalled and demultiplexed using ONT's Albacore 2.3.4 followed by Porechop 2.3. Quality filtering was performed using Filtlong 2.0. Hybrid Assembly was performed using Unicycler 4.7. Plasmids were compared with reference sequences in plasmid-RefSeq and pATLAS. RESULTS: A total of 35 mcr-1 positive Enterobacteriaceae were investigated, which resulted in 34 qualitatively robust hybrid assemblies of 2 Klebsiella pneumoniae and 32 Escherichia coli. mcr-1.1 was present in 33/34 isolates. One isolate contained an mcr-1.1-like resistance gene, due to a deletion of one codon. Two mcr-1.1 genes were located on the chromosome, while the majority of the mcr-1 genes were found on IncX4 type plasmids (n = 19). Almost all plasmids identified in this study were highly similar to plasmids found in human-derived strains. CONCLUSIONS: The mcr-1.1-containing plasmids from retail chicken show high sequence similarity to human mcr-1.1 plasmids, suggesting that this may be a contributor to the presence of colistin resistance in humans.

Centre-specific bacterial pathogen typing affects infection-control decision making.

Whole-genome sequencing is becoming the de facto standard for bacterial outbreak surveillance and infection prevention. This is accompanied by a variety of bioinformatic tools and needs bioinformatics expertise for implementation. However, little is known about the concordance of reported outbreaks when using different bioinformatic workflows. In this multi-centre proficiency testing among 13 major Dutch healthcare-affiliated centres, bacterial whole-genome outbreak analysis was assessed. Centres who participated obtained two randomized bacterial datasets of Illumina sequences, a Klebsiella pneumoniae and a Vancomycin-resistant Enterococcus faecium, and were asked to apply their bioinformatic workflows. Centres reported back on antimicrobial resistance, multi-locus sequence typing (MLST), and outbreak clusters. The reported clusters were analysed using a method to compare landscapes of phylogenetic trees and calculating Kendall-Colijn distances. Furthermore, fasta files were analysed by state-of-the-art single nucleotide polymorphism (SNP) analysis to mitigate the differences introduced by each centre and determine standardized SNP cut-offs. Thirteen centres participated in this study. The reported outbreak clusters revealed discrepancies between centres, even when almost identical bioinformatic workflows were used. Due to stringent filtering, some centres failed to detect extended-spectrum beta-lactamase genes and MLST loci. Applying a standardized method to determine outbreak clusters on the reported de novo assemblies, did not result in uniformity of outbreak-cluster composition among centres.

Appearance of vanD-positive Enterococcus faecium in a tertiary hospital in the Netherlands: prevalence of vanC and vanD in hospitalized patients.

Vancomycin-resistant enterococci (VRE) can rapidly spread through hospitals. Therefore, our hospital employs a screening program whereby rectal swabs are screened for the presence of vanA and vanB, and only PCR-positive broths are cultured on VRE selection agar. Early November 2016, a clinical vanA-/vanB-negative VRE isolate was detected in a vanA/vanB-screening-negative patient, giving the possibility that an undetected VRE might be spreading within our hospital. Whole-genome-sequencing of the isolate showed that resistance was vanD-mediated and core genome multilocus sequence typing showed it was a rare type: ST17/CT154. To determine the prevalence of vanA/B/C/D-carrying enterococci, we designed a real-time PCR for vanC1/2/3 and vanD and screened rectal swabs from 360 patients. vanD was found in 27.8% of the patients, yet culture demonstrated only E. faecium from vanA-positive broths and E. gallinarum from vanC1-positive broths. No vanD-positive VRE were found, limiting the possibility of nosocomial spread of this VRE. Moreover, the high prevalence of non-VRE vanD in rectal swabs makes it unfeasible to include the vanD PCR in our VRE screening. However, having validated the vanC1/2/3 and vanD PCRs allows us to rapidly check future vanA/B-negative VRE for the presence of vanC and vanD genes.

Citrate, low pH and amino acid limitation induce citrate utilization in Lactococcus lactis biovar diacetylactis.

In Lactococcus lactis subsp. lactis biovar diacetylactis, citrate transport is facilitated by the plasmid-encoded citrate permease (CitP). In this work, we analysed regulation of citrate utilization by pH, nutrient limitation and the presence of citrate at four different levels: (i) plasmid copy number, (ii) citP transcription, (iii) citP mRNA processing and (iv) citrate utilization capacity. Citrate was supplied as cosubstrate together with lactose. The citP gene is known to be induced in cells grown at low pH. However, we demonstrated that transcription of citP was even higher in the presence of citrate (3.8-fold compared with 2.0-fold). The effect of citrate has been overlooked by other researchers because they determined the effect of citrate using M17 medium, which already contains 0.80 ± 0.07 mM citrate. The plasmid copy number increased in cells grown under amino acid limitation (1.6-fold) and/or at low pH (1.4-fold). No significant differences in citP mRNA processing were found. Citrate utilization rates increased from approximately 1 to 65 µmol min-1  gDW-1 from lowest to highest citP expression. Acetoin formation increased during growth in an acidic environment due to induction of the acetoin pathway. Quantification of the relative contributions allowed us to construct a model for regulation of citrate utilization in L. lactis biovar diacetylactis. This knowledge will help to select conditions to improve flavour formation from citrate.