The Accessory Genome of Pseudomonas aeruginosa

Summary: Pseudomonas aeruginosa strains exhibit significant variability in pathogenicity and ecological flexibility. Such interstrain differences reflect the dynamic nature of the P. aeruginosa genome, which is composed of a relatively invariable “core genome” and a highly variable “accessory genome.” Here we review the major classes of genetic elements comprising the P. aeruginosa accessory genome and highlight emerging themes in the acquisition and functional importance of these elements. Although the precise phenotypes endowed by the majority of the P. aeruginosa accessory genome have yet to be determined, rapid progress is being made, and a clearer understanding of the role of the P. aeruginosa accessory genome in ecology and infection is emerging.     

Characterization of the core and accessory genomes of Pseudomonas aeruginosa using bioinformatic tools Spine and AGEnt

Background

Pseudomonas aeruginosa is an important opportunistic pathogen responsible for many infections in hospitalized and immunocompromised patients. Previous reports estimated that approximately 10% of its 6.6 Mbp genome varies from strain to strain and is therefore referred to as “accessory genome”. Elements within the accessory genome of P. aeruginosa have been associated with differences in virulence and antibiotic resistance. As whole genome sequencing of bacterial strains becomes more widespread and cost-effective, methods to quickly and reliably identify accessory genomic elements in newly sequenced P. aeruginosa genomes will be needed.

Results

We developed a bioinformatic method for identifying the accessory genome of P. aeruginosa. First, the core genome was determined based on sequence conserved among the completed genomes of twelve reference strains using Spine, a software program developed for this purpose. The core genome was 5.84 Mbp in size and contained 5,316 coding sequences. We then developed an in silico genome subtraction program named AGEnt to filter out core genomic sequences from P. aeruginosa whole genomes to identify accessory genomic sequences of these reference strains. This analysis determined that the accessory genome of P. aeruginosa ranged from 6.9-18.0% of the total genome, was enriched for genes associated with mobile elements, and was comprised of a majority of genes with unknown or unclear function. Using these genomes, we showed that AGEnt performed well compared to other publically available programs designed to detect accessory genomic elements. We then demonstrated the utility of the AGEnt program by applying it to the draft genomes of two previously unsequenced P. aeruginosa strains, PA99 and PA103.

Conclusions

The P. aeruginosa genome is rich in accessory genetic material. The AGEnt program accurately identified the accessory genomes of newly sequenced P. aeruginosa strains, even when draft genomes were used. As P. aeruginosa genomes become available at an increasingly rapid pace, this program will be useful in cataloging the expanding accessory genome of this bacterium and in discerning correlations between phenotype and accessory genome makeup. The combination of Spine and AGEnt should be useful in defining the accessory genomes of other bacterial species as well.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-737) contains supplementary material, which is available to authorized users.     

Complete genome sequence of Pseudomonas aeruginosa lytic bacteriophage PA1O which resembles temperate bacteriophage D3112

A novel Pseudomonas aeruginosa lytic bacteriophage (phage), PA1Ø, was isolated, and its genome was sequenced completely. This phage is able to lyse not only P. aeruginosa but also Staphylococcus aureus. Genome analysis of PA1Ø showed that it is similar to a P. aeruginosa temperate phage, D3112, with the exception of the absence of a c repressor-encoding gene, which is known to play a critical role in the maintenance of the lysogenic state of D3112 in P. aeruginosa.     

Complete Genome Sequence of the Bacteriophage YMC01/01/P52 PAE BP,Which Causes Lysis of Verona Integron-Encoded Metallo-β-Lactamase-Producing,Carbapenem-Resistant Pseudomonas aeruginosa

Multidrug-resistant Pseudomonas aeruginosa commonly causes serious nosocomial infections. In this study, a novel lytic bacteriophage belonging to a member of the family Podoviridae, YMC01/01/P52 PAE BP, which infects carbapenem-resistant Pseudomonas aeruginosa, was isolated and characterized. YMC01/01/P52 PAE BP genome was analyzed by whole-genome sequencing and putative function identification. The bacteriophage genome consists of a double-stranded linear DNA genome of 49,381 bp with a GC content of 62.16%.     

