Framing Bacterial Genomics for Public Health(care)

Development and implementation of a core genome multilocus sequence typing scheme for Yersinia enterocolitica: a tool for surveillance and outbreak detection

Yersinia enterocolitica (Y. enterocolitica) is the most frequent etiological agent of yersiniosis and has been responsible for several national outbreaks in Norway and elsewhere. A standardized high-resolution method, such as core genome Multilocus Sequence Typing (cgMLST), is needed for pathogen traceability at the national and international levels. In this study, we developed and implemented a cgMLST scheme for Y. enterocolitica. We designed a cgMLST scheme in SeqSphere + using high-quality genomes from different Y. enterocolitica biotype sublineages. The scheme was validated if more than 95% of targets were found across all tested Y. enterocolitica: 563 Norwegian genomes collected between 2012 and 2022 and 327 genomes from public data sets. We applied the scheme to known outbreaks to establish a threshold for identifying major complex types (CTs) based on the number of allelic differences. The final cgMLST scheme included 2,582 genes with a median of 97.9% (interquartile range 97.6%-98.8%) targets found across all tested genomes. Analysis of outbreaks identified all outbreak strains using single linkage clustering at four allelic differences. This threshold identified 311 unique CTs in Norway, of which CT18, CT12, and CT5 were identified as the most frequently associated with outbreaks. The cgMLST scheme showed a very good performance in typing Y. enterocolitica using diverse data sources and was able to identify outbreak clusters. We recommend the implementation of this scheme nationally and internationally to facilitate Y. enterocolitica surveillance and improve outbreak response in national and cross-border outbreaks.

Genomics for public health and international surveillance of antimicrobial resistance

Historically, epidemiological investigation and surveillance for bacterial antimicrobial resistance (AMR) has relied on low-resolution isolate-based phenotypic analyses undertaken at local and national reference laboratories. Genomic sequencing has the potential to provide a far more high-resolution picture of AMR evolution and transmission, and is already beginning to revolutionise how public health surveillance networks monitor and tackle bacterial AMR. However, the routine integration of genomics in surveillance pipelines still has considerable barriers to overcome. In 2022, a workshop series and online consultation brought together international experts in AMR and pathogen genomics to assess the status of genomic applications for AMR surveillance in a range of settings. Here we focus on discussions around the use of genomics for public health and international AMR surveillance, noting the potential advantages of, and barriers to, implementation, and proposing recommendations from the working group to help to drive the adoption of genomics in public health AMR surveillance. These recommendations include the need to build capacity for genome sequencing and analysis, harmonising and standardising surveillance systems, developing equitable data sharing and governance frameworks, and strengthening interactions and relationships among stakeholders at multiple levels.

Comparison of carbapenem-susceptible and carbapenem-resistant Enterobacterales at nine sites in the USA, 2013-2016: a resource for antimicrobial resistance investigators

Carbapenem-resistant Enterobacterales (CRE) are an urgent public health threat. Genomic sequencing is an important tool for investigating CRE. Through the Division of Healthcare Quality Promotion Sentinel Surveillance system, we collected CRE and carbapenem-susceptible Enterobacterales (CSE) from nine clinical laboratories in the USA from 2013 to 2016 and analysed both phenotypic and genomic sequencing data for 680 isolates. We describe the molecular epidemiology and antimicrobial susceptibility testing (AST) data of this collection of isolates. We also performed a phenotype-genotype correlation for the carbapenems and evaluated the presence of virulence genes in Klebsiella pneumoniae complex isolates. These AST and genomic sequencing data can be used to compare and contrast CRE and CSE at these sites and serve as a resource for the antimicrobial resistance research community.

Antibiotic resistomes in face-mask biofilm along an urban river: Multiple drivers and co-occurrence with human opportunistic pathogens

Discarded face masks from the global COVID-19 pandemic have contributed significantly to plastic pollution in surface water, whereas their potential as a reservoir for aquatic pollutants is not well understood. Herein, we conducted a field experiment along a human-impacted urban river, investigating the variations of antibiotic resistance genes (ARGs), pathogens, and water-borne contaminants in commonly-used face masks. Results showed that high-biomass biofilms formed on face masks selectively enriched more ARGs than stone biofilm (0.08-0.22 vs 0.07-0.15 copies/16 S rRNA gene copies) from bulk water, which mainly due to unique microbial communities, enhanced horizontal gene transfer, and selective pressure of accumulated contaminants based on redundancy analysis and variation partitioning analysis. Several human opportunistic pathogens (e.g., Acinetobacter, Escherichia-Shigella, Bacillus, and Klebsiella), which are considered potential ARG carriers, were also greatly concentrated in face-mask biofilms, imposing a potential threat to aquatic ecological environment and human health. Moreover, wastewater treatment plant effluents, as an important source of pollutants to urban rivers, further aggravated the abundances of ARGs and opportunistic pathogens in face-mask biofilms. Our findings demonstrated that discarded face masks provide a hotspot for the proliferation and spread of ARGs and pathogens in urban water, highlighting the urgent requirement for implementing stricter regulations in face mask disposal.

