A core-genome multilocus sequence typing scheme for the detection of genetically related Streptococcus pyogenes clusters

The recently observed increase in invasive Streptococcus pyogenes infections causes concern in Europe. However, conventional molecular typing methods lack discriminatory power to aid investigations of outbreaks caused by S. pyogenes. Therefore, there is an urgent need for high-resolution molecular typing methods to assess genetic relatedness between S. pyogenes isolates. In the current study, we aimed to develop a novel high-resolution core-genome multilocus sequence typing (cgMLST) scheme for S. pyogenes and compared its discriminatory power to conventional molecular typing methods. The cgMLST scheme was designed with the commercial Ridom SeqSphere+ software package. To define a cluster threshold, the scheme was evaluated using publicly available data from nine defined S. pyogenes outbreaks in the United Kingdom. The cgMLST scheme was then applied to 23 isolates from a suspected S. pyogenes outbreak and 117 S. pyogenes surveillance isolates both from the Netherlands. MLST and emm-typing results were used for comparison to cgMLST results. The allelic differences between isolates from defined outbreaks ranged between 6 and 31 for isolates with the same emm-type, resulting in a proposed cluster threshold of <5 allelic differences out of 1,095 target loci. Seven out of twenty-three (30%) isolates from the suspected outbreak had an allelic difference of <2, thereby identifying a potential cluster that could not be linked to other isolates. The proposed cgMLST scheme shows a higher discriminatory ability when compared to conventional typing methods. The rapid and simple analysis workflow allows for extended detection of clusters of potential outbreak isolates and surveillance and may facilitate the sharing of sequencing results between (inter)national laboratories.

A comparison of five Illumina, Ion Torrent, and nanopore sequencing technology-based approaches for whole genome sequencing of SARS-CoV-2

Rapid identification of the rise and spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern remains critical for monitoring of the efficacy of diagnostics, therapeutics, vaccines, and control strategies. A wide range of SARS-CoV-2 next-generation sequencing (NGS) methods have been developed over the last years, but cross-sequence technology benchmarking studies have been scarce. In the current study, 26 clinical samples were sequenced using five protocols: AmpliSeq SARS-CoV-2 (Illumina), EasySeq RC-PCR SARS-CoV-2 (Illumina/NimaGen), Ion AmpliSeq SARS-CoV-2 (Thermo Fisher), custom primer sets (Oxford Nanopore Technologies (ONT)), and capture probe-based viral metagenomics (Roche/Illumina). Studied parameters included genome coverage, depth of coverage, amplicon distribution, and variant calling. The median SARS-CoV-2 genome coverage of samples with cycle threshold (Ct) values of 30 and lower ranged from 81.6 to 99.8% for, respectively, the ONT protocol and Illumina AmpliSeq protocol. Correlation of coverage with PCR Ct values varied per protocol. Amplicon distribution signatures differed across the methods, with peak differences of up to 4 log₁₀ at disbalanced positions in samples with high viral loads (Ct values ≤ 23). Phylogenetic analyses of consensus sequences showed clustering independent of the workflow used. The proportion of SARS-CoV-2 reads in relation to background sequences, as a (cost-)efficiency metric, was the highest for the EasySeq protocol. The hands-on time was the lowest when using EasySeq and ONT protocols, with the latter additionally having the shortest sequence runtime. In conclusion, the studied protocols differed on a variety of the studied metrics. This study provides data that assist laboratories when selecting protocols for their specific setting.

Panton-Valentine Leukocidin-Positive CC398 MRSA in Urban Clinical Settings, the Netherlands

We report detection of Panton-Valentine leukocidin-positive clonal complex 398 human-origin methicillin-resistant Staphylococcus aureus L2 in the Netherlands. This hypervirulent lineage originated in the Asia-Pacific Region and could become community-acquired in Europe after recurrent travel-related introductions. Genomic surveillance enables early detection to guide control measures and help limit pathogen spread in urban settings.

Longitudinal Monitoring of DNA Viral Loads in Transplant Patients Using Quantitative Metagenomic Next-Generation Sequencing

