Bacterial Whole-genome Sequencing Research Articles

BakRep - a searchable large-scale web repository for bacterial genomes, characterizations and metadata.

Bacteria are fascinating research objects in many disciplines for countless reasons, and whole-genome sequencing (WGS) has become the paramount methodology to advance our microbiological understanding. Meanwhile, access to cost-effective sequencing platforms has accelerated bacterial WGS to unprecedented levels, introducing new challenges in terms of data accessibility, computational demands, heterogeneity of analysis workflows and, thus, ultimately its scientific usability. To this end, a previous study released a uniformly processed set of 661 405 bacterial genome assemblies obtained from the European Nucleotide Archive as of November 2018. Building on these accomplishments, we conducted further genome-based analyses like taxonomic classification, multilocus sequence typing and annotation of all genomes. Here, we present BakRep, a searchable large-scale web repository of these genomes enriched with consistent genome characterizations and original metadata. The platform provides a flexible search engine combining taxonomic, genomic and metadata information, as well as interactive elements to visualize genomic features. Furthermore, all results can be downloaded for offline analyses via an accompanying command line tool. The web repository is accessible via https://bakrep.computational.bio.

Easy analysis of bacterial whole-genome sequencing data for clinical microbiologists using open-source Galaxy platform: Characterization of ESBL-producing Enterobacterales from bloodstream infections

ObjectivesClinical microbiologists require easy-to-use open access tools with graphical interfaces to perform bacterial whole-genome sequencing (WGS) in routine practice. This study aimed to build a bioinformatics pipeline on the open-source Galaxy platform, facilitating comprehensive and reproducible analysis of bacterial WGS data in a few steps. We then used it to characterize our local epidemiology of ESBL-producing Enterobacterales isolated from patients with bacteremia. MethodsWe built a bioinformatics pipeline consisting of the following sequential tools: Fastp (input data trimming); FastQC (read quality control); SPAdes (genome assembly); Quast (quality control of genome assembly); Prokka (gene annotation); Staramr (ResFinder database) and ABRicate (CARD database) for antimicrobial resistance (AMR) gene screening and molecular strain typing. Paired-end short read WGS data from all ESBL-producing Enterobacterales strains isolated from patients with bacteremia over one year were analysed. ResultsThe Galaxy platform does not require command line tools. The bioinformatics pipeline was constructed within one hour. It only required uploading fastq files and facilitated systematization of de novo assembly of genomes, MLST typing, and AMR gene screening in one step. Among the 66 ESBL-producing strains analysed, the two most frequent ESBL genes were blaCTX−M-15 (62.1%) and blaCTX−M-27 (13.6%). ConclusionsThe open-access Galaxy platform provides a graphical interface and easy-to-use tools suitable for routine use in clinical microbiology laboratories without bioinformatics specialists. We believe that this platform will facilitate fast and low-cost analysis of bacterial WGS data, especially in resource-limited settings.

Whole-genome sequencing of multidrug-resistant Escherichia coli causing urinary tract infection in an immunocompromised patient: a case report

BackgroundEscherichia coli is a major human pathogen responsible for a broad range of clinical illnesses. It has been linked to endemic and epidemic nosocomial diseases caused by multidrug-resistant pathogens in Sudan as well as throughout the globe.Case presentationA 76-year-old African woman arrived at Saad Rashwan Medical Centre complaining of backaches and discomfort during urination. Throughout the preceding 5 years, the patient had recurrent urinary tract infections. Following overnight incubation at 37 °C, Escherichia coli was found in her midstream urine specimen on cysteine lactose electrolyte deficient agar media. Minimum inhibitory concentration (colorimetric/turbidimetric method) was employed to test a wide range of antimicrobial drugs against this bacterial strain, and the results revealed significant multidrug resistance. QIAamp® DNA Mini Kit was used to obtain DNA Template from the purified Escherichia coli (Qiagen, Hilden, Germany). The bacterial whole-genome sequence was done by Novogene company (Hong Kong) using Illumina HiSeq 2500 (Illumina, San Diego, CA, USA), followed by whole genome reconstructions, and identification of antibiotic-resistant genes. Phylogenetic analysis revealed that our strain was related to the Escherichia coli DSM30083 ( genome sequence ID: CP033092.2) from the USA. Our strain possessed the following antimicrobial-resistant genes: aminoglycoside (kdpE, baeR, cpxA, aadA5), nitroimidazole (msbA), phosphonic acid (mdtG), tetracycline (emrY), macrolide, penam, tetracycline, (evgA, TolC, H-NS), fluoroquinolone, cephalosporin, glycylcycline, penam, tetracycline, rifamycin, phenicol antibiotic, disinfecting agents and antiseptics (acrB; marA), sulfonamide (sul1), macrolide (Mrx), cephalosporin, penam (CTX-M-15), carbapenem, cephalosporin, and penam (OXA-1).ConclusionThis study found that the isolated Escherichia coli strain had varied antimicrobial resistance genes on the basis of whole-genome sequencing and phenotypic resistance analyses. Whole-genome sequencing is critical for control and preventative methods to battle the growing threat of antimicrobial resistance. A larger investigation is recommended for improved generalization of results.

