Pseudomonas aeruginosa Population Structure Revisited

Abstract

At present there are strong indications that Pseudomonas aeruginosa exhibits an epidemic population structure; clinical isolates are indistinguishable from environmental isolates, and they do not exhibit a specific (disease) habitat selection. However, some important issues, such as the worldwide emergence of highly transmissible P. aeruginosa clones among cystic fibrosis (CF) patients and the spread and persistence of multidrug resistant (MDR) strains in hospital wards with high antibiotic pressure, remain contentious. To further investigate the population structure of P. aeruginosa, eight parameters were analyzed and combined for 328 unrelated isolates, collected over the last 125 years from 69 localities in 30 countries on five continents, from diverse clinical (human and animal) and environmental habitats. The analysed parameters were: i) O serotype, ii) Fluorescent Amplified-Fragment Length Polymorphism (FALFP) pattern, nucleotide sequences of outer membrane protein genes, iii) oprI, iv) oprL, v) oprD, vi) pyoverdine receptor gene profile (fpvA type and fpvB prevalence), and prevalence of vii) exoenzyme genes exoS and exoU and viii) group I pilin glycosyltransferase gene tfpO. These traits were combined and analysed using biological data analysis software and visualized in the form of a minimum spanning tree (MST). We revealed a network of relationships between all analyzed parameters and non-congruence between experiments. At the same time we observed several conserved clones, characterized by an almost identical data set. These observations confirm the nonclonal epidemic population structure of P. aeruginosa, a superficially clonal structure with frequent recombinations, in which occasionally highly successful epidemic clones arise. One of these clones is the renown and widespread MDR serotype O12 clone. On the other hand, we found no evidence for a widespread CF transmissible clone. All but one of the 43 analysed CF strains belonged to a ubiquitous P. aeruginosa “core lineage” and typically exhibited the exoS+/exoU− genotype and group B oprL and oprD alleles. This is to our knowledge the first report of an MST analysis conducted on a polyphasic data set.

About the Authors:

Jean-Paul Pirnay

* E-mail: [email protected]

Affiliation: Laboratory for Molecular and Cellular Technology, Burn Centre, Queen Astrid Military Hospital, Brussel, Belgium

Florence Bilocq

Affiliation: Laboratory for Molecular and Cellular Technology, Burn Centre, Queen Astrid Military Hospital, Brussel, Belgium

Bruno Pot

Affiliation: Applied Maths, Sint-Martens-Latem, Belgium

Pierre Cornelis

Affiliation: Laboratory of Microbial Interactions, Department of Molecular and Cellular Interactions, Flanders Interuniversity Institute of Biotechnology, Vrije Universiteit Brussel, Brussel, Belgium

Martin Zizi

Affiliation: Department of Physiology, Vrije Universiteit Brussel, Brussel, Belgium

Johan Van Eldere

Affiliation: Rega Institute and Laboratory of Microbiology, Gasthuisberg University Hospital, Catholic University of Leuven, Leuven, Belgium

Pieter Deschaght

Affiliation: Laboratory for Bacteriology Research, Department of Chemistry, Microbiology and Immunology, Ghent University, Ghent, Belgium

Mario Vaneechoutte

Affiliation: Laboratory for Bacteriology Research, Department of Chemistry, Microbiology and Immunology, Ghent University, Ghent, Belgium

Serge Jennes

Affiliation: Laboratory for Molecular and Cellular Technology, Burn Centre, Queen Astrid Military Hospital, Brussel, Belgium

Tyrone Pitt

Affiliation: Laboratory of Healthcare Associated Infection, Specialist and Reference Microbiology Division, Health Protection Agency, London, United Kingdom

Daniel De Vos

Affiliations Laboratory for Molecular and Cellular Technology, Burn Centre, Queen Astrid Military Hospital, Brussel, Belgium, Laboratory of Microbial Interactions, Department of Molecular and Cellular Interactions, Flanders Interuniversity Institute of Biotechnology, Vrije Universiteit Brussel, Brussel, Belgium

Introduction

In his 1882 paper, “Sur les colorations bleue et verte des linges à pansements”, introduced by Louis Pasteur, Carle Gessard describes the isolation of an organism causing a blue-green coloration of wound dressings [1]. He describes this ‘accidental’ organism as colourless, globular, 1 to 1.5 thousandths of a millimetre in length, aerobic and very motile. The bacterium was named Bacillus (rod) pyocyaneus. Today we refer to this organism as Pseudomonas aeruginosa. This species is ubiquitous in the biosphere, has wide metabolic versatility and high intrinsic and acquired resistance to antimicrobials. It can be found in a wide variety of ecological environments ranging from fresh and salt water to the rhizosphere in which they colonize the endemic fauna (e.g. nematodes), flora and fungi (e.g. Pythium spp.) [2]. The opportunistic bacterium P. aeruginosa occasionally migrates from its natural environment and causes disease in animals (wild, domestic and livestock) and humans. In the latter it has emerged, partly due to its intrinsic antibiotic resistance, as a major pathogen in the airways of cystic fibrosis (CF) patients, causing often-fatal chronic respiratory infections, and as one of the most clinically significant opportunist nosocomial agents. Immunosuppressed patients such as those with severe burns, cancer or AIDS are particularly at risk.

Numerous research groups have demonstrated that P. aeruginosa clinical isolates are genotypically, chemotaxonomically, and functionally indistinguishable from environmental isolates. Römling et al. observed that the most frequently identified clone in CF patients was also detected at a relatively high frequency in aquatic environments [3] and Rahme et al. demonstrated the infectivity of a P. aeruginosa strain in both plant and animal models [4]. Similarly, P. aeruginosa strains isolated from a gasoline-contaminated aquifer were indistinguishable from clinical isolates [5] and both oil-contaminated soil isolates and clinical isolates showed pathogenic and biodegradative properties [6].

Population structure

Using multilocus enzyme electrophoresis, Maynard Smith and colleagues demonstrated that bacterial population structures could range from panmictic or fully sexual, with random association between alleles, to clonal, with nonrandom association of alleles, the latter resulting in the frequent recovery of relatively few of the many possible multilocus genotypes [7]. An intermediate type of population structure that is predominantly sexual, but harbours some epidemic clones, which show significant association between loci, was called ‘epidemic’.

The population structure of P. aeruginosa has been the subject of numerous investigations, we present an overview. Both Denamur et al. in 1993, and Picard et al. in 1994, suggested a panmictic population structure for the species but highlighted the need for caution in inferring the population structure from any single class of genetic marker [8], [9]. In 2000, comparative sequencing of 19 environmental and clinical isolates revealed a net-like population with a high frequency of recombination between isolates [10]. Using randomly amplified polymorphic DNA typing, Ruimy et al. demonstrated that bacteremia and pneumonia were not caused by specific P. aeruginosa clones [11]. In 2001 Lomholt and colleagues suggested an epidemic population structure for a P. aeruginosa population isolated mainly from patients with keratitis and their environment [12]. They found evidence for an epidemic clone that is pathogenic to the eye and is characterized by a distinct combination of virulence factors. In 2002, we combined the data obtained by 4 different typing methods, performed on a batch of 73 unrelated clinical and environmental P. aeruginosa isolates, collected across the world and observed a clear mosaicism in the results and a non-congruence between experiments, features of a panmictic population structure [13]. But, in this network we also observed some clonal complexes characterized by an almost identical data set. There was no obvious correlation between these dominant clones and habitat or, with the exception of some recent clones, their geographical origin. Therefore, we suggested an epidemic population structure for P. aeruginosa. Using multi locus sequence typing (MLST), Curran et al. confirmed in 2004 that P. aeruginosa exhibits a nonclonal epidemic population structure [14]. The P. aeruginosa population in the River Woluwe in Brussels was found to be almost as diverse as the global population, harbouring members of nearly all successful clonal complexes [15].

Several groups found that P. aeruginosa possessed a highly conserved genome, which encoded genes important for survival in numerous environments including humans and evolved through the acquisition, loss, and reorganisation of genome islands and genome islets [16]–[20]. Horizontal gene transfer (HGT) might play a more important role than point mutation in the adaptation of P. aeruginosa to different habitats. Despite not believed to be naturally competent, P. aeruginosa displays a high level of interstrain genomic plasticity and contains a high number of unfixed genes. Shen et al. put forward the idea of a population-based supra-genome that is substantially larger than the genome size of any of the component strains [21]. No two strains would be identical in terms of their genetic content and HGT continuously creates new strains with unique genetic characteristics.

Environmentally endemic bacteriophages are probably responsible for a fair amount of HGT, as they were shown to be formidable transducers of naturally occurring microbial communities of P. aeruginosa [22].

In 2006 Lee and colleagues tested the pathogenicity of diverse P. aeruginosa strains in a Caenorhabditis elegans pathogenicity model and showed that genes required for pathogenicity in one strain of P. aeruginosa were neither required for, nor predictive of virulence in other strains [23]. They concluded that virulence in P. aeruginosa is both multifactorial and combinatorial, the result of a pool of pathogenicity-related genes that interact in various combinations in different backgrounds.

In 2007 Wiehlmann and colleagues analysed 240 P. aeruginosa strains with a DNA array tube assay and reported the segregation of strains from diverse habitats and geographic origin into two large nonoverlapping clusters and 45 isolated clonal complexes composed of a few or even single strains [24]. The majority of strains belonged to a few dominant clones widespread in disease and environmental habitats.

In conclusion, there appears to be a consensus that the P. aeruginosa population structure is nonclonal epidemic, that clinical isolates are indistinguishable from environmental isolates, and that there are no specific clones with a specific (disease) habitat selection. The P. aeruginosa genome consists of a highly conserved core spiked with mobile islands and elements, which are exchanged between strains through intensive and basically phage-mediated HGT, thus creating the striking diversity of this ubiquitous opportunistic pathogen.

Despite the above-mentioned studies, some important contentious issues remain. First, since the 1980s several studies have reported the emergence, spread and persistence of multidrug resistant (MDR) clones in hospitals, mainly in intensive care wards with high antibiotic pressure. Two serotypes, O11 and O12, are highly associated with these epidemic strains [25]–[47]. Typing of these strains supported a heterogeneous population in serotype O11 but those of serotype O12 often appeared to lack significant diversity.

Second, since the second half of the 1990s, an increasing number of P. aeruginosa ‘transmissible’ CF clones have been reported worldwide [48]–[60], suggestive of an emergence of specific clones that have adapted to the CF airway environment and are spreading within CF populations.

This study

To provide a reference evolutionary framework and to position these emergent P. aeruginosa clones in the global population structure, we decided to expand our earlier study [13] both in terms of number and range of isolates and of characters investigated. Our starting material consisted of a collection of 328 unrelated isolates, collected over the last 125 years from 69 localities in 30 countries on 5 continents, including isolates from diverse clinical (human and animal) and environmental habitats (Table 1).

[Figure omitted. See PDF.]

Table 1. Origin of the P. aeruginosa strains (summary).

https://doi.org/10.1371/journal.pone.0007740.t001

Since different (genetic) markers have been shown to measure different evolutionary forces, confirming the importance of a polyphasic approach to population analysis [8], [9], [19], [61], [62], we decided to analyse and combine data from eight parameters that are equally dispersed over the P. aeruginosa genome (Table 2). The parameters investigated were i) O-serotype, ii) total genome profile by fluorescent amplified-fragment length polymorphism (FAFLP) analysis, nucleotide sequence of the outer membrane protein genes iii) oprI, iv) oprL, and v) oprD, vi) pyoverdine receptor gene profile (fpvA type and fpvB prevalence), and the prevalence of vii) exoenzyme genes exoS and exoU and viii) group I pilin glycosyltransferase gene tfpO.

[Figure omitted. See PDF.]

Table 2. Genomic localisation of the parameters investigated in this study.

https://doi.org/10.1371/journal.pone.0007740.t002

Serotyping only allows for a crude discrimination between different P. aeruginosa isolates, but because it has been performed all over the world for more than 80 years [63] it forms a bridge between old and new epidemiological studies.

FAFLP is a highly discriminatory and reproducible genotyping method based on the selective amplification of a subset of DNA fragments generated by restriction enzyme digestion [64]–[66].

Although it is generally assumed that the best means of indexing natural variation in a population structure is to sequence housekeeping genes [67] we previously showed that the DNA sequence of the oprI, oprL and oprD genes generated equally discriminative data [13]. The oprI and oprL genes, which code for outer membrane lipoproteins [68]–[71], showed sequence diversity comparable to that of housekeeping genes [13] and have been included in SNP schemes [19].

The P. aeruginosa oprD gene codes for a specialized pore protein, OprD, which allows selective permeation of basic amino acids and their structural analogs like the carbapenem antibiotics imipenem and meropenem [72], [73]. It exhibits important sequence variability with multiple non-silent mutations and a microscale mosaic structure resulting from multiple recombinational events [74]. The oprD sequence data have proven to be an extremely interesting genetic marker, for the following reasons: (i) resistance to carbapenems is often achieved by defective oprD mutations (DOMs), (ii) the mosaic structure of the oprD gene exposes evidence of recombination events between P. aeruginosa strains, (iii) the virtually unlimited number of oprD alleles provides high discriminatory power, (iv) despite this extremely high sequence variability, members of narrow clonal complexes often show identical oprD sequences, thus illustrating the stability of these complexes [13].

