Accessory genome of the multi-drug resistant ocular isolate of Pseudomonas aeruginosa PA34

Model-driven characterization of functional diversity of Pseudomonas aeruginosa clinical isolates with broadly representative phenotypes

Pseudomonas aeruginosa is a leading cause of infections in immunocompromised individuals and in healthcare settings. This study aims to understand the relationships between phenotypic diversity and the functional metabolic landscape of P. aeruginosa clinical isolates. To better understand the metabolic repertoire of P. aeruginosa in infection, we deeply profiled a representative set from a library of 971 clinical P. aeruginosa isolates with corresponding patient metadata and bacterial phenotypes. The genotypic clustering based on whole-genome sequencing of the isolates, multilocus sequence types, and the phenotypic clustering generated from a multi-parametric analysis were compared to each other to assess the genotype-phenotype correlation. Genome-scale metabolic network reconstructions were developed for each isolate through amendments to an existing PA14 network reconstruction. These network reconstructions show diverse metabolic functionalities and enhance the collective P. aeruginosa pangenome metabolic repertoire. Characterizing this rich set of clinical P. aeruginosa isolates allows for a deeper understanding of the genotypic and metabolic diversity of the pathogen in a clinical setting and lays a foundation for further investigation of the metabolic landscape of this pathogen and host-associated metabolic differences during infection.

Pseudomonas aeruginosa and related antibiotic resistance genes as indicators for wastewater treatment

This review aims to examine the existence of Pseudomonas aeruginosa (P. aeruginosa) and their antibiotic resistance genes (ARGs) in aquatic settings and the alternative treatment ways. P. aeruginosa in a various aquatic environment have been identified as contaminants with impacts on human health and the environment. P. aeruginosa resistance to multiple antibiotics, such as sulfamethoxazole, ciprofloxacin, quinolone, trimethoprim, tetracycline, vancomycin, as well as specific antibiotic resistance genes including sul1, qnrs, blaVIM, blaTEM, blaCTX, blaAIM-1, tetA, ampC, blaVIM. The development of resistance can occur naturally, through mutations, or via horizontal gene transfer facilitated by sterilizing agents. In addition, an overview of the current knowledge on inactivation of Pseudomonas aeruginosa and ARG and the mechanisms of action of various disinfection processes in water and wastewater (UV chlorine processes, catalytic oxidation, Fenton reaction, and ozonation) is given. An overview of the effects of nanotechnology and the resulting wetlands is also given.

Prioritization of Critical Factors for Surveillance of the Dissemination of Antibiotic Resistance in Pseudomonas aeruginosa: A Systematic Review

Pseudomonas aeruginosa is the primary opportunistic human pathogen responsible for a range of acute and chronic infections; it poses a significant threat to immunocompromised patients and is the leading cause of morbidity and mortality for nosocomial infections. Its high resistance to a diverse array of antimicrobial agents presents an urgent health concern. Among the mechanisms contributing to resistance in P. aeruginosa, the horizontal acquisition of antibiotic resistance genes (ARGs) via mobile genetic elements (MGEs) has gained recognition as a substantial concern in clinical settings, thus indicating that a comprehensive understanding of ARG dissemination within the species is strongly required for surveillance. Here, two approaches, including a systematic literature analysis and a genome database survey, were employed to gain insights into ARG dissemination. The genome database enabled scrutinizing of all the available sequence information and various attributes of P. aeruginosa isolates, thus providing an extensive understanding of ARG dissemination within the species. By integrating both approaches, with a primary focus on the genome database survey, mobile ARGs that were linked or correlated with MGEs, important sequence types (STs) carrying diverse ARGs, and MGEs responsible for ARG dissemination were identified as critical factors requiring strict surveillance. Although human isolates play a primary role in dissemination, the importance of animal and environmental isolates has also been suggested. In this study, 25 critical mobile ARGs, 45 critical STs, and associated MGEs involved in ARG dissemination within the species, are suggested as critical factors. Surveillance and management of these prioritized factors across the One Health sectors are essential to mitigate the emergence of multidrug-resistant (MDR) and extensively resistant (XDR) P. aeruginosa in clinical settings.

Model-driven characterization of functional diversity of Pseudomonas aeruginosa clinical isolates with broadly representative phenotypes

