A Putative ABC Transporter,HatABCDE, Is among Molecular Determinants of Pyomelanin Production in Pseudomonas aeruginosa

Pyomelanin overproduction is a common phenotype among Pseudomonas aeruginosa isolates recovered from cystic fibrosis and urinary tract infections. Its prevalence suggests that it contributes to the persistence of the producing microbial community, yet little is known about the mechanisms of its production. Using transposon mutagenesis, we identified factors that contribute to melanogenesis in a clinical isolate of P. aeruginosa. In addition to two enzymes already known to be involved in its biosynthesis (homogentisate dioxygenase and hydroxyphenylpyruvate dioxygenase), we identified 26 genes that encode regulatory, metabolic, transport, and hypothetical proteins that contribute to the production of homogentisic acid (HGA), the monomeric precursor of pyomelanin. One of these, PA14_57880, was independently identified four times and is predicted to encode the ATP-binding cassette of an ABC transporter homologous to proteins in Pseudomonas putida responsible for the extrusion of organic solvents from the cytosol. Quantification of HGA production by P. aeruginosa PA14 strains missing the predicted subcomponents of this transporter confirmed its role in HGA production: mutants unable to produce the ATP-binding cassette (PA14_57880) or the permease (PA14_57870) produced substantially less extracellular HGA after growth for 20 h than the parental strain. In these mutants, concurrent accumulation of intracellular HGA was observed. In addition, quantitative real-time PCR revealed that intracellular accumulation of HGA elicits upregulation of these transport genes. Based on their involvement in homogentisic acid transport, we rename the genes of this operon hatABCDE.Pseudomonas aeruginosa is a metabolically versatile, opportunistic pathogen that is a major cause of life-threatening infections in patients with burn wounds, compromised immunity, chronic obstructive pulmonary disease (COPD), and cystic fibrosis (CF) (23, 41). A major contributor to P. aeruginosa''s pathogenicity is its remarkable genomic plasticity, which often results is a wide range of phenotypic variation among isolates obtained from both acute and chronic infections. These phenotypes include small colony variant formation (24), alginate overproduction (36), hyperpigmentation (22), autoaggregation (13), and autolysis (64). Many of these phenotypes evolve as infections progress, and most have been ascribed to “loss-of-function” genome diversification that promotes long-term survival in the host environment (54). In this regard, recent studies have stimulated interest in another example of a loss-of-function phenotype, the mutation or deletion of hmgA, which encodes the homogentisate 1,2-dioxygenase enzyme. The absence of this protein leads to the accumulation and subsequent export of homogentisic acid (HGA), which ultimately aggregates into the pyomelanin polymer that manifests as a reddish brown pigmentation of P. aeruginosa colonies and their surrounding milieu (Fig. ​(Fig.1A)1A) (5, 49).Open in a separate windowFIG. 1.Pyomelanin production by the PA14 ΔhmgA and DKN343 strains. (A) Homogentisate pathway for the catabolism of chorismate and aromatic amino acids. Enzyme names are shown above the arrows for each step. A mutation or deletion of the hmgA gene (encoding homogentisate 1,2-dioxygenase) leads to the accumulation of pyomelanin. (B) Pyomelanin overproduction by the PA14 ΔhmgA mutant is abolished when complemented with an intact hmgA gene. Complementation of a melanogenic clinical P. aeruginosa isolate, DKN343, with hmgA results in no phenotypic change, indicating that other factors contribute to its pigmentation.Production of pyomelanin (and other forms of melanin) has been described to occur in a wide range of bacterial species, including Aeromonas (4), Azotobacter (51), Azospirillum (50), Bacillus (3), Legionella (8), Marinomonas (33), Micrococcus (40), Mycobacterium (45), Proteus (1), Rhizobium (12), Shewanella (61), Sinorhizobium (38), Streptomyces (67), and Vibrio (63) species. Notably, isolates of other bacterial species associated with chronic infections of the CF lung, Burkholderia cenocepacia and Stenotrophomonas maltophilia, can also be melanogenic (28, 58), suggesting a possible role for this pigment in the establishment and/or persistence of infection. Some genera produce melanin under normal conditions via polyphenol oxidases or laccases, while others synthesize the pigment only in response to specific environmental conditions (17, 35). Many species, however, including P. aeruginosa, overproduce pyomelanin as a result of a point mutation in hmgA or large chromosomal deletions of loci containing the homogentisate operon (2, 19). While these genetic variations have been frequently reported, there is little understanding of the competitive advantage, if any, that this pigment confers to the producing bacterium.Proposed roles for pyomelanin include the enhancement of bacterial surface attachment (20), extracellular electron transfer (61), iron reduction/acquisition (8), induction of virulence factor expression (63), heavy metal binding (21), and protection from environmental stress (11, 28, 32, 44, 53, 65). A protective role has also been proposed to occur in P. aeruginosa PA14, where pyomelanin was shown to contribute to the persistence of an overproducing strain in a chronic CF infection model in mice (49). However, given that melanogenic isolates have been recovered from laboratory-grown communities of P. aeruginosa PAO1 (5, 56), it is probable that pyomelanin plays other roles in addition to protection against host defense mechanisms.As a first step toward gaining a better understanding of pyomelanin function, we sought to identify the molecular determinants of its production in P. aeruginosa. By screening a library of pTnTet/minimariner transposon mutants of a pyomelanin-overproducing clinical isolate for alterations in pigmentation, we identified several genes whose products are involved in tyrosine catabolism, central metabolic pathways, nucleotide biosynthesis, regulation, and membrane transport, in addition to hypothetical proteins of unknown function. We chose to further characterize the gene identified most frequently in our screen, one annotated as encoding a putative ATP-binding cassette of an ABC-type transporter. Here, we demonstrate that this transporter is involved in HGA transport and the subsequent extracellular formation of pyomelanin.     

