Development and Application of a Core Genome Multilocus Sequence Typing Scheme for the Health Care-Associated Pathogen Pseudomonas aeruginosa

A comparison of bacterial colonization between nasogastric and orogastric enteral feeding tubes in infants in the neonatal intensive care unit

OBJECTIVE

Feeding tubes harbor microbial contaminants; studies to date have not explored differences between orogastric (OG) and nasogastric (NG) tube biofilms. We sought to extend a previous analysis by comparing bacterial colonization by location (OG v NG) and by evaluating clinical factors that may affect tube bacterial populations.

STUDY DESIGN

The pharyngeal segments of 41 infant feeding tubes (14 OG and 27 NG) from 41 infants were analyzed by next generation 16 S rRNA sequencing on the MiSeq platform.

RESULTS

At the phylum level, Proteobacteria had the highest relative abundance of both OG and NG tubes. At the genus/species level, nine taxa differed significantly between OG and NG tubes. Alpha and beta diversity analyses showed significant differences between OG and NG tubes with relatively little contribution from clinical factors.

CONCLUSION

The route of feeding tube insertion (oral vs nasal) had a greater impact on bacterial colonization than the assessed clinical factors.

Whole-Genome Sequencing Reveals Diversity of Carbapenem-Resistant Pseudomonas aeruginosa Collected through CDC's Emerging Infections Program, United States, 2016-2018

The CDC's Emerging Infections Program (EIP) conducted population- and laboratory-based surveillance of US carbapenem-resistant Pseudomonas aeruginosa (CRPA) from 2016 through 2018. To characterize the pathotype, 1,019 isolates collected through this project underwent antimicrobial susceptibility testing and whole-genome sequencing. Sequenced genomes were classified using the seven-gene multilocus sequence typing (MLST) scheme and a core genome (cg)MLST scheme was used to determine phylogeny. Both chromosomal and horizontally transmitted mechanisms of carbapenem resistance were assessed. There were 336 sequence types (STs) among the 1,019 sequenced genomes, and the genomes varied by an average of 84.7% of the cgMLST alleles used. Mutations associated with dysfunction of the porin OprD were found in 888 (87.1%) of the genomes and were correlated with carbapenem resistance, and a machine learning model incorporating hundreds of genetic variations among the chromosomal mechanisms of resistance was able to classify resistant genomes. While only 7 (0.1%) isolates harbored carbapenemase genes, 66 (6.5%) had acquired non-carbapenemase β-lactamase genes, and these were more likely to have OprD dysfunction and be resistant to all carbapenems tested. The genetic diversity demonstrates that the pathotype includes a variety of strains, and clones previously identified as high-risk make up only a minority of CRPA strains in the United States. The increased carbapenem resistance in isolates with acquired non-carbapenemase β-lactamase genes suggests that horizontally transmitted mechanisms aside from carbapenemases themselves may be important drivers of the spread of carbapenem resistance in P. aeruginosa.

Whole Genome Multi-Locus Sequence Typing and Genomic Single Nucleotide Polymorphism Analysis for Epidemiological Typing of Pseudomonas aeruginosa From Indonesian Intensive Care Units