Draft Genome Sequence for Pseudomonas aeruginosa Strain PAO579, a Mucoid Derivative of PAO381

Pseudomonas aeruginosa is an opportunistic pathogen that establishes a chronic lung infection in individuals afflicted with cystic fibrosis. Here, we announce the draft genome of P. aeruginosa strain PAO579, an alginate-overproducing derivative of strain PAO381.     

Two copies of bla NDM-1 gene are present in NDM-1 producing Pseudomonas aeruginosa isolates from Serbia

New Delhi metallo-β-lactamase producing Pseudomonas aeruginosa isolates are of special interest since P. aeruginosa is a major cause of nosocomial infections, the treatment of which could now be jeopardized, especially in developing countries. Six additional NDM-1 positive P. aeruginosa clinical isolates belonging to two different genotypes were shown to be plasmid-free. PFGE-hybridization experiments revealed the chromosomal location of the bla _NDM-1 gene. Restriction analysis and hybridization revealed that two copies of the bla _NDM-1 gene are present in the genomes of all tested isolates, as in previously characterized P. aeruginosa MMA83. Moreover, it was shown that increasing imipenem concentration did not have the effect on copy number of the bla _NDM-1 gene in the genome of P. aeruginosa MMA83.     

A Novel Virulence Strategy for Pseudomonas aeruginosa Mediated by an Autotransporter with Arginine-Specific Aminopeptidase Activity

The opportunistic human pathogen, Pseudomonas aeruginosa, is a major cause of infections in chronic wounds, burns and the lungs of cystic fibrosis patients. The P. aeruginosa genome encodes at least three proteins exhibiting the characteristic three domain structure of autotransporters, but much remains to be understood about the functions of these three proteins and their role in pathogenicity. Autotransporters are the largest family of secreted proteins in Gram-negative bacteria, and those characterised are virulence factors. Here, we demonstrate that the PA0328 autotransporter is a cell-surface tethered, arginine-specific aminopeptidase, and have defined its active site by site directed mutagenesis. Hence, we have assigned PA0328 with the name AaaA, for arginine-specific autotransporter of P. aeruginosa. We show that AaaA provides a fitness advantage in environments where the sole source of nitrogen is peptides with an aminoterminal arginine, and that this could be important for establishing an infection, as the lack of AaaA led to attenuation in a mouse chronic wound infection which correlated with lower levels of the cytokines TNFα, IL-1α, KC and COX-2. Consequently AaaA is an important virulence factor playing a significant role in the successful establishment of P. aeruginosa infections.     

Genome Sequence of Pseudomonas aeruginosa Strain AH16, Isolated from a Patient with Chronic Pneumonia in China

Pseudomonas aeruginosa AH16 is a virulent strain isolated from a patient with chronic pneumonia in China. Here, we present a 6.8-Mb (G+C content, 66.13%) assembly of its genome with 6,332 putative coding sequences, which may provide insights into the genomic basis of activity of the clinical P. aeruginosa strain in China.     

Comprehensive MALDI-TOF Biotyping of the Non-Redundant Harvard Pseudomonas aeruginosa PA14 Transposon Insertion Mutant Library

