Clustering of Pseudomonas aeruginosa transcriptomes from planktonic cultures, developing and mature biofilms reveals distinct expression profiles

Background Pseudomonas aeruginosa is a genetically complex bacterium which can adopt and switch between a free-living or biofilm lifestyle, a versatility that enables it to thrive in many different environments and contributes to its success as a human pathogen. Results Transcriptomes derived from growth states relevant to the lifestyle of P. aeruginosa were clustered using three different methods (K-means, K-means spectral and hierarchical clustering). The culture conditions used for this study were; biofilms incubated for 8, 14, 24 and 48 hrs, and planktonic culture (logarithmic and stationary phase). This cluster analysis revealed the existence and provided a clear illustration of distinct expression profiles present in the dataset. Moreover, it gave an insight into which genes are up-regulated in planktonic, developing biofilm and confluent biofilm states. In addition, this analysis confirmed the contribution of quorum sensing (QS) and RpoS regulated genes to the biofilm mode of growth, and enabled the identification of a 60.69 Kbp region of the genome associated with stationary phase growth (stationary phase planktonic culture and confluent biofilms). Conclusion This is the first study to use clustering to separate a large P. aeruginosa microarray dataset consisting of transcriptomes obtained from diverse conditions relevant to its growth, into different expression profiles. These distinct expression profiles not only reveal novel aspects of P. aeruginosa gene expression but also provide a growth specific transcriptomic reference dataset for the research community.


Background
Pseudomonas aeruginosa is an opportunistic pathogen which causes a variety of infections including bacteremia in burns patients, urinary tract catheter infections and chronic colonisation of the lungs of Cystic Fibrosis patients [1]. This bacterium is a common nosocomial contaminant and its persistence as a major cause of human disease is linked to its intrinsic resistance to many antibiotics. In addition to this natural resistance, in many infections P. aeruginosa survives as a biofilm -adhered communities which themselves are more resistant to antibiotic therapy than free-living or planktonic organisms [2,3]. Given the decreased therapeutic options for this organism, there is a clear need for novel approaches to treat or prevent P. aeruginosa infection, and the key to this is to understand more about its biology. A greater knowledge of the genetics of growth of P. aeruginosa will help unravel the molecular basis of infection and in particular the mechanics of survival in the biofilm mode of growth.
The success of P. aeruginosa as a human pathogen, and its ability to thrive in diverse environmental habitats, is thought to be due to its large and complex genome (6.3 million base-pairs, 5570 open reading frames) which includes a variety of virulence factors and also its ability to alternate between a free-living and biofilm state [4]. In order to understand this flexibility we have used wholegenome microarrays to obtain transcriptomes from six conditions relevant to the growth of this organism. These conditions were two planktonic phases (logarithmic and stationary) and multiple biofilm time points (8,14,24 and 48 hrs). Using confocal microscopy we have previously characterised the static biofilm model used for this study, and found that confluent biofilms develop in a sequential manner from microcolonies [5]. In a preliminary analysis of this transcriptomic dataset, we found major differences in gene expression when comparing actively growing growth states with stationary phase growth states; 19.4% of the PAO1 genome was differentially expressed between log phase and stationary phase planktonic cultures and >15.5% of the PAO1 genome was differentially expressed when developing biofilms (8 hrs) was compared to confluent biofilm time-points (14,24 and 48 hrs) [5]. In contrast we determined that developing and confluent biofilm transcriptomes were related to those of logarithmic and stationary phase planktonic cultures, respectively and that gene expression was conserved between confluent biofilms [5]. We also identified a unique set of 20 and 26 novel genes that were ≥ 2.5 fold up-regulated (P < 0.05) in developing biofilms and confluent biofilms respectively, when compared to all other conditions.
Although our previous analysis gave information regarding the relatedness of these different P. aeruginosa growth states, in this study we provide a more comprehensive analysis of this data using clustering. We performed such analysis using three different clustering methods: Kmeans, K-means spectral and hierarchical clustering. This analysis provides an overview of P. aeruginosa gene expression in conditions relevant to its lifestyle and allowed us to reveal the existence of distinct expression profiles in the transcriptomic data from these six different conditions. We were then able to observe gene expression that is unique to single conditions, or common to two or more culture conditions.

