Genome sequence of Pseudomonas aeruginosa PAO1161, a PAO1 derivative with the ICEPae1161 integrative and conjugative element

Genomic and phylogenomic characterization of carbapenem-resistant Pseudomonas aeruginosa 'high-risk' clone O4/ExoS+/ST654 circulating in Chilean hospitals

INTRODUCTION

Carbapenem-resistant Pseudomonas aeruginosa (CRPA) is a serious threat to public health. Globally, carbapenemases-producing CRPA isolates mainly belong to 'high-risk' clones; however, the molecular epidemiology of CRPA isolates circulating in Chile are scarce, where this pathogen is the main aetiological agent of ventilator-associated pneumonia.

OBJECTIVES

To characterize the phylogenomics and molecular features of ST654 CRPA isolates collected in Chile between 2016 - 2022.

METHODS

89 CRPA isolates collected in different Chilean hospitals from clinical specimens between 2005 and 2022 were analyzed. Antibiotic susceptibility tests and carbapenemases production were carried out on the CRPA ST654 isolates. Also, they were subjected to whole-genome sequencing (WGS) from which in silico analyses were performed.

RESULTS

Thirty-four strains (38.2%) belonged to the ST654 'high risk' clone, being the most predominant lineage of the collection. Most of these isolates belonged to a sub-clade including KPC-producers that also clustered with strains from Argentina and the USA, whereas few VIM and NDM co-producers clustered in two different smaller sub-clades. The isolates exhibited a broad resistome encompassing genes mediating resistance to several other clinically relevant drugs. Additionally, all the 34 ST654 isolates were ExoS+ as a virulence factor and associated to the O4-serotype.

CONCLUSIONS

Our report represents the most comprehensive phylogenomic study of CRPA 'high risk' clone ST654 to date. Our analyses suggest that this lineage is undergoing a divergent evolutionary path in Chile, since most of the isolates were KPC-producers and were O4-serotype, differing from previous descriptions, which underline the relevance of performing molecular surveillance on this pathogen.

Isolation and characterization of novel plasmid-dependent phages infecting bacteria carrying diverse conjugative plasmids

IMPORTANCE

This work was undertaken because plasmid-dependent phages can reduce the prevalence of conjugative plasmids and can be leveraged to prevent the acquisition and dissemination of ARGs by bacteria. The two novel phages described in this study, Lu221 and Hi226, can infect Escherichia coli, Salmonella enterica, Kluyvera sp. and Enterobacter sp. carrying conjugative plasmids. This was verified with plasmids carrying resistance determinants and belonging to the most common plasmid families among Gram-negative pathogens. Therefore, the newly isolated phages could have the potential to help control the spread of ARGs and thus help combat the antimicrobial resistance crisis.

Robust ParB Binding to Half-parS Sites in Pseudomonas aeruginosa-A Mechanism for Retaining ParB on the Nucleoid?

Chromosome segregation in Pseudomonas aeruginosa is assisted by the tripartite ParAB-parS system, composed of an ATPase (ParA), a DNA-binding protein (ParB) and its target parS sequence(s). ParB forms a nucleoprotein complex around four parSs (parS1-parS4) that overlaps oriC and facilitates relocation of newly synthesized ori domains inside the cells by ParA. Remarkably, ParB of P. aeruginosa also binds to numerous heptanucleotides (half-parSs) scattered in the genome. Here, using chromatin immunoprecipitation-sequencing (ChIP-seq), we analyzed patterns of ParB genome occupancy in cells growing under conditions of coupling or uncoupling between replication and cell division processes. Interestingly, a dissipation of ParB-parS complexes and a shift of ParB to half-parSs were observed during the transition from the exponential to stationary phase of growth on rich medium, suggesting the role of half-parSs in retaining ParB on the nucleoid within non-dividing P. aeruginosa cells. The ChIP-seq analysis of strains expressing ParB variants unable to dislocate from parSs showed that the ParB spreading ability is not required for ParB binding to half-parSs. Finally, a P. aeruginosa strain with mutated 25 half-parSs of the highest affinity towards ParB was constructed and analyzed. It showed altered ParB coverage of the oriC region and moderate changes in gene expression. Overall, this study characterizes a novel aspect of conserved bacterial chromosome segregation machinery.