We have previously studied carbapenem non-susceptible Pseudomonas aeruginosa (CNPA) strains from intensive care units (ICUs) in a referral hospital in Jakarta, Indonesia (Pelegrin et al., 2019). We documented that CNPA transmissions and acquisitions among patients were variable over time and that these were not significantly reduced by a set of infection control measures. Three high risk international CNPA clones (sequence type (ST)235, ST823, ST357) dominated, and carbapenem resistance was due to carbapenemase-encoding genes and mutations in the porin OprD. Pelegrin et al. (2019) reported core genome analysis of these strains. We present a more refined and detailed whole genome-based analysis of major clones represented in the same dataset. As per our knowledge, this is the first study reporting Single Nucleotide Polymorphisms (wgSNP) analysis of Pseudomonas strains. With whole genome-based Multi Locus Sequence Typing (wgMLST) of the 3 CNPA clones (ST235, ST357 and ST823), three to eleven subgroups with up to 200 allelic variants were observed for each of the CNPA clones. Furthermore, we analyzed these CNPA clone clusters for the presence of wgSNP to redefine CNPA transmission events during hospitalization. A maximum number 35350 SNPs (including non-informative wgSNPs) and 398 SNPs (ST-specific_informative-wgSNPs) were found in ST235, 34,570 SNPs (including non-informative wgSNPs) and 111 SNPs (ST-specific_informative-wgSNPs) in ST357 and 26,443 SNPs (including non-informative SNPs) and 61 SNPs (ST-specific_informative-wgSNPs) in ST823. ST-specific_Informative-wgSNPs were commonly noticed in sensor-response regulator genes. However, the majority of non-informative wgSNPs was found in conserved hypothetical proteins or in uncharacterized proteins. Of note, antibiotic resistance and virulence genes segregated according to the wgSNP analyses. A total of 8 transmission chains for ST235 strains followed by 9 and 4 possible transmission chains for ST357 and ST823 were traceable on the basis of pairwise distances of informative-wgSNPs (0 to 4 SNPs) among the strains. The present study demonstrates the value of detailed whole genome sequence analysis for highly refined epidemiological analysis of P. aeruginosa.

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.

Extensively Drug-Resistant Carbapenemase-Producing Pseudomonas aeruginosa and Medical Tourism from the United States to Mexico, 2018-2019

Carbapenem-resistant Pseudomonas aeruginosa (CRPA) producing the Verona integron‒encoded metallo-β-lactamase (VIM) are highly antimicrobial drug-resistant pathogens that are uncommon in the United States. We investigated the source of VIM-CRPA among US medical tourists who underwent bariatric surgery in Tijuana, Mexico. Cases were defined as isolation of VIM-CRPA or CRPA from a patient who had an elective invasive medical procedure in Mexico during January 2018‒December 2019 and within 45 days before specimen collection. Whole-genome sequencing of isolates was performed. Thirty-eight case-patients were identified in 18 states; 31 were operated on by surgeon 1, most frequently at facility A (27/31 patients). Whole-genome sequencing identified isolates linked to surgeon 1 were closely related and distinct from isolates linked to other surgeons in Tijuana. Facility A closed in March 2019. US patients and providers should acknowledge the risk for colonization or infection after medical tourism with highly drug-resistant pathogens uncommon in the United States.

OUP accepted manuscript

Objectives

To describe the propensity of carbapenem-resistant Pseudomonas aeruginosa to spread within a hospital critical care setting.

Methods

The study was conducted in a 700-bed tertiary centre in Cologne, Germany. P. aeruginosa resistant to piperacillin, ceftazidime, cefepime, imipenem, meropenem and ciprofloxacin, isolated from clinical and screening specimens from four critical care units from 2015 to 2020 were analysed. Genotyping was carried out by WGS (Illumina and MinION). MLST, core genome MLST (cgMLST) and resistome analysis was performed and merged with epidemiological data.

Results

Fifty-five out of 79 non-duplicate P. aeruginosa isolates were available, of which 20 were carbapenemase producers as follows: bla_VIM-1 (n = 1), bla_VIM-2 (n = 17), bla_VIM-4 (n = 1), and bla_NDM-1/bla_GES-5 (n = 1). Forty-two of 55 isolates were hospital-acquired. cgMLST revealed three clusters: Cluster 1 (n = 15, ST111, bla_VIM-2, recovered between 2015 and 2020); Cluster 2 (n = 4, ST970, carbapenemase negative); and Cluster 3 (n = 2, ST357, carbapenemase negative). The bla_VIM-2 gene of Cluster 1 was integrated on the chromosome in a class 1 integron (type In59). Using conventional epidemiology, we were only able to confirm two patient-to-patient transmissions and one room-to-patient transmission on three different ICUs within Cluster 1. Isolates from Cluster 2 represented an outbreak occurring in 2019.