An enterococcal phage-derived enzyme suppresses graft-versus-host disease

Changes in the gut microbiome have pivotal roles in the pathogenesis of acute graft-versus-host disease (aGVHD) after allogenic haematopoietic cell transplantation (allo-HCT)1–6. However, effective methods for safely resolving gut dysbiosis have not yet been established. An expansion of the pathogen Enterococcus faecalis in the intestine, associated with dysbiosis, has been shown to be a risk factor for aGVHD7–10. Here we analyse the intestinal microbiome of patients with allo-HCT, and find that E. faecalis escapes elimination and proliferates in the intestine by forming biofilms, rather than by acquiring drug-resistance genes. We isolated cytolysin-positive highly pathogenic E. faecalis from faecal samples and identified an anti-E. faecalis enzyme derived from E. faecalis-specific bacteriophages by analysing bacterial whole-genome sequencing data. The antibacterial enzyme had lytic activity against the biofilm of E. faecalis in vitro and in vivo. Furthermore, in aGVHD-induced gnotobiotic mice that were colonized with E. faecalis or with patient faecal samples characterized by the domination of Enterococcus, levels of intestinal cytolysin-positive E.faecalis were decreased and survival was significantly increased in the group that was treated with the E. faecalis-specific enzyme, compared with controls. Thus, administration of a phage-derived antibacterial enzyme that is specific to biofilm-forming pathogenic E. faecalis—which is difficult to eliminate with existing antibiotics—might provide an approach to protect against aGVHD.

2023 Neer Award for Basic Science: Genetics of Cutibacterium acnes in revision shoulder arthroplasty: a large-scale bacterial whole-genome sequencing study

BackgroundCutibacterium acnes is the bacterium most commonly responsible for shoulder periprosthetic joint infection (PJI) and is often cultured from samples obtained at the time of revision for failed shoulder arthroplasty. We sought to determine whether these bacteria originate from the patient or from exogenous sources. We also sought to identify which C. acnes genetic traits were associated with the development of shoulder PJI. MethodsWe performed bacterial whole-genome sequencing of C. acnes from a single-institution repository of cultures obtained before or during primary and revision shoulder arthroplasty and correlated the molecular epidemiology and genetic content of strains with clinical features of infection. ResultsA total of 341 isolates collected over a four-year period from 88 patients were sequenced. C. acnes cultured from surgical specimens demonstrated significant similarity to the strains colonizing the skin of the same patient (p<0.001). Infrequently, there was evidence of strains shared across unrelated patients, suggesting that exogenous sources of C. acnes culture-positivity were uncommon. Phylotypes IB and II were modestly associated with clinical features of PJI, but all phylotypes appeared inherently capable of causing disease. Chronic shoulder PJI was associated with the absence of common C. acnes genes involved in bacterial quorum-sensing (luxS, tqsA). ConclusionC. acnes strains cultured from deep intraoperative sources during revision shoulder arthroplasty demonstrate strong genetic similarity to the strains colonizing a patient’s skin. Some phylotypes of C. acnes commonly colonizing human skin are modestly more virulent than others, but all phylotypes have a capacity for PJI. C. acnes cultured from cases of PJI commonly demonstrated genetic hallmarks associated with adaptation from acute to chronic phases of infection. This is the strongest evidence to date supporting the role of the patient’s own, cutaneous C. acnes strains in the pathogenesis of shoulder arthroplasty infection. Our findings support the importance of further research focused on perioperative decolonization and management of endogenous bacteria that are likely to be introduced into the arthroplasty wound at the time of skin incision.