Introduction: Immunocompromised patients are prone to reactivations and (re-)infections of multiple DNA viruses. Viral load monitoring by single-target quantitative PCRs (qPCR) is the current cornerstone for virus quantification. In this study, a metagenomic next-generation sequencing (mNGS) approach was used for the identification and load monitoring of transplantation-related DNA viruses. Methods: Longitudinal plasma samples from six patients that were qPCR-positive for cytomegalovirus (CMV), Epstein-Barr virus (EBV), BK polyomavirus (BKV), adenovirus (ADV), parvovirus B19 (B19V), and torque teno-virus (TTV) were sequenced using the quantitative metagenomic Galileo Viral Panel Solution (Arc Bio, LLC, Cambridge, MA, USA) reagents and bioinformatics pipeline combination. Qualitative and quantitative performance was analysed with a focus on viral load ranges relevant for clinical decision making. Results: All pathogens identified by qPCR were also identified by mNGS. BKV, CMV, and HHV6B were additionally detected by mNGS, and could be confirmed by qPCR or auxiliary bioinformatic analysis. Viral loads determined by mNGS correlated with the qPCR results, with inter-method differences in viral load per virus ranging from 0.19 log₁₀ IU/mL for EBV to 0.90 log₁₀ copies/mL for ADV. TTV, analysed by mNGS in a semi-quantitative way, demonstrated a mean difference of 3.0 log₁₀ copies/mL. Trends over time in viral load determined by mNGS and qPCR were comparable, and clinical thresholds for initiation of treatment were equally identified by mNGS. Conclusions: The Galileo Viral Panel for quantitative mNGS performed comparably to qPCR concerning detection and viral load determination, within clinically relevant ranges of patient management algorithms.

Viral metagenomic sequencing in a cohort of international travellers returning with febrile illness

BACKGROUND

Diagnosis of infections in returning international travellers can be challenging because of the broad spectrum of potential infectious etiologies potentially involved. Viral metagenomic next-generation sequencing (mNGS) has the potential to detect any virus present in a patient sample and is increasingly being used for difficult to diagnose cases. The aim of this study was to analyze the performance of mNGS for viral pathogen detection in the clinical setting of international travellers returning with febrile illness.

METHODS

Thirty-eight serum samples from international travellers returning with febrile illness and presenting at the outpatient clinic of the Leiden University Medical Center in the Netherlands in the time period 2015-2016 were selected retrospectively. Samples were processed for viral metagenomic sequencing using a probe panel capturing all known vertebrate viruses. Bioinformatic analysis was performed using Genome Detective software for metagenomic virus detection. Metagenomic virus findings were compared with viral pathogen detection using conventional methods.

RESULTS

In 8 out of the 38 patients (21%), a pathogenic virus was detected by mNGS. All viral pathogens detected by conventional assays were also detected by mNGS: dengue virus (n=4 patients), Epstein-Barr virus (n=2), hepatitis B virus (n=1). In addition, mNGS resulted in additional pathogenic findings in 2 patients (5%): dengue virus (n=1), and hepatitis C virus (n=1). Non-pathogenic viruses detected were: GB virus C (n=1) and torque teno viruses (n=3). High genome coverage and depth using capture probes enabled typing of the dengue viruses detected.

CONCLUSIONS

Viral metagenomics has the potential to assist the detection of viral pathogens and co-infections in one step in international travellers with a febrile syndrome. Furthermore, viral enrichment by probes resulted in high genome coverage and depth which enabled dengue virus typing.

Hospital transmission of borderline oxacillin-resistant Staphylococcus aureus evaluated by whole-genome sequencing

Introduction

Staphylococcus aureus is a major cause of hospital infections worldwide. Awareness towards methicillin-resistant S. aureus (MRSA) infections is high but attention towards borderline oxacillin-resistant S. aureus (BORSA) is limited, possibly due to an underestimated clinical relevance, presumption of low incidence and diagnostic limitations.

Gap statement

BORSA surveillance has not been routinely implemented, and thus consensus with regard to a definition and infection control measures is lacking.

Aim

Our goals were to investigate the occurrence, molecular characteristics and clinical manifestations of BORSA infections in the hospital setting.

Methodology

Following an increased incidence in 2016, BORSA cases in 2014/2016 (in our institution) were more specifically evaluated. Medical records were reviewed to investigate epidemiological links, clinical characteristics and outcomes. Resistance and virulence markers were assessed by whole genome sequencing (WGS). Conventional methods: amplified fragment length polymorphism (AFLP) ; multilocus sequence typing (MLST) and multiple locus variable-number tandem repeat analysis (MLVA) were compared with core genome MLST (cgMLST) and whole-genome single nucleotide polymorphism (wgSNP) analysis to confirm genetic clusters.

