Combining Standard Molecular Typing and Whole Genome Sequencing to Investigate Pseudomonas aeruginosa Epidemiology in Intensive Care Units

Retrospective genome-oriented analysis reveals low transmission rate of multidrug-resistant Pseudomonas aeruginosa from contaminated toilets at a bone marrow transplant unit

BACKGROUND

Prevention of toilet-to-patient transmission of multidrug-resistant Pseudomonas aeruginosa (MDR PA) poses management-related challenges at many bone marrow transplant units (BMTU). Using whole-genome sequencing (WGS), we conducted a longitudinal retrospective analysis of the toilet-to-patient transmission rate for MDR PA under existing infection control (IC) measures at a BMTU with persistent MDR PA toilet colonization.

METHODS

The local IC bundle comprised 1.) patient education regarding IC, 2.) routine patient screening, 3.) toilet flushing volume of 9L, 4.) bromination of toilet water tanks, and 5.) toilet decontamination using hydrogen peroxide. Toilet water was sampled periodically between 2016-2021 (minimum every three months - 26 intervals). Upon MDR PA detection, disinfection and re-sampling were repeated until ≤ 3 cfu/100ml was reached. WGS was performed retrospectively on all available MDR PA isolates (90 of 117 positive environmental samples, 10 of 14 patients - including nine nosocomial).

RESULTS

WGS of patient isolates identified six sequence types (STs), with ST235/CT1352/FIM-1 and ST309/CT3049/no-carbapenemase being predominant (three isolates each). Environmental sampling consistently identified MDR PA ST235 (65.5% ST235/CT1352/FIM-1), showing low genetic diversity (difference of ≤29 alleles by cgMLST). This indicates that direct toilet-to-patient transmission was infrequent although MDR PA was widespread (detection on 79 occasions, detection in every toilet). Only three MDR PA patient isolates can be attributed to the ST235/CT1352/FIM-1 toilet MRD PA population over six years.

CONCLUSIONS

Genome-oriented environmental and patient surveillance suggests that the persistent presence of MDR PA poses a potential risk for acquisition, but with stringent targeted toilet disinfection, only three highly vulnerable patients experienced nosocomial transmission.

Epidemiology of healthcare-associated Pseudomonas aeruginosa in intensive care units: are sink drains to blame?

BACKGROUND

Pseudomonas aeruginosa (PA) is a common cause of healthcare-associated infection (PA-HAI) in the intensive care unit (ICU).

AIM

To describe the epidemiology of PA-HAI in ICUs in Ontario, Canada, and to identify episodes of sink-to-patient PA transmission.

METHODS

This was a prospective cohort study of patients in six ICUs from 2018 to 2019, with retrieval of PA clinical isolates, and PA-screening of antimicrobial-resistant organism surveillance rectal swabs, and of sink drain, air, and faucet samples. All PA isolates underwent whole-genome sequencing. PA-HAI was defined using US National Healthcare Safety Network criteria. ICU-acquired PA was defined as PA isolated from specimens obtained ≥48 h after ICU admission in those with prior negative rectal swabs. Sink-to-patient PA transmission was defined as ICU-acquired PA with close genomic relationship to isolate(s) previously recovered from sinks in a room/bedspace occupied 3-14 days prior to collection date of the relevant patient specimen.

FINDINGS

Over ten months, 72 PA-HAIs occurred among 60/4263 admissions. The rate of PA-HAI was 2.40 per 1000 patient-ICU-days; higher in patients who were PA-colonized on admission. PA-HAI was associated with longer stay (median: 26 vs 3 days uninfected; P < 0.001) and contributed to death in 22/60 cases (36.7%). Fifty-eight admissions with ICU-acquired PA were identified, contributing 35/72 (48.6%) PA-HAIs. Four patients with five PA-HAIs (6.9%) had closely related isolates previously recovered from their room/bedspace sinks.

CONCLUSION

Nearly half of PA causing HAI appeared to be acquired in ICUs, and 7% of PA-HAIs were associated with sink-to-patient transmission. Sinks may be an under-recognized reservoir for HAIs.

