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.     

Identification of Pseudomonas aeruginosa genes required for epithelial cell injury

We have developed a simple, reproducible and rapid genetic screen for Pseudomonas aeruginosa -induced epithelial cell cytotoxicity in cultures of MDCK cells. This screen was used to isolate isogenic transposon-tagged non-cytotoxic mutants of a cytotoxic and lung-virulent strain of P. aeruginosa (PA103). The transposon-insertion site was determined by using an inverse polymerase chain reaction followed by DNA-sequence analysis. On the basis of phenotype and sequence analysis, these mutants fell into four classes. One class had absent or defective pili, based on their resistance to phage PO4 and/or loss of twitching motility (twt⁻). A second class exhibited decreased adherence. A third class of mutants exhibited probable defects in the machinery or targets of type III protein secretion. A final class of mutants exhibited decreased but not absent cytotoxicity. This class included members of the first three classes as well as other mutants. These results suggest that localized cytotoxicity is likely to require several steps and several components, including pili and other (unidentified) extracellular proteins. The type III protein-secretion apparatus appears to be involved in this process.     

铜绿假单胞菌泳动能力相关新基因的筛选及鉴定

从Mu转座突变子文库中经过表型筛选,得到12株泳动(Swimming motility)能力缺陷的突变子,经Mu转座子插入位点的确认、基因克隆及测序分析发现其中10个突变子中Mu转座子分别插入到10个不同的与鞭毛运动和功能相关的基因中,2个突变子中Mu转座子插入到功能未知的新基因(PA2950和PA5022)中,电镜观察结果表明这2个突变株均具有完整的鞭毛,初步推测这2个基因可能是参与鞭毛泳动的能量代谢、趋化作用或信息传递的新基因。     

Identification of in vivo essential genes from Pseudomonas aeruginosa by PCR-based signature-tagged mutagenesis

We adapted PCR-based signature-tagged mutagenesis (STM) to Pseudomonas aeruginosa. A collection of 1056 mutants was screened in a chronic lung infection rat model. Thirteen mutants were confirmed to be attenuated. Analysis revealed that these STM mutants represented transposon insertions into eight genes previously described in databases, three genes encoding proteins sharing identity with hypothetical proteins and two genes that shared no significant identity with sequences in databases. Five strains mutated in genes involved in protein degradation, stress tolerance, cation transport, ABC transporter, and an unknown protein were shown to be highly attenuated when tested individually in the rat chronic lung infection model.     

Identification of RegA protein from Pseudomonas aeruginosa using anti-RegA antibody

The product of Pseudomonas aeruginosa regA gene acts as a positive regulator of exotoxin A expression. The protein, RegA, was overproduced in E. coli transformed with an expression vector containing the regA gene. The overproduced RegA accumulated in E. coli in the form of inclusion bodies. The latter were isolated and served as a source of antigen for raising polyclonal antibodies. The antibodies reacted specifically with a P. aeruginosa protein whose molecular weight corresponded to that predicted for RegA from its known DNA sequence, and whose response to modulating factors matched that expected for RegA. The immunodetectable RegA was localized in the membrane fraction of P. aeruginosa strain PA103.     

铜绿假单胞菌必需基因的研究进展

铜绿假单胞菌为专性需氧非发酵革兰氏阴性杆菌,是医院感染的常见条件致病菌之一,可引起呼吸道、泌尿道、烧伤创面和菌血症等严重感染.铜绿假单胞菌耐药形势日益严峻,给临床治疗带来困难.必需基因是生长过程中必不可少的看家基因,对铜绿假单胞菌必需基因进行深入研究,不仅有助于了解细菌的生长、毒力等基本特性,也有助于筛选新的抗菌药物靶...     

Identification of genes involved in swarming motility using a Pseudomonas aeruginosa PAO1 mini-Tn5-lux mutant library

During a screening of a mini-Tn5-luxCDABE transposon mutant library of Pseudomonas aeruginosa PAO1 for alterations in swarming motility, 36 mutants were identified with Tn5 insertions in genes for the synthesis or function of flagellin and type IV pilus, in genes for the Xcp-related type II secretion system, and in regulatory, metabolic, chemosensory, and hypothetical genes with unknown functions. These mutants were differentially affected in swimming and twitching motility but in most cases had only a minor additional motility defect. Our data provide evidence that swarming is a more complex type of motility, since it is influenced by a large number of different genes in P. aeruginosa. Conversely, many of the swarming-negative mutants also showed an impairment in biofilm formation, indicating a strong relationship between these types of growth states.     

Cloning and characterization of chemotaxis genes in Pseudomonas aeruginosa