Conclusions

These data give insight into the epidemiology of carbapenem-resistant P. aeruginosa. Transmission dynamics differed between carbapenemase- and non-carbapenemase-producing isolates. A continuous acquisition of clonally related ST111 VIM-2 P. aeruginosa, being the main carbapenemase-producing strain, was observed over the whole study period, as well as an overall higher genomic diversity among non-carbapenemase-producing P. aeruginosa.

A panel of diverse Pseudomonas aeruginosa clinical isolates for research and development

Objectives

Pseudomonas aeruginosa is a leading cause of community- and hospital-acquired infections. Successful treatment is hampered by its remarkable ability to rapidly develop resistance to antimicrobial agents, primarily through mutation. In response, WHO listed carbapenem-resistant P. aeruginosa as a Priority 1 (Critical) pathogen for research and development of new treatments. A key resource in developing effective countermeasures is access to diverse and clinically relevant strains for testing. Herein we describe a panel of 100 diverse P. aeruginosa strains to support this endeavour.

Methods

WGS was performed on 3785 P. aeruginosa isolates in our repository. Isolates were cultured from clinical samples collected from healthcare facilities around the world between 2003 and 2017. Core-genome MLST and high-resolution SNP-based phylogenetic analyses were used to select a panel of 100 strains that captured the genetic diversity of this collection. Antibiotic susceptibility testing was also performed using 14 clinically relevant antibiotics.

Results

This 100-strain diversity panel contained representative strains from 91 different STs, including genetically distinct strains from major epidemic clones ST-111, ST-235, ST-244 and ST-253. Seventy-one distinct antibiotic susceptibility profiles were identified ranging from pan-susceptible to pan-resistant. Known resistance alleles as well as the most prevalent mutations underlying the antibiotic susceptibilities were characterized for all isolates.

Conclusions

This panel provides a diverse and comprehensive set of P. aeruginosa strains for use in developing solutions to antibiotic resistance. The isolates and available metadata, including genome sequences, are available to industry, academia, federal and other laboratories at no additional cost.

National surveillance pilot study unveils a multicenter, clonal outbreak of VIM-2-producing Pseudomonas aeruginosa ST111 in the Netherlands between 2015 and 2017

Verona Integron-encoded Metallo-beta-lactamase (VIM) is the most frequently-encountered carbapenemase in the healthcare-related pathogen Pseudomonas aeruginosa. In the Netherlands, a low-endemic country for antibiotic-resistant bacteria, no national surveillance data on the prevalence of carbapenemase-producing P. aeruginosa (CPPA) was available. Therefore, in 2016, a national surveillance pilot study was initiated to investigate the occurrence, molecular epidemiology, genetic characterization, and resistomes of CPPA among P. aeruginosa isolates submitted by medical microbiology laboratories (MMLs) throughout the country. From 1221 isolates included in the study, 124 (10%) produced carbapenemase (CIM-positive); of these, the majority (95, 77%) were positive for the bla_VIM gene using PCR. Sequencing was performed on 112 CIM-positive and 56 CIM-negative isolates (n = 168), and genetic clustering revealed that 75/168 (45%) isolates were highly similar. This genetic cluster, designated Group 1, comprised isolates that belonged to high-risk sequence type ST111/serotype O12, had similar resistomes, and all but two carried the bla_VIM-2 allele on an identical class 1 integron. Additionally, Group 1 isolates originated from around the country (i.e. seven provinces) and from multiple MMLs. In conclusion, the Netherlands had experienced a nationwide, inter-institutional, clonal outbreak of VIM-2-producing P. aeruginosa for at least three years, which this pilot study was crucial in identifying. A structured, national surveillance program is strongly advised to monitor the spread of Group 1 CPPA, to identify emerging clones/carbapenemase genes, and to detect transmission in and especially between hospitals in order to control current and future outbreaks.