Diverse Partners of the Partitioning ParB Protein in Pseudomonas aeruginosa

In the majority of bacterial species, the tripartite ParAB-parS system, composed of an ATPase (ParA), a DNA-binding protein (ParB), and its target parS sequence(s), assists in the chromosome partitioning. ParB forms large nucleoprotein complexes at parS(s), located in the vicinity of origin of chromosomal replication (oriC), which after replication are subsequently positioned by ParA in cell poles. Remarkably, ParA and ParB participate not only in the chromosome segregation but through interactions with various cellular partners they are also involved in other cell cycle-related processes, in a species-specific manner. In this work, we characterized Pseudomonas aeruginosa ParB interactions with the cognate ParA, showing that the N-terminal motif of ParB is required for these interactions, and demonstrated that ParAB-parS-mediated rapid segregation of newly replicated ori domains prevented structural maintenance of chromosome (SMC)-mediated cohesion of sister chromosomes. Furthermore, using proteome-wide techniques, we have identified other ParB partners in P. aeruginosa, which encompass a number of proteins, including the nucleoid-associated proteins NdpA(PA3849) and NdpA2, MinE (PA3245) of Min system, and transcriptional regulators and various enzymes, e.g., CTP synthetase (PA3637). Among them are also NTPases PA4465, PA5028, PA3481, and FleN (PA1454), three of them displaying polar localization in bacterial cells. Overall, this work presents the spectrum of P. aeruginosa ParB partners and implicates the role of this protein in the cross-talk between chromosome segregation and other cellular processes. IMPORTANCE In Pseudomonas aeruginosa, a Gram-negative pathogen causing life-threatening infections in immunocompromised patients, the ParAB-parS system is involved in the precise separation of newly replicated bacterial chromosomes. In this work, we identified and characterized proteins interacting with partitioning protein ParB. We mapped the domain of interactions with its cognate ParA partner and showed that ParB-ParA interactions are crucial for the chromosome segregation and for proper SMC action on DNA. We also demonstrated ParB interactions with other DNA binding proteins, metabolic enzymes, and NTPases displaying polar localization in the cells. Overall, this study uncovers novel players cooperating with the chromosome partition system in P. aeruginosa, supporting its important regulatory role in the bacterial cell cycle.

Whole-genome sequencing reveals high-risk clones of Pseudomonas aeruginosa in Guangdong, China

The ever-increasing prevalence of infections produced by multidrug-resistant or extensively drug-resistant Pseudomonas aeruginosa is commonly linked to a limited number of aptly-named epidemical 'high-risk clones' that are widespread among and within hospitals worldwide. The emergence of new potential high-risk clone strains in hospitals highlights the need to better and further understand the underlying genetic mechanisms for their emergence and success. P. aeruginosa related high-risk clones have been sporadically found in China, their genome sequences have rarely been described. Therefore, the large-scale sequencing of multidrug-resistance high-risk clone strains will help us to understand the emergence and transmission of antibiotic resistances in P. aeruginosa high-risk clones. In this study, 212 P. aeruginosa strains were isolated from 2 tertiary hospitals within 3 years (2018-2020) in Guangdong Province, China. Whole-genome sequencing, multi-locus sequence typing (MLST) and antimicrobial susceptibility testing were applied to analyze the genomic epidemiology of P. aeruginosa in this region. We found that up to 130 (61.32%) of the isolates were shown to be multidrug resistant, and 196 (92.45%) isolates were Carbapenem-Resistant Pseudomonas aeruginosa. MLST analysis demonstrated high diversity of sequence types, and 18 reported international high-risk clones were identified. Furthermore, we discovered the co-presence of exoU and exoS genes in 5 collected strains. This study enhances insight into the regional research of molecular epidemiology and antimicrobial resistance of P. aeruginosa in China. The high diversity of clone types and regional genome characteristics can serve as a theoretical reference for public health policies and help guide measures for the prevention and control of P. aeruginosa resistance.

Functional Characterization of TetR-like Transcriptional Regulator PA3973 from Pseudomonas aeruginosa