Growth conditions used for microarray analysis
The biofilm microarray data used in this study was derived from static biofilms grown on nitrocellulose filters placed on 20% LB agar [5]. Reduced strength LB (20%) was used in order to reflect a low nutrient environment which would be encountered in vivo. Visual monitoring using strain PAO1 tagged with the green fluorescent protein (GFP) and confocal scanning laser microscope (CSLM) imaging demonstrated that after filters were inoculated with single P. aeruginosa cells, large microcolonies had formed after 8 hours which then further developed into confluent biofilms after 14 hrs [5]. No further change in biofilm architecture was observed after 24 and 48 hrs incubation. Microarray analysis using P. aeruginosa Affymetrix GeneChip arrays was carried out using RNA derived from four biofilm time points; 8, 14, 24, and 48 hours [5].
The planktonic culture microarray data used in this study was derived from logarithmic and stationary phase planktonic cultures grown in 250 ml Erlenmeyer flasks. Briefly, 100 ml LB broth (20%) was inoculated with 1 ml PAO1 overnight LB culture (diluted 10 fold), and the resulting cultures were grown at 37°C with agitation (200 rpm). Microarray analysis was carried out on RNA extracted after 4 hours (OD 600 0.1) (Logarithmic phase -LP) and 24 hours growth (OD 600 0.3) (Stationary phase -SP).

Transcriptome analysis: clustering of microarray data reveals different expression profiles
On analysis of the microarray data, 1429 genes (26% of genome) [see Additional file 2 - Supplementary Table S1] were found to show low expression (under 50 signal units for all three replicates) for all six different conditions. The most abundant functional classes in this list were: 'transport of small molecules', 'hypothetical, unclassified, unknown', 'putative enzymes' and 'transcriptional regulators'. This list contained two characterised multidrug efflux resistance systems (MexC-MexD-OprJ and MexE-MexF-OprN) and a putative multidrug efflux resistance system (PA3521-PA3523). Also, in this group was the gene encoding the β-lactamase precursor (ampC) and katB, one of the four genes encoding possible catalases in the P. aeruginosa genome. Both, ampC and katB have been shown previously to be induced in biofilms by the presence of sub-inhibitory doses of imipenem [6,7] and by hydrogen peroxide [7], respectively.
We began the transcriptome analysis by filtering out these 1429 genes showing very low expression for all six different conditions. A 1-way ANOVA (assuming unequal variances) was then performed which removed a further 805 genes that did not vary significantly between different conditions [see Additional file 2 -Supplementary Table  S2]. As we are interested in the behaviour of genes between the six conditions and not within the 6 conditions, for the remaining 3315 genes we averaged the values of the three replicates for each of the six conditions, and then normalized these averaged profiles to zero mean and unit variance. We used the correlation coefficient as distance measure between gene profiles, and we clustered them using K-means clustering [8], K-means spectral clustering [9,10] and hierarchical clustering [8]. The first two methods were applied using a value of K (i.e. number of clusters) of 4, 7, 8, 9, 10, 12, 13, 15, 16 and 19. For the hierarchical clustering we tried the complete, centroid and ward linkage metrics for building the hierarchical tree and different values of the threshold for cutting it. An important finding is that for equivalent settings of the parameters of these different algorithms we obtained broadly the same groupings of the genes [see Additional file 1 -Supplementary Figure S1]. This revealed the existence of distinct expression profiles in the transcriptomic data from these six conditions. In the following we shall discuss the results obtained using K-means clustering with K = 10 ( Figure 1a and 1b). These expression profiles are defined in Table 1.

Confirmation of microarray transcription profiles by quantitative reverse transcriptase PCR (qRT-PCR)
In order to verify the results obtained using the P. aeruginosa GeneChip microarrays (Affymetrix), genes from four different clusters (PA0020, cluster 8; fliE, cluster 6; PA2184, cluster 1; PA5555, cluster 4) representing different transcriptional profiles were examined by qRT-PCR. This was performed on RNA samples taken from 4 conditions (LP planktonic culture, SP planktonic culture, 8 hr biofilm, and 48 hr biofilm) and the results are displayed relative to expression under LP planktonic culture conditions ( Figure 2). The gene expression profiles were consistent with those obtained from the microarray profiling experiments. Therefore, data from qRT-PCR provides independent verification of the microarray results.