BackgroundPseudomonas aeruginosa is a gram-negative bacterium that is ubiquitously present in the aerobic biosphere. As an antibiotic-resistant facultative pathogen, it is a major cause of hospital-acquired infections. Its rapid and accurate identification is crucial in clinical and therapeutic environments.MethodsIn a large-scale MALDI-TOF mass spectrometry-based screen of the Harvard transposon insertion mutant library of P. aeruginosa strain PA14, intact-cell proteome profile spectra of 5547 PA14 transposon mutants exhibiting a plethora of different phenotypes were acquired and analyzed.ResultsOf all P. aeruginosa PA14 mutant profiles 99.7% were correctly identified as P. aeruginosa with the Biotyper software on the species level. On the strain level, 99.99% of the profiles were mapped to five different individual P. aeruginosa Biotyper database entries. A principal component analysis-based approach was used to determine the most important discriminatory mass features between these Biotyper groups. Although technical replicas were consistently categorized to specific Biotyper groups in 94.2% of the mutant profiles, biological replicas were not, indicating that the distinct proteotypes are affected by growth conditions.ConclusionsThe PA14 mutant profile collection presented here constitutes the largest coherent P. aeruginosa MALDI-TOF spectral dataset publicly available today. Transposon insertions in thousands of different P. aeruginosa genes did not affect species identification from MALDI-TOF mass spectra, clearly demonstrating the robustness of the approach. However, the assignment of the individual spectra to sub-groups proved to be non-consistent in biological replicas, indicating that the differentiation between biotyper groups in this nosocomial pathogen is unassured.     

New temperate Pseudomonas aeruginosa phage, phi297: Specific features of genome structure

The genome structure and some specific features of temperate Pseudomonas aeruginosa phage phi297 are considered. Analysis of sequencing data and genome annotation suggest that the phi297 genome displays a mosaic structure, which has formed through combining gene blocks from bacteria of taxonomically remote groups and/or their phages. The results of a comparison of the phi297 DNA homology level and pattern with the genome sequences of the currently known related P. aeruginosa bacteriophages are interpreted from the perspective of assumed active migration of these phages between different bacterial species.     

Hybrid pathogenicity island PAGI-5 contributes to the highly virulent phenotype of a Pseudomonas aeruginosa isolate in mammals

Most known virulence determinants of Pseudomonas aeruginosa are remarkably conserved in this bacterium's core genome, yet individual strains differ significantly in virulence. One explanation for this discrepancy is that pathogenicity islands, regions of DNA found in some strains but not in others, contribute to the overall virulence of P. aeruginosa. Here we employed a strategy in which the virulence of a panel of P. aeruginosa isolates was tested in mouse and plant models of disease, and a highly virulent isolate, PSE9, was chosen for comparison by subtractive hybridization to a less virulent strain, PAO1. The resulting subtractive hybridization sequences were used as tags to identify genomic islands found in PSE9 but absent in PAO1. One 99-kb island, designated P. aeruginosa genomic island 5 (PAGI-5), was a hybrid of the known P. aeruginosa island PAPI-1 and novel sequences. Whereas the PAPI-1-like sequences were found in most tested isolates, the novel sequences were found only in the most virulent isolates. Deletional analysis confirmed that some of these novel sequences contributed to the highly virulent phenotype of PSE9. These results indicate that targeting highly virulent strains of P. aeruginosa may be a useful strategy for identifying pathogenicity islands and novel virulence determinants.     

Analysis of the draft genome of Pseudomonas fluorescens ATCC17400 indicates a capacity to take up iron from a wide range of sources,including different exogenous pyoverdines

All fluorescent pseudomonads (Pseudomonas aeruginosa, P. putida, P. fluorescens, P. syringae and others) are known to produce the high-affinity peptidic yellow-green fluorescent siderophore pyoverdine. These siderophores have peptide chains that are quite diverse and more than 50 pyoverdine structures have been elucidated. In the majority of the cases, a Pseudomonas species is also able to produce a second siderophore of lower affinity for iron. Pseudomonas fluorescens ATCC 17400 has been shown to produce a unique second siderophore, (thio)quinolobactin, which has an antimicrobial activity against the phytopathogenic Oomycete Pythium debaryanum. We show that this strain has the capacity to utilize 16 different pyoverdines, suggesting the presence of several ferripyoverdine receptors. Analysis of the draft genome of P. fluorescens ATCC 17400 confirmed the presence of 55 TonB-dependent receptors, the largest so far for Pseudomonas, among which 15 are predicted to be ferripyoverdine receptors (Fpv). Phylogenetic analysis revealed the presence of two different clades containing ferripyoverdine receptors, with sequences similar to the P. aeruginosa type II FpvA forming a separate cluster. Among the other receptors we confirmed the presence of the QbsI (thio)quinolobactin receptor, an ferri-achromobactin and an ornicorrugatin receptor, several catecholate and four putative heme receptors. Twenty five of the receptors genes were found to be associated with genes encoding extracytoplasmic sigma factors (ECF σ) and transmembrane anti-σ sensors.     