Two chemotaxis-defective mutants of Pseudomonas aeruginosa, designated PC3 and PC4, were selected by the swarm plate method after N-methyl-N'-nitro-N-nitrosoguanidine mutagenesis. These mutants were not complemented by the P. aeruginosa cheY and cheZ genes, which had been previously cloned (Masduki et al., J. Bacteriol., 177, 948-952, 1995). DNA sequences downstream of the cheY and cheZ genes were able to complement PC3 but not PC4. Sequence analysis of a 9.7-kb region directly downstream of the cheZ gene found three chemotaxis genes, cheA, cheB, and cheW, and seven unknown open reading frames (ORFs). The predicted translation products of the cheA, cheB, and cheW genes showed 33, 36, and 31% amino acid identity with Escherichia coli CheA, CheB, and CheW, respectively. Two of the unknown ORFs, ORF1 and ORF2, encoded putative polypeptides that resembled Bacillus subtilis MotA (40% amino acid identity) and MotB (34% amino acid identity) proteins, respectively. Although P. aeruginosa was found to have proteins similar to the enteric chemotaxis proteins CheA, CheB, CheW, CheY, and CheZ, the gene encoding a CheR homologue did not reside in the chemotaxis gene cluster. The P. aeruginosa cheR gene could be cloned by phenotypic complementation of the PC4 mutant. This gene was located at least 1,800 kb away from the chemotaxis gene cluster and encoded a putative polypeptide that had 32% amino acid identity with E. coli CheR.     

耐亚胺培南铜绿假单胞菌相关耐药基因的检测

目的对耐亚胺培南（IMP）的铜绿假单胞菌（IRPa）相关耐药基因进行检测。方法 2003年至2009年从临床标本中分离到（P.aeruginosa）共220株,采用三维试验筛选产β-内酰胺酶的铜绿假单胞菌,应用普通PCR和多重PCR分别检测碳青霉烯酶基因和质粒携带的C类头孢菌素酶（AmpC酶）耐药基因,应用荧光定量RT-PCR的方法检测oprD2基因的表达情况。结果共检出43株产β-内酰胺酶的菌株,其中产AmpC酶、超广谱β-内酰胺酶（ESBLs）、金属β-内酰胺酶（MBLs）和未知酶菌株的构成比分别58.14%（25/43）、18.60%（8/43）、4.65%（2/43）和16.28%（7/43）。74株耐亚胺培南的铜绿假单胞菌中,有2株菌携带IMP-9基因,1株菌携带DHA质粒型AmpC酶基因,其他碳青霉烯酶基因检测为阴性。40株菌株oprD2基因表达蛋白量降低,34株oprD2基因表达蛋白量正常。结论 oprD2基因的突变或蛋白表达量降低是IRPa对亚胺培南耐药的主要原因,AmpC酶可水解亚胺培南可能与铜绿假单胞菌对亚胺培南的耐药有一定的关系,而KPC-1酶和MBLs在铜绿假单胞菌对亚胺培南耐药机制中不是主要因素。     

Transductional analysis of the flagellar genes in Pseudomonas aeruginosa.

Complementation in bacteriophage E79 tv-l-mediated transduction and the phenotypic properties of the flagellar genes in Pseudomonas aeruginosa PAO were investigated by using 195 flagellar mutants of this organism. A total of 15 fla. 1 mot, and 2 che cistrons were identified. At least 5 fla cistrons (fla V to flaZ) and one mot cistron resided in one region, and at least 10 fla cistrons (flaA to flaJ) and two che cistrons (cheA and cheB) resided in another. The flaC mutants exhibited cistron-specific leakiness on motility agar plates. The flaE cistron may be the structural gene for the component protein of the flagellar filament. The cheA mutations, which resulted in pleiotropic phenotypes for flagellar formation, motility, and taxis, belonged to the same complementation group as the flaF mutations; that is, we inferred that cheA and flaF are synonymous.     

Identification and characterization of porins in Pseudomonas aeruginosa.

Earlier studies have shown that the major porin species in Pseudomonas aeruginosa outer membrane is protein F (OprF), which produces channels wider than those produced by Escherichia coli porins. In contrast, Yoshihara and Nakae ((1989) J. Biol. Chem. 264, 6297-6301) reported that protein F has no pore-forming activity as measured by the flux of L-arabinose, and that the channels in P. aeruginosa outer membrane, being produced by proteins C, "D," and "E," are much narrower than E. coli porin channels. In this study, we followed the protein purification scheme of Yoshihara and Nakae as closely as possible, and found that protein F had a specific activity for pore formation similar to that of proteins D1, D2, and E2. Furthermore, proteoliposome reconstitution assays showed conclusively that the channels formed by protein F, as well as by unfractionated outer membranes, allowed the diffusion of a tetrasaccharide, stachyose, at a significant rate, indicating that these channels are much larger than E. coli porin channels. It appears likely that in the study of Yoshihara and Nakae protein F was inadvertently inactivated during purification. We further suggest a hypothesis that resolves the apparent conflict between the presence of large diameter channels and the low permeability of the outer membrane in P. aeruginosa.     

Cloning of genes specifying carbohydrate catabolism in Pseudomonas aeruginosa and Pseudomonas putida.