Analysis of Antimicrobial Resistance in Bacterial Pathogens Recovered from Food and Human Sources: Insights from 639,087 Bacterial Whole-Genome Sequences in the NCBI Pathogen Detection Database.

Understanding the role of foods in the emergence and spread of antimicrobial resistance necessitates the initial documentation of antibiotic resistance genes within bacterial species found in foods. Here, the NCBI Pathogen Detection database was used to query antimicrobial resistance gene prevalence in foodborne and human clinical bacterial isolates. Of the 1,843,630 sequence entries, 639,087 (34.7%) were assigned to foodborne or human clinical sources with 147,788 (23.14%) from food and 427,614 (76.88%) from humans. The majority of foodborne isolates were either Salmonella (47.88%), Campylobacter (23.03%), Escherichia (11.79%), or Listeria (11.3%), and the remaining 6% belonged to 20 other genera. Most foodborne isolates were from meat/poultry (95,251 or 64.45%), followed by multi-product mixed food sources (29,892 or 20.23%) and fish/seafood (6503 or 4.4%); however, the most prominent isolation source varied depending on the genus/species. Resistance gene carriage also varied depending on isolation source and genus/species. Of note, Klebsiella pneumoniae and Enterobacter spp. carried larger proportions of the quinolone resistance gene qnrS and some clinically relevant beta-lactam resistance genes in comparison to Salmonella and Escherichia coli. The prevalence of mec in S. aureus did not significantly differ between meat/poultry and multi-product sources relative to clinical sources, whereas this resistance was rare in isolates from dairy sources. The proportion of biocide resistance in Bacillus and Escherichia was significantly higher in clinical isolates compared to many foodborne sources but significantly lower in clinical Listeria compared to foodborne Listeria. This work exposes the gaps in current publicly available sequence data repositories, which are largely composed of clinical isolates and are biased towards specific highly abundant pathogenic species. We also highlight the importance of requiring and curating metadata on sequence submission to not only ensure correct information and data interpretation but also foster efficient analysis, sharing, and collaboration. To effectively monitor resistance carriage in food production, additional work on sequencing and characterizing AMR carriage in common commensal foodborne bacteria is critical.

Hound: a novel tool for automated mapping of genotype to phenotype in bacterial genomes assembled de novo.

Increasing evidence suggests that microbial species have a strong within species genetic heterogeneity. This can be problematic for the analysis of prokaryote genomes, which commonly relies on a reference genome to guide the assembly process. Differences between reference and sample genomes will therefore introduce errors in final assembly, jeopardizing the detection from structural variations to point mutations-critical for genomic surveillance of antibiotic resistance. Here we present Hound, a pipeline that integrates publicly available tools to assemble prokaryote genomes de novo, detect user-given genes by similarity to report mutations found in the coding sequence, promoter, as well as relative gene copy number within the assembly. Importantly, Hound can use the query sequence as a guide to merge contigs, and reconstruct genes that were fragmented by the assembler. To showcase Hound, we screened through 5032 bacterial whole-genome sequences isolated from farmed animals and human infections, using the amino acid sequence encoded by blaTEM-1, to detect and predict resistance to amoxicillin/clavulanate which is driven by over-expression of this gene. We believe this tool can facilitate the analysis of prokaryote species that currently lack a reference genome, and can be scaled either up to build automated systems for genomic surveillance or down to integrate into antibiotic susceptibility point-of-care diagnostics.

Edwardsiella tarda Causing Septicemia in a Wild Crested Ibis (Nipponia nippon).

An adult Crested Ibis (Nipponia nippon) was found moribund in the Qinling area of China. Postmortem examination and histopathological analysis revealed lung inflammation and multi-organ hemorrhage. Bacterial isolation and whole-genome sequencing confirmed Edwardsiella tarda infection.

A prospective bacterial whole-genome-sequencing-based surveillance programme for comprehensive early detection of healthcare-associated infection transmission in paediatric oncology patients

