Previous Article |
Table of Contents
| Next Article
From the Cover
BIOLOGICAL SCIENCES / MICROBIOLOGY
Genetic adaptation by Pseudomonas aeruginosa to the airways of cystic fibrosis patients








,

,||,
,

Genome Center, *Program in Molecular and Cellular Biology, and Departments of ¶Microbiology,
Pediatrics, 
Medicine, and ||Genome Sciences, University of Washington, Seattle, WA 98195; and **Division of Infectious and Immunological Diseases, Department of Pediatrics, British Columbia Research Institute for Childrens and Womens Health, 950 West 28th Avenue, Vancouver, BC, Canada V5Z 4H4
Contributed by Maynard V. Olson, March 20, 2006
| Abstract |
|---|
|
|
|---|
In many human infections, hosts and pathogens coexist for years or decades. Important examples include HIV, herpes viruses, tuberculosis, leprosy, and malaria. With the exception of intensively studied viral infections such as HIV/AIDs, little is known about the extent to which the clonal expansion that occurs during long-term infection by pathogens involves important genetic adaptations. We report here a detailed, whole-genome analysis of one such infection, that of a cystic fibrosis (CF) patient by the opportunistic bacterial pathogen Pseudomonas aeruginosa. The bacteria underwent numerous genetic adaptations during 8 years of infection, as evidenced by a positive-selection signal across the genome and an overwhelming signal in specific genes, several of which are mutated during the course of most CF infections. Of particular interest is our finding that virulence factors that are required for the initiation of acute infections are often selected against during chronic infections. It is apparent that the genotypes of the P. aeruginosa strains present in advanced CF infections differ systematically from those of "wild-type" P. aeruginosa and that these differences may offer new opportunities for treatment of this chronic disease.
chronic infection | positive selection | virulence | antibiotic resistance
The overall picture is reminiscent of typical cancers: a clone of cells, albeit in this instance one of exogenous origin, experiences selection for an accumulation of genetic variants that promote long-term survival and clonal expansion. Our data validate this model for P. aeruginosa infections in CF and provide strong evidence for the role of selection in shaping the genotypes of the bacteria that are present during the late, life-threatening phase of the infections. Our data also focus attention on particular aspects of P. aeruginosa metabolism that are premier targets of selection, both in the patient we studied in most detail and in other, independently evolving, P. aeruginosa infections in additional CF patients.
| Results and Discussion |
|---|
|
|
|---|
97% of the composite genome in high-quality sequence data derived from each isolate. Hence, we expected to detect nearly all mutational differences between the two isolates. As is summarized in Table 1, we detected 68 mutations.
|
139 genes. Many of the mutations are likely to cause partial or complete loss of function in the altered genes because
25% of the mutations cause a premature translation stop or a shift in the reading frame, or involve either a transposon insertion or outright deletion of a gene. A standard computational method for predicting the effect of each nonsynonymous mutation on protein function suggests that approximately half (33 of 68) of all mutations in the 96-month isolate are likely to affect protein function (see Tables 4 and 5, which are published as supporting information on the PNAS web site). Many Genes Functioning in Virulence and Multidrug Efflux Were Mutated During the Patient 1 Infection. The functions of mutated genes cluster in specific categories (see Tables 2, 4, and 5). Particularly striking is the observation that a whole array of "virulence" factors, and their regulators, are mutated in the 96-month isolate. Virulence is classically defined in terms of establishment-of-infection assays, in which the loss of a virulence factor results in a decrease in the ability to cause disease. The virulence factors and regulators mutated in the 96-month isolate include genes functioning in O-antigen biosynthesis, type III secretion, twitching motility, exotoxin A regulation, multidrug efflux, osmotic balance, phenazine biosynthesis, quorum sensing, and iron acquisition. Phenotypes of the 96-month isolate, when compared with those of the 6-month isolate, are consistent with the expected effects of these mutations, including loss of serotype-specific antigenicity, loss of twitching motility, loss of pyoverdine production, loss of secreted protease activities (including elastase), and reduced biofilm formation (see Fig. 2, which is published as supporting information on the PNAS web site). Each of these phenotypes is associated with decreased virulence in models of acute infection (911).
|
-lactam antibiotic ticarcillan, perhaps as a result of mutation that occurred before our collection of this isolate. The resistance in these isolates may be the result of antibiotic treatments that patient 1 received, which included both aminoglycoside and
-lactam antibiotics. The 96-month isolate has a nonsynonymous mutation in mutS, a DNA mismatch repair gene that, when disrupted, results in an increased rate of mutation. Mutators have been reported in isolates from numerous infectious diseases, including P. aeruginosa isolates from CF airways (8), and in cells isolated from many cancers. An increased mutation rate seems beneficial when there is a need for rapid accumulation of genetic adaptations. For example, when a bacterial culture is exposed to two antibiotics successively, the proportion of mutators in the surviving population is significantly higher than in the original culture (12). Analogously, the presence of a mutator mutation in the 96-month isolate makes our observation of a strong positive-selection signal particularly significant. We were unable to determine whether the 96-month isolate had an increased rate of mutation compared with the 6-month isolate by using the standard assay for a mutator phenotype. This assay measures the frequency of mutation to resistance to the antibiotic rifampicin, to which the 96-month isolate is already resistant. Because the mutS mutation occurred relatively late in the infection (see Fig. 1 and Table 6, which is published as supporting information on the PNAS web site), it is uncertain to what extent a mutator phenotype might have influenced the existing mutations in the 96-month isolate. Despite the expectation that a mutator may produce a large number of background mutations, we did not observe this effect. As discussed above, a strikingly high proportion of the mutations we did observe appear to affect gene function.
|
This retrospective tracing of mutations present in the 96-month isolate provides no information about the heterogeneity of the infection at any given time point. However, we were able to gain some insight into the degree of heterogeneity between isolates by finding mutations in intermediate isolates that were not present in the 96-month isolate. During the PCR-based genotyping of intermediate isolates, we discovered nine additional mutations that are present only in intermediate isolates. We also found three mutations specific to the 60-month isolate by comparing its partial genome sequence, published previously (13), to the composite genome. Multiple, independently acquired mutations in mucA have also been reported (13). Viewed together, the known mutations in patient 1 isolates indicate that parallel lineages existed during the early years of the infection and coexisted, in some cases, for several years (see Fig. 1 and Table 6). For example, mutations in mucA commonly occur during CF infections, and mucA mutations arose independently in three different lineages during the patient 1 infection. Isolates from a single time point are also heterogeneous, e.g., five of the six isolates collected at 36 months had detectably different genotypes.
Signals of Positive Selection and Loss of Function Are Typical in Isolates from Chronic CF Airway Infections. A major issue is whether the mutations observed in the patient 1 infection are, to some extent, representative of the mutational change associated with P. aeruginosa infections in other CF patients with independently acquired infections. To address this issue, we tested other patients isolates for similar mutations. We assembled a collection of longitudinally collected isolates from 29 different CF patients by using the following criteria: (i) at least two clonally related isolates were available from each patient, (ii) the earliest collected isolate from each patient was collected before 11 years of age, and (iii) the latest collected isolate from each patient was collected >5 years (and in some cases 20 years) after the earliest collected isolate. In the resulting collection of 91 isolates, we could sequence a gene in each isolate and, by comparing isolates that were collected from the same patient, find all mutational differences in that gene that had occurred during 29 different CF infections. In each isolate of the collection, we chose to sequence the entire gene and regulatory region of 24 genes that had been mutated early in the patient 1 infection (genes from mutations m1, m3m15, m17m22, m24m26, and s15, as shown in Fig. 1). We also sequenced the entire gene and regulatory region of 10 genes that were candidates for mutation in many CF infections, as determined from other studies (3339). Table 3 shows the results of this sequencing, reporting a comprehensive view of all mutational differences in these genes that exist between isolates collected from the same patient.
|
The Multidrug Efflux Gene mexZ and the Quorum-Sensing Regulator lasR Are Common Mutational Targets During Chronic CF Airway Infections. The most frequently mutated genes are multidrug-efflux pumps (see Tables 3, 7, and 8). Of the genes we sequenced, the most frequently mutated gene is mexZ, a negative regulator of mexX and mexY, which are components of the MexXY-OprM multidrug-efflux pump (14, 15). Isolates with resistance to aminoglycosides often have increased expression of mexX and mexY (16), and over half of the mutations we identified in mexZ are conspicuous loss-of-function mutations that lead to increased expression of mexX and mexY. CF patients are routinely treated with antibiotics such as the aminoglycoside tobramycin, and these treatments may select for mutations in multidrug-efflux pumps and their regulators.
Virulence-related genes are also mutated during most infections, most notably the primary regulator of quorum sensing, lasR (see Tables 3, 7, and 8). Bacteria use quorum sensing to communicate, by using the concentration of a secreted small molecule to measure the prevailing cell density. On reaching high cell density, the cells dramatically alter their pattern of gene expression. At high cell density, wild-type P. aeruginosa produces virulence factors, forms biofilms, and becomes more antibiotic resistant (17, 18). All of these phenotypes are disrupted in lasR mutants. P. aeruginosa is thought to exist in a biofilm state during CF airway infections (19); if so, our results are paradoxical because they indicate that there is intense selection to mutate lasR and thus impair biofilm formation during chronic airways infections. Other commonly mutated genes that contribute to virulence are vfr, a virulence regulator that also regulates lasR; exsA, the regulator of type III secretion, a system that injects virulence factors into eukaryotic cells (20); and mexA, a component of a multidrug-efflux pump (21).
Virulence Factors Are Selected Against During Chronic CF Airway Infections. Virulence factors are, by definition, required for acute infections, but our results show that they are selected against during chronic infection of CF airways. The presumed explanation is immune evasion, because the immune system recognizes numerous virulence factors and attempts to kill those cells that express them (22), thereby selecting for a population of cells with mutations in virulence factors and virulence regulators. Other bacteria appear to reduce virulence-factor expression during chronic infections by regulation rather than by mutational mechanisms. An example is the regulated virulence system of some Bordetella species that appear to shut off virulence-factor expression and reduce biofilm development after sensing appropriate environmental conditions (23).
By sequencing 10 candidate genes in our panel of CF isolates from 29 patients, we could test whether we had enriched for commonly mutated genes by focusing on 24 genes that were mutated in the patient 1 infection. Genes from both of these groups show a pattern of mutation consistent with positive selection acting at many different genes, but only in the genes that were mutated during the patient 1 infection did we find genes such as mexA and lasR that are common mutational targets (see Table 3). Overall, we conclude that a surprisingly large number of genes in the genome are targets for mutation during adaptation to CF airways, although most of these genes are mutated in only a fraction of infections. The effect of the specific pattern of mutation in particular patients on clinical outcomes requires further study.
One concern with this study is whether the infection in patient 1 is typical of other CF infections. At the outset of this study, we did not have access to clinical information about patient 1, such as the severity of the patients disease, antibiotics used during treatment, and the CF transmembrane conductance regulator (CFTR) genotype of the patient. However, late in the study we were granted Institutional Review Board approval to learn this information. Patient 1 is homozygous for the most common CF disease allele, CFTR
F508, and was hospitalized for acute respiratory illness on multiple occasions in each of the first 7 years of life. Despite these infections, patient 1 had normal lung function at age 5 (measured by forced expiratory volume for 1 sec) but had decreased lung function by age 7, suggestive of the typical decline in lung function from chronic airways infection. Patient 1 has been treated with inhaled corticosteroids and numerous different i.v. and oral antibiotics, including
-lactam and aminoglycoside antibiotics, trimethoprim/sulfamethoxazole, and azithromycin.
Concluding Remarks. Using a detailed, whole-genome analysis, we have shown that P. aeruginosa adapts genetically to CF airways. As the cost of large-scale DNA sequencing continues to drop, this analytical approach will become increasingly practical and, therefore, applicable to many types of chronic infection. Common loss-of-function mutations in genes such as lasR and mexZ may offer new therapeutic opportunities. Whereas mutations in these genes appear to favor the clonal expansion of P. aeruginosa in the airways of CF patients, the mutations may also create new vulnerabilities that can be used for treatment. This hope arises by analogy with the situation in typical cancers, in which mutational loss of DNA-repair checkpoints favors clonal expansion but also confers sensitivity to DNA-damaging agents. What is already clear is that the genetic properties of the bacterial cells present late in P. aeruginosa infections of CF airways differ greatly from the properties of the cells that initiated the infections many years before decline in lung function became life-threatening.
| Materials and Methods |
|---|
|
|
|---|
Isolate Phenotyping. To assay for mutation rate, overnight cultures were serially diluted and plated on solid agar media containing LB or LB + 300 µg/ml rifampicin; each assay was performed in triplicate. To assay for twitching motility, single colonies were stabbed into thin (3-mm) LB agar plates and incubated at room temperature for 72 hr (25). To assay for pyoverdine production, strains were streaked onto Casamino acid-containing medium (low-iron medium), incubated at room temperature for 48 hr, and visually scored under UV light for fluorescent pyoverdine production (26). Serotyping with O1-specific antibodies was performed as described in ref. 27. Antibiotic sensitivities of isolates were determined at the time of initial isolation by broth microdilution (28). The panel of antibiotics tested included amikacin, gentamicin, tobramycin, ceftazidime, ticarcillan, aztreonam, ciprofloxacin, meropenem, and piperacillin.
Secreted protease production was assayed as described in ref. 29. Briefly, single colonies of each strain were patched onto brainheart infusion/milk agar [1.5% agar, 3.7% brainheart infusion (Difco), and 1.5% nonfat powdered milk] and incubated at 37°C for 28 hr or 45 hr. Zones of clearing surrounding bacterial growth reflected degradation of casein by secreted proteases, including the lasB protease.
Biofilm formation was assayed as described in ref. 30. Briefly, overnight cultures were diluted 1:100 in MuellerHinton broth (Difco) and adjusted for equivalent cell densities. From each culture, 100 µl was inoculated into 6 different wells of a sterile polystyrene 96-well culture plate (Nunc), covered, and incubated at 37°C for 24 hr without shaking. The culture medium was then removed, the wells were washed with sterile water three times, and the wells were stained for 15 min with 125 µl of 0.1% crystal violet (Fisher). The wells were again washed with water three times and air-dried at room temperature. The crystal violet was dissolved from the wells with 125 µl of ethanol for 15 min at room temperature, and the solution was then removed to a fresh plate for spectrophotometry. Results are representative of three separate experiments with six replicates each.
Sequencing. The genomes of isolates 1 and 35 from patient 1 were sequenced by the whole-genome shotgun method, by using standard protocols. Each genome was sequenced to 5.7x coverage, and a combined assembly was built by using reads from both isolates. Mutations were found by identifying high-quality discrepancies between the isolates. Each mutation was confirmed by whole-cell PCR and sequencing of the mutated site.
We used whole-cell PCR and sequencing to sequence the entire gene and regulatory region of 24 genes mutated early in the patient 1 infection, as well as 10 candidate genes in a panel of 91 isolates from 29 CF patients. We generated high-quality sequence data for
96% of isolates in each of these regions. In a few cases, one isolate produced no sequence for any PCR products in a given gene, suggesting that the gene was deleted; this result occurred for multiple isolates in mexZ, one isolate in lasR, and several isolates in both mexT and PA2312. When a mutation was found between clonally related isolates from the same patient, we again performed PCR and sequencing to confirm the mutation. wbpA (mutation m2) was not included in this survey because the O-antigen locus is highly divergent between strains, and dnaX (mutation m23) was not included because we identified the mutation late in the study.
Computational Tools. GENEMARK was used to make gene predictions (31). P1-001 is the temporary name for the gene in a 45-kb genomic island that is integrated at tRNA-Pro (PA2736.1) in patient 1 isolates. Gene annotations were obtained from www.pseudomonas.com.
The program SIFT was used to predict the effect of nonsynonymous mutations on protein function (32). As a control for neutral nonsynonymous mutations, we compared the single-nucleotide polymorphisms in 100 genes between the PAO1 and PA14 genome sequences. Approximately 1015% of nonsynonymous mutations in these genes were predicted to lose function; this percentage does not differ from the estimated false-positive rate in the SIFT analysis. We used the default parameters of SIFT (median conservation of sequences = 3.00 and removal of sequences >90% identical to the query) and used the National Center for Biotechnology Information nonredundant database (www.ncbi.nlm.nih.gov) (October 2004) to find homologs. A mutation was scored as "not tolerated" when SIFT predicted that the probability of observing this mutation in other functional homologs was 0.06 or less.
| Acknowledgements |
|---|
|
|
|---|
| Footnotes |
|---|
Abbreviations: CF, cystic fibrosis.
To whom correspondence may be addressed at: University of Washington, 1959 Northeast Pacific Street, Box 357730, Seattle, WA 98195. E-mail: smithee{at}u.washington.edu