A Long-Chain Flavodoxin Protects Pseudomonas aeruginosa from Oxidative Stress and Host Bacterial Clearance

Long-chain flavodoxins, ubiquitous electron shuttles containing flavin mononucleotide (FMN) as prosthetic group, play an important protective role against reactive oxygen species (ROS) in various microorganisms. Pseudomonas aeruginosa is an opportunistic pathogen which frequently has to face ROS toxicity in the environment as well as within the host. We identified a single ORF, hereafter referred to as fldP (for flavodoxin from P . aeruginosa), displaying the highest similarity in length, sequence identity and predicted secondary structure with typical long-chain flavodoxins. The gene was cloned and expressed in Escherichia coli. The recombinant product (FldP) could bind FMN and exhibited flavodoxin activity in vitro. Expression of fldP in P. aeruginosa was induced by oxidative stress conditions through an OxyR-independent mechanism, and an fldP-null mutant accumulated higher intracellular ROS levels and exhibited decreased tolerance to H₂O₂ toxicity compared to wild-type siblings. The mutant phenotype could be complemented by expression of a cyanobacterial flavodoxin. Overexpression of FldP in a mutT-deficient P. aeruginosa strain decreased H₂O₂-induced cell death and the hypermutability caused by DNA oxidative damage. FldP contributed to the survival of P. aeruginosa within cultured mammalian macrophages and in infected Drosophila melanogaster, which led in turn to accelerated death of the flies. Interestingly, the fldP gene is present in some but not all P. aeruginosa strains, constituting a component of the P. aeruginosa accessory genome. It is located in a genomic island as part of a self-regulated polycistronic operon containing a suite of stress-associated genes. The collected results indicate that the fldP gene encodes a long-chain flavodoxin, which protects the cell from oxidative stress, thereby expanding the capabilities of P. aeruginosa to thrive in hostile environments.     

Isolation and Characterization of a “phiKMV-Like” Bacteriophage and Its Therapeutic Effect on Mink Hemorrhagic Pneumonia

The objective of this study was to investigate the potential of using phages as a therapy against hemorrhagic pneumonia in mink both in vitro and in vivo. Five Pseudomonas aeruginosa (P. aeruginosa) strains were isolated from lungs of mink with suspected hemorrhagic pneumonia and their identity was confirmed by morphological observation and 16S rDNA sequence analysis. Compared to P. aeruginosa strains isolated from mink with hemorrhagic pneumonia in 2002, these isolates were more resistant to antibiotics selected. A lytic phage vB_PaeP_PPA-ABTNL (PPA-ABTNL) of the Podoviridae family was isolated from hospital sewage using a P. aeruginosa isolate as host, showing broad host range against P. aeruginosa. A one-step growth curve analysis of PPA-ABTNL revealed eclipse and latent periods of 20 and 35 min, respectively, with a burst size of about 110 PFU per infected cell. Phage PPA-ABTNL significantly reduced the growth of P. aeruginosa isolates in vitro. The genome of PPA-ABTNL was 43,227 bp (62.4% G+C) containing 54 open reading frames and lacked regions encoding known virulence factors, integration-related proteins and antibiotic resistance determinants. Genome architecture analysis showed that PPA-ABTNL belonged to the “phiKMV-like Viruses” group. A repeated dose inhalational toxicity study using PPA-ABTNL crude preparation was conducted in mice and no significantly abnormal histological changes, morbidity or mortality were observed. There was no indication of any potential risk associated with using PPA-ABTNL as a therapeutic agent. The results of a curative treatment experiment demonstrated that atomization by ultrasonic treatment could efficiently deliver phage to the lungs of mink and a dose of 10 multiplicity of infection was optimal for treating mink hemorrhagic pneumonia. Our work demonstrated the potential for phage to fight P. aeruginosa involved in mink lung infections when administered by means of ultrasonic nebulization.     

