Pseudomonas aeruginosa clinical and environmental isolates constitute a single population with high phenotypic diversity
- María-Victoria Grosso-Becerra†1,
- Christian Santos-Medellín†1,
- Abigail González-Valdez1,
- José-Luis Méndez2,
- Gabriela Delgado2,
- Rosario Morales-Espinosa2,
- Luis Servín-González1,
- Luis-David Alcaraz3 and
- Gloria Soberón-Chávez1Email author
© Grosso-Becerra et al.; licensee BioMed Central Ltd. 2014
Received: 24 July 2013
Accepted: 24 March 2014
Published: 28 April 2014
Pseudomonas aeruginosa is an opportunistic pathogen with a high incidence of hospital infections that represents a threat to immune compromised patients. Genomic studies have shown that, in contrast to other pathogenic bacteria, clinical and environmental isolates do not show particular genomic differences. In addition, genetic variability of all the P. aeruginosa strains whose genomes have been sequenced is extremely low. This low genomic variability might be explained if clinical strains constitute a subpopulation of this bacterial species present in environments that are close to human populations, which preferentially produce virulence associated traits.
In this work, we sequenced the genomes and performed phenotypic descriptions for four non-human P. aeruginosa isolates collected from a plant, the ocean, a water-spring, and from dolphin stomach. We show that the four strains are phenotypically diverse and that this is not reflected in genomic variability, since their genomes are almost identical. Furthermore, we performed a detailed comparative genomic analysis of the four strains studied in this work with the thirteen previously reported P. aeruginosa genomes by means of describing their core and pan-genomes.
Contrary to what has been described for other bacteria we have found that the P. aeruginosa core genome is constituted by a high proportion of genes and that its pan-genome is thus relatively small. Considering the high degree of genomic conservation between isolates of P. aeruginosa from diverse environments, including human tissues, some implications for the treatment of infections are discussed. This work also represents a methodological contribution for the genomic study of P. aeruginosa, since we provide a database of the comparison of all the proteins encoded by the seventeen strains analyzed.
KeywordsP. aeruginosa genomics and proteomics Core genome Pangenomics Phenotypic diversity
Pseudomonas aeruginosa is an environmental, ubiquitous γ-proteobacterium that is also an important opportunistic human pathogen . In contrast to other bacterial pathogens, P. aeruginosa genomes of clinical and environmental isolates are highly conserved, and all isolates are able to produce virulence-associated traits and are thus potential pathogens [2, 3].
The pathogenicity of this bacterium depends on the production and secretion of multiple virulence factors that are regulated at the transcriptional level by the so-called quorum sensing (QS) response . This genetic response is mediated by the bacterial production of acyl-homoserine lactones (autoinducers) that act as signal molecules interacting with transcriptional regulators of the LuxR family. Pseudomonas aeruginosa QS is a hierarchical regulatory cascade: LasR interacts with 3-oxo-dodecanoyl-homoserine lactone (3O-C12-HSL), produced by the LasI enzyme, and activates the transcription of lasI, and of several genes coding for virulence factors. It also activates transcription of rhlR, which encodes the second QS transcriptional regulator, and of rhlI, which encodes the enzyme that produces butanoyl-homoserine lactone (C4-HSL), which is the autoinducer that interacts with RhlR . RhlR/C4-HSL in turn promotes the expression of genes responsible for the production of several virulence factors. These include among others, pyocyanin, lectin PA-IL (encoded by lecA), and biosurfactant rhamnolipids (synthesized by the products of the rhlAB operon).
This bacterium represents an important public health problem, as one of the leading causes of nosocomial infections  and as an infectious agent for cystic fibrosis patients . Additionally, the high level of antibiotic resistance shown by this bacterium  makes it very difficult to treat P. aeruginosa infections.
The first P. aeruginosa genome to be completely sequenced was that of the type strain PAO1 , which was described in 1955, as an isolate from the wound of a patient in Australia . At present, the genome sequences of twelve clinical isolates and of a strain isolated from watermelon rhizosphere (M18) have been reported [NCBI data base, [8, 10–18]]. Strain M18 produces higher levels of phenazine-carboxylic acid (PCA), the precursor of pyocyanin, at 30°C than at 37°C , and the hierarchy of the autoinducer-based transcriptional regulation of virulence associated traits, i.e. the QS response, is different to that of the type strain PAO1 . However, the genomes of clinical and environmental isolates, including strain M18 are very highly conserved . This high sequence identity at the genome level is not observed in other bacterial species, such as Escherichia coli, or even in other Pseudomonas species . Furthermore, it has been reported that the sequenced genomes of P. aeruginosa strains, including PAO1, show more than 95% sequence identity in the genes of their core genomes. The only sequenced genome that shows a slightly lower sequence identity is that of the multiresistant strain PA7 isolated in Argentina, which belongs to a different clade .
General features of the P. aeruginosa strains used in this work
Genome Size (Mb)
Environment of isolation
Dolphin, gastric juice
Cystic fibrosis patient
Cornea from a patient with ulcerative keratitis
Isolate from an infant with community acquired diarrhea
Cystic fibrosis patient
Cystic fibrosis patient
Epidemic strain from Manchester England
Sweet melon rhizosphere
Urinary tract infection patient
Cystic fibrosis patient, nosocomial infection
Burn wound isolate
Non-respiratory human isolate
Human burn patient
Persistent isolate from a Cystic Fibrosis patient
The characterized environmental isolates are: strain ID4365, a marine isolate from the Indian Ocean which produces very high levels of pyoverdins and phenazines ; strain M10, that was isolated from the water spring of the Churince system at Cuatro Cienegas, Coahuila, Mexico (a system of ponds in the middle of the Chihuahuan Desert with features that resemble those of an ancient ocean ); and strain IGB83, which is a highly lipolytic strain isolated in the rainforest of the Mexican State of Chiapas from a rotten coconut . We also studied strain 148 (Table 1), since it was isolated from a dolphin kept in captivity in a marine aquarium in Cancún, Quintana Roo, Mexico and it enabled us to determine whether pathogenic interactions with mammals were possibly being selected in the P. aeruginosa genomic content, and therefore, whether strain 148 was probably more closely related to clinical strains than to those isolated from the ocean water column or other wild environments.
The genomic similarity of strain M18 to clinical P. aeruginosa strains (http://www.pseudomonas.com) suggests that clinical isolates are not a subpopulation of P. aeruginosa, which is the main question of our research, but it is still possible that strain M18 is a particular case that does not reflect the general situation of environmental P. aeruginosa strains.
In this work we show that the three environmental strains studied, and the dolphin-associated strain, are all pathogenic and phenotypically diverse with respect to production of virulence factors and motility. However, this phenotypic diversity is not reflected in genomic variability, since their genomes have a very high degree of similarity. Furthermore, detailed comparison of all the orthologous proteins encoded in the genomes of the seventeen P. aeruginosa strains studied in this work (Table 1), shows that the genomes of such a diverse group of P. aeruginosa strains (twelve isolates from human, two water inhabitants, two plant-associated strains and one strain isolated from a dolphin) are surprisingly similar, thus having a large proportion of their genome conserved between them (core genome). We also found that the genomes of the four isolates described in this study (Table 1) harbor sequences that have been reported to be part of genomic islands (GIs) that have been related to pathogenic phenotypes [18, 27–29], suggesting that these strains share genetic information with clinical isolates. However the four studied strains do not show antibiotic resistance, suggesting that the high incidence of multi-resistant strains among clinical isolates is due to selection pressures for these traits caused by the use of antibiotics for the treatment of P. aeruginosa infections.
The database of all the protein families, their abundance and the known molecular functions encoded by the genomes of the seventeen P. aeruginosa strains analyzed in this work, is provided (Additional file 1: Table S1 http://figshare.com/articles/Pan_genome_Pseudomonas_aeruginosa_Supplementary_Table_S1/760583). This constitutes a contribution for researchers working with this opportunistic bacterium.
Results and discussion
Phenotypic characterization of three environmental and a dolphin-associated P. aeruginosastrains
We tested the resistance-pattern of strains ID4365, M10, IGB83 and 148 to 20 antibiotics (see Methods section) and found that all four were sensitive to all antibiotics tested. The multisensitive pattern of these strains is completely different to the patterns shown by a collection of 100 P. aeruginosa clinical isolates, which are predominantly multi-resistant . The result obtained is expected since neither the environmental strains nor the dolphin-associated strain have been subjected to selection by the presence of antibiotics used in the treatment of P. aeruginosa infections.
Determination of strain virulence in the mouse model*
% lethality (108)
% lethality (109)
% lethality (1010)
Phenotypic characterization of P. aeruginosa strains analyzed compared to the type strain PAO1 a
141.1 ± 15.7
178.2 ± 16.8
2.34 ± 0.46
5.29 ± 0.23
0.23 ± 0.06
0.34 ± 0.04
170.2 ± 9.9
255 ± 11.1
197.2 ± 6
24.5 ± 1.5
234.5 ± 17.8
152.8 ± 5.3
1.6 ± 0.6
0.36 ± 0.1
0.41 ± 0.04
0.53 ± 0.11
297.6 ± 5
257.6 ± 31.9
0.32 ± 0.02
0.43 ± 0.03
87.9 ± 10.4
158.6 ± 5.2
14.05 ± 0.44
8.95 ± 0.28
In contrast to what has been reported for the clinical strains PAO1 and PA14 strains (Grosso-Becerra MV, Croda-García G, Merino E, Servín-González L, Mojica-Espinosa R, Soberón-Chávez G: Regulation of Pseudomonas aeruginosa virulence factors by two novel RNA thermometers, submitted, ), we found that the Indian Ocean isolate, strain ID4365 produces high levels of pyocyanin at 37°C  and even higher levels at 30°C (Table 3). This pattern of thermoregulation is similar to that reported for strain M18 , but the marine strain ID4365 lacks the sequence present in the 3’ region of the operon phzA1B1C1D1E1F1G1 that has been implicated in thermoregulation in strain M18 .
Strain M10, the water-spring isolate, produces similar levels of rhamnolipids, pyocyanin, elastase (Table 3), 3O-C12-HSL autoinducer (Figure 1) and LasR (Figure 2) to those produced by PAO1, but its production of rhamnolipids and pyocyanin is higher at 30°C than at 37°C (Table 3). Another particular feature of this strain is the low concentration of RhlR and C4-HSL produced at both temperatures tested (Figures 1 and 2).
Pseudomonas aeruginosa has been found to establish pathogenic interactions with different hosts including plants . Strain IGB83 can be considered not only an environmental strain, but also a potential plant pathogen, since it was isolated from a rotten fruit. We found that this strain produces very low levels of pyocyanin at either temperature tested (Table 3), and that rhamnolipids production is not thermoregulated, as it is in strains PAO1, ID4365 and 148 (Table 3).
Strain 148, isolated from a dolphin stomach shows the same pattern of thermoregulation of rhamnolipids, C4-HSL and RhlR production as PAO1 type strain (Table 3; Figures 1 and 2), but its pyocyanin production is higher at 30°C than at 37°C (Table 3). This strain completely lacks elastase activity and production of 3O-C12-HSL. We were unable to detect the presence of LasR (Table 3; Figures 1 and 2).
Contrary to what has been found for two P. aeruginosa clinical isolates, PAO1 and PA14 (Grosso-Becerra MV, Croda-García G, Merino E, Servín-González L, Mojica-Espinosa R, Soberón-Chávez G: Regulation of Pseudomonas aeruginosa virulence factors by two novel RNA thermometers, submitted, ), three of the strains studied in this work (ID4365, M10 and 148), synthesize higher levels of pyocyanin at an environmental temperature than at the human body temperature (Table 3).
Strain 148 is unable to swim or swarm (Table 3; Additional file 2: Figure S1). The non-motility of the two strains isolated from marine-related habitats (ID4375 and 148) suggests that this might be a common phenotype among P. aeruginosa living in the sea. The non-motile phenotype of strain 148 is due to multiple mutations in fliC and fliD genes (Additional file 3: Figure S2). However the genomic basis of the strain ID4365 inability to swim or swarm was not detected, since there are no-apparent mutations in the different genes involved in swimming or swarming that were analyzed (Adittional file 3: Figure S2).
Analysis of QS-dependent genes
Detailed analysis of the sequence of QS-dependent genes such as rhlA, rhlB, rhlR, rhlI, lecA, lasB, and lasI showed that they are identical in the four P. aeruginosa strains studied (ID4365, M10, IGB83, and 148), and that is also the case for lasR in the three environmental isolates. Genome assembly of these four strains does not enable the full-length sequence annotation of the duplicated phzA1B1C1D1E1F1G1 and phzA2B2C2D2E2F2G2 operons since they are nearly identical and all of the sequences collapse into a single contig. This is a well-known problem for duplicated sequences in genome assembly. In order to determine whether both operons are present in the four strains, we amplified the 5’ regions of both operons using PAO1 DNA as a control. We found that both phz operons are present in strains 148, ID4365, IGB83 and M10 (Additional file 2: Figure S3).
The comparisons, using RNA-Seq, of transcriptomes from three environmental P. aeruginosa strains showed that the expression of quorum-sensing related genes is widely different, but that the master quorum-sensing regulators in all of them are well conserved . Recently it was reported that two P. aeruginosa isolates from cystic fibrosis patients showed a wide phenotypic variability, even though both strains have a highly conserved genome sequence . These observations are completely in accordance with the results presented in this work.
We have been unable to detect the genomic characteristics responsible for the low level of expression of the Las regulon in strain ID4365, nor the very high levels of pyocyanin that it produces, particularly at 30°C (Table 3; Figures 1 and 2). The genomic basis responsible for particular phenotypic characteristics of the four studied strains, which are different to those of clinical isolates, are also not clear.
The high genomic conservation of P. aeruginosa strains, however, is not reflected in a homogeneous pattern of QS-regulated virulence-traits expression. In this work we detected one strain that does not produce pyocyanin (IGB83, Table 3), two other strains that do not produce elastase (ID4365 and 148), and one strain (148) that does not produce the autoinducer 3O-C12-HSL, nor the LasR protein (Table 3; Figures 1 and 2). We detected that, contrary to the case of PAO1 and PA14 strains [31, 32], higher levels of pyocyanin are produced at 30°C rather than at 37°C by the three studied strains that synthesize this phenazine, even though RhlR and the autoinducer C4-HSL tend to be produced at higher levels at 37°C by all the four strains studied. Despite these differences, all strains are virulent on the mice model, but show a slightly lower virulence than PAO1 strain (Table 2). It is our interest to determine the molecular mechanisms that underlie the different patterns of QS-regulated virulence trait expression on strains 148, ID4365, IGB83 and M10.
Presence of genomic islands in strains 148, ID4365, IGB83 and M10
All sequenced P. aeruginosa strains sequenced to date show DNA sequences inserted in the core genome that are strain-specific or shared only by a fraction of isolates, which appeared to be transferred horizontally between strains [18, 27, 28]. Among these insertions, called genomic islands (GIs), there are some that are implicated in pathogenic interactions, and are thus assumed to be a characteristic of clinical isolates. The resistance to multiple antibiotics is also assumed to be only present in clinical isolates.
It was recently reported that a collection of P. aeruginosa clinical isolates has a high incidence of GIs . The best characterized P. aeruginosa GIs related to pathogenicity are PAGI-1 to PAGI-4, PAPI-1, PAPI-2 and pKLC102 , but it has not been determined whether these elements can also be found in environmental isolates. We therefore searched for the presence of GIs in the chromosome of strains 148, ID4365, IGB83 and M10, and found that all four contain at least one ORF that has been reported to be part of GIs (Figure 3). The marine strain ID4365 has the highest number of GIs-related genes (Figure 3). This result reinforces the observation that environmental and clinical P. aeruginosa isolates do not constitute different subpopulations, and that GIs sequences are transferred horizontally between strains present in environments that are not close to humans and clinical strains.
P. aeruginosacore and pan-genome analysis
Global analysis of seventeen P. aeruginosa strains (twelve isolated from humans, two from water, two from plants and one from a dolphin, Table 1) was performed by best reciprocal blast hits (BBH) of all the proteins encoded in their genomes. Using this analysis we found that the P. aeruginosa core genome consists of 4,455 orthologous encoding genes present in all the seventeen strains analyzed. Considering that some of the genes in the core genome are part of paralogous gene families we determined a total of 74,731 core genes in all the seventeen P. aeruginosa strains analyzed.
This conserved genomic structure of P. aeruginosa is completely different to that reported for other bacterial species. A unique feature is that the core genome of this species constitutes more than 80% of its genetic repertoire and that its pan-genome is considerably smaller than the ones previously defined for other bacteria. For example the analysis of 61 sequenced Escherichia coli genomes showed that only 993 genes were present in all the genomes analyzed (core genome), and that this represented only 6% of the total number of gene families identified .
The analysis of the core genome and of the pan-genome by GSS provides the opportunity to get a broader view of the shared features across all the strains and to highlight the differences that can help us to understand particular niche adaptation strategies. The complete list of predicted proteins of the P. aeruginosa pan-genome and their distribution among strains is a valuable tool for the identification of genes that might code for traits that are important for niche adaptation (Additional file 1: Table S1). Particularly, this is a valuable resource to target human-related traits for diagnosis and direct therapy efforts. In this respect, the group of human-related P. aeruginosa strains is the one with the largest number of unique genes (Figure 4). The high number of unique genes in this group might be due to the over-representation of these strains among our sample (they constitute more than two thirds of the strains analyzed). Among the isolates from humans, there are several examples of genes that might be related to virulence, like those coding for a transcription regulator of the AraC family with high similarity to CdhR transcription regulators , fimbria related proteins, type IV secretion system, and adhesion related proteins, for example.
In conclusion, the analysis of the genome sequence of three environmental and one dolphin-associated P. aeruginosa strains shows that their genomic content is extremely well conserved among them and with respect to previously sequenced strains (Figure 5). The seventeen strains analyzed (Table 1), which have representatives of isolates from diverse environments (human-, water-, plant- and dolphin-associated strains), grouped as a single clade with very low genetic diversity. The inset in Figure 5 shows that there is no such thing as a subgroup of the species that clusters apart based on the isolation environment. Strains isolated from humans distribute evenly and across the P. aeruginosa cluster and they can be more closely related to water, dolphin or plant isolated strains than to each other. Furthermore we found that the core genome defined as the number of proteins that are encoded in all the seventeen genomes sequences analyzed in detail represents a high proportion (more than 80%) of the total number of proteins encoded in them. These results rule out the possibility that clinical strains constitute a subpopulation of P. aeruginosa and that the genetic variability of this bacterial species is presented in environmental isolates. Our results reinforce previously reported results that suggest that genetic variability of P. aeruginosa is extremely low with high levels of recombination within the species [2, 3, 10, 21, 39].
The extremely high sequence similarity among different P. aeruginosa isolates is difficult to explain. One possibility to account for this high degree of genetic conservation between P. aeruginosa strains isolated from such diverse habitats would be that these bacteria have an inefficient mechanism of DNA exchange, but this does not seem to be the case as judged by the prevalence of GIs among clinical  and environmental isolates (Figure 3), and the population studies suggesting high levels of homologous recombination . Other possibilities for explaining the low genetic variability of P. aeruginosa are that the population size of this bacterium is very low, or that a mechanism of gene conversion exists among different strains. These possibilities remain to be evaluated.
The high degree of sequence conservation of P. aeruginosa strains makes it a difficult task to device strategies for the treatment of infections caused by these bacteria that are based on vaccination or hygiene measures, and strengthens the need for designing strategies that inhibit the expression of virulence associated traits to contend with P. aeruginosa infections. Novel mechanisms of virulence inhibition, which take into account all the phenotypic diversity of this ubiquitous opportunistic pathogen, remain to be discovered.
Strains, culture conditions and microbiological procedures
P. aeruginosa strains (Table 1) were routinely propagated at 37°C on Luria-Bertani (LB) agar or LB broth  and all liquid cultures were grown with shaking (225 rpm). For experiments assessing the effect of temperature PPGAS media (phosphate-limited-peptone-glucose-ammonium) , was inoculated at a starting OD600 of 0.1, with overnight cultures in LB supplemented with the appropriate antibiotics and incubated at 30°C or at 37°C.
SDS gel electrophoresis and western blot analysis
For the preparation of crude cell extracts, PPGAS cultures were grown at 30°C and 37°C to an OD600 of 1.5. Cells were harvested and re-suspended in PA buffer (10 mM sodium phosphate buffer, 30 mM NaCl, 0.25% Tween-20, 10 mM EDTA, 10 mM β-mercaptoethanol, pH 7.5) prior to cell disruption through sonication. Cellular debris was removed by centrifugation (13,000 g for 15 min at 4°C) and the supernatant (containing soluble protein fraction) was mixed 1:1 v/v with Laemmli loading buffer . Total protein concentration was determined by protein assay kit (Bio-Rad) with bovine serum albumin as standard. Equal amounts of proteins were separated by 12% sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and gels were electro-transferred onto a Hybond-C Extra nitrocellulose membrane (Amersham Biosciences). The membrane was blocked by 5% nonfat milk and incubated with a 1:1000 dilution of rabbit polyclonal antibody of either anti-RhlR or anti-LasR. Goat anti-rabbit immunoglobulin G (Santa Cruz Biotechnology) secondary antibody conjugated to horseradish peroxidase was used at a 1:10,000 of dilution. Detection was performed with a chemiluminescence-based system SuperSignal West Femto (Pierce) followed by exposure to an X-ray film (Amersham Biosciences) for autoradiography. PAO1 derived rhlR and lasR mutants were used as negative controls and as positive controls we use Escherichia coli strains harboring plasmids expressing rhlR (pECP61.5) and lasR (pECP64) .
Virulence factors production
Pyocyanin was extracted from culture supernatants and measured as previously described . The pyocyanin assay is based on the absorbance of pyocyanin at 520 nm in acidic solution. A 5 ml sample of culture supernatant was extracted with 3 ml of chloroform and then re-extracted into 1 ml of 0.2 N HCl to give a pink to deep red solution. The absorbance of this solution was measured at 520 nm. Concentrations, expressed as micrograms of pyocyanin produced per milliliter of culture supernatant, were determined by multiplying the optical density at 520 nm (OD520) by 17.072.
The concentration of rhamnolipids in the sample was estimated by the Orcinol method . A 333 μl portion of the filtered supernatant was extracted twice with 1 ml of diethyl ether. The diethyl ether was evaporated to dryness and dissolved in 1 ml of deionized water. To 100 μl of each sample, 900 μl of a solution containing 0.19% orcinol (in 53% sulfuric acid) was added. The samples were heated at 80°C in a water-bath for 30 min and cooled for 15 min at room temperature and the A421 was measured. Concentrations of rhamnolipids were determined by comparing the data with those obtained with L-rhamnose standards between 0 and 50 μg/ml.
The elastolytic activity of LasB elastase was determinate using elastin Congo red as substrate and the procedure was modified from that described previously . Briefly, the cells were grown in LB broth at 30° and 37°C, respectively for 24 h. Samples of the filter-sterilized supernatants were diluted 1:10 with LB and 1 mL was added 10 mg of Elastin Congo Red (Sigma) in glass tubes. The mixture was incubated at 37°C for 16 h with constant rotation (225 rpm), insoluble substrate was pelleted with centrifugation (1300 g for 10 min at 4°C) and absorbance of the supernatant was measured at 495 nm with a spectrophotometer using as blank elastin Congo red sample incubated with medium alone. The experiment was performed three times in triplicate with supernatants from three individual growth experiments. The values from the triplicate experiments were averaged and used as one value to represent each of the three experiments.
Acyl-homoserine lactone (AHL) extraction and analytical thin-layer chromatography (TLC)
To evaluate profiles of AHLs, cells were grown in PPGAS media at 30°C and 37°C, respectively for 24 h. A 10-ml sample of culture supernatant was extracted twice with equal volumes of acidified ethyl acetate and then dried in a fume hood. The residues of extraction were then dissolved in 100 μl ethyl acetate and 5 μl were analyzed by thin layer chromatography (TLC). Analytical TLC was performed on reverse phase aluminium-backed RP18 F254S TLC plates (20 cm X 20 cm; Merck) . Chromatograms were developed with methanol: water (60:40, v/v), then air-dried in a fume hood. The TLC plate was then overlaid with a thin film of agar seeded with the AHL reporter strain C. violaceum CV026 that produces the purple colour violacein in response to AHLs with N-acyl side chains between 4 and 8 carbons in length . After incubation of the plate at 30°C for 24 h, AHLs were located as purple spots on a white background. Alternatively, TLC plates were overlaid with a thin layer of agar seeded with a culture of E. coli biosensor strain containing lux-based bioluminescence AHL reporter plasmid (pSB1075). This reporter contains the P. aeruginosa lasR gene and lasI promoter fused to luxCDABE from P. luminescens and detects 3-oxo-substituted AHL derivatives with acyl chain length from 4 to 12 carbons . Light emission was detected using a system for chemiluminescence detection Gel Logic 112 (Kodak). All the experiments were performed at least twice.
Determination of antibiotic resistance
To assess the susceptibility profiles to 20 antimicrobial agents of the strains of Pseudomonas aeruginosa the agar dilution method was used according to the guidelines established by the National Committee for Clinical Laboratory Standards (NCCLS), as previously described . ATCC 27853 Pseudomonas aeruginosa, ATCC 25922 Escherichia coli, ATCC 35218 Escherichia coli, ATCC 29213 Staphylococcus aureus, and ATCC 29212 Enterococcus faecalis were used as controls in the susceptibility tests. All the strains were grown in Muller Hinton agar and harvested in sterile saline solution to achieve a turbidity equivalent to that of a No. 0.5 McFarland opacity standard. The antimicrobial agents tested against P. aeruginosa were: carbenicillin (16-64 μg/mL), ticarcillin (8-32 μg/mL), piperacillin (1-8 μg/mL) ticarcillin/clavulanic acid (8/2-32/2 μg/mL), piperacillin/tazobactam (1/4-8/4 μg/mL) ceftazidime (1-4 μg/mL) ceftriaxone (8-64 μg/mL) cefotaxime (8-32 μg/mL) cefepime (1-8 μg/mLl) imipenem (1-4 μg/mL) meropenem (0.25-1 μg/mlL) aztreonam (2-8 μg/mL) amikacin (1-4 μg/mL) gentamicin (0.5-2 μg/mL) tobramycin (0.25-1 μg/mL) polymyxin b (0.25-2 μg/mL) ciprofloxacin (0.25-1 μg/mL) norfloxacin (1-4 μg/mL) and levofloxacin (0.5-4 μg/mL). Agar dilution was performed using two-fold increments (across a range of 0.125 to 512 μg/mL) of each antimicrobial agent incorporated into Muller-Hinton agar. The concentration range for susceptibility and resistance are indicated in parenthesis with a MIC value lower than the cut-off to indicate susceptibility and two-fold dilutions above the cut-off to determine resistance. The criterion for intermediate susceptibility was based on isolates growing within one-fold dilution higher than the MIC value.
All mouse studies were conducted in accord with Guide for the Care and Use of Laboratory Animals (Committee for the Update of the Guide for the Care and Use of Laboratory Animals and Institute for Laboratory Animal Research, Washington D. C., 2011) and the Comité para el Cuidado y Uso de Animales de Laboratorio (CCUAL) and were approved by the Ethics committees of Instituto de Investigaciones Biomédicas - UNAM (approval No. ID201 09/02-2010).
Bacteria to be injected into mice were cultured as described previously and quantitated by plate serial dilution method. For virulence studies, groups of five female BALB/cAnNHsd mice 6–8 weeks old (Harlan, Indianapolis, IN, United States) were inoculated intraperitoneally at day 0 with 100 μl of bacterial suspension at 1X107, 1X108, 1X109 and 1X1010 colony forming units (CFU) of the each P. aeruginosa strain, respectively. Bacterial suspensions were prepared in sterile water for injection. Intraperitoneal injection without bacterial cells was used as a control group. Before and after experimental use, the animals were maintained in cages at 5 mice per cage with food and water available ad libitum. Mice from all treatment groups were monitored at 24 h after infection and died mice were incinerated.
Genomes assembly and annotation
The genomes of isolates 148, ID4365, IGB83, and M10 were sequenced using the Illumina Genome Analyzer II. The libraries were prepared using a pair-ended protocol with an insertion average length of 500 bp. Quality control analysis and trimmi  and Celera assemblers . The read average length goes to 35.99 bp. The estimated coverage goes from 39.49 to 73.48X fold. The draft genomes assemblies were done first by quality control check by using FastX (http://hannonlab.cshl.edu/fastx_toolkit), and then inputted into Velvet genome assembler . Multiple assemblies were conducted to calibrate the optimum kmers parameters. We chose a global optimum of 23 kmers for de novo pair-ended assembly. The genomic alignments against the P. aeruginosa reference genomes were done with the MUMmeralignment suite  using the maxmatch option for nucleotide alignments. The MUMmer alignments were used to aid in the further assembly and scaffolding of the sequenced genomes, additional processing of the assembly was done with the Bambus program of the AMOS suite . A second independent assembly and processing was done by the CLC Genomics Workbench version 5.5.1, and scaffolding and gap closure using SSPACE  and GapFiller version 1.10 . Automated gene calling and annotation was conducted using the RAST Server .
The Whole Genome Shotgun (WGS) projects for P. aerugionsa isolates M10, IGB83, 148, ID4365 have been deposited at DDBJ/EMBL/Genbank under the project accessions ATAG00000000, ATAH00000000, ATAI00000000, ATAJ00000000, respectively. The versions described in this paper are the first versions.
Genome Similarity Score
Genome Similarity Score (GSS)  was built with seven of the public released Pseudomonas strains. This is done by getting all the pairwise orthologs by means of reciprocal best blast hits (BBH) and getting individual bit-scores, comparing the sum of the comparison bit-scores and normalizing it against the self GSS bit-scores sum. This score has been successfully used before when comparing the overall genomic similarity of related species [59, 60]. The range of GSS goes from 1 (identical strains) to 0 (non-related strains). The strength of this score relies on having genetic distances by means of comparing all the pairwise shared orthologs. The genetic distances are then sorted into a distance matrix and then plotted as a dendrogram. We also calculated the GSS distance with two out-groups strains of the genus Xanthomonas. To plot the GSS as dendrograms we got the inverse (1-GSS) distance value and then plot it using Mega (Neighbor-Joining) , and Figtree v 1.4.0 (http://tree.bio.ed.ac.uk/software/figtree/). The accessions for the genomes used in the GSS calculations are: NC_002516.2, NC_002947.3, NC_004129.6, NC_004578.1, NC_004632.1, NC_004633.1, NC_005244.2, NC_005773.3, NC_007005.1, NC_007274.1, NC_007275.1, NC_007492.2, NC_008027.1, NC_008463.1, NC_009434.1, NC_009439.1, NC_009512.1, NC_009656.1, NC_010322.1, NC_010501.1, NC_011770.1, NC_012660.1, NC_015379.1, NC_015410.1, NC_015556.1, NC_015733.1, NC_015740.1, NC_016830.1, NC_017530.1, NC_017532.1, NC_017548.1, NC_017549.1, NC_017911.1, NC_017986.1, NC_018028.1, NC_018080.1, NC_018177.1, NC_018220.1, NC_018746.1, NC_019670.1, NC_019905.1, NC_019906.1, NC_019936.1, NC_019937.1, NC_019938.1, NC_019939.1, NC_020209.1, NC_020829.1, NC_020912.1, and the out-groups NC_003902, NC_20800.1. The P. aeruginosa genomes 39016, 2192, B136, RP73, C3719, and PACS2 where downloaded in its last version from http://www.pseudomonas.org.
In order to assess the presence, absence and variation of certain GIs in the environmental strains’ genomes, we performed tblastn comparisons. As database, we used a Multi-FASTA file with the nucleotide sequences of different genomic islands from the PAGI and PAPI families, as well as from the integrative plasmid pKLC102; as query, we used the protein sequences from all the annotated genes in strains 148, ID4365, IGB83 and M10. A ring representation was then created with the BLAST Ring Image Generator . The accessions for the GIs used in this analysis are: AF241171, AF440523, AF440524, AY258138, EF611301, EF611302, EF611303, EF611304, EF611305, EF611306, EF611307, NC_019202, AY273869, and AY273870.
The core genome sensu stricto was obtained by means of BBH, as described before , of all the available P. aeruginosa complete genomes available up-to date (May 2013). We define the core genome as the orthologs that are present in each and every one of the analyzed genomes. This is done by parsing the BBH results into tables and sorting them out by the presence/absence of a BBH match. The complete list of the isolates used in this study is shown in Table 1. The core genome sensu latu was done by protein clustering (80% identity) and present in at least 11 representatives in each cluster.
The pan-genome was described by means of comparing and clustering all the predicted proteins of the eleven analyzed strains (Table 1). The clustering was done using CD-HIT version 4.6 , using a 0.80 identity threshold, with longest-sequence-first list removal algorithm. Results were parsed into a table (Additional file 1: Table S1) where a representative sequence was selected to represent every single protein family, each P. aeruginosa isolate has a column indicating the protein family representatives copy numbers within each of the analyzed genomes. Next, some basic calculations are shown displaying the total number of the members of each protein family. Each orthologous identified by the core genome analysis is shown in the column titled “Core Orthologous”. A second way to identify the core genome would be to make protein clustering and a cut-off of abundance, this is shown in the column named “Core.Prot.Fam”, using a threshold of 11 for the minimum abundance. The genomes where clustered based on its environmental isolation types defining clusters of Water (ID4365, M10), Plant (IGB83, M18), Dolphin (148), and human-associated isolates (NCGM2, LESB58, PAO1, PA14, DK2, PA7). Additional file 1: Table S1 shows the count of the predicted proteins that are exclusive to one kind of these environmental types (Water.only, Plant.only, Dolphin.only, Human.only). Finally the predicted protein sequence is available to perform further analysis, at the last field of the table.
All the statistical analyses were done on R (ver 2.14.1) . Using the gplots , and vennerable (https://r-forge.r-project.org/projects/vennerable/) packages.
Availability of supporting data
The supporting data of our paper is now available. The reference is: Pan-genome Pseudomonas aeruginosa Additional file 1: Table S1. (http://figshare.com/articles/Pan_genome_Pseudomonas_aeruginosa_Supplementary_Table_S1/760583).
We are grateful to Dra. Jeiry Toribio and Dr. Rigoberto Hernández for providing strains M10 and 148, respectively. We acknowledge the technical assistance in a polyclonal antibodies production to Alejandro Olvera Rodríguez and Hector Cardoso Torres. We greatly acknowledge Valeria Souza for critical discussion of this work’s fundament. This work was supported in part by grant PAPIIT IN202510 (DGAPA-UNAM) and CONACYT 100343. Victoria Grosso-Becerra and Abigail González-Valdez are PhD students of the Programa de Ciencias Bioquímicas, UNAM.
- Hardalo HC, Edberg SC: Pseudomonas aeruginosa: Assessment of risk from drinking water. Crit Rev Microbiol. 1997, 23: 47-75. 10.3109/10408419709115130.PubMedView ArticleGoogle Scholar
- Pirnay JP, Bilocq F, Pot B, Cornelis P, Zizi M, Van Eldere J, Deschaght P, Vaneechoutte M, Jennes S, Pitt T, De Vos D: Pseudomonas aeruginosa population structure revisited. PLoS One. 2009, 4: e7740-10.1371/journal.pone.0007740.PubMed CentralPubMedView ArticleGoogle Scholar
- Silby MW, Winstanley C, Godfrey SAC, Levy SB, Jackson RW: Pseudomonas genomes: diverse and adaptable. FEMS Microbiol Rev. 2001, 35: 652-680.View ArticleGoogle Scholar
- Williams P, Cámara M: Quorum sensing and environmental adaptation in Pseudomonas aeruginosa: a tale of regulatory networks and multifunctional signal molecules. Curr Opin Microbiol. 2009, 12: 182-191. 10.1016/j.mib.2009.01.005.PubMedView ArticleGoogle Scholar
- Costerton JW: 1980. Pseudomonas aeruginosa. The organism, diseases it causes and their treatment. Edited by: Sabath CD. 1980, Switzerland: Hans Huber Publishers, 15-24.Google Scholar
- Govan JR, Deretic V: Microbial pathogenesis in cystic fibrosis: mucoid Pseudomonas aeruginosa and Burkhordelia cepacia. Microbiol Rev. 1996, 60: 539-574.PubMed CentralPubMedGoogle Scholar
- Arruda EA, Marinho IS, Boulos M, Sinto SI, Caiaffa HH, Mendes CM, Oplustil CP, Sader H, Levy CE, Levin AS: Nosocomial infections caused by multiresistant Pseudomonas aeruginosa. Infect Control Hosp Epidemiol. 1999, 20: 620-623. 10.1086/501683.PubMedView ArticleGoogle Scholar
- Stover CK, Pham XQ, Erwin AL, Mizoguchi SD, Warrener P, Hickey MJ, Brinkman FS, Hufnagle WO, Kowalik DJ, Lagrou M, Garber RL, Goltry L, Tolentino E, Westbrock-Wadman S, Yuan Y, Brody LL, Coulter SN, Folger KR, Kas A, Larbig K, Lim R, Smith K, Spencer D, Wong GK, Wu Z, Paulsen IT, Reizer J, Saier MH, Hancock RE, Lory S, Olson MV: Complete genome sequence of Pseudomonas aeruginosaPAO1, an opportunistic pathogen. Nature. 2000, 406: 959-964. 10.1038/35023079.PubMedView ArticleGoogle Scholar
- Holloway BW: Genetic recombination in Pseudomonas aeruginosa. J Gen Microbiol. 1955, 13: 572-581. 10.1099/00221287-13-3-572.PubMedGoogle Scholar
- Mathee K, Narasimhan G, Valdes C, Qiu X, Matewish JM, Koehrsen M, Rokas A, Yandava CN, Engels R, Zeng E, Olavarietta R, Doud M, Smith RS, Montgomery P, White JR, Godfrey PA, Kodira C, Birren B, Galagan JE, Lory S: Dynamics of Pseudomonas aeruginosagenome evolution. Proc Natl Acad Sci USA. 2008, 105: 3100-3105. 10.1073/pnas.0711982105.PubMed CentralPubMedView ArticleGoogle Scholar
- Stewart RMK, Wiehlmann L, Ashelford KE, Preston SJ, Frimmersdorf E, Campbell BJ, Neal TJ, Hall N, Tuft S, Kaye SB, Winstanley C: Genetic characterization indicates that a specific subpopulation of Pseudomonas aeruginosais associated with keratitis infections. J Clin Microbiol. 2011, 49: 993-1003. 10.1128/JCM.02036-10.PubMed CentralPubMedView ArticleGoogle Scholar
- Rau MH, Marvig RL, Ehrlich GD, Molin S, Jelsbak L: Deletion and acquisition of genomic content during early stage adaptation of Pseudomonas aeruginosa to a human host environment. 2012, 14: 2200-2211.Google Scholar
- Winstanley C, Langille MG, Fothergill JL, Kukavica-Ibrulj I, Paradis-Bleau C, Sanschagrin F, Thomson NR, Winsor GL, Quail MA, Lennard N, Bignell A, Clarke L, Seeger K, Saunders D, Harris D, Parkhill J, Hancock RE, Brinkman FS, Levesque RC: Newly introduced genomic prophage islands are critical determinants of in vivo competitiveness in the Liverpool Epidemic Strain of Pseudomonas aeruginosa. Genome Res. 2009, 19: 12-23.PubMed CentralPubMedView ArticleGoogle Scholar
- Wu DQ, Ye J, Ou HY, Wei X, Huang X, He YW, Xu Y: Genomic analysis and temperature-dependent transcriptome profiles of the rhizosphere originating strain Pseudomonas aeruginosa M18. BMC Genomics. 2011, 12: 438-10.1186/1471-2164-12-438.PubMed CentralPubMedView ArticleGoogle Scholar
- Miyoshi-Akiyama T, Kuwahara T, Tada T, Kitao T, Kirikae T: Complete Genome Sequence of Highly Multidrug-Resistant Pseudomonas aeruginosa NCGM2.S1, a Representative Strain of a Cluster Endemic to Japan. J Bacteriol. 2011, 193: 7010-10.1128/JB.06312-11.PubMed CentralPubMedView ArticleGoogle Scholar
- Roy PH, Tetu SG, Larouche A, Elbourne L, Tremblay S, Ren Q, Dodson R, Harkins D, Shay R, Watkins K, Mahamoud Y, Paulsen IT: Complete genome sequence of the multiresistant taxonomic outlier Pseudomonas aeruginosaPA7. PLoS ONE. 2010, 5: e8842-10.1371/journal.pone.0008842.PubMed CentralPubMedView ArticleGoogle Scholar
- Lee DG, Urbach JM, Wu G, Liberati NT, Feinbaum RL, Miyata S, Diggins LT, He J, Saucier M, Deziel E, Friedman L, Li L, Grills G, Montgomery K, Kucherlapati R, Rahme LG, Ausubel FM: Genomic analysis reveals that Pseudomonas aeruginosa virulence is combinatorial. Genome Biol. 2006, 7: R90-10.1186/gb-2006-7-10-r90.PubMed CentralPubMedView ArticleGoogle Scholar
- Jeukens J, Boyle B, Bianconi I, Kukavica-Ibrulj I, Tümmler B, Bragonzi A, Levesque RC: Complete genome sequence of persistent Cystic Fibrosis isolate Pseudomonas aeruginosa strain RP73. Genome Announc. 2013, 1 (4): e00568-13. http://genomea.asm.org/content/1/4/e00568-13,PubMed CentralPubMedView ArticleGoogle Scholar
- Huang J, Xu Y, Zhang H, Li Y, Huang X, Ren B, Zhang X: Temperature-dependent expression of phzMand its regulatory genes lasIand ptsPin rhizosphereisolate Pseudomonassp. strain M18. Appl Environ Microbiol. 2009, 75: 6568-6580. 10.1128/AEM.01148-09.PubMed CentralPubMedView ArticleGoogle Scholar
- Lu J, Huang X, Zhang M, Li S, Jiang H, Xu Y: The distinct quorum sensing hierarchy of lasand rhlin Pseudomonassp. M18. Curr Microbiol. 2009, 59: 621-627. 10.1007/s00284-009-9483-y.PubMedView ArticleGoogle Scholar
- Wolfgang MC, Kulasekara BR, Liang X, Boyd D, Wu K, Yang Q, Miyada CG, Lory S: Conservation of genome content and virulence determinants among clinical and environmental isolates of Pseudomonas aeruginosa. Proc Natl Acad Sci U S A. 2003, 100: 8484-8489. 10.1073/pnas.0832438100.PubMed CentralPubMedView ArticleGoogle Scholar
- Welch RA, Burland V, Plunkett G, Redford P, Roesch P, Rasko D, Buckles EL, Liou SR, Boutin A, Hackett J, Stroud D, Mayhew GF, Rose DJ, Zhou S, Schwartz DC, Perna NT, Mobley HL, Donnenberg MS, Blattner FR: Extensive mosaic structure revealed by the complete genome sequence of uropathogenic Escherichia coli. 2002, 99: 17020-17024.Google Scholar
- Silby MW, Cerdeño-Tárraga AM, Vernikos GS, Giddens SR, Jackson RW, Preston GM, Zhang XX, Moon CD, Gehrig SM, Godfrey SA, Knight CG, Malone JG, Robinson Z, Spiers AJ, Harris S, Challis GL, Yaxley AM, Harris D, Seeger K, Murphy L, Rutter S, Squares R, Quail MA, Saunders E, Mavromatis K, Brettin TS, Bentley SD, Hothersall J, Stephens E, Thomas CM, et al: Genomic and genetic analyses of diversity and plant interactions of Pseudomonas fluorescens. Genome Biol. 2009, 10: R51-10.1186/gb-2009-10-5-r51.PubMed CentralPubMedView ArticleGoogle Scholar
- Manwar AV, Khandelwal SR, Chaudhari BL, Meyer JM, Chincholkar SB: Siderophore production by a marine Pseudomonas aeruginosaand its antagonistic action against phytopathogenic fungi. Appl Biochem Biotechnol. 2004, 118: 243-251. 10.1385/ABAB:118:1-3:243.PubMedView ArticleGoogle Scholar
- Souza V, Espinosa-Asuar L, Escalante AE, Eguiarte LE, Farmer J, Forney L, Lloret L, Rodríguez-Martínez JM, Soberón X, Dirzo R, Elser JJ: An endangered oasis of aquatic microbial biodiversity in the Chihuahuan desert. Proc Natl Acad Sci USA. 2006, 103: 6565-6570. 10.1073/pnas.0601434103.PubMed CentralPubMedView ArticleGoogle Scholar
- Palmeros B, Güereca L, Alagón A, Soberón-Chávez G: Biochemical characterization of the lipolytic activity of Pseudomonas aeruginosaIGB 83. Process Biochem. 1994, 29: 207-212. 10.1016/0032-9592(94)85005-4.View ArticleGoogle Scholar
- Liang X, Pham XQ, Olson M, Lory S: Identification of a genomic island present in the majority of pathogenic isolates of Pseudomonas aeruginosa. J Bacteriol. 2000, 183: 843-853.View ArticleGoogle Scholar
- Juhas M, van der Meer JR, Gaillard M, Harding RM, Hood DW, Crook DW: Genomic Island: tools of bacterial horizontal gene transfer and evolution. FEMS Microbiol Rev. 2008, 33: 376-393.PubMed CentralPubMedView ArticleGoogle Scholar
- Morales-Espinosa R, Soberón-Chávez G, Delgado-Sapién G, Sandner-Miranda L, Méndez JL, González-Valencia G, Cravioto A: Genetic and phenotypic characterization of a Pseudomonas aeruginosa populationwith high frequency of genomic islands. PLoS ONE. 2012, 7: e37459-10.1371/journal.pone.0037459.PubMed CentralPubMedView ArticleGoogle Scholar
- Kumar A, Munder A, Aravind R, Eapen SJ, Tümmler B, Raaijmakers JM: Friend or foe: genetic and functional characterization of plant endophytic Pseudomonas aeruginosa. Environ Microbiol. 2012, 15: 764-779.PubMedView ArticleGoogle Scholar
- Wurtzel O, Yoder-Himes DR, Han K, Dandekar AA, Edelheit S, Greenberg EP, Sorek R, Lory S: The single-nucleotide resolution transcriptome of Pseudomonas aeruginosagrown in body temperature. PLoS Pathog. 2012, 8: e1002945-10.1371/journal.ppat.1002945.PubMed CentralPubMedView ArticleGoogle Scholar
- Rahme LG, Le MWL, Wong SM, Tompkins RG, Calderwood SB, Ausubel F: Use of model plant hosts to identify Pseudomonas aeruginosavirulence factors. Proc Natl Acad Sci USA. 1997, 94: 13245-13250. 10.1073/pnas.94.24.13245.PubMed CentralPubMedView ArticleGoogle Scholar
- Ni N, Li M, Wang J, Wang B: Inhibitors and antagonist of bacterial quorum sensing. Med Res Rev. 2009, 29: 65-124. 10.1002/med.20145.PubMedView ArticleGoogle Scholar
- Chugani S, Kim BS, Pattarosokol S, Brittnacher MJ, Choi SH, Harwood CS, Geenberg EP: Strain-dependent diversity in the Pseudomonas aeruginosaquorum-sensing regulon. Proc. Natl. Acad. Sci. USA. 2012, 109: E2823-E2831. 10.1073/pnas.1214128109.PubMed CentralPubMedView ArticleGoogle Scholar
- Klockgether J, Miethke N, Kubesch P, Bohn Y-S, Brockhausen I, Cramer N, Eberl L, Greipel J, Herrmann C, Herrmann S, Horatzek S, Lingner M, Luciano L, Salunkhe L, Schomburg D, Wehsling M, Wiehlmann L, Davenport CF, Tümmler B: Intraclonal diversity of the Pseudomonas aeruginosa cystic fibrosis airway isolates TBCF10839 and TBCF121838: distinct signatures of transcriptome, proteome, metabolome, adherence and pathogenicity despite an almost identical genome sequence. Environ Microbiol. 2013, 15: 191-210. 10.1111/j.1462-2920.2012.02842.x.PubMedView ArticleGoogle Scholar
- Lukjancenko O, Wassenaar TM, Ussery DW: Comparison of 61 sequenced Escherichia coligenomes. Microb Ecol. 2010, 60: 708-720. 10.1007/s00248-010-9717-3.PubMed CentralPubMedView ArticleGoogle Scholar
- Moreno-Hagelsieb G, Janga SC: Operons and the effect of genome redundancy in deciphering functional relationships using phylogenetic profiles. Proteins. 2008, 70: 344-352.PubMedView ArticleGoogle Scholar
- Chawla A, Hirano T, Bainbridge BW, Demuth DR, Xie H, Lamont RJ: Community signalling between Streptococcus gordoniiand Porphyromonas gingivalisis controlled by the transcriptional regulator CdhR. Mol Microbiol. 2010, 78: 1510-1522. 10.1111/j.1365-2958.2010.07420.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Van Mansfeld R, Jongerden I, Bootsma M, Buiting A, Bonten M, Willems R: The population genetics of Pseudomonas aeruginosaisolates from different patient populations exhibits high-level host specificity. PloS ONE. 2010, 5: e13482-10.1371/journal.pone.0013482.PubMed CentralPubMedView ArticleGoogle Scholar
- Miller J: In Experiments in Molecular Genetics. 1972, New York: Cold Spring Harbor Laboratory, 352-355.Google Scholar
- Zhang Y, Miller RM: Enhancement of octadecane dispersion and biodegradation by a Pseudomonasrhamnolipid surfactant (biosurfactant). Appl Environ Microbiol. 1992, 58: 3276-3282.PubMed CentralPubMedGoogle Scholar
- Masduki A, Nakamura J, Ohga T, Umezaki R, Kato J, Ohtake H: Isolation and characterization of chemotaxis mutants and genes of Pseudomonas aeruginosa. J Bacteriol. 1995, 177: 948-952.PubMed CentralPubMedGoogle Scholar
- Tremblay J, Déziel E: Improving the reproducibility of Pseudomonas aeruginosaswarming motility assays. J Basic Microbiol. 2008, 48: 509-515. 10.1002/jobm.200800030.PubMedView ArticleGoogle Scholar
- Laemmli UK: Cleavage of structural proteins during the assembly of the head of bacteriophage T4. Nature- London. 1970, 277: 680-685.View ArticleGoogle Scholar
- Pearson JP, Pesci EC, Iglewski BH: Roles of Pseudomonas aeruginosa lasand rhlquorum sensing systems in control of elastase and rhamnolipid biosynthesis genes. J Bacteriol. 1997, 179: 3127-31321.PubMed CentralPubMedGoogle Scholar
- Essar DW, Eberly L, Crawford IP: Evolutionary differences in chromosomal locations of four early genes of tryptophan pathway in fluorescent Pseudomonas: DNA sequences and characterization of Pseudomonas putida trpEand trpGDC. J Bacteriol. 1990, 172: 867-883.PubMed CentralPubMedGoogle Scholar
- Chandrasekaran EV, Bemiller JN: Constituent analyses of glycosaminoglycans. Methods Carbohydr Chem. 1980, 8: 89-96.Google Scholar
- Beatson SA, Whitchurch CB, Sargent JL, Levesque RC, Mattick JS: Differential Regulation of Twitching Motility and Elastase Production by Vfr in Pseudomonas aeruginosa. J Bacteriol. 2002, 184: 3605-3613. 10.1128/JB.184.13.3605-3613.2002.PubMed CentralPubMedView ArticleGoogle Scholar
- Middleton B, Rodgers HC, Cámara M, Knox AJ, Williams P, Hardman A: Direct detection of N-acylhomoserine lactones in cystic fibrosis sputum. FEMS Microbiol Lett. 2002, 207: 1-7. 10.1111/j.1574-6968.2002.tb11019.x.PubMedView ArticleGoogle Scholar
- Shaw PD, Ping G, Daly SL, Cha C, Cronan JE, Rinehart KL, Farrand SK: Detecting and characterizing N-acyl-homoserine lactone signal molecules by thin-layer chromatography. Proc Natl Acad Sci USA. 1997, 94: 6036-6041. 10.1073/pnas.94.12.6036.PubMed CentralPubMedView ArticleGoogle Scholar
- Winson MK, Swift S, Fish L, Throup JP, JØrgensen F, Chhabra SR, Bycroft EW, Williams P, Stewart G: Construction and analysis of luxCDABE-based plasmid sensors for investigating N-acyl homoserine lactone-mediated quorum sensing. FEMS Microbiol Lett. 1998, 163: 185-192. 10.1111/j.1574-6968.1998.tb13044.x.PubMedView ArticleGoogle Scholar
- Zerbino DR, Birney E: Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Research. 2004, 18: 821-829.View ArticleGoogle Scholar
- Miller JR, Delcher AL, Koren S, Venter E, Walenz BP, Brownley A, Johnson J, Li K, Mobarry C, Sutton G: Aggressive assembly of pyrosequencing reads with mates. Bioinformatics. 2008, 24: 2818-2824. 10.1093/bioinformatics/btn548.PubMed CentralPubMedView ArticleGoogle Scholar
- Kurtz S, Phillippy A, Delcher AL, Smoot M, Shumway M, Antonescu C, Salzberg SL: Versatile and open software for comparing large genomes. Genome Biol. 2004, 5: R12-10.1186/gb-2004-5-2-r12.PubMed CentralPubMedView ArticleGoogle Scholar
- Sommer DD, Delcher AL, Salzberg SL, Pop M: Minimus: a fast, lightweight genome assembler. BMC Bioinformatics. 2007, 8: 64-10.1186/1471-2105-8-64.PubMed CentralPubMedView ArticleGoogle Scholar
- Boetzer M, Henkel CV, Jansen HJ, Butler D, Pirovano W: Scaffolding pre-assembled contigs using SSPACE. Bioinformatics. 2011, 27: 578-579. 10.1093/bioinformatics/btq683.PubMedView ArticleGoogle Scholar
- Boetzer M, Pirovano W: Toward almost closed genomes with GapFiller. Genome Biol. 2012, 13: R56-10.1186/gb-2012-13-6-r56.PubMed CentralPubMedView ArticleGoogle Scholar
- Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, Formsma K, Gerdes S, Glass EM, Kubal M, Meyer F, Olsen GJ, Olson R, Osterman AL, Overbeek RA, McNeil LK, Paarmann D, Paczian T, Parrello B, Pusch GD, Reich C, Stevens R, Vassieva O, Vonstein V, Wilke A, Zagnitko O: The RAST Server: rapid annotations using subsystems technology. BMC Genomics. 2008, 9: 75-10.1186/1471-2164-9-75.PubMed CentralPubMedView ArticleGoogle Scholar
- Alcaraz LD, Moreno-Hagelsieb G, Eguiarte LE, Souza V, Herrera-Estrella L, Olmedo G: Understanding the evolutionary relationships and major traits of Bacillus through comparative genomics. BMC Genomics. 2010, 11: 332-10.1186/1471-2164-11-332.PubMed CentralPubMedView ArticleGoogle Scholar
- Kubicek CP, Herrera-Estrella A, Seidl-Seiboth V, Martinez DA, Druzhinina IS, Thon M, Zeilinger S, Casas-Flores S, Horwitz BA, Mukherjee PK, Mukherjee M, Kredics L, Alcaraz LD, Aerts A, Antal Z, Atanasova L, Cervantes-Badillo MG, Challacombe J, Chertkov O, McCluskey K, Coulpier F, Deshpande N, von Döhren H, Ebbole DJ, Esquivel-Naranjo EU, Fekete E, Flipphi M, Glaser F, Gómez-Rodríguez EY, Gruber S, et al: Comparative genome sequence analysis underscores mycoparasitism as the ancestral life style of Trichoderma. Genome Biol. 2011, 12: R40-10.1186/gb-2011-12-4-r40.PubMed CentralPubMedView ArticleGoogle Scholar
- Kumar S, Tamura K, Nei M: MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Brief Bioinform. 2004, 5: 150-163. 10.1093/bib/5.2.150.PubMedView ArticleGoogle Scholar
- Alikhan NF, Petty NK, Zakour NLB, Beatson SA: BLAST Ring Image Generator (BRIG): simple prokaryote genome comparisons. BMC Genomics. 2011, 12: 402-10.1186/1471-2164-12-402.PubMed CentralPubMedView ArticleGoogle Scholar
- Huang Y, Niu B, Gao Y, Fu L, Li W: CD-HIT Suite: a web server for clustering and comparing biological sequences. Bioinformatics. 2010, 26: 680-682. 10.1093/bioinformatics/btq003.PubMed CentralPubMedView ArticleGoogle Scholar
- Team TRDC: R: A language and environment for statistical computing. 2008, R Foundation for Statistical Computing: Vienna, Austria, Retrieved from http://cran.r-project.org/Google Scholar
- Warnes G, Bolker B, Lumley T: gplots: Various R programming tools for 10978 Q12 plotting data. 2014, R package, http://CRAN.R-project.org/,Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.