Genomic diversity of β-lactamase producing Pseudomonas aeruginosa in Iran; the impact of global high-risk clones

BACKGROUND

Hospital-acquired infections caused by multidrug-resistant Pseudomonas aeruginosa incline hospital stay and costs of treatment that resulted in an increased mortality rate. The frequency of P. aeruginosa high-risk clones producing carbapenemases was investigated in our clinical samples.

METHODS

In this cross-sectional study, 155 non-repetitive P. aeruginosa isolates were included from different medical centers of Iran. Antibiotic susceptibility testing was determined, and the presence of β-lactamases were sought by phenotypic and genotypic methods. The clonal relationship of all isolates was investigated, and multi-locus sequence typing (MLST) was used for finding the sequence types of carbapenemase-producers.

RESULTS

The agent with highest percent susceptibility rate was recorded for colistin (94.9%). MOX and FOX were found both as low as 1.95% (3/155). The most frequent narrow spectrum β-lactamase was SHV with 7.7% (12/155) followed by PER, OXA-1, and TEM with the frequency of 7.1% (11/155), 3.2% (5/155), and 1.3% (2/155), respectively. Carbapenemases were detected in 28 isolates (18%). The most frequent carbapenemase was IMP with 9% (14/155) followed by NDM, 8.4% (13/155). OXA-48 and VIM were also detected both per one isolate (0.65%). MLST of carbapenem resistant P. aeruginosa isolates revealed that ST244, ST664, ST235, and ST357 were spread in subjected clinical settings. REP-PCR uncovered high genomic diversity in our clinical setting.

CONCLUSION

Clonal proliferation of ST235 strain plays a key role in the propagation of MDR pattern in P. aeruginosa. Our data showed that high-risk clones has distributed in Iran, and programs are required to limit spreading of these clones.

High-Throughput Short Sequence Typing Schemes for Pseudomonas aeruginosa and Stenotrophomonas maltophilia Pure Culture and Environmental DNA

Molecular typing techniques are utilized to determine genetic similarities between bacterial isolates. However, the use of environmental DNA profiling to assess epidemiologic links between patients and their environment has not been fully explored. This work reports the development and validation of two high-throughput short sequence typing (HiSST) schemes targeting the opportunistic pathogens Pseudomonas aeruginosa and Stenotrophomonas maltophilia, along with a modified SM2I selective medium for the specific isolation of S. maltophilia. These HiSST schemes are based on four discriminative loci for each species and demonstrate high discriminating power, comparable to pairwise whole-genome comparisons. Each scheme includes species-specific PCR primers for precise differentiation from closely related taxa, without the need for upstream culture-dependent methods. For example, the primers targeting the bvgS locus make it possible to distinguish P. aeruginosa from the very closely related Pseudomonas paraeruginosa sp. nov. The selected loci included in the schemes are adapted to massive parallel amplicon sequencing technology. An R-based script implemented in the DADA2 pipeline was assembled to facilitate HiSST analyses for efficient and accurate genotyping of P. aeruginosa and S. maltophilia. We demonstrate the performance of both schemes through in silico validations, assessments against reference culture collections, and a case study involving environmental samples.

Diversity and evolution of an abundant ICEclc family of integrative and conjugative elements in Pseudomonas aeruginosa

IMPORTANCE

Microbial populations swiftly adapt to changing environments through horizontal gene transfer. While the mechanisms of gene transfer are well known, the impact of environmental conditions on the selection of transferred gene functions remains less clear. We investigated ICEs, specifically the ICEclc-type, in Pseudomonas aeruginosa clinical isolates. Our findings revealed co-evolution between ICEs and their hosts, with ICE transfers occurring within strains. Gene functions carried by ICEs are positively selected, including potential virulence factors and heavy metal resistance. Comparison to publicly available P. aeruginosa genomes unveiled widespread antibiotic-resistance determinants within ICEclc clades. Thus, the ubiquitous ICEclc family significantly contributes to P. aeruginosa's adaptation and fitness in diverse environments.