Results

From 2009 to 2013, BORSA comprised 0.1 % of all clinical S. aureus strains. In 2016, the incidence was six-fold higher in comparison to the baseline. Whole-genome SNP and cgMLST confirmed two BORSA clusters among patients with dermatological conditions. Patients with BORSA presented with skin infections, and one case developed a severe invasive infection with a fatal outcome. Infection control measures successfully prevented further transmission in both clusters. WGS findings showed that BORSA strains carried multiple resistance and virulence genes with increased pathogenic potential.

Conclusion

WGS and cgMLST effectively characterized and confirmed BORSA clusters among at-risk patients with clinical manifestations ranging from mild skin infections to life-threatening bacteraemia. Clinical awareness and active monitoring are therefore warranted for the timely implementation of infection control measures to prevent BORSA transmission in high-risk patients.

Ceftazidime-avibactam resistance and restoration of carbapenem susceptibility in KPC-producing Klebsiella pneumoniae infections: A case series

OBJECTIVES

Since the introduction of the β-lactam/β-lactamase inhibitor ceftazidime-avibactam (CZA), rapid evolution of resistance has been reported in different KPC-producing Klebsiella pneumoniae isolates. In this multicenter retrospective study, we describe the emergence of CZA resistance and evaluate the mutations that might be responsible for the restoration of carbapenem susceptibility.

METHODS

During a study period of 18 months, KPC-producing K. pneumoniae isolates of five hospitalized patients were collected with phenotypic development of CZA resistance.

RESULTS

In vitro restoration of carbapenem susceptibility during treatment was observed in 3 isolates. Whole genome sequencing of these isolates showed a D179Y mutation in the KPC gene of 2 variants and a KPC-2 with a Δ242-GT-243 deletion (KPC-14). Two KPC-3 variants showed CZA resistance with sustained carbapenemase activity without genomic adaptations in the KPC gene.

CONCLUSIONS

This study confirms the emergence of CZA resistance in KPC K. pneumoniae. The role of carbapenems in treating patients with these variants is unclear and combination therapies warrant further investigation.

Borderline oxacillin-resistant Staphylococcus aureus carriage among healthcare workers at neonatal intensive care unit and paediatric ward

BACKGROUND

During a meticillin-resistant Staphylococcus aureus contact tracing and screening investigation, two borderline oxacillin-resistant Staphylococcus aureus (BORSA)-positive screening cultures were encountered among neonatal intensive care unit (NICU) healthcare workers (HCWs). This finding led to further investigations.

AIM

To assess the likelihood of an outbreak with direct transmission among HCWs.

METHODS

An infection control team was initiated after the discovery. The team initiated additional infection control measures and evaluated new findings. All NICUs and paediatric ward HCWs were screened for BORSA carriage, and a prospective BORSA seven-week monitoring period for patients was observed. To assess the likelihood of an outbreak with direct transmission among HCWs, the BORSA isolates were analysed using augmented fragment length polymorphism and whole-genome sequencing (WGS).

FINDINGS

Positive HCWs were prohibited from clinical work while awaiting the results from the screening programme. In all, 127 NICU and 77 general paediatric ward HCWs were screened for BORSA carriage; five HCWs were BORSA positive. Seventy-two patients were screened during the seven-week period yielding a total of 138 cultures, ranging from one to nine cultures per patient. No spread from HCWs to patients occurred, and the BORSA screening programme was discontinued. WGS analysis with core genome multi-locus sequence typing of all five BORSA strains showed relatedness between two NICU strains.

CONCLUSION

During a seven-week period, no transmission from BORSA-positive HCWs to neonates was observed in either screening or clinical cultures. More vigilance and experience is needed to design adequate evidence-based interventions in the future for this vulnerable population.

Taxonomic classification and abundance estimation using 16S and WGS-A comparison using controlled reference samples