Pseudomonas aeruginosa is a leading cause of infections in immunocompromised individuals and in healthcare settings. This study aims to understand the relationships between phenotypic diversity and the functional metabolic landscape of P. aeruginosa clinical isolates. To better understand the metabolic repertoire of P. aeruginosa in infection, we deeply profiled a representative set from a library of 971 clinical P. aeruginosa isolates with corresponding patient metadata and bacterial phenotypes. The genotypic clustering based on whole-genome sequencing of the isolates, multi-locus sequence types, and the phenotypic clustering generated from a multi-parametric analysis were compared to each other to assess the genotype-phenotype correlation. Genome-scale metabolic network reconstructions were developed for each isolate through amendments to an existing PA14 network reconstruction. These network reconstructions show diverse metabolic functionalities and enhance the collective P. aeruginosa pangenome metabolic repertoire. Characterizing this rich set of clinical P. aeruginosa isolates allows for a deeper understanding of the genotypic and metabolic diversity of the pathogen in a clinical setting and lays a foundation for further investigation of the metabolic landscape of this pathogen and host-associated metabolic differences during infection.

Comparative Genomic Analyses of Virulence and Antimicrobial Resistance in Citrobacter werkmanii, an Emerging Opportunistic Pathogen

Citrobacter werkmanii is an emerging and opportunistic human pathogen found in developing countries and is a causative agent of wound, urinary tract, and blood infections. The present study conducted comparative genomic analyses of a C. werkmanii strain collection from diverse geographical locations and sources to identify the relevant virulence and antimicrobial resistance genes. Pangenome analyses divided the examined C. werkmanii strains into five distinct clades; the subsequent classification identified genes with functional roles in carbohydrate and general metabolism for the core genome and genes with a role in secretion, adherence, and the mobilome for the shell and cloud genomes. A maximum-likelihood phylogenetic tree with a heatmap, showing the virulence and antimicrobial genes' presence or absence, demonstrated the presence of genes with functional roles in secretion systems, adherence, enterobactin, and siderophore among the strains belonging to the different clades. C. werkmanii strains in clade V, predominantly from clinical sources, harbored genes implicated in type II and type Vb secretion systems as well as multidrug resistance to aminoglycoside, beta-lactamase, fluoroquinolone, phenicol, trimethoprim, macrolides, sulfonamide, and tetracycline. In summary, these comparative genomic analyses have demonstrated highly pathogenic and multidrug-resistant genetic profiles in C. werkmanii strains, indicating a virulence potential for this commensal and opportunistic human pathogen.

Comparative Genomics of Pseudomonas aeruginosa Strains Isolated from Different Ecological Niches

The Pseudomonas aeruginosa genome can change to adapt to different ecological niches. We compared four genomes from a Mexican hospital and 59 genomes from GenBank from different niches, such as urine, sputum, and environmental. The ST analysis showed that high-risk STs (ST235, ST773, and ST27) were present in the genomes of the three niches from GenBank, and the STs of Mexican genomes (ST167, ST2731, and ST549) differed from the GenBank genomes. Phylogenetic analysis showed that the genomes were clustering according to their ST and not their niche. When analyzing the genomic content, we observed that environmental genomes had genes involved in adapting to the environment not found in the clinics and that their mechanisms of resistance were mutations in antibiotic resistance-related genes. In contrast, clinical genomes from GenBank had resistance genes, in mobile/mobilizable genetic elements in the chromosome, except for the Mexican genomes that carried them mostly in plasmids. This was related to the presence of CRISPR-Cas and anti-CRISPR; however, Mexican strains only had plasmids and CRISPR-Cas. bla_OXA-488 (a variant of bla_OXA50) with higher activity against carbapenems was more prevalent in sputum genomes. The virulome analysis showed that exoS was most prevalent in the genomes of urinary samples and exoU and pldA in sputum samples. This study provides evidence regarding the genetic variability among P. aeruginosa isolated from different niches.

Application of prophage sequence analysis to investigate a disease outbreak involving Salmonella Adjame, a rare serovar and implications for the population structure

Introduction

Outbreak investigation of foodborne salmonellosis is hindered when the food source is contaminated by multiple strains of Salmonella, creating difficulties matching an incriminated organism recovered from patients with the specific strain in the suspect food. An outbreak of the rare Salmonella Adjame was caused by multiple strains of the organism as revealed by single-nucleotide polymorphism (SNP) variation. The use of highly discriminatory prophage analysis to characterize strains of Salmonella should enable a more precise strain characterization and aid the investigation of foodborne salmonellosis.

Methods

We have carried out genomic analysis of S. Adjame strains recovered during the course of a recent outbreak and compared them with other strains of the organism (n = 38 strains), using SNPs to evaluate strain differences present in the core genome, and prophage sequence typing (PST) to evaluate the accessory genome. Phylogenetic analyses were performed using both total prophage content and conserved prophages.

Results