Engineered live bacteria suppress Pseudomonas aeruginosa infection in mouse lung and dissolve endotracheal-tube biofilms

Engineered live bacteria could provide a new modality for treating lung infections, a major cause of mortality worldwide. In the present study, we engineered a genome-reduced human lung bacterium, Mycoplasma pneumoniae, to treat ventilator-associated pneumonia, a disease with high hospital mortality when associated with Pseudomonas aeruginosa biofilms. After validating the biosafety of an attenuated M. pneumoniae chassis in mice, we introduced four transgenes into the chromosome by transposition to implement bactericidal and biofilm degradation activities. We show that this engineered strain has high efficacy against an acute P. aeruginosa lung infection in a mouse model. In addition, we demonstrated that the engineered strain could dissolve biofilms formed in endotracheal tubes of patients with ventilator-associated pneumonia and be combined with antibiotics targeting the peptidoglycan layer to increase efficacy against Gram-positive and Gram-negative bacteria. We expect our M. pneumoniae-engineered strain to be able to treat biofilm-associated infections in the respiratory tract.

Resistance genomics and molecular epidemiology of high-risk clones of ESBL-producing Pseudomonas aeruginosa in young children

Introduction

The emergence of multidrug-resistant Pseudomonas aeruginosa poses a global threat, but the distribution and resistance profiling are unclear, especially in young children. Infections due to P. aeruginosa are common, associated with high mortality, and increasingly β-lactam drug resistant.

Methods

We studied the molecular epidemiology and antibiotic resistance mechanisms in 294 clinicalisolates of P. aeruginosa from a pediatric hospital in China. Non-duplicate isolates were recovered from clinical cases and were identified using an API-20 kit followed by antimicrobial susceptibility testing using the VITEK®2 compact system (BioMerieux, France) and also by broth dilution method. In addition, a double-disc synergy test for the ESBL/E-test for MBL was performed. The presence of beta-lactamases, plasmid types, and sequence types was determined by PCR and sequencing.

Results

Fifty-six percent (n = 164) of the isolates were resistant to piperacillin-tazobactam, followed by cefepime (40%; n = 117), ceftazidime (39%; n = 115), imipenem (36%; n = 106), meropenem (33%; n = 97), and ciprofloxacin (32%; n = 94). Forty-two percent (n = 126) of the isolates were positive for ESBL according to the double-disc synergy test. The blaCTX-M-15 cephalosporinase was observed in 32% (n = 40/126), while 26% (n = 33/126) werepositive for blaNDM-1 carbapenemase. Aminoglycoside resistance gene aac(3)IIIawas observed in 16% (n = 20/126), and glycylcyclines resistance gene tet(A) was observed in 12% (n = 15/126) of the isolates. A total of 23 sequence types were detected, including ST1963 (12%; n = 16), followed by ST381 (11%; n = 14), ST234 (10%; n = 13), ST145 (58%; n = 10), ST304 (57%; n = 9), ST663 (5%; n = 7), and a novel strain. In ESBL-producing P. aeruginosa, 12 different Incompatibility groups (Inc) were observed, the most common being IncFI, IncFIS, and IncA/C. The MOBP was the most common plasmid type, followed by MOBH, MOBF, and MOBQ.

Discussion

Our data suggest that the spread of antibiotic resistance is likely due toclonal spread and dissemination of different clinical strains of P. aeruginosa harbouring different plasmids. This is a growing threat in hospitals particularly in young children which needs robust prevention strategies.

Molecular epidemiology of antimicrobial-resistant Pseudomonas aeruginosa in a veterinary teaching hospital environment