Key functional classes in each cluster
The genes within each cluster were analysed and the percentage of twenty-six functional classes http://www.pseu domonas.com/ in each cluster was determined [see Additional file 1 -Supplementary Figure. S2]. Figure 3 shows the distribution of ten functional classes throughout the clusters.
Of the functional classes, 'hypothetical, unclassified, unknown' was the most abundant, comprising at least 25% of the total number of genes in each cluster. As expected, clusters 4 and 5 (up-regulation in LP planktonic culture and developing biofilms) contain many genes in functional classes associated with growth (ie. 'transcription, RNA processing and degradation'; 'translation, posttranslational modification, degradation', 'nucleotide biosynthesis and metabolism', 'DNA replication, recombination, modification and repair', 'cell division', 'amino acid biosynthesis and metabolism') ( Figure 3). There is also a large percentage of genes from the classes 'transcription, RNA processing and degradation', 'cell wall/LPS/capsule' and 'cell division' in clusters 7 (up-regulation in LP planktonic culture) and 9 (up-regulation in LP and SP planktonic culture) ( Figure 3). The 'motility and attachment' class features most prominently in clusters 6 (up-regulation in SP planktonic) and 8 (up-regulation in developing biofilms) ( Figure 3); many genes in cluster 6 are involved in the biosynthesis and assembly of flagellar (essential for swimming and swarming motility) whereas cluster 8, contains many genes involved in biosynthesis and assembly of type-IV pili (essential for twitching motility). The regulators of biosynthesis of these motility structures are also present in these groups. Perhaps, the most interesting results are the appearance of a large percentage of genes involved in 'chemotaxis' and genes from the functional class 'phage, transposon, or plasmid' in clusters 1 (up-regulation in the confluent biofilm time points -expression peak 14 hrs) and 5 (up-regulation in LP planktonic culture and developing biofilms), respectively (discussed later) ( Figure 3).
Gene lists for each cluster were compiled, together with their microarray expression results [see Additional file 3, Tables S3 -S12]. The following section describes some of the key genes found in the different expression profiles.  Table S5].

Genes up-regulated in confluent-biofilms, and genes upregulated in confluent-biofilms and in SP planktonic
Although the previously identified QS and RpoS regulated genes were identified under planktonic growth conditions [11][12][13], it is known that these regulatory systems play an important role in biofilm maturation [14][15][16]. Therefore, we utilized the available QS and RpoS transcriptomic datasets [11][12][13] to determine the distribution of these genes throughout the 10 clusters. The majority of genes previously found to be regulated by QS and RpoS in planktonic culture were also expressed in confluent biofilms. 423 of the 504 genes found to be up-regulated by RpoS [12] were in the gene list used for the cluster analysis and 80.1% of these genes were found distributed between clusters 1, 2 and 3 ( Figure 4a). 434 of the 534 genes (in total) shown to be up-regulated in two QS studies [11,13] were in the gene list used for the cluster analysis and 71% of these genes were found distributed amongst clusters 1, 2 and 3 ( Figure 4b). This percentage would probably have been higher had full strength LB been used rather than the reduced strength LB broth (20%) used to culture planktonic and biofilm cells in this study. Taking into account the very different conditions we used when compared to previous QS and RpoS transcriptomic studies [11][12][13] and also considering the cell densities of the different growth states we used, the distribution of the QS and RpoS regulated genes over the 10 clusters that we obtained would be expected.
As many P. aeruginosa virulence determinants are regulated by QS, it is not surprising that the following are found in clusters 1, 2 and 3; elastase precursor (lasB, cluster 1), phenazine biosynthetic genes (PA1901 -PA1905, PA4209, PA4210-PA4211 cluster 1), and rhamnosyltransferases involved in rhamnolipid biosurfactant synthesis (rhlAB cluster 3, rhlC cluster 2). Many of the genes in clusters 1, 2 and 3 are activated in response to stress in P. aeruginosa or are orthologues of stress induced genes in other Results of K-means clustering of the expression data (K = 10) Figure 1 Results of K-means clustering of the expression data (K = 10). The three replicates for each of the six conditions were averaged, and then normalized to zero mean and unit variance. A) The Eisen diagram [56] of the expression profiles grouped according to the results of the clustering. Red denotes that the value observed is above the mean of the observation across the dataset. The y-axis shows the number of genes in each cluster. The labels on the X-axis denote the different expression conditions: LP, LP planktonic culture; SP, SP planktonic culture; 8, 14, 24 and 48, biofilm time points. B) The plot of the profiles in each cluster. The X-axis shows the conditions in the same order as A). The thick yellow lines represent the learned cluster centres.