The PST analysis of the S. Adjame isolates showed a high degree of strain heterogeneity. We observed small clusters made up of 2-6 isolates (n = 27) and singletons (n = 11) in stark contrast with the three clusters observed by SNP analysis. In total, we detected 24 prophages of which only four were highly prevalent, namely: Entero_p88 (36/38 strains), Salmon_SEN34 (35/38 strains), Burkho_phiE255 (33/38 strains) and Edward_GF (28/38 strains). Despite the marked strain diversity seen with prophage analysis, the distribution of the four most common prophages matched the clustering observed using core genome.

Discussion

Mutations in the core and accessory genomes of S. Adjame have shed light on the evolutionary relationships among the Adjame strains and demonstrated a convergence of the variations observed in both fractions of the genome. We conclude that core and accessory genomes analyses should be adopted in foodborne bacteria outbreak investigations to provide a more accurate strain description and facilitate reliable matching of isolates from patients and incriminated food sources. The outcomes should translate to a better understanding of the microbial population structure and an 46 improved source attribution in foodborne illnesses.

Distribution and spread of the mobilised RND efflux pump gene cluster tmexCD-toprJ in clinical Gram-negative bacteria: a molecular epidemiological study

BACKGROUND

TMexCD1-TOprJ1, which is associated with phenotypic resistance to multiple classes of antibiotics, is a transmissible resistance-nodulation-division (RND) family efflux pump. However, the prevalence and genomic and phenotypic characteristics of clinical isolates with this important resistance determinant are poorly understood. In this study, we aimed to survey tmexCD-toprJ among clinical Gram-negative isolates collected from hospitals in China between 1991 and 2020 and characterise tmexCD-toprJ-positive clinical isolates.

METHODS

We conducted online data retrieval and active nationwide surveillance in China to screen tmexCD-toprJ-positive strains. We characterised tmexCD-toprJ-positive clinical strains for their antimicrobial susceptibility, genetic and functional characteristics, and the potential inter-species transmission route of tmexCD-toprJ with whole genome sequencing and bioinformatics analyses. The function of tmexCD-toprJ in Pseudomonas sp and Proteus sp was investigated by tmexD gene knockdown using an isopropylthio-β-galactoside-inducible CRISPR interference system.

FINDINGS

Data retrieval obtained 53 strains carrying tmexCD-toprJ, comprising 32 Pseudomonas spp, 11 Klebsiella pneumoniae, one Aeromonas spp, one Citrobacter freundii, and one uncultured bacterium from diverse niches. 48 (0·64%) of 7517 clinical isolates from China, including seven Klebsiella spp, one Proteus mirabilis, and 40 Pseudomonas spp, carried tmexCD-toprJ. These isolates exhibited multidrug resistance phenotypes and co-harboured resistance genes, such as mcr and carbapenemases genes. tmexCD-toprJ was encoded on both plasmids and chromosomes in all Klebsiella spp that carried plasmid-borne tmexCD-toprJ (n=7), P mirabilis carried chromosome-borne tmexCD-toprJ, and Pseudomonas spp carried either plasmid-borne (n=19) or chromosome-borne (n=21) ones. tmexCD-toprJ had undergone clonal and horizontal transmission among clinical pathogens. Eight different types of genetic context of tmexCD-toprJ were identified, each of which was associated with different mobile elements, including IntI, IS6100, TnAs1-like, ISRor5, ISVsa3, ISCfr-like, Tn5393, and IS222-like, which might facilitate its transmission. Knockdown of tmexD led to a four times decrease in tigecycline minimum inhibitory concentrations in both Pseudomonas spp and Proteus spp.

INTERPRETATION

Our study provides evidence to suggest that tmexCD-toprJ contributes to the antimicrobial resistance phenotypes in different bacterial species. tmexCD-toprJ has disseminated among diverse species of clinical pathogens, which warrants timely monitoring in clinical pathogens.

FUNDING

National Natural Science Foundation of China, Guangdong Major Project of Basic and Applied Basic Research, Natural Science Foundation of Jiangsu Province.

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.

CRISPR-Cas systems restrict horizontal gene transfer in Pseudomonas aeruginosa

CRISPR-Cas systems provide bacteria and archaea with an adaptive immune system that targets foreign DNA. However, the xenogenic nature of immunity provided by CRISPR-Cas raises the possibility that these systems may constrain horizontal gene transfer. Here we test this hypothesis in the opportunistic pathogen Pseudomonas aeruginosa, which has emerged as an important model system for understanding CRISPR-Cas function. Across the diversity of P. aeruginosa, active CRISPR-Cas systems are associated with smaller genomes and higher GC content, suggesting that CRISPR-Cas inhibits the acquisition of foreign DNA. Although phage is the major target of CRISPR-Cas spacers, more than 80% of isolates with an active CRISPR-Cas system have spacers that target integrative conjugative elements (ICE) or the conserved conjugative transfer machinery used by plasmids and ICE. Consistent with these results, genomes containing active CRISPR-Cas systems harbour a lower abundance of both prophage and ICE. Crucially, spacers in genomes with active CRISPR-Cas systems map to ICE and phage that are integrated into the chromosomes of closely related genomes lacking CRISPR-Cas immunity. We propose that CRISPR-Cas acts as an important constraint to horizontal gene transfer, and the evolutionary mechanisms that ensure its maintenance or drive its loss are key to the ability of this pathogen to adapt to new niches and stressors.