Validation and Application of Long-Read Whole-Genome Sequencing for Antimicrobial Resistance Gene Detection and Antimicrobial Susceptibility Testing

Next-generation sequencing applications are increasingly used for detection and characterization of antimicrobial-resistant pathogens in clinical settings. Oxford Nanopore Technologies (ONT) sequencing offers advantages for clinical use compared with other sequencing methodologies because it enables real-time basecalling, produces long sequencing reads that increase the ability to correctly assemble DNA fragments, provides short turnaround times, and requires relatively uncomplicated sample preparation. A drawback of ONT sequencing, however, is its lower per-read accuracy than short-read sequencing. We sought to identify best practices in ONT sequencing protocols. As some variability in sequencing results may be introduced by the DNA extraction methodology, we tested three DNA extraction kits across three independent laboratories using a representative set of six bacterial isolates to investigate accuracy and reproducibility of ONT technology. All DNA extraction techniques showed comparable performance; however, the DNeasy PowerSoil Pro kit had the highest sequencing yield. This kit was subsequently applied to 42 sequentially collected bacterial isolates from blood cultures to assess Ares Genetics's pipelines for predictive whole-genome sequencing antimicrobial susceptibility testing (WGS-AST) performance compared to phenotypic triplicate broth microdilution results. WGS-AST results ranged across the organisms and resulted in an overall categorical agreement of 95% for penicillins, 82.4% for cephalosporins, 76.7% for carbapenems, 86.9% for fluoroquinolones, and 96.2% for aminoglycosides. Very major errors/major errors were 0%/16.7% (penicillins), 11.7%/3.6% (cephalosporins), 0%/24.4% (carbapenems), 2.5%/7.7% (fluoroquinolones), and 0%/4.1% (aminoglycosides), respectively. This work showed that, although additional refinements are necessary, ONT sequencing demonstrates potential as a method to perform WGS-AST on cultured isolates for patient care.

PRAWNS: compact pan-genomic features for whole-genome population genomics

MOTIVATION

Scientists seeking to understand the genomic basis of bacterial phenotypes, such as antibiotic resistance, today have access to an unprecedented number of complete and nearly complete genomes. Making sense of these data requires computational tools able to perform multiple-genome comparisons efficiently, yet currently available tools cannot scale beyond several tens of genomes.

RESULTS

We describe PRAWNS, an efficient and scalable tool for multiple-genome analysis. PRAWNS defines a concise set of genomic features (metablocks), as well as pairwise relationships between them, which can be used as a basis for large-scale genotype-phenotype association studies. We demonstrate the effectiveness of PRAWNS by identifying genomic regions associated with antibiotic resistance in Acinetobacter baumannii.

AVAILABILITY AND IMPLEMENTATION

PRAWNS is implemented in C++ and Python3, licensed under the GPLv3 license, and freely downloadable from GitHub (https://github.com/KiranJavkar/PRAWNS.git).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Molecular Epidemiology of Carbapenem-Resistant Acinetobacter baumannii in the United States, 2013-2017

Healthcare-associated carbapenem-resistant Acinetobacter baumannii (CRAB) infections are a serious threat associated with global epidemic clones and a variety of carbapenemase gene classes. In this study, we describe the molecular epidemiology, including whole-genome sequencing analysis and antimicrobial susceptibility profiles of 92 selected, nonredundant CRAB collected through public health efforts in the United States from 2013 to 2017. Among the 92 isolates, the Oxford (OX) multilocus sequence typing scheme identified 30 sequence types (STs); the majority of isolates (n = 59, 64%) represented STs belonging to the international clonal complex 92 (CC92_OX). Among these, ST208_OX (n = 21) and ST281_OX (n = 20) were the most common. All isolates carried an OXA-type carbapenemase gene, comprising 20 alleles. Ninety isolates (98%) encoded an intrinsic OXA-51-like enzyme; 67 (73%) harbored an additional acquired bla_OXA gene, most commonly bla_OXA-23 (n = 45; 49%). Compared with isolates harboring only intrinsic oxacillinase genes, acquired bla_OXA gene presence was associated with higher prevalence of resistance and a higher median minimum inhibitory concentration to the carbapenem imipenem (64 μg/mL vs. 8 μg/mL), and antibiotics from other drug classes, including penicillin, aminoglycosides, cephalosporins, and polymyxins. These data illustrate the wide distribution of CC92_OX and high prevalence of acquired bla_OXA carbapenemase genes among CRAB in the United States.