Organoselenium Coating on Cellulose Inhibits the Formation of Biofilms by Pseudomonas aeruginosa and Staphylococcus aureus

Among the most difficult bacterial infections encountered in treating patients are wound infections, which may occur in burn victims, patients with traumatic wounds, necrotic lesions in people with diabetes, and patients with surgical wounds. Within a wound, infecting bacteria frequently develop biofilms. Many current wound dressings are impregnated with antimicrobial agents, such as silver or antibiotics. Diffusion of the agent(s) from the dressing may damage or destroy nearby healthy tissue as well as compromise the effectiveness of the dressing. In contrast, the antimicrobial agent selenium can be covalently attached to the surfaces of a dressing, prolonging its effectiveness. We examined the effectiveness of an organoselenium coating on cellulose discs in inhibiting Pseudomonas aeruginosa and Staphylococcus aureus biofilm formation. Colony biofilm assays revealed that cellulose discs coated with organoselenium completely inhibited P. aeruginosa and S. aureus biofilm formation. Scanning electron microscopy of the cellulose discs confirmed these results. Additionally, the coating on the cellulose discs was stable and effective after a week of incubation in phosphate-buffered saline. These results demonstrate that 0.2% selenium in a coating on cellulose discs effectively inhibits bacterial attachment and biofilm formation and that, unlike other antimicrobial agents, longer periods of exposure to an aqueous environment do not compromise the effectiveness of the coating.Among the most difficult bacterial infections encountered in treating patients are wound infections, which may occur in burn victims (10), patients with traumatic wounds (33), people with diabetes (27), and patients with surgical wounds (29, 31). Two of the more common causative agents of wound infections are Staphylococcus aureus and Pseudomonas aeruginosa (10, 27, 29, 31, 33). Such infections often lead to fatality; the mortality rate among patients infected with P. aeruginosa ranges from 26% to 55% (9, 49), while mortality from S. aureus infection ranges from 19% to 38% (28, 46, 50). As opportunistic pathogens, S. aureus and P. aeruginosa cause few infections in healthy individuals but readily cause infection once host defenses are compromised, such as with the removal of skin from burns (10). S. aureus infection originates from the normal flora of either the patient or health care workers (48), while P. aeruginosa is acquired from the environment surrounding the patient (41). Once established on the skin, S. aureus and P. aeruginosa are then able to colonize the wound. Infection results if the organisms proliferate in the wound environment.Both P. aeruginosa and S. aureus often exist within burn wounds as biofilms (43, 47). A biofilm is presently defined as a sessile microbial community characterized by cells that are irreversibly attached either to a substratum or to each other (16). Biofilms, which can attain over 100 μm in thickness, are made up of multiple layers of bacteria in an exopolysaccharide matrix (12, 16, 42). Sauer et al. showed that P. aeruginosa biofilms form in distinct developmental stages: reversible attachment, irreversible attachment, two stages of maturation, and a dispersion phase (42). Clinically, biofilms present serious medical management problems through their association with different chronic infections (37). During vascular catheter-related infections and sepsis, biofilms serve as a reservoir of bacteria from which planktonic cells detach and spread throughout the tissue and/or enter the circulatory system, resulting in bacteremia or septicemia (15). Factors specific to the bacterium may influence the formation of bacterial biofilms at different infection sites or surfaces. For example, during the initial attachment stage the flagellum, lipopolysaccharide, and possibly outer membrane proteins play a major role in bringing P. aeruginosa into proximity with the surface as well as mediating the interaction with the substratum (12). Using the murine model of thermal injury, we recently showed that P. aeruginosa forms a biofilm within the thermally injured tissues (43). Clinically, the surgeons debride the infected or dead tissues; however, a few microorganisms may remain on the tissue surface and reinitiate biofilm formation.Antibiotics, silver, or chitosan, attached to or embedded in gauze, have been shown to be efficacious in preventing wound infection (21, 24, 26, 36). However, the resistance of P. aeruginosa and S. aureus to available antibiotics severely limits the choices for antibiotic treatment (13, 40). Additionally, silver compounds, such as silver nitrate and silver sulfadiazine, leaching from dressings are toxic to human fibroblasts even at low concentrations (20, 25). Thus, effective alternative antimicrobial agents that contact the thermally injured/infected tissues and prevent the development of bacterial biofilms are required. Previous studies have shown that selenium (Se) can be covalently bound to a solid matrix and retain its ability to catalyze the formation of superoxide radicals (O₂^·−) (30). These superoxide radicals inhibit bacterial attachment to the solid surface (30). In this study, we examined the ability of a newly synthesized organoselenium-methacrylate polymer (Se-MAP) to block biofilm formation by both S. aureus and P. aeruginosa. These bacteria were chosen since they cause a major share of wound infections and because drug-resistant forms of these bacteria have become a serious problem in the treatment and management of these wound infections (6, 13, 17, 18, 38). Results of the study show that 0.2% (wt/wt) Se in Se-MAP covalently attached to cellulose discs inhibited P. aeruginosa and S. aureus biofilm formation. This could lead to the development of a selenium-based antimicrobial coating for cotton materials that will prevent the bacterial attachment and colonization that can ultimately lead to bacterial biofilm formation during chronic infections.     