Pseudomonas aeruginosa cytochrome P450 CYP168A1 is a fatty acid hydroxylase that metabolizes arachidonic acid to the vasodilator 19-HETE

Pseudomonas aeruginosa is a Gram-negative opportunistic human pathogen that is highly prevalent in individuals with cystic fibrosis (CF). A major problem in treating CF patients infected with P. aeruginosa is the development of antibiotic resistance. Therefore, the identification of novel P. aeruginosa antibiotic drug targets is of the utmost urgency. The genome of P. aeruginosa contains four putative cytochrome P450 enzymes (CYPs) of unknown function that have never before been characterized. Analogous to some of the CYPs from Mycobacterium tuberculosis, these P. aeruginosa CYPs may be important for growth and colonization of CF patients’ lungs. In this study, we cloned, expressed, and characterized CYP168A1 from P. aeruginosa and identified it as a subterminal fatty acid hydroxylase. Spectral binding data and computational modeling of substrates and inhibitors suggest that CYP168A1 has a large, expansive active site and preferentially binds long chain fatty acids and large hydrophobic inhibitors. Furthermore, metabolic experiments confirm that the enzyme is capable of hydroxylating arachidonic acid, an important inflammatory signaling molecule present in abundance in the CF lung, to 19-hydroxyeicosatetraenoic acid (19-HETE; K_m = 41 μM, V_max = 220 pmol/min/nmol P450), a potent vasodilator, which may play a role in the pathogen’s ability to colonize the lung. Additionally, we found that the in vitro metabolism of arachidonic acid is subject to substrate inhibition and is also inhibited by the presence of the antifungal agent ketoconazole. This study identifies a new metabolic pathway in this important human pathogen that may be of utility in treating P. aeruginosa infections.     

Simple Sequence Repeats Together with Mismatch Repair Deficiency Can Bias Mutagenic Pathways in Pseudomonas aeruginosa during Chronic Lung Infection

Pseudomonas aeruginosa is an opportunistic pathogen that chronically infects the airways of cystic fibrosis (CF) patients and undergoes a process of genetic adaptation based on mutagenesis. We evaluated the role of mononucleotide G:C and A:T simple sequence repeats (SSRs) in this adaptive process. An in silico survey of the genome sequences of 7 P. aeruginosa strains showed that mononucleotide G:C SSRs but not A:T SSRs were greatly under-represented in coding regions, suggesting a strong counterselection process for G:C SSRs with lengths >5 bp but not for A:T SSRs. A meta-analysis of published whole genome sequence data for a P. aeruginosa strain from a CF patient with chronic airway infection showed that G:C SSRs but not A:T SSRs were frequently mutated during the infection process through the insertion or deletion of one or more SSR subunits. The mutation tendency of G:C SSRs was length-dependent and increased exponentially as a function of SSR length. When this strain naturally became a stable Mismatch Repair System (MRS)-deficient mutator, the degree of increase of G:C SSRs mutations (5-fold) was much higher than that of other types of mutation (2.2-fold or less). Sequence analysis of several mutated genes reported for two different collections, both containing mutator and non-mutator strains of P. aeruginosa from CF chronic infections, showed that the proportion of G:C SSR mutations was significantly higher in mutators than in non-mutators, whereas no such difference was observed for A:T SSR mutations. Our findings, taken together, provide genome-scale evidences that under a MRS-deficient background, long G:C SSRs are able to stochastically bias mutagenic pathways by making the genes in which they are harbored more prone to mutation. The combination of MRS deficiency and virulence-related genes that contain long G:C SSRs is therefore a matter of concern in P. aeruginosa CF chronic infection.     

Core and accessory genome architecture in a group of Pseudomonas aeruginosa Mu-like phages