Pseudomonas aeruginosa, a human opportunistic pathogen, is a common cause of nosocomial infections. Its ability to survive under different conditions relies on a complex regulatory network engaging transcriptional regulators controlling metabolic pathways and capabilities to efficiently use the available resources. P. aeruginosa PA3973 encodes an uncharacterized TetR family transcriptional regulator. In this study, we applied a transcriptome profiling (RNA-seq), genome-wide identification of binding sites using ChIP-seq, as well as the phenotype analyses to unravel the biological role of PA3973. Transcriptional profiling of P. aeruginosa PAO1161 overexpressing PA3973 showed changes in the mRNA level of 648 genes. Concomitantly, ChIP-seq analysis identified more than 300 PA3973 binding sites in the P. aeruginosa genome. A 13 bp sequence motif was indicated as the binding site of PA3973. The PA3973 regulon encompasses the PA3972-PA3971 genes encoding a probable acyl-CoA dehydrogenase and a thioesterase. In vitro analysis showed PA3973 binding to PA3973p. Accordingly, the lack of PA3973 triggered increased expression of PA3972 and PA3971. The ∆PA3972-71 PAO1161 strain demonstrated impaired growth in the presence of stress-inducing agents hydroxylamine or hydroxyurea, thus suggesting the role of PA3972-71 in pathogen survival upon stress. Overall our results showed that TetR-type transcriptional regulator PA3973 has multiple binding sites in the P. aeruginosa genome and influences the expression of diverse genes, including PA3972-PA3971, encoding proteins with a proposed role in stress response.

The Lrp/AsnC-Type Regulator PA2577 Controls the EamA-like Transporter Gene PA2576 in Pseudomonas aeruginosa

The regulatory network of gene expression in Pseudomonas aeruginosa, an opportunistic human pathogen, is very complex. In the PAO1 reference strain, about 10% of genes encode transcriptional regulators, many of which have undefined regulons and unknown functions. The aim of this study is the characterization of PA2577 protein, a representative of the Lrp/AsnC family of transcriptional regulators. This family encompasses proteins involved in the amino acid metabolism, regulation of transport processes or cell morphogenesis. The transcriptome profiling of P. aeruginosa cells with mild PA2577 overproduction revealed a decreased expression of the PA2576 gene oriented divergently to PA2577 and encoding an EamA-like transporter. A gene expression analysis showed a higher mRNA level of PA2576 in P. aeruginosa ΔPA2577, indicating that PA2577 acts as a repressor. Concomitantly, ChIP-seq and EMSA assays confirmed strong interactions of PA2577 with the PA2577/PA2576 intergenic region. Additionally, phenotype microarray analyses indicated an impaired metabolism of ΔPA2576 and ΔPA2577 mutants in the presence of polymyxin B, which suggests disturbances of membrane functions in these mutants. We show that PA2576 interacts with two proteins, PA5006 and PA3694, with a predicted role in lipopolysaccharide (LPS) and membrane biogenesis. Overall, our results indicate that PA2577 acts as a repressor of the PA2576 gene coding for the EamA-like transporter and may play a role in the modulation of the cellular response to stress conditions, including antimicrobial peptides, e.g., polymyxin B.

carP, encoding a Ca2+-regulated putative phytase, is evolutionarily conserved in Pseudomonas aeruginosa and has potential as a biomarker

Pseudomonas aeruginosa infects patients with cystic fibrosis, burns, wounds and implants. Previously, our group showed that elevated Ca²⁺ positively regulates the production of several virulence factors in P. aeruginosa, such as biofilm formation, production of pyocyanin and secreted proteases. We have identified a Ca²⁺-regulated β-propeller putative phytase, CarP, which is required for Ca²⁺ tolerance, regulation of the intracellular Ca²⁺ levels, and plays a role in Ca²⁺ regulation of P. aeruginosa virulence. Here, we studied the conservation of carP sequence and its occurrence in diverse phylogenetic groups of bacteria. In silico analysis revealed that carP and its two paralogues PA2017 and PA0319 are primarily present in P. aeruginosa and belong to the core genome of the species. We identified 155 single nucleotide alterations within carP, 42 of which lead to missense mutations with only three that affected the predicted 3D structure of the protein. PCR analyses with carP-specific primers detected P. aeruginosa specifically in 70 clinical and environmental samples. Sequence comparison demonstrated that carP is overall highly conserved in P. aeruginosa isolated from diverse environments. Such evolutionary preservation of carP illustrates its importance for P. aeruginosa adaptations to diverse environments and demonstrates its potential as a biomarker.