Bacterial whole-genome sequencing (WGS) and determination of genetic relatedness is an important tool for investigation of epidemiologically suspected outbreaks. This prospective cohort study evaluated a comprehensive, prospective bacterial WGS-based surveillance programme for early detection of transmission of most bacterial pathogens among patients at a paediatric oncology hospital. Cultured bacterial isolates from clinical diagnostic specimens collected prospectively from both inpatient and outpatient encounters between January 2019 and December 2021 underwent routine WGS and core genome multi-locus sequence typing to determine isolates' relatedness. Previously collected isolates from January to December 2018 were retrospectively analysed for identification of prior or ongoing transmission. Multi-patient clusters were investigated to identify potential transmission events based on temporal and spatial epidemiological links and interventions were introduced. A total of 1497 bacterial isolates from 1025 patients underwent WGS. A total of 259 genetically related clusters were detected, of which 18 (6.9%) multi-patient clusters involving 38 (3.7%) patients were identified. Sixteen clusters involved two patients each, and two clusters involved three patients. Following investigation, epidemiologically plausible transmission links were identified in five (27.8%) multi-patient clusters. None of the multi-patient clusters were suspected by conventional epidemiological surveillance. Bacterial WGS-based surveillance for early detection of hospital transmission detected several limited multi-patient clusters that were unrecognized by conventional epidemiological methods. Genomic surveillance helped efficiently focus interventions while reducing unnecessary investigations.

Application of Next-Generation Sequencing to Enterobacter Hormaechei Subspecies Analysis during a Neonatal Intensive Care Unit Outbreak.

The Enterobacter cloacae complex (ECC) species are potential neonatal pathogens, and ECC strains are among the most commonly encountered Enterobacter spp. associated with nosocomial bloodstream infections. Outbreaks caused by ECC can lead to significant morbidity and mortality in susceptible neonates. At the molecular level, ECC exhibits genomic heterogeneity, with six closely related species and subspecies. Genetic variability poses a challenge in accurately identifying outbreaks by determining the clonality of ECC isolates. This difficulty is further compounded by the limitations of the commonly used molecular typing methods, such as pulsed field gel electrophoresis, which do not provide reliable accuracy in distinguishing between ECC strains and can lead to incorrect conclusions. Next-generation sequencing (NGS) offers superior resolution in determining strain relatedness. Therefore, we investigated the clinical pertinence of incorporating NGS into existing bundle measures to enhance patient management during an outbreak of ECC in a level-3 neonatal intensive care unit (NICU) in Germany. As the standard of care, all neonates on the NICU received weekly microbiological swabs (nasopharyngeal and rectal) and analysis of endotracheal secretion, where feasible. During the 2.5-month outbreak, colonisation with ECC was detected in n = 10 neonates. The phylogenetic relationship and potential antimicrobial resistance genes as well as mobile genetic elements were identified via bacterial whole-genome sequencing (WGS) using Illumina MiSeq followed by in silico data analysis. Although all ECC isolates exhibited almost identical antimicrobial susceptibility patterns, the WGS data revealed the involvement of four different ECC clones. The isolates could be characterised as Enterobacter hormaechei subspecies steigerwaltii (n = 6, clonal), subsp. hoffmannii (n = 3, two clones) and subsp. oharae (n = 1). Despite the collection of environmental samples, no source of this diffuse outbreak could be identified. A new standardised operating procedure was implemented to enhance the management of neonates colonised with MRGN. This collaborative approach involved both parents and medical professionals and successfully prevented further transmission of ECC. Initially, it was believed that the NICU outbreak was caused by a single ECC clone due to the similarity in antibiotic resistance. However, our findings show that antibiotic susceptibility patterns can be misleading in investigating outbreaks of multi-drug-resistant ECC. In contrast, bacterial WGS accurately identified ECC at the clonal level, which significantly helped to delineate the nature of the observed outbreak.

CamPype: an open-source workflow for automated bacterial whole-genome sequencing analysis focused on Campylobacter

BackgroundThe rapid expansion of Whole-Genome Sequencing has revolutionized the fields of clinical and food microbiology. However, its implementation as a routine laboratory technique remains challenging due to the growth of data at a faster rate than can be effectively analyzed and critical gaps in bioinformatics knowledge.ResultsTo address both issues, CamPype was developed as a new bioinformatics workflow for the genomics analysis of sequencing data of bacteria, especially Campylobacter, which is the main cause of gastroenteritis worldwide making a negative impact on the economy of the public health systems. CamPype allows fully customization of stages to run and tools to use, including read quality control filtering, read contamination, reads extension and assembly, bacterial typing, genome annotation, searching for antibiotic resistance genes, virulence genes and plasmids, pangenome construction and identification of nucleotide variants. All results are processed and resumed in an interactive HTML report for best data visualization and interpretation.ConclusionsThe minimal user intervention of CamPype makes of this workflow an attractive resource for microbiology laboratories with no expertise in bioinformatics as a first line method for bacterial typing and epidemiological analyses, that would help to reduce the costs of disease outbreaks, or for comparative genomic analyses. CamPype is publicly available at https://github.com/JoseBarbero/CamPype.