Activation of ExoU Phospholipase Activity Requires Specific C-Terminal Regions

Pseudomonas aeruginosa is a Gram-negative opportunistic pathogen that utilizes a type III secretion system to subvert host innate immunity. Of the 4 known effector proteins injected into eukaryotic cells, ExoS and ExoU are cytotoxic. The cytotoxic phenotype of ExoU depends on the enzymatic activity of the patatin-like phospholipase A₂ domain localized to the N-terminal half of the protein. Amino acid residues located within the C-terminal region of ExoU are postulated to be required for trafficking or localization to the plasma membrane of eukaryotic cells. This report describes the characterization of a transposon-based linker insertion library in ExoU. Utilizing an unbiased screening approach and sensitive methods for measuring enzymatic activity, we identified regions of ExoU that are critical for activation of the phospholipase activity by the only known cofactor, SOD1. Insertions at D572 and L618 reduced the rate of substrate cleavage. Enzymatic activity could be restored to almost parental levels when SOD1 concentrations were increased, suggesting that the linker insertion disrupted the interaction between ExoU and SOD1. An enzyme-linked immunosorbent assay (ELISA)-based binding test was developed to measure ExoU-SOD1 binding. These experiments suggest that ExoU activation by SOD1 is hampered by linker insertion. ExoU derivatives harboring minimal phospholipase activity retained biological activity in tissue culture assays. These proteins affected primarily cellular architecture in a manner similar to that of ExoT. Our studies suggest that conformational changes in ExoU are facilitated by SOD1. Importantly, the level of phospholipase activity influences the biological outcome of ExoU intoxication.Pseudomonas aeruginosa is a Gram-negative bacterium responsible for severe and potentially fatal opportunistic infections. As a contributor to nosocomial infections, P. aeruginosa is a leading cause of hospital-acquired and ventilator-associated pneumonias (40). Furthermore, P. aeruginosa is responsible for ulcerative keratitis and ocular disease found in conjunction with the use of soft contact lenses (2, 10, 54). Infections with this pathogen are of critical concern for individuals admitted with severe burns, due to the bacterium''s ability to colonize and persist in damaged tissues (35). Patients suffering from cystic fibrosis often succumb to severe lung infections and inflammation due to colonization with antibi otic-resistant, mucoid strains of P. aeruginosa (3). The expression of multiple efflux pumps and the ability to inactivate and modify antibiotics make P. aeruginosa dangerous and difficult to treat (27). Several investigators are exploring ways, as adjuncts or alternatives to antibiotic treatment, to neutralize virulence factors that contribute to the ability of P. aeruginosa to suppress host innate and adaptive immune responses (17, 21, 22, 52).Many Gram-negative bacteria, including P. aeruginosa, encode one or more type III secretion systems (T3SS), which are thought to aid in pathogenesis and increase disease severity (19, 32, 39). Four effectors are translocated by the T3SS of P. aeruginosa and include ExoS, ExoT, ExoU, and ExoY (8, 23, 56, 57). The activity of each effector is dependent upon interaction with a cofactor present in eukaryotic but not prokaryotic cells. ExoS and ExoT are bifunctional enzymes that possess both Rho GTPase-activating protein and ADP-ribosyltransferase activities (23, 25, 51). The ADP ribosylation of eukaryotic proteins by ExoS and ExoT requires activation by members of the 14-3-3 family of scaffolding proteins (13). ExoY is an adenylyl cyclase that causes the accumulation of cyclic AMP (cAMP) in intoxicated cells. The eukaryotic cofactor required for ExoY activity has not been identified (57). ExoU, a potent A₂ phospholipase responsible for membrane disruption and cellular lysis, requires superoxide dismutase 1 (SOD1) for the detection of enzymatic activity (43, 46).ExoU is an important virulence factor of P. aeruginosa, as it causes rapid cell death during in vitro infections and is associated with poor clinical outcomes (19, 39, 44). Several studies have used truncation analyses, linker mutagenesis, and site-specific amino acid substitutions to define regions of ExoU important for various functions (7, 36). ExoU is a 74-kDa, hydrophilic, and slightly acidic protein with a pI of 5.9 (8). The first 52 amino acids are required for interaction with the chaperone SpcU and may be important for translocation through the T3SS (7, 9). Enzymatic activity is attributed to the patatin-like phospholipase domain located between residues 107 and 357 (34, 46). Two catalytic residues, S142 and D344, and a sequence encoding an oxyanion hole (₁₁₂GGAK₁₁₅) are located within this domain (34, 46). The oxyanion hole is thought to stabilize the negative charge of the intermediate structure during substrate cleavage (5). C-terminal residues of ExoU, specifically the last 137 amino acids, have been implicated in membrane localization after translocation into mammalian cells (37). The domain or region(s) required for the activation of ExoU by SOD1 have not been identified.In this study, linker-scanning mutagenesis (the insertion of 15 nucleotides randomly throughout the coding sequence) was used to identify regions of exoU that impair activation of phospholipase activity by SOD1. Our data support the model that SOD1 may be facilitating the activation of ExoU by altering the conformational properties of the enzyme. Understanding the molecular mechanisms mediating SOD1 and ExoU interaction may contribute to the design of therapeutics for the treatment of acute P. aeruginosa infections.     