The LysR-Type Transcriptional Regulator BsrA (PA2121) Controls Vital Metabolic Pathways in Pseudomonas aeruginosa

Pseudomonas aeruginosa, a facultative human pathogen causing nosocomial infections, has complex regulatory systems involving many transcriptional regulators. LTTR (LysR-Type Transcriptional Regulator) family proteins are involved in the regulation of various processes, including stress responses, motility, virulence, and amino acid metabolism. The aim of this study was to characterize the LysR-type protein BsrA (PA2121), previously described as a negative regulator of biofilm formation in P. aeruginosa. Genome wide identification of BsrA binding sites using chromatin immunoprecipitation and sequencing analysis revealed 765 BsrA-bound regions in the P. aeruginosa PAO1161 genome, including 367 sites in intergenic regions. The motif T-N₁₁-A was identified within sequences bound by BsrA. Transcriptomic analysis showed altered expression of 157 genes in response to BsrA excess; of these, 35 had a BsrA binding site within their promoter regions, suggesting a direct influence of BsrA on the transcription of these genes. BsrA-repressed loci included genes encoding proteins engaged in key metabolic pathways such as the tricarboxylic acid cycle. The panel of loci possibly directly activated by BsrA included genes involved in pilus/fimbria assembly, as well as secretion and transport systems. In addition, DNA pull-down and regulatory analyses showed the involvement of PA2551, PA3398, and PA5189 in regulation of bsrA expression, indicating that this gene is part of an intricate regulatory network. Taken together, these findings reveal the existence of a BsrA regulon, which performs important functions in P. aeruginosa.

IMPORTANCE This study shows that BsrA, a LysR-type transcriptional regulator from Pseudomonas aeruginosa, previously identified as a repressor of biofilm synthesis, is part of an intricate global regulatory network. BsrA acts directly and/or indirectly as the repressor and/or activator of genes from vital metabolic pathways (e.g., pyruvate, acetate, and tricarboxylic acid cycle) and is involved in control of transport functions and the formation of surface appendages. Expression of the bsrA gene is increased in the presence of antibiotics, which suggests its induction in response to stress, possibly reflecting the need to redirect metabolism under stressful conditions. This is particularly relevant for the treatment of infections caused by P. aeruginosa. In summary, the findings of this study demonstrate that the BsrA regulator performs important roles in carbon metabolism, biofilm formation, and antibiotic resistance in P. aeruginosa.

The AraC-Type Transcriptional Regulator GliR (PA3027) Activates Genes of Glycerolipid Metabolism in Pseudomonas aeruginosa

Pseudomonas aeruginosa encodes a large set of transcriptional regulators (TRs) that modulate and manage cellular metabolism to survive in variable environmental conditions including that of the human body. The AraC family regulators are an abundant group of TRs in bacteria, mostly acting as gene expression activators, controlling diverse cellular functions (e.g., carbon metabolism, stress response, and virulence). The PA3027 protein from P. aeruginosa has been classified in silico as a putative AraC-type TR. Transcriptional profiling of P. aeruginosa PAO1161 overexpressing PA3027 revealed a spectacular increase in the mRNA levels of PA3026-PA3024 (divergent to PA3027), PA3464, and PA3342 genes encoding proteins potentially involved in glycerolipid metabolism. Concomitantly, chromatin immunoprecipitation-sequencing (ChIP-seq) analysis revealed that at least 22 regions are bound by PA3027 in the PAO1161 genome. These encompass promoter regions of PA3026, PA3464, and PA3342, showing the major increase in expression in response to PA3027 excess. In Vitro DNA binding assay confirmed interactions of PA3027 with these regions. Furthermore, promoter-reporter assays in a heterologous host showed the PA3027-dependent activation of the promoter of the PA3026-PA3024 operon. Two motifs representing the preferred binding sites for PA3027, one localized upstream and one overlapping with the -35 promoter sequence, were identified in PA3026p and our data indicate that both motifs are required for full activation of this promoter by PA3027. Overall, the presented data show that PA3027 acts as a transcriptional regulator in P. aeruginosa, activating genes likely engaged in glycerolipid metabolism. The GliR name, from a glycerolipid metabolism regulator, is proposed for PA3027 of P. aeruginosa.