This study aimed to investigate sites for colonization and molecular epidemiology of antimicrobial-resistant Pseudomonas aeruginosa in a veterinary teaching hospital. Bacterial specimens from surface and liquid samples (n = 165) located in five rooms were collected three times every 2 months, and antimicrobial susceptibility was subsequently determined by minimum inhibitory concentrations. The genomes of resistant strains were further analyzed using whole-genome sequencing. Among 19 P. aeruginosa isolates (11.5%, 19/165), sinks were the most frequent colonization site (53.3%), followed by rubber tubes (44.4%), and anesthesia-breathing circuit (33.3%). The highest resistance to gentamicin (47.4%), followed by piperacillin/tazobactam (36.8%), levofloxacin (36.8%), and ciprofloxacin (36.8%), was observed from 19 P. aeruginosa isolates, of which 10 were resistant strains. Of these 10 antimicrobial-resistant isolates, five were multidrug-resistant isolates, including carbapenem. From the multilocus sequence typing (MLST) analysis, five sequence types (STs), including a high-risk clone of human ST235 (n = 3), and ST244 (n = 3), ST606 (n = 2), ST485 (n = 1), and ST3405 (n = 1) were identified in resistant strains. Multiresistant genes were identified consistent with STs, except ST235. The MLST approach and single nucleotide polymorphism analysis revealed a link between resistant strains from ward rooms and those from examination, wound care, and operating rooms. The improvement of routine cleaning, especially of sink environments, and the continued monitoring of antimicrobial resistance of P. aeruginosa in veterinary hospitals are necessary to prevent the spread of resistant clones and ensure infection control.

Parallel Evolution of Pseudomonas aeruginosa during a Prolonged ICU-Infection Outbreak

Most knowledge about Pseudomonas aeruginosa pathoadaptation is derived from studies on airway colonization in cystic fibrosis; little is known about adaptation in acute settings. P. aeruginosa frequently affects burned patients and the burn wound niche has distinct properties that likely influence pathoadaptation. This study aimed to genetically and phenotypically characterize P. aeruginosa isolates collected during an outbreak of infection in a burn intensive care unit (ICU). Sequencing reads from 58 isolates of ST1076 P. aeruginosa taken from 23 patients were independently mapped to a complete reference genome for the lineage (H25338); genetic differences were identified and were used to define the population structure. Comparative genomic analysis at single-nucleotide resolution identified pathoadaptive genes that evolved multiple, independent mutations. Three key phenotypic assays (growth performance, motility, carbapenem resistance) were performed to complement the genetic analysis for 47 unique isolates. Population structure for the ST1076 lineage revealed 11 evolutionary sublineages. Fifteen pathoadaptive genes evolved mutations in at least two sublineages. The most prominent functional classes affected were transcription/two-component regulatory systems, and chemotaxis/motility and attachment. The most frequently mutated gene was oprD, which codes for outer membrane porin involved in uptake of carbapenems. Reduced growth performance and motility were found to be adaptive phenotypic traits, as was high level of carbapenem resistance, which correlated with higher carbapenem consumption during the outbreak. Multiple prominent linages evolved each of the three traits in parallel providing evidence that they afford a fitness advantage for P. aeruginosa in the context of human burn infection. IMPORTANCE Pseudomonas aeruginosa is a Gram-negative pathogen causing infections in acutely burned patients. The precise mechanisms required for the establishment of infection in the burn setting, and adaptive traits underpinning prolonged outbreaks are not known. We have assessed genotypic data from 58 independent P. aeruginosa isolates taken from a single lineage that was responsible for an outbreak of infection in a burn ICU that lasted for almost 2.5 years and affected 23 patients. We identified a core set of 15 genes that we predict to control pathoadaptive traits in the burn infection based on the frequency with which independent mutations evolved. We combined the genotypic data with phenotypic data (growth performance, motility, antibiotic resistance) and clinical data (antibiotic consumption) to identify adaptive phenotypes that emerged in parallel. High-level carbapenem resistance evolved rapidly, and frequently, in response to high clinical demand for this antibiotic class during the outbreak.

A High-Throughput Short Sequence Typing Scheme for Serratia marcescens Pure Culture and Environmental DNA