Abundance of Novel and Diverse tfdA-Like Genes,Encoding Putative Phenoxyalkanoic Acid Herbicide-Degrading Dioxygenases,in Soil

Phenoxyalkanoic acid (PAA) herbicides are widely used in agriculture. Biotic degradation of such herbicides occurs in soils and is initiated by α-ketoglutarate- and Fe²⁺-dependent dioxygenases encoded by tfdA-like genes (i.e., tfdA and tfdAα). Novel primers and quantitative kinetic PCR (qPCR) assays were developed to analyze the diversity and abundance of tfdA-like genes in soil. Five primer sets targeting tfdA-like genes were designed and evaluated. Primer sets 3 to 5 specifically amplified tfdA-like genes from soil, and a total of 437 sequences were retrieved. Coverages of gene libraries were 62 to 100%, up to 122 genotypes were detected, and up to 389 genotypes were predicted to occur in the gene libraries as indicated by the richness estimator Chao1. Phylogenetic analysis of in silico-translated tfdA-like genes indicated that soil tfdA-like genes were related to those of group 2 and 3 Bradyrhizobium spp., Sphingomonas spp., and uncultured soil bacteria. Soil-derived tfdA-like genes were assigned to 11 clusters, 4 of which were composed of novel sequences from this study, indicating that soil harbors novel and diverse tfdA-like genes. Correlation analysis of 16S rRNA and tfdA-like gene similarity indicated that any two bacteria with D > 20% of group 2 tfdA-like gene-derived protein sequences belong to different species. Thus, data indicate that the soil analyzed harbors at least 48 novel bacterial species containing group 2 tfdA-like genes. Novel qPCR assays were established to quantify such new tfdA-like genes. Copy numbers of tfdA-like genes were 1.0 × 10⁶ to 65 × 10⁶ per gram (dry weight) soil in four different soils, indicating that hitherto-unknown, diverse tfdA-like genes are abundant in soils.Phenoxyalkanoic acid (PAA) herbicides such as MCPA (4-chloro-2-methyl-phenoxyacetic acid) and 2,4-D (2,4-dichlorophenoxyacetic acid) are widely used to control broad-leaf weeds in agricultural as well as nonagricultural areas (19, 77). Degradation occurs primarily under oxic conditions in soil, and microorganisms play a key role in the degradation of such herbicides in soil (62, 64). Although relatively rapidly degraded in soil (32, 45), both MCPA and 2,4-D are potential groundwater contaminants (10, 56, 70), accentuating the importance of bacterial PAA herbicide-degrading bacteria in soils (e.g., references 3, 5, 6, 20, 41, 59, and 78).Degradation can occur cometabolically or be associated with energy conservation (15, 54). The first step in the degradation of 2,4-D and MCPA is initiated by the product of cadAB or tfdA-like genes (29, 30, 35, 67), which constitutes an α-ketoglutarate (α-KG)- and Fe²⁺-dependent dioxygenase. TfdA removes the acetate side chain of 2,4-D and MCPA to produce 2,4-dichlorophenol and 4-chloro-2-methylphenol, respectively, and glyoxylate while oxidizing α-ketoglutarate to CO₂ and succinate (16, 17).Organisms capable of PAA herbicide degradation are phylogenetically diverse and belong to the Alpha-, Beta-, and Gammproteobacteria and the Bacteroidetes/Chlorobi group (e.g., references 2, 14, 29-34, 39, 60, 68, and 71). These bacteria harbor tfdA-like genes (i.e., tfdA or tfdAα) and are categorized into three groups on an evolutionary and physiological basis (34). The first group consists of beta- and gammaproteobacteria and can be further divided into three distinct classes based on their tfdA genes (30, 46). Class I tfdA genes are closely related to those of Cupriavidus necator JMP134 (formerly Ralstonia eutropha). Class II tfdA genes consist of those of Burkholderia sp. strain RASC and a few strains that are 76% identical to class I tfdA genes. Class III tfdA genes are 77% identical to class I and 80% identical to class II tfdA genes and linked to MCPA degradation in soil (3). The second group consists of alphaproteobacteria, which are closely related to Bradyrhizobium spp. with tfdAα genes having 60% identity to tfdA of group 1 (18, 29, 34). The third group also harbors the tfdAα genes and consists of Sphingomonas spp. within the alphaproteobacteria (30).Diverse PAA herbicide degraders of all three groups were identified in soil by cultivation-dependent studies (32, 34, 41, 78). Besides CadAB, TfdA and certain TfdAα proteins catalyze the conversion of PAA herbicides (29, 30, 35). All groups of tfdA-like genes are potentially linked to the degradation of PAA herbicides, although alternative primary functions of group 2 and 3 TfdAs have been proposed (30, 35). However, recent cultivation-independent studies focused on 16S rRNA genes or solely on group 1 tfdA sequences in soil (e.g., references 3-5, 13, and 41). Whether group 2 and 3 tfdA-like genes are also quantitatively linked to the degradation of PAA herbicides in soils is unknown. Thus, tools to target a broad range of tfdA-like genes are needed to resolve such an issue. Primers used to assess the diversity of tfdA-like sequences used in previous studies were based on the alignment of approximately 50% or less of available sequences to date (3, 20, 29, 32, 39, 47, 58, 73). Primers specifically targeting all major groups of tfdA-like genes to assess and quantify a broad diversity of potential PAA degraders in soil are unavailable. Thus, the objectives of this study were (i) to develop primers specific for all three groups of tfdA-like genes, (ii) to establish quantitative kinetic PCR (qPCR) assays based on such primers for different soil samples, and (iii) to assess the diversity and abundance of tfdA-like genes in soil.     