In Silico Evaluation of Variant Calling Methods for Bacterial Whole-Genome Sequencing Assays.

Identification and analysis of clinically relevant strains of bacteria increasingly relies on whole-genome sequencing. The downstream bioinformatics steps necessary for calling variants from short-read sequences are well-established but seldom validated against haploid genomes. We devised an in silico workflow to introduce single nucleotide polymorphisms (SNP) and indels into bacterial reference genomes, and computationally generate sequencing reads based on the mutated genomes. We then applied the method to Mycobacterium tuberculosis H37Rv, Staphylococcus aureus NCTC 8325, and Klebsiella pneumoniae HS11286, and used the synthetic reads as truth sets for evaluating several popular variant callers. Insertions proved especially challenging for most variant callers to correctly identify, relative to deletions and single nucleotide polymorphisms. With adequate read depth, however, variant callers that use high quality soft-clipped reads and base mismatches to perform local realignment consistently had the highest precision and recall in identifying insertions and deletions ranging from1 to 50 bp. The remaining variant callers had lower recall values associated with identification of insertions greater than 20 bp.

Machine-Learning Model for Prediction of Cefepime Susceptibility in Escherichia coli from Whole-Genome Sequencing Data.

The declining cost of performing bacterial whole-genome sequencing (WGS) coupled with the availability of large libraries of sequence data for well-characterized isolates have enabled the application of machine-learning (ML) methods to the development of nonlinear sequence-based predictive models. We tested the ML-based model developed by Next Gen Diagnostics for prediction of cefepime phenotypic susceptibility results in Escherichia coli. A cohort of 100 isolates of E. coli recovered from urine (n = 77) and blood (n = 23) cultures were used. The cefepime MIC was determined in triplicate by reference broth microdilution and classified as susceptible (MIC of ≤2 μg/mL) or not susceptible (MIC of ≥4 μg/mL) using the 2022 Clinical and Laboratory Standards Institute breakpoints. Five isolates generated both susceptible and not susceptible MIC results, yielding categorical agreement of 95% for the reference method to itself. Categorical agreement of ML to MIC interpretations was 97%, with 2 very major (false, susceptible) and 1 major (false, not susceptible) errors. One very major error occurred for an isolate with blaCTX-M-27 (MIC mode, ≥32 μg/mL) and one for an isolate with blaTEM-34 for which the MIC cefepime mode was 4 μg/mL. One major error was for an isolate with blaCTX-M-27 but with a MIC mode of 2 μg/mL. These preliminary data demonstrated performance of ML for a clinically important antimicrobial-species pair at a caliber similar to phenotypic methods, encouraging wider development of sequence-based susceptibility prediction and its validation and use in clinical practice.

Benchmarking of two bioinformatic workflows for the analysis of whole-genome sequenced Staphylococcus aureus collected from patients with suspected sepsis

BackgroundThe rapidly growing area of sequencing technologies, and more specifically bacterial whole-genome sequencing, could offer applications in clinical microbiology, including species identification of bacteria, prediction of genetic antibiotic susceptibility and virulence genes simultaneously. To accomplish the aforementioned points, the commercial cloud-based platform, 1928 platform (1928 Diagnostics, Gothenburg, Sweden) was benchmarked against an in-house developed bioinformatic pipeline as well as to reference methods in the clinical laboratory.MethodsWhole-genome sequencing data retrieved from 264 Staphylococcus aureus isolates using the Illumina HiSeq X next-generation sequencing technology was used. The S. aureus isolates were collected during a prospective observational study of community-onset severe sepsis and septic shock in adults at Skaraborg Hospital, in the western region of Sweden. The collected isolates were characterized according to accredited laboratory methods i.e., species identification by MALDI-TOF MS analysis and phenotypic antibiotic susceptibility testing (AST) by following the EUCAST guidelines. Concordance between laboratory methods and bioinformatic tools, as well as concordance between the bioinformatic tools was assessed by calculating the percent of agreement.ResultsThere was an overall high agreement between predicted genotypic AST and phenotypic AST results, 98.0% (989/1006, 95% CI 97.3–99.0). Nevertheless, the 1928 platform delivered predicted genotypic AST results with lower very major error rates but somewhat higher major error rates compared to the in-house pipeline. There were differences in processing times i.e., minutes versus hours, where the 1928 platform delivered the results faster. Furthermore, the bioinformatic workflows showed overall 99.4% (1267/1275, 95% CI 98.7–99.7) agreement in genetic prediction of the virulence gene characteristics and overall 97.9% (231/236, 95% CI 95.0–99.2%) agreement in predicting the sequence types (ST) of the S. aureus isolates.ConclusionsAltogether, the benchmarking disclosed that both bioinformatic workflows are able to deliver results with high accuracy aiding diagnostics of severe infections caused by S. aureus. It also illustrates the need of international agreement on quality control and metrics to facilitate standardization of analytical approaches for whole-genome sequencing based predictions.