A) B)
Cluster No bacteria; general stress (PA1753, PA2190, cspD, rpoS, and rmf), nutrient stress (glgA, glgB, glgP -discussed in the next section), oxidative stress (katE -discussed in the next section), anaerobicity (anr), and osmolarity (osmC, osmE and osmY). Cluster 2 contains the gene encoding the transcriptional regulator -Anr, the anaerobic regulator of the arginine deiminase (ADI) and nitrate reductase pathways. Anr induces the ADI pathway in P. aeruginosa during oxygen limitation [17] enabling this organism to use arginine as an energy source. Three of the four genes of the ADI pathway (arcD, arcB, arcC) are present in clusters 1 and 3.
Interestingly, cluster 1 includes a six ORF cluster napEFD-ABC (PA1172 -PA1177) and PA4130 which are thought to be involved in nitrate reduction, and nitrite or sulfite reduction, respectively. P. aeruginosa is known to use nitrate or nitrite as a terminal electron acceptor under anaerobic conditions [18]. Biofilms are known to be heterogeneous environments and it is known that even if biofilms are cultured under aerobic conditions there is little oxygen deep within mature biofilms [19] and in the CF lung [20]. Therefore, we can speculate that in our study cells at the surface of the biofilm are in an aerobic envi- Comparison of microarray and quantitative reverse transcriptase PCR (qRT-PCR) transcription profiles

C) D)
Fold change -condition / LP -(microarray) , condition / LP -(qRT-PCR) Distribution of ten functional classes throughout the clusters created by K-means clustering (K = 10) Figure 3 Distribution of ten functional classes throughout the clusters created by K-means clustering (K = 10). Percentages were obtained by dividing the number of genes of a functional class in each cluster, by the total number of genes in that functional class. For a representation of all twenty-six functional classes in all ten clusters see Additional File 1 ( Figure. S2). Cluster 10

DNA replication, recombination, modification and repair
Pe rce ntage ronment whereas cells nearest the filter surface after 48 hours biofilm growth are in a localised anaerobic environment. It is therefore not surprising that the gene encoding the anaerobic regulator Anr was found in cluster 2 with expression peaking after 48 hours biofilm growth.
Cluster 1 contains PA1753, a homologue of an E. coli universal stress protein (UspA). UspA is thought to have a general protective function, is one of most abundant proteins in growth arrested cells and is stimulated by numerous conditions such as SP conditions, nutrient limitation, and exposure to oxidants and antibiotics [21]. Cluster 2 contains PA2622 which is 71% similar to the cspD gene product of Escherichia coli, a member of the CspA cold shock stress adaptation family of proteins [22]. Although E. coli CspD is not induced by cold shock it is dramatically induced by SP growth and glucose starvation [22].
Other general stress genes in cluster 2 are rpoS, PA2190 (discussed in the next section) and rmf. rmf encodes a putative protein which is 66% similar to the characterised ribosome modulation factor of E. coli and was very highly expressed in the confluent biofilm time points. In E. coli this protein has been shown to be expressed during the transition from LP to SP growth and in slow growing cells and its role is to associate with 70S ribosomes, converting them into a dimeric form and reducing protein synthesis [23,24].
Cluster 1 contains PA0059, PA4876 and PA4739, homologues of three E. coli proteins (osmC, osmE, osmY) that are osmotically induced. Interestingly, in E. coli it has been previously found that bacteria within biofilms encounter higher osmolarity conditions than bacteria in a planktonic state [25].
In cluster 2 are two genes that encode lectins, lecA (PA2570 -galactophilic lectin) and lecB (fucose-binding lectin) which are very highly expressed in the three confluent biofilm time points (>2.5 fold up-regulated, when compared to the three other conditions, P < 0.05). Interestingly, a lecB deficient mutant has previously been found to be significantly impaired in biofilm formation [26].
Three genes of the mexGHI-opmD efflux system were found in cluster 2. This efflux pump is thought to be involved in acyl-homoserine lactone (AHL) homeostatis and resistance to vanadium, but mutants defective in this efflux system do not display decreased antibiotic resistance [27].
P. aeruginosa has five (I to V) 'clusters' (which will be called 'sets' to avoid confusion) of chemotaxis genes. Genes from sets I, II and III are found in cluster 1 and genes from set V is found in cluster 3. Why genes from these chemotaxis sets are up-regulated under biofilm conditions in which motility biosynthetic genes (flagellar and type IV pili) are down-regulated is not clear. One explanation could be that these chemotactic sets are involved in environmental sensing and responding to stress under biofilm conditions [12].