A New Borrelia Species Defined by Multilocus Sequence Analysis of Housekeeping Genes

Analysis of Lyme borreliosis (LB) spirochetes, using a novel multilocus sequence analysis scheme, revealed that OspA serotype 4 strains (a rodent-associated ecotype) of Borrelia garinii were sufficiently genetically distinct from bird-associated B. garinii strains to deserve species status. We suggest that OspA serotype 4 strains be raised to species status and named Borrelia bavariensis sp. nov. The rooted phylogenetic trees provide novel insights into the evolutionary history of LB spirochetes.Multilocus sequence typing (MLST) and multilocus sequence analysis (MLSA) have been shown to be powerful and pragmatic molecular methods for typing large numbers of microbial strains for population genetics studies, delineation of species, and assignment of strains to defined bacterial species (4, 13, 27, 40, 44). To date, MLST/MLSA schemes have been applied only to a few vector-borne microbial populations (1, 6, 30, 37, 40, 41, 47).Lyme borreliosis (LB) spirochetes comprise a diverse group of zoonotic bacteria which are transmitted among vertebrate hosts by ixodid (hard) ticks. The most common agents of human LB are Borrelia burgdorferi (sensu stricto), Borrelia afzelii, Borrelia garinii, Borrelia lusitaniae, and Borrelia spielmanii (7, 8, 12, 35). To date, 15 species have been named within the group of LB spirochetes (6, 31, 32, 37, 38, 41). While several of these LB species have been delineated using whole DNA-DNA hybridization (3, 20, 33), most ecological or epidemiological studies have been using single loci (5, 9-11, 29, 34, 36, 38, 42, 51, 53). Although some of these loci have been convenient for species assignment of strains or to address particular epidemiological questions, they may be unsuitable to resolve evolutionary relationships among LB species, because it is not possible to define any outgroup. For example, both the 5S-23S intergenic spacer (5S-23S IGS) and the gene encoding the outer surface protein A (ospA) are present only in LB spirochete genomes (36, 43). The advantage of using appropriate housekeeping genes of LB group spirochetes is that phylogenetic trees can be rooted with sequences of relapsing fever spirochetes. This renders the data amenable to detailed evolutionary studies of LB spirochetes.LB group spirochetes differ remarkably in their patterns and levels of host association, which are likely to affect their population structures (22, 24, 46, 48). Of the three main Eurasian Borrelia species, B. afzelii is adapted to rodents, whereas B. valaisiana and most strains of B. garinii are maintained by birds (12, 15, 16, 23, 26, 45). However, B. garinii OspA serotype 4 strains in Europe have been shown to be transmitted by rodents (17, 18) and, therefore, constitute a distinct ecotype within B. garinii. These strains have also been associated with high pathogenicity in humans, and their finer-scale geographical distribution seems highly focal (10, 34, 52, 53).In this study, we analyzed the intra- and interspecific phylogenetic relationships of B. burgdorferi, B. afzelii, B. garinii, B. valaisiana, B. lusitaniae, B. bissettii, and B. spielmanii by means of a novel MLSA scheme based on chromosomal housekeeping genes (30, 48).     