Comparative genome analysis of multidrug-resistant Pseudomonas aeruginosa JNQH-PA57, a clinically isolated mucoid strain with comprehensive carbapenem resistance mechanisms

BACKGROUND

The prevalence of clinical multidrug-resistant (MDR) Pseudomonas aeruginosa has been increasing rapidly worldwide over the years and responsible for a wide range of acute and chronic infections with high mortalities. Although hundreds of complete genomes of clinical P. aeruginosa isolates have been sequenced, only a few complete genomes of mucoid strains are available, limiting a comprehensive understanding of this important group of opportunistic pathogens. Herein, the complete genome of a clinically isolated mucoid strain P. aeruginosa JNQH-PA57 was sequenced and assembled using Illumina and Oxford nanopore sequencing technologies. Genomic features, phylogenetic relationships, and comparative genomics of this pathogen were comprehensively analyzed using various bioinformatics tools. A series of phenotypic and molecular-genetic tests were conducted to investigate the mechanisms of carbapenem resistance in this strain.

RESULTS

Several genomic features of MDR P. aeruginosa JNQH-PA57 were identified based on the whole-genome sequencing. We found that the accessory genome of JNQH-PA57 including several prophages, genomic islands, as well as a PAPI-1 family integrative and conjugative element (ICE), mainly contributed to the larger genome of this strain (6,747,067 bp) compared to other popular P. aeruginosa strains (with an average genome size of 6,445,223 bp) listed in Pseudomonas Genome Database. Colony morphology analysis and biofilm crystal staining assay respectively demonstrated an enhanced alginate production and a thicker biofilm formation capability of JNQH-PA57. A deleted mutation at nt 424 presented in mucA gene, resulted in the upregulated expression of a sigma-factor AlgU and a GDP mannose dehydrogenase AlgD, which might explain the mucoid phenotype of this strain. As for the carbapenem resistance mechanisms, our results revealed that the interplay between impaired OprD porin, chromosomal β-lactamase OXA-488 expression, MexAB-OprM and MexXY-OprM efflux pumps overexpression, synergistically with the alginates-overproducing protective biofilm, conferred the high carbapenem resistance to P. aeruginosa JNQH-PA57.

CONCLUSION

Based on the genome analysis, we could demonstrate that the upregulated expression of algU and algD, which due to the truncation variant of MucA, might account for the mucoid phenotype of JNQH-PA57. Moreover, the resistance to carbapenem in P. aeruginosa JNQH-PA57 is multifactorial. The dataset presented in this study provided an essential genetic basis for the comprehensive cognition of the physiology, pathogenicity, and carbapenem resistance mechanisms of this clinical mucoid strain.

The MarR-Type Regulator PA3458 Is Involved in Osmoadaptation Control in Pseudomonas aeruginosa

Pseudomonas aeruginosa is a facultative human pathogen, causing acute and chronic infections that are especially dangerous for immunocompromised patients. The eradication of P. aeruginosa is difficult due to its intrinsic antibiotic resistance mechanisms, high adaptability, and genetic plasticity. The bacterium possesses multilevel regulatory systems engaging a huge repertoire of transcriptional regulators (TRs). Among these, the MarR family encompasses a number of proteins, mainly acting as repressors, which are involved in response to various environmental signals. In this work, we aimed to decipher the role of PA3458, a putative MarR-type TR from P. aeruginosa. Transcriptional profiling of P. aeruginosa PAO1161 overexpressing PA3458 showed changes in the mRNA level of 133 genes; among them, 100 were down-regulated, suggesting the repressor function of PA3458. Concomitantly, ChIP-seq analysis identified more than 300 PA3458 binding sites in P. aeruginosa. The PA3458 regulon encompasses genes involved in stress response, including the PA3459–PA3461 operon, which is divergent to PA3458. This operon encodes an asparagine synthase, a GNAT-family acetyltransferase, and a glutamyl aminopeptidase engaged in the production of N-acetylglutaminylglutamine amide (NAGGN), which is a potent bacterial osmoprotectant. We showed that PA3458-mediated control of PA3459–PA3461 expression is required for the adaptation of P. aeruginosa growth in high osmolarity. Overall, our data indicate that PA3458 plays a role in osmoadaptation control in P. aeruginosa.

