- Open Access
Genomic analysis reveals that Pseudomonas aeruginosa virulence is combinatorial
Genome Biology volume 7, Article number: R90 (2006)
Pseudomonas aeruginosa is a ubiquitous environmental bacterium and an important opportunistic human pathogen. Generally, the acquisition of genes in the form of pathogenicity islands distinguishes pathogenic isolates from nonpathogens. We therefore sequenced a highly virulent strain of P. aeruginosa, PA14, and compared it with a previously sequenced (and less pathogenic) strain, PAO1, to identify novel virulence genes.
The PA14 and PAO1 genomes are remarkably similar, although PA14 has a slightly larger genome (6.5 megabses [Mb]) than does PAO1 (6.3 Mb). We identified 58 PA14 gene clusters that are absent in PAO1 to determine which of these genes, if any, contribute to its enhanced virulence in a Caenorhabditis elegans pathogenicity model. First, we tested 18 additional diverse strains in the C. elegans model and observed a wide range of pathogenic potential; however, genotyping these strains using a custom microarray showed that the presence of PA14 genes that are absent in PAO1 did not correlate with the virulence of these strains. Second, we utilized a full-genome nonredundant mutant library of PA14 to identify five genes (absent in PAO1) required for C. elegans killing. Surprisingly, although these five genes are present in many other P. aeruginosa strains, they do not correlate with virulence in C. elegans.
Genes required for pathogenicity in one strain of P. aeruginosa are neither required for nor predictive of virulence in other strains. We therefore propose that virulence in this organism is both multifactorial and combinatorial, the result of a pool of pathogenicity-related genes that interact in various combinations in different genetic backgrounds.
The potential virulence of bacterial pathogens is significantly modulated by the presence of pathogenicity islands [1, 2], which are clusters of one or more virulence-related genes that are often acquired by horizontal gene transfer. The introduction of these virulence islands can allow a previously nonvirulent isolate to infect a particular host. Commonly, this switch to a simpler and more stable environment within a host (as opposed to the more complex outside environment) is followed by gene loss and genome reduction that improve the ability of the pathogen to survive in the host but also restrict the range of hosts available to the bacterium [3, 4]. In contrast, free-living bacteria that dominate in complex environments (such as soil) tend to have genomes that continue to acquire DNA and undergo expansion rather than reduction.
Pseudomonas aeruginosa, a ubiquitous Gram-negative soil organism, is an important opportunistic human pathogen that infects injured, burned, immunodeficient, and immunocompromised patients, and causes persistent respiratory infections in individuals suffering from cystic fibrosis (CF) [5, 6]. The genome sequence of the widely studied P. aeruginosa strain PAO1 (originally a wound isolate) revealed that it possesses a large number of genes that are involved in regulation, catabolism, transport, and efflux of organic compounds, as well as several putative chemotaxis systems , all of which potentially contribute to the remarkable ability of this bacterium to adapt to a wide range of environmental niches.
Different P. aeruginosa isolates share a remarkable amount of similarity in their genomes. When DNA derived from several P. aeruginosa strains was hybridized to a PAO1 microarray, between 89% and 98% of the PAO1 sequences were detected [8, 9]. Whole-genome shotgun sequencing of two CF isolates and one environmental strain revealed that, aside from this apparent highly conserved core set of P. aeruginosa genes, differences were largely due to strain-specific islands of genes, consisting either of genes with similar or related function but divergent DNA sequences (such as genes for biosynthesis of the O-antigen component of lipopolysaccharide, genes for flagellar biosynthesis, or alternate forms of genes for the bacteriocidal pyocins) or genes that are entirely absent in some strains .
Despite the overall genome similarity among diverse P. aeruginosa strains, differences in complex phenotypes such as pathogenicity can be striking. For example, the clinical isolate PA14 is significantly more virulent than PAO1 in a wide range of hosts, including mice, the nematode Caenorhabditis elegans, the insect Galleria mellonella, and the plant Arabidopsis thaliana [11–13]. PA14 genes required for full virulence include genes common to many if not all P. aeruginosa strains, including global transcriptional regulators such as gacA; genes that are involved in pathogenesis-related processes such as motility, quorum sensing, and phenazine biosynthesis; and genes that encode secreted cellulytic factors and toxins such as ExoU, exotoxin A, phospholipase C, and elastase [12–18]. On the other hand, novel PA14 genes that are absent in PAO1 (and potentially absent in other isolates) have also been identified as being required for pathogenicity in model hosts and mice [11, 15, 17, 18], and at least some of these genes appear to reside on large pathogenicity islands . Taken together, these studies suggest that PA14 pathogenicity is multifactorial, requiring the cumulative (and potentially coordinated) action of multiple virulence factors, some of which are components of the basic core genome, whereas others are located on classically defined virulence islands.
The experiments described in this paper were designed to test the hypothesis that the enhanced virulence of PA14 compared with PAO1 is mostly a consequence of recognizable pathogenicity (virulence) islands that are present in PA14 but absent in PAO1. To expand our tools for dissecting P. aeruginosa virulence and to test the hypothesis that strain differences in virulence are due to the acquisition of strain-specific genes, we sequenced the PA14 genome and performed a functional analysis of genes that are present in PA14 but absent in PAO1 to assess their contribution to pathogenicity.
Results and discussion
PA14 genome sequence
To identify all of the putative pathogenicity islands that distinguish PA14 from PAO1, we sequenced the PA14 genome and found that it contains a slightly larger chromosome (6.5 megabases [Mb] versus 6.3 Mb for PAO1; genome sequence and annotations are available at the Ausubel lab PA14 sequencing website  and have also been deposited in GenBank [GenBank: CP000438]). Consistent with previous observations that overall strain similarity is high, we found that approximately 91.7% of the PA14 genome is present in PAO1, and that 95.8% of the PAO1 genome is present in PA14 (Table 1).
The PA14 and PAO1 genomes are largely colinear with the exception of one major rearrangement that has previously been described: an inversion between two of four dispersed copies of a large ribosomal RNA cluster (Figure 1). We used long-range polymerase chain reaction (PCR) spanning the rRNA clusters to verify the inversion in the sequenced PA14 and PAO1 genomes; PCR products predicted to be indicative of one orientation or the other were only generated when either PAO1 or PA14 DNA was used as the template (Figure 2a,b). When 18 diverse P. aeruginosa strains were subjected to this same PCR analysis, all 18 of these isolates were the same or resembled PA14 genome structure with respect to the large inversion, whereas none of the 18 strains had or resembled the PAO1 genome structure (data not shown). Of note, additional PAO1 clones (distinct from the isolate that was sequenced) also appear to contain the same inversion found in PA14 .
We also compared the GC skews of the PA14 and PAO1 genomes. GC skew is defined as the value of [G-C]/[G+C] where G and C represent the local base frequencies of G and C, respectively. In prokaryotic genomes there is a bias toward G over C on the leading strand of DNA synthesis . Therefore, when measured at regular intervals along the chromosome (1 kilobase [kb] segments) GC skew tends to have a positive value on the leading strand of DNA synthesis and a negative value on the lagging strand, resulting in polarity changes at the origin and terminus of replication. The putative position of the replicative terminus is therefore operationally defined as the peak of the cumulative GC skew and it typically resides opposite the origin of replication in bacterial genomes . For PA14, the peak GC skew was indeed mapped opposite the replication origin (at 49.2% of the genome; Figure 2a,c). In contrast, the position of the terminus in the sequenced PAO1 chromosome is shifted relative to the origin (at 38.8% of the genome). These observations argue that the PA14 chromosome with respect to this large inversion is more representative of the canonical or ancestral P. aeruginosa genome than that of the sequenced PAO1 isolate. This interpretation that is further supported by the finding that the inversion in the sequenced strain occurred sometime after 1990 (a PAO1 cosmid library originally described in that year was shown not to contain the inversion ). However, the physiologic consequences of such an inversion and asymmetric replication cycle are not clear.
PA14 gene annotation
We annotated the PA14 genome and identified 5973 predicted open reading frames (ORFs; 322 more than in PAO1; Table 1). Summaries of all predicted PA14 genes and their distribution among different functional categories are presented in Additional data files 2 and 3, respectively. We performed a reciprocal NCBI-BLAST (basic local alignment search tool) search of each gene in either PA14 or PAO1 against the total collection of genes in the other strain (specifically, each PA14 gene was BLASTed against the complete set of PAO1 genes, and vice versa). The position of each best BLAST hit was plotted (with the PA14 genomic location on the x-axis, and the PAO1 coordinate on the y-axis; Figure 3), confirming that the two genomes are largely colinear, with the exception of the large inversion described above. Individual PA14 genes that had no counterpart in PAO1 are shown as individual (green) data points on the x-axis, and PAO1 genes absent in PA14 are represented as (pink) data points on the y-axis. Large clusters of these strain-specific genes often correlated with regions of the respective genomes in which the local GC content was lower than that of the total genome (Figure 3).
Using a combination of raw sequence and ORF-based global alignments of the two genomes, we compiled a list of gene clusters present in one strain but absent in the other (58 PA14 regions absent in PAO1 containing 478 genes, and 54 PAO1 regions absent in PA14 containing 234 genes; Table 1 and Additional data file 4). We refer to these as PA14-specific or PAO1-specific regions for the purposes of this discussion (recognizing that these genes may be present in other isolates of P. aeruginosa and are not strictly strain specific). Many of these gene clusters have hallmarks of horizontally acquired DNA, including direct repeats, insertion sequences, tRNA genes at their boundaries (data not shown), and/or anomalous GC contents (PA14-specific clusters have an average GC content of 59.6%, which is more than one standard deviation below the whole genome average of 66.3%; Table 1).
These strain-specific regions contain a high proportion of genes of unknown function as compared with the whole genome. Each PA14 (and PAO1) annotation has a confidence score associated with the described gene function [7, 22]. A class 1 designation refers to genes whose functions have been experimentally validated in P. aeruginosa, class 2 genes are highly similar to genes whose functions have been validated in other organisms, class 3 genes have hypothesized functions based on limited similarity to other genes or structural/functional domains, and class 4 genes are ORFs of unknown function. For the PA14-specific regions, 63% of the predicted ORFs have no known function, whereas only 28.9% of genes in the whole genome have class 4 annotations (Table 1). Therefore, gene identity alone could not indicate whether virulence-related genes were enriched in PA14-specific regions.
Conservation of PA14-specific genes and their putative role in virulence
If a PA14 strain-specific region were functioning as a canonical pathogenicity island, then it should be more prevalent among other pathogenic isolates (and less prevalent among avirulent isolates). To test this hypothesis we needed to establish an objective measure of pathogenicity to compare different P. aeruginosa strains as well as to develop a high throughput method for determining the genomic content of the different strains. We used a model host infection system (the nematode Caenorhabditis elegans) to rank order the virulence of 18 diverse P. aeruginosa strains, using PA14 and PAO1 as reference strains. The 18 P. aeruginosa strains included 13 clinical isolates from a variety of infection types (CF lung infections, urinary tract infections, ocular infections, and blood isolates), one laboratory strain, and four environmental isolates. This same set of strains had been used in a previous study of P. aeruginosa strain diversity . We used a custom microarray-based system (described below) to determine relatedness of the 20 (18 plus PAO1 and PA14) P. aeruginosa strains.
In the C. elegans model pathogenicity system , an age-synchronized population of nematodes is fed the pathogen to be tested (instead of Escherichia coli, its traditional laboratory food source) and the longevity of the nematodes is determined. As shown previously, the longevity of C. elegans feeding on a particular P. aeruginosa strain correlates with virulence of the strain in mice [13, 15, 18]. The set of P. aeruginosa strains tested exhibited a full range of virulence phenotypes, including both the upper and lower limits that the assay system is capable of measuring (Figure 4a). PA14 (dark blue diamonds, second curve from the left) is extremely efficient at killing nematodes, whereas the less pathogenic PAO1 is intermediate (pink squares), killing more slowly than PA14 but more quickly than E. coli (negative control; yellow squares, second curve from the right).
The observation that PAO1 is more virulent than many other tested strains suggests that it has not become attenuated because of extensive passaging in the laboratory; rather, its virulence is probably that of a moderately pathogenic strain, with other strains more representative of truly 'avirulent' isolates. When comparing strains derived from the same type of infection, there was no consistent clustering with respect to their phenotype in C. elegans. For example, both the most and the least virulent strains tested were isolates from CF infections. Similarly, the five urinary tract infection strains exhibited a wide range of virulent to avirulent phenotypes. Importantly, two closely related environmental isolates (MSH3 and MSH10) were the fourth and fifth most virulent strains tested, indicating that nonclinically isolated strains can also have the potential to be infectious.
To test the hypothesis that the virulence of the 20 P. aeruginosa strains correlates with the presence of particular virulence islands, we performed a microarray-based analysis of the genomic content (a process described as genomotyping ) of the P. aeruginosa strains. We arrayed 285 synthetic oligonucleotides (70 mers) corresponding to PA14 genes that are absent in PAO1 and 130 oligonucleotides corresponding to PAO1 genes that are absent in PA14, along with additional sequences serving as positive and negative controls (see Materials and methods, below, and Additional data files 5 and 6). We first used the array data to generate a hierarchic clustering dendrogram showing the relatedness of the 20 P. aeruginosa strains (Figure 4b). Next, we used the C. elegans virulence data in Figure 4a to rank order the 20 strains and then examine the microarray data to determine whether the presence or absence of PA14-specific (or PAO1-specific) genes correlated with virulence (also see Table 2 and Additional data file 6). Strains were ordered from the most virulent (rank order 1) to the least virulent (rank order 20), with strains that had equivalent virulence assigned a tied rank. When the virulence rank order of each strain was superimposed on the dendrogram describing the relatedness of each strain based on genomic content (Figure 4b), no strong correlation was observed with the relative virulence of each strain in C. elegans. There are small clusters of similarly avirulent or virulent strains; for example, strains 6077, U2504, JJ692, S54485, and X13273 grouped together and possess similar virulence rankings (10, 6, 12, 10, and 10, respectively). However, these small clusters were occasionally punctuated by exceptions: a cluster of weakly virulent or avirulent strains (CF5, E2, PAK and CF27, with rank orders of 20, 17, 13, and 18, respectively) also includes strain UDL, which is the third most pathogenic strain tested.
In a previous study, Wolfgang and coworkers  used Affymetrix GeneChips to survey these same 18 strains (with the exception of PA14) for the presence or absence of PAO1 genes and found no obvious pattern that would correlate genomic content with disease (or the type of infection from which the isolate was derived). Similarly, examination of the distribution of PA14-specific genes among the tested isolates did not reveal any obvious clustering of genomic content with respect to the source of the strain (see dendrogram in Figure 4b).
Identification of PA14-specific virulence genes and their conservation in other strains
The experiments described above showed that there was no correlation between the PA14-specific sequences in general and virulence in the C. elegans killing assay. We therefore performed a functional analysis of PA14-specific ORFs, identifying genes that were specifically required for pathogenicity and subsequently assessing their distribution among the other strains. Our laboratory has constructed a genome-wide, nonredundant transposon insertion mutant library in PA14 [24, 25]. Using this library, we conducted a screen for mutants in PA14-specific genes that had reduced virulence in C. elegans. Nine genes were identified, which fell into six distinct clusters (Table 3). Three of these clusters are present in all of the other 19 P. aeruginosa strains but contain highly divergent sequences, including the O-antigen biosynthetic cluster (region PA14R38) and two groups of type 4 fimbrial biogenesis genes (regions PA14R77 and PA14R79). Oligonucleotides corresponding to the PA14-version of the O-antigen biosynthesis genes were included on the array and demonstrated that none of the other strains contain the PA14-versions of these genes; the genes for the two fimbrial synthesis clusters were not included on the array (because the sequences were not sufficiently divergent to meet our criteria for oligonucleotide design as outlined in Additional data file 1).
The remaining three clusters contain genes that are absent in PAO1. As indicated in Table 3, these three clusters contain five virulence genes, including one gene with putative transcriptional regulator activity and four ORFs of unknown function. Figure 4c summarizes the array data for these three regions. Each of the 20 strains tested are represented in columns arranged from left to right in order of decreasing pathogenicity (see Table 2 for rank order of virulence). The PA14 genes tested (shown in rows) are described as present (blue), absent (yellow), or indeterminate (red, indicating that the hybridization intensity was intermediate and a present or absent call could not be made with confidence; see Additional data file 1 for further details). The positions of the genes that result in an avirulent phenotype when mutated in PA14 are indicated on the right (also see Table 3). If a cluster of genes was required for or predictive of virulence in nematodes then we would expect to see a bias for present (blue) calls toward the left and absent (yellow) calls toward the right.
For these five virulence genes (and their associated three clusters), Spearman's rank correlation coefficients (relating their presence, absence, or indeterminate status with the virulence of the strain) were found to vary dramatically. The best correlation was found for two genes in region PA14R41 (correlation coefficients of 0.63). This is a cluster of eight genes known as the clone C-specific region common to clone C isolates (members of a clone family associated with CF infections ). Intermediate to no correlations were observed for the remaining two regions, namely PA14R78 and PAR09, which had correlation coefficients of 0.44 and 0.02, respectively. Of note, PA14R78 is a previously described pathogenicity island (PAPI-1) shown to contain genes required for pathogenicity in plants and mammals , although the two mutations in genes of unknown function that we identified in this cluster (PA14_59010 and PA14_59070; positions 6 and 7 in Figure 4c) were not among those previously examined. Regardless of the magnitude of the correlation, each virulence gene had exceptions to the expected trend, being present in attenuated strains and/or absent in virulent strains. Taken together, the functionally defined PA14 genes required for C. elegans killing are neither required for nor necessarily predictive of another strain's ability to be pathogenic.
In this study we combined a traditional comparative genomic analysis with functional analyses, including genomotyping of other isolates, use of a genome-wide mutant library, and a model host infection system amenable to high-throughput screens, to ascertain the relationship between genomic content and virulence. We have found that PA14 and PAO1 are remarkably similar in total content and are largely colinear (with the exception of a previously described inversion). Because PA14 is significantly more virulent than PAO1 in most hosts, and previous studies have identified PA14 pathogenicity genes in regions that are absent in PAO1, we extended these observations to examine all of the PA14-specific genes for potential contributions to virulence. We began by examining 18 additional diverse isolates to assess the conservation of PA14-specific genes in other strains and whether they correlated with pathogenicity. In general, there was no obvious relationship between the presence of PA14-specific or PAO1-specific genes in other isolates and their virulence in C. elegans or the source of the strain. A specific examination of genes in these PA14-specific regions that were experimentally shown to contribute to virulence also failed to show a strong correlation between the presence of these genes and the pathogenicity of other isolates.
Our results amend the general view of pathogenicity islands, in which the acquisition of an island leads to the addition of gene products whose contribution to virulence is apparent based on the known function of the virulence gene and does not depend (necessarily) on genes outside of the island. For example, an island may contain a complete complement of genes required for the synthesis of secreted toxins, adhesins, invasion systems, iron uptake systems, or secretion systems (such as the type III and type IV systems) [1, 2]. In these 'classic' cases, the presence or absence of individual pathogenicity islands correlates directly with a given gene product or process known to be required for virulence (and, therefore, with overall virulence). P. aeruginosa strains contain genes that fall into this category; for example, the presence of a pathogenicity island containing the type III secretion effector ExoU (present on PA14 region PA14R72) makes strains more cytotoxic to mammalian cells  and is required for pathogenicity in Galleria and the amoeba Dictyostelium discoideum [16, 28]. However, we have also identified PA14 pathogenicity-related genes whose presence or absence does not correlate directly with degree of virulence, suggesting that these genes do not function autonomously to affect virulence. Our genomic analysis of PA14 virulence has demonstrated that pathogenicity in this organism is both multifactorial and combinatorial. Within a given isolate, virulence is multifactorial in that several factors combine to result in an overall virulence phenotype. Additionally, when comparing different strains, virulence is combinatorial in that pathogenicity factors may behave differently and that distinct combinations or groupings of these determinants may result in comparable virulence phenotypes.
What might account for the apparent complexity of PA14-specific virulence factors with respect to their conservation and role in other strains? First, it is clear that our analysis of genomic content in other strains is only an initial step in addressing similarities or differences in gene function among isolates; the apparent absence of a gene cannot exclude the existence of a functionally similar gene with significant sequence divergence, and the apparent presence of a gene cannot determine whether a gene is expressed. Transcriptional profiling of strain-specific genes (particularly in the presence of a putative host) will be a critical next step in clarifying the genes that contribute to virulence. Additional sequence information from other strains will be required to determine whether other functionally similar genes exist or whether orthologous genes are likely to be fully functional (intact ORFs with no polymorphisms that might alter protein function). Second, given the multifactorial and combinatorial nature of P. aeruginosa virulence, a full understanding of pathogenicity will require elucidation of how strain-specific genes (potentially responsible for differences in severity of disease among isolates) interact with core genome genes required for a base level of virulence. Our laboratory is currently extending the screen for PA14 mutants attenuated in C. elegans killing beyond the PA14-specific genes to include the entire genome (Liberati NT, Feinbaum RL, Ausubel FM, unpublished observations). Third, an important future direction will be to determine how generalizable our observations are in other model hosts. We are currently conducting pilot experiments to assess the viability of screening both the set of 20 P. aeruginosa strains and the PA14 mutant library in insects and plants to determine whether the presence or absence of identified PA14 virulence genes in other isolates correlates with their overall pathogenicity. A thorough understanding of networks of genes that are necessary for virulence in many or the majority of P. aeruginosa isolates (as opposed to genes that contribute to pathogenicity only in a subset of strains in which they are present) will be crucial for the design of effective therapeutics to combat the wide variety of human infections observed in clinical settings.
The evolution of virulence for a dedicated human pathogen generally involves the acquisition of discrete virulence functions required for specific interactions with the host, followed by gene loss related to specialization and potential restriction to the new environmental niche. In contrast, ubiquitous environmental micro-organisms continuously encounter dramatic changes in their ecosystem and the maintenance of genome complexity is preferable to optimization for a single niche. For environmental pathogenic fungi such as Cryptococcus neoformans that do not require animal hosts for replication or survival, the phenomenon of 'ready made' virulence has been described in which the selection and maintenance of virulence factors occurs during infection of environmental predators such as nematodes and ameba [29–31]. We propose that this multi-host view of pathogenic evolution also applies to environmental bacteria such as P. aeruginosa that infects nematodes, insects, plants, and ameba in the laboratory and probably encounter a similar range of potential hosts in the wild. The recent discovery that C. elegans pre-exposed to PA14 can modify its olfactory preferences to avoid the pathogen, and that this 'learned' avoidance behavior does not occur with nonpathogenic mutants of PA14, suggests that an interaction between these two organisms is biologically relevant . Therefore, selection for pathogenicity may occur constantly in the environment; indeed, among the 20 isolates surveyed in this study, two environmental strains (MSH3 and MSH10) were among the most virulent, demonstrating that pathogenic potential exists even in nonclinical isolates. Furthermore, as an environmental organism, P. aeruginosa is not likely to undergo genome reduction associated with host restriction as in dedicated human pathogens. A notable exception is the case of CF infections, in which a clonal population remains isolated in a defined environment over a long period of time and tends not to spread to other patients. Indeed, a recent comparison of an early and late CF isolate from the same patient (at 6 months and 8 years of age) revealed an accumulation of multiple mutations leading to gene inactivation or loss, including a large 188 kb deletion . When these mutated genes were examined in early and late infections of 29 additional CF patients, most of the observed mutations arose relatively late in the infections.
Our laboratory and others have previously shown that aspects of PA14 pathogenicity are conserved in model hosts and in mammalian hosts; PA14 virulence factors have been shown to be required for disease in both model hosts and mammals [13, 14, 16, 17, 19, 28, 34], and components of the host defense response are conserved in model hosts and humans [34–36]. Although these observations are critical in validating the use of simple, genetically tractable organisms as surrogate hosts to study human disease, perhaps the more appropriate perspective is to view nematodes, insects, plants, and ameba as the relevant natural hosts in which the selection for and evolution of pathogenic traits occurs, and the ability to subsequently infect humans is a secondary effect of these interactions.
Materials and methods
Shotgun genome sequencing
Shotgun sequencing of PA14 was performed using 65,800 plasmids containing 2-4 kb fragments of genomic PA14 DNA, resulting in over 10-fold coverage. Sequence reads were assembled using the Phred, Phrap, and Consed tools [37–39]. Details of finishing methods are described in Additional data file 1.
A long-range PCR-based method was used to assess whether particular P. aeruginosa genomes contain a PA14-like or a PAO1-like arrangement of genome sequence between the two most distant copies of a large ribosomal RNA cluster (Figures 1 and 2a,b; also see Stover and coworkers ). Primer pairs used to amplify two products specific for the PA14 chromosomal arrangement are as follows: 807.LL_rev (CGAACTGGAGGAAGTCTTCG) + 2-5'_2 (CGAGGCTTTCGTCTATCCAG), and 837.LL (AACTGGTGGAGGGAGAAGGAT) + 803.RR.rev (TAGCCTTCAATTCCACCTGG). Primer combinations used to amplify two PAO1-specific products are as follows: 807.LL_rev + 803.RR.rev, and 837.LL + 2-5'_2. These diagnostic primer pairs were also used to survey 18 additional strains (Table 2) to determine whether the PA14 or the PAO1 orientation was more representative of the ancestral P. aeruginosa chromosome. Of the 18 strains tested, none generated a PAO1-specific amplification product. Using the primers for the PA14 arrangement, seven of the 18 strains yielded products with both primer pairs, nine of 18 had a strong PCR product with one of the two primer pairs and a weak product with the second primer pair, and two of 18 exhibited a weak product with one primer pair and no product with the second.
Sequence analysis and gene annotation
In determining GC content for PA14 and PAO1, standard deviations were determined using a continuous sliding window of 1 kb, and comparing the GC content for these regions with that of the whole genome. GC skew analysis was performed (using the formula GC skew = [n(G) - n(C)]/[n(G) + n(C)], where n [i] is the number of nucleotides, applied to 1 kb segments of each genome) to identify the peak cumulative GC skew as the likely position of the replication terminus . Global alignments of the PA14 and PAO1 genomes were performed using the MUMmer 3.0 software package  to identify strain-specific regions. Automated ORF predictions were made using a combination of the BLAST and Glimmer2 algorithms [41, 42]. Each predicted ORF was assigned a PA14 LocusName, beginning with 'PA14_' followed by five numerals. ORFs were numbered starting with PA14_00010 (dnaA), increasing in increments of 10 to allow for future insertions of additional genes or functional RNAs. Further details of ORF prediction and annotation are provided in Additional data file 1.
Microarray genomotyping of P. aeruginosastrains
70-mer oligonucleotides for microarrays were designed as described previously  for 285 PA14-specific sequences, 130 PAO1-specific sequences, 11 genes common to both strains (to serve as positive controls), and additional genes and negative controls (see Additional data file 1). Chromosomal DNA from PA14, PAO1 (the sequenced isolate), and 18 additional P. aeruginosa strains described previously  was digested with HaeIII and labeled with Cyanine-3 or Cyanine-5 using the MICROMAX ASAP labeling kit (part #MPS544001KT; PerkinElmer, Wellesley, MA, USA). Labeled samples were combined in random pairs and hybridized to the arrays. Observed intensities for replicate hybridizations for each strain were averaged and log2 ratios were computed to determine the presence/absence of each gene based on a cut-off determined independently for each sample (see Additional data file 1). For each gene, a Spearman's rank correlation coefficient was calculated to describe the relationship between the spectrum of present/absent calls and the determined rank order virulence in C. elegans (see below). Heirarchical clustering analysis of strain relationships was performed using Cluster 3.0  and Java Treeview . All microarray data are included in Additional data files 5 and 6 and have also been deposited in ArrayExpress (ArrayExpress E-MEXP-824).
C. eleganspathogenicity assays
The 20 P. aeruginosa strains analyzed by genomotyping were cultured overnight in Luria broth (LB) and assayed in the C. elegans killing system as described previously , with the exception that assays were performed using fer-15(b26)ts;fem-1(hc17) temperature-sensitive sterile mutant nematodes. The assays were carried out at the restrictive temperature (25°C) to prevent progeny formation and to allow the experiments to continue long enough to examine less pathogenic strains (for which progeny would typically overwhelm the assay plate by the end of the experiment). Bacterial strains are rank ordered from most to least virulent based on the position of the LT50 (the time at which 50% lethality was observed) in Figure 4a. In cases in which the LT50 was overlapping, rank orders were represented as a tied value for the strains in question. Strain CF27 is not included in the data shown in Figure 4a; however, its rank order in a similar experiment places it between strains E2 and S35004. The relative rank orders shown in Figure 4a are consistent with those observed in two additional experiments; although the absolute LT50 values increased or decreased for given strains between experiments, the relative rank orders remained consistent.
To screen for avirulent PA14 mutants, we utilized a nonredundant transposon insertion library . All available mutants in PA14-specific genes (349 in total) were grown in 150 ml of LB media in 0.5 ml 96-well Masterblocks (part #786261; Greiner, Monroe, NC, USA). Cultures were agitated at 225 rpm and incubated for 16 hours at 37°C. A volume of 10 μl of each culture was spotted onto slow killing agar in 2 wells of a six-well plate [13, 18], allowed to grow at 37°C for 24 hours, and left at room temperature for an additional 24 hours. Five L4 stage N2 nematodes were transferred manually to each well, and the number and age of progeny were recorded after 4 days at 25°C. The 349 mutants in PA14-specific genes were screened on two separate occasions, and any that were highly attenuated in either screen or were attenuated in both screens were re-examined in a secondary screen as described using wild-type N2 nematodes . Six mutants were tested for growth in minimal (M63) media and all were found to have wild-type growth curves (Table 3). These included four of the five mutants in PA14 genes completely absent in PAO1 (examined in Figure 4c).
Additional data files
The following additional data are available with the online version of this paper. Additional data file 1 is a document summarizing supplementary methods. Additional data file 2 is a table of PA14 gene annotations: the first worksheet contains the data, and the second worksheet contains a detailed description of the contents of the table and serves as a legend for the first worksheet. Additional data file 3 is a table showing the percentage of PA14 and PAO1 genes corresponding to each of 28 functional categories (shown for both the whole genome and strain-specific genes as described in Additional data file 4): the first worksheet contains the data, and the second worksheet contains a detailed description of the contents of the table and serves as a legend for the first worksheet. Additional data file 4 is a table summarizing properties of PA14-specific and PAO1-specific regions: the first worksheet contains the data, and the second worksheet contains a detailed description of the contents of the table and serves as a legend for the first worksheet. Additional data file 5 is a table describing the properties of oligos found on the genotyping microarray: the first worksheet contains the data, and the second worksheet contains a detailed description of the contents of the table and serves as a legend for the first worksheet. Additional data file 6 is a table of the microarray genotyping data for the 20 strains examined in this study (normalized average data and log2 ratio of normalized averages): the first worksheet contains the data, and the second worksheet contains a detailed description of the contents of the table and serves as a legend for the first worksheet. Additional data file 7 contains the coordinates of regions in PA14 or PAO1 that are common to four sequenced P. aeruginosa strains (PA14, PAO1, and two CF isolates: 2192 and C3719) and therefore comprise a putative P. aeruginosa core genome: the first worksheet contains the data, and the second worksheet contains a detailed description of the contents of the table and serves as a legend for the first worksheet.
Hacker J, Blum-Oehler G, Muhldorfer I, Tschape H: Pathogenicity islands of virulent bacteria: structure, function and impact on microbial evolution. Mol Microbiol. 1997, 23: 1089-1097. 10.1046/j.1365-2958.1997.3101672.x.
Oelschlaeger TA, Hacker J: Impact of pathogenicity islands in bacterial diagnostics. Apmis. 2004, 112: 930-936. 10.1111/j.1600-0463.2004.apm11211-1214.x.
Bentley SD, Parkhill J: Comparative genomic structure of prokaryotes. Annu Rev Genet. 2004, 38: 771-792. 10.1146/annurev.genet.38.072902.094318.
Ochman H, Davalos LM: The nature and dynamics of bacterial genomes. Science. 2006, 311: 1730-1733. 10.1126/science.1119966.
Doring D: Chronic Pseudomonas aeruginosa lung infection in cystic fibrosis patients. Pseudomonas aeruginosa as an opportunistic pathogen. 1993, New York: Plenum Press, 245-273.
Wood RE: Pseudomonas: the compromised host. Hosp Pract. 1976, 11: 91-100.
Stover CK, Pham XQ, Erwin AL, Mizoguchi SD, Warrener P, Hickey MJ, Brinkman FS, Hufnagle WO, Kowalik DJ, Lagrou M, et al: Complete genome sequence of Pseudomonas aeruginosa PA01, an opportunistic pathogen. Nature. 2000, 406: 959-964. 10.1038/35023079.
Ernst RK, D'Argenio DA, Ichikawa JK, Bangera MG, Selgrade S, Burns JL, Hiatt P, McCoy K, Brittnacher M, Kas A, et al: Genome mosaicism is conserved but not unique in Pseudomonas aeruginosa isolates from the airways of young children with cystic fibrosis. Environ Microbiol. 2003, 5: 1341-1349. 10.1111/j.1462-2920.2003.00518.x.
Wolfgang MC, Kulasekara BR, Liang X, Boyd D, Wu K, Yang Q, Miyada CG, Lory S: Conservation of genome content and virulence determinants among clinical and environmental isolates of Pseudomonas aeruginosa. Proc Natl Acad Sci USA. 2003, 100: 8484-8489. 10.1073/pnas.0832438100.
Spencer DH, Kas A, Smith EE, Raymond CK, Sims EH, Hastings M, Burns JL, Kaul R, Olson MV: Whole-genome sequence variation among multiple isolates of Pseudomonas aeruginosa. J Bacteriol. 2003, 185: 1316-1325. 10.1128/JB.185.4.1316-1325.2003.
Choi JY, Sifri CD, Goumnerov BC, Rahme LG, Ausubel FM, Calderwood SB: Identification of virulence genes in a pathogenic strain of Pseudomonas aeruginosa by representational difference analysis. J Bacteriol. 2002, 184: 952-961. 10.1128/jb.184.4.952-961.2002.
Rahme LG, Stevens EJ, Wolfort SF, Shao J, Tompkins RG, Ausubel FM: Common virulence factors for bacterial pathogenicity in plants and animals. Science. 1995, 268: 1899-1902. 10.1126/science.7604262.
Tan MW, Mahajan-Miklos S, Ausubel FM: Killing of Caenorhabditis elegans by Pseudomonas aeruginosa used to model mammalian bacterial pathogenesis. Proc Natl Acad Sci USA. 1999, 96: 715-720. 10.1073/pnas.96.2.715.
Jander G, Rahme LG, Ausubel FM: Positive correlation between virulence of Pseudomonas aeruginosa mutants in mice and insects. J Bacteriol. 2000, 182: 3843-3845. 10.1128/JB.182.13.3843-3845.2000.
Mahajan-Miklos S, Tan MW, Rahme LG, Ausubel FM: Molecular mechanisms of bacterial virulence elucidated using a Pseudomonas aeruginosa-Caenorhabditis elegans pathogenesis model. Cell. 1999, 96: 47-56. 10.1016/S0092-8674(00)80958-7.
Miyata S, Casey M, Frank DW, Ausubel FM, Drenkard E: Use of the Galleria mellonella caterpillar as a model host to study the role of the type III secretion system in Pseudomonas aeruginosa pathogenesis. Infect Immun. 2003, 71: 2404-2413. 10.1128/IAI.71.5.2404-2413.2003.
Rahme LG, Tan MW, Le L, Wong SM, Tompkins RG, Calderwood SB, Ausubel FM: Use of model plant hosts to identify Pseudomonas aeruginosa virulence factors. Proc Natl Acad Sci USA. 1997, 94: 13245-13250. 10.1073/pnas.94.24.13245.
Tan MW, Rahme LG, Sternberg JA, Tompkins RG, Ausubel FM: Pseudomonas aeruginosa killing of Caenorhabditis elegans used to identify P. aeruginosa virulence factors. Proc Natl Acad Sci USA. 1999, 96: 2408-2413. 10.1073/pnas.96.5.2408.
He J, Baldini RL, Deziel E, Saucier M, Zhang Q, Liberati NT, Lee D, Urbach J, Goodman HM, Rahme LG: The broad host range pathogen Pseudomonas aeruginosa strain PA14 carries two pathogenicity islands harboring plant and animal virulence genes. Proc Natl Acad Sci USA. 2004, 101: 2530-2535. 10.1073/pnas.0304622101.
PA14 Gene Annotations. [http://ausubellab.mgh.harvard.edu/cgi-bin/pa14/annotation/start.cgi]
Huang B, Whitchurch CB, Croft L, Beatson SA, Mattick JS: A minimal tiling path cosmid library for functional analysis of the Pseudomonas aeruginosa PAO1 genome. Microb Comp Genomics. 2000, 5: 189-203.
Pseudomonas Genome Project and PseudoCAP. [http://www.pseudomonas.com/]
Kim CC, Joyce EA, Chan K, Falkow S: Improved analytical methods for microarray-based genome-composition analysis. Genome Biol. 2002, 3: RESEARCH0065-10.1186/gb-2002-3-11-research0065.
PA14 Transposon Insertion Mutant Library. [http://ausubellab.mgh.harvard.edu/cgi-bin/pa14/home.cgi]
Liberati NT, Urbach JM, Miyata S, Lee DG, Drenkard E, Wu G, Villanueva J, Wei T, Ausubel FM: An ordered, nonredundant library of Pseudomonas aeruginosa strain PA14 transposon insertion mutants. Proc Natl Acad Sci USA. 2006, 103: 2833-2838. 10.1073/pnas.0511100103.
Romling U, Kader A, Sriramulu DD, Simm R, Kronvall G: Worldwide distribution of Pseudomonas aeruginosa clone C strains in the aquatic environment and cystic fibrosis patients. Environ Microbiol. 2005, 7: 1029-1038. 10.1111/j.1462-2920.2005.00780.x.
Sato H, Frank DW: ExoU is a potent intracellular phospholipase. Mol Microbiol. 2004, 53: 1279-1290. 10.1111/j.1365-2958.2004.04194.x.
Pukatzki S, Kessin RH, Mekalanos JJ: The human pathogen Pseudomonas aeruginosa utilizes conserved virulence pathways to infect the social amoeba Dictyostelium discoideum. Proc Natl Acad Sci USA. 2002, 99: 3159-3164. 10.1073/pnas.052704399.
Casadevall A, Steenbergen JN, Nosanchuk JD: 'Ready made' virulence and 'dual use' virulence factors in pathogenic environmental fungi: the Cryptococcus neoformans paradigm. Curr Opin Microbiol. 2003, 6: 332-337. 10.1016/S1369-5274(03)00082-1.
Levitz SM: Does amoeboid reasoning explain the evolution and maintenance of virulence factors in Cryptococcus neoformans?. Proc Natl Acad Sci USA. 2001, 98: 14760-14762. 10.1073/pnas.261612398.
Steenbergen JN, Casadevall A: The origin and maintenance of virulence for the human pathogenic fungus Cryptococcus neoformans. Microbes Infect. 2003, 5: 667-675. 10.1016/S1286-4579(03)00092-3.
Zhang Y, Lu H, Bargmann CI: Pathogenic bacteria induce aversive olfactory learning in Caenorhabditis elegans. Nature. 2005, 438: 179-184. 10.1038/nature04216.
Smith EE, Buckley DG, Wu Z, Saenphimmachak C, Hoffman LR, D'Argenio DA, Miller SI, Ramsey BW, Speert DP, Moskowitz SM, et al: From the cover: genetic adaptation by Pseudomonas aeruginosa to the airways of cystic fibrosis patients. Proc Natl Acad Sci USA. 2006, 103: 8487-8492. 10.1073/pnas.0602138103.
Lau GW, Goumnerov BC, Walendziewicz CL, Hewitson J, Xiao W, Mahajan-Miklos S, Tompkins RG, Perkins LA, Rahme LG: The Drosophila melanogaster toll pathway participates in resistance to infection by the gram-negative human pathogen Pseudomonas aeruginosa. Infect Immun. 2003, 71: 4059-4066. 10.1128/IAI.71.7.4059-4066.2003.
Kim DH, Feinbaum R, Alloing G, Emerson FE, Garsin DA, Inoue H, Tanaka-Hino M, Hisamoto N, Matsumoto K, Tan MW, et al: A conserved p38 MAP kinase pathway in Caenorhabditis elegans innate immunity. Science. 2002, 297: 623-626. 10.1126/science.1073759.
Liberati NT, Fitzgerald KA, Kim DH, Feinbaum R, Golenbock DT, Ausubel FM: Requirement for a conserved Toll/interleukin-1 resistance domain protein in the Caenorhabditis elegans immune response. Proc Natl Acad Sci USA. 2004, 101: 6593-6598. 10.1073/pnas.0308625101.
Gordon D, Abajian C, Green P: Consed: a graphical tool for sequence finishing. Genome Res. 1998, 8: 195-202.
Gordon D, Desmarais C, Green P: Automated finishing with autofinish. Genome Res. 2001, 11: 614-625. 10.1101/gr.171401.
Phred, Phrap, and Consed. [http://www.phrap.org/phredphrapconsed.html]
Kurtz S, Phillippy A, Delcher AL, Smoot M, Shumway M, Antonescu C, Salzberg SL: Versatile and open software for comparing large genomes. Genome Biol. 2004, 5: R12-10.1186/gb-2004-5-2-r12.
Delcher AL, Harmon D, Kasif S, White O, Salzberg SL: Improved microbial gene identification with GLIMMER. Nucleic Acids Res. 1999, 27: 4636-4641. 10.1093/nar/27.23.4636.
Salzberg SL, Delcher AL, Kasif S, White O: Microbial gene identification using interpolated Markov models. Nucleic Acids Res. 1998, 26: 544-548. 10.1093/nar/26.2.544.
Wang X, Seed B: Selection of oligonucleotide probes for protein coding sequences. Bioinformatics. 2003, 19: 796-802. 10.1093/bioinformatics/btg086.
Open Source Clustering Software. [http://bonsai.ims.u-tokyo.ac.jp/~mdehoon/software/cluster/software.htm#ctv]
Saldanha AJ: Java Treeview: extensible visualization of microarray data. Bioinformatics. 2004, 20: 3246-3248. 10.1093/bioinformatics/bth349.
Broad Institute Pseudomonas aeruginosa Database. [http://www.broad.mit.edu/annotation/microbes/pseudomonas_aeruginosa_c3719/]
Genomic DNA Labelling Protocol_Version B. [http://derisilab.ucsf.edu/pdfs/GenomicDNALabel_B.pdf]
We thank J Decker, W Brown, K Osborn, A Perera, R Elliott, and L Gendal for technical assistance during the sequencing of PA14; R Jackson for suggestions on genomic DNA preparations for microarray analysis; N El Massadi and J Frietas for suggestions on sample labelling, hybridizations, and scanning for microarray experiments; D Park for assistance with microarray data analysis; and E Mylonakis for discussions regarding the evolution of fungal pathogens. We especially acknowledge S Lory for encouragement and for many helpful discussions and suggestions. This work was funded by grants to FMA from the National Institutes of Health (Grant Numbers U01 HL66678 and R01 AI064332) and the Department of Energy (Grant Number DE-FG02-ER63445), a postdoctoral fellowship from the Jane Coffin Childs Memorial fund for Medical Research to DGL, and by support from the Harvard-Partners Center for Genetics and Genomics.
DGL and FMA were responsible for strategic planning and managing the overall project. LL, GG, KM, and RK were responsible for overseeing the raw sequencing, finishing, and initial assembly of PA14. DGL, NTL, SM, JH, and MS contributed to the raw sequencing effort. DGL performed manual edits to the assembly and additional finishing reactions. JMU and GW developed bioinformatic tools to analyze the sequence data. JMU, GW, and DGL performed the bulk of the bioinformatic analyses. DGL, JMU, LTD, NTL, JH, ED, and LF contributed to the manual annotation of PA14 genes. DGL designed the oligonucleotide array, performed the microarray experiments, and DGL and JMU analyzed the microarray data. NTL and RF performed the C. elegans assays. DGL and FMA wrote the paper. All authors discussed the results and commented on the manuscript.
Electronic supplementary material
Additional data file 3: The percentage of PA14 and PAO1 genes corresponding to each of 28 functional categories (shown for both the whole genome and strain-specific genes as described in Additional data file 4). (XLS 27 KB)
Additional data file 6: Microarray genotyping data for the 20 strains examined in this study (normalized average data and log2 ratio of normalized averages). (XLS 307 KB)
Additional data file 7: Coordinates of regions in PA14 or PAO1 that are common to four sequenced P. aeruginosa strains (PA14, PAO1, and two CF isolates: 2192 and C3719) and therefore comprise a putative P. aeruginosa core genome. (XLS 46 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Lee, D.G., Urbach, J.M., Wu, G. et al. Genomic analysis reveals that Pseudomonas aeruginosa virulence is combinatorial. Genome Biol 7, R90 (2006). https://doi.org/10.1186/gb-2006-7-10-r90
- Cystic Fibrosis
- Additional Data File
- Pathogenicity Island
- Genomic Content
- Cystic Fibrosis Isolate