Metabolic Network Analysis of Pseudomonas aeruginosa during Chronic Cystic Fibrosis Lung Infection

System-level modeling is beginning to be used to decipher high throughput data in the context of disease. In this study, we present an integration of expression microarray data with a genome-scale metabolic reconstruction of Pseudomonas aeruginosa in the context of a chronic cystic fibrosis (CF) lung infection. A genome-scale reconstruction of P. aeruginosa metabolism was tailored to represent the metabolic states of two clonally related lineages of P. aeruginosa isolated from the lungs of a CF patient at different points over a 44-month time course, giving a mechanistic glimpse into how the bacterial metabolism adapts over time in the CF lung. Metabolic capacities were analyzed to determine how tradeoffs between growth and other important cellular processes shift during disease progression. Genes whose knockouts were either significantly growth reducing or lethal in silico were also identified for each time point and serve as hypotheses for future drug targeting efforts specific to the stages of disease progression.The last decade has witnessed an explosion in both the quantity and the pace of biological discovery. High throughput methods have been developed and leveraged at an expanding rate, with the accumulation of high throughput data outstripping the capacity for analysis using conventional methods (16, 21). To face these new challenges, systems-focused methods have come to the forefront of biological discovery, enabling a synergistic merging of network analysis with the existing reductionist paradigms that have fueled biology for the past half-century (25, 40).One of the most pressing applications of systems analysis is unraveling the myriad factors that combine to form human disease. This ambitious goal has motivated a surge of interest in the collection and analysis of microarray data, which has emerged as a dominant technology for gathering genome-scale data due to its relatively low cost, ubiquity, ease, and increasingly high resolution and reproducibility (42). In particular, microarrays for gene expression profiling have been used in longitudinal studies of disease, as it enables a glimpse at the internal changes cells undergo as a disease progresses. While many such studies have been published, very little model-driven analysis has been leveraged toward interpreting these data at the network level. There is a tremendous need for this next level of analysis, as a network approach promises a deeper mechanistic understanding of whole-cell phenotypes that will be crucial for determining better therapies in the future.With the increase in life span of cystic fibrosis (CF) patients over the last several decades, bacterial infections of the thickened mucus of the lung have become the primary disease burden that must be managed in these patients today (23). The peculiarities of the CF lung mucosal environment render it a ripe environment for growth of Pseudomonas aeruginosa in particular, a notorious opportunistic pathogen that chronically infects the lungs of nearly every CF patient by an early age (32). Due to the ability of P. aeruginosa to thrive in many varied environments and its possession of a large number of regulators, it has been hypothesized that an important determinant of the virulence of this pathogen is its exceptional metabolic versatility and adaptability (37).CF lung infections involve many adaptive stages as the bacteria respond to the host lung environment and as the lungs contemporaneously remodel based on the stresses of infection (18, 20, 35). Long-term bacterial adaptations have been studied in part through gene expression profiling, and it has been noted that a significant percentage of genes differentially expressed during chronic infection encode physiological or metabolic functions (12, 36). This finding reinforces the hypothesis that the metabolic versatility of P. aeruginosa is a large factor in its pathogenicity. As a tool in studying the metabolism of this opportunistic human pathogen, we have previously published a genome-scale reconstruction of the P. aeruginosa PAO1 strain (26). This reconstruction accounts for the functions of 1,056 genes, 883 reactions, and 760 metabolites, incorporating the functions of approximately 20% of the genes in the genome into a functional computational model that is amenable to metabolic flux-level analysis (9, 17, 31). Methods for integrating high-throughput data, including gene expression array data, with genome-scale models of metabolism in order to study tissue- or condition-dependent metabolic phenotypes are developing (1, 4-6, 22, 34). By integrating gene expression data from a longitudinal study of P. aeruginosa growth (12) with our model of P. aeruginosa metabolism (26), we are providing the first network-driven analysis of metabolic changes in P. aeruginosa growing in the CF lung. By evaluating the metabolic changes that occur in this environment, we offer a deeper understanding of how the metabolism of this pathogen adapts during a chronic CF infection and present a new way to view its evolving metabolic lifestyle.