A genome region associated with stationary phase growth (planktonic and biofilm)
Whilst studying the genes in clusters 1 and 2 it became apparent that there was a region of the genome (PA2134 -PA2190; 60.69 Kbp from position 2349488 to 2410181) that had the same expression profile; up-regulation in the confluent biofilm time points (14, 24 and 48 hrs) ( Figure 5). 32 of these 57 sequential genes have previously been found to be QS and/or RpoS regulated [11][12][13] using SP planktonic cultures of higher cell density than that of the SP planktonic cultures used in this study. Therefore, it is possible that this region could be associated with stationary phase growth (SP planktonic culture and/or confluent biofilms). There is no putative function assigned to the majority (40) of these genes. Interestingly this region contains genes with possible roles in the accumulation and breakdown of storage materials (glycogen, trehalose), oxidative stress protection (two catalases; katE, katN) and general stress (PA2190) ( Table 2). glgA, glgB, glgP encode enzymes with putative roles in glycogen synthesis and degradation (GlgA -a probable glycogen synthase, GlgB -a 1, 4-alpha-glucan branching enzyme and GlgP -a putative glycogen phosphorylase). An E. coli glgA mutant was impaired in biofilm formation [28] and GlgP has been suggested to have a role in the slow degradation of endogenous glycogen during long stationary conditions [29]. This region also contains three probable glycosyl hydrolases (PA2160, PA2162, PA2164) which may also be involved in the hydrolysis of storage polysaccharides and PA2190 which is 62% similar to the glucose or phos-phate starvation inducible gsiB general stress gene of Bacillus subtilis [30]. The two catalase genes present in this region; katE which encodes catalase HPII and katN which encodes a non-haem catalase, could aid the major housekeeping catalase KatA (PA4236) in SP conditions and in particular in the protection of biofilms against hydrogen peroxide.
Given that amongst this region there are genes that appear to be involved in general stress, glycogen accumulation and breakdown, and oxidative stress survival, it would be interesting to determine whether there are any genomic or expression differences in this region between different isolates (clinical and environmental) and if so, whether there is any difference in survival of these isolates under long term stationary conditions (biofilm and planktonic). A previous P. aeruginosa strain PAO1 microarray study suggested phage induction to be important for gene transfer in biofilms, as a region of genes (PA0718 -PA0727) of the functional class 'phage, transposon, or plasmid' were found to be highly up-regulated in mature biofilms when compared to planktonic culture [15]. Therefore, it would be expected that these genes would appear in clusters 1 and 2 (ie. the late biofilm time points). However, only one gene from this region was found in cluster 1 (PA0723). In contrast, we found cluster 5 to contain the highest percentage of 'phage, transposon, or plasmid' functional class genes, all of which were in a different region (PA0616-PA0624, PA0632-PA0636, PA0643-PA0646) to the previous study. Many of these genes have homology to P. aeruginosa R pyocin-related phage family genes (e.g. phi CTX, PS17). Phi CTX is a temperate phage that produces a pore-forming cytotoxin [33] and R pyocin is a type of bacteriocin that has a structure similar to contractile phage tails [34,35]. PA0622 and PA0623 encode proteins with high homology (91 and 84%, respectively) to contractile tail proteins of phage PS17. Therefore, it is possible that these phage genes which have recently been shown to be induced in response to hydrogen peroxide stress [36] may play a role in virulence under actively growing conditions (LP planktonic culture and 8 hr biofilm) which could be encountered in vivo.

