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 log10 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 spread of pathogens 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 log10 IU/mL for EBV to 0.90 log10 copies/mL for ADV. TTV, analysed by mNGS in a semi-quantitative way, demonstrated a mean difference of 3.0 log10 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.

Rhinovirus Detection in the Nasopharynx of Children Undergoing Cardiac Surgery Is Not Associated With Longer PICU Length of Stay: Results of the Impact of Rhinovirus Infection After Cardiac Surgery in Kids (RISK) Study.

OBJECTIVES: To determine whether children with asymptomatic carriage of rhinovirus in the nasopharynx before elective cardiac surgery have an increased risk of prolonged PICU length of stay. STUDY DESIGN: Prospective, single-center, blinded observational cohort study. SETTING: PICU in a tertiary hospital in The Netherlands. PATIENTS: Children under 12 years old undergoing elective cardiac surgery were enrolled in the study after informed consent of the parents/guardians. INTERVENTIONS: The parents/guardians filled out a questionnaire regarding respiratory symptoms. On the day of the operation, a nasopharyngeal swab was obtained. Clinical data were collected during PICU admission, and PICU/hospital length of stay were reported. If a patient was still intubated 3 days after operation, an additional nasopharyngeal swab was collected. Nasopharyngeal swabs were tested for rhinovirus and other respiratory viruses with polymerase chain reaction. MEASUREMENTS AND MAIN RESULTS: Of the 163 included children, 74 (45%) tested rhinovirus positive. Rhinovirus-positive patients did not have a prolonged PICU length of stay (median 2 d each; p = 0.257). Rhinovirus-positive patients had a significantly shorter median hospital length of stay compared with rhinovirus-negative patients (8 vs 9 d, respectively; p = 0.006). Overall, 97 of the patients (60%) tested positive for one or more respiratory virus. Virus-positive patients had significantly shorter PICU and hospital length of stay, ventilatory support, and nonmechanical ventilation. Virus-negative patients had respiratory symptoms suspected for a respiratory infection more often. In 31% of the children, the parents reported mild upper respiratory complaints a day prior to the cardiac surgery, this was associated with postextubation stridor, but no other clinical outcome measures. CONCLUSIONS: Preoperative rhinovirus polymerase chain reaction positivity is not associated with prolonged PICU length of stay. Our findings do not support the use of routine polymerase chain reaction testing for respiratory viruses in asymptomatic children admitted for elective cardiac surgery.

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.

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.

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.

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.

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.