<title>Background</title> <p>Nosocomial infections caused by opportunistic pathogens are one of the most important health problems in developed countries. Depending on the geographic location, <italic>P. aeruginosa </italic>is the first or second causative agent of nosocomial infections [##REF##11714042##1##,##REF##12146499##2##]. <italic>P. aeruginosa </italic>infects patients suffering from AIDS, people at intensive care units, and burned people among others, and is the major cause of morbidity and mortality in patients with cystic fibrosis, the most prevalent hereditary disease in Caucasian populations [##REF##14555458##3##]. A successful infection by this type of pathogens depends on the interplay of multiple factors including the susceptibility of the host, the virulence of the strain and its resistance to antibiotics [##REF##12364374##4##]. Previous work has shown that the physiological fitness and the virulence of <italic>P. aeruginosa </italic>and other opportunists are affected by the expression of antibiotic resistance mechanisms such as MDR-pumping systems [##REF##14615205##5##, ####REF##12407121##6##, ##REF##14739147##7##, ##REF##15687203##8####15687203##8##].</p> <p>The pathogenicity of <italic>Pseudomonas aeruginosa </italic>involves various components operating at different levels. The flagella and <italic>pili </italic>facilitate contact with the bacterium's cell target and play a role in its adhesion, which is a critical step in the infection [##REF##9423837##9##,##REF##14607814##10##]. After contact, the type III secretion system is able to inject into the cytoplasm of the target cell a series of cytotoxic molecules that act at various levels. The mechanism of action involves, in many cases, the presence of host cofactors still unidentified [##REF##16427231##11##]. Other virulence factors involve products secreted into the extracellular medium by systems I and II such as elastase, alkaline phosphatase and exotoxin A among others. The expression of many of these virulence factors is regulated by a mechanism of bacteria-to-bacteria cell signalling known as quorum-sensing [##REF##9866731##12##]. Despite the functional and genomic similarity among different <italic>P. aeruginosa </italic>strains [##REF##11207762##13##,##REF##15205425##14##], some differences in their pathogenicity have been observed [##REF##17038190##15##]. For example, the clinical isolate PA14 is more virulent than PAO1 in a wide range of hosts [##REF##17038190##15##, ####REF##7604262##16##, ##REF##9892699##17####9892699##17##]. It has been shown that the genome of PA14 contains two pathogenicity islands that are not present in PAO1 and it has been proposed that the virulence in this organism (and the difference between PA14 and PAO1) is the result of a pool of pathogenicity genes interacting in various combinations in different genetic backgrounds [##REF##17038190##15##]. In spite of these suggestions, the cause of the different virulence behavior of PAO1 and PA14 is not yet fully understood.</p> <p>Although most of the work on pathogenesis has been focused on understanding the bacterial factors that render a virulence phenotype, increasing attention is being paid to the host and those aspects connected to the susceptibility or resistance to infection. Understanding the host-pathogen relationship, at both the cellular and molecular level, is essential to identify new targets and develop new strategies to fight infection. Molecular analysis of host-pathogen interactions would benefit from the use of model systems allowing a systematic study of the factors involved. In this regard the social amoeba <italic>D. discoideum </italic>has proven particularly useful for its ease of handling, genetic tractability [##REF##18042255##18##, ####REF##17178462##19##, ##REF##11867744##20##, ##REF##12003944##21##, ##UREF##0##22####0##22##] and fully sequenced genome [##REF##15875012##23##].</p> <p><italic>D. discoideum </italic>is a soil microorganism that feeds on bacteria by phagocytosis. The interaction between bacteria and their natural predators (<italic>Dictyostelium</italic>, other protists and worms) is believed to have shaped both predators [##REF##17699615##24##] and bacterial evolution. As a consequence, some of the mechanisms developed by bacteria to avoid the activity of their natural predators in the environment might have been adapted later in evolution to allow the infection of higher organisms such as humans [##REF##15679834##25##]. Specifically, it was found that the quorum-sensing mechanisms and type III secretion, which are essential factors in the infectivity to humans are also responsible for the infectivity of <italic>P. aeruginosa </italic>in <italic>D. discoideum </italic>[##REF##18042255##18##,##REF##11867744##20##,##REF##12003944##21##].</p> <p>Our previous studies have shown the utility of this model system of infection to analyze the virulence of other opportunistic pathogens like <italic>Stenotrophomonas maltophilia </italic>[##REF##14739147##7##]. It has been also demonstrated the validity of <italic>D. discoideum </italic>as a model of infection by intracellular pathogens such as <italic>Legionella, Cryptococcus and Mycobacterium </italic>[##REF##17178462##19##,##UREF##0##22##]. Consequently, the conservation of the mechanisms of infection needed to infect mammals and <italic>D. discoideum </italic>in a wide variety of pathogens reinforces the use of this system as a valid model to study host-pathogen relations. We have used whole-genome <italic>D. discoideum </italic>microarrays to study global host transcription upon infection with <italic>Pseudomonas aeruginosa </italic>PAO1 and PA14 to determine whether or not transcription is being affected as a defense mechanism or altered in the benefit of the pathogen.</p>
<title>Methods</title> <title><italic>D. discoideum </italic>growth and plating assays</title> <p>Dictyostelium AX4 cells were grown axenically in HL5 medium (14 g/l tryptone, 7 g/l yeast extract, 0.35 g/l Na<sub>2</sub>HPO<sub>4</sub>, 1.2 g/l KH<sub>2</sub>PO<sub>4</sub>, 14 g/l glucose, pH 6.5) or in association with <italic>Klebsiella aerogenes </italic>on SM plates (10 g/l glucose, 1 g/l yeast extract, 10 g/l peptone, 1 g/l MgSO<sub>4</sub>·6H<sub>2</sub>0, 1.9 g/l KH<sub>2</sub>PO<sub>4</sub>, 0.6 g/l K<sub>2</sub>HPO<sub>4</sub>, 20 g/l agar, pH 6.5) [##REF##3298997##39##].</p> <p>For the nutrient SM-plating assay <italic>Pseudomonas aeruginosa </italic>(PAO1 and PA14) and <italic>Klebsiella aerogenes (KA) </italic>were grown overnight in LB. After washing, bacteria were resuspended with PDF (20 mM KCl; 9 mM K<sub>2</sub>HPO<sub>4</sub>, 13 mM KH<sub>2</sub>PO<sub>4</sub>, 1 mM CaCl<sub>2</sub>; 1 mM MgSO<sub>4</sub>; pH: 6.4) and the optical density (OD) determined at 600 nm. After adjusting to 0.5 OD units, 300 μl of <italic>Klebsiella </italic>and <italic>Pseudomonas </italic>at the indicated proportions were plated in SM-agar plates with approximately 100 <italic>D. discoideum </italic>cells.</p> <p>For non-nutrient plating assay bacteria were grown as before, washed and resuspended in PDF (at an OD of 15 units at 600 nm. 100 μl of bacteria at the indicated proportions were mixed with <italic>D. discoideum </italic>cells and deposited in a drop over a PDF-based non-nutrient agar plates.</p> <title>Microarrays</title> <p><italic>Dictyostelium </italic>cells (5 × 10<sup>7 </sup>cells) were deposited in 10 ml of HL5 (without antibiotics) in shaking culture and exposed during 4 hours to 1.0 OD (approximate multiplicity of infection: 1000) of <italic>Klebsiella aerogenes </italic>as a control or to a mixture of <italic>Pseudomonas aeruginosa </italic>(PAO1 and PA14) and <italic>Klebsiella </italic>(used in excess to provide similar conditions of food supply). The proportion of <italic>Klebsiella </italic>to <italic>Pseudomonas </italic>(either PAO1 or PA14) was 7:3. Previous experiments showed that these proportions are adequate for a clear inhibition of <italic>Dictyostelium </italic>growth in plating assays (see the results section). <italic>Dictyostelium </italic>cells were separated from the bacteria by gentle centrifugation (twice at 1000 rpm, 5 minutes) and RNA isolated by Trizol (Life Technologies) according to manufacturer's instructions. Three independent biological experiments were performed making a total of three independent treatments for each condition (<italic>Klebsiella </italic>control, PAO1, and PA14). RNA from the three treatments were paired in all three combinations (Control/PAO1; Control/PA14 and PAO1/PA14) and hybridized to three different arrays. The same was performed for the other two biological replicates making a total of 9 microarrays hybridized. One of the three biological replicates was hybridized in the opposite dye orientation to the other two. The arrays, and protocols for labelling, hybridisation and scanning were as previously described [##REF##18430225##40##]. Background fluorescence was subtracted [##REF##11911795##41##], linear models were fitted and the significance of apparent changes in expression was assessed using limma [##REF##15657102##42##,##UREF##1##43##]. The data were normalised within each array with the printtip loess algorithm to counteract scanning and spatial biases, and further between each array to normalise mean absolute deviations using the 'scale' algorithm [##REF##14597310##44##].</p> <p>Preliminary analysis using ANOVA methods indicated that many genes had significant differences in expression between treatments, so we proceeded to examine each pairwise contrast in turn. We filtered out less reliable data by selecting genes with p-values adjusted for multiple testing [##UREF##2##45##] less than 0.05, making use of the moderated t statistics calculated by the eBayes function of limma. Since many genes passing this cutoff showed small changes in expression, we filtered further by absolute log2ratio. The commonly-used cutoff of greater than 2-fold change would have excluded a large number of genes that appeared to change in expression quite consistently, so we used the less stringent criterion of absolute log2ratio &gt;0.5. A lower log2ratio would have included genes with differences in expression too small to be corroborated by other methods. The array design is available from ArrayExpress [##UREF##3##46##] under the accession A-SGRP-3. The array experiment was deposited in the ArrayExpress database under the accession E-TABM-464.</p> <title>Quantitative PCR</title> <p>The same RNA samples subjected to microarray study were used as templates for retrotranscription with High Capacity cDNA Reverse transcription kit (Applied Biosystems) using 250 ng of each RNA in a final volume of 20 μl. For each sample, a triplicated blank was used. This cDNAs served as template in the PCR reaction carried out in 7900 HT Fast Real- Time PCR System using Power Sybrgreen PCR Master Mix 2× with 300 nM oligonucleotides concentration in a final volume of 10 μl. The results were acquired with SDS 2.3 software by Applied Biosystems and handled with EXCEL software by Microsoft. A total of seven genes were studied for each sample and their amount were related to one control gene, DDB0217951, whose expression is not affected by the treatments with the different strains.</p>
<title>Results</title> <title><italic>Pseudomonas aeruginosa </italic>strains PAO1 and PA14 show a different virulence behavior in <italic>D. discoideum</italic></title> <p>PAO1 and PA14 are two clinical isolates of <italic>P. aeruginosa </italic>frequently used as model strains to analyze the virulence of this bacterial pathogen. Since they behave differently in some aspects dealing with the expression of virulence determinants, we wanted to compare the differential response of the host to these strains. For this purpose, we made use of <italic>D. discoideum </italic>as a model for virulence. As a first step a plating assay of virulence was set up. Figure ##FIG##0##1## shows a representative experiment of three independent assays in which <italic>D. discoideum </italic>cells were grown in association with bacteria on nutrient SM plates. <italic>Klebsiella aerogenes</italic>, a non -pathogenic bacteria, was used as an appropriate food supply and <italic>P. aeruginosa </italic>mixed at the indicated proportions. An effect in the size of the clearing plaques could already be seen when only 3.5% of <italic>P. aeruginosa </italic>cells were mixed with 96.5% of <italic>K. aerogenes </italic>cells and this effect was even clearer using 17% of <italic>P. aeruginosa </italic>cells. When the behavior of the strains was analyzed in more detail, it was found that PAO1 is reproducibly more permissive than PA14 as observed by the higher growth of <italic>D. discoideum </italic>on PAO1. The differences in the area of the cleared bacterial lawn between PAO1 and PA14 were measured for the condition corresponding to the 3.5 % mixture. The average area and the standard deviation were 1.65 ± 1.2 mm<sup>2 </sup>for PAO1 and 0.11 ± 0.07 mm<sup>2 </sup>for PA14 (the number of clear plaques measured in each condition was 50). The significance of differences between groups as determined by Student's <italic>t</italic>-test was p &lt; 10<sup>-8</sup>. To further confirm these results a different plating assay was performed on non-nutrient agar. PAO1, PA14 and <italic>K. aerogenes </italic>were previously grown in LB overnight, washed out of the media by centrifugation and deposited with <italic>D. discoideum </italic>cells in agar plates at the indicated proportions. Under these conditions the difference in the virulence between PAO1 and PA14 was even more evident as shown in a representative experiment in Figure ##FIG##1##2##. Interestingly PAO1 is permissive to <italic>D. discoideum </italic>growth under these non-nutrient conditions. However, PA14 still shows a strong virulence against <italic>D. discoideum</italic>. All together these results suggest that PA14 is more virulent than PAO1 in the <italic>D. discoideum </italic>model of virulence.</p> <title><italic>Pseudomonas aeruginosa </italic>induces a specific gene expression response in <italic>Dictyostelium</italic></title> <p>Little is known about the interplay between the host and the pathogen in terms of gene expression responses. We wanted to determine if there is a specific gene expression response of <italic>D. discoideum </italic>to their interaction with <italic>P. aeruginosa</italic>. <italic>D. discoideum </italic>cells were exposed to <italic>P. aeruginosa </italic>strains PAO1 and PA14 mixed with an excess of <italic>K. aerogenes </italic>in HL5 for 4 hours. <italic>K. aerogenes </italic>alone was used as a control to which the gene expression levels were compared. RNA was extracted from <italic>D. discoideum </italic>and used to study the global pattern of gene expression using whole-genome <italic>D. discoideum </italic>microarrays (see Additional file ##SUPPL##0##1## for the complete data). Using a P &lt; 0.05 cutoff, there were 752 genes whose expression was significantly different between the PAO1-treated cells and the controls and 624 genes between PA14-treated cells and controls (Table ##TAB##0##1## summarizes the results at different P values and log-ratios). The heat map shown in Figure ##FIG##2##3## indicates that the responses were broadly comparable between the two strains with very few genes oppositely altered in the two conditions. The differences in the gene expression are approximately in the range between four-fold repression and three-fold induction (log-ratios between of -2 to +1.5 as shown in the histogram of Figure ##FIG##2##3##). These results were validated by real time PCR of the same samples used for the transcriptomic assays, measuring the expression of 7 representative genes that were up-regulated or down-regulated in the different conditions. Figure ##FIG##3##4## shows a good correlation between the data obtained from the microarray transcriptomic experiment as compared with that obtained by quantitative RT-PCR. Although the log-ratio changes in the gene expression showed some differences the overall trend were consistent, supporting the reliability of our data.</p> <title>Common and specific responses of <italic>D. discoideum </italic>to the infection with PAO1 and PA14 strains</title> <p>As shown in Table ##TAB##0##1## there were 364 genes that showed similar differential regulation with both bacterial strains compared with the controls (labeled as PAO1+PA14 vs control). We have considered in the analysis those genes showing differences in log-ratios that are higher than +0.5 or lower than -0.5. Interestingly the expression of another group of 169 genes (labeled as PAO1 vs PA14) was different depending on whether the infecting strain was PAO1 or PA14. We have studied in detail both groups by manual annotation and categorizing using the extended categorization for <italic>D. discoideum </italic>previously described [##REF##15010511##26##]. Genes of unknown function and those showing weak homologies were not included in the list. Table ##TAB##1##2## contains the genes that were similarly regulated upon infection with any of both strains, and in Figure ##FIG##4##5## the genes are categorized by function (see Additional file ##SUPPL##1##2## for the complete data). The first interesting conclusion from this experiment is the existence of a common transcriptional response that affects many different genes that are involved in a wide range of functions. The proportion of the genes that were downregulated by the treatment with both strains of <italic>P. aeruginosa </italic>is higher (258 genes) compared to those upregulated (106 genes). This difference is more evident in categories such as stress response and transport (Figure ##FIG##4##5##).</p> <p>Table ##TAB##2##3## displays the <italic>D. discoideum </italic>genes whose expression changed differentially between PAO1 and PA14 infection and Figure ##FIG##5##6## shows the number of genes in each category (see Additional file ##SUPPL##2##3## for complete data). In general a higher proportion of the genes showed a higher level of expression by PAO1 infection (109 genes) when compared with the levels observed by PA14. On the other hand, 60 genes behaved oppositely showing lower levels of expression upon exposure to PAO1 compared to those levels obtained after PA14 infection. Interestingly, all the genes represented in the categories stress-response and protein targeting had a higher level of expression in the cells exposed to PAO1 compared to PA14. The behavior of these genes in comparison with the control is also displayed in Table ##TAB##2##3##.</p>
<title>Discussion</title> <p><italic>P. aeruginosa </italic>is able to infect <italic>D. discoideum </italic>cells using several virulence traits that are similar to those used to infect mammalian cells and other hosts [##REF##18042255##18##]. The clinical <italic>P. aeruginosa </italic>isolates PA14 and PAO1 have been used independently to study the infection of <italic>Dictyostelium </italic>by <italic>Pseudomonas </italic>in two different laboratories [##REF##11867744##20##,##REF##12003944##21##]. However, no direct comparison had been reported so far between these strains in this pathogenicity model. We now report that PA14 is indeed more virulent in <italic>D. discoideum </italic>using different plating assays. Since <italic>P. aeruginosa </italic>is phagocytosed at much lower rate than the non-pathogenic <italic>K. aerogenes </italic>(commonly used to grow <italic>Dictyostelium</italic>) [##REF##11867744##20##], the assays were designed to provide sufficient food to <italic>D. discoideum </italic>to avoid cell starvation. Thus, <italic>K. aerogenes </italic>was always used in excess together with the pathogenic strains. In the first assay (Figure ##FIG##0##1##) a nutrient plate was used to allow the growth of bacteria and <italic>D. discoideum </italic>simultaneously. Under these conditions the presence of PA14 inhibits <italic>D. discoideum </italic>growth to a greater extent than PAO1. To avoid differences in the growth rates between bacteria that might alter their final proportions, a non-nutrient assay was performed (Figure ##FIG##1##2##). In these experiments, <italic>D. discoideum </italic>feed on bacteria that have been previously grown and deposited at different proportions in non-nutrient agar. Interestingly, PAO1 is not virulent in this condition suggesting that bacterial growth is necessary for the expression of the virulence in this strain. However, even in these conditions PA14 is capable of inhibiting <italic>D. discoideum </italic>growth. Some studies have suggested that PA14 pathogenicity is multifactorial and required the action of multiple virulence mechanisms [##REF##17038190##15##,##REF##11807055##27##]. These differences between strains prompted us to study the transcriptional profile of <italic>D. discoideum </italic>upon infection with PAO1 and PA14 to gain insights not only into the possible common transcriptional response but also into any specific response that could explain the observed differences in their virulence. Pilot experiments showed that 4 hours of exposure of <italic>D. discoideum </italic>cells to either <italic>Pseudomonas aeruginosa </italic>strains did not result in any apparent cell death or change in cell morphology (data not shown). Since we wanted to study the early transcriptional response we chose this short time of exposure to avoid changes due to cell death. The existence of a rapid gene expression response between 1–6 hours upon exposure of <italic>D. discoideum </italic>cells to <italic>Legionella</italic>, an intracellular pathogen, has also been described [##REF##16469056##28##].</p> <p>Our results show the existence of a common transcriptional response to the infection with <italic>P. aeruginoa </italic>PAO1 or PA14 that affects 364 genes grouped in many different cellular functions. The complexity of the observed transcriptional changes could be the result of the induction of <italic>D. discoideum </italic>defensive responses or triggered by <italic>P. aeruginosa </italic>to make a less hostile cell environment that would support a better survival of the pathogen. In this scenario downregulation of genes involved in stress response might be beneficial for a successful infection. Interestingly, we have observed a clear down-regulation of genes dedicated to stress in the common response to PAO1 and PA14 but also in the specific response to the more virulent strain PA14. For example the gene coding for Strictosidine synthase (DDB0185428) which is involved in the synthesis of alkaloids related to defense mechanisms in plants [##REF##17401749##29##], Trap1 (DDB0169033) that plays a central role in cell cycle regulation and differentiation [##REF##15507488##30##] or the genes coding for lysozymes involved in bacterial degradation (DDB0167491) [##REF##15640146##31##], to mention just a few.</p> <p>Besides stress response other categories are affected by an overall downregulation such as metabolism, translation and transport facilitation. A subcategorization of the genes included in metabolism (see supplementary table ##TAB##1##2##) showed that all the genes coding for proteins involved in nucleotide metabolism were downregulated in the common response suggesting an effect on cell proliferation. Moreover, a number of other downregulated genes are directly involved in cell growth as demonstrated functionally in previous studies. This is the case for example of DDB0192001 (ppkA, polyphosphate kinase), whose disruption leads to reduced growth on bacteria [##REF##17940044##32##,##REF##15701689##33##], DDB0186120 (gcsA, glutamylcysteine synthetase), which is essential for cell growth as mutants in the gene are not viable in the absence of glutathione [##REF##15993406##34##]. DDB0168860 (sgkA, sphingosine kinase) that is involved in cell proliferation [##REF##15643073##35##], among others that have been annotated in Supplementary Table ##TAB##1##2##, ##TAB##2##3##.</p> <p>Two different expression microarray analysis in mammals upon infection with <italic>Pseudomonas aeruginosa </italic>have been reported. In the first report epithelial cells were exposed to the pathogen for 3 hours, a short exposure similar to our experimental design. Unfortunately the number of genes represented in the array was very limited (1500 cDNAs) [##REF##10931941##36##]. Only 22 genes were differentially regulated and we have not found any homologous gene in common. The other work reported the analysis of <italic>Pseudomonas aeruginosa </italic>corneal infection using an oligonucleotide microarray [##REF##12882789##37##]. This experiment is not directly comparable to ours since a long exposure to the pathogen (1 day) was performed to assure an infection process. As a consequence most of the regulated genes were associated with the immune response and apoptosis, aspects that are not present in <italic>Dictyostelium</italic>.</p> <p><italic>D. discoideum </italic>is also susceptible to the infection by <italic>Legionella pneumophila</italic>, a facultative intracellular parasite, which uses different infective mechanisms from <italic>P. aeruginosa</italic>. It is important to note here that the transcriptional response of <italic>D. discoideum </italic>upon infection with <italic>Legionella </italic>[##REF##16469056##28##] was essentially different to the one we report for <italic>P. aeruginosa</italic>. Only 8 genes were found to be altered in both experiments (DDB0186332, DDB0219578, DDB0167879, DDB0205386, DDB0185740, DDB0167345, DDB0201617, DDB0202615). This indicates that the host response is rather specific of the type of infection and the bacterial pathogen involved. Nevertheless, some responses can be also common. For instance, DDB0202615 (nramp1, natural resistance-associated macrophage protein) whose expression is downregulated in PAO1 and PA14, plays an important role in Legionella infection since the null mutant has increased sensitivity to the infection [##REF##16445684##38##]. Nramp1 transports metal cations out of the phagolysosome in an ATP-dependent process. This activity is believed to be necessary to avoid the growth of intracellular pathogens and might also contribute to the efficient killing of other bacterial pathogens.</p> <p>The variety of genes whose expression is altered by <italic>P. aeruginosa </italic>infection suggests a complex scenario in which a combined downregulation of the expression of some of the mentioned genes might affect <italic>D. discoideum </italic>fitness thus favoring the infection. The precise role of these genes in the pathogenesis and the mechanisms that regulates their expression will promote further investigation.</p>
<title>Conclusion</title> <p>Our results showed that <italic>P. aeruginosa </italic>PA14 is more virulent than PAO1 in the <italic>D. discoideum </italic>model using different plating assays. The transcriptional responses of <italic>D. discoideum </italic>infected by either <italic>P. aeruginosa </italic>PAO1 or PA14 were analyzed by whole-genome microarrays and the expression of 364 genes changed similarly upon infection with any of both strains as compared with the control. Effects on metabolism, signaling, stress response and cell cycle can be inferred from the genes affected. Interestingly there were 169 genes differentially regulated between PAO1 and PA14, and this differential response might contribute to the different virulence behavior displayed by these two model strains. This is a starting point to begin to understand the complex relationships between environmental opportunistic pathogens and their natural hosts. Besides, our data support the idea that the host responses to different isolates of the same bacterial pathogen are largely different, thus indicating that the crosstalk between the pathogen and its host is more specific and more complex than previously thought.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p><italic>Pseudomonas aeruginosa </italic>is one of the most relevant human opportunistic bacterial pathogens. Two strains (PAO1 and PA14) have been mainly used as models for studying virulence of <italic>P. aeruginosa</italic>. The strain PA14 is more virulent than PAO1 in a wide range of hosts including insects, nematodes and plants. Whereas some of the differences might be attributable to concerted action of determinants encoded in pathogenicity islands present in the genome of PA14, a global analysis of the differential host responses to these <italic>P. aeruginosa </italic>strains has not been addressed. Little is known about the host response to infection with <italic>P. aeruginosa </italic>and whether or not the global host transcription is being affected as a defense mechanism or altered in the benefit of the pathogen. Since the social amoeba <italic>Dictyostelium discoideum </italic>is a suitable host to study virulence of <italic>P. aeruginosa </italic>and other pathogens, we used available genomic tools in this model system to study the transcriptional host response to <italic>P. aeruginosa </italic>infection.</p> <title>Results</title> <p>We have compared the virulence of the <italic>P. aeruginosa </italic>PAO1 and PA14 using <italic>D. discoideum </italic>and studied the transcriptional response of the amoeba upon infection. Our results showed that PA14 is more virulent in <italic>Dictyostelium </italic>than PA01using different plating assays. For studying the differential response of the host to infection by these model strains, <italic>D. discoideum </italic>cells were exposed to either <italic>P. aeruginosa </italic>PAO1 or <italic>P. aeruginosa </italic>PA14 (mixed with an excess of the non-pathogenic bacterium <italic>Klebsiella aerogenes </italic>as food supply) and after 4 hours, cellular RNA extracted. A three-way comparison was made using whole-genome <italic>D. discoideum </italic>microarrays between RNA samples from cells treated with the two different strains and control cells exposed only to <italic>K. aerogenes</italic>. The transcriptomic analyses have shown the existence of common and specific responses to infection. The expression of 364 genes changed in a similar way upon infection with one or another strain, whereas 169 genes were differentially regulated depending on whether the infecting strain was either <italic>P. aeruginosa </italic>PAO1 or PA14. Effects on metabolism, signalling, stress response and cell cycle can be inferred from the genes affected.</p> <title>Conclusion</title> <p>Our results show that pathogenic <italic>Pseudomonas </italic>strains invoke both a common transcriptional response from <italic>Dictyostelium </italic>and a strain specific one, indicating that the infective process of bacterial pathogens can be strain-specific and is more complex than previously thought.</p>
<title>Authors' contributions</title> <p>SC and JC performed the biological experiments and analyzed the data. GB, JS and AI designed and constructed the microarray, and provided informatics and analysis tools; GB contributed to the array experimental design, carried out the array experiments, and contributed to the analysis of results; RK was array project manager. JM contributed to the experimental design of the array and biological experiments. RE designed, coordinated the experiments and drafted the manuscript. All the authors read and approved the final manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>This work was supported by grants BMC2006-00394 from the Spanish Ministerio de Educación y Ciencia (R.E.), CCG07-CSIC/SAL-1959 from Comunidad de Madrid/CSIC (R.E.) and BIO2005-04278, LSHM-CT-2005-518152 and LSHM-CT-2005-018705 to J.M. We would like to thank Theresa Feltwell for technical assistance with the microarray experiment. GB, JS and AI were supported by the Wellcome Trust. JC is recipient of a pre-doctoral fellowship from Spanish Ministerio de Educación y Ciencia.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>PA14 is more virulent than PAO1 in SM-plating assay</bold>. Approximately 100 <italic>D. discoideum </italic>cells were cultivated in SM-plates with the indicated proportion of <italic>Klebsiella </italic>and <italic>Pseudomonas </italic>strains (PAO1 or PA14) previously grown and adjusted to the same optical density. Plates were maintained at 22°C for 5 days. Growth of <italic>D. discoideum </italic>is severely affected by the presence of <italic>Pseudomonas </italic>but the inhibition is stronger when PA14 is used.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>PA14 is more virulent than PAO1 in PDF-agar plating assay</bold>. <italic>D. discoideum </italic>cells were cultivated in non-nutrient agar on a lawn of <italic>Klebsiella </italic>and <italic>Pseudomonas </italic>(PAO1 and PA14) at the indicated proportion. Under these conditions PA14 maintain a high virulence as seen by the strong inhibition of <italic>D. discoideum </italic>growth. Heat inactivated PA14 is used as a control.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold><italic>Pseudomonas aeruginosa </italic>induces gene expression changes in <italic>Dictyostelium</italic></bold>. Heat map comparing the genes significantly altered (p &lt; 0.05) between PAO1-treated cells versus control (975) and PA14-treated cells versus control (838). Each row of the plot is a gene and was colored according to the log2ratio of expression with red meaning up-regulation in relation to the controls and blue downregulation. The histogram shows the range of changes in a log2 scale. The data presented are for the three independent experiments combined. The heatmap was generated using the heatmap.2 function of the gplots package in R [##UREF##4##47##]. The dendrogram was generated using Euclidean distance and the "complete" agglomeration method.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Correlation of microaray and real-time PCR</bold>. Real-time PCR measurements of the mRNA levels for seven representative genes whose expression were affected in the array. Upper panel shows a direct comparison of the changes in a log2 scale for PAO1 versus control and the lower panel shows the same genes for PA14 versus control. Blue bars corresponded to quantitative real time PCR and the purple bars to the array data. The array data and the real time PCR displayed are the combination of three independent biological experiments. The correlation coefficients were: R<sup>2 </sup>= 0.87 for PAO1 and R<sup>2 </sup>= 0.91 for PA14.</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p><bold>Functional categories of the genes affected similarly by the exposure with PAO1 and PA14</bold>. The genes whose expression was altered by PAO1 and PA14 were manually annotated (see Table 2) and grouped in functional categories. The size of the blue bars indicates the number of genes upregulated in each category related to the control and purple bars the number of genes downregulated.</p></caption></fig> <fig position="float" id="F6"><label>Figure 6</label><caption><p><bold>Functional categories of the genes differentially affected by the exposure to PAO1 and PA14</bold>. The genes whose expression was differentially altered in PAO1 versus PA14 were manually annotated (see Table 3) and grouped in functional categories. The size of the blue bars indicates the number of genes upregulated in PAO1 versus PA14 in each category and purple bars the number of genes downregulated in PAO1 versus PA14.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Differential genes at p &lt; 0.05 and different Log 2 ratios</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td align="left">PAO1 vs Control</td><td align="left">PA14 vs Control</td><td align="left">PAO1+PA14 vs Control</td><td align="left">PAO1 vs PA14</td></tr></thead><tbody><tr><td align="left">Log 2 ratio</td><td align="left"><bold>752</bold></td><td align="left"><bold>623</bold></td><td align="left"><bold>364 (Table 2)</bold></td><td align="left"><bold>169 (Table 3)</bold></td></tr><tr><td align="left">(&gt;+0.5 or &lt;-0.5)</td><td align="left">461 down 291 up</td><td align="left">396 down 227 up</td><td align="left">258 down 106 up</td><td align="left">60 down 109 up</td></tr><tr><td align="left">Log 2 ratio</td><td align="left"><bold>150</bold></td><td align="left"><bold>125</bold></td><td align="left"><bold>70</bold></td><td align="left"><bold>35</bold></td></tr><tr><td align="left">(&gt;+1 or &lt;-1)</td><td align="left">126 down 24 up</td><td align="left">105 down 20 up</td><td align="left">66 down 4 up</td><td align="left">14 down 21 up</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Genes differentially expressed upon infection with PAO1 and PA14 versus <italic>Klebsiella</italic></p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Gene ID</td><td align="center">Gene function and name</td><td align="left">Log 2 ratio: PAO1-C</td><td align="left">PA14-C</td></tr></thead><tbody><tr><td align="left"><bold>Cell proliferation</bold></td><td/><td/><td/></tr><tr><td align="left">DDB0216882</td><td align="left">Cyclin-dependent kinase regulatory</td><td align="right">-0,89</td><td align="right">-0,97</td></tr><tr><td align="left">DDB0188449</td><td align="left">cdc40, conserved splicing factor</td><td align="right">0,68</td><td align="right">0,78</td></tr><tr><td align="left">DDB0205486</td><td align="left">CDK family protein kinase</td><td align="right">0,81</td><td align="right">0,89</td></tr><tr><td align="left">DDB0168249</td><td align="left">cdk1, "cyclin-dependent_kinase, p34-cdc2_protein"</td><td align="right">-0,78</td><td align="right">-0,83</td></tr><tr><td align="left">DDB0216532</td><td align="left">cdk10;"putative_CDK_family_protein_kinase</td><td align="right">0,91</td><td align="right">1,08</td></tr><tr><td align="left">DDB0185341</td><td align="left">PP-loop family</td><td align="right">-0,60</td><td align="right">-0,70</td></tr><tr><td align="left">DDB0205486</td><td align="left">putative protein serine/threonine kinase, CDK family protein kinase</td><td align="right">0,68</td><td align="right">0,64</td></tr><tr><td align="left">DDB0218360</td><td align="left">PhoPQ-activated pathogenicity-related protein</td><td align="right">-0,59</td><td align="right">-0,57</td></tr><tr><td align="left" colspan="2"><bold>Cellular biogenesis and organization</bold></td><td/><td/></tr><tr><td align="left">DDB0189693</td><td align="left">copA, coatomer protein complex alpha subunit</td><td align="right">0,52</td><td align="right">0,68</td></tr><tr><td align="left">DDB0216892</td><td align="left">lvsB;"BEACH_domain-containing_protein"</td><td align="right">1,05</td><td align="right">0,91</td></tr><tr><td align="left">DDB0202609</td><td align="left">Transport protein particle (TRAPP)</td><td align="right">-0,77</td><td align="right">-0,61</td></tr><tr><td align="left">DDB0187558</td><td align="left">putative mitochondrial import inner membrane translocase</td><td align="right">-0,65</td><td align="right">-0,65</td></tr><tr><td align="left">DDB0186481</td><td align="left">atg9, APG9, "autophagy_protein_9"</td><td align="right">-0,88</td><td align="right">-0,77</td></tr><tr><td align="left">DDB0217942</td><td align="left">Putative Mpv17/PMP22 family</td><td align="right">-0,66</td><td align="right">-0,93</td></tr><tr><td align="left" colspan="2"><bold>Energy</bold></td><td/><td/></tr><tr><td align="left">DDB0217090</td><td align="left">Isocitrate lyase family</td><td align="right">-0,80</td><td align="right">-0,84</td></tr><tr><td align="left">DDB0192001</td><td align="left">ppkA, "poly_P_kinase, polyphosphate_kinase"</td><td align="right">-0,65</td><td align="right">-0,51</td></tr><tr><td align="left">DDB0190821</td><td align="left">sdhB, complex II, iron-sulfur protein (IP) subunit</td><td align="right">0,74</td><td align="right">0,53</td></tr><tr><td align="left">DDB0190821</td><td align="left">sdhB;"complex_II,(ubiquinone), succinic_dehydrogenase"</td><td align="right">0,59</td><td align="right">0,61</td></tr><tr><td align="left">DDB0167662</td><td align="left">similar to Coenzyme Q9</td><td align="right">-0,71</td><td align="right">-0,72</td></tr><tr><td align="left">DDB0204006</td><td align="left">AMPK beta-2 chain</td><td align="right">0,52</td><td align="right">0,73</td></tr><tr><td align="left" colspan="2"><bold>Metabolism</bold></td><td/><td/></tr><tr><td align="left">DDB0187528</td><td align="left">cysteine dioxygenase</td><td align="right">0,69</td><td align="right">0,66</td></tr><tr><td align="left">DDB0185702</td><td align="left">hgd;"homogentisate_1,2-dioxygenase"</td><td align="right">1,05</td><td align="right">1,10</td></tr><tr><td align="left">DDB0184361</td><td align="left">Thiamine pyrophosphate enzyme</td><td align="right">-0,57</td><td align="right">-0,88</td></tr><tr><td align="left">DDB0186120</td><td align="left">gshA, "gamma_glutamylcysteine_synthetase, glutamate-cysteine_ligase"</td><td align="right">-1,72</td><td align="right">-0,94</td></tr><tr><td align="left">DDB0167249</td><td align="left">Aldehyde dehydrogenase</td><td align="right">-0,77</td><td align="right">-0,94</td></tr><tr><td align="left">DDB0192169</td><td align="left">alrA;"aldehyde_reductase, aldo-keto_reductase"</td><td align="right">-0,54</td><td align="right">-0,88</td></tr><tr><td align="left">DDB0218652</td><td align="left">alrB;"aldo-keto_reductase"</td><td align="right">-0,60</td><td align="right">-0,69</td></tr><tr><td align="left">DDB0189745</td><td align="left">alrC;"aldo-keto_reductase"</td><td align="right">-0,84</td><td align="right">-0,72</td></tr><tr><td align="left">DDB0186332</td><td align="left">alrE;"aldo-keto_reductase"</td><td align="right">-1,23</td><td align="right">-0,99</td></tr><tr><td align="left">DDB0204015</td><td align="left">D-Lactate dehydrogenase</td><td align="right">0,92</td><td align="right">0,65</td></tr><tr><td align="left">DDB0187572</td><td align="left">Endoglucanase_E_like</td><td align="right">0,64</td><td align="right">0,78</td></tr><tr><td align="left">DDB0203268</td><td align="left">glgB;"1,4-alpha-glucan_branching_enzyme, branching_enzyme"</td><td align="right">1,02</td><td align="right">0,85</td></tr><tr><td align="left">DDB0187562</td><td align="left">glk;glucokinase</td><td align="right">0,69</td><td align="right">0,56</td></tr><tr><td align="left">DDB0217973</td><td align="left">Gluconolactonase</td><td align="right">0,85</td><td align="right">0,82</td></tr><tr><td align="left">DDB0202855</td><td align="left">Glycoside hydrolase</td><td align="right">-1,63</td><td align="right">-1,45</td></tr><tr><td align="left">DDB0202233</td><td align="left">Glycosyl hydrolase family 7</td><td align="right">-1,26</td><td align="right">-1,38</td></tr><tr><td align="left">DDB0167594</td><td align="left">Glycosyl hydrolases family</td><td align="right">0,60</td><td align="right">0,56</td></tr><tr><td align="left">DDB0204016</td><td align="left">gpt10;"putative_glycophosphotransferase"</td><td align="right">-0,69</td><td align="right">-1,11</td></tr><tr><td align="left">DDB0204037</td><td align="left">Legume lectins beta-chain signature</td><td align="right">-0,86</td><td align="right">-1,16</td></tr><tr><td align="left">DDB0206405</td><td align="left">Mannosyl oligosaccharide glucosidase</td><td align="right">-0,80</td><td align="right">-0,68</td></tr><tr><td align="left">DDB0205896</td><td align="left">NAD-dependent epimerase/dehydratase family protein</td><td align="right">-2,12</td><td align="right">-1,73</td></tr><tr><td align="left">DDB0204752</td><td align="left">Phosphoglycerate mutase family</td><td align="right">0,64</td><td align="right">0,61</td></tr><tr><td align="left">DDB0190464</td><td align="left">Predicted kinase related to galactokinase and mevalonate kinase</td><td align="right">0,52</td><td align="right">1,02</td></tr><tr><td align="left">DDB0186919</td><td align="left">zinc-containing alcohol dehydrogenase</td><td align="right">-2,09</td><td align="right">-1,99</td></tr><tr><td align="left">DDB0168737</td><td align="left">zinc-containing alcohol dehydrogenase (ADH)</td><td align="right">-1,35</td><td align="right">-1,03</td></tr><tr><td align="left">DDB0169356</td><td align="left">carboxylic ester hydrolase</td><td align="right">-1,03</td><td align="right">-0,72</td></tr><tr><td align="left">DDB0190523</td><td align="left">CAS1;"cycloartenol_synthase"</td><td align="right">0,85</td><td align="right">0,88</td></tr><tr><td align="left">DDB0185601</td><td align="left">cutA;"fatty_acid_elongase_3-ketoacyl-CoA_synthase, long_chain_fatty_acid_elongase"</td><td align="right">-0,93</td><td align="right">-0,62</td></tr><tr><td align="left">DDB0188166</td><td align="left">delta-24-sterol methyltransferase</td><td align="right">-1,42</td><td align="right">-1,74</td></tr><tr><td align="left">DDB0186908</td><td align="left">eapA;"alkyl-dihydroxyacetonephosphate_synthase"</td><td align="right">0,72</td><td align="right">0,81</td></tr><tr><td align="left">DDB0187604</td><td align="left">enoyl-CoA hydratase/isomerase domain-containing protei</td><td align="right">-0,65</td><td align="right">-0,63</td></tr><tr><td align="left">DDB0205302</td><td align="left">Enoyl-CoA hydratase/isomerase family</td><td align="right">-0,64</td><td align="right">-0,68</td></tr><tr><td align="left">DDB0205157</td><td align="left">fadB, des5-2, "delta_5_fatty_acid_desaturase"</td><td align="right">0,89</td><td align="right">1,11</td></tr><tr><td align="left">DDB0190288</td><td align="left">fcsB, fatty acyl-CoA synthetase, long-chain-fatty-acid-CoA ligase</td><td align="right">-1,45</td><td align="right">-1,57</td></tr><tr><td align="left">DDB0191679</td><td align="left">GNS1/SUR4 family</td><td align="right">-1,03</td><td align="right">-0,59</td></tr><tr><td align="left">DDB0205505</td><td align="left">mfeB;"hypothetical_peroxisomal_multifunctional_enzyme_2"</td><td align="right">-1,06</td><td align="right">-0,65</td></tr><tr><td align="left">DDB0191653</td><td align="left">patatin family protein</td><td align="right">-2,07</td><td align="right">-1,77</td></tr><tr><td align="left">DDB0218187</td><td align="left">Perilipin family</td><td align="right">-1,35</td><td align="right">-1,04</td></tr><tr><td align="left">DDB0184443</td><td align="left">Saposin (B) Domains</td><td align="right">0,59</td><td align="right">0,59</td></tr><tr><td align="left">DDB0190553</td><td align="left">Similar to sterol-C4-methyl oxidase-like</td><td align="right">-0,71</td><td align="right">-0,91</td></tr><tr><td align="left">DDB0217332</td><td align="left">stearoyl-CoA desaturase</td><td align="right">0,90</td><td align="right">1,02</td></tr><tr><td align="left">DDB0206478</td><td align="left">allC;"allantoate_amidinohydrolase, allantoicase"</td><td align="right">-1,12</td><td align="right">-1,25</td></tr><tr><td align="left">DDB0205700</td><td align="left">cysteine desulfurase</td><td align="right">-1,31</td><td align="right">-1,26</td></tr><tr><td align="left">DDB0169540</td><td align="left">MOSC domain</td><td align="right">-1,00</td><td align="right">-0,67</td></tr><tr><td align="left">DDB0187599</td><td align="left">5 prime nucleotidase family</td><td align="right">-0,90</td><td align="right">-0,70</td></tr><tr><td align="left">DDB0187063</td><td align="left">3-methyl-2-oxobutanoate hydroxymethyltransferase</td><td align="right">-1,83</td><td align="right">-1,41</td></tr><tr><td align="left">DDB0190860</td><td align="left">adenine phosphoribosyltransferase</td><td align="right">-0,98</td><td align="right">-0,75</td></tr><tr><td align="left">DDB0203073</td><td align="left">adenosine deaminase-related growth factor</td><td align="right">-1,35</td><td align="right">-0,74</td></tr><tr><td align="left">DDB0215237</td><td align="left">ATP:D-ribose_5-phosphotransferase, ribokinase</td><td align="right">-0,56</td><td align="right">-0,51</td></tr><tr><td align="left">DDB0206047</td><td align="left">CTP synthase</td><td align="right">-1,01</td><td align="right">-1,20</td></tr><tr><td align="left">DDB0187738</td><td align="left">Cytidine and deoxycytidylate deaminase zinc-binding region</td><td align="right">-0,76</td><td align="right">-0,73</td></tr><tr><td align="left">DDB0191911</td><td align="left">putative RNA methylase</td><td align="right">-1,01</td><td align="right">-1,21</td></tr><tr><td align="left">DDB0185785</td><td align="left">Putative RNA methylase family</td><td align="right">-0,59</td><td align="right">-0,63</td></tr><tr><td align="left">DDB0219236</td><td align="left">pyrK, cytidylate kinase</td><td align="right">-0,98</td><td align="right">-0,71</td></tr><tr><td align="left">DDB0217423</td><td align="left">rnrB_2, ribonucleotide reductase small subunit</td><td align="right">-0,79</td><td align="right">-0,84</td></tr><tr><td align="left">DDB0215284</td><td align="left">tRNA/rRNA methyltransferase SpoU family protein</td><td align="right">-0,77</td><td align="right">-0,93</td></tr><tr><td align="left">DDB0186269</td><td align="left">Thioredoxin family</td><td align="right">-1,41</td><td align="right">-1,91</td></tr><tr><td align="left">DDB0206431</td><td align="left">FAD binding domain</td><td align="right">-1,63</td><td align="right">-1,13</td></tr><tr><td align="left">DDB0187958</td><td align="left">gchA, "GTP_cyclohydrolase_I"</td><td align="right">-0,91</td><td align="right">-0,89</td></tr><tr><td align="left">DDB0185963</td><td align="left">Oxysterol-binding protein</td><td align="right">-0,55</td><td align="right">-0,75</td></tr><tr><td align="left">DDB0205608</td><td align="left">pks18, putative polyketide synthase</td><td align="right">0,54</td><td align="right">0,60</td></tr><tr><td align="left">DDB0168380</td><td align="left">pks5, putative polyketide synthase</td><td align="right">0,57</td><td align="right">0,99</td></tr><tr><td align="left">DDB0219613</td><td align="left">stlB, putative polyketide synthase</td><td align="right">0,68</td><td align="right">1,05</td></tr><tr><td align="left">DDB0186173</td><td align="left">Histidine acid phosphatase</td><td align="right">-0,79</td><td align="right">-0,80</td></tr><tr><td align="left">DDB0184156</td><td align="left">Acetyltransferase (GNAT) family</td><td align="right">0,81</td><td align="right">0,78</td></tr><tr><td align="left">DDB0205937</td><td align="left">Aldehyde dehydrogenase family</td><td align="right">1,02</td><td align="right">0,81</td></tr><tr><td align="left">DDB0183800</td><td align="left">dihydrolipoamide_dehydrogenase</td><td align="right">-0,60</td><td align="right">-0,57</td></tr><tr><td align="left">DDB0188526</td><td align="left">FAD binding domain</td><td align="right">-1,14</td><td align="right">-0,93</td></tr><tr><td align="left">DDB0169374</td><td align="left">haloacid dehalogenase-like hydrolase</td><td align="right">-1,26</td><td align="right">-1,45</td></tr><tr><td align="left">DDB0204714</td><td align="left">hemA, ALAS, "5-aminolevulinate_synthase, ALA_synthase"</td><td align="right">1,00</td><td align="right">0,99</td></tr><tr><td align="left">DDB0187575</td><td align="left">monooxygenase, FAD-binding</td><td align="right">-1,30</td><td align="right">-1,14</td></tr><tr><td align="left">DDB0203608</td><td align="left">NADH:flavin oxidoreductase/NADH oxidase domain-containing protein</td><td align="right">-2,99</td><td align="right">-2,49</td></tr><tr><td align="left">DDB0186877</td><td align="left">Predicted hydrolases or acyltransferases</td><td align="right">-0,70</td><td align="right">-0,83</td></tr><tr><td align="left">DDB0203708</td><td align="left">Putative dehydrogenase domain</td><td align="right">0,60</td><td align="right">0,66</td></tr><tr><td align="left">DDB0186921</td><td align="left">Putative quinone oxidoreductase</td><td align="right">-0,68</td><td align="right">-0,80</td></tr><tr><td align="left">DDB0218378</td><td align="left">selD, selenophosphate synthase</td><td align="right">0,83</td><td align="right">0,85</td></tr><tr><td align="left">DDB0219578</td><td align="left">short chain dehydrogenase</td><td align="right">-0,56</td><td align="right">-0,76</td></tr><tr><td align="left">DDB0168766</td><td align="left">short chain dehydrogenase</td><td align="right">0,75</td><td align="right">0,90</td></tr><tr><td align="left">DDB0201995</td><td align="left">Short-chain alcohol dehydrogenase</td><td align="right">-0,64</td><td align="right">-0,69</td></tr><tr><td align="left">DDB0191047</td><td align="left">Sucrolytic enzyme/ferredoxin homolog protei</td><td align="right">-0,76</td><td align="right">-0,55</td></tr><tr><td align="left">DDB0205223</td><td align="left">Ubiquinone biosynthesis protein</td><td align="right">1,02</td><td align="right">0,99</td></tr><tr><td align="left" colspan="2"><bold>Movement</bold></td><td/><td/></tr><tr><td align="left">DDB0190345</td><td align="left">Actin</td><td align="right">-0,62</td><td align="right">-0,64</td></tr><tr><td align="left">DDB0216677</td><td align="left">tubB;"beta_tubulin"</td><td align="right">0,78</td><td align="right">0,66</td></tr><tr><td align="left" colspan="2"><bold>Protein targeting</bold></td><td/><td/></tr><tr><td align="left">DDB0189280</td><td align="left">CLN3 protein;Major Facilitator Superfamily</td><td align="right">-1,38</td><td align="right">-0,81</td></tr><tr><td align="left">DDB0186130</td><td align="left">pigF, phosphatidylinositol glycan, class Fphosphoethanolamine N- methyltransferase family</td><td align="right">-1,13</td><td align="right">-1,33</td></tr><tr><td align="left">DDB0187271</td><td align="left">Protein prenyltransferase, alpha subunit</td><td align="right">-0,69</td><td align="right">-0,66</td></tr><tr><td align="left">DDB0187195</td><td align="left">Ubiquitin family protein</td><td align="right">-0,68</td><td align="right">-0,56</td></tr><tr><td align="left">DDB0217546</td><td align="left">vps13B, vacuolar protein sorting-associated protein</td><td align="right">0,80</td><td align="right">0,57</td></tr><tr><td align="left">DDB0205767</td><td align="left">Importin-beta N-terminal domain</td><td align="right">-1,39</td><td align="right">-0,91</td></tr><tr><td align="left">DDB0219696</td><td align="left">CSN3, COP9 signalosome complex subunit 3</td><td align="right">0,86</td><td align="right">0,56</td></tr><tr><td align="left">DDB0188097</td><td align="left">mppB;"mitochondrial_processing_peptidase_beta_subunit"</td><td align="right">0,85</td><td align="right">0,93</td></tr><tr><td align="left">DDB0188792</td><td align="left">npl4, nuclear protein localization 4</td><td align="right">-1,46</td><td align="right">-1,26</td></tr><tr><td align="left">DDB0189322</td><td align="left">Peptidase family M41</td><td align="right">0,66</td><td align="right">0,82</td></tr><tr><td align="left">DDB0204548</td><td align="left">putative E3 ubiquitin ligase</td><td align="right">0,56</td><td align="right">0,88</td></tr><tr><td align="left">DDB0203213</td><td align="left">Putative serine protease</td><td align="right">0,57</td><td align="right">0,64</td></tr><tr><td align="left">DDB0190305</td><td align="left">RING-finger-containing ubiquitin ligase</td><td align="right">-1,71</td><td align="right">-2,23</td></tr><tr><td align="left">DDB0191910</td><td align="left">sigB, GP63, orfGP63", "leishmanolysin_family_protein, peptidase</td><td align="right">-1,44</td><td align="right">-1,19</td></tr><tr><td align="left">DDB0219558</td><td align="left">usp12, putative ubiquitin carboxyl-terminal hydrolase (UCH)</td><td align="right">0,81</td><td align="right">0,50</td></tr><tr><td align="left">DDB0188490</td><td align="left">usp40, putative ubiquitin carboxyl-terminal hydrolase (UCH)</td><td align="right">0,83</td><td align="right">0,60</td></tr><tr><td align="left" colspan="2"><bold>Signal transduction</bold></td><td/><td/></tr><tr><td align="left">DDB0167328</td><td align="left">ArfGAP, Arf GTPase activating protein</td><td align="right">0,58</td><td align="right">0,67</td></tr><tr><td align="left">DDB0187828</td><td align="left">gacX, RacGAP</td><td align="right">0,62</td><td align="right">0,52</td></tr><tr><td align="left">DDB0217797</td><td align="left">gpaG, "G-protein_subunit_alpha_7"</td><td align="right">-0,93</td><td align="right">-0,77</td></tr><tr><td align="left">DDB0205484</td><td align="left">GTPase-activator protein for Ras-like</td><td align="right">0,58</td><td align="right">0,69</td></tr><tr><td align="left">DDB0202545</td><td align="left">rabX;"Rab_GTPase"</td><td align="right">-0,59</td><td align="right">-0,76</td></tr><tr><td align="left">DDB0186244</td><td align="left">abkD, AdckB2, "putative_ABC1_family_protein_kinase"</td><td align="right">0,80</td><td align="right">0,88</td></tr><tr><td align="left">DDB0217600</td><td align="left">nek3, putative protein serine/threonine kinase</td><td align="right">0,55</td><td align="right">0,56</td></tr><tr><td align="left">DDB0203684</td><td align="left">tyrosine kinase-like</td><td align="right">0,53</td><td align="right">0,56</td></tr><tr><td align="left">DDB0189806</td><td align="left">vwkA;"protein_serine/threonine_kinase"</td><td align="right">-0,77</td><td align="right">-0,67</td></tr><tr><td align="left">DDB0205355</td><td align="left">Calcineurin-like phosphoesterase</td><td align="right">-1,27</td><td align="right">-0,86</td></tr><tr><td align="left">DDB0218779</td><td align="left">pdsA, "PDE, pde1, pdeA", "cAMP_phosphodiesterase"</td><td align="right">0,79</td><td align="right">1,40</td></tr><tr><td align="left">DDB0204820</td><td align="left">rabS;"Rab_GTPase"</td><td align="right">-0,80</td><td align="right">-0,87</td></tr><tr><td align="left">DDB0190872</td><td align="left">NLI interacting factor-like phosphatase</td><td align="right">-0,55</td><td align="right">-0,60</td></tr><tr><td align="left">DDB0185382</td><td align="left">Protein phosphatase 5, catalytic subunit</td><td align="right">-0,68</td><td align="right">-0,62</td></tr><tr><td align="left">DDB0186390</td><td align="left">protein tyrosine phosphatase</td><td align="right">-0,85</td><td align="right">-0,74</td></tr><tr><td align="left">DDB0218065</td><td align="left">ptpB, DdPTPa, "phosphotyrosine_phosphatase_ptp2</td><td align="right">-1,20</td><td align="right">-1,01</td></tr><tr><td align="left">DDB0189698</td><td align="left">Tyrosine specific protein phosphatases family</td><td align="right">-0,76</td><td align="right">-0,89</td></tr><tr><td align="left">DDB0203756</td><td align="left">G-protein-coupled receptor (GPCR) family protein</td><td align="right">0,69</td><td align="right">0,99</td></tr><tr><td align="left">DDB0189216</td><td align="left">gacD, RacGAP</td><td align="right">-0,63</td><td align="right">-0,58</td></tr><tr><td align="left">DDB0217433</td><td align="left">Regulator of G protein signaling</td><td align="right">0,90</td><td align="right">0,72</td></tr><tr><td align="left">DDB0169375</td><td align="left">cGMP-specific phosphodiesterase</td><td align="right">-0,83</td><td align="right">-0,79</td></tr><tr><td align="left">DDB0167494</td><td align="left">plbF, PLB, "phospholipase_B-like"</td><td align="right">-1,17</td><td align="right">-1,01</td></tr><tr><td align="left">DDB0168860</td><td align="left">sgkA, "SK, SPHK", "sphingosine_kinase"</td><td align="right">-0,97</td><td align="right">-0,79</td></tr><tr><td align="left">DDB0167227</td><td align="left">Cytochrome b5-like Heme/Steroid binding domain</td><td align="right">-0,71</td><td align="right">-1,57</td></tr><tr><td align="left">DDB0216720</td><td align="left">Tetraspanin family</td><td align="right">-0,73</td><td align="right">-0,68</td></tr><tr><td align="left" colspan="2"><bold>Stress response</bold></td><td/><td/></tr><tr><td align="left">DDB0202483</td><td align="left">AhpC/TSA family</td><td align="right">-0,51</td><td align="right">-0,68</td></tr><tr><td align="left">DDB0203727</td><td align="left">AhpC/TSA family protein. Thioredoxin-like</td><td align="right">-1,25</td><td align="right">-1,89</td></tr><tr><td align="left">DDB0203727</td><td align="left">AhpC/TSA family protein. Thioredoxin-like</td><td align="right">-1,11</td><td align="right">-2,00</td></tr><tr><td align="left">DDB0205904</td><td align="left">AhpC/TSA family protein. Thioredoxin-like</td><td align="right">-1,10</td><td align="right">-1,63</td></tr><tr><td align="left">DDB0168230</td><td align="left">Cytochrome P450</td><td align="right">0,61</td><td align="right">0,67</td></tr><tr><td align="left">DDB0217979</td><td align="left">cytochrome P450 family protein</td><td align="right">-1,11</td><td align="right">-1,24</td></tr><tr><td align="left">DDB0187276</td><td align="left">cytochrome P450 family protein</td><td align="right">-0,73</td><td align="right">-0,88</td></tr><tr><td align="left">DDB0186118</td><td align="left">cytochrome P450 family protein</td><td align="right">-0,54</td><td align="right">-0,51</td></tr><tr><td align="left">DDB0167587</td><td align="left">Glutathione S-transferase</td><td align="right">-2,24</td><td align="right">-2,02</td></tr><tr><td align="left">DDB0218804</td><td align="left">Glutathione S-transferase</td><td align="right">-1,07</td><td align="right">-0,86</td></tr><tr><td align="left">DDB0185602</td><td align="left">putative FMN-dependent NAD(P)H:quinone reductase</td><td align="right">-0,66</td><td align="right">-0,93</td></tr><tr><td align="left">DDB0168563</td><td align="left">putative glutathione S-transferase</td><td align="right">-0,64</td><td align="right">-1,23</td></tr><tr><td align="left">DDB0201962</td><td align="left">Ku70-binding family protein</td><td align="right">-0,60</td><td align="right">-0,51</td></tr><tr><td align="left">DDB0191833</td><td align="left">TFIIH4, "TFIIH_subunit, general_transcription_factor_IIH,_polypeptide_4"</td><td align="right">-0,81</td><td align="right">-0,90</td></tr><tr><td align="left">DDB0204089</td><td align="left">NUDIX hydrolase family</td><td align="right">-0,89</td><td align="right">-1,14</td></tr><tr><td align="left">DDB0185428</td><td align="left">Strictosidine synthase</td><td align="right">-1,74</td><td align="right">-1,01</td></tr><tr><td align="left">DDB0169033</td><td align="left">trap1, Dd-trap1, "TNF_receptor-associated_protein", member of the HSP90 fam</td><td align="right">-1,93</td><td align="right">-1,56</td></tr><tr><td align="left">DDB0188234</td><td align="left">Chaperone clpB</td><td align="right">-1,00</td><td align="right">-1,25</td></tr><tr><td align="left">DDB0192088</td><td align="left">heat shock cognate protein</td><td align="right">-0,82</td><td align="right">-0,99</td></tr><tr><td align="left">DDB0192086</td><td align="left">heat shock protein, 70 kDa heat shock protein</td><td align="right">-0,94</td><td align="right">-1,18</td></tr><tr><td align="left">DDB0169051</td><td align="left">Hsp20/alpha crystallin family</td><td align="right">-1,13</td><td align="right">-1,30</td></tr><tr><td align="left">DDB0169044</td><td align="left">Hsp20/alpha crystallin family</td><td align="right">-1,07</td><td align="right">-1,12</td></tr><tr><td align="left">DDB0169207</td><td align="left">hspG12, heat shock protein Hsp20 domain-containing protein</td><td align="right">-0,72</td><td align="right">-0,86</td></tr><tr><td align="left" colspan="2"><bold>Transcription</bold></td><td/><td/></tr><tr><td align="left">DDB0204405</td><td align="left">CRTF;"transcription_factor"</td><td align="right">0,97</td><td align="right">1,15</td></tr><tr><td align="left">DDB0167879</td><td align="left">IWS1 C-terminus</td><td align="right">0,70</td><td align="right">0,69</td></tr><tr><td align="left">DDB0205969</td><td align="left">snd1, tudor domain-containing protein</td><td align="right">0,66</td><td align="right">0,79</td></tr><tr><td align="left">DDB0188840</td><td align="left">TFIIAL, "transcription_factor_IIA"</td><td align="right">-0,79</td><td align="right">-0,55</td></tr><tr><td align="left">DDB0167865</td><td align="left">ddx52, DEAD/DEAH box helicase</td><td align="right">-0,55</td><td align="right">-0,55</td></tr><tr><td align="left">DDB0184074</td><td align="left">ddx6, DEAD/DEAH box helicase</td><td align="right">0,62</td><td align="right">0,66</td></tr><tr><td align="left">DDB0184228</td><td align="left">DEAD/DEAH box helicase</td><td align="right">-0,63</td><td align="right">-0,60</td></tr><tr><td align="left">DDB0206136</td><td align="left">myb domain-containing protein</td><td align="right">0,59</td><td align="right">0,58</td></tr><tr><td align="left">DDB0189583</td><td align="left">rpmA;"DNA-dependent_RNA_polymerase"</td><td align="right">-1,33</td><td align="right">-1,17</td></tr><tr><td align="left">DDB0186101</td><td align="left">pwp2, ortholog of H. sapiens and S. cerevisiae PWP2</td><td align="right">-0,80</td><td align="right">-0,88</td></tr><tr><td align="left">DDB0192008</td><td align="left">rpa2, RNA polymerase I, second largest subunit</td><td align="right">-0,98</td><td align="right">-1,13</td></tr><tr><td align="left">DDB0218008</td><td align="left">rpc4;"putative_RNA_polymerase_III_subunit"</td><td align="right">-0,58</td><td align="right">-0,55</td></tr><tr><td align="left">DDB0216877</td><td align="left">tRNA pseudouridine synthase</td><td align="right">-0,56</td><td align="right">-0,65</td></tr><tr><td align="left">DDB0204724</td><td align="left">DNA helicase TIP49, TBP-interacting protei</td><td align="right">0,63</td><td align="right">0,64</td></tr><tr><td align="left">DDB0219410</td><td align="left">pirin-like protein</td><td align="right">-0,99</td><td align="right">-0,81</td></tr><tr><td align="left" colspan="2"><bold>Translation</bold></td><td/><td/></tr><tr><td align="left">DDB0167043</td><td align="left">Ribosomal protein L10</td><td align="right">-0,86</td><td align="right">-1,37</td></tr><tr><td align="left">DDB0190639</td><td align="left">Fibrillarin</td><td align="right">-0,65</td><td align="right">-0,84</td></tr><tr><td align="left">DDB0184302</td><td align="left">Mitochondrial small ribosomal subunit Rsm22</td><td align="right">-0,70</td><td align="right">-0,72</td></tr><tr><td align="left">DDB0205674</td><td align="left">MPP10, U3 small nucleolar ribonucleoprotein</td><td align="right">-0,77</td><td align="right">-0,66</td></tr><tr><td align="left">DDB0201601</td><td align="left">mrpl11, S60 ribosomal protein L11, mitochondrial</td><td align="right">-0,91</td><td align="right">-1,26</td></tr><tr><td align="left">DDB0204554</td><td align="left">Ribosomal protein L28</td><td align="right">-0,50</td><td align="right">-0,69</td></tr><tr><td align="left">DDB0188692</td><td align="left">Ribosomal protein S8e</td><td align="right">-0,94</td><td align="right">-0,79</td></tr><tr><td align="left">DDB0183814</td><td align="left">Ribosomal RNA processing protein 4</td><td align="right">0,58</td><td align="right">0,56</td></tr><tr><td align="left">DDB0188661</td><td align="left">rps9, "rp1024, v12", "40S_ribosomal_protein_S9</td><td align="right">0,60</td><td align="right">0,70</td></tr><tr><td align="left">DDB0219852</td><td align="left">u3 small nucleolar RNA interacting protein 2, putative</td><td align="right">-0,75</td><td align="right">-0,63</td></tr><tr><td align="left">DDB0191852</td><td align="left">eukaryotic translation initiation factor 3 subunit 5</td><td align="right">0,59</td><td align="right">0,54</td></tr><tr><td align="left">DDB0203843</td><td align="left">Eukaryotic translation initiation factor 6 (EIF-6)-like protein</td><td align="right">-0,96</td><td align="right">-0,64</td></tr><tr><td align="left">DDB0189529</td><td align="left">gfm2, mitochondrial translation elongation factor G</td><td align="right">-0,59</td><td align="right">-0,50</td></tr><tr><td align="left">DDB0202851</td><td align="left">NMD3 family</td><td align="right">-0,95</td><td align="right">-0,97</td></tr><tr><td align="left">DDB0168814</td><td align="left">aspartyl-tRNA_synthetase</td><td align="right">0,70</td><td align="right">0,98</td></tr><tr><td align="left" colspan="2"><bold>Transport</bold></td><td/><td/></tr><tr><td align="left">DDB0217304</td><td align="left">ABC transporter AbcG17</td><td align="right">-1,58</td><td align="right">-1,44</td></tr><tr><td align="left">DDB0167281</td><td align="left">ABC transporter mdrA2</td><td align="right">-2,65</td><td align="right">-1,44</td></tr><tr><td align="left">DDB0167281</td><td align="left">ABC transporter mdrA2</td><td align="right">-2,26</td><td align="right">-1,45</td></tr><tr><td align="left">DDB0191940</td><td align="left">abcB2;"ABC_transporter_B_family_protein"</td><td align="right">0,87</td><td align="right">0,93</td></tr><tr><td align="left">DDB0188931</td><td align="left">abcE1;"RNaseL_inhibitor-like_protein, non-transporter_ABC_protein"</td><td align="right">0,61</td><td align="right">0,74</td></tr><tr><td align="left">DDB0189332</td><td align="left">amino acid permease family protein</td><td align="right">-0,72</td><td align="right">-0,87</td></tr><tr><td align="left">DDB0168564</td><td align="left">Amino acid/polyamine transporter</td><td align="right">-1,28</td><td align="right">-1,06</td></tr><tr><td align="left">DDB0190286</td><td align="left">mcfF, Mitochondrial carrier protein</td><td align="right">-0,99</td><td align="right">-0,86</td></tr><tr><td align="left">DDB0216936</td><td align="left">mftA;"carrier_protein_RIM"</td><td align="right">-0,91</td><td align="right">-1,05</td></tr><tr><td align="left">DDB0188529</td><td align="left">nucleoporin family protein</td><td align="right">-0,94</td><td align="right">-0,72</td></tr><tr><td align="left">DDB0189222</td><td align="left">ccsA, copper chaperone for superoxide dismutase</td><td align="right">-0,69</td><td align="right">-0,62</td></tr><tr><td align="left">DDB0205129</td><td align="left">Co/Zn/Cd efflux system component</td><td align="right">-0,56</td><td align="right">-0,71</td></tr><tr><td align="left">DDB0202441</td><td align="left">nhe1, DdNHE1, "Na-H_exchanger, sodium/hydrogen_exchanger"</td><td align="right">-0,77</td><td align="right">-0,59</td></tr><tr><td align="left">DDB0168533</td><td align="left">porA;porin</td><td align="right">0,68</td><td align="right">0,59</td></tr><tr><td align="left">DDB0218156</td><td align="left">P-type cation-transporting ATPase</td><td align="right">-0,69</td><td align="right">-0,57</td></tr><tr><td align="left">DDB0189480</td><td align="left">mcfT, mitochondrial substrate carrier family protein</td><td align="right">-0,74</td><td align="right">-0,86</td></tr><tr><td align="left">DDB0202337</td><td align="left">Nodulin, Major Facilitator Superfamily</td><td align="right">-0,60</td><td align="right">-0,68</td></tr><tr><td align="left">DDB0185520</td><td align="left">Nucleoside transporter</td><td align="right">-0,53</td><td align="right">-0,56</td></tr><tr><td align="left">DDB0203447</td><td align="left">Sugar (and other) transporter</td><td align="right">1,01</td><td align="right">0,71</td></tr><tr><td align="left">DDB0168979</td><td align="left">Sugar transport proteins signature 1</td><td align="right">-0,78</td><td align="right">-1,28</td></tr><tr><td align="left">DDB0189650</td><td align="left">sodium/potassium-transporting ATPase alpha chain 2</td><td align="right">-1,48</td><td align="right">-0,81</td></tr><tr><td align="left">DDB0190036</td><td align="left">Major Facilitator Superfamily</td><td align="right">-1,13</td><td align="right">-1,05</td></tr><tr><td align="left">DDB0205693</td><td align="left">Major Facilitator Superfamily</td><td align="right">-1,33</td><td align="right">-1,08</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Genes differentially expressed upon infection with PAO1 versus PA14</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Gene ID</td><td align="center">Gene function and name</td><td align="left">Log2 ratio: PAO1-14</td><td align="left">PAO1-C</td><td align="left">PA14-C</td></tr></thead><tbody><tr><td align="left" colspan="2"><bold>Cellular Biogenesis and organization</bold></td><td/><td/><td/></tr><tr><td align="left">DDB0187116</td><td align="left">vps13A, vacuolar protein sorting-associated protein</td><td align="right">0,81</td><td align="right">0,67</td><td align="right">-0,13</td></tr><tr><td align="left">DDB0189855</td><td align="left">vps46, Vacuolar Protein Sorting</td><td align="right">0,83</td><td align="right">0,59</td><td align="right">-0,24</td></tr><tr><td align="left" colspan="2"><bold>Energy</bold></td><td/><td/><td/></tr><tr><td align="left">DDB0204335</td><td align="left">cxgE, "cox7E, coxVIIe", "cytochrome_c_oxidase_subunit_VII_E"</td><td align="right">1, 20</td><td align="right">0,80</td><td align="right">-0,40</td></tr><tr><td align="left" colspan="2"><bold>Metabolism</bold></td><td/><td/><td/></tr><tr><td align="left">DDB0190752</td><td align="left">diaminopimelate epimerase</td><td align="right">1,08</td><td align="right">0,02</td><td align="right">-1,06</td></tr><tr><td align="left">DDB0204319</td><td align="left">Hydroxymethyltransferase</td><td align="right">0,63</td><td align="right">0,43</td><td align="right">-0,20</td></tr><tr><td align="left">DDB0168738</td><td align="left">putative arginine deiminase</td><td align="right">-0,62</td><td align="right">-0,78</td><td align="right">-0,16</td></tr><tr><td align="left">DDB0205389</td><td align="left">acly, ATP citrate lyase</td><td align="right">0,89</td><td align="right">0,75</td><td align="right">-0,15</td></tr><tr><td align="left">DDB0169357</td><td align="left">methylenetetrahydrofolate dehydrogenase</td><td align="right">0,78</td><td align="right">0,21</td><td align="right">-0,57</td></tr><tr><td align="left">DDB0187393</td><td align="left">NAD-dependent epimerase/dehydratase family protein</td><td align="right">-0,91</td><td align="right">-1,18</td><td align="right">-0,28</td></tr><tr><td align="left">DDB0205386</td><td align="left">putative ATP citrate synthase</td><td align="right">0,96</td><td align="right">0,79</td><td align="right">-0,17</td></tr><tr><td align="left">DDB0205339</td><td align="left">rpe;"ribulose_phosphate_3-epimerase"</td><td align="right">0,78</td><td align="right">0,43</td><td align="right">-0,35</td></tr><tr><td align="left">DDB0187942</td><td align="left">Short-chain alcohol dehydrogenase of unknown specificit</td><td align="right">0,70</td><td align="right">0,43</td><td align="right">-0,27</td></tr><tr><td align="left">DDB0187544</td><td align="left">smlA</td><td align="right">1,43</td><td align="right">1,00</td><td align="right">-0,43</td></tr><tr><td align="left">DDB0217455</td><td align="left">zinc-containing alcohol dehydrogenase (ADH)</td><td align="right">-1,21</td><td align="right">-1,16</td><td align="right">0,06</td></tr><tr><td align="left">DDB0217374</td><td align="left">zinc-containing alcohol dehydrogenase (ADH)</td><td align="right">-1,02</td><td align="right">-1,37</td><td align="right">-0,35</td></tr><tr><td align="left">DDB0190948</td><td align="left">acid ceramidase-like protein</td><td align="right">-0,87</td><td align="right">-0,31</td><td align="right">0,56</td></tr><tr><td align="left">DDB0188248</td><td align="left">Acyltransferase</td><td align="right">0,64</td><td align="right">0,16</td><td align="right">-0,48</td></tr><tr><td align="left">DDB0219652</td><td align="left">cinB, "esterase/lipase/thioesterase_domain-containing_protein</td><td align="right">0,92</td><td align="right">1,28</td><td align="right">0,36</td></tr><tr><td align="left">DDB0189754</td><td align="left">esterase/lipase/thioesterase domain-containing protein</td><td align="right">1,07</td><td align="right">1,07</td><td align="right">0,00</td></tr><tr><td align="left">DDB0185740</td><td align="left">esterase/lipase/thioesterase domain-containing protein</td><td align="right">1,08</td><td align="right">1,27</td><td align="right">0,20</td></tr><tr><td align="left">DDB0184141</td><td align="left">Phosphate acyltransferases</td><td align="right">0,55</td><td align="right">-0,26</td><td align="right">-0,82</td></tr><tr><td align="left">DDB0167446</td><td align="left">pks16, putative fatty acid synthase</td><td align="right">1,19</td><td align="right">1,81</td><td align="right">0,63</td></tr><tr><td align="left">DDB0189182</td><td align="left">Putative esterase/lipase/thioesterase</td><td align="right">0,99</td><td align="right">1,43</td><td align="right">0,44</td></tr><tr><td align="left">DDB0191907</td><td align="left">dUTP diphosphatase</td><td align="right">0,61</td><td align="right">0,14</td><td align="right">-0,47</td></tr><tr><td align="left">DDB0217901</td><td align="left">purH, AICAR transformylase/IMP cyclohydrolase</td><td align="right">0,89</td><td align="right">0,34</td><td align="right">-0,55</td></tr><tr><td align="left">DDB0167009</td><td align="left">pyr4;"dihydroorotate_dehydrogenase, dihydroorotate_oxidase"</td><td align="right">0,68</td><td align="right">0,18</td><td align="right">-0,50</td></tr><tr><td align="left">DDB0217842</td><td align="left">rpiA;"ribose-5-phosphate_isomerase"</td><td align="right">0,63</td><td align="right">0,42</td><td align="right">-0,21</td></tr><tr><td align="left">DDB0189571</td><td align="left">Sulfite reductase, alpha subunit (flavoprotein)</td><td align="right">-0,89</td><td align="right">-0,09</td><td align="right">0,79</td></tr><tr><td align="left">DDB0202318</td><td align="left">Cyclopropyl sterol isomerase</td><td align="right">0,73</td><td align="right">0,31</td><td align="right">-0,42</td></tr><tr><td align="left">DDB0167227</td><td align="left">Cytochrome b5-like Heme/Steroid binding domain</td><td align="right">0,86</td><td align="right">-0,71</td><td align="right">-1,57</td></tr><tr><td align="left">DDB0168923</td><td align="left">dihydropteridine reductase</td><td align="right">0,60</td><td align="right">0,14</td><td align="right">-0,47</td></tr><tr><td align="left">DDB0185998</td><td align="left">ERG24, Ergosterol biosynthesis</td><td align="right">0,59</td><td align="right">-0,02</td><td align="right">-0,61</td></tr><tr><td align="left">DDB0219255</td><td align="left">Fe(II) oxygenase superfamily</td><td align="right">1,14</td><td align="right">0,91</td><td align="right">-0,23</td></tr><tr><td align="left">DDB0217308</td><td align="left">Predicted iron-dependent peroxidase</td><td align="right">-0,55</td><td align="right">-0,12</td><td align="right">0,43</td></tr><tr><td align="left">DDB0192180</td><td align="left">putative O-methyltransferase</td><td align="right">0,56</td><td align="right">0,71</td><td align="right">0,15</td></tr><tr><td align="left">DDB0202301</td><td align="left">putative SAM dependent methyltransferase</td><td align="right">0,88</td><td align="right">0,94</td><td align="right">0,05</td></tr><tr><td align="left">DDB0167345</td><td align="left">short-chain dehydrogenase/reductase (SDR) family protein</td><td align="right">0,73</td><td align="right">0,56</td><td align="right">-0,17</td></tr><tr><td align="left" colspan="2"><bold>Movement</bold></td><td/><td/><td/></tr><tr><td align="left">DDB0188280</td><td align="left">myoB, "DMIB, abmB", "myosin_IB"</td><td align="right">-0,85</td><td align="right">-0,30</td><td align="right">0,55</td></tr><tr><td align="left">DDB0167337</td><td align="left">myoD, DMID, "myosin_ID_heavy_chain"</td><td align="right">-0,81</td><td align="right">-0,76</td><td align="right">0,05</td></tr><tr><td align="left" colspan="2"><bold>Multicellular organization</bold></td><td/><td/><td/></tr><tr><td align="left">DDB0216906</td><td align="left">comC;"FIBROSURFIN_PRECURSOR"</td><td align="right">-0,65</td><td align="right">-0,39</td><td align="right">0,26</td></tr><tr><td align="left" colspan="2"><bold>Protein destination</bold></td><td/><td/><td/></tr><tr><td align="left">DDB0189735</td><td align="left">homolog to co-chaperone p23</td><td align="right">0,56</td><td align="right">0,45</td><td align="right">-0,11</td></tr><tr><td align="left">DDB0187409</td><td align="left">Acetyltransferase (GNAT) family</td><td align="right">0,57</td><td align="right">0,33</td><td align="right">-0,24</td></tr><tr><td align="left">DDB0189799</td><td align="left">SET domain-containing protein</td><td align="right">0,52</td><td align="right">0,33</td><td align="right">-0,19</td></tr><tr><td align="left">DDB0202482</td><td align="left">26S proteasome non-ATPase regulatory subunit 9</td><td align="right">0,78</td><td align="right">1,25</td><td align="right">0,47</td></tr><tr><td align="left">DDB0219654</td><td align="left">Cysteine proteinase 1 precursor</td><td align="right">0,70</td><td align="right">0,61</td><td align="right">-0,09</td></tr><tr><td align="left">DDB0167298</td><td align="left">Dipeptidyl aminopeptidase</td><td align="right">0,51</td><td align="right">-0,07</td><td align="right">-0,59</td></tr><tr><td align="left">DDB0190542</td><td align="left">Probable proteasome subunit beta type 3</td><td align="right">0,57</td><td align="right">-0,06</td><td align="right">-0,63</td></tr><tr><td align="left">DDB0216902</td><td align="left">prtA, M3L, "proteosomal_alpha-subunit_M3"</td><td align="right">0,93</td><td align="right">0,94</td><td align="right">0,01</td></tr><tr><td align="left">DDB0216901</td><td align="left">prtB, M3R, "proteosomal_alpha-subunit_7-1"</td><td align="right">1,37</td><td align="right">1,30</td><td align="right">-0,06</td></tr><tr><td align="left">DDB0186869</td><td align="left">small ubiquitin-like protein</td><td align="right">0,61</td><td align="right">0,33</td><td align="right">-0,28</td></tr><tr><td align="left">DDB0217344</td><td align="left">ubiquitin-like domain containing CTD phosphatase</td><td align="right">0,52</td><td align="right">0,29</td><td align="right">-0,23</td></tr><tr><td align="left" colspan="2"><bold>Signal transduction</bold></td><td/><td/><td/></tr><tr><td align="left">DDB0169410</td><td align="left">gpaB, "Ga2, Galpha2, gpa2", "G-protein_subunit_alpha_2"</td><td align="right">-0,82</td><td align="right">-0,55</td><td align="right">0,27</td></tr><tr><td align="left">DDB0190318</td><td align="left">rab1C;"Rab_GTPase"</td><td align="right">0,87</td><td align="right">0,76</td><td align="right">-0,10</td></tr><tr><td align="left">DDB0202066</td><td align="left">pakC, STE20 family protein kinase</td><td align="right">-0,61</td><td align="right">-0,18</td><td align="right">0,43</td></tr><tr><td align="left">DDB0169250</td><td align="left">putative protein serine/threonine kinase</td><td align="right">-0,87</td><td align="right">-0,58</td><td align="right">0,29</td></tr><tr><td align="left">DDB0205782</td><td align="left">roco6;"putative_protein_serine/threonine_kinase</td><td align="right">-0,61</td><td align="right">-0,48</td><td align="right">0,14</td></tr><tr><td align="left">DDB0167076</td><td align="left">sepA, putative protein serine/threonine kinase</td><td align="right">-0,57</td><td align="right">-0,31</td><td align="right">0,26</td></tr><tr><td align="left">DDB0217465</td><td align="left">fslH, G-protein-coupled receptor (GPCR) family protein</td><td align="right">0,71</td><td align="right">0,24</td><td align="right">-0,46</td></tr><tr><td align="left">DDB0205174</td><td align="left">grlF, GABA-B receptor-like protein</td><td align="right">-1,41</td><td align="right">-1,54</td><td align="right">-0,12</td></tr><tr><td align="left">DDB0168770</td><td align="left">grlJ, GABA-B receptor-like protein</td><td align="right">-0,63</td><td align="right">-0,96</td><td align="right">-0,32</td></tr><tr><td align="left">DDB0204083</td><td align="left">grlL, GABA-B receptor-like protein</td><td align="right">-1,84</td><td align="right">-1,25</td><td align="right">0,60</td></tr><tr><td align="left">DDB0229801</td><td align="left">grlQ, G-protein-coupled receptor (GPCR) family protein</td><td align="right">0,77</td><td align="right">0,98</td><td align="right">0,21</td></tr><tr><td align="left">DDB0167432</td><td align="left">gacFF, RacGAP</td><td align="right">-0,57</td><td align="right">0,08</td><td align="right">0,66</td></tr><tr><td align="left">DDB0167384</td><td align="left">putative guanine nucleotide exchange factor (GEF)</td><td align="right">-0,70</td><td align="right">-0,01</td><td align="right">0,70</td></tr><tr><td align="left">DDB0167541</td><td align="left">dpoA;"prolyl_oligopeptidase"</td><td align="right">-0,66</td><td align="right">-0,11</td><td align="right">0,55</td></tr><tr><td align="left" colspan="2"><bold>Stress response and cell rescue</bold></td><td/><td/><td/></tr><tr><td align="left">DDB0218719</td><td align="left">AhpC/TSA family protein</td><td align="right">0,82</td><td align="right">0,76</td><td align="right">-0,05</td></tr><tr><td align="left">DDB0205882</td><td align="left">AhpC/TSA family protein. Thioredoxin-like</td><td align="right">0,75</td><td align="right">0,57</td><td align="right">-0,18</td></tr><tr><td align="left">DDB0217383</td><td align="left">Cytochrome P450</td><td align="right">1,16</td><td align="right">0,42</td><td align="right">-0,74</td></tr><tr><td align="left">DDB0168563</td><td align="left">putative glutathione S-transferase</td><td align="right">0,59</td><td align="right">-0,64</td><td align="right">-1,23</td></tr><tr><td align="left">DDB0167491</td><td align="left">alyA;lysozyme</td><td align="right">0,63</td><td align="right">0,73</td><td align="right">0,10</td></tr><tr><td align="left">DDB0167489</td><td align="left">alyB;lysozyme</td><td align="right">0,69</td><td align="right">0,75</td><td align="right">0,07</td></tr><tr><td align="left">DDB0167490</td><td align="left">alyC;lysozyme</td><td align="right">0,71</td><td align="right">0,87</td><td align="right">0,16</td></tr><tr><td align="left" colspan="2"><bold>Transcription</bold></td><td/><td/><td/></tr><tr><td align="left">DDB0206051</td><td align="left">member of NOD protein family</td><td align="right">1,05</td><td align="right">0,90</td><td align="right">-0,15</td></tr><tr><td align="left">DDB0202276</td><td align="left">srfC;"putative_MADS-box_transcription_factor"</td><td align="right">-0,90</td><td align="right">-0,30</td><td align="right">0,60</td></tr><tr><td align="left">DDB0217613</td><td align="left">wrky1;"putative_WRKY_transcription_factor"</td><td align="right">1,04</td><td align="right">1,11</td><td align="right">0,07</td></tr><tr><td align="left">DDB0167422</td><td align="left">putative histone acetyltransferase</td><td align="right">-0,64</td><td align="right">-0,55</td><td align="right">0,09</td></tr><tr><td align="left" colspan="2"><bold>Translation</bold></td><td/><td/><td/></tr><tr><td align="left">DDB0201621</td><td align="left">mrps2, ribosomal protein S2, mitochondrial</td><td align="right">0,75</td><td align="right">0,49</td><td align="right">-0,27</td></tr><tr><td align="left">DDB0218535</td><td align="left">Eukaryotic elongation factor 1 (EF1) alpha subfamily</td><td align="right">-0,83</td><td align="right">-0,18</td><td align="right">0,65</td></tr><tr><td align="left">DDB0167339</td><td align="left">tRNA-ribosyltransferase</td><td align="right">-0,52</td><td align="right">-0,56</td><td align="right">-0,04</td></tr><tr><td align="left" colspan="2"><bold>Transport facilitation</bold></td><td/><td/><td/></tr><tr><td align="left">DDB0187089</td><td align="left">abcC5;"ABC_transporter_C_family_protein"</td><td align="right">-1,12</td><td align="right">-0,68</td><td align="right">0,44</td></tr><tr><td align="left">DDB0218568</td><td align="left">Copper-transporting P-type ATPase</td><td align="right">-1,14</td><td align="right">-0,60</td><td align="right">0,54</td></tr><tr><td align="left">DDB0204460</td><td align="left">patA, PAT1, "Ca2+-ATPase, P-type_ATPase"</td><td align="right">-0,93</td><td align="right">-0,15</td><td align="right">0,78</td></tr><tr><td align="left">DDB0205031</td><td align="left">potassium channel tetramerization domain-containing protein</td><td align="right">0,92</td><td align="right">0,84</td><td align="right">-0,08</td></tr><tr><td align="left">DDB0217251</td><td align="left">P-type cation-transporting ATPase</td><td align="right">-2,09</td><td align="right">-1,46</td><td align="right">0,63</td></tr><tr><td align="left">DDB0186223</td><td align="left">P-type cation-transporting ATPase</td><td align="right">-0,69</td><td align="right">-0,40</td><td align="right">0,29</td></tr><tr><td align="left">DDB0167361</td><td align="left">vatE;"vacuolar_H+-ATPase_E_subunit"</td><td align="right">0,54</td><td align="right">0,47</td><td align="right">-0,08</td></tr><tr><td align="left">DDB0217699</td><td align="left">phospholipid-translocating P-type ATPase family protein</td><td align="right">-0,67</td><td align="right">-0,63</td><td align="right">0,04</td></tr><tr><td align="left">DDB0187945</td><td align="left">Sugar (and other) transporter</td><td align="right">-0,90</td><td align="right">-0,78</td><td align="right">0,11</td></tr><tr><td align="left">DDB0206579</td><td align="left">mcfQ, mitochondrial substrate carrier family protein</td><td align="right">0,65</td><td align="right">0,78</td><td align="right">0,13</td></tr><tr><td align="left">DDB0192172</td><td align="left">mcfZ, mitochondrial substrate carrier family protein</td><td align="right">0,73</td><td align="right">0,92</td><td align="right">0,18</td></tr><tr><td align="left">DDB0183815</td><td align="left">Mitochondrial carrier protein</td><td align="right">0,55</td><td align="right">0,81</td><td align="right">0,26</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>Microsoft excel document containing all the results and genes contained in the microarray.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p>Microsoft excel document containing the array data for the genes that showed differences in both PAO1 and PA14. The genes were manually annotated and filtered for p &lt; 0.05 and log2 ratio &gt;0.5 or &lt;-0.5. Table ##TAB##1##2## was obtained from these data.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional file 3</title><p>Microsoft excel document containing the array data and manual annotation of the genes that showed differential regulation between PAO1 and PA14. The genes were filtered for p &lt; 0.05 and log2 ratio &gt;0.5 or &lt;-0.5. Table ##TAB##2##3## was obtained from these data.</p></caption></supplementary-material>
<graphic xlink:href="1471-2180-8-109-1"/> <graphic xlink:href="1471-2180-8-109-2"/> <graphic xlink:href="1471-2180-8-109-3"/> <graphic xlink:href="1471-2180-8-109-4"/> <graphic xlink:href="1471-2180-8-109-5"/> <graphic xlink:href="1471-2180-8-109-6"/>
<media xlink:href="1471-2180-8-109-S1.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2180-8-109-S2.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2180-8-109-S3.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Unal", "Steinert"], "given-names": ["C", "M"], "article-title": ["Dictyostelium discoideum as a model to study host-pathogen interactions"], "source": ["Methods in molecular biology (Clifton, NJ)"], "year": ["2006"], "volume": ["346"], "fpage": ["507"], "lpage": ["515"]}, {"surname": ["Smyth"], "given-names": ["GK"], "article-title": ["Linear models and empirical bayes methods for assessing differential expression in microarray experiments"], "source": ["Stat Appl Genet Mol Biol"], "year": ["2004"], "volume": ["3"]}, {"surname": ["Benjamini", "Hochberg"], "given-names": ["Y", "Y"], "article-title": ["Controlling the flase discovery rate: a practical and powerful approach to multiple testing"], "source": ["J R Statist Soc B"], "year": ["1995"], "volume": ["57"], "fpage": ["289"], "lpage": ["300"]}, {"article-title": ["ArrayExpress"]}, {"surname": ["Warnes"], "given-names": ["GR"], "article-title": ["Includes R source code and/or documentation contributed by Ben Bolker and Thomas Lumley. gplots: Various R programming tools for plotting data. R package version 2.3.2"], "year": ["2007"]}]
BMC Microbiol. 2008 Jun 30; 8:109