The assessment of microbiome biodiversity is the most common application of metagenomics. While 16S sequencing remains standard procedure for taxonomic profiling of metagenomic data, a growing number of studies have clearly demonstrated biases associated with this method. By using Whole Genome Shotgun sequencing (WGS) metagenomics, most of the known restrictions associated with 16S data are alleviated. However, due to the computationally intensive data analyses and higher sequencing costs, WGS based metagenomics remains a less popular option. Selecting the experiment type that provides a comprehensive, yet manageable amount of information is a challenge encountered in many metagenomics studies. In this work, we created a series of artificial bacterial mixes, each with a different distribution of skin-associated microbial species. These mixes were used to estimate the resolution of two different metagenomic experiments - 16S and WGS - and to evaluate several different bioinformatics approaches for taxonomic read classification. In all test cases, WGS approaches provide much more accurate results, in terms of taxa prediction and abundance estimation, in comparison to those of 16S. Furthermore, we demonstrate that a 16S dataset, analysed using different state of the art techniques and reference databases, can produce widely different results. In light of the fact that most forensic metagenomic analysis are still performed using 16S data, our results are especially important.

Spread of ESBL-producing Escherichia coli in nursing home residents in Ireland and the Netherlands may reflect infrastructural differences

A prevalence study in two nursing homes (one each in the Netherlands and Ireland) found four (11%) Dutch and six (9%) Irish residents colonized with 11 extended-spectrum beta-lactamase-producing Escherichia coli, 10 of which contained CTX-M-15. Four Dutch isolates, from three residents of the same ward, belonged to E. coli O25:H4, sequence type (ST) 131 and were part of the same cluster type by whole-genome sequencing. Four Irish residents on three different wards were colonized with an identical E. coli O89:H9, ST131, complex type 1478. Cross-transmission between three Irish wards may reflect differences in nursing home infrastructure, specifically communal areas and multi-bedded resident rooms.

BacTag - a pipeline for fast and accurate gene and allele typing in bacterial sequencing data based on database preprocessing

BACKGROUND

Bacteria carry a wide array of genes, some of which have multiple alleles. These different alleles are often responsible for distinct types of virulence and can determine the classification at the subspecies levels (e.g., housekeeping genes for Multi Locus Sequence Typing, MLST). Therefore, it is important to rapidly detect not only the gene of interest, but also the relevant allele. Current sequencing-based methods are limited to mapping reads to each of the known allele reference, which is a time-consuming procedure.

RESULTS

To address this limitation, we developed BacTag - a pipeline that rapidly and accurately detects which genes are present in a sequencing dataset and reports the allele of each of the identified genes. We exploit the fact that different alleles of the same gene have a high similarity. Instead of mapping the reads to each of the allele reference sequences, we preprocess the database prior to the analysis, which makes the subsequent gene and allele identification efficient. During the preprocessing, we determine a representative reference sequence for each gene and store the differences between all alleles and this chosen reference. Throughout the analysis we estimate whether the gene is present in the sequencing data by mapping the reads to this reference sequence; if the gene is found, we compare the variants to those in the preprocessed database. This allows to detect which specific allele is present in the sequencing data. Our pipeline was successfully tested on artificial WGS E. coli, S. pseudintermedius, P. gingivalis, M. bovis, Borrelia spp. and Streptomyces spp. data and real WGS E. coli and K. pneumoniae data in order to report alleles of MLST house-keeping genes.

CONCLUSIONS

We developed a new pipeline for fast and accurate gene and allele recognition based on database preprocessing and parallel computing and performed better or comparable to the current popular tools. We believe that our approach can be useful for a wide range of projects, including bacterial subspecies classification, clinical diagnostics of bacterial infections, and epidemiological studies.

Direct detection of extended-spectrum beta-lactamases (CTX-M) from blood cultures by LC-MS/MS bottom-up proteomics

Rapid bacterial species identification and antibiotic susceptibility testing in positive blood cultures have an important impact on the antibiotic treatment for patients. To identify extended-spectrum beta-lactamases (ESBL) directly in positive blood culture bottles, we developed a workflow of saponin extraction followed by a bottom-up proteomics approach using liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS). The workflow was applied to positive blood cultures with Escherichia coli and Klebsiella pneumoniae collected prospectively in two academic hospitals over a 4-month period. Of 170 positive blood cultures, 22 (12.9%) contained ESBL-positive isolates based on standard susceptibility testing. Proteomic analysis identified CTX-M ESBLs in 95% of these isolates directly in positive blood cultures, whereas no false positives were found in the non-ESBL producing positive blood cultures. The results were confirmed by molecular characterisation of beta-lactamase genes. Based on this proof-of-concept study, we conclude that LC-MS/MS-based protein analysis can directly identify extended-spectrum beta lactamases in E. coli and K. pneumoniae positive blood cultures, and could be further developed for application in routine diagnostics.