Genes up-regulated in SP planktonic culture
Cluster 6 contains genes up-regulated in SP planktonic culture [see Additional file 3, Supplementary Table S8]. As discussed above this cluster is unique, as many of the genes displaying this expression profile are involved in the biosynthesis, assembly and regulation of flagellar motility.
Other genes that display this profile suggest magnesium (phoQ, oprH, mgtA, mgtE) and potassium (KdpFABCD) may be limiting in the SP planktonic culture conditions used in this study. Genes known to be stimulated in response to magnesium starvation are; the oprH-phoPQ operon which encodes an outer-membrane protein (OprH) and a two component regulatory system (PhoP/ Q) [37], mgtA (a ATP-dependent Mg 2+ transporter homologue) [38], and mgtE (a probable Mg transporter). Whilst the kdp gene cluster (PA1632 -1637) encodes a putative high affinity K + ion transport system. However, the expression profile of these genes differs greatly to the kdpABC operon of Staphylococcus aureus, which has been found to be up-regulated in biofilms compared to planktonic culture [39]. In this study the general stress gene recA was found to be up-regulated in SP planktonic culture, unlike in E. coli in which recA has been found to be up-regulated in biofilms [40]. This cluster also contains algU the P. aeruginosa orthologue of E. coli extreme stress sigma factor σ E .
Genes up-regulated in planktonic culture; LP, or LP and SP planktonic culture Cluster 7 contains genes up-regulated in LP planktonic culture whereas cluster 9 contains genes up-regulated in LP and SP planktonic culture [see Additional file 3, Supplementary Tables S9 and S11]. These two clusters contain genes encoding two characterised multidrug efflux systems (MexAB-OprM, and MexXY) between them. Genes from two putative multidrug efflux systems (PA2527, PA2528, PA1541) were also present in these two clusters together with a two-component regulatory system (PmrA-PmrB) that regulates resistance to polymyxin B and cationic antimicrobial peptides. This two-component system is involved in the regulation of a putative LPS modification operon (PA3552 -PA3559) [41] which is also present in cluster 7. Biofilms are well known for their decreased susceptibility to antibiotics when compared to planktonic cells. However, no previously characterised antibiotic efflux system was found to be up-regulated in the static biofilm system used in this study. It is well known that other factors such as slow bacterial growth, oxygen limitation and decreased penetration of antibiotic through the exopolysaccharide matrix play a role in the decreased susceptibility of biofilms to antibiotic therapy [42,43]. However, recent studies have shown biofilm specific induction of the MexCD-OprJ efflux system in the presence of azithromycin [44] and biofilm specific induction of genes from two probable efflux systems upon exposure to tobramycin [15].

Conclusion
The three clustering methods employed in this study (Kmeans, K-means spectral and hierarchical clustering) allowed us to sort this large transcriptomic dataset into genes with similar expression profiles. This enabled direct visualisation of gene expression common to two or more culture conditions and unique to single conditions. Confidence in our results was gained from the fact that for equivalent settings of the parameters all three clustering methods gave broadly the same grouping of expression profiles. Thus, the distinct expression clusters obtained together with the genes with low expression and consistent expression [see Additional file 2, Table S1 and S2] over all conditions, gives a complete picture of gene expression for this organism under conditions relevant to its lifestyle.
Here we have discussed a selection of genes representative of each cluster but given the large number of hypothetical genes and un-characterised orthologues in each cluster, there is still much to understand about the molecular and genetic basis of planktonic and biofilm growth of Pseudomonas aeruginosa. However, we believe that the availability of this dataset as a resource for the research community will accelerate the determination of the roles of these genes.
Structuring this data in clusters provided a clear illustration of the main expression profiles of many novel genes, such as genes from the functional classes 'chemotaxis' and 'phage, transposon, or plasmid' which are found in clusters 1 and 5 respectively. The cluster analysis also enabled us to observe the up-regulation of many individual genes in confluent biofilms, for example; 1) un-characterised orthologues associated with general stress (PA1753, PA2190, cspD, rmf), nutrient stress (glgA, glgB, glgP), oxidative stress (katE), and osmolarity (osmC, osmE and osmY), 2) the characterised anaerobic transcriptional regulator Anr (encoded by anr). It also facilitated the identification of a region of 57 sequential genes (PA2134 -PA2190) which demonstrate the same expression profile (up-regulation in the confluent biofilm time points). 32 of these genes have previously been found to be QS and/ or RpoS regulated [11][12][13].
Given the rise in antibiotic resistance, its ability to survive in the less antibiotic susceptible state -as a biofilm and the natural recalcitrance of this important human pathogen to many commonly used antibiotics, this dataset should aid in the pursuit of novel therapies against P.