Background

Bacteriophages that infect the opportunistic pathogen Pseudomonas aeruginosa have been classified into several groups. One of them, which includes temperate phage particles with icosahedral heads and long flexible tails, bears genomes whose architecture and replication mechanism, but not their nucleotide sequences, are like those of coliphage Mu. By comparing the genomic sequences of this group of P. aeruginosa phages one could draw conclusions about their ontogeny and evolution.

Results

Two newly isolated Mu-like phages of P. aeruginosa are described and their genomes sequenced and compared with those available in the public data banks. The genome sequences of the two phages are similar to each other and to those of a group of P. aeruginosa transposable phages. Comparing twelve of these genomes revealed a common genomic architecture in the group. Each phage genome had numerous genes with homologues in all the other genomes and a set of variable genes specific for each genome. The first group, which comprised most of the genes with assigned functions, was named “core genome”, and the second group, containing mostly short ORFs without assigned functions was called “accessory genome”. Like in other phage groups, variable genes are confined to specific regions in the genome.

Conclusion

Based on the known and inferred functions for some of the variable genes of the phages analyzed here, they appear to confer selective advantages for the phage survival under particular host conditions. We speculate that phages have developed a mechanism for horizontally acquiring genes to incorporate them at specific loci in the genome that help phage adaptation to the selective pressures imposed by the host.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1146) contains supplementary material, which is available to authorized users.     

Site-specific integration of the phage ΦCTX genome into the Pseudomonas aeruginosa chromosome: characterization of the functional integrase gene located close to and upstream of attP

The site-specific integration of the phage ?CTX genome, which carries the gene for a pore-forming cytotoxin, into the Pseudomonas aeruginosa chromosome was analysed. The 1,167 by integrase gene, int, located immediately upstream of the attachment site, attP, was characterized using plasmid constructs, harbouring the integration functions, and serving as an integration probe in both P. aeruginosa and Escherichia coli. The attP plasmids p1000/p400 in the presence of the int plasmid pIBH and attP-int plasmids pINT/pINTS can be stably integrated into the P. aeruginosa chromosome. Successful recombination between the attP plasmid p1000 and the attB plasmid p5.1, in the presence of the int plasmid pIBH in E. coli HB101 showed that the int gene is active in trans in E. coli. The int gene product was detected as a 43 kDa protein in E. coli maxicells harbouring pINT. Proposed integration arm regions downstream of attP are not necessary for the integration process. pINT and phage ?CTX could be integrated together into P. aeruginosa chromosomal DNA, yielding double integrates.     

A novel Pseudomonas aeruginosa Bacteriophage,Ab31, a Chimera Formed from Temperate Phage PAJU2 and P. putida Lytic Phage AF: Characteristics and Mechanism of Bacterial Resistance

A novel temperate bacteriophage of Pseudomonas aeruginosa, phage vB_PaeP_Tr60_Ab31 (alias Ab31) is described. Its genome is composed of structural genes related to those of lytic P. putida phage AF, and regulatory genes similar to those of temperate phage PAJU2. The virion structure resembles that of phage AF and other lytic Podoviridae (S. enterica Epsilon 15 and E. coli phiv10) with similar tail spikes. Ab31 was able to infect P. aeruginosa strain PA14 and two genetically related strains called Tr60 and Tr162, out of 35 diverse strains from cystic fibrosis patients. Analysis of resistant host variants revealed different phenotypes, including induction of pigment and alginate overproduction. Whole genome sequencing of resistant variants highlighted the existence of a large deletion of 234 kbp in two strains, encompassing a cluster of genes required for the production of CupA fimbriae. Stable lysogens formed by Ab31 in strain Tr60, permitted the identification of the insertion site. During colonization of the lung in cystic fibrosis patients, P. aeruginosa adapts by modifying its genome. We suggest that bacteriophages such as Ab31 may play an important role in this adaptation by selecting for bacterial characteristics that favor persistence of bacteria in the lung.