A Practical Bioinformatics Workflow for Routine Analysis of Bacterial WGS Data.

The use of whole-genome sequencing (WGS) for bacterial characterisation has increased substantially in the last decade. Its high throughput and decreasing cost have led to significant changes in outbreak investigations and surveillance of a wide variety of microbial pathogens. Despite the innumerable advantages of WGS, several drawbacks concerning data analysis and management, as well as a general lack of standardisation, hinder its integration in routine use. In this work, a bioinformatics workflow for (Illumina) WGS data is presented for bacterial characterisation including genome annotation, species identification, serotype prediction, antimicrobial resistance prediction, virulence-related genes and plasmid replicon detection, core-genome-based or single nucleotide polymorphism (SNP)-based phylogenetic clustering and sequence typing. Workflow was tested using a collection of 22 in-house sequences of Salmonella enterica isolates belonging to a local outbreak, coupled with a collection of 182 Salmonella genomes publicly available. No errors were reported during the execution period, and all genomes were analysed. The bioinformatics workflow can be tailored to other pathogens of interest and is freely available for academic and non-profit use as an uploadable file to the Galaxy platform.

A low-cost genomics workflow enables isolate screening and strain-level analyses within microbiomes

Earth’s environments harbor complex consortia of microbes that affect processes ranging from host health to biogeochemical cycles. Understanding their evolution and function is limited by an inability to isolate genomes in a high-throughput manner. Here, we present a workflow for bacterial whole-genome sequencing using open-source labware and the OpenTrons robotics platform, reducing costs to approximately $10 per genome. We assess genomic diversity within 45 gut bacterial species from wild-living chimpanzees and bonobos. We quantify intraspecific genomic diversity and reveal divergence of homologous plasmids between hosts. This enables population genetic analyses of bacterial strains not currently possible with metagenomic data alone.

Accurate Identification of Klebsiella variicola by MALDI-TOF Mass Spectrometry in Clinical Microbiology Laboratories.

ABSTRACTKlebsiella variicola is a pathogen that is increasingly recognized as being associated with human infections, but the methods available to clinical microbiology laboratories for accurate identification are limited. In this study, we assessed the accuracy of identification of K. variicola by matrix-assisted laser desorption ionization–time of flight (MALDI-TOF) mass spectrometry using genetic identification with multiplex PCR as the reference method. Antimicrobial susceptibilities and virulence of K. variicola strains were also investigated. Fifty-five Klebsiella pneumoniae, 26 K. variicola, and 2 Klebsiella quasipneumoniae clinical strains were used for evaluation. Both MALDI Biotyper with library version 9 and Klebsiella MALDI TypeR, a web-based species identification tool using MALDI-TOF data, accurately identified all K. variicola strains. In addition, two strains of K. quasipneumoniae were accurately identified with Klebsiella MALDI TypeR. Whole-genome sequencing confirmed the accurate identification to the subspecies level by Klebsiella MALDI TypeR for four strains (two strains each of K. variicola subsp. variicola and K. quasipneumoniae subsp. similipneumoniae). While 13 strains, 3 strains, and 1 strain of K. pneumoniae showed nonsusceptibility to ampicillin-sulbactam, ceftriaxone, and meropenem, respectively, all strains of K. variicola were susceptible to all tested antimicrobial agents. Although two K. variicola strains were positive for the string test, no K. variicola strains harbored any of the genes associated with hypervirulence of K. pneumoniae. Accurate identification of the K. pneumoniae complex, including K. variicola, by MALDI-TOF in clinical microbiology laboratories is expected to clarify the clinical characteristics of each species in the future.IMPORTANCE Recent widespread use of bacterial whole-genome sequencing analysis has resulted in the proposal of novel bacterial species and reclassification of taxonomy. Accurate methods for identification of bacterial species in clinical microbiology laboratories are essential to accumulate information on the clinical characteristics of each bacterial species. Klebsiella variicola is a member of the Klebsiella pneumoniae complex, and its association with human infections has been increasingly recognized, but accurate identification methods approved for use in clinical microbiology laboratories have been limited thus far. The findings of the present study suggest that K. variicola can be accurately identified by matrix-assisted laser desorption ionization–time of flight mass spectrometry using updated library or web-based identification tools. Accurate identification will promote exploration of clinical characteristics of K. variicola.