Pyoverdines are high-affinity fluorescent peptidic siderophores secreted by P. aeruginosa in order to scavenge Fe(III) in the extracellular environment and shuttle it into the cell [75]. Uptake of the pyoverdine-Fe(III) complex is mediated by FpvA, a specific outer membrane receptor protein. Three P. aeruginosa siderovars can be distinguished, each producing a different pyoverdine (type I, II and III) and a matching cognate FpvA receptor [76]–[78]. The type II pyoverdine receptors are more diverse and it has been suggested that they are under positive selection [79]–[80]. Two distinct type II pyoverdine receptor gene clusters were observed: IIa and IIb [81]. In 2004, an additional pyoverdine receptor, FpvB, was discovered [82]. It was found to confer, in pyoverdine type II and III producing P. aeruginosa strains, the capacity to utilize type I pyoverdine as a source of iron. The majority of P. aeruginosa strains were shown to possess the fpvB gene.

ExoS and ExoU are effector molecules (exoenzymes) that can be injected directly into the host cell by the type III secretion system. There are indications that ExoS is the major cytotoxin required for colonization and dissemination during infection, while secretion of ExoU has been associated with increased virulence [83], [84].

The pilin glycosyltransferase TfpO (also called PilO) is an inner membrane protein that captures O antigen subunits and attaches them to a serine residue at the carboxy terminus of the group I pilins [85]. The group I pilin-containing strains can be divided into subgroups: TfpOa (pilin group Ia) strains and TfpOb (pilin group Ib) strains. Analysis of pilin allele distribution among isolates from various sources revealed a striking bias in the prevalence of isolates with group I pilin genes from CF compared with non-CF human sources, suggesting that this particular pilin type may confer a colonization or persistence advantage in the CF host [86].

The above-described traits were combined and analysed using biological data analysis software. The results were visualised using a minimum spanning tree (MST).

Finally, the minimum inhibitory concentrations (MIC) of 21 antimicrobials were determined for the 328 isolates.

Results and Discussion

Serotype

Only 215 (65%) out of the 328 strains could be serotyped (Table 3). This surprisingly low percentage is partially due to the nonagglutinability of 33 out of 43 CF isolates. Additionally, we suspect that the commercially available monoclonal antibody suspensions are not as potent as some of the homemade antisera that were used in past studies. Eleven strains, including 5 CF isolates, were polyagglutinable. Nonagglutinable strains have lost most or all of their lipopolysaccharide (LPS) and polyagglutinable strains have lost part or all of their O-repeating saccharide units, which determine serotype specificity, due to a defective LPS side chain synthesis [87], [88]. Cross-reactions in agglutination are due to core LPS epitopes, which are conserved in all the serotypes. Loss of O serotype reaction was described as one of the distinctive features for P. aeruginosa strains isolated from CF patients with chronic bronchopulmonary infection [89]. Already in 1975 Zierdt and Williams reported that isolates from CF patients were frequently polyagglutinable [90]. The predominant serotypes in our collection were O11 (20.1%), O6 (14.2%), O1 (11.9%) and O12 (7.9%) (Table 3). This is in congruence with the findings of Bert and Lambert-Zechowsky, who determined the O-serotypes of 2952 P. aeruginosa isolates and found serotypes O11, O6 and O1 to be predominant [91]. The incidence of O12, however, was low. The higher prevalence of serotype O12 in our collection is due to an overrepresentation of MDR strains. As could be expected, most MDR clinical isolates exhibited serotypes O11 and O12 (Figure 1, Figure 2, Figure 3, Figure 4, and Table 3). Finally, we would like to stress that the occasional clustering of isolates with different serotypes is not necessarily the result of recombinational events. It was demonstrated that anti-pseudomonal drugs [92] and bacteriophages [93] were able to induce serotype conversion in P. aeruginosa.

[Figure omitted. See PDF.]

Figure 1. Overview of the characteristics and test results of P. aeruginosa strains 1-82/328.

ABR: antibiotic resistance; ATCC: American Type Culture Collection; BWC: burn wound centre; CC: clonal complex; CF: cystic fibrosis; COPD: chronic obstructive pulmonary disease; CPHL: Central Public Health Laboratory, London; DOM: defective oprD mutation; ESP: Ecole de Santé Publique, Brussels; FAFLP: fluorescent amplified fragment length polymorphism; HPA: Health protection Agency, Colindale; HSV: high sequence variability; IDEXX: laboratory for veterinary, food and water testing; ICU: intensive care unit; IG: insufficient growth; LMG: Laboratorium voor Microbiologie Gent, public bacteria collection; NA: no amplification; NIMR: National Institute for Medical Research, London; NT: not typable; POS: positive; TRM: reaction terminated. * ABR: antibiotic resistance, expressed as the number of antibiotic classes to which resistance was observed. ** Detected using degenerate primers.

https://doi.org/10.1371/journal.pone.0007740.g001

[Figure omitted. See PDF.]

Figure 2. Overview of the characteristics and test results of P. aeruginosa strains 83-164/328.

https://doi.org/10.1371/journal.pone.0007740.g002

[Figure omitted. See PDF.]

Figure 3. Overview of the characteristics and test results of P. aeruginosa strains 165-246/328.

https://doi.org/10.1371/journal.pone.0007740.g003

[Figure omitted. See PDF.]

Figure 4. Overview of the characteristics and test results of P. aeruginosa strains 247-328/328.

https://doi.org/10.1371/journal.pone.0007740.g004

[Figure omitted. See PDF.]

Table 3. Differential prevalence (%) of strain characteristics.

https://doi.org/10.1371/journal.pone.0007740.t003

FAFLP

The FAFLP patterns of the P. aeruginosa strains were normalised and clustered using the Unweighted Pair Group Method with Arithmetic mean (UPGMA). By applying the criteria for differentiation of P. aeruginosa by FAFLP [94], which were based on the criteria for pulsed-field gel electrophoresis [95], 44 clusters of related isolates (with ≥80% homology) were identified and numbered (Figure 1, Figure 2, Figure 3, Figure 4, and Figure S1). The close genetic relationship between some isolates, illustrated by an almost identical data set (Figure 1, Figure 2, Figure 3, Figure 4), also resulted in very similar FAFLP patterns (Figure S1). This shows that FAFLP can be used, in clinical settings for example, to recognize epidemic P. aeruginosa clones during short time spans. In contrast, the relationship between the different clonal complexes, and sometimes even between distinct clones within a complex, was not always supported by FAFLP (Figure 1, Figure 2, Figure 3, Figure 4). This is illustrated by a congruence of only 54% (lineair correlation) between the similarity matrix of FAFLP and the matrix derived from a combination of all the methods (Figure 5). FAFLP is useful to discriminate between isolates, when investigating local epidemics, but on its own it is not capable to identify clonal complexes and elucidate the population structure of P. aeruginosa.

[Figure omitted. See PDF.]

Figure 5. Congruence between experiments, as calculated using the Pearson product-moment correlation coefficient.

https://doi.org/10.1371/journal.pone.0007740.g005

ExoS and exoU

Seventy-three percent of all isolates harboured the exoS gene and 23% the exoU gene (Table 3). With the exception of three strains, the carriage of exoU and exoS was mutually exclusive and in 10 isolates neither of the genes could be amplified by PCR (Table 3). Interestingly, 42 of the 43 CF isolates exhibited the exoS+/exoU- genotype (Table 3). This could mean that the presence of exoS, which is indicative of an invasive phenotype [96], and/or the absence of exoU, which has been associated with virulence [83]–[84] and severity of disease [97], is mandatory for successful colonisation of the CF lung. These results are in congruence with earlier reports. Feltman and colleagues observed that 72% of P. aeruginosa isolates contained the exoS gene and 28% the exoU gene [98]. The presence of the exoS and exoU genes appeared to be mutually exclusive and they also observed that CF isolates harboured more frequently the exoS gene and less frequently the exoU gene than did isolates from other sites of infection, including the respiratory tract of patients without CF. Wareham and Curtis also observed an association of the exoS+/exoU− genotype with chronic infection in CF patients, whilst the exoS−/exoU+ genotype was associated with strains isolated from blood [99]. The mutual exclusion of exoS or exoU indicates that selective pressures contributed to the evolution of these genomes in different environmental niches [17]. Because the type III secretion system secretes both ExoS and ExoU, the adaptation to either one of these exoenzymes almost certainly involved interaction with different target eukaryotic organisms. Accordingly, Ferguson et al. suggested that in the transition of P. aeruginosa from the soil to certain clinical settings, the loss of ExoS expression is favoured [100]. In clinical settings the inactivation of host cell function [101] and the antiphagocytic properties [102] of ExoS should aid in the infectious process, but its limited cytotoxicity, combined with its inefficient targeting of cells of lymphoid origin, may favour the production of more cytotoxic factors, such as ExoU and exotoxin A [103], at certain sites of P. aeruginosa infection.

Kulasekara et al. suggested that the evolutionary history of the exoU locus more than likely involved transposition of the ExoU determinant onto a transmissible plasmid, followed by transfer of this plasmid into different P. aeruginosa strains [104]. This is in accordance with our results and would explain the three strains that harbour both exoS and exoU. The acquisition of novel genetic material, such as the exoU genomic island, through HGT, may enhance colonisation and survival in different host environments [17].

OprI, oprL, and oprD