Typing Pseudomonas aeruginosa Isolates with Ultrahigh Resolution MALDI-FTICR Mass Spectrometry

The introduction of standardized matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) platforms in the medical microbiological practice has revolutionized the way microbial species identification is performed on a daily basis. To a large extent, this is due to the ease of operation. Acquired spectra are compared to profiles obtained from cultured colonies present in a reference spectra database. It is fast and reliable, and costs are low compared to previous diagnostic approaches. However, the low resolution and dynamic range of the MALDI-TOF profiles have shown limited applicability for the discrimination of different bacterial strains, as achieved with typing based on genetic markers. This is pivotal in cases where certain strains are associated with, e.g., virulence or antibiotic resistance. Ultrahigh resolution MALDI-FTICR MS allows the measurement of small proteins at isotopic resolution and can be used to analyze complex mixtures with increased dynamic range and higher precision than MALDI-TOF MS, while still generating results in a similar time frame. Here, we propose to use ultrahigh resolution 15T MALDI-Fourier transform ion cyclotron resonance (FTICR) MS to discriminate clinically relevant bacterial strains after species identification performed by MALDI-TOF MS. We used a collection of well characterized Pseudomonas aeruginosa strains, featuring distinct antibiotic resistance profiles, and isolates obtained during hospital outbreaks. Following cluster analysis based on amplification fragment length polymorphism (AFLP), these strains were grouped into three different clusters. The same clusters were obtained using protein profiles generated by MALDI-FTICR MS. Subsequent intact protein analysis by electrospray ionization (ESI)-collision-induced dissociation (CID)-FTICR MS was applied to identify protein isoforms that contribute to the separation of the different clusters, illustrating the additional advantage of this analytical platform.

Extended-spectrum β-lactamase-producing enterobacteriaceae among travelers from the Netherlands

A prospective cohort study was performed among travelers from the Netherlands to investigate the acquisition of carbapenemase-producing Enterobacteriaceae (CP-E) and extended-spectrum β-lactamase–producing Enterobacteriaceae (ESBL-E) and associated risk factors. Questionnaires were administered and rectal swab samples were collected and tested before and after traveler return. Of 370 travelers, 32 (8.6%) were colonized with ESBL-E before trave,; 113 (30.5%) acquired an ESBL-E during travel, and 26 were still colonized 6 months after return. No CP-E were found. Independent risk factors for ESBL-E acquisition were travel to South and East Asia. Multilocus sequence typing showed extensive genetic diversity among Escherichia coli. Predominant ESBLs were CTX-M enzymes. The acquisition rate, 30.5%, of ESBL-E in travelers from the Netherlands to all destinations studied was high. Active surveillance for ESBL-E and CP-E and contact isolation precautions may be recommended at admission to medical facilities for patients who traveled to Asia during the previous 6 months.

Exploring the contribution of efflux on the resistance to fluoroquinolones in clinical isolates of Escherichia coli

Resistance to ciprofloxacin in Escherichia coli is increasing parallel to increased use of fluoroquinolones both in The Netherlands and in other European countries. The objective was to investigate the contribution of active efflux and expression of outer membrane proteins (OMPs) in a collection of clinical E. coli isolates collected at a clinical microbiology department in a Dutch hospital. Forty-seven E. coli isolates a wide range of ciprofloxacin minimum inhibitory concentrations and known mutations in the quinolone resistance determining region were included. A fluorometric determination of bisbenzimide efflux was used two different efflux pump inhibitors and compared to quantitative reverse transcription-polymerase chain reaction (qRT-PCR) for the expression levels of acrA, acrB, tolC, yhiV, and mdfA efflux pump genes and the OMPs ompF and ompX. Six isolates (12.7%) showed increased efflux. Although in 35 isolates (76%), overexpression of ≥1 efflux pump genes using qRT-PCR was present. Only the combined overexpression of acrAB-TolC and mdfA correlated with the phenotypic efflux assay using glucose/carbonyl cyanide m-chlorophenylhydrazone with glucose. Thus, efflux was involved in ciprofloxacin resistance in a limited number of E. coli isolates collected at a clinical microbiology department in a Dutch hospital complementing other resistance mechanisms.