The advent of high-throughput sequencing techniques has made it possible to follow the genomic evolution of pathogenic bacteria by comparing longitudinally collected bacteria sampled from human hosts. Such studies in the context of chronic airway infections by Pseudomonas aeruginosa in cystic fibrosis (CF) patients have indicated high bacterial population diversity. Such diversity may be driven by hypermutability resulting from DNA mismatch repair system (MRS) deficiency, a common trait evolved by P. aeruginosa strains in CF infections. No studies to date have utilized whole-genome sequencing to investigate within-host population diversity or long-term evolution of mutators in CF airways. We sequenced the genomes of 13 and 14 isolates of P. aeruginosa mutator populations from an Argentinian and a Danish CF patient, respectively. Our collection of isolates spanned 6 and 20 years of patient infection history, respectively. We sequenced 11 isolates from a single sample from each patient to allow in-depth analysis of population diversity. Each patient was infected by clonal populations of bacteria that were dominated by mutators. The in vivo mutation rate of the populations was ∼100 SNPs/year–∼40-fold higher than rates in normo-mutable populations. Comparison of the genomes of 11 isolates from the same sample showed extensive within-patient genomic diversification; the populations were composed of different sub-lineages that had coexisted for many years since the initial colonization of the patient. Analysis of the mutations identified genes that underwent convergent evolution across lineages and sub-lineages, suggesting that the genes were targeted by mutation to optimize pathogenic fitness. Parallel evolution was observed in reduction of overall catabolic capacity of the populations. These findings are useful for understanding the evolution of pathogen populations and identifying new targets for control of chronic infections.

Survival of Pseudomonas aeruginosa in cystic fibrosis (CF) chronic infections is based on a genetic adaptation process consisting of mutations in specific genes, which can produce advantageous phenotypic switches and ensure its persistence in the lung. Among these, mutations inactivating the regulators MucA (alginate biosynthesis), LasR (quorum sensing) and MexZ (multidrug-efflux pump MexXY) are the most frequently observed, with those inactivating the DNA mismatch repair system (MRS) being also highly prevalent in P. aeruginosa CF isolates, leading to hypermutator phenotypes that could contribute to this adaptive mutagenesis by virtue of an increased mutation rate. Here, we characterized the mutations found in the mucA, lasR, mexZ and MRS genes in P. aeruginosa isolates obtained from Argentinean CF patients, and analyzed the potential association of mucA, lasR and mexZ mutagenesis with MRS-deficiency and antibiotic resistance. Thus, 38 isolates from 26 chronically infected CF patients were characterized for their phenotypic traits, PFGE genotypic patterns, mutations in the mucA, lasR, mexZ, mutS and mutL gene coding sequences and antibiotic resistance profiles. The most frequently mutated gene was mexZ (79%), followed by mucA (63%) and lasR (39%) as well as a high prevalence (42%) of hypermutators being observed due to loss-of-function mutations in mutL (60%) followed by mutS (40%). Interestingly, mutational spectra were particular to each gene, suggesting that several mechanisms are responsible for mutations during chronic infection. However, no link could be established between hypermutability and mutagenesis in mucA, lasR and mexZ, indicating that MRS-deficiency was not involved in the acquisition of these mutations. Finally, although inactivation of mucA, lasR and mexZ has been previously shown to confer resistance/tolerance to antibiotics, only mutations in MRS genes could be related to an antibiotic resistance increase. These results help to unravel the mutational dynamics that lead to the adaptation of P. aeruginosa to the CF lung.