To whom correspondence may be addressed at: University of Washington, Fluke Hall 316, Box 352145, Seattle, WA 98195. E-mail: mvo{at}u.washington.edu
See Commentary on page 8305.
Author contributions: E.E.S. and M.V.O. designed research; E.E.S., D.G.B., Z.W., C.S., L.R.H., D.A.D., and R.K. performed research; S.I.M., B.W.R., D.P.S., S.M.M., and J.L.B. contributed new reagents/analytic tools; E.E.S. analyzed data; and E.E.S. wrote the paper.
Conflict of interest statement: No conflicts declared.
Data deposition: The sequence reported in this paper has been deposited in the GenBank database (accession no. DQ470842).
© 2006 by The National Academy of Sciences of the USA
| References |
|---|
|
|
|---|
Related articles in PNAS:
This article has been cited by other articles in HighWire Press-hosted journals:
![]() |
E. I. Lutter, M. M. P. Faria, H. R. Rabin, and D. G. Storey Pseudomonas aeruginosa Cystic Fibrosis Isolates from Individual Patients Demonstrate a Range of Levels of Lethality in Two Drosophila melanogaster Infection Models Infect. Immun., May 1, 2008; 76(5): 1877 - 1888. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. McAlester, F. O'Gara, and J. P. Morrissey Signal-mediated interactions between Pseudomonas aeruginosa and Candida albicans J. Med. Microbiol., May 1, 2008; 57(5): 563 - 569. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Hernandez, P. Francois, L. Farinelli, M. Osteras, and J. Schrenzel De novo bacterial genome sequencing: Millions of very short reads assembled on a desktop computer Genome Res., May 1, 2008; 18(5): 802 - 809. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Yang, J. A. J. Haagensen, L. Jelsbak, H. K. Johansen, C. Sternberg, N. Hoiby, and S. Molin In Situ Growth Rates and Biofilm Development of Pseudomonas aeruginosa Populations in Chronic Lung Infections J. Bacteriol., April 15, 2008; 190(8): 2767 - 2776. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. B. Goldberg, R. E. W. Hancock, R. E. Parales, J. Loper, and P. Cornelis Pseudomonas 2007 J. Bacteriol., April 15, 2008; 190(8): 2649 - 2662. [Full Text] [PDF] |
||||
![]() |
M. A. Oberhardt, J. Puchalka, K. E. Fryer, V. A. P. Martins dos Santos, and J. A. Papin Genome-Scale Metabolic Network Analysis of the Opportunistic Pathogen Pseudomonas aeruginosa PAO1 J. Bacteriol., April 15, 2008; 190(8): 2790 - 2803. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. E. A. Zlosnik, T. J. Hird, M. C. Fraenkel, L. M. Moreira, D. A. Henry, and D. P. Speert Differential Mucoid Exopolysaccharide Production by Members of the Burkholderia cepacia Complex J. Clin. Microbiol., April 1, 2008; 46(4): 1470 - 1473. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Fox, D. Haas, C. Reimmann, S. Heeb, A. Filloux, and R. Voulhoux Emergence of Secretion-Defective Sublines of Pseudomonas aeruginosa PAO1 Resulting from Spontaneous Mutations in the vfr Global Regulatory Gene Appl. Envir. Microbiol., March 15, 2008; 74(6): 1902 - 1908. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Mathee, G. Narasimhan, C. Valdes, X. Qiu, J. M. Matewish, M. Koehrsen, A. Rokas, C. N. Yandava, R. Engels, E. Zeng, et al. Dynamics of Pseudomonas aeruginosa genome evolution PNAS, February 26, 2008; 105(8): 3100 - 3105. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Rao, A. DiGiandomenico, J. Unger, Y. Bao, R. K. Polanowska-Grabowska, and J. B. Goldberg A novel oxidized low-density lipoprotein-binding protein from Pseudomonas aeruginosa Microbiology, February 1, 2008; 154(2): 654 - 665. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. P. Bernier, D. T. Nguyen, and P. A. Sokol A LysR-Type Transcriptional Regulator in Burkholderia cenocepacia Influences Colony Morphology and Virulence Infect. Immun., January 1, 2008; 76(1): 38 - 47. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Kohler, J.-L. Dumas, and C. Van Delden Ribosome Protection Prevents Azithromycin-Mediated Quorum-Sensing Modulation and Stationary-Phase Killing of Pseudomonas aeruginosa Antimicrob. Agents Chemother., December 1, 2007; 51(12): 4243 - 4248. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. D. Macia, A. Mena, N. Borrell, J. L. Perez, and A. Oliver Increased Susceptibility to Colistin in Hypermutable Pseudomonas aeruginosa Strains from Chronic Respiratory Infections Antimicrob. Agents Chemother., December 1, 2007; 51(12): 4531 - 4532. [Full Text] [PDF] |
||||
![]() |
S. Y. Kassim, S. A. Gharib, B. H. Mecham, T. P. Birkland, W. C. Parks, and J. K. McGuire Individual Matrix Metalloproteinases Control Distinct Transcriptional Responses in Airway Epithelial Cells Infected with Pseudomonas aeruginosa Infect. Immun., December 1, 2007; 75(12): 5640 - 5650. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. S. Son, W. J. Matthews Jr., Y. Kang, D. T. Nguyen, and T. T. Hoang In Vivo Evidence of Pseudomonas aeruginosa Nutrient Acquisition and Pathogenesis in the Lungs of Cystic Fibrosis Patients Infect. Immun., November 1, 2007; 75(11): 5313 - 5324. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. M. Sandoz, S. M. Mitzimberg, and M. Schuster From the Cover: Social cheating in Pseudomonas aeruginosa quorum sensing PNAS, October 2, 2007; 104(40): 15876 - 15881. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. M. Kirov, J. S. Webb, C. Y. O'May, D. W. Reid, J. K. K. Woo, S. A. Rice, and S. Kjelleberg Biofilm differentiation and dispersal in mucoid Pseudomonas aeruginosa isolates from patients with cystic fibrosis Microbiology, October 1, 2007; 153(10): 3264 - 3274. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Hoffmann, B. Lee, M. Hentzer, T. B. Rasmussen, Z. Song, H. K. Johansen, M. Givskov, and N. Hoiby Azithromycin Blocks Quorum Sensing and Alginate Polymer Formation and Increases the Sensitivity to Serum and Stationary-Growth-Phase Killing of Pseudomonas aeruginosa and Attenuates Chronic P. aeruginosa Lung Infection in Cftr / Mice Antimicrob. Agents Chemother., October 1, 2007; 51(10): 3677 - 3687. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Duan and M. G. Surette Environmental Regulation of Pseudomonas aeruginosa PAO1 Las and Rhl Quorum-Sensing Systems J. Bacteriol., July 1, 2007; 189(13): 4827 - 4836. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. T. Kenna, C. J. Doherty, J. Foweraker, L. Macaskill, V. A. Barcus, and J. R. W. Govan Hypermutability in environmental Pseudomonas aeruginosa and in populations causing pulmonary infection in individuals with cystic fibrosis Microbiology, |