Tn6603, a Carrier of Tn5053 Family Transposons, Occurs in the Chromosome and in a Genomic Island of Pseudomonas aeruginosa Clinical Strains

Transposons of the Pseudomonasaeruginosa accessory gene pool contribute to phenotype and to genome plasticity. We studied local P. aeruginosa strains to ascertain the encroachment of mer-type res site hunter transposons into clinical settings and their associations with other functional modules. Five different Tn5053 family transposons were detected, all chromosomal. Some were solitary elements; one was in res of Tn1013^#, a relative of a reported carrier of int-type res site hunters (class 1 integrons), but most were in res of Tn6603, a new Tn501-related transposon of unknown phenotype. Most of the Tn6603::Tn elements, and some Tn6603 and Tn6603::Tn elements found in GenBank sequences, were at identical sites in an hypothetical gene of P. aeruginosa genomic island PAGI-5v. The island in clonally differing strains was at either of two tRNA^Lys loci, suggesting lateral transfer to these sites. This observation is consistent with the membership of the prototype PAGI-5 island to the ICE family of mobile genetic elements. Additionally, the res site hunters in the nested transposons occupied different positions in the Tn6603 carrier. This suggested independent insertion events on five occasions at least. Tn5053 family members that were mer-/tni-defective were found in Tn6603- and Tn501-like carriers in GenBank sequences of non-clinical Pseudomonas spp. The transposition events in these cases presumably utilized tni functions in trans, as can occur with class 1 integrons. We suggest that in the clinical context, P. aeruginosa strains that carry Tn6603 alone or in PAGI-5v can serve to disseminate functional res site hunters; these in turn can provide the requisite trans-acting tni functions to assist in the dissemination of class 1 integrons, and hence of their associated antibiotic resistance determinants.

Pseudomonas aeruginosa mexT is an indicator of PAO1 strain integrity

Laboratory research with Pseudomonas aeruginosa commonly involves the prototype strain PAO1. There is continued concern that PAO1 sublines maintained and propagated in the same laboratory or different laboratories exhibit genetic and phenotypic variability that may affect the reproducibility and validity of research. Whole-genome sequencing and other research identified the mexT locus as a mutational hotspot, but the explication of the diverse mutations present in the various sublines and consequences remained rather cursory. Here we present evidence that MexT sequence diversity is a predictor of PAO1 lineage integrity and define the protein's prototype sequence.

Interaction of ArmZ with the DNA-binding domain of MexZ induces expression of mexXY multidrug efflux pump genes and antimicrobial resistance in Pseudomonas aeruginosa

Multidrug efflux pumps play an important role in antibiotic resistance in bacteria. In Pseudomonas aeruginosa, MexXY pump provides intrinsic resistance to many antimicrobials, including aminoglycosides. The expression of mexXY operon is negatively regulated by MexZ repressor. The repression is alleviated in response to the antibiotic-induced ribosome stress, which results in increased synthesis of anti-repressor ArmZ, interacting with MexZ. The molecular mechanism of MexZ inactivation by ArmZ is not known. Here, we showed that the N-terminal part of MexZ, encompassing the DNA-binding domain, is required for interaction with ArmZ. Using the bacterial two hybrid system based mutant screening and pull-down analyses we identified substitutions in MexZ that diminished (R3S, K6E, R13H) or completely impaired (K53E) the interaction with ArmZ without blocking MexZ activity as a transcriptional repressor. Introduction of corresponding mexZ missense mutations into P aeruginosa PAO1161 chromosome impaired (mexZ K6E, mexZ R13H) or blocked (mexZ K53E) tetracycline mediated induction of mexY expression. Concomitantly, PAO1161 mexZ K53E strain was more susceptible to aminoglycosides. The identified residues are highly conserved in MexZ-like transcriptional regulators found in bacterial genomes encoding both MexX/MexY/MexZ and ArmZ/PA5470 orthologs, suggesting that a similar mechanism may contribute to induction of efflux mediated resistance in other bacterial species. Overall, our data shed light on the molecular mechanism of ArmZ mediated induction of intrinsic antimicrobial resistance in P. aeruginosa.