Previous Article | Next Article ![]()
Infection and Immunity, November 2007, p. 5313-5324, Vol. 75, No. 11
0019-9567/07/$08.00+0 doi:10.1128/IAI.01807-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
,
Department of Microbiology, College of Natural Sciences,1 Department of Pediatrics, John A. Burns School of Medicine, University of Hawaii at Manoa, Honolulu, Hawaii2
Received 13 November 2006/ Returned for modification 23 February 2007/ Accepted 15 August 2007
|
|
|---|
|
|
|---|
Recovery of an HCD P. aeruginosa population from the lungs of CF patients presents a unique opportunity to analyze the transcriptome of this organism in vivo. The primary focus of our approach was to reveal bacterial in vivo gene expression prior to collection of expectorate from a patient. Specifically, we examined bacteria immediately after they were collected from fresh sputum samples in order to decipher their metabolic activities in vivo. This was achieved by isolating bacterial mRNA immediately after a single expectoration from a 42-year-old CF patient (>109 bacteria/ml of sputum). The sample was used for microarray and real-time reverse transcription (RT)-PCR, and the results of gene expression profiling experiments were compared to the results for a clinical isolate pool from the same patient grown in vitro in 1x M9 medium containing citrate. This comparison, performed using Affymetrix P. aeruginosa GeneChips, revealed hundreds of genes that were induced
2-fold (P
0.05) in vivo compared to in vitro growth. Constitutively expressed genes were also observed by comparing the expression profile of the clinical isolate pool to that of PAO1, both grown in 1x M9 medium containing citrate. In this report, our data provide evidence showing that P. aeruginosa as a population exhibited a specific repertoire of highly expressed genes related to virulence, drug resistance, and nutrient utilization in vivo.
|
|
|---|
0.6), cells were harvested for RNA isolation as described below.
Bacterial RNA isolation from sputum and bacterial culture.
Harvested cell cultures described above and 30 ml of a fresh sputum sample (single expectoration with >109 CFU/ml sputum) collected from a 42-year-old CF patient were kept on ice and processed immediately as follows. The approach and method used for sputum collection were approved by the Institutional Review Board Committee on Human Studies at the University of Hawaii at Manoa, and full informed consent was obtained from the patient. The cell pellets from cultures and the fresh sputum sample were treated on ice for 20 min with an equal volume of Sputolysin reagent (Calbiochem, La Jolla, CA) plus Sigma DNase I (10 U/ml of cell suspension) and proteinase K (600 µg/ml of cell suspension) with intermittent mixing to remove extracellular chromosomal DNA and proteins. Samples were then centrifuged at 8,000 x g, and the cell pellets were washed twice with ice-cold sterile double-distilled water (DDW) to lyse eukaryotic cells and remove soluble cell debris in the sputum. For the sputum sample, an aliquot of the washed cells was serially diluted and plated on no-salt Luria-Bertani medium (Teknova, Hollister, CA) and Pseudomonas isolation agar (Difco) to obtain a clinical isolate pool and for determining bacterial counts. The use of a clinical isolate pool, as opposed to a single clinical isolate, eliminated the possibility of any bias towards a particular isolate, and the pool was more representative of the P. aeruginosa population in vivo. Total RNA was isolated from P. aeruginosa using the established protocol of Stephen Lory (http://cfgenomics.unc.edu/protocols_rna_prep.htm), with a few minor modifications. Briefly, cells were harvested at 4°C and then resuspended and sonicated in Trizol reagent (Invitrogen, Carlsbad, CA) and treated with chloroform. An additional phenol-chloroform extraction was performed, and total nucleic acids were precipitated with ethanol. DNA was digested with DNase I (Promega, Madison, WI), and the RNA was phenol-chloroform extracted, ethanol precipitated, and resuspended in diethyl pyrocarbonate (DEPC)-treated DDW. tRNA was removed using QIAgen RNeasy purification columns as recommended by the manufacturer (QIAGEN, Valencia, CA), with an additional on-column DNase I digestion. The final RNA concentrations and purities were determined with a Beckman DU7500 spectrophotometer (A260/A280, between 1.8 and 2.0). The RNA yield for the sputum sample was
55 µg and was sufficient for performing up to three GeneChip analyses and several real-time PCR runs. In our experience, the volume of sputum was not as critical as the P. aeruginosa cell counts in the sputum for successful microarray and real-time RT-PCR experiments. In particular, we ensured that processed sputum had an HCD (
109 CFU/ml of sputum).
cDNA synthesis, labeling, and hybridization and microarray data analyses.
cDNA was synthesized, fragmented, labeled, and processed as recommended by Affymetrix (Affymetrix, Santa Clara, CA). Hybridization to P. aeruginosa GeneChips was performed at the Greenwood Biotechnology Facility, University of Hawaii. Raw data were obtained using the Affymetrix GeneChip Operating System 1.4 software (Affymetrix, Santa Clara, CA). The tab-delimited files were then imported into the GeneSpring 7.0 software (Agilent Technologies, Redwood City, CA) for further analysis. Analysis was done by conducting pairwise comparisons between duplicate or triplicate GeneChips for two conditions at a time (i.e., PAO1 grown in PC versus PAO1 grown in citrate, PAO1 grown in C16:0 versus PAO1 grown in citrate, in vivo sputum versus in vitro-grown clinical pool, in vitro-grown clinical pool versus PAO1 grown in citrate). The list of genes was subjected to statistical analysis using analysis of variance (ANOVA) and including the two-tailed Student t test, and only significant expression data (P
0.05) were kept. The gene list was further analyzed based on fold changes, where only genes showing a change of twofold or greater were kept, and all hypothetical genes were removed from the final list. Fold change values were averaged for nine independent pairwise comparisons (three GeneChips per in vitro condition) with P values of
0.05. However, for in vivo conditions, fold change values were averaged for six independent pairwise comparisons (two GeneChips for in vivo conditions and three GeneChips for in vitro conditions) with P values of
0.05. Average expression levels of a few selected genes (glp genes, mexY, and plcH) that did not pass the statistical stringency test (P
0.05) were analyzed individually by determining the relative expression level for each pairwise comparison, and the average expression levels were calculated (see Fig. 2).
![]() View larger version (32K): [in a new window] |
FIG. 2. Average fold changes for pairwise comparisons of glp genes, plcH, and mexY, which were not detected due to the statistical stringency (P 0.05) imposed on the GeneChip data for glp genes (Table 4) and plcH and mexY (see Table S1 in the supplemental material). (A to E) Pairwise comparisons of three data sets for PAO1 grown in 1x M9 medium containing 0.4% PC (indicated by the numbers 1 to 3 on the x axis) to three data sets for PAO1 grown in 1x M9 medium containing 20 mM citrate (indicated by the letters A to C on the x axis). (F and G) Pairwise comparisons of two data sets for clinical in vivo samples (indicated by the numbers 1 and 2 on the x axis) to three data sets for the clinical isolate pool grown in 1x M9 medium containing 20 mM citrate (indicated by the letters A to C on the x axis). The glycerol uptake facilitator gene glpF (A), the glycerol metabolism regulator gene glpR (B), the glycerol kinase gene glpK (C), the glycerol-3-phosphate transporter gene glpT (D), and the glycerol-3-phosphate dehydrogenase gene glpD (E) all demonstrated induction when PAO1 was grown on PC compared to when PAO1 was grown on citrate, yielding average changes of 4.4-, 2.4-, 7.0-, 3.5-, and 4.0-fold, respectively. The hemolytic phospholipase C precursor gene plcH (F) and the RND multidrug efflux transporter gene mexY (G) were induced more on average in vivo than in the clinical isolate pool grown in 1x M9 containing 20 mM citrate, with average changes of 2.8- and 2.2-fold, respectively. Real-time RT-PCR confirmation data for glpK and glpD are shown in Table 4, and confirmation data for plcH and mexY are shown in Table 3.
|
cDNA synthesis for real-time RT-PCR. Three micrograms of the same purified mRNA for each condition used for microarrays was used for cDNA synthesis. An iScript cDNA synthesis kit was used as recommended by the manufacturer (Bio-Rad, Hercules, CA). A "no reverse-transcriptase" control was included for each sample during cDNA synthesis to ensure that there was no DNA carryover contamination from the RNA isolation. Final volumes were brought to 1,400 µl with DDW (with no DEPC), and 10 µl was used for each real-time PCR as described below.
Bacterial mRNA processing of a sample from a second CF patient for real-time RT-PCR. A sputum sample was obtained from a young adult (an 18-year-old patient). The total volume of sputum was approximately 2 ml, and the P. aeruginosa plate count was 2.3 x 108 CFU/ml. Total RNA was initially isolated as described above, but the subsequent cDNA synthesis procedure was slightly modified as follows. The total RNA on a QIAGEN RNeasy column was eluted in 30 µl of DEPC-treated DDW and protected with 40 U of rRNasin (Promega, Madison, WI). A 2.75-µg aliquot of RNA was treated with 2 U of RQ1 RNase-free DNase I (Promega) for a second time in 2 µl of iScript cDNA synthesis buffer (Bio-Rad) for 10 min at 37°C. The DNase I was inactivated at 70°C for 10 min, and the sample was chilled on ice. Then 14 µl of iScript cDNA synthesis buffer, 4 µl of iScript reverse transcriptase (Bio-Rad), and 23 µl of DEPC-treated DDW were added to the same tube. cDNA was synthesized according to the recommendations of Bio-Rad. The volume of the final product was adjusted to 1,200 µl with DDW (with no DEPC) rather than 1,400 µl because there was less input RNA, and 10 µl was used as a template for each quantitative real-time TaqMan PCR as described below.
Primers and TaqMan probe design. Primers and probes (see Table S3 in the supplemental material) for each gene in the real-time PCRs were designed using Integrated DNA Technologies Primer Quest software (http://www.idtdna.com). Briefly, the amplicon sizes ranged from 67 to 86 bp and the primer melting temperatures were designed for 62°C, with a melting temperature difference of less than 4°C for each primer pair. The probe melting temperatures were designed to be 5 to 10°C higher than the melting temperatures for the corresponding primer pairs. Primer and probe sequences were also subjected to BLAST analysis against the P. aeruginosa PAO1 genome to eliminate the possibility of nonspecific binding. It was particularly ensured that the primers had no significant complementarity at the 3' end and the probe had no significant complementarity at the 5' end to other nonspecific locations on the PAO1 genome.
Real-time PCR. Supermixtures for all reactions were made and aliquoted into subsupermixtures for each gene assayed. Essentially, each real-time PCR mixture (final volume, 25 µl) contained 10 µl of cDNA, 12.5 µl of iQ Supermix (Bio-Rad), 120 nM of each forward and reverse primer, and 12 nM probe. Real-time PCR was performed with an iCycler iQ (Bio-Rad) using the following protocol: denaturation at 95°C for 10 min and then 55 cycles of amplification and quantification at 95°C for 20 s and at 65°C for 45 s. To control for variations between runs, all three housekeeping genes and the various target genes for each individual condition were amplified at the same time on a 96-well plate.
Real-time PCR data analyses.
For analysis we used the method of Peirson et al. (29) as previously reported, which gives more accurate quantitative real-time PCR data. Because normalization by geometric averaging using multiple housekeeping genes has been shown to yield more accurate fold changes than normalization using a single housekeeping gene (47), we opted to use three housekeeping genes for our normalization analyses. Real-time PCR was conducted for each of the 11 target genes and the three housekeeping genes in eight replicates. Real-time PCR fold change values were averaged for eight replicates for each gene and were determined by comparison with the geometric mean of three housekeeping genes. A list of housekeeping genes was determined based on in vivo and in vitro microarray data from over 40 GeneChips for 16 different growth conditions in various growth phases, from mid-log phase to stationary phase (data not shown). These housekeeping genes were analyzed using GeneSpring 7.0 and were found to be expressed consistently across all conditions (fold change, <2; P
0.05). Three housekeeping genes (PA1769, PA1795, and PA1805) which showed the least amount of variability and were expressed consistently across all conditions (as determined using Spearman's correlation coefficient calculations) were selected. Real-time PCR fold changes were calculated using the amplification plot method and the available macro for data analyses of real-time PCR (DART-PCR) (29). Our requirement for real-time PCR is that the efficiencies of each gene (including housekeeping genes) are within 5% from one condition to the next, but similar efficiencies are not necessarily required for the different genes. Accordingly, the average efficiencies of each gene in this study were very similar for the conditions compared (<4.8% efficiency differences), allowing accurate analysis. We recommend not using the standard curve method for clinical samples, since it requires at least three serial dilution points in order to generate a reasonable standard curve. Due to the initial limited amount of input RNA isolated directly from sputum and the possibility of rare transcripts and low transcript abundance within bacterial cells, we were unable to use the dilution-standard curve method to obtain consistent standard curves and efficiencies (<5% difference) for different conditions for many genes. Using the amplification plot analysis method solves this problem, and this method has been demonstrated (29) to yield very good results comparable to those of standard curve analysis and the absolute quantification method, and it is much more accurate than the 
CT and the efficiency correction method.
Assay for constitutive alginate production. Starter cultures of PAO1 and the clinical isolate pool from the 42-year-old CF patient were grown overnight in Pseudomonas isolation broth and then subcultured (1:200 dilution) into no-salt Luria-Bertani medium (Teknova) and incubated at 250 rpm and 37°C. Cell density (OD540) was measured at regular intervals to monitor cell growth and the growth phases, and aliquots were taken for a simultaneous alginate assay. The amount of alginate produced for each strain was normalized to the corresponding OD540. Cells were harvested at early, mid-, and late-log and early stationary phases, and the alginate assays with carbazole reagent were performed as previously described (13, 16). Extracellular alginate production was quantified in triplicate (± standard error of the mean) at various growth phases. The amount of alginate produced by the clinical isolate pool at early stationary phase was defined as 100%, and all other alginate quantities were normalized to this value.
Microarray data accession number. Our microarray data have been deposited in the NCBI Gene Expression Omnibus (GEO) repository under accession number GSE7704.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Summary of genes induced in vivoa
|
This observation prompted our efforts to analyze the GeneChip expression profiles of the in vitro-grown clinical pool to those of prototype strain PAO1, both grown in 1x M9 medium containing citrate. Our results revealed hundreds of genes constitutively expressed compared to the expression in PAO1, where the expression of these genes in PAO1 is more controlled (Table 2; see Table S2 in the supplemental material). As suspected and summarized in Table 2 (see Table S2 in the supplemental material), the clinical isolate pool obtained from the 42-year-old patient had higher constitutive expression levels of several classical virulence factors (muc and alg genes for biofilm synthesis, lipase, phospholipase, and several different proteases, rhl for rhamnolipid hemolysin, hcn for hydrogen cyanide production, and pch for pyochelin). The significance of biofilms in the lungs of CF patients cannot be emphasized enough, since the levels of expression of many alg genes were several-hundred-fold; the results for one of these genes (algK) were confirmed by real-time RT-PCR (Table 3). This constitutive expression of biofilm synthesis genes in the CF patient isolate pool relative to the basal biofilm production in PAO1 translated into greater alginate production at all growth phases (Fig. 1). The lack of virulence gene (plcH, encoding phospholipase C) expression, despite the presence of plcR expression in vivo, was due to larger variations in the microarray data (P > 0.05) for this gene, since PlcR is the accessory protein that is required for PlcH secretion (see Table S2 in the supplemental material) (9). Although not detected by microarrays at high stringency (P
0.05), plcH was expressed in vivo, and this was confirmed by real-time RT-PCR (Table 3 and Fig. 2F). There were nonidentical fold changes in some genes that may be organized in an operon (e.g., bet genes and PA5372 to PA5375; and alg genes and PA3540 to PA3551) (see Tables S1 and S2 in the supplemental material), which are normal and typical of many microarray data sets described previously (41, 48, 50). However, the general trends for genes expressed in the same direction (positive increase) were observed for genes arranged in operons. In a few cases (e.g., mexX, arcC, and arcD) (see below), gene expression was observed both in vitro and in vivo (see Tables S1 and S2 in the supplemental material). We reasoned that these genes were expressed in a gradient fashion, with the expression highest in vivo, lower in the in vitro clinical isolate pool, and lowest or repressed in PAO1. Therefore, Tables 1 and 2 represent a summary of in vivo expressed genes for in vivo survival and virulence; some genes are regulated (Table 1; see Table S1 in the supplemental material), while others have evolved to be constitutive (Table 2; see Table S2 in the supplemental material).
|
View this table: [in a new window] |
TABLE 2. Summary of constitutively expressed genesa
|
|
View this table: [in a new window] |
TABLE 3. Verification of microarray data by real-time RT-PCR for genes required for in vivo PC degradation, drug resistance (mexY), and biofilm biosynthesis (algK) by P. aeruginosa
|
![]() View larger version (30K): [in a new window] |
FIG. 1. Constitutive expression of P. aeruginosa biofilm/alginate biosynthetic genes results in high production of alginate at all growth phases by the clinical isolate pool from the CF patient relative to the production by strain PAO1. Although the growth of the pooled isolates from the CF patient and the growth of PAO1 are comparable, the pooled isolates from the CF patient showed excessive alginate production at all growth phases.
|
In vivo expression of drug resistance genes. In terms of mechanisms of resistance to antibacterial agents, P. aeruginosa is second to none, having several such mechanisms in its huge arsenal, including biofilm production, target modification and overexpression, drug modification, and more than 12 multidrug efflux pumps (1, 18, 22, 33). All possible mechanisms of drug resistance were observed in our data (see Tables S1 and S2 in the supplemental material), with the exception of target modification since our approach cannot detect such changes. At the time of sputum collection, the patient was being treated with tobramycin (an aminoglycoside) by inhalation, tobramycin and cefepime (a fourth-generation cephalosporin) intravenously, and Bactrim (trimethoprim-sulfamethoxazole) orally. Certainly, the highly expressed biofilm biosynthesis genes could present biofilm as the first barrier for drug resistance. The target of tobramycin is at the level of the ribosome to inhibit protein translation, and we observed significant overexpression of 38 different ribosomal proteins (target overexpression) (see Table S1 in the supplemental material). In addition to the induction of two aminoglycoside drug modification proteins, PA1409 (aphA) and PA4119 (aph), we also observed induction of genes for the aminoglycoside efflux pump, MexXY-OprM (19) (Table 3 and Fig. 2G; see Table S1 in the supplemental material). Mutations in mexZ, a negative regulator of mexX and mexY, have been found in the majority of isolates from several CF patients (45). It has been previously demonstrated that MexXY-OprM also effluxes trimethoprim (7). Therefore, expression of this pump may also provide resistance against administered trimethoprim. Increased expression of five folate biosynthetic genes was also observed (cofactor [see Table S1 in the supplemental material]), likely due to the inhibition of folate synthesis by the trimethoprim-sulfamethoxazole, and one of the proteins overexpressed was dihydropteroate synthase (FolP, PA4750), which is the target of sulfamethoxazole (target overexpression). For potential resistance to the ß-lactam cefepime, overexpression of ß-lactamases was observed in vivo (ampC, 87.9-fold [see Table S1 in the supplemental material]; PA5514, 13.9-fold [see Table S2 in the supplemental material]). A previous study showed that antibiotic resistance mechanisms arise from hypermutable mutS variants (44), and isolates from CF patients do harbor mutS mutations (45). Some of the antibiotic resistance mechanisms discussed here in the population may have arisen from such mutS variants. Regardless, our important findings suggest that P. aeruginosa utilizes several mechanisms in vivo to deal with antibacterial agents in a natural infection, by inducing many genes to counteract the negative effects of the administered drugs (see Tables S1 and S2 in the supplemental material).
General metabolism. Various central pathways expressed by P. aeruginosa (Tables 1 and 2) indicate active metabolism in vivo. Gene expression in vivo for central pathways, including those of saccharide metabolism, the Krebs and glyoxylate cycles, the pentose phosphate pathway, and oxidative phosphorylation, along with vitamin/cofactor biosyntheses and nutrient transport, indicates that P. aeruginosa is very metabolically active in vivo. The carbon sources utilized by P. aeruginosa seem to be both lipids of lung surfactant and amino acids, as discussed below. All these activities may contribute to replication and division in vivo, as evident from the expression of genes for cell wall metabolism and division. The data also indicate that P. aeruginosa actively senses and interacts with the in vivo environment by expressing many genes involved in two-component systems, secretion, and transport. In addition, various processes are controlled by induction of a large number of transcriptional regulators and sigma factors. Finally, in an animal model, the expression of the napEFDABC and arcDABC gene clusters has been monitored to indicate anaerobic respiration by reflecting an increase in nitrate and arginine turnover (4). Our data demonstrate that there is significant expression of one of these two gene clusters, arcDABC (13.1-, 28.0-, 8.8-, and 4.2-fold changes, respectively) (see Table S1 in the supplemental material). Together with expression of genes involved in oxidative respiration, this suggests that anaerobic respiration also occurs in the population, which may result from the diverse microenvironment within the lung or biofilm (21). As alluded to by Nguyen and Singh (21), we believe that in a sputum biofilm in a CF patient there may be subpopulations of P. aeruginosa undergoing oxidative respiration on the biofilm surface, while other subpopulations are anaerobically respiring within the biofilm or oxygen-deficient areas, as indicated by the expression of both aerobic and anaerobic respiration genes in the population.
Active amino acid degradation contributing to HCD replication in vivo. Our data indicate that P. aeruginosa metabolizes amino acids as a nutrient source in vivo. Previous studies utilizing sterilized mucopurulent respiratory liquid or sputum from CF patients to grow P. aeruginosa in vitro have discovered several sputum-induced genes (25, 53). One of these in vitro studies (25), in which P. aeruginosa was grown on a pool of several lyophilized sputum samples, showed that there was induction of a subset of 20 genes for amino acid degradation, suggesting that P. aeruginosa metabolizes amino acids as nutrient sources in vivo. However, this hypothesis requires in vivo confirmation. Our in vivo data (see Tables S1 and S2 in the supplemental material), obtained by isolating mRNA directly from a single expectoration source of sputum, indicate that there is induction of a much larger set of genes for amino acid transport and degradation (36 genes induced and 26 genes constitutively expressed). This set represents a total of 60 different genes (since arcC and arcD are present in both Table S1 and Table S2 in the supplemental material) that are involved in amino acid degradation, which suggests that amino acid degradation in vivo could be much more important than previously thought (25). Since proteases were expressed in vivo (LasA, LasB, and others) (see Table S2 in the supplemental material) to cleave proteins into peptides and amino acids, it was significant to find that many amino acid and peptide transporters were also expressed in our data. However, these data indicate that P. aeruginosa may utilize more than amino acids and that lung surfactant lipids may also serve as a great nutrient source in vivo.
Degradation of lung surfactant lipids by P. aeruginosa.
Lung surfactant components, especially lipids, are essential for proper lung functions, and the lack of these components can lead to respiratory distress syndromes in adults and prematurely born infants. Pulmonary surfactant consists of
10% surfactant proteins (SP-A, SP-B, SP-C, and SP-D) and
90% lipids (PC, phosphatidylglycerol, phosphatidylethanolamine, phosphatidylinositol, phosphatidylserine, and small amounts of other lipids and free fatty acids), and PC accounts for
80% of the lipids (2, 14). Thus, the most plentiful sources of potential bacterial nutrients in lung surfactant are lipids, especially PC.
We initially performed GeneChip studies to determine which genes were induced by P. aeruginosa strain PAO1 during growth on PC in vitro as the sole carbon source (Table 4). By identifying the genes involved in PC degradation in vitro, we could relate and search for the expression of the same set of genes in vivo to determine if P. aeruginosa metabolizes PC in vivo. Through the action of lipases (Fig. 3), the richest nutrient obtained from the lung surfactant PC molecule in vivo is fatty acids, mostly palmitic acid (C16:0) (50 to 60%) but also approximately 10 to 20% each C14:0, C16:1, C18:1, and C18:2 (34). Accordingly, the PC utilized in our experiments (Sigma, St. Louis, MO) contained mostly long-chain fatty acids (LCFA) (C16:0, 33%; C18:1, 30%; C18:2, 14%; C18:0, 14%). These LCFA can be metabolized via ß-oxidation in the fatty acid degradation (Fad) pathway (Fig. 3C). As observed for the relevant genes in Table 4, several key genes signal the degradation of PC and LCFA in vitro. Specifically, lipases (LipA and LipC) and phospholipases (PlcH and the accessory protein PlcR) were expressed to break down PC into smaller constituents: fatty acid, glycerol, and phosphorylcholine. These constituents could be further metabolized in vitro, via the expression of Bet enzymes (choline head group metabolism) (49), Glp enzymes (glycerol metabolism) (42), and Fad enzymes (fatty acid degradation) (8). Fold changes for the glycerol degradation genes in Table 4 were determined with less stringency, because of the larger variations (P > 0.05) between the pairwise comparisons in our microarray data for these genes, and all showed positive induction (Fig. 2). Expression of the fadBA5 operon was observed only with C16:0 and not with PC relative to citrate (Table 4), although these genes were involved in PC degradation in vitro (Fig. 4) and were expressed in vivo (Table 3). FadB5 has 73% homology to the well-characterized Escherichia coli FadB protein (35). To demonstrate that FabBA5 is involved in PC degradation, we created a
fadBA5 deletion mutant. This mutant showed a decrease in the ability to grow on LCFA and PC (Fig. 4). The
fadBA5 mutant could still grow on LCFA, because there are at least two other fadBA operons in P. aeruginosa (unpublished data). Of course, there was less of a defect for this mutant to grow on PC than for it to grow on palmitate, since the PC molecule still contains two other nutrient sources (i.e., glycerol and choline) (Fig. 3).
|
View this table: [in a new window] |
TABLE 4. In vitro microarray analysis of relevant genes expressed twofold or more in P. aeruginosa strain PAO1grown on 1x M9 medium containing 0.4% PC or 0.4% palmitate (C16:0) compared to the expression with 20 mM citrate
|
![]() View larger version (43K): [in a new window] |
FIG. 3. (A) Lung surfactant, made of 90% lipids and 10% proteins, coats the trachea, bronchioles, and alveoli of the lung. (B) Of the 90% lipids, 80% consists of PC, which can be cleaved by P. aeruginosa lipases and phospholipase C into three constituents, including fatty acids (FA), glycerol, and phosphorylcholine. (C) Glycerol and phosphorylcholine constituents are further metabolized by known Glp and Bet enzymes of P. aeruginosa (42, 49). The regulators GlpR and BetR control the expression of these enzymes. GlpT and GlpF of the cytoplasmic membrane facilitate glycerol-3-phosphate and glycerol transport, respectively. Based on the E. coli model (8), predicted steps in the Fad pathway of P. aeruginosa, which has not been characterized, are shown. In addition, the regulation of the fad genes in P. aeruginosa is an enigma. FadL, an outer membrane protein, is involved in fatty acid transport. FAD, flavin adenine dinucleotide; FADH2, reduced flavin adenine dinucleotide; CoA, coenzyme A.
|
![]() View larger version (15K): [in a new window] |
FIG. 4. fadBA5 operon is involved in PC degradation. (A) P. aeruginosa fusion strain PAO1-PfadBA5-lacZ was grown in LB, PC, citrate, and palmitate (C16:0), and the ß-galactosidase (ß-Gal) assay was performed (20). The results indicate the reason for no apparent difference in the fold change in the fadBA5 operon expression shown in Table 4 for growth on PC relative to growth on citrate. The lack of fadBA5 induction when P. aeruginosa strain PAO1 was grown on 1x M9 medium plus PC compared to growth on citrate seems to be due to the high expression of this operon when bacteria were grown in 1x M9 medium plus citrate compared to the expression during growth on PC. (B and C) fadBA5 mutant showed a defect in growth on palmitate compared to wild-type strain PAO1 (WT) (B), and this mutant had a defect in PC degradation compared to wild-type strain PAO1 (C), although the defect was not as dramatic as the defect with palmitate, probably due to growth on other carbon sources (glycerol and phosphorylcholine of PC). Both panel B and panel C show that the fadBA5 mutant had a lower growth rate and overall lower final cell density, which were due to a partial defect in the ability to degrade fatty acid as one of the components of PC as a nutrient. Data for construction of the mutant and fusion are not shown.
|
Expression of PC degradation genes in a second CF patient as determined by real-time RT-PCR. To further demonstrate the importance of potential PC degradation in vivo through gene expression, we analyzed the expression of lipA, plcH, fadB5, fadD1, fadD2, fadD4, betA, glpD, and glpK from P. aeruginosa isolated directly from an 18-year-old CF patient by real-time RT-PCR. However, not armed with microarray data for this patient, we wanted to perform more thorough and careful analyses to clearly show gene expression as it relates to PC degradation. Table 5 indicates that genes required for PC degradation in vitro were also expressed in vivo, compared to both the same clinical isolate pool or PAO1 grown in vitro in 1x M9 medium with 20 mM citrate. As Table 5 shows, using gene expression in PAO1 as a reference is essential to get a baseline for quantifying the relative expression of many genes. In this second CF patient, as in the first patient described above, P. aeruginosa clearly expressed genes essential for the degradation of PC. It is interesting that glp genes, fadD2, and fadD4 in the first patient were deregulated (Table 3). Likewise, these same genes were also constitutively expressed in the second patient (Table 5). These gene expression data further support our overall hypothesis that P. aeruginosa may utilize PC as one of the major nutrient sources in vivo.
|
View this table: [in a new window] |
TABLE 5. Expression of genes involved in PC degradation in a second CF patient by real-time RT-PCR
|
In summary, our study revealed several important aspects of P. aeruginosa pathogenesis, drug resistance, and metabolism during infection in CF patients in relation to population diversity. P. aeruginosa has a huge arsenal of metabolic capabilities, which allows it to exploit many possible nutrients in the lung environment. This is the first study that defines the metabolic capability of P. aeruginosa in vivo and the host pulmonary nutrient factors that may contribute to bacterial replication. Our objective in this study was to identify the expression of metabolic pathways which may contribute to nutrient acquisition, leading to replication and maintenance of P. aeruginosa in the lung at an HCD. Our data suggested that P. aeruginosa degrades amino acids as a nutrient source in the lung. In addition, P. aeruginosa also induces genes in vivo to potentially metabolize lung surfactant lipids, adapting to the many possible nutrient sources.
The gene expression of P. aeruginosa in vivo may be influenced by years of evolution during chronic lung infection and by the microenvironment of each bacterium, which leads to a diverse infectious population. Our study revealed a different level of evolution (i.e., deregulation) for several genes in the transcriptome rather than the deletion or acquisition of genes previously demonstrated in the genome (12, 54). We suggest that the population diversity and microenvironment might influence an individual bacterium's behavior, perhaps contributing to a cooperative and diverse infectious population. The shift in thinking of the bacteria in a chronic infection in a CF patient as a diverse infectious population and not as individual isolates may yield novel insight into future treatment. For example, the existence of hypermutable mutS variants coupled to the bacterial population diversity found in the lungs of CF patients with chronic lung infections as opposed to acute infections has important ramifications for future testing for effective antibiotic combinations and concentrations. Thus, future tests of antibiotics on planktonic cells or biofilm would be performed more effectively with a clinical population or a pool of isolates than with a single isolate, in order to obtain the true resistance potential of the infectious P. aeruginosa population. Therefore, both the diversity of the population and the microenvironment may contribute to the full potential of P. aeruginosa within the lungs of CF patients, such that each member of the population plays a distinct role in order to contribute to the pathogenesis of the infectious population. However, this view awaits novel technology to dissect the metabolism and virulence of individual bacteria in the population relative to the microenvironment, which should answer some questions concerning the extent of bacterial cooperation.
The results of this study should serve as an initial model for understanding (i) how P. aeruginosa reacts to antibacterial treatment in vivo, (ii) nutrient acquisition for extracellular pathogens of the lung, and (iii) the detrimental effects of nutrient metabolism in the lung. Finally, this study opens opportunities for further and more thorough analyses of these important pathogenesis aspects with clinical samples from other CF patients in order to determine the transcriptome conservation and heterogeneity of P. aeruginosa populations in vivo.
The contents of this paper are solely the responsibility of the authors and do not necessarily represent the official views of NCRR or NIH.
We are very grateful to E. P Greenberg and M. J. Schurr for their microarray data, which allowed us to search for appropriate housekeeping genes with our data. We thank Joseph Lam and Chad Walton for critical reading of the manuscript.
Published ahead of print on 27 August 2007. ![]()
Supplemental material for this article may be found at http://iai.asm.org/. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»