|Home | About | Journals | Submit | Contact Us | Français|
In previous work (E. E. Smith, D. G. Buckley, Z. Wu, C. Saenphimmachack, L. R. Hoffman, D. A. D'Argenio, S. I. Miller, B. W. Ramsey, D. P. Speert, S. M. Moskowitz, J. L. Burns, R. Kaul, and M. V. Olson, Proc. Natl. Acad. Sci. USA 103:8487-8492, 2006) it was shown that Pseudomonas aeruginosa undergoes intense genetic adaptation during chronic respiratory infection (CRI) in cystic fibrosis (CF) patients. We used the same collection of isolates to explore the role of hypermutation in this process, since one of the hallmarks of CRI is the high prevalence of DNA mismatch repair (MMR) system-deficient mutator strains. The presence of mutations in 34 genes (many of them positively linked to adaptation in CF patients) in the study collection of 90 P. aeruginosa isolates obtained longitudinally from 29 CF patients was not homogeneous; on the contrary, mutations were significantly concentrated in the mutator lineages, which represented 17% of the isolates (87% MMR deficient). While sequential nonmutator lineages acquired a median of only 0.25 mutation per year of infection, mutator lineages accumulated more than 3 mutations per year. On the whole-genome scale, data for the first fully sequenced late CF isolate, which was also shown to be an MMR-deficient mutator, also support these findings. Moreover, for the first time the predicted amplification of mutator populations due to hitchhiking with adaptive mutations in the course of natural human infections is clearly documented. Interestingly, increased accumulation of mutations in mutator lineages was not a consequence of overrepresentation of mutations in genes involved in antimicrobial resistance, the only adaptive trait linked so far to hypermutation in CF patients, demonstrating that hypermutation also plays a major role in P. aeruginosa genome evolution and adaptation during CRI.
Chronic respiratory infection (CRI) with Pseudomonas aeruginosa is the main driver of morbidity and mortality in cystic fibrosis (CF) patients (4, 15). The establishment of P. aeruginosa CRI is mediated by a complex adaptive process that includes physiological changes produced by the activation of specific regulatory pathways, including the induction of the biofilm mode of growth or the differential expression of virulence genes (39), and genetic changes leading to selection of an important number of adaptive mutations required for long-term persistence (19, 27, 34).
Although CF patients with CRI are generally infected by a single P. aeruginosa strain that in most cases persists throughout the patient's life (32), one of the hallmarks of such infections is the emergence and fixation over time of multiple phenotypic variants of the underlying clonal populations (24), a process known as adaptive radiation (30). Many of the selected phenotypes have been clearly linked to adaptation to the lung environment that favors lifelong persistence of CRI (27). Indeed, once the adaptation stage is reached and the CRI is fully established, eradication is generally no longer possible. At this point, with resignation, the therapeutic goals change from attempting to cure the infection to slowing the decline of lung function and improving the patient's quality of life.
The intense genetic adaptation process that takes place during the establishment of CRI has recently been quantified (34). Whole-genome sequencing of early (6 months of age) and late (96 months of age) P. aeruginosa isolates from a CF patient revealed that up to 68 mutations were acquired during establishment of the CRI. A clear signal for positive selection was demonstrated, since an extraordinarily high ratio of nonsynonymous mutations per site to synonymous mutations per site was documented. Many of the mutations detected led to loss of function of the affected genes, which were frequently involved in bacterial virulence. The work was completed by analysis of a collection of sequential P. aeruginosa isolates from 29 additional CF patients; 24 of the genes that had been mutated early in patient 1 and 10 additional genes that had been shown to be candidates for mutation in many CF infections were sequenced. Again, the signal for positive selection of mutations in these genes was overwhelming, since only 5 synonymous mutations and 103 nonsynonymous mutations were found.
One of the hallmarks of CRI with P. aeruginosa is the high prevalence of DNA mismatch repair (MMR) system-deficient hypermutable (or mutator) strains, which are detected in 30 to 60% of the patients, in contrast to acute infections, for which a prevalence of less than 1% has been reported (1, 6, 16, 21, 24, 25). While the presence of mutator strains has been found to be linked to the high rates of antibiotic resistance of P. aeruginosa clinical isolates recovered from patients with CRI (1, 7, 16, 24), the role of these strains in genetic adaptation to the CRI setting remains uncertain. Furthermore, despite the fact that the effect of hypermutation on genetic adaptation and bacterial evolution has been explored for a number of pathogens using in vivo models of colonization or infection (5, 12, 20, 23), there is still little information for natural human infectious diseases. Therefore, the goal of this work was to explore the contribution of hypermutation to the genetic adaptation of P. aeruginosa to the airways of CF patients, using the collection of clinical strains described above that were analyzed previously by Smith et al. (34).
Laboratory strains and plasmids used or constructed in this work are listed in Table Table1.1. The collection of P. aeruginosa CF isolates used was the collection previously analyzed by Smith et al. (34). This collection includes 35 longitudinal isolates from one CF patient (patient 1). The complete genomes (6.5 Mb) of the earliest available isolate (isolate 1, collected at 6 months of age) and the latest available isolate (isolate 35, collected at 96 months of age) have been sequenced, and the mutations have been detected and traced by PCR-based genotyping in the remaining intermediate isolates (34). The collection also included 90 P. aeruginosa clinical isolates obtained longitudinally from respiratory samples from 29 additional CF patients (patients 2 to 30). There were at least two clonal isolates from each patient; the earliest isolate was collected from each patient before he or she was 11 years old, and the latest isolate was collected >5 years after the recovery of the earliest isolate. The isolates were designated by using a number for identification of a patient (2 to 30), followed by the age of the patient (in years) when the isolate was collected. For this set of 90 CF isolates a total of 34 genes and their regulatory regions have been fully sequenced (34). The 34 genes sequenced include 24 of the genes that were found to be mutated early in the isolates from patient 1 and 10 additional genes that have previously been shown to be candidates for mutation in many CF infections. The list of sequenced genes includes mexZ (PA2020), lasR (PA1430), PA0313, mexA (PA0425), accC (PA4848), vfr (PA0652), mexS (PA2491), exsA (PA1713), PA0506, wspF (PA3703), rpoN (PA4462), fleQ (PA1097), mexT (PA2492), nalD (PA3574), ampD (PA4522), PA0366, cyaB (PA3217), P1-001 (in a 45-kb genomic island), PA3817, pilB (PA4526), PA1333, PA3565, PA2312, anr (PA1544), rhlR (PA3476), phoP (PA1179), rhlI (PA3476), PA4796, pqsB (PA0997), toxR (PA0707), PA2435, PA4420, ORF C/SG114, and PA2121. For each isolate, we recorded the total number of mutations for the 34 genes. We also recorded the number of new mutations that accumulated per year of infection in the sequential isolates, which was determined by dividing the number of mutations detected in an isolate that were not present in the preceding isogenic isolate by the time period (in years) between the collection times for the two isolates.
The frequencies of mutation to rifampin (300 μg/ml) resistance and/or streptomycin (1,000 μg/ml) resistance were determined by using a modified protocol based on protocols described previously (24, 25). For each strain, independent aliquots containing approximately 103 cells were inoculated into five flasks containing 10 ml of Mueller-Hinton broth and incubated at 37°C and 180 rpm for 16 to 18 h. Cultures were then collected by centrifugation and resuspended in 1 ml of a saline solution, and serial 1:10 dilutions were plated on Mueller-Hinton agar plates plates and Mueller-Hinton agar plates supplemented with 300 μg/ml of rifampin and/or 1,000 μg/ml of streptomycin. Mutation frequencies were then calculated by dividing the median numbers of mutants by the median numbers of total cells. The breakpoint used to define hypermutable strains was a frequency of mutation to rifampin resistance of >2 × 10−7, as recommended previously (24). Additionally, as a control, frequencies of mutation to rifampin and streptomycin resistance were determined with reference strain PAO1, as well as mutant derivatives with mutations in mucA (PAOMA), mutS (PAOMS), and mucA mutS (PAOMSA), which were constructed in another study (submitted for publication). At least three independent mutation frequency experiments (each with five cultures) were performed for the control strains, and the results were expressed as means ± standard deviations.
To explore the genetic basis for the mutator phenotypes, complementation studies were performed with all hypermutable strains. Plasmid pUCPMS harboring PAO1 wild-type mutS, plasmid pUCPML harboring PAO1 wild-type mutL, and plasmid pUCP24, a control cloning vector, were electroporated into the hypermutable isolates as previously described (33). The MICs of gentamicin (a resistance marker for these plasmids) were determined by a conventional broth microdilution method for the hypermutable isolates before electroporation, and depending on the level of resistance, transformants were selected on Luria-Bertani agar plates containing 50 or 250 μg/ml of the antibiotic.
Complementation was demonstrated by reversion of the increased rate of mutation to rifampin or streptomycin resistance in two independent transformant colonies for each strain. Complementation was considered positive if the mutation frequency after pUCPMS or pUCPML transformation decreased ≥2 logs, negative if the mutation frequency decreased <1 log, and partial if the reduction was at least 1 log but was not as high as the wild-type level.
To construct plasmid pUCPML, the PAO1 wild-type mutL gene was obtained from plasmid pJMML by BamHI restriction and was ligated to plasmid pUCP24 opened with the same enzyme. The ligation product was then transformed into Escherichia coli XL-1 Blue made competent by CaCl2, and the resulting transformants were selected on MacConkey agar plates containing 5 μg/ml of gentamicin and confirmed by PCR and sequencing.
The mutS and mutL genes from all the mutator strains were PCR amplified using primer pairs mutSF1-mutSR1 and mutLF-mutLR, respectively (Table (Table1).1). After purification (QIAquick PCR purification kit; Qiagen, Germany), the PCR products were fully sequenced. A BigDye Terminator kit (PE-Applied Biosystems, United States) was used for the sequencing reactions, which were analyzed with an ABI Prism 3100 DNA sequencer (PE-Applied Biosystems, United States). Additionally, the mutS and/or mutL alleles of the nonmutator isolates that preceded and followed the mutator isolates were PCR amplified and sequenced for all patients. In all cases, the mutations detected were confirmed by sequencing a new PCR product. The sequences were compared with the sequences of PAO1 using the BLAST computer program (http://www.ncbi.nlm.nih.gov/BLAST).
The epidemiological relatedness of the sequential isolates obtained from the patients infected by hypermutable strains was studied by using pulsed-field gel electrophoresis (PFGE). Bacterial DNA embedded in agarose plugs was digested with SpeI. DNA was separated in a contour-clamped homogeneous electric field DRIII apparatus (Bio-Rad, La Jolla, CA) under the following conditions: 6 V/cm2 for 26 h with pulse times of 5 to 40 s. DNA macrorestriction patterns were interpreted by using the criteria established by Tenover et al. (38).
A statistical analysis was conducted using Fischer's exact test, Student's t test, or the Mann-Whitney U test, as appropriate. A P value of <0.05 was considered significant.
In this work we analyzed the collection of P. aeruginosa CF isolates previously characterized by Smith et al. (34). This collection includes 35 longitudinal isolates from one CF patient (patient 1). The complete genomes (6.5 Mb) of the earliest available isolate (isolate 1, collected at 6 months of age) and the latest available isolate (isolate 35, collected at 96 months of age) were sequenced previously, and the mutations in the remaining intermediate isolates were detected and traced by PCR-based genotyping (34). In that work it was shown that isolate 35 (96-month isolate) had a nonsynonymous mutation in mutS (R490L), but hypermutability was not assessed because this strain was rifampin resistant. In order to determine whether this fully sequenced late CF isolate was indeed a mutator, we determined the frequencies of mutation to streptomycin resistance for isolates 1 and 35 and compared the results with the results for the control strains. As shown in Table Table2,2, isolate 35 showed a 2-log-higher mutation frequency than isolate 1. Moreover, the frequencies of mutation to streptomycin resistance for isolates 1 and 35 were very similar to those for the wild-type and mutS-deficient control strains, respectively (Table (Table2).2). Furthermore, the increased mutation frequency was restored by complementation with the plasmid harboring the wild-type mutS gene, demonstrating that the first fully sequenced late CF isolate was indeed a mutator and also that the R490L substitution, affecting DNA binding domain IV of MutS (13), results in a nonfunctional MutS protein. As expected, isolates 32 to 34, which also harbored the R490L substitution, showed increased mutation frequencies, whereas the latest isolate that did not harbor the mutation (isolate 31, 60-month isolate) showed wild-type mutation frequencies (not shown).
We then continued to characterize the collection of Smith et al. (34), which also included 90 P. aeruginosa clinical isolates obtained longitudinally from sputum samples from 29 additional CF patients (patients 2 to 30). There were at least two clonal isolates from each patient; the earliest isolate from each patient was collected before he or she was 11 years old, and the latest isolate was collected >5 years after the recovery of the earliest isolate. For this set of 90 CF isolates a total of 34 genes and their regulatory regions were fully sequenced by Smith et al. (34). Because all the isolates showed wild-type rifampin susceptibility, the frequencies of mutation to rifampin resistance were assessed for all 90 isolates. As shown in Fig. Fig.1,1, the analysis of the mutation frequencies revealed that there were 15 (16.7%) hypermutable isolates, based on the defined breakpoint (>2 × 10−7). Mutator isolates were detected for 9 of the 29 patients (31.0%); for 2 of these patients (6.9%) hypermutable strains were detected among the earliest isolates at the onset of lung infection, whereas for the 7 remaining patients (24.1%) mutators emerged during the course of the CRI.
The genetic relatedness of the sequential isolates was explored by PFGE. In all seven cases in which mutators emerged during the course of the chronic infection, the mutator isolates and the corresponding preceding nonmutator isolates showed the same macrorestriction pattern. For one additional patient (patient 8) the first isolate and the last isolate had mutator phenotypes and produced the same PFGE pattern. Finally, for the single patient in whom nonmutator isolates followed a mutator, the same PFGE pattern was also found. These results suggest that the transition from the nonmutator phenotype to the mutator phenotype (and the other way around) was a consequence of mutational events and not of strain replacement.
The genetic basis for the hypermutability of all the P. aeruginosa mutator CF isolates was investigated by performing complementation studies and by sequencing the MMR system genes. The main results are shown in Table Table3,3, and the complete allelic profiles of MMR system genes are shown in Table S1 in the supplemental material. Of the 15 hypermutable strains, 13 (86.6%) were found to be defective in either mutS (7) or mutL (6). As shown in Table Table3,3, in most cases (11 of the 13 strains), the MMR system-deficient mutator isolates contained directly inactivating mutations (frameshifts or stop codons) in mutS or mutL. As expected, in all cases, sequencing of the nonmutator isolates that preceded or followed the mutator isolates revealed the presence of identical mutS or mutL alleles except for the inactivating mutation. In the two remaining MMR-deficient isolates (22-20 and 21-19.5), missense mutations were responsible for the mutator phenotypes (Table (Table3).3). The mutS-deficient isolate 22-20 contained two nonsynonymous substitutions in mutS, A187S and R306C, although the former was also detected in the preceding nonhypermutable isogenic isolates. These results show that the previously undescribed R306C substitution, involving a highly conserved residue located in DNA binding domain IV (13), must be primarily responsible for the inactivation of MutS in this strain. The A187S substitution was also found in the preceding nonhypermutable isogenic isolates of mutS-deficient (4-bp insertion) isolate 13-19 (Table (Table3).3). Nevertheless, this polymorphism is apparently frequent in mutator isolates in several other different collections (8, 11, 21). This suggests that it may have a potential effect on MutS activity that should be explored further. The results for patient 21 are particularly interesting. Early nonhypermutable isolates (21-0.5 and 21-10.7) were followed by two sequential mutS-deficient isolates (21-14.3 and 21-15.2) and a final mutL-deficient isolate (21-19.5) (Table (Table3).3). All the isolates were found to be isogenic by PFGE, suggesting that there was independent emergence of two different mutator lineages during the chronic infection. Intriguingly, isolate 21-19.5 harbored two missense mutations in mutL, neither of which was present in the preceding isolates, despite the fact that the allelic profiles were otherwise conserved (see Table S1 in the supplemental material). It seems unlikely that acquisition of both missense mutations is required for the inactivation of MutL, since there are hundreds of possible mutations leading to one-step gene inactivation. A more likely explanation is that one of the mutations leading to MutL inactivation was selected first and that the second mutation was selected afterwards in the mutator background as a compensatory mutation to restore a potential deleterious effect of the first mutation. The emergence of two independent mutator genotypes apparently also occurred for patient 8, in which an early weak mutator (not MMR deficient) was followed by a mutS-deficient isogenic lineage (Table (Table33).
In order to analyze the effect of hypermutation during the process of genetic adaptation to the airways of CF patients, we recorded the mutations that accumulated in mutator and nonmutator isolates. Initially, we analyzed the data for patient 1. A total of 68 mutations were detected in the latest isolate (isolate 35, 96 month of age) that were not present in the earliest isolate (isolate 1, 6 months of age) by whole-genome (6.5 Mb) sequencing, and only 8 of these mutations were synonymous mutations (34). Interesting, the vast majority (42 of 68) of the detected mutations were present only in the mutator isolates (isolates 32 to 35) selected during the last 36 months of the 90-month period evaluated. On the other hand, the last nonhypermutable isolate (isolate 31, 60 months of age) had accumulated only 26 of the 68 final mutations in a longer period (54 months, from 6 to 60 months of age). According to these data, the established rate of accumulation of mutations per year of infection was notably higher in the period in which mutators were selected (14.0 mutations per year) than in the preceding period in which only nonmutators were present (5.7 mutations per year).
We then analyzed the collection of 90 clinical isolates obtained longitudinally from respiratory samples from 29 additional CF patients (patients 2 to 30) in which 34 genes and their regulatory regions had been fully sequenced (34). We recorded the total number of mutations in the 34 genes, as well as the number of new mutations accumulated per year of infection in the sequential isolates, for each mutator and nonmutator strain. As shown in Fig. Fig.2,2, there were dramatic differences in the median numbers of all mutations in the 34 genes; nonmutator isolates had a median of only 1 mutation, whereas the median number of mutations for mutators was 6.5 (P = 0.0002). Similarly, the proportion of mutators was dramatically increased for CF isolates containing a growing number of total mutations (Fig. (Fig.3).3). For instance, the proportion of mutators increased from 16.7% when all isolates were considered to 87.5% when only isolates containing at least seven mutations were considered (Fig. (Fig.3).3). Obviously, the total number of mutations in a particular isolate is expected to be influenced by the time frame of the chronic infection (late isolates should have more mutations than early isolates), and therefore the number of mutations should be higher in mutator isolates because they are overrepresented in late isolates. For instance, the mean age of patients at the time of isolation of mutator isolates (15.5 ± 5.8 years) was significantly higher (P = 0.009) than the mean age of patients at the time of isolation of nonmutator isolates (11.3 ± 5.5 years). To normalize the time frame of chronic infection, we analyzed the number of new mutations that accumulated per year of infection for the sequential isolates. As shown in Fig. Fig.4,4, notable differences were observed in the median numbers of mutations in the 34 genes; nonmutator isolates had a median of only 0.26 mutation per year of infection, whereas the median value for mutators was 0.91 mutation per year (P = 0.0002). Furthermore, these data are expected to underestimate the real rate of accumulation of new mutations in the mutator isolates, since in an important proportion of the cases these isolates were preceded by a nonmutator isolate. Therefore, in the time between the times of acquisition of the isolates there had to be a period when the isolate was nonmutator and a period when the isolate was a mutator. Consistent with this observation, the number of new mutations that accumulated per year of infection by mutators whose preceding isolate was also a mutator (3.33 new mutations per year) was significantly (P = 0.04) higher than the number of new mutations that accumulated per year of infection by the mutator isolates whose preceding isolate was a nonmutator (0.81 new mutation per year). Using the same reasoning, nonmutator isolates whose preceding isolate was a mutator should also have had a higher rate of accumulation of new mutations; indeed, the median rate for these isolates was 0.90 new mutation per year, in contrast to the rate for nonmutator isolates preceded by nonmutator isolates (0.25 new mutation per year).
Overall, these results showed that mutators dramatically enhance the process of genetic adaptation to the airways of CF patients. Subsequently, we explored whether, in addition to quantitative differences in the accumulation of new mutations, there were also qualitative differences between mutator and nonmutator isolates in the types of new mutations accumulated and in the genes affected. Table Table44 shows the distribution of the types of new mutations detected in mutator and nonmutator isolates. Overall, the observed differences were not statistically significant, although nonmutator isolates proportionally tended to have more base substitutions (particularly nonsense mutations and single-nucleotide polymorphisms in intergenic regions) and mutator isolates tended to have more insertions and deletions, specifically due to the contribution of frameshift mutations produced by short (1- to 4-bp) insertions or deletions, as would be expected for MMR-system deficient strains (10).
Table Table55 shows the distribution of the new mutations among the sequenced genes in mutator and nonmutator isolates. Since a clear link between hypermutation and development of antimicrobial resistance has been well established in previous studies (1, 7, 16, 24), we examined whether the increased numbers of mutated genes in mutator isolates were a consequence of overrepresentation of mutations in genes involved in antimicrobial resistance. As shown in Table Table5,5, the proportion of mutations that affect genes involved in antimicrobial resistance for mutator isolates (35.4%) was not higher than the proportion for nonmutator isolates (38.6%). Furthermore, significant differences were not observed for any of the genes studied, demonstrating that hypermutation has a generalized effect on genetic adaptation not linked to a specific adaptive trait. For instance, of the two genes showing a higher number of mutations in CF isolates, one (mexZ) is involved in antibiotic resistance and the other (lasR) is involved in quorum sensing regulation and virulence (34). As shown in Table Table5,5, mutations in these two genes were increased equally in mutator isolates (0.60 new mutation per mutator isolate versus 0.27 mutations per nonmutator isolate [P = 0.015] for mexZ and 0.53 new mutation per mutator isolate versus 0.19 mutation per nonmutator isolate [P = 0.008] for lasR).
Over the last decade, several in vitro and in vivo experiments have shown that mutator phenotypes may confer an evolutionary advantage during bacterial exposure to new, stressful, or fluctuating environments by increasing the rate of appearance of adaptive mutations (3, 5, 35, 37). Adaptation to the host immune system (23), antibiotic treatments (17), or viral parasites (28) has been shown to be accelerated in mutator backgrounds in various experimental models. Mutator strains have also been found in diverse natural bacterial populations with variable prevalence, although they appear to be overrepresented in chronic infections, while they are rather uncommon in acute settings (16, 24, 29, 31). The case of P. aeruginosa is highly illustrative; the prevalence of mutators has been reported to be <1% in acute infections (6), 5 to 10% in environmental and early CF isolates (11), and 30 to 60% in long-term CRI (1, 16, 24). Not surprisingly, our results obtained with the collection of CF isolates of Smith et al. (34), the first collection reported for North America, agree well within these previous findings. In most cases, mutator lineages emerged during CRI from otherwise clonal strains, although in two of the patients (6.9%) mutators were detected in the earliest isolates available. Furthermore, the vast majority of the mutator strains detected, including the first fully sequenced late CF isolate, were found to be defective in the MMR system gene mutS or mutL, also confirming previous observations (21, 25). Moreover, the documented independent selection of mutS and mutL mutants in a single patient provides further clear evidence of positive selection of MMR system-inactivating mutations during CRI.
CRI with P. aeruginosa in CF patients is one of the most illustrative examples of extreme natural bacterial evolution during the course of human infections. Indeed, the genetic properties of the bacterial cells present late in P. aeruginosa infections of CF patients differ greatly from the genetic properties of the cells that initiated the infections many years before the decline in lung function became life threatening (19, 34). Long-term persistence in the airways of CF patients seems to be driven by the selection of multiple, frequently loss-of-function mutations leading to an overall pattern of increased antimicrobial resistance, reduced virulence, and specific metabolic adaptation. Mutational inactivation of the transcriptional regulator lasR can be considered one of the paradigms of such a pattern of evolution, since it is extraordinarily frequent in the CF setting and meets all three criteria mentioned above (2). Several other frequently selected mutations, such as those leading to alginate hyperproduction, the loss of flagellum expression, or a modified lipopolysaccharide, have been demonstrated experimentally to significantly contribute to bacterial adaptation during CRI (27).
Consistent with previous hypotheses and in vitro and in vivo models, we showed that hypermutation has a dramatic effect on the intense evolution of P. aeruginosa during CRI. As shown here, the presence of adaptive traits is not homogeneous among the CF isolates; on the contrary, it is significantly concentrated in MMR system-deficient lineages. Indeed, the presence of mutations in the 34 genes studied was rather infrequent for nonhypermutable isolates (median, one mutation per isolate), whereas it was up to 6.5-fold higher in mutator backgrounds. Moreover, for the first time the predicted amplification of mutator populations by hitchhiking with adaptive mutations (18) is clearly documented in the course of natural human infections, as shown in Fig. Fig.4.4. The documented differences in the number of mutations that accumulated per year of infection are certainly overwhelming; while sequential nonmutator lineages acquire a median of only 0.25 mutation per year of infection, sequential mutator lineages accumulate over 3 mutations per year, a 13-fold-higher rate. Interestingly, increased accumulation of mutations in mutator isolates was not a consequence of overrepresentation of mutations in genes involved in antimicrobial resistance, the only adaptive trait linked so far to hypermutation in CF patients (7, 24), demonstrating that hypermutation has a generalized effect on P. aeruginosa evolution and adaptation during CRI. Our results are consistent with recent findings for the CF mouse model of chronic colonization, in which P. aeruginosa mutator strains were shown to favor long-term persistence even in the absence of antimicrobial therapy (20). They also support the results of recent studies showing that inactivation of the MMR system favors the emergence of multiple P. aeruginosa phenotypic variants in vitro, including the typical markers of CF lung infection lasR and mucA (mucA22) mutants (14, 22). Nevertheless, it should be noted that while the effect of hypermutation in bacterial evolution during CRI is clear, its implications for disease progression are still uncertain. In a very recent study, the presence of mutators was indeed associated with poorer lung function, although hypermutation could not be demonstrated to be an independent risk factor in this negative outcome (40). It should also be noted that the shortcut to genetic adaptation to the airways of CF patients provided by mutators is not expected to be cost free. While the acquisition of mutations under positive selective pressure should speed up adaptation, it may also lead to accumulation of mutations deleterious for secondary environments. Indeed, adapted mutator lineages recovered from CF mice showed reduced transmissibility in a previous study (20). Moreover, mutator isolates recovered from late CRI in CF patients have recently been shown to have reduced fitness and virulence once they are retrieved from the primary lung environment (9, 21).
This work was supported by the Ministerio de Educación y Ciencia of Spain (grant SAF2003-02851), by the Ministerio de Sanidad y Consumo, Instituto de Salud Carlos III through the Spanish Network for Research in Infectious Diseases (grants REIPI C03/14 and RD06/0008), and by the Canadian Cystic Fibrosis Foundation (D.P.S.).
Published ahead of print on 10 October 2008.
†Supplemental material for this article may be found at http://jb.asm.org/.