Molecular typing methods are used to characterize the relatedness between bacterial isolates involved in infections. These approaches rely mostly on discrete loci or whole-genome sequencing (WGS) analyses of pure cultures. On the other hand, their application to environmental DNA profiling to evaluate epidemiological relatedness among patients and environments has received less attention. We developed a specific, high-throughput short sequence typing (HiSST) method for the opportunistic human pathogen Serratia marcescens. Genes displaying the highest polymorphism were retrieved from the core genome of 60 S. marcescens strains. Bioinformatics analyses showed that use of only three loci (within bssA, gabR, and dhaM) distinguished strains with a high level of efficiency. This HiSST scheme was applied to an epidemiological survey of S. marcescens in a neonatal intensive care unit (NICU). In a first case study, a strain responsible for an outbreak in the NICU was found in a sink drain of this unit, by using HiSST scheme and confirmed by WGS. The HiSST scheme was also applied to environmental DNA extracted from sink-environment samples. Diversity of S. marcescens was modest, with 11, 6, and 4 different sequence types (ST) of gabR, bssA, and dhaM loci among 19 sink drains, respectively. Epidemiological relationships among sinks were inferred on the basis of pairwise comparisons of ST profiles. Further research aimed at relating ST distribution patterns to environmental features encompassing sink location, utilization, and microbial diversity is needed to improve the surveillance and management of opportunistic pathogens. IMPORTANCE Serratia marcescens is an important opportunistic human pathogen, often multidrug resistant and involved in outbreaks of nosocomial infections in neonatal intensive care units. Here, we propose a quick and user-friendly method to select the best typing scheme for nosocomial outbreaks in relating environmental and clinical sources. This method, named high-throughput short sequence typing (HiSST), allows to distinguish strains and to explore the diversity profile of nonculturable S. marcescens. The application of HiSST profile analysis for environmental DNA offers new possibilities to track opportunistic pathogens, identify their origin, and relate their distribution pattern with environmental features encompassing sink location, utilization, and microbial diversity. Adaptation of the method to other opportunistic pathogens is expected to improve knowledge regarding their ecology, which is of significant interest for epidemiological risk assessment and elaborate outbreak mitigation strategies.

Outbreak of Pseudomonas aeruginosa producing VIM carbapenemase in an intensive care unit and its termination by implementation of waterless patient care

BACKGROUND

Long-term outbreaks of multidrug-resistant Gram-negative bacilli related to hospital-building water systems have been described. However, successful mitigation strategies have rarely been reported. In particular, environmental disinfection or replacement of contaminated equipment usually failed to eradicate environmental sources of Pseudomonas aeruginosa.

METHODS

We report the investigation and termination of an outbreak of P. aeruginosa producing VIM carbapenemase (PA-VIM) in the adult intensive care unit (ICU) of a Swiss tertiary care hospital with active case finding, environmental sampling and whole genome sequencing (WGS) of patient and environmental strains. We also describe the implemented control strategies and their effectiveness on eradication of the environmental reservoir.

RESULTS

Between April 2018 and September 2020, 21 patients became either infected or colonized with a PA-VIM strain. For 16 of them, an acquisition in the ICU was suspected. Among 131 environmental samples collected in the ICU, 13 grew PA-VIM in sink traps and drains. WGS confirmed the epidemiological link between clinical and environmental strains and the monoclonal pattern of the outbreak. After removing sinks from patient rooms and implementation of waterless patient care, no new acquisition was detected in the ICU within 8 months after the intervention.

DISCUSSION

Implementation of waterless patient care with removal of the sinks in patient rooms was successful for termination of a PA-VIM ICU outbreak linked to multiple environmental water sources. WGS provides highly discriminatory accuracy to investigate environment-related outbreaks.

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.

It's Not Easy Being Green: A Narrative Review on the Microbiology, Virulence and Therapeutic Prospects of Multidrug-Resistant Pseudomonas aeruginosa