A 6.0-kilobase EcoRI fragment of the Pseudomonas aeruginosa PAO chromosome containing a cluster of genes specifying carbohydrate catabolism was cloned into the multicopy plasmid pRO1769. The vector contains a unique EcoRI site for cloning within a streptomycin resistance determinant and a selectable gene encoding gentamicin resistance. Mutants of P. aeruginosa PAO transformed with the chimeric plasmid pRO1816 regained the ability to grow on glucose, and the following deficiencies in enzyme or transport activities corresponding to the specific mutations were complemented: glcT1, glucose transport and periplasmic glucose-binding protein; glcK1, glucokinase; and edd-1, 6-phosphogluconate dehydratase. Two other carbohydrate catabolic markers that are cotransducible with glcT1 and edd-1 were not complemented by plasmid pRO1816: zwf-1, glucose-6-phosphate dehydrogenase; and eda-9001, 2-keto-3-deoxy-6-phosphogluconate aldolase. However, all five of these normally inducible activities were expressed at markedly elevated basal levels when transformed cells of prototrophic strain PAO1 were grown without carbohydrate inducer. Vector plasmid pRO1769 had no effect on the expression of these activities in transformed mutant or wild-type cells. Thus, the chromosomal insert in pRO1816 contains the edd and glcK structural genes, at least one gene (glcT) that is essential for expression of the glucose active transport system, and other loci that regulate the expression of the five clustered carbohydrate catabolic genes. The insert in pRO1816 also complemented the edd-1 mutation in a glucose-negative Pseudomonas putida mutant but not the eda-1 defect in another mutant. Moreover, pRO1816 caused the expression of high specific activities of glucokinase, an enzyme that is naturally lacking in these strains of Pseudomonas putida.     

Phenotypic mixing of pyocin R2 and bacteriophage PS17 in Pseudomonas aeruginosa PAO.

Previous results indicate that a group of bacteriocins in Pseudomonas aeruginosa, named R-type pyocins, have a structure resembling bacteriophage tails and share some serological homology with certain bacteriophages. This paper presents genetic evidence which strongly suggests that components of pyocin R2, an R-type pyocin of P. aeruginosa PAO, and tail components of bacteriophage PS17 are interchangeable. Complementation tests with pyocin R2-deficient mutants of PAO and ts mutants of PS17 revealed that various phenotypic interactions occur between the pyocin and bacteriophage in PAO cells lysogenized or infected with PS17. (i) Certain pyocin R2-deficient mutations were phenotypically suppressed in cells carrying PS17 prophage. (ii) A temperature-sensitive mutant of PS17, tsQ31, was phenotypically suppressed in PAO cells treated with mitomycin C. (iii) Phenotypically mixed phages with receptor and serological specificities of pyocin R2 were formed in PS17 lysogens of certain pyocin R2-deficient mutants.     

Identification of genes and proteins necessary for catabolism of acyclic terpenes and leucine/isovalerate in Pseudomonas aeruginosa

Geranyl-coenzyme A (CoA)-carboxylase (GCase; AtuC/AtuF) and methylcrotonyl-CoA-carboxylase (MCase; LiuB/LiuD) are characteristic enzymes of the catabolic pathway of acyclic terpenes (citronellol and geraniol) and of saturated methyl-branched compounds, such as leucine or isovalerate, respectively. Proteins encoded by two gene clusters (atuABCDEFGH and liuRABCDE) of Pseudomonas aeruginosa PAO1 were essential for acyclic terpene utilization (Atu) and for leucine and isovalerate utilization (Liu), respectively, as revealed by phenotype analysis of 10 insertion mutants, two-dimensional gel electrophoresis, determination of GCase and MCase activities, and Western blot analysis of wild-type and mutant strains. Analysis of the genome sequences of other pseudomonads (P. putida KT2440 and P. fluorescens Pf-5) revealed candidate genes for Liu proteins for both species and candidate genes for Atu proteins in P. fluorescens. This result concurred with the finding that P. fluorescens, but not P. putida, could grow on acyclic terpenes (citronellol and citronellate), while both species were able to utilize leucine and isovalerate. A regulatory gene, atuR, was identified upstream of atuABCDEFGH and negatively regulated expression of the atu gene cluster.     

Identification of the Major Ubiquitin-binding Domain of the Pseudomonas aeruginosa ExoU A2 Phospholipase

Numerous Gram-negative bacterial pathogens use type III secretion systems to deliver effector molecules into the cytoplasm of a host cell. Many of these effectors have evolved to manipulate the host ubiquitin system to alter host cell physiology or the location, stability, or function of the effector itself. ExoU is a potent A₂ phospholipase used by Pseudomonas aeruginosa to destroy membranes of infected cells. The enzyme is held in an inactive state inside of the bacterium due to the absence of a required eukaryotic activator, which was recently identified as ubiquitin. This study sought to identify the region of ExoU required to mediate this interaction and determine the properties of ubiquitin important for binding, ExoU activation, or both. Biochemical and biophysical approaches were used to map the ubiquitin-binding domain to a C-terminal four-helix bundle of ExoU. The hydrophobic patch of ubiquitin is required for full binding affinity and activation. Binding and activation were uncoupled by introducing an L8R substitution in ubiquitin. Purified L8R demonstrated a parental binding phenotype to ExoU but did not activate the phospholipase in vitro. Utilizing these new biochemical data and intermolecular distance measurements by double electron-electron resonance, we propose a model for an ExoU-monoubiquitin complex.