Estimation of pathogenic potential of an environmental Pseudomonas aeruginosa isolate using comparative genomics

The isolation and sequencing of new strains of Pseudomonas aeruginosa created an extensive dataset of closed genomes. Many of the publicly available genomes are only used in their original publication while additional in silico information, based on comparison to previously published genomes, is not being explored. In this study, we defined and investigated the genome of the environmental isolate P. aeruginosa KRP1 and compared it to more than 100 publicly available closed P. aeruginosa genomes. By using different genomic island prediction programs, we could identify a total of 17 genomic islands and 8 genomic islets, marking the majority of the accessory genome that covers ~ 12% of the total genome. Based on intra-strain comparisons, we are able to predict the pathogenic potential of this environmental isolate. It shares a substantial amount of genomic information with the highly virulent PSE9 and LESB58 strains. For both of these, the increased virulence has been directly linked to their accessory genome before. Hence, the integrated use of previously published data can help to minimize expensive and time consuming wetlab work to determine the pathogenetic potential.

The Basis for Natural Multiresistance to Phage in Pseudomonas aeruginosa

Pseudomonas aeruginosa is responsible for long-term infections and is particularly resistant to treatments when hiding inside the extracellular matrix or biofilms. Phage therapy might represent an alternative to antibiotic treatment, but up to 10% of clinical strains appear to resist multiple phages. We investigated the characteristics of P. aeruginosa clinical strains naturally resistant to phages and compared them to highly susceptible strains. The phage-resistant strains were defective in lipopolysaccharide (LPS) biosynthesis, were nonmotile and displayed an important degree of autolysis, releasing phages and pyocins. Complete genome sequencing of three resistant strains showed the existence of a large accessory genome made of multiple insertion elements, genomic islands, pyocins and prophages, including two phages performing lateral transduction. Mutations were found in genes responsible for the synthesis of LPS and/or type IV pilus, the major receptors for most phages. CRISPR-Cas systems appeared to be absent or inactive in phage-resistant strains, confirming that they do not play a role in the resistance to lytic phages but control the insertion of exogenous sequences. We show that, despite their apparent weakness, the multiphage-resistant strains described in this study displayed selective advantages through the possession of various functions, including weapons to eliminate other strains of the same or closely related species.

A Core Genome Multilocus Sequence Typing Scheme for Pseudomonas aeruginosa

Pseudomonas aeruginosa is a ubiquitous microorganism and an important opportunistic pathogen responsible for a broad spectrum of infections mainly in immunosuppressed and critically ill patients. Molecular investigations traditionally rely on pulsed field gel electrophoresis (PFGE) and multilocus sequence typing (MLST). In this work we propose a core genome multilocus sequence typing (cgMLST) scheme for P. aeruginosa, a methodology that combines traditional MLST principles with whole genome sequencing data. All publicly available complete P. aeruginosa genomes, representing the diversity of this species, were used to establish a cgMLST scheme targeting 2,653 genes. The scheme was then tested using genomes available at contig, chromosome and scaffold levels. The proposed cgMLST scheme for P. aeruginosa typed over 99% (2,314/2,325) of the genomes available for this study considering at least 95% of the cgMLST target genes present. The absence of a certain number gene targets at the threshold considered for both the creation and validation steps due to low genome sequence quality is possibly the main reason for this result. The cgMLST scheme was compared with previously published whole genome single nucleotide polymorphism analysis for the characterization of the population structure of the epidemic clone ST235 and results were highly similar. In order to evaluate the typing resolution of the proposed scheme, collections of isolates belonging to two important STs associated with cystic fibrosis, ST146 and ST274, were typed using this scheme, and ST235 isolates associated with an outbreak were evaluated. Besides confirming the relatedness of all the isolates, earlier determined by MLST, the higher resolution of cgMLST denotes that it may be suitable for surveillance programs, overcoming possible shortcomings of classical MLST. The proposed scheme is publicly available at: https://github.com/BioinformaticsHIAEMolecularMicrobiology/cgMLST-Pseudomonas-aeruginosa.