Imaging sensitive and drug-resistant bacterial infection with [11C]-trimethoprim.

BACKGROUNDSeveral molecular imaging strategies can identify bacterial infections in humans. PET affords the potential for sensitive infection detection deep within the body. Among PET-based approaches, antibiotic-based radiotracers, which often target key bacterial-specific enzymes, have considerable promise. One question for antibiotic radiotracers is whether antimicrobial resistance (AMR) reduces specific accumulation within bacteria, diminishing the predictive value of the diagnostic test.METHODSUsing a PET radiotracer based on the antibiotic trimethoprim (TMP), [11C]-TMP, we performed in vitro uptake studies in susceptible and drug-resistant bacterial strains and whole-genome sequencing (WGS) in selected strains to identify TMP resistance mechanisms. Next, we queried the NCBI database of annotated bacterial genomes for WT and resistant dihydrofolate reductase (DHFR) genes. Finally, we initiated a first-in-human protocol of [11C]-TMP in patients infected with both TMP-sensitive and TMP-resistant organisms to demonstrate the clinical feasibility of the tool.RESULTSWe observed robust [11C]-TMP uptake in our panel of TMP-sensitive and -resistant bacteria, noting relatively variable and decreased uptake in a few strains of P. aeruginosa and E. coli. WGS showed that the vast majority of clinically relevant bacteria harbor a WT copy of DHFR, targetable by [11C]-TMP, and that despite the AMR, these strains should be “imageable.” Clinical imaging of patients with [11C]-TMP demonstrated focal radiotracer uptake in areas of infectious lesions.CONCLUSIONThis work highlights an approach to imaging bacterial infection in patients, which could affect our understanding of bacterial pathogenesis as well as our ability to better diagnose infections and monitor response to therapy.TRIAL REGISTRATIONClinicalTrials.gov NCT03424525.FUNDINGInstitute for Translational Medicine and Therapeutics, Burroughs Wellcome Fund, NIH Office of the Director Early Independence Award (DP5-OD26386), and University of Pennsylvania NIH T32 Radiology Research Training Grant (5T32EB004311-12).

Deep clustering of bacterial tree images.

The field of genomic epidemiology is rapidly growing as many jurisdictions begin to deploy whole-genome sequencing (WGS) in their national or regional pathogen surveillance programmes. WGS data offer a rich view of the shared ancestry of a set of taxa, typically visualized with phylogenetic trees illustrating the clusters or subtypes present in a group of taxa, their relatedness and the extent of diversification within and between them. When methicillin-resistant Staphylococcus aureus (MRSA) arose and disseminated widely, phylogenetic trees of MRSA-containing types of S. aureus had a distinctive ‘comet’ shape, with a ‘comet head’ of recently adapted drug-resistant isolates in the context of a ‘comet tail’ that was predominantly drug-sensitive. Placing an S. aureus isolate in the context of such a ‘comet’ helped public health laboratories interpret local data within the broader setting of S. aureus evolution. In this work, we ask what other tree shapes, analogous to the MRSA comet, are present in bacterial WGS datasets. We extract trees from large bacterial genomic datasets, visualize them as images and cluster the images. We find nine major groups of tree images, including the ‘comets’, star-like phylogenies, ‘barbell’ phylogenies and other shapes, and comment on the evolutionary and epidemiological stories these shapes might illustrate.This article is part of a discussion meeting issue ‘Genomic population structures of microbial pathogens’.

Use of whole-genome sequencing to detect transmission of group A Streptococcus in Houston, TX.

We used a combination of local, comprehensive strain surveillance and bacterial whole-genome sequencing to identify potential transmission events of group A streptococcus (GAS) in Houston, TX, USA. We identified pharyngeal and skin and soft tissue sources of infection as having important roles in community GAS transmission, including invasive diseases.