Pseudomonas aeruginosa is the most frequent cause of infection among non-fermenting Gram-negative bacteria, predominantly affecting immunocompromised patients, but its pathogenic role should not be disregarded in immunocompetent patients. These pathogens present a concerning therapeutic challenge to clinicians, both in community and in hospital settings, due to their increasing prevalence of resistance, and this may lead to prolonged therapy, sequelae, and excess mortality in the affected patient population. The resistance mechanisms of P. aeruginosa may be classified into intrinsic and acquired resistance mechanisms. These mechanisms lead to occurrence of resistant strains against important antibiotics-relevant in the treatment of P. aeruginosa infections-such as β-lactams, quinolones, aminoglycosides, and colistin. The occurrence of a specific resistotype of P. aeruginosa, namely the emergence of carbapenem-resistant but cephalosporin-susceptible (Car-R/Ceph-S) strains, has received substantial attention from clinical microbiologists and infection control specialists; nevertheless, the available literature on this topic is still scarce. The aim of this present review paper is to provide a concise summary on the adaptability, virulence, and antibiotic resistance of P. aeruginosa to a readership of basic scientists and clinicians.

Novel Molecular Markers Linked to Pseudomonas aeruginosa Epidemic High-Risk Clones

The population structure of Pseudomonas aeruginosa is panmictic-epidemic in nature, with the prevalence of some high-risk clones. These clones are often linked to virulence, antibiotic resistance, and more morbidity. The clonal success of these lineages has been linked to acquisition and spread of mobile genetic elements. The main aim of the study was to explore other molecular markers that explain their global success. A comprehensive set of 528 completely sequenced P. aeruginosa genomes was analyzed. The population structure was examined using Multilocus Sequence Typing (MLST). Strain relationships analysis and diversity analysis were performed using the geoBURST Full Minimum Spanning Tree (MST) algorithm and hierarchical clustering. A phylogenetic tree was constructed using the Unweighted Pair Group Method with Arithmetic mean (UPGMA) algorithm. A panel of previously investigated resistance markers were examined for their link to high-risk clones. A novel panel of molecular markers has been identified in relation to risky clones including armR, ampR, nalC, nalD, mexZ, mexS, gyrAT83I, gyrAD87N, nalCE153Q, nalCS46A, parCS87W, parCS87L, ampRG283E, ampRM288R, pmrALeu71Arg, pmrBGly423Cys, nuoGA890T, pstBE89Q, phoQY85F, arnAA170T, arnDG206C, and gidBE186A. In addition to mobile genetic elements, chromosomal variants in membrane proteins and efflux pump regulators can play an important role in the success of high-risk clones. Finding risk-associated markers during molecular surveillance necessitates applying more infection-control precautions.

NGS-Based S. aureus Typing and Outbreak Analysis in Clinical Microbiology Laboratories: Lessons Learned From a Swiss-Wide Proficiency Test

Whole genome sequencing (WGS) enables high resolution typing of bacteria up to the single nucleotide polymorphism (SNP) level. WGS is used in clinical microbiology laboratories for infection control, molecular surveillance and outbreak analyses. Given the large palette of WGS reagents and bioinformatics tools, the Swiss clinical bacteriology community decided to conduct a ring trial (RT) to foster harmonization of NGS-based bacterial typing. The RT aimed at assessing methicillin-susceptible Staphylococcus aureus strain relatedness from WGS and epidemiological data. The RT was designed to disentangle the variability arising from differences in sample preparation, SNP calling and phylogenetic methods. Nine laboratories participated. The resulting phylogenetic tree and cluster identification were highly reproducible across the laboratories. Cluster interpretation was, however, more laboratory dependent, suggesting that an increased sharing of expertise across laboratories would contribute to further harmonization of practices. More detailed bioinformatic analyses unveiled that while similar clusters were found across laboratories, these were actually based on different sets of SNPs, differentially retained after sample preparation and SNP calling procedures. Despite this, the observed number of SNP differences between pairs of strains, an important criterion to determine strain relatedness given epidemiological information, was similar across pipelines for closely related strains when restricting SNP calls to a common core genome defined by S. aureus cgMLST schema. The lessons learned from this pilot study will serve the implementation of larger-scale RT, as a mean to have regular external quality assessments for laboratories performing WGS analyses in a clinical setting.

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.