Bacterial strains
Pseudomonas aeruginosa wild-type strain PAO1 [55] was used for both biofilm and planktonic culture expression analysis.
Biofilm modelling, confocal imaging, planktonic phase culturing, RNA isolation, cDNA labelling and microarray scanning Biofilm modelling, confocal imaging, planktonic phase culturing, RNA isolation, cDNA labelling and microarray scanning were performed as previously described [5]. The RNA samples used in this analysis were biological replicates, as they were derived from three independent cultures for each condition.

Quantitative reverse transcriptase PCR (qRT-PCR)
Planktonic and biofilm gene expression profiles generated using Affymetrix GeneChip arrays were verified using qRT-PCR. The same RNA samples used for the array analysis were also used for the qRT-PCR. 500 ng of RNA was reverse transcribed using the iScript cDNA Synthesis Kit as described by the manufacturer (BioRad Laboratories) and as a control the reaction was performed without reverse transcriptase (RT). Spectrophotometry was used to quantitate total cDNA and 67.5 ng of cDNA was used for each reaction using universal master mix (Applied Biosystems). Gene specific primers and probes were designed using Primer Express Software (Applied Biosystems), the sequences of which are shown in Table 3. Probes (purchased from Sigma-Genosys) were labeled 5' with 6-carboxyfluorescein (FAM) as the reporter and 3' with 6carboxytetramethylrhodamine as the quencher (TAMRA).
Each reaction was carried out in triplicate (in the presence of a no template control) using a model ABI 7900 HT thermal cycler (Applied Biosystems); thermal cycling parameters were 2 minutes at 50°C, 10 minutes at 95°C, and 40 cycles of 95°C for 15 seconds and 60°C for 1 minute. Calculation of expression value's for each gene (where cycle threshold is abbreviated as Ct): Ct (+RT) -Ct (-RT) = cCt (corrected CT). ΔCt = cCt (sample gene) -cCt (control gene, clpP). Expression value = 2 -ΔCt . clpP was used as a control as it was found to be consistently expressed under all six conditions using microarray profiling [see Additional file 2, Supplementary Table S2].

Cluster analysis of microarray data
Before data processing, all microarray data was globally scaled to set the average signal intensity of each array to a target signal of 100.
We began the transcriptome analysis by filtering out 1429 genes showing very low expression for all six different conditions (under 50 signal units for all three replicates). A 1-way ANOVA (assuming unequal variances) was then performed which removed a further 805 genes that did not vary significantly between different conditions [see Additional file 2, Supplementary Table S2]. For the remaining 3315 genes, we averaged the values of the three replicates for each of the six conditions, and then normalized these averaged profiles to zero mean and unit variance. We used the correlation coefficient as distance measure between gene profiles, and we clustered them using K-means clustering [8], K-means spectral clustering [9,10] and hierarchical clustering [8]. The first two methods were applied using a value of K (i.e. number of clusters) of 4, 7, 8, 9, 10, 12, 13, 15, 16 and 19. For the hierarchical clustering we tried the complete, centroid and ward linkage metrics for building the hierarchical tree and different values of the threshold for cutting it. An important finding is that for equivalent settings of the parameters of these different algorithms we obtained broadly the same groupings of genes. This revealed the existence of distinct expression profiles in the transcriptomic data from these six conditions. As within the paper we discuss the results obtained using K-means clustering with K = 10, in the following section we discuss some details of Kmeans clustering. K-means is a standard algorithm for partitioning N data points into K disjoint subsets S j containing N j data points. The algorithm consists of a reestimation procedure that attempts to minimize the sumof-squares criterion: where x n is a vector representing the nth data point and μ j is the mean of the data points in S j .
When using this algorithm for clustering gene expression profiles, each profile was considered a data point, and the distance in (A) was the correlation coefficient between x n and μ j . This amounts to say that we are learning the expression profiles which are the centers of the clusters (the thick yellow line Figure 1b).
The clearest representation of the expression profiles within the dataset was obtained with a K = 10 (Figure 1a and 1b). The 10 clusters are defined in Table 1. The clusters found were quite stable with respect to different random initialization of the algorithm. Gene lists for each cluster were compiled, together with their microarray expression results [see Additional file 3, Tables S3 -S12].
We also performed the analysis directly on the experimental data, without averaging the replicates. As expected, the results obtained [see Additional file 1, Supplementary Figure S3] were similar to the ones presented in Figure 1 although more noisy.