Techniques in bacterial strain typing: past, present, and future

PURPOSE OF REVIEW

The advancement of molecular techniques such as whole-genome sequencing (WGS) has revolutionized the field of bacterial strain typing, with important implications for epidemiological surveillance and outbreak investigations. This review summarizes state-of-the-art techniques in strain typing and examines barriers faced by clinical and public health laboratories in implementing these new methodologies.

RECENT FINDINGS

WGS-based methodologies are on track to become the new 'gold standards' in bacterial strain typing, replacing traditional methods like pulsed-field gel electrophoresis and multilocus sequence typing. These new techniques have an improved ability to identify genetic relationships among organisms of interest. Further, advances in long-read sequencing approaches will likely provide a highly discriminatory tool to perform pangenome analyses and characterize relevant accessory genome elements, including mobile genetic elements carrying antibiotic resistance determinants in real time. Barriers to widespread integration of these approaches include a lack of standardized workflows and technical training.

SUMMARY

Genomic bacterial strain typing has facilitated a paradigm shift in clinical and molecular epidemiology. The increased resolution that these new techniques provide, along with epidemiological data, will facilitate the rapid identification of transmission routes with high confidence, leading to timely and effective deployment of infection control and public health interventions in outbreak settings.

Detection of CTX-M-27 β-Lactamase Genes on Two Distinct Plasmid Types in ST38 Escherichia coli from Three U.S. States

Infections caused by extended-spectrum-β-lactamase (ESBL)-producing Escherichia coli are a significant cause of morbidity and health care costs. Globally, the prevailing clonal type is ST131 in association with the bla_CTX-M-15 β-lactamase gene. However, other ESBLs, such as bla_CTX-M-14 and bla_CTX-M-27, can also be prevalent in some regions. We identified ST38 ESBL-producing E. coli from different regions in the United States which carry bla_CTX-M-27 embedded on two distinct plasmid types, suggesting the potential emergence of new ESBL lineages.

Framing Bacterial Genomics for Public Health(care)

Advancements in comparative genomics have generated significant interest in defining applications for healthcare-associated pathogens. Clinical microbiology, however, relies on increasingly automated platforms to quickly identify pathogens, resistance mechanisms, and therapy options within CLIA- and FDA-approved frameworks. Additionally, and most notably, healthcare-associated pathogens, especially those that are resistant to antibiotics, represent a diverse spectrum of genera harboring complex genetic targets including antibiotic, biocide, and virulence determinants that can be highly transmissible and, at least for antibiotic resistance, serve as potential targets for containment efforts. U.S. public health investments have focused on rapidly detecting outbreaks and emerging resistance in healthcare-associated pathogens using reference, culture-based, and molecular methods that are distributed, for example, across national laboratory network infrastructures. Herein we describe the public health applications of genomic science that are built from the top-down for broad surveillance, as well as the bottom-up, starting with identification of infections and infectious clusters. For healthcare-associated, including antimicrobial-resistant, pathogens, we propose a combination of top-down and bottom-up genomic approaches leveraged across the public health spectrum, from local infection control, to regional and national containment efforts, to national surveillance for understanding emerging strain ecology and fitness of healthcare pathogens.

Establishment and Evaluation of a Core Genome Multilocus Sequence Typing Scheme for Whole-Genome Sequence-Based Typing of Pseudomonas aeruginosa