The oprI, oprL, and oprD sequences of the 328 studied P. aeruginosa strains were aligned and clustered using UPGMA. Allele codes were arbitrarily assigned and consisted of a capital letter for the allele group and a number, according to their position in the alignment (Figure 6). The oprI and oprL genes showed moderate sequence variability comparable to that of housekeeping genes, as could be expected since both genes code for a structural outer membrane lipoprotein (Table 4). The oprI, oprL and oprD sequences of strains LMG 5031, Br680, CPHL 11451 and PA7 diverged considerably (Figure 6). With the exception of one isolate, all mutations in oprI and oprL were silent (http://www.pseudomonas.com/related_links.jsp#alleles). All CF isolates but one possessed the group B oprL allele (Figure 1, Figure 2, Figure 3, Figure 4, and Table 3). Since non-silent mutations are extremely rare in oprL, the conservation of distinct alleles within a clonal complex or clone is likely the result of a genetic linkage or co-selection.

[Figure omitted. See PDF.]

Figure 6. UPGMA dendrograms, with assignment of allele codes, for the oprI (a), oprL (b), and oprD (c) alleles detected in 328 P. aeruginosa strains.

https://doi.org/10.1371/journal.pone.0007740.g006

[Figure omitted. See PDF.]

Table 4. Analysis of the oprI, L, and D nucleotide sequence data in 328 P. aeruginosa strains.

https://doi.org/10.1371/journal.pone.0007740.t004

The oprD gene showed the expected high sequence variability (Figure 6C and Table 4), typical for a gene that is under strong selection for diversity (http://www.pseudomonas.com/related_links.jsp#alleles). The oprD genes of strains US376 and W15 Oct 31 could not be amplified by PCR. The oprD gene of these strains is probably not present or exhibits an aberrant nucleotide sequence (at least at the primer annealing sites). With the exception of three isolates (Li004, 5BR2 and MC086), all CF isolates exhibited a group B oprD allele (Figure 1, Figure 2, Figure 3, Figure 4, and Table 3). Although genetic linkage or co-selection between group B oprD alleles and parameters that have a significant impact on selection in the CF lung are likely, the important intergroup amino acid differences, especially in the external loops of the OprD porin [74], could be indicative of a more decisive role for OprD in the selection for strains in the CF niche. The actual weight of OprD as a selection force in the CF niche cannot be determined from our data and needs further investigation. Twenty-one different defective oprD mutations (DOMs), conferring resistance to carbapenem antibiotics, were observed (Table 3 and Table 5). Seven (16.3%) of the 43 CF isolates exhibited a DOM and as a consequence are expected to express no OprD porin. This could mean that if OprD is truly a selective force in the CF niche, it is likely that it will only be of importance in the early stages of colonisation. Additionally, this relatively high percentage of DOMs suggests that carbapenems have an impact on P. aeruginosa in the CF lung. Finally, from the congruence chart (Figure 5) we learn that the similarity matrix obtained from the oprD data alone is almost identical to that obtained from oprD, oprL and oprI combined, indicating that oprL and oprI add little or nothing to the discriminatory power of oprD.

[Figure omitted. See PDF.]

Table 5. Defective oprD mutations (DOMs) in 328 P. aeruginosa isolates.

https://doi.org/10.1371/journal.pone.0007740.t005

Pyoverdine receptors

No significant correlation could be established between the fpvA pyoverdine receptor gene type and habitat (Table 3). De Vos et al. reported a prevalence of pyoverdine type II isolates in CF patients and suggested that there might be a correlation between fpvA type and the (clinical) origin of the P. aeruginosa isolates [105]. We did observe a higher prevalence of pyoverdine type II, and more specifically type IIb (46.5%), in the CF isolates as compared to the total collection (31.0%) (Table 3), but it seems unlikely that the pyoverdine receptor is in itself a selective force in the CF niche.

The relatively unordered distribution of the different pyoverdine receptor types over the different clonal complexes is suggestive for multiple recombinatorial events involving pyoverdine receptors (Figure 1, Figure 2, Figure 3, Figure 4) and a complex evolutionary history. Tümmler and Cornelis reviewed the evolution of the pyoverdine receptor in P. aeruginosa and claimed that the pyoverdine region is the most divergent locus of the core genome because it is subject to speciation and coevolution, encodes a trait of altruistic cooperation (the production of siderophores), and encodes a receptor that is both a major fitness allele and a major deleterious allele [79]. Indeed, the mosaic dispersal of fpvA types among the different clonal complexes (Figure 1, Figure 2, Figure 3, Figure 4) is possibly the result of the selection pressure caused by bacteriocins, which use the pyoverdine receptors to enter the bacteria. Pyocin S3 was shown to use the type II FpvA receptor, while pyocin S2 was found to kill strains harbouring the type I FpvA receptor [106]–[107].

As expected, the fpvB gene was present in the majority (93.4%) of P. aeruginosa isolates, including all 43 CF isolates.

TfpO

The tfpO gene, indicative for group I pilins, was detected in 48.2% of isolates (Table 3). The tfpOa allele was very rare; it was only detected in four isolates (Table 3). The tfpO gene was present in 55.8% of CF isolates, which is only slightly higher than the average (48.2%). Thus, in contrast to Kus et al. [86] who detected the tfpO gene in 69.7% of CF isolates, we did not find a strong association of tfpO with CF. The tfpO data were found to have only very limited value and discriminatory power and were therefore not included in the combined analysis.

Population structure

MSTs have long been used in the context of mathematical topology. When a set of distances is given between entries (strains in this case), a minimum spanning tree connects all entries in such a way that the summed distance of all branches of the tree is the shortest possible [108]. In a biological context, this principle adheres to the idea that evolution should be explained in as few events as possible. MST suffers from a serious degree of degeneracy as it generates a large number of solutions, many of which have no biological relevance. Hence, priority rules are applied in order to find or assign the biologically most relevant solution amongst the many solutions. MST analysis was originally developed to link MLVA-derived sequence types (STs) [109], but technically it can be used for any data type, as long as a true distance matrix can be calculated. The MST principle, however, requires that all samples are present in the data set to construct the tree. Internal branches are normally also based upon existing samples. This means that, when an MST is calculated for evolutionary studies, there are two important conditions that have to be met: (1) the study must focus on a short time-frame, assuming that all forms or states are still present, and (2) the sampled data set must be sufficiently complete to enable the method to construct a valid tree, i.e. representing the full biodiversity of forms or states as closely as possible [108]. A major advantage of the MST approach is that the algorithm may result in trees with star-like branches, which allows for a correct classification of population systems that have a strong mutational or recombinational rate, such as P. aeruginosa, and where a large number of single locus variants (SLVs) may evolve from one common type [7]. As mentioned above, MSTs can be calculated from a true distance matrix. A distance matrix based upon a data matrix (in the case of fingerprint type data, derived after a global band matching), whether derived from one or multiple data sources, can be used. In theory, every distance coefficient applied on a data matrix produces a distance matrix suitable for analysis with the MST method [108]. Recently MST was used to determine the phylogenetic framework of Listeria monocytogenes [110]. In this study MST was used, for the first time, to link the Polyphasic Profiles (PPs) of 328 unrelated P. aeruginosa strains in such a way that the sum of the distances (number of differences between two distinct PPs) is minimized.

In our previous P. aeruginosa population structure study a UPGMA dendrogram, based on the comparison of the composite data set consisting of 4 markers in 73 strains, revealed 7 distinct clonal complexes, arbitrarily labelled CC A to CC G [13]. In the present MST, based on the composite similarity matrix derived from the combination of 7 markers in 328 strains (Figure 7), we identified 4 additional clonal complexes (CC H to CC K). The former CC C was renamed ‘clone O12’ to avoid confusion with the worldwide CF and aquatic clone C [3] and the former CC B disappeared as its members no longer clustered into a distinct clonal complex. We also observed several distinct isolates with a unique PP, some of which diverged considerably from the rest of the population (e.g. reference strains PA7 and UCBPP-PA14). Strains isolated from inanimate environments, animals and humans, separated by thousands of miles, often clustered into the same clonal complex, confirming that, in general, there is no clear correlation between the clonal complexes and geographical origin or (disease) habitat. As in our previous study, there was again strong evidence that the relation among the isolates was distorted by recombination. We observed a network of relationships between all analysed parameters (Figure 1, Figure 2, Figure 3, Figure 4) and a relatively low congruence between experiments (Figure 5). Evidence of recombination is additionally supported by the mosaic structure of the oprD gene (http://www.pseudomonas.com/related_links.jsp#alleles), which is the result of a history of intra and possibly inter species recombinational exchanges of DNA blocks [74]. We also observed several conserved clones, characterized by an almost identical data set (Figure 1, Figure 2, Figure 3, Figure 4) and represented by relatively large numbers of isolates (circles) in the MST (Figure 7). The results of this polyphasic characterization confirm the nonclonal epidemic population structure of P. aeruginosa, i.e. a superficially clonal structure with frequent recombinations, in which occasionally highly successful epidemic clones arise.

[Figure omitted. See PDF.]

Figure 7. Minimum spanning tree of the similarity matrix of the composite data set consisting of the FAFLP pattern, serotype, oprI, oprL, and oprD gene sequences, pyoverdine receptor profile and prevalence of exoS/U for 328 P. aeruginosa strains.

Each circle corresponds to a polyphasic profile (PP). The circles are scaled with member count. Branch lengths are logarithmic. Coloured zones surround PPs that belong to the same clonal complex. These complexes are also indicated with a capital letter. The lines between PPs indicate inferred phylogenetic relationships and are represented as bold, plain, discontinuous and light discontinuous depending on the number of differences between profile types. Discontinuous links are only indicative. Two bold black indent lines delimit the P. aeruginosa “core lineage”; the MDR serotype O12 clone is encircled by a red dotted line.

https://doi.org/10.1371/journal.pone.0007740.g007

A conventional UPGMA dendrogram based on the composite similarity matrix is shown in Figure S2.

CF “transmissible” clones

According to this study, a typical CF strain shows the following profile: non- or polyagglutinable (76.5%), oprD goup B (93.0%), oprL group B (97.7%), exoS+ (97.7%) and fpvB+ (100%) (Table 3). Although CF isolates exhibited a genetic diversity that was comparable to that observed in other habitats, all of them, with the exception of Li004, clustered in, or were located at the border of what appears to be a large ‘core lineage’ (Figure 7). This ‘core lineage’ seems to be predominant in disease and environmental habitats across the world and is composed of CCs A, D and J (Figure 7). Li004 was isolated from a CF patient in Lisbon (Figure 3), but it remains unclear whether it is an ‘early’ sporadic strain or a ‘late’ persistent strain. All characteristics that were associated with the CF niche (exoS, group B oprL, group B oprD, fpvA type II and presence of fpvB) in this study were also prevalent in this ‘core lineage’. This supports the argument that not one parameter in itself, but rather a multitude of linked characteristics are responsible for the selection of particular strains in the CF niche.

Although CF strains isolated in different locations across the world were shown to be genotypically non-identical and thus probably not directly related (Figure 1, Figure 2, Figure 3, Figure 4, and Figure 7), they all clustered into the ‘core lineage’. It is thus quite understandable that CF strains isolated in distant places show some level of relatedness, which should however not be confused with clonality. Lanotte et al. determined some genetic features of 162 isolates from different ecological origins [111] and found that 3 major genogroups of P. aeruginosa isolates were able to colonize CF patients. Unfortunately, due to different choices of typing techniques and strains between studies, we are not able to match these genogroups to our clonal complexes.

Occasional transmission of CF strains in CF clinics and holiday or rehabilitation camps has been reported [37], [112], [113]. Our results indicate that a widespread or global transmission of successful P. aeruginosa CF strains is unlikely to have occurred. Our data suggest that strains belonging to the successful ‘core lineage’ are ubiquitous in the natural environment and are therefore more likely to infect CF patients. In 1994, Römling et al. observed that clone C, the major clone in the CF population in Germany, was also overrepresented in soil and aquatic habitats, and suggested that the isolation frequency in CF patients simply reflected the distribution of clones in the environment [3].

MDR serotype O11 and O12 strains

All confirmed MDR O12 strains, showing resistance to one or more representatives of at least 3 antibiotic classes (Figure 1, Figure 2, Figure 3, Figure 4), clustered into a very conserved clone (Figure 7). These strains typically exhibited the following profile: serotype O12 (100%), oprD goup A (100%), oprL group B (100%), oprI group B (100%), exoS+ (100%), fpvB+ (100%) and tfpOb+ (100%) (Figure 1, Figure 2, Figure 3, Figure 4, and Table 3). MDR serotype O11 strains, in contrast, clustered into CCs F, G, H and I (Figure 7) and showed an overall higher genetic divergence (Figure 1, Figure 2, Figure 3, Figure 4, and Table 3). Serotyping of 7089 P. aeruginosa strains, isolated in 16 Belgian hospitals in the period from 1977 to 1986, revealed a steady increase of P. aeruginosa O12 isolates from 2% in 1982 to 22% in 1986 [114]. The majority of these O12 isolates showed the same distinctive pyocin and phage types, suggesting a high degree of homogeneity within the O12 strains in Belgium. A multicentre European study provided evidence for a common O12 P. aeruginosa strain in Europe [41]. In the present study, all MDR O12 strains, isolated between 1985 and 2006 in 9 countries, some of them separated by thousands of miles, were shown to cluster into a very conserved clone exhibiting virtually no divergence after more than 20 years of ‘evolution’ (Figure 1, Figure 2, Figure 3, Figure 4, and Figure 7). This MDR O12 clone consists exclusively of clinical isolates; absolutely no environmental, animal or CF isolates were part of this clone (Figure 7). Furthermore, only recent strains, isolated post 1980, clustered into this clone (Figure 1, Figure 2, Figure 3, Figure 4, and Table 3). These observations could be indicative of a recent, rapid and widespread dissemination. Natural forces are likely to sustain global dispersal of organisms as small and abundant as bacteria [115], but the increased mobility of humans and the simultaneous worldwide increase of high care facilities is likely to have accelerated the dispersal of these MDR epidemic strains. One could state that the MDR O12 clone is a genuine global epidemic clone. Strains can acquire characteristics (e.g. antibiotic resistance determinants), which are advantageous in a specific niche (e.g. an intensive care unit) and this can lead to a rapid clonal expansion. The O12 clone was, to the best of our knowledge, never isolated from the natural environment and it has been suggested that colonised or infected patients might be the primary reservoirs of the prevalent O12 clone [54], [116].

We feel that the emerging MDR O12 clone is an example of a rapid and sustained adaptation of P. aeruginosa to a novel environment. Man-made changes to the (hospital) environment, like the introduction of antimicrobials, are affecting the P. aeruginosa population structure.

Conclusion

This present study is to our knowledge the first report of an MST analysis conducted on a polyphasic data set. The population structure of P. aeruginosa was determined by means of a combination of seven valuable experiments. Analysis and clustering based on a single experiment broadly conserved the clonal complexes and clones designated in the MST based on the combined experiments (Figure 1, Figure 2, Figure 3, Figure 4). The relationship between these groups of strains, however, varied according to the considered experiment, which is visualized as a mosaic pattern in Figure 1, Figure 2, Figure 3, Figure 4. Therefore, we are convinced that the ultimate or ‘true’ population structure is most faithfully approached combining as many experiments as feasible, which are then again performed on as many unrelated and diverse strains as feasible.

This polyphasic characterization of 328 diverse and unrelated P. aeruginosa strains confirmed the nonclonal epidemic population structure of P. aeruginosa. Our results also indicate that there are no widespread CF epidemic clones. CF strains are part of a successful and ubiquitous ‘core lineage’ that have infected CF patients from the natural environment and spread through short to medium range transmission between patients in CF clinics and holiday and rehabilitation camps, possibly helped by breaches in basic infection control measures. In contrast, we report the worldwide spread and persistence of MDR clone O12. The excessive use of antibiotics has caused a worldwide ‘unnatural’ selection for multiply resistant or even panresistant P. aeruginosa strains.

We hope that the evolutionary framework presented in this study will serve as the basis for more specific studies that will prove helpful in designing public health policies (e.g. segregation of CF patients or not). Additionally, the exchange of standardized data between laboratories and the creation of international reference databases of typed microorganisms should be encouraged. It will enable an efficient monitoring of changes in microbial populations and consequently allow more adequate infection control measures. Knowing a species population structure and evolutionary paths is the cornerstone of strategies aiming to control it. Specialised follow-up papers, based on the evolutionary framework presented here and dealing with some clinically relevant issues, are in preparation.

Materials and Methods

P. aeruginosa isolates

A total of 328 P. aeruginosa clinical and environmental isolates, collected worldwide (69 localities, 30 countries and 5 continents) were examined (Table 1).

Most of them were isolated in the late eighties and nineties, but 49 were isolated before 1980, including 14 P. aeruginosa strains isolated at the ‘Institut Pasteur’ in Paris by Carle Gessard [1] and his colleagues in the late 19th century. The studied collection contained 185 human clinical isolates (including 43 CF, 33 burn, 32 wound, 18 urine, 15 sputum, 6 faeces and 5 blood isolates), 63 animal clinical isolates (39 dogs, 6 turtles, 4 horses, 3 parrots, 3 dolphins, 2 cats, 2 cows, 1 kangaroo, 1 goat, 1 rabbit and 1 seal) and 55 environmental isolates (17 sea water, 16 river water, 6 lake water, 5 turtle egg, 4 plant, 3 tap water, 2 swimming pool water and 2 drink water isolates). Geographical origin, isolation site and time and other relevant characteristics of all P. aeruginosa isolates can be found in Figure 1, Figure 2, Figure 3, Figure 4.

The P. aeruginosa strains were kindly provided by: Dr. A. T. McManus, US Army Institute of Surgical Research, Texas, USA; Dr. L. Ménesi, General Hospital St. Istvan, Budapest, Hungary; Dr. A. Vanderkelen, Dr. S. Jennes, G. Verbeken and D. Schoeters, Queen Astrid Military Hospital, Neder-Over-Heembeek (Brussels), Belgium; Dr. J. A. Clark, Queen Mary's University Hospital, London, England; Dr. A. F. Vloemans, Rode Kruis Ziekenhuis, Beverwijk, The Netherlands; Dr. T. Taddonio, University of Michigan, Michigan, USA; Dr. A. Radke, Klinik für Verbrennungs- und Plastische Wiederherstellungschirurgie, Aachen, Germany; Prof. R. Konigova, Charles University Hospital, Prague, Czech Republic; Dr. R. G. Tompkins, Burns Institute, Shriners Hospital for Children, Boston, USA; Prof. B. Tümmler, Medizinische Hochschule, Hannover, Germany; Dr. M. Caneira, Hospital de Santa Maria, Lisboa, Portugal; Prof. A. Boudabous, Science Faculty, Tunis, Tunisia; Dr. M. Mergeay, Environmental Technology Expertise Centre, Mol, Belgium; Dr. A. E. Lim Jr., St. Scholastica's College of Health Sciences, Tacloban City, Philippines; Prof. O. Hadjiiski, Scientific Institute of Emergency Medicine Pirogov, Sofia, Bulgaria; Prof. K. Taviloglu, University of Istanbul, Istanbul, Turkey; Dr. W. D. H. Hendriks, Zuiderziekenhuis, Rotterdam, The Netherlands; Dr. G. Wauters, University of Louvain, Brussels, Belgium; Dr. O. Vandenberg, Universitair Ziekenhuis St.-Pierre, Brussels, Belgium; Prof. M. Vaneechoutte, University Hospital Ghent, Gent, Belgium; Prof. J. Van Eldere, Catholic University of Leuven, Leuven, Belgium; Dr. U. Römling, Karolinska Institute, Stockholm, Sweden; Dr. L. Roddam and R. Bradbury, University of Tasmania, Hobart, Australia; Dr. T. L. Pitt, Health Protection Agency, London, UK; Dr. A. Leitão, Faculty of Veterinary Medicine, Lisboa, Portugal; Dr. R. W. Brimicombe, Haga Ziekenhuis, Den Haag, The Netherlands; Prof. N. J. Legakis and Dr. P. T. Tassios, National and Kapodastrian University of Athens, Athens, Greece; Prof. J.-M. Meyer, University Louis Pasteur, Strasbourg, France and Dr. M. P. Crespo, Universidad Santiago de Cali, Cali, Colombia; Dr. M. Merabishvili and Dr. Nina Chanishvili, Eliava Institute, Tbilisi, Georgia; Dr. L. Griffiths, Dr. K. Craven and J. Awong-Taylor, Armstrong Atlantic State University, Savannah, US; Prof. M. de Chial, University of Panama, Panama City, Panama; Dr. N. H. Khan, Dr. N. Kimata and Prof. K. Kogure, University of Tokyo, Tokyo, Japan; Dr. D. Armstrong, Monash Medical Center, Melbourne, Australia; A. Catrijsse, Vlaams Instituut voor de Zee, Oostende, Belgium.

Strains LMG 2107, 5031, 10643, and 14083-5 were purchased from the BCCM™/LMG bacteria collection. The 20 ‘CPHL (Central Public health Laboratory) strains’ were purchased from the National Collection of Type Cultures in London (UK). Strain PAO-1 was kindly provided by Dr. C. K. Stover (PathoGenesis Corporation, Seattle, USA). Strain ATCC 27853 was purchased from Gibson Laboratories (USA).

All isolates were grown overnight in Luria-Bertani broth medium (Gibco-BRL-Life Technologies, Belgium) at 37°C on a rotary shaker (150 rpm). The overnight cultures were mixed with equal amounts of sterile 50% (vol/vol) glycerol (Sigma Aldrich, Belgium) in PBS buffer (Sigma Aldrich, Belgium) and stored in duplicate at −80°C.

FAFLP

FAFLP utilized an ABI 377 automated fluorescence sequencer (Applied Biosystems, Belgium), and the AFLP™ Microbial Fingerprinting Kit (Applied Biosystems) as detailed in the manufacturer's protocols. The enzymes used were T4 DNA ligase, EcoRI, and Tru9I (all purchased from Roche Diagnostics, Belgium). The primer pair used was EcoRI-0[FAM]/MseI-C. GeneScan-500[ROX] internal standard (Applied Biosystems) was co-electrophoresed with each sample in order to allow an accurate calculation of fragment lengths and correction for variation rates and gel distortions. Normalization and fragment sizing were carried out using GeneScan software (Applied Biosystems, Belgium). Band patterns were imported into the BioNumerics v5.1 software (Applied Maths, Belgium) and normalised; parameters used: background subtraction (10% disc diameter), filtering (arithmetic average), band search (minimum profiling 0.5% relative to the maximum value). Cluster analysis was performed by pairwise calculation of the Pearson correlation; the similarity matrix was clustered using the UPGMA algorithm with optimisation: 0%, position tolerance: 1%; uncertain bands were ignored.

Serotyping

Strains were grown overnight on Luria-Bertani agar medium (Gibco-BRL-Life Technologies) at 37°C. Isolates were serotyped by slide agglutination according to the International Antigenic Typing Scheme (IATS) for P. aeruginosa [117], using a panel of 16 type O monovalent antisera (Bio-Rad, Belgium). Some strains had already been serotyped by the strain providers (e.g. isolates LMG 14084 and CPHL 12447, which were shown to exhibit the provisionary O17 and O18 serotypes).

PCR and sequencing

Strains were grown overnight in Luria-Bertani broth medium (Gibco-BRL-Life Technologies) at 37°C on a rotary shaker (150 rpm). DNA was extracted from the overnight cultures using the High Pure™ PCR Template Preparation Kit (Roche Diagnostics) according to the manufacturer's guidelines. The complete oprI, oprL, and oprD genes and a fragment of the exoS, exoU, fpvA, fpvB and tfpO genes were amplified by PCR, using the primers described in Table 6. PCR was performed in 200-µl microcentrifuge tubes. The PCR mixture (50 µl final volume) contained the following: 25.5 µl sterile distilled water, 5 µl 10×PCR buffer (500 mmol/l KCl and 100 mmol/l Tris-HCl: pH 8.3), 4 µl of a deoxynucleotide mixture (dGTP, dTTP, dATP, and dCTP; 2 mmol/l each), 5 µl MgCl2 (2.5 mmol/l), 5 µl of a primer mixture (10 µmol/l each), 5 µl template DNA, and 0.5 µl AmpliTaq DNA polymerase (5 U/µl). All PCR-reagents and primers were ordered from PE-Applied Biosystems. The amplification was performed in a GeneAmp® 9700 thermocycler (Applied Biosystems). The amplification program was set at 50 cycles of denaturation at 94°C for 30 s, annealing at a temperature in accordance to the primers (Table 6), for 30 s, and elongation at 72°C for 1 min. In a few strains the amplification of the fpvA gene required degenerate primers. The reaction mixture was put on an agarose gel of 1.5 % (wt/vol) for electrophoresis and visualization of the PCR-product after staining with ethidium bromide on a transilluminator. Prior to the sequencing of the oprD, oprL and oprI genes, the respective PCR-products were purified, using centricon® 100 micro-concentrators (Millipore, Brussels, Belgium) according to the manufacturer's instructions. Five µl of the purified PCR fragment was used as a template in the sequencing reaction. PCR primers were used for sequencing. Sequencing of the coding and anti-coding strand of the oprD PCR products necessitated two additional internal primers (Table 6). DNA sequencing utilized an ABI 377 automated fluorescence sequencer (Applied Biosystems), and the ABI Prism® BigDye™ Terminator cycle sequencing kit (Applied Biosystems) as detailed in the manufacturer's protocols. The oprD gene of isolate Be128 was sequenced directly from genomic DNA. Some genes were sequenced in the VIB Genetic Service Facility (Belgium) using a capillary Applied Biosystems 3730 DNA Analyzer. PCR and sequencing were performed in duplicate in order to be able to detect eventual PCR mistakes.

[Figure omitted. See PDF.]

Table 6. Primers for PCR and sequencing.

https://doi.org/10.1371/journal.pone.0007740.t006

Using the BioNumerics v5.1 software, sequences were grouped via a pairwise clustering (pairwise alignment parameters: open gap penalty: 100%, unit gap penalty: 0%, min. match sequence: 2, max n° of gaps: 9, fast algorithm). The obtained UPGMA tree was used to seed a multiple alignment (multiple alignment parameters: open gap penalty: 100%, unit gap penalty: 0%, min. match sequence: 2, max n° of gaps: 98). Finally, multiple aligned sequences were clustered using the same parameters as used in the initial pairwise clustering, resulting in the final UPGMA tree.

Conbined data analysis

A data set consisting of the serotype, FAFLP pattern, oprI, oprL, and oprD gene sequences, pyoverdine receptor profile (fpvA and fpvB) and prevalence of the genes exoS and exoU of 328 P. aeruginosa isolates was analyzed using the biological data analysis software BioNumerics v5.1. The settings used for the comparison of the FAFLP fingerprints and the gene sequences are described in the respective paragraphs. Serotype, pyoverdine receptor profile and presence of exoS/U were compared using the Pearson correlation. These individual comparisons resulted in individual similarity matrices, which were averaged into the similarity matrix of the composite data set. No correction for internal weights was applied. Each isolate was thus assigned a ‘polyphasic profile’ (PP) contributing to the composite similarity matrix. Grouping of the averaged composite similarity matrix was achieved by MST analysis using BioNumerics v5.1 software. The MST coefficient was taken from the composite similarity matrix. The Degeneracy of the MST was reduced through the use of a priority rule by which types that had a maximum number of entries were linked first, confirming a biological meaning that these clones are most likely older. For visual purposes, isolates were further grouped into clonal complexes. For the creation of the clonal complexes, the similarity bin size (1 change) was set to 2.5%; the maximal neighbour distance between two complexes was 5 changes (12.5%) and the minimum size of a complex was 5 types. Originally a clonal complex was defined as a cluster of STs in a burstdiagram in which all STs are linked as SLVs to at least one other ST. In our case a clonal complex is a cluster of PPs, after MST analysis, in which all PPs with less than 5 changes ( =  less than 12.5% distance in the similarity matrix) are linked. Congruence between experiments was calculated using the Pearson product-moment correlation coefficient between the respective similarity matrices.

Antimicrobial susceptibility tests

Strains were grown 18–24 h at 37°C on Columbia agar containing 5% horse blood (bioMérieux). Suspensions of these cultures were made in 0.45% saline, adjusted to the turbidity of a 0.6 McFarland standard, and used to load the test cards for VITEK 2 (bioMérieux), which was used in accordance with the manufacturer's directions. The following antibiotics were tested using the AST-N020 antimicrobial susceptibility cards: AMP, ampicilin; AMC, amoxicillin + clavulanic acid; PIP, piperacillin; TZP, piperacillin + tazobactam; CEF, cephalothin; CXM, cefuroxime; CTX, cefotaxime; CAZ, ceftazidime; CPD, cefpodoxime; FOX, cefoxitin; FEP, cefepime; MEM, meropenem; GEN, gentamicin; TOB, tobramycin; AMK, amikacin; NOR, norfloxacin; OFX, ofloxacin; CIP, ciprofloxacin; NIT, nitrofurantoin; SXT, trimethoprim + sulfamethoxazole. Antibiotic resistance phenotypes, represented by the minimum inhibitory concentrations (MICs) for the above-mentioned antibiotics, were determined using VITEK 2 Advanced Expert System (AES) [118]. P. aeruginosa ATCC 27853 was used as control strain. For some isolates the MIC was determined by the broth microdilution method [119].

Nucleotide sequences

The nucleotide sequences generated in this study have been deposited in the Pseudomonas aeruginosa Genome database (http://www.pseudomonas.com/related_links.jsp#alleles).

Strain collection

All studied P. aeruginosa strains were deposited in the Belgian Coordinated Collections of Microorganisms (BCCM) of the Laboratorium voor Microbiologie (LMG) of the Ghent University. Strains were assigned a BCCM/LMG number (LMG 24881 - 25202). Strains that were obtained from a culture collection (BCCM/LMG or ATCC) maintained their original reference number.

Strains can be obtained from the LMG bacteria collection for research use only and with the consent of the strain donors.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

UPGMA dendrogram of the FAFLP patterns of the 328 studied P. aeruginosa strains.

https://doi.org/10.1371/journal.pone.0007740.s001

(0.42 MB PDF)

Figure S2.

UPGMA dendrogram of the similarity matrix of the composite data set consisting of the serotype, FAFLP pattern, oprI, L, and D gene sequences, pyoverdine receptor profile and prevalence of exoS/U genes for the 328 studied P. aeruginosa strains.

https://doi.org/10.1371/journal.pone.0007740.s002

(0.02 MB PDF)

Acknowledgments

We thank all the collectors and providers of P. aeruginosa isolates.

Author Contributions

Conceived and designed the experiments: JPP PC DDV. Performed the experiments: JPP FB. Analyzed the data: JPP FB BP PC MZ JVE PD MV SJ TP DDV. Contributed reagents/materials/analysis tools: JPP BP PC MZ JVE PD MV SJ TP DDV. Wrote the paper: JPP BP PC MZ JVE MV SJ TP DDV.

Citation: Pirnay J-P, Bilocq F, Pot B, Cornelis P, Zizi M, Van Eldere J, et al. (2009) Pseudomonas aeruginosa Population Structure Revisited. PLoS ONE4(11): e7740. https://doi.org/10.1371/journal.pone.0007740

References

1. Gessard C (1882) Sur les colorations bleue et verte des linges à pansements. Comptes-rendus hebdomadaire des séances de l'Académie des Sciences 94: 536–538.C. Gessard1882Sur les colorations bleue et verte des linges à pansements.Comptes-rendus hebdomadaire des séances de l'Académie des Sciences94536538

2. Goldberg JB (2000) Pseudomonas: global bacteria. Trends Microbiol 8: 55–57.JB Goldberg2000Pseudomonas: global bacteria.Trends Microbiol85557

3. Römling U, Wingender J, Müller H, Tümmler B (1994) A major Pseudomonas aeruginosa clone common to patients and aquatic habitats. Appl Environ Microbiol 60: 1734–1738.U. RömlingJ. WingenderH. MüllerB. Tümmler1994A major Pseudomonas aeruginosa clone common to patients and aquatic habitats.Appl Environ Microbiol6017341738

4. Rahme LG, Stevens EJ, Wolfort SF, Shao J, Tompkins RG, et al. (1995) Common virulence factors for bacterial pathogenicity in plants and animals. Science 268: 1899–1902.LG RahmeEJ StevensSF WolfortJ. ShaoRG Tompkins1995Common virulence factors for bacterial pathogenicity in plants and animals.Science26818991902

5. Foght JM, Westlake DWS, Johnson WM, Ridgway HF (1996) Environmental gasoline-utilizing isolates of Pseudomonas aeruginosa are taxonomically indistinguishable by chemotaxonomic and molecular techniques. Microbiology 142: 2333–2340.JM FoghtDWS WestlakeWM JohnsonHF Ridgway1996Environmental gasoline-utilizing isolates of Pseudomonas aeruginosa are taxonomically indistinguishable by chemotaxonomic and molecular techniques.Microbiology14223332340

6. Alonso A, Rojo F, Martinez JL (1999) Environmental and clinical isolates of Pseudomonas aeruginosa show pathogenic and biodegradative properties irrespective of their origin. Environ Microbiol 1: 421–430.A. AlonsoF. RojoJL Martinez1999Environmental and clinical isolates of Pseudomonas aeruginosa show pathogenic and biodegradative properties irrespective of their origin.Environ Microbiol1421430

7. Maynard Smith J, Smith NH, O'Rourke M, Spratt BG (1993) How clonal are bacteria? Proc Natl Acad Sci U S A 90: 4384–4388.J. Maynard SmithNH SmithM. O'RourkeBG Spratt1993How clonal are bacteria?Proc Natl Acad Sci U S A9043844388

8. Denamur E, Picard B, Decoux G, Denis JB, Elion J (1993) The absence of correlation between allozyme and rrn RFLP analysis indicates a high gene flow rate within human clinical Pseudomonas aeruginosa isolates. FEMS Microbiol Lett 110: 275–280.E. DenamurB. PicardG. DecouxJB DenisJ. Elion1993The absence of correlation between allozyme and rrn RFLP analysis indicates a high gene flow rate within human clinical Pseudomonas aeruginosa isolates.FEMS Microbiol Lett110275280

9. Picard B, Denamur E, Barakat A, Elion J, Goullet P (1994) Genetic heterogenicity of Pseudomonas aeruginosa clinical isolates revealed by esterase electrophoretic polymorphism and restriction fragment length polymorphism of the ribosomal RNA gene region. J Med Microbiol 40: 313–322.B. PicardE. DenamurA. BarakatJ. ElionP. Goullet1994Genetic heterogenicity of Pseudomonas aeruginosa clinical isolates revealed by esterase electrophoretic polymorphism and restriction fragment length polymorphism of the ribosomal RNA gene region.J Med Microbiol40313322

10. Kiewitz C, Tümmler B (2000) Sequence diversity of Pseudomonas aeruginosa: impact on population structure and genome evolution. J Bacteriol 182: 3125–3135.C. KiewitzB. Tümmler2000Sequence diversity of Pseudomonas aeruginosa: impact on population structure and genome evolution.J Bacteriol18231253135

11. Ruimy R, Genauzeau E, Barnabe C, Beaulieu A, Tibayrenc M, et al. (2001) Genetic diversity of Pseudomonas aeruginosa strains isolated from ventilated patients with nosocomial pneumonia, cancer patients with bacteremia, and environmental water. Infect Immun 69: 584–588.R. RuimyE. GenauzeauC. BarnabeA. BeaulieuM. Tibayrenc2001Genetic diversity of Pseudomonas aeruginosa strains isolated from ventilated patients with nosocomial pneumonia, cancer patients with bacteremia, and environmental water.Infect Immun69584588

12. Lomholt JA, Poulsen K, Kilian M (2001) Epidemic population structure of Pseudomonas aeruginosa: evidence for a clone that is pathogenic to the eye and that has a distinct combination of virulence factors. Infect Immun 69: 6284–6295.JA LomholtK. PoulsenM. Kilian2001Epidemic population structure of Pseudomonas aeruginosa: evidence for a clone that is pathogenic to the eye and that has a distinct combination of virulence factors.Infect Immun6962846295

13. Pirnay JP, De Vos D, Cochez C, Bilocq F, Vanderkelen A, et al. (2002) Pseudomonas aeruginosa displays an epidemic population structure. Environ Microbiol 4: 898–911.JP PirnayD. De VosC. CochezF. BilocqA. Vanderkelen2002Pseudomonas aeruginosa displays an epidemic population structure.Environ Microbiol4898911

14. Curran B, Jonas D, Grundmann H, Pitt T, Dowson CG (2004) Development of a multilocus sequence typing scheme for the opportunistic pathogen Pseudomonas aeruginosa. J Clin Microbiol 42: 5644–5649.B. CurranD. JonasH. GrundmannT. PittCG Dowson2004Development of a multilocus sequence typing scheme for the opportunistic pathogen Pseudomonas aeruginosa.J Clin Microbiol4256445649

15. Pirnay JP, Matthijs S, Colak H, Chablain P, Bilocq F, et al. (2005) Global Pseudomonas aeruginosa biodiversity as reflected in a Belgian river. Environ Microbiol 7: 969–980.JP PirnayS. MatthijsH. ColakP. ChablainF. Bilocq2005Global Pseudomonas aeruginosa biodiversity as reflected in a Belgian river.Environ Microbiol7969980

16. Spencer DH, Kas A, Smith EE, Raymond CK, Sims EH, et al. (2003) Whole-genome sequence variation among multiple isolates of Pseudomonas aeruginosa. J Bacteriol 185: 1316–1325.DH SpencerA. KasEE SmithCK RaymondEH Sims2003Whole-genome sequence variation among multiple isolates of Pseudomonas aeruginosa.J Bacteriol18513161325

17. Wolfgang MC, Kulasekara BR, Liang X, Boyd D, Wu K, et al. (2003) Conservation of genome content and virulence determinants among clinical and environmental isolates of Pseudomonas aeruginosa. Proc Natl Acad Sci U S A 100: 8484–8489.MC WolfgangBR KulasekaraX. LiangD. BoydK. Wu2003Conservation of genome content and virulence determinants among clinical and environmental isolates of Pseudomonas aeruginosa.Proc Natl Acad Sci U S A10084848489

18. He J, Baldini RL, Déziel E, Saucier M, Zhang Q, et al. (2004) The broad host range pathogen Pseudomonas aeruginosa strain PA14 carries two pathogenicity islands harboring plant and animal virulence genes. Proc Natl Acad Sci U S A 101: 2530–2535.J. HeRL BaldiniE. DézielM. SaucierQ. Zhang2004The broad host range pathogen Pseudomonas aeruginosa strain PA14 carries two pathogenicity islands harboring plant and animal virulence genes.Proc Natl Acad Sci U S A10125302535

19. Morales G, Wiehlmann L, Gudowius P, van Delden C, Tümmler B, et al. (2004) Structure of Pseudomonas aeruginosa populations analyzed by single nucleotide polymorphism and pulsed-field gel electrophoresis genotyping. J Bacteriol 186: 4228–4237.G. MoralesL. WiehlmannP. GudowiusC. van DeldenB. Tümmler2004Structure of Pseudomonas aeruginosa populations analyzed by single nucleotide polymorphism and pulsed-field gel electrophoresis genotyping.J Bacteriol18642284237

20. Mathee K, Narasimhan G, Valdes C, Qiu X, Matewish JM, et al. (2008) Dynamics of Pseudomonas aeruginosa genome evolution. Proc Natl Acad Sci U S A 105: 3100–3105.K. MatheeG. NarasimhanC. ValdesX. QiuJM Matewish2008Dynamics of Pseudomonas aeruginosa genome evolution.Proc Natl Acad Sci U S A10531003105

21. Shen K, Sayeed S, Antalis P, Gladitz J, Ahmed A, et al. (2006) Extensive genomic plasticity in Pseudomonas aeruginosa revealed by identification and distribution studies of novel genes among clinical isolates. Infect Immun 74: 5272–5283.K. ShenS. SayeedP. AntalisJ. GladitzA. Ahmed2006Extensive genomic plasticity in Pseudomonas aeruginosa revealed by identification and distribution studies of novel genes among clinical isolates.Infect Immun7452725283

22. Ripp S, Ogunseitan OA, Miller RV (1994) Transduction of a freshwater microbial community by a new Pseudomonas aeruginosa generalized transducing phage, UT1. Mol Ecol 3: 121–126.S. RippOA OgunseitanRV Miller1994Transduction of a freshwater microbial community by a new Pseudomonas aeruginosa generalized transducing phage, UT1.Mol Ecol3121126

23. Lee DG, Urbach JM, Wu G, Liberati NT, Feinbaum RL, et al. (2006) Genomic analysis reveals that Pseudomonas aeruginosa virulence is combinatorial. Genome Biol 7: R90.DG LeeJM UrbachG. WuNT LiberatiRL Feinbaum2006Genomic analysis reveals that Pseudomonas aeruginosa virulence is combinatorial.Genome Biol7R90

24. Wiehlmann L, Wagner G, Cramer N, Siebert B, Gudowius P, et al. (2007) Population structure of Pseudomonas aeruginosa. Proc Natl Acad Sci U S A 104: 8101–8106.L. WiehlmannG. WagnerN. CramerB. SiebertP. Gudowius2007Population structure of Pseudomonas aeruginosa.Proc Natl Acad Sci U S A10481018106

25. Farmer JJ 3rd, Weinstein RA, Zierdt CH, Brokopp CD (1982) Hospital outbreaks caused by Pseudomonas aeruginosa: Importance of serogroup O11. J Clin Microbiol 16: 266–270.JJ Farmer 3rdRA WeinsteinCH ZierdtCD Brokopp1982Hospital outbreaks caused by Pseudomonas aeruginosa: Importance of serogroup O11.J Clin Microbiol16266270

26. Giammanco A, Di Stefano R, Arista S, Sinatra A, Chiarini A (1985) Infections caused by Pseudomonas aeruginosa: relatively frequent isolation of serogroup 12 from clinical specimens. Eur J Epidemiol 1: 104–109.A. GiammancoR. Di StefanoS. AristaA. SinatraA. Chiarini1985Infections caused by Pseudomonas aeruginosa: relatively frequent isolation of serogroup 12 from clinical specimens.Eur J Epidemiol1104109

27. Legakis NJ, Koukoubanis N, Malliara K, Michalitsianos D, Papavassiliou J (1987) Importance of carbenicillin and gentamicin cross-resistant serotype O:12 Pseudomonas aeruginosa in six Athens hospitals. Eur J Clin Microbiol 6: 300–303.NJ LegakisN. KoukoubanisK. MalliaraD. MichalitsianosJ. Papavassiliou1987Importance of carbenicillin and gentamicin cross-resistant serotype O:12 Pseudomonas aeruginosa in six Athens hospitals.Eur J Clin Microbiol6300303

28. Pitt TL, Livermore DM, Pitcher D, Vatopoulos AC, Legakis NJ (1989) Multiresistant serotype O12 Pseudomonas aeruginosa: evidence for a common strain in Europe. Epidemiol Infect 103: 565–576.TL PittDM LivermoreD. PitcherAC VatopoulosNJ Legakis1989Multiresistant serotype O12 Pseudomonas aeruginosa: evidence for a common strain in Europe.Epidemiol Infect103565576

29. Grattard F, Gaudin OG, Pozzetto B, Ros A, Mbida AD (1993) Genotypic homogeneity of nosocomial Pseudomonas aeruginosa O12 strains demostrated by analysis of protein profiles, DNA fingerprints and rRNA gene restriction patterns. Eur J Clin Microbiol Infect Dis 12: 57–61.F. GrattardOG GaudinB. PozzettoA. RosAD Mbida1993Genotypic homogeneity of nosocomial Pseudomonas aeruginosa O12 strains demostrated by analysis of protein profiles, DNA fingerprints and rRNA gene restriction patterns.Eur J Clin Microbiol Infect Dis125761

30. Tsakris A, Vatopoulos AC, Tzouvelekis LS, Legakis NJ (1992) Diversity of resistance phenotypes and plasmid analysis in multi-resistant O:12 Pseudomonas aeruginosa. Eur J Epidemiol 8: 865–870.A. TsakrisAC VatopoulosLS TzouvelekisNJ Legakis1992Diversity of resistance phenotypes and plasmid analysis in multi-resistant O:12 Pseudomonas aeruginosa.Eur J Epidemiol8865870

31. Elaichouni A, Verschraegen G, Claeys G, Devleeschouwer M, Godard C, et al. (1994) Pseudomonas aeruginosa serotype O12 outbreak studied by arbitrary primer PCR. J Clin Microbiol 32: 666–671.A. ElaichouniG. VerschraegenG. ClaeysM. DevleeschouwerC. Godard1994Pseudomonas aeruginosa serotype O12 outbreak studied by arbitrary primer PCR.J Clin Microbiol32666671

32. Richard P, Le Floch R, Chamoux C, Pannier M, Espaze E, et al. (1994) Pseudomonas aeruginosa outbreak in a burn unit: role of antimicrobials in the emergence of multiply resistant strains. J Infect Dis 170: 377–383.P. RichardR. Le FlochC. ChamouxM. PannierE. Espaze1994Pseudomonas aeruginosa outbreak in a burn unit: role of antimicrobials in the emergence of multiply resistant strains.J Infect Dis170377383

33. Kettner M, Milosovic P, Hletková M, Kallová J (1995) Incidence and mechanisms of aminoglycoside resistance in Pseudomonas aeruginosa serotype O11 isolates. Infection 23: 380–383.M. KettnerP. MilosovicM. HletkováJ. Kallová1995Incidence and mechanisms of aminoglycoside resistance in Pseudomonas aeruginosa serotype O11 isolates.Infection23380383

34. Bingen E, Bonacorsi S, Rohrlich P, Duval M, Lhopital S, et al. (1996) Molecular evidence of genotypic heterogeneity of multidrug-resistant Pseudomonas aeruginosa serotype O:12 outbreak isolates from a pediatric hospital. J Clin Microbiol 34: 3226–3229.E. BingenS. BonacorsiP. RohrlichM. DuvalS. Lhopital1996Molecular evidence of genotypic heterogeneity of multidrug-resistant Pseudomonas aeruginosa serotype O:12 outbreak isolates from a pediatric hospital.J Clin Microbiol3432263229

35. Mifsud AJ, Watine J, Picard B, Charet JC, Solignac-Bourrel C, et al. (1997) Epidemiologically related and unrelated strains of Pseudomonas aeruginosa serotype O12 cannot be distinguished by phenotypic and genotypic typing. J Hosp Infect 36: 105–116.AJ MifsudJ. WatineB. PicardJC CharetC. Solignac-Bourrel1997Epidemiologically related and unrelated strains of Pseudomonas aeruginosa serotype O12 cannot be distinguished by phenotypic and genotypic typing.J Hosp Infect36105116

36. Tassios PT, Gennimata V, Maniatis AN, Fock C, Legakis NJ, et al. (1998) Emergence of multidrug resistance in ubiquitous and dominant Pseudomonas aeruginosa serogroup O:11. J Clin Microbiol 36: 897–901.PT TassiosV. GennimataAN ManiatisC. FockNJ Legakis1998Emergence of multidrug resistance in ubiquitous and dominant Pseudomonas aeruginosa serogroup O:11.J Clin Microbiol36897901

37. Watine J, Hacini J, Vidal I (1999) Is there a connection between prolonged carriage and clonal hospital-to-hospital clonal spread of multiresistant Pseudomonas aeruginosa of the O12 serotype? Are the specific habits of the hospitals involved the cause? Pathol Biol (Paris) 47: 457–461.J. WatineJ. HaciniI. Vidal1999Is there a connection between prolonged carriage and clonal hospital-to-hospital clonal spread of multiresistant Pseudomonas aeruginosa of the O12 serotype? Are the specific habits of the hospitals involved the cause?Pathol Biol (Paris)47457461

38. Müller-Premru M, Gubina M (2000) Serotype, antimicrobial susceptibility and clone distribution of Pseudomonas aeruginosa in a university hospital. Zentralbl Bakteriol 289: 857–867.M. Müller-PremruM. Gubina2000Serotype, antimicrobial susceptibility and clone distribution of Pseudomonas aeruginosa in a university hospital.Zentralbl Bakteriol289857867

39. Dubois V, Arpin C, Melon M, Melon B, Andre C, et al. (2001) Nosocomial outbreak due to a multiresistant strain of Pseudomonas aeruginosa P12: efficacy of cefepime-amikacin therapy and analysis of β-lactam resistance. J Clin Microbiol 39: 2072–2078.V. DuboisC. ArpinM. MelonB. MelonC. Andre2001Nosocomial outbreak due to a multiresistant strain of Pseudomonas aeruginosa P12: efficacy of cefepime-amikacin therapy and analysis of β-lactam resistance.J Clin Microbiol3920722078

40. Loureiro MM, de Moraes BA, Mendonca VL, Quadra MR, Pinheiro GS, et al. (2002) Pseudomonas aeruginosa: study of antibiotic resistance and molecular typing in hospital infection cases in a neonatal intensive care unit from Rio de Janeiro City, Brazil. Mem Inst Oswaldo Cruz 97: 387–394.MM LoureiroBA de MoraesVL MendoncaMR QuadraGS Pinheiro2002Pseudomonas aeruginosa: study of antibiotic resistance and molecular typing in hospital infection cases in a neonatal intensive care unit from Rio de Janeiro City, Brazil.Mem Inst Oswaldo Cruz97387394

41. Pirnay JP, De Vos D, Cochez C, Bilocq F, Pirson J, et al. (2003) Molecular epidemiology of Pseudomonas aeruginosa colonization in a burn unit: persistence of a multidrug-resistant clone and a silver sulfadiazine-resistant clone. J Clin Microbiol 41: 1192–1202.JP PirnayD. De VosC. CochezF. BilocqJ. Pirson2003Molecular epidemiology of Pseudomonas aeruginosa colonization in a burn unit: persistence of a multidrug-resistant clone and a silver sulfadiazine-resistant clone.J Clin Microbiol4111921202

42. Crespo MP, Woodford N, Sinclair A, Kaufmann ME, Turton J, et al. (2004) Outbreak of carbapenem-resistant Pseudomonas aeruginosa producing VIM-8, a novel metallo-beta-lactamase, in a tertiary care center in Cali, Colombia. J Clin Microbiol 42: 5094–5101.MP CrespoN. WoodfordA. SinclairME KaufmannJ. Turton2004Outbreak of carbapenem-resistant Pseudomonas aeruginosa producing VIM-8, a novel metallo-beta-lactamase, in a tertiary care center in Cali, Colombia.J Clin Microbiol4250945101

43. Deplano A, Denis O, Poirel L, Hocquet D, Nonhoff C, et al. (2005) Molecular characterization of an epidemic clone of panantibiotic-resistant Pseudomonas aeruginosa. J Clin Microbiol 43: 1198–1204.A. DeplanoO. DenisL. PoirelD. HocquetC. Nonhoff2005Molecular characterization of an epidemic clone of panantibiotic-resistant Pseudomonas aeruginosa.J Clin Microbiol4311981204

44. Riccio ML, Pallecchi L, Docquier JD, Cresti S, Catania MR, et al. (2005) Clonal relatedness and conserved integron structures in epidemiologically unrelated Pseudomonas aeruginosa strains producing the VIM-1 metallo-β-lactamase from different Italian hospitals. Antimicrob Agents Chemother 49: 104–110.ML RiccioL. PallecchiJD DocquierS. CrestiMR Catania2005Clonal relatedness and conserved integron structures in epidemiologically unrelated Pseudomonas aeruginosa strains producing the VIM-1 metallo-β-lactamase from different Italian hospitals.Antimicrob Agents Chemother49104110

45. Sekiguchi J, Asagi T, Miyoshi-Akiyama T, Kasai A, Mizuguchi Y, et al. (2007) Outbreaks of multidrug-resistant Pseudomonas aeruginosa in community hospitals in Japan. J Clin Microbiol 45: 979–989.J. SekiguchiT. AsagiT. Miyoshi-AkiyamaA. KasaiY. Mizuguchi2007Outbreaks of multidrug-resistant Pseudomonas aeruginosa in community hospitals in Japan.J Clin Microbiol45979989

46. Hammami S, Ghozzi R, Burghoffer B, Arlet G, Redjeb S (2008) Mechanisms of carbapenem resistance in non-metallo-beta-lactamase-producing clinical isolates of Pseudomonas aeruginosa from a Tunisian hospital. Pathol Biol (Paris) Oct 30 [Epub ahead of print]. S. HammamiR. GhozziB. BurghofferG. ArletS. Redjeb2008Mechanisms of carbapenem resistance in non-metallo-beta-lactamase-producing clinical isolates of Pseudomonas aeruginosa from a Tunisian hospital.Pathol Biol (Paris) Oct 30 [Epub ahead of print]

47. Libisch B, Watine J, Balogh B, Gacs M, Muzslay M, et al. (2008) Molecular typing indicates an important role for two international clonal complexes in dissemination of VIM-producing Pseudomonas aeruginosa clinical isolates in Hungary. Res Microbiol 159: 162–168.B. LibischJ. WatineB. BaloghM. GacsM. Muzslay2008Molecular typing indicates an important role for two international clonal complexes in dissemination of VIM-producing Pseudomonas aeruginosa clinical isolates in Hungary.Res Microbiol159162168

48. Cheng K, Smyth RL, Govan JRW, Doherty C, Winstanley C, et al. (1996) Spread of β-lactam-resistant Pseudomonas aeruginosa in a cystic fibrosis clinic. Lancet 348: 639–642.K. ChengRL SmythJRW GovanC. DohertyC. Winstanley1996Spread of β-lactam-resistant Pseudomonas aeruginosa in a cystic fibrosis clinic.Lancet348639642

49. Fluge G, Ojeniyi B, Høiby N, Digranes A, Ciofu O, et al. (2001) Typing of Pseudomonas aeruginosa strains in Norwegian cystic fibrosis patients. Clin Microbiol Infect 7: 238–243.G. FlugeB. OjeniyiN. HøibyA. DigranesO. Ciofu2001Typing of Pseudomonas aeruginosa strains in Norwegian cystic fibrosis patients.Clin Microbiol Infect7238243

50. Jones AM, Govan JRW, Doherty CJ, Dodd ME, Isalska BJ, et al. (2001) Spread of a multiresistant strain of Pseudomonas aeruginosa in an adult cystic fibrosis clinic. Lancet 358: 557–558.AM JonesJRW GovanCJ DohertyME DoddBJ Isalska2001Spread of a multiresistant strain of Pseudomonas aeruginosa in an adult cystic fibrosis clinic.Lancet358557558

51. McCallum SJ, Corkill J, Gallagher M, Ledson MJ, Hart CA, Walshaw MJ (2001) Superinfection with a transmissible strain of Pseudomonas aeruginosa in adults with cystic fibrosis chronically colonised by P. aeruginosa. Lancet 358: 558–560.SJ McCallumJ. CorkillM. GallagherMJ LedsonCA HartMJ Walshaw2001Superinfection with a transmissible strain of Pseudomonas aeruginosa in adults with cystic fibrosis chronically colonised by P. aeruginosa.Lancet358558560

52. Anthony M, Rose B, Pegler MB, Elkins M, Service H, et al. (2002) Genetic analysis of Pseudomonas aeruginosa isolates from the sputa of Australian adult cystic fibrosis patients. J Clin Microbiol 40: 2772–2778.M. AnthonyB. RoseMB PeglerM. ElkinsService H, et al.2002Genetic analysis of Pseudomonas aeruginosa isolates from the sputa of Australian adult cystic fibrosis patients.J Clin Microbiol4027722778

53. Armstrong D, Bell S, Robinson M, Bye P, Rose B, et al. (2003) Evidence for spread of a clonal strain of Pseudomonas aeruginosa among cystic fibrosis clinics. J Clin Microbiol 41: 2266–2267.D. ArmstrongS. BellM. RobinsonP. ByeB. Rose2003Evidence for spread of a clonal strain of Pseudomonas aeruginosa among cystic fibrosis clinics.J Clin Microbiol4122662267

54. O'Carroll MR, Syrmis MW, Wainwright CE, Greer RM, Mitchell P, et al. (2004) Clonal strains of Pseudomonas aeruginosa in paediatric and adult cystic fibrosis units. Eur Respir J 24: 101–106.MR O'CarrollMW SyrmisCE WainwrightRM GreerP. Mitchell2004Clonal strains of Pseudomonas aeruginosa in paediatric and adult cystic fibrosis units.Eur Respir J24101106

55. Scott F, Pitt TL (2004) Identification and characterization of transmissible Pseudomonas aeruginosa strains in cystic fibrosis patients in England and Wales. J Med Microbiol 53: 609–615.F. ScottTL Pitt2004Identification and characterization of transmissible Pseudomonas aeruginosa strains in cystic fibrosis patients in England and Wales.J Med Microbiol53609615

56. Chambers D, Scott F, Bangur R, Davies R, Lim A, et al. (2005) Factors associated with infection by Pseudomonas aeruginosa in adult cystic fibrosis. Eur Respir J 26: 651–656.D. ChambersF. ScottR. BangurR. DaviesA. Lim2005Factors associated with infection by Pseudomonas aeruginosa in adult cystic fibrosis.Eur Respir J26651656

57. Lewis DA, Jones A, Parkhill J, Speert DP, Govan JRW, et al. (2005) Identification of DNA markers for a transmissible Pseudomonas aeruginosa cystic fibrosis strain. Am J Respir Cell Mol Biol 33: 56–64.DA LewisA. JonesJ. ParkhillDP SpeertJRW Govan2005Identification of DNA markers for a transmissible Pseudomonas aeruginosa cystic fibrosis strain.Am J Respir Cell Mol Biol335664

58. Jelsbak L, Johansen HK, Frost AL, Thøgersen R, Thomsen LE, et al. (2007) Molecular epidemiology and dynamics of Pseudomonas aeruginosa populations in lungs of cystic fibrosis patients. Infect Immun 75: 2214–2224.L. JelsbakHK JohansenAL FrostR. ThøgersenLE Thomsen2007Molecular epidemiology and dynamics of Pseudomonas aeruginosa populations in lungs of cystic fibrosis patients.Infect Immun7522142224

59. Bradbury R, Champion A, Reid DW (2008) Poor clinical outcomes associated with a multi-drug resistant clonal strain of Pseudomonas aeruginosa in the Tasmanian cystic fibrosis population. Respirology 13: 886–892.R. BradburyA. ChampionDW Reid2008Poor clinical outcomes associated with a multi-drug resistant clonal strain of Pseudomonas aeruginosa in the Tasmanian cystic fibrosis population.Respirology13886892

60. Brimicombe RW, Dijkshoorn L, van der Reijden TJ, Kardoes I, Pitt TL, et al. (2008) Transmission of Pseudomonas aeruginosa in children with cystic fibrosis attending summer camps in The Netherlands. J Cyst Fibros 7: 30–36.RW BrimicombeL. DijkshoornTJ van der ReijdenI. KardoesTL Pitt2008Transmission of Pseudomonas aeruginosa in children with cystic fibrosis attending summer camps in The Netherlands.J Cyst Fibros73036

61. van Belkum A (1996) Current trends in typing of bacterial strains for medical purposes. Zentralbl Bakteriol 1045: 249–252.A. van Belkum1996Current trends in typing of bacterial strains for medical purposes.Zentralbl Bakteriol1045249252

62. Vandamme P, Pot B, Gillis M, de Vos P, Kersters K, et al. (1996) Polyphasic taxonomy, a consensus approach to bacterial systematics. Microbiol Rev 60: 407–438.P. VandammeB. PotM. GillisP. de VosK. Kersters1996Polyphasic taxonomy, a consensus approach to bacterial systematics.Microbiol Rev60407438

63. Aoki K (1926) Agglutinatorische einteilung von Pyocyaneus-bazillen welche bei verschiedenen menschenerkrangungen nach gewiesen wurden. Zentralbl Bakteriol Parasitenkd Infectionskr Hyg Abt 1 Orig 98: 186–195.K. Aoki1926Agglutinatorische einteilung von Pyocyaneus-bazillen welche bei verschiedenen menschenerkrangungen nach gewiesen wurden.Zentralbl Bakteriol Parasitenkd Infectionskr Hyg Abt 1 Orig98186195

64. Vos P, Hogers R, Bleeker M, Reijans M, van de Lee T, et al. (1995) AFLP: a new technique for DNA fingerprinting. Nucleic Acids Res 23: 4407–4414.P. VosR. HogersM. BleekerM. ReijansT. van de Lee1995AFLP: a new technique for DNA fingerprinting.Nucleic Acids Res2344074414

65. Janssen P, Coopman R, Huys G, Swings J, Bleeker M, et al. (1996) Evaluation of the DNA fingerprinting method AFLP as a new tool in bacterial taxonomy. Microbiology 142: 1881–1893.P. JanssenR. CoopmanG. HuysJ. SwingsM. Bleeker1996Evaluation of the DNA fingerprinting method AFLP as a new tool in bacterial taxonomy.Microbiology14218811893

66. Savelkoul PHM, Aarts HJM, de Haas J, Dijkshoorn L, Duim B, et al. (1999) Amplified-fragment length polymorphism analysis: the state of an art. J Clin Microbiol 37: 3083–3091.PHM SavelkoulHJM AartsJ. de HaasL. DijkshoornB. Duim1999Amplified-fragment length polymorphism analysis: the state of an art.J Clin Microbiol3730833091

67. Spratt BG, Maiden MCJ (1999) Bacterial population genetics, evolution and epidemiology. Phil Trans R Soc Lond 354: 701–710.BG SprattMCJ Maiden1999Bacterial population genetics, evolution and epidemiology.Phil Trans R Soc Lond354701710

68. Cornelis P, Bouia A, Belarbi A, Guyonvarch A, Kammerer B, et al. (1989) Cloning and analysis of the gene for the major outer membrane lipoprotein from Pseudomonas aeruginosa. Mol Microbiol 3: 421–428.P. CornelisA. BouiaA. BelarbiA. GuyonvarchB. Kammerer1989Cloning and analysis of the gene for the major outer membrane lipoprotein from Pseudomonas aeruginosa.Mol Microbiol3421428

69. Lim A Jr, De Vos D, Brauns M, Mossialos D, Gaballa A, et al. (1997) Molecular and immunological characterization of OprL, the 18 kDa outer-membrane peptidoglycan-associated lipoprotein (PAL) of Pseudomonas aeruginosa. Microbiology 143: 1709–1716.A. Lim JrD. De VosM. BraunsD. MossialosA. Gaballa1997Molecular and immunological characterization of OprL, the 18 kDa outer-membrane peptidoglycan-associated lipoprotein (PAL) of Pseudomonas aeruginosa.Microbiology14317091716

70. De Vos D, Lim A Jr, Pirnay JP, Struelens M, Vandenvelde C, et al. (1997) Direct detection and identification of Pseudomonas aeruginosa in clinical samples such as skin biopsy specimens and expectorations by multiplex PCR based on two outer membrane lipoprotein genes, oprI and oprL. J Clin Microbiol 35: 1295–1299.D. De VosA. Lim JrJP PirnayM. StruelensC. Vandenvelde1997Direct detection and identification of Pseudomonas aeruginosa in clinical samples such as skin biopsy specimens and expectorations by multiplex PCR based on two outer membrane lipoprotein genes, oprI and oprL.J Clin Microbiol3512951299

71. De Vos D, Bouton C, Sarniguet A, De Vos P, Vauterin M, et al. (1998) Sequence diversity of the oprI gene, coding for major outer membrane lipoprotein I, among rRNA group I pseudomonads. J Bacteriol 180: 6551–6556.D. De VosC. BoutonA. SarniguetP. De VosM. Vauterin1998Sequence diversity of the oprI gene, coding for major outer membrane lipoprotein I, among rRNA group I pseudomonads.J Bacteriol18065516556

72. Trias J, Nikaido H (1990) Outer membrane protein D2 catalyzes facilitated diffusion of carbapenems and penems through the outer membrane of Pseudomonas aeruginosa. Antimicrob Agents Chemother 34: 52–57.J. TriasH. Nikaido1990Outer membrane protein D2 catalyzes facilitated diffusion of carbapenems and penems through the outer membrane of Pseudomonas aeruginosa.Antimicrob Agents Chemother345257

73. Trias J, Nikaido H (1990) Protein D2 channel of the Pseudomonas aeruginosa outer membrane has a binding site for basic amino acids and peptides. J Biol Chem 265: 15680–15684.J. TriasH. Nikaido1990Protein D2 channel of the Pseudomonas aeruginosa outer membrane has a binding site for basic amino acids and peptides.J Biol Chem2651568015684

74. Pirnay JP, De Vos D, Mossialos D, Vanderkelen A, Cornelis P, et al. (2002) Analysis of the Pseudomonas aeruginosa oprD gene from clinical and environmental isolates. Environ Microbiol 4: 872–882.JP PirnayD. De VosD. MossialosA. VanderkelenP. Cornelis2002Analysis of the Pseudomonas aeruginosa oprD gene from clinical and environmental isolates.Environ Microbiol4872882

75. Meyer JM (2000) Pyoverdines: pigments, siderophores and potential taxonomic markers of fluorescent Pseudomonas species. Arch Microbiol 174: 135–142.JM Meyer2000Pyoverdines: pigments, siderophores and potential taxonomic markers of fluorescent Pseudomonas species.Arch Microbiol174135142

76. Cornelis P, Hohnadel D, Meyer JM (1989) Evidence for different pyoverdine-mediated iron uptake systems among Pseudomonas aeruginosa strains. Infect Immun 57: 3491–3497.P. CornelisD. HohnadelJM Meyer1989Evidence for different pyoverdine-mediated iron uptake systems among Pseudomonas aeruginosa strains.Infect Immun5734913497

77. Meyer JM, Stintzi A, De Vos D, Cornelis P, Tappe R, et al. (1997) Use of siderophores to type pseudomonads: the three Pseudomonas aeruginosa pyoverdine systems. Microbiology 143: 35–43.JM MeyerA. StintziD. De VosP. CornelisR. Tappe1997Use of siderophores to type pseudomonads: the three Pseudomonas aeruginosa pyoverdine systems.Microbiology1433543

78. de Chial M, Ghysels B, Beatson SA, Geoffroy V, Meyer JM, et al. (2003) Identification of type II and type III pyoverdine receptors from Pseudomonas aeruginosa. Microbiology 149: 821–831.M. de ChialB. GhyselsSA BeatsonV. GeoffroyJM Meyer2003Identification of type II and type III pyoverdine receptors from Pseudomonas aeruginosa.Microbiology149821831

79. Tümmler B, Cornelis P (2005) Pyoverdine receptor: a case of positive Darwinian selection in Pseudomonas aeruginosa. J Bacteriol 187: 3289–3292.B. TümmlerP. Cornelis2005Pyoverdine receptor: a case of positive Darwinian selection in Pseudomonas aeruginosa.J Bacteriol18732893292

80. Smith EE, Sims EH, Spencer DH, Kaul R, Olson MV (2005) Evidence for diversifying selection at the pyoverdine locus of Pseudomonas aeruginosa. J Bacteriol 187: 2138–2147.EE SmithEH SimsDH SpencerR. KaulMV Olson2005Evidence for diversifying selection at the pyoverdine locus of Pseudomonas aeruginosa.J Bacteriol18721382147

81. Bodilis J, Ghysels B, Osayande J, Matthijs S, Pirnay JP, et al. (2009) Distribution and evolution of ferripyoverdine receptors in Pseudomonas aeruginosa. Environ Microbiol. 2009 Apr 21. [Epub ahead of print]. J. BodilisB. GhyselsJ. OsayandeS. MatthijsJP Pirnay2009Distribution and evolution of ferripyoverdine receptors in Pseudomonas aeruginosa.Environ Microbiol. 2009 Apr 21. [Epub ahead of print]

82. Ghysels B, Dieu BT, Beatson SA, Pirnay JP, Ochsner UA, et al. (2004) FpvB, an alternative type I ferripyoverdine receptor of Pseudomonas aeruginosa. Microbiology 150: 1671–1680.B. GhyselsBT DieuSA BeatsonJP PirnayUA Ochsner2004FpvB, an alternative type I ferripyoverdine receptor of Pseudomonas aeruginosa.Microbiology15016711680

83. Schulert GS, Feltman H, Rabin SD, Martin CG, Battle SE, et al. (2003) Secretion of the toxin ExoU is a marker for highly virulent Pseudomonas aeruginosa isolates obtained from patients with hospital-acquired pneumonia. J Infect Dis 188: 1695–1706.GS SchulertH. FeltmanSD RabinCG MartinSE Battle2003Secretion of the toxin ExoU is a marker for highly virulent Pseudomonas aeruginosa isolates obtained from patients with hospital-acquired pneumonia.J Infect Dis18816951706

84. Lin HH, Huang SP, Teng HC, Ji DD, Chen YS, et al. (2006) Presence of the exoU gene of Pseudomonas aeruginosa is correlated with cytotoxicity in MDCK cells but not with colonization in BALB/c mice. J Clin Microbiol 44: 4596–4597.HH LinSP HuangHC TengDD JiYS Chen2006Presence of the exoU gene of Pseudomonas aeruginosa is correlated with cytotoxicity in MDCK cells but not with colonization in BALB/c mice.J Clin Microbiol4445964597

85. Castric P (1995) pilO, a gene required for glycosylation of Pseudomonas aeruginosa 1244 pilin. Microbiology 141: 1247–1254.P. Castric1995pilO, a gene required for glycosylation of Pseudomonas aeruginosa 1244 pilin.Microbiology14112471254

86. Kus JV, Tullis E, Cvitkovitch DG, Burrows LL (2004) Significant differences in type IV pilin allele distribution among Pseudomonas aeruginosa isolates from cystic fibrosis (CF) versus non-CF patients. Microbiology 150: 1315–1326.JV KusE. TullisDG CvitkovitchLL Burrows2004Significant differences in type IV pilin allele distribution among Pseudomonas aeruginosa isolates from cystic fibrosis (CF) versus non-CF patients.Microbiology15013151326

87. Hancock RE, Mutharia LM, Chan L, Darveau RP, Speert DP, et al. (1983) Pseudomonas aeruginosa isolates from patients with cystic fibrosis: a class of serum-sensitive, nontypable strains deficient in lipopolysaccharide O side chains. Infect Immun 42: 170–177.RE HancockLM MuthariaL. ChanRP DarveauDP Speert1983Pseudomonas aeruginosa isolates from patients with cystic fibrosis: a class of serum-sensitive, nontypable strains deficient in lipopolysaccharide O side chains.Infect Immun42170177

88. Pitt TL, MacDougall J, Penketh AR, Cooke EM (1986) Polyagglutinating and non-typable strains of Pseudomonas aeruginosa in cystic fibrosis. J Med Microbiol 21: 179–186.TL PittJ. MacDougallAR PenkethEM Cooke1986Polyagglutinating and non-typable strains of Pseudomonas aeruginosa in cystic fibrosis.J Med Microbiol21179186

89. Penketh A, Pitt T, Roberts D, Hodson ME, Batten JC (1983) The relationship of phenotype changes in Pseudomonas aeruginosa to the clinical condition of patients with cystic fibrosis. Am Rev Respir Dis 127: 605–608.A. PenkethT. PittD. RobertsME HodsonJC Batten1983The relationship of phenotype changes in Pseudomonas aeruginosa to the clinical condition of patients with cystic fibrosis.Am Rev Respir Dis127605608

90. Zierdt CH, Williams RL (1975) Serotyping of Pseudomonas aeruginosa isolates from patients with cystic fibrosis of the pancreas. J Clin Microbiol 1: 521–526.CH ZierdtRL Williams1975Serotyping of Pseudomonas aeruginosa isolates from patients with cystic fibrosis of the pancreas.J Clin Microbiol1521526

91. Bert F, Lambert-Zechovsky N (1996) Comparative distribution of resistance patterns and serotypes in Pseudomonas aeruginosa isolates from intensive care units and other wards. J Antimicrob Chemother 37: 809–813.F. BertN. Lambert-Zechovsky1996Comparative distribution of resistance patterns and serotypes in Pseudomonas aeruginosa isolates from intensive care units and other wards.J Antimicrob Chemother37809813

92. Kobayashi L, Hasegawa M, Miyazaki S, Nishida M, Goto S (1994) In vitro and in vivo changes of serotype in Pseudomonas aeruginosa isolates by anti-pseudomonal drugs. J Antibiot (Tokyo) 47: 72–79.L. KobayashiM. HasegawaS. MiyazakiM. NishidaS. Goto1994In vitro and in vivo changes of serotype in Pseudomonas aeruginosa isolates by anti-pseudomonal drugs.J Antibiot (Tokyo)477279

93. Kuzio J, Kropinski AM (1983) O-antigen conversion in Pseudomonas aeruginosa PAO1 by bacteriophage D3. J Bacteriol 155: 203–212.J. KuzioAM Kropinski1983O-antigen conversion in Pseudomonas aeruginosa PAO1 by bacteriophage D3.J Bacteriol155203212

94. Speijer H, Savelkoul PH, Bonten MJ, Stobberingh EE, Tjhie JH (1999) Application of different genotyping methods for Pseudomonas aeruginosa in a setting of endemicity in an intensive care unit. J Clin Microbiol 37: 3654–3661.H. SpeijerPH SavelkoulMJ BontenEE StobberinghJH Tjhie1999Application of different genotyping methods for Pseudomonas aeruginosa in a setting of endemicity in an intensive care unit.J Clin Microbiol3736543661

95. Tenover FC, Arbeit RD, Goering RV, Micklesen PA, Murray BE, et al. (1995) Interpreting chromosomal DNA restriction patterns produced by pulsed-field gel electrophoresis: criteria for bacterial strain typing. J Clin Microbiol 33: 2233–2239.FC TenoverRD ArbeitRV GoeringPA MicklesenBE Murray1995Interpreting chromosomal DNA restriction patterns produced by pulsed-field gel electrophoresis: criteria for bacterial strain typing.J Clin Microbiol3322332239

96. Fleiszig SM, Wiener-Kronish JP, Miyazaki H, Vallas V, Mostov KE, et al. (1997) Pseudomonas aeruginosa-mediated cytotoxicity and invasion correlate with distinct genotypes at the loci encoding exoenzyme S. Infect Immun 65: 579–586.SM FleiszigJP Wiener-KronishH. MiyazakiV. VallasKE Mostov1997Pseudomonas aeruginosa-mediated cytotoxicity and invasion correlate with distinct genotypes at the loci encoding exoenzyme S.Infect Immun65579586

97. Hauser AR, Cobb E, Bodi M, Mariscal D, Vallés J, et al. (2002) Type III protein secretion is associated with poor clinical outcomes in patients with ventilator-associated pneumonia caused by Pseudomonas aeruginosa. Crit Care Med 30: 521–528.AR HauserE. CobbM. BodiD. MariscalJ. Vallés2002Type III protein secretion is associated with poor clinical outcomes in patients with ventilator-associated pneumonia caused by Pseudomonas aeruginosa.Crit Care Med30521528

98. Feltman H, Schulert G, Khan S, Jain M, Peterson L, et al. (2001) Prevalence of type III secretion genes in clinical and environmental isolates of Pseudomonas aeruginosa. Microbiology 147: 2659–2669.H. FeltmanG. SchulertS. KhanM. JainL. Peterson2001Prevalence of type III secretion genes in clinical and environmental isolates of Pseudomonas aeruginosa.Microbiology14726592669

99. Wareham DW, Curtis MA (2007) A genotypic and phenotypic comparison of type III secretion profiles of Pseudomonas aeruginosa cystic fibrosis and bacteremia isolates. Int J Med Microbiol 297: 227–234.DW WarehamMA Curtis2007A genotypic and phenotypic comparison of type III secretion profiles of Pseudomonas aeruginosa cystic fibrosis and bacteremia isolates.Int J Med Microbiol297227234

100. Ferguson MW, Maxwell JA, Vincent TS, da Silva J, Olson JC (2001) Comparison of the exoS gene and protein expression in soil and clinical isolates of Pseudomonas aeruginosa. Infect Immun 69: 2198–2210.MW FergusonJA MaxwellTS VincentJ. da SilvaJC Olson2001Comparison of the exoS gene and protein expression in soil and clinical isolates of Pseudomonas aeruginosa.Infect Immun6921982210

101. Olson JC, Fraylick JE, McGuffie EM, Dolan KM, Yahr TL, et al. (1999) Interruption of multiple cellular processes in HT-29 epithelial cells by Pseudomonas aeruginosa exoenzyme S. Infect Immun 67: 2847–2854.JC OlsonJE FraylickEM McGuffieKM DolanTL Yahr1999Interruption of multiple cellular processes in HT-29 epithelial cells by Pseudomonas aeruginosa exoenzyme S.Infect Immun6728472854

102. Rocha CL, Coburn J, Rucks EA, Olson JC (2003) Characterization of Pseudomonas aeruginosa exoenzyme S as a bifunctional enzyme in J774A.1 macrophages. Infect Immun 71: 5296–5305.CL RochaJ. CoburnEA RucksJC Olson2003Characterization of Pseudomonas aeruginosa exoenzyme S as a bifunctional enzyme in J774A.1 macrophages.Infect Immun7152965305

103. Gallant CV, Raivio TL, Olson JC, Woods DE, Storey DG (2000) Pseudomonas aeruginosa cystic fibrosis clinical isolates produce exotoxin A with altered ADP-ribosyltransferase activity and cytotoxicity. Microbiology 146: 1891–1899.CV GallantTL RaivioJC OlsonDE WoodsDG Storey2000Pseudomonas aeruginosa cystic fibrosis clinical isolates produce exotoxin A with altered ADP-ribosyltransferase activity and cytotoxicity.Microbiology14618911899

104. Kulasekara BR, Kulasekara HD, Wolfgang MC, Stevens L, Frank DW, et al. (2006) Acquisition and evolution of the exoU locus in Pseudomonas aeruginosa. J Bacteriol 188: 4037–4050.BR KulasekaraHD KulasekaraMC WolfgangL. StevensDW Frank2006Acquisition and evolution of the exoU locus in Pseudomonas aeruginosa.J Bacteriol18840374050

105. De Vos D, de Chial M, Cochez C, Jansen S, Tümmler B, et al. (2001) Study of pyoverdine type and production by Pseudomonas aeruginosa isolated from cystic fibrosis patients: prevalence of type II pyoverdine isolates and accumulation pyoverdine-negative mutations. Arch Microbiol 175: 384–388.D. De VosM. de ChialC. CochezS. JansenB. Tümmler2001Study of pyoverdine type and production by Pseudomonas aeruginosa isolated from cystic fibrosis patients: prevalence of type II pyoverdine isolates and accumulation pyoverdine-negative mutations.Arch Microbiol175384388

106. Baysse C, Meyer JM, Plesiat P, Geoffroy V, Michel-Briand Y, et al. (1999) Uptake of pyocin S3 occurs through the outer membrane ferripyoverdine type II receptor of Pseudomonas aeruginosa. J Bacteriol 181: 3849–3851.C. BaysseJM MeyerP. PlesiatV. GeoffroyY. Michel-Briand1999Uptake of pyocin S3 occurs through the outer membrane ferripyoverdine type II receptor of Pseudomonas aeruginosa.J Bacteriol18138493851

107. Denayer S, Matthijs S, Cornelis P (2007) Pyocin S2 (Sa) kills Pseudomonas aeruginosa strains via the FpvA type I ferripyoverdine receptor. J Bacteriol 189: 7663–7668.S. DenayerS. MatthijsP. Cornelis2007Pyocin S2 (Sa) kills Pseudomonas aeruginosa strains via the FpvA type I ferripyoverdine receptor.J Bacteriol18976637668

108. Vauterin L, Vauterin P (2006) Integrated databasing and analysis. In: Stackebrandt E, editor. Molecular identification, systematics, and population structure of prokaryotes. Berlin: Springer-Verlag. pp. 141–217.L. VauterinP. Vauterin2006Integrated databasing and analysis.E. StackebrandtMolecular identification, systematics, and population structure of prokaryotesBerlinSpringer-Verlag141217

109. Schouls LM, van der Heide HG, Vauterin L, Vauterin P, Mooi FR (2004) Multiple-locus variable-number tandem repeat analysis of Dutch Bordetella pertussis strains reveals rapid genetic changes with clonal expansion during the late 1990s. J Bacteriol 186: 5496–5505.LM SchoulsHG van der HeideL. VauterinP. VauterinFR Mooi2004Multiple-locus variable-number tandem repeat analysis of Dutch Bordetella pertussis strains reveals rapid genetic changes with clonal expansion during the late 1990s.J Bacteriol18654965505

110. Ragon M, Wirth T, Hollandt F, Lavenir R, Lecuit M, et al. (2008) A new perspective on Listeria monocytogenes evolution. PLoS Pathog 4: e1000146.M. RagonT. WirthF. HollandtR. LavenirM. Lecuit2008A new perspective on Listeria monocytogenes evolution.PLoS Pathog4e1000146

111. Lanotte P, Watt S, Mereghetti L, Dartiguelongue N, Rastegar-Lari A, et al. (2004) Genetic features of Pseudomonas aeruginosa isolates from cystic fibrosis patients compared with those of isolates from other origins. J Med Microbiol 53: 73–81.P. LanotteS. WattL. MereghettiN. DartiguelongueA. Rastegar-Lari2004Genetic features of Pseudomonas aeruginosa isolates from cystic fibrosis patients compared with those of isolates from other origins.J Med Microbiol537381

112. Pitt TL (2002) Cross infection of cystic fibrosis patients with Pseudomonas aeruginosa. Thorax 57: 921.TL Pitt2002Cross infection of cystic fibrosis patients with Pseudomonas aeruginosa.Thorax57921

113. Van Daele SG, Franckx H, Verhelst R, Schelstraete P, Haerynck F, et al. (2005) Epidemiology of Pseudomonas aeruginosa in a cystic fibrosis rehabilitation centre. Eur Respir J 25: 474–481.SG Van DaeleH. FranckxR. VerhelstP. SchelstraeteF. Haerynck2005Epidemiology of Pseudomonas aeruginosa in a cystic fibrosis rehabilitation centre.Eur Respir J25474481

114. Allemeersch D, Beumer J, Devleeschouwer M, De Maeyer S, Dony J, et al. (1988) Marked increase of Pseudomonas aeruginosa serotype O12 in Belgium since 1982. Eur J Clin Microbiol Infect Dis 7: 265–269.D. AllemeerschJ. BeumerM. DevleeschouwerS. De MaeyerJ. Dony1988Marked increase of Pseudomonas aeruginosa serotype O12 in Belgium since 1982.Eur J Clin Microbiol Infect Dis7265269

115. Finlay BJ (2002) Global dispersal of free-living microbial eukaryote species. Science 296: 1061–1063.BJ Finlay2002Global dispersal of free-living microbial eukaryote species.Science29610611063

116. Réseau épidémiologique des utilisateurs du systeme SIR (2001) Are the prolonged carrier and interhospital clonal diffusion of serotype O12 multiresistant Pseudomonas aeruginosa connected? Pathol Biol (Paris) 49: 620–623.Réseau épidémiologique des utilisateurs du systeme SIR2001Are the prolonged carrier and interhospital clonal diffusion of serotype O12 multiresistant Pseudomonas aeruginosa connected?Pathol Biol (Paris)49620623

117. Liu PV, Matsumoto H, Kusama H, Bergan T (1983) Survey of heat-stable, major somatic antigens of Pseudomonas aeruginosa. Int J Syst Bacteriol 33: 256–264.PV LiuH. MatsumotoH. KusamaT. Bergan1983Survey of heat-stable, major somatic antigens of Pseudomonas aeruginosa.Int J Syst Bacteriol33256264

118. Livermore DM, Struelens M, Amorim J, Baquero F, Bille J, et al. (2002) Multicentre evaluation of the VITEK 2 Advanced Expert System for interpretive reading of antimicrobial resistance tests. J Antimicrob Chemother 49: 289–300.DM LivermoreM. StruelensJ. AmorimF. BaqueroJ. Bille2002Multicentre evaluation of the VITEK 2 Advanced Expert System for interpretive reading of antimicrobial resistance tests.J Antimicrob Chemother49289300

119. Amsterdam D (1991) Susceptibility testing of antimicrobials in liquid media. In: Lorian V, editor. Antibiotics in laboratory medicine. Baltimore: Williams & Wilkins. pp. 72–78.D. Amsterdam1991Susceptibility testing of antimicrobials in liquid media.V. LorianAntibiotics in laboratory medicineBaltimoreWilliams & Wilkins7278