The environmental bacterium Pseudomonas aeruginosa, particularly multidrug-resistant clones, is often associated with nosocomial infections and outbreaks. Today, core genome multilocus sequence typing (cgMLST) is frequently applied to delineate sporadic cases from nosocomial transmissions. However, until recently, no cgMLST scheme for a standardized typing of P. aeruginosa was available. To establish a novel cgMLST scheme for P. aeruginosa, we initially determined the breadth of the P. aeruginosa population based on MLST data with a Bayesian approach (BAPS). Using genomic data of representative isolates for the whole population and all 12 serogroups, we extracted target genes and further refined them using a random data set of 1,000 P. aeruginosa genomes. Subsequently, we investigated reproducibility and discriminatory ability with repeatedly sequenced isolates and isolates from well-defined outbreak scenarios, respectively, and compared clustering applying two recently published cgMLST schemes. BAPS generated seven P. aeruginosa groups. To cover these and all serogroups, 15 reference strains were used to determine genes common in all strains. After refinement with the data set of 1,000 genomes, the cgMLST scheme consisted of 3,867 target genes, which are representative of the P. aeruginosa population and highly reproducible using biological replicates. We finally evaluated the scheme by reanalyzing two published outbreaks where the authors used single-nucleotide polymorphism (SNP) typing. In both cases, cgMLST was concordant with the previous SNP results and the results of the two other cgMLST schemes. In conclusion, the highly reproducible novel P. aeruginosa cgMLST scheme facilitates outbreak investigations due to the publicly available cgMLST nomenclature.

Molecular epidemiology in current times

Motivated to find options for prevention or intervention, molecular epidemiology aims to identify the host and microbial factors that determine the transmission, manifestation and progression of infectious disease. The genotyping of cultivatable bacterial strains is performed by either anonymous fingerprinting techniques or sequence-based exploration of variable genomic sites. Multilocus sequence typing of housekeeping genes and allele profiling of the core genome have become standard techniques of bacterial strain typing that may be supplemented by whole genome sequencing to explore all single nucleotide variants and/or the composition of the accessory genome. Next, novel protocols to investigate host and microbiome based upon smart third generation sequencing technologies are being developed for an effective surveillance, rapid diagnosis and real-time tracking of infectious diseases.

Comparison of Whole Genome (wg-) and Core Genome (cg-) MLST (BioNumericsTM) Versus SNP Variant Calling for Epidemiological Investigation of Pseudomonas aeruginosa

Whole genome sequencing (WGS) is increasingly used for epidemiological investigations of pathogens. While SNP variant calling is currently considered as the most suitable method, the choice of a representative reference genome and the isolate dependency of results limit standardization and affect resolution in an unknown manner. Whole or core genome Multi Locus Sequence Typing (wg-, cg-MLST) represents an attractive alternative. Here, we assess the accuracy of wg- and cg-MLST by comparing results of four Pseudomonas aeruginosa datasets for which epidemiological and genomic data were previously described. Three datasets included 155 isolates from three different sequence types (ST) of P. aeruginosa collected in our ICUs over a 5-year period. The fourth dataset consisted of 10 isolates from an investigation of P. aeruginosa contaminated hand soap. All isolates were previously analyzed by a core SNP approach. In this study, wg- and cg-MLST were performed in BioNumerics^TM using a scheme developed by Applied-Maths. Correlation between SNP calling and wg- or cg-MLST results were evaluated by calculating linear regressions and their coefficient of correlations (R²) between the number of SNPs and the number of allele differences in pairwise comparison of isolates. The number of SNPs and allele difference between isolates with close epidemiological linkage varies between 0–26 and 0–13, respectively. When compared to core-SNP calling, a higher coefficient of correlation was obtained with cgMLST (R² of 0.92–0.99) than with wgMLST (0.78–0.99). In one dataset, a putative homologous recombination of a large DNA fragment (202 loci) was identified among these isolates, affecting its phylogeny, but with no impact on the epidemiological analysis of outbreak isolates. In conclusion, we showed that the P. aeruginosa wgMLST scheme in BioNumerics^TM is as discriminatory as the core-SNP calling approach and apparently useful for outbreak investigations. We also showed that epidemiological linked isolates showed less than 26 SNPs or 13 allele differences. These are important figures for the distinction between outbreak and non-outbreak isolates when interpreting WGS results. However, as P. aeruginosa is highly recombinant, a cgMLST approach is preferable and caution should be addressed to possible recombination of large DNA fragments.