Whole-genome sequence variation among multiple isolates of Pseudomonas aeruginosa

Whole-genome shotgun sequencing was used to study the sequence variation of three Pseudomonas aeruginosa isolates, two from clonal infections of cystic fibrosis patients and one from an aquatic environment, relative to the genomic sequence of reference strain PAO1. The majority of the PAO1 genome is represented in these strains; however, at least three prominent islands of PAO1-specific sequence are apparent. Conversely, approximately 10% of the sequencing reads derived from each isolate fail to align with the PAO1 backbone. While average sequence variation among all strains is roughly 0.5%, regions of pronounced differences were evident in whole-genome scans of nucleotide diversity. We analyzed two such divergent loci, the pyoverdine and O-antigen biosynthesis regions, by complete resequencing. A thorough analysis of isolates collected over time from one of the cystic fibrosis patients revealed independent mutations resulting in the loss of O-antigen synthesis alternating with a mucoid phenotype. Overall, we conclude that most of the PAO1 genome represents a core P. aeruginosa backbone sequence while the strains addressed in this study possess additional genetic material that accounts for at least 10% of their genomes. Approximately half of these additional sequences are novel.     

Genomic Variation among Contemporary Pseudomonas aeruginosa Isolates from Chronically Infected Cystic Fibrosis Patients

The airways of individuals with cystic fibrosis (CF) often become chronically infected with unique strains of the opportunistic pathogen Pseudomonas aeruginosa. Several lines of evidence suggest that the infecting P. aeruginosa lineage diversifies in the CF lung niche, yet so far this contemporary diversity has not been investigated at a genomic level. In this work, we sequenced the genomes of pairs of randomly selected contemporary isolates sampled from the expectorated sputum of three chronically infected adult CF patients. Each patient was infected by a distinct strain of P. aeruginosa. Single nucleotide polymorphisms (SNPs) and insertions/deletions (indels) were identified in the DNA common to the paired isolates from different patients. The paired isolates from one patient differed due to just 1 SNP and 8 indels. The paired isolates from a second patient differed due to 54 SNPs and 38 indels. The pair of isolates from the third patient both contained a mutS mutation, which conferred a hypermutator phenotype; these isolates cumulatively differed due to 344 SNPs and 93 indels. In two of the pairs of isolates, a different accessory genome composition, specifically integrated prophage, was identified in one but not the other isolate of each pair. We conclude that contemporary isolates from a single sputum sample can differ at the SNP, indel, and accessory genome levels and that the cross-sectional genomic variation among coeval pairs of P. aeruginosa CF isolates can be comparable to the variation previously reported to differentiate between paired longitudinally sampled isolates.     

Genome mosaicism is conserved but not unique in Pseudomonas aeruginosa isolates from the airways of young children with cystic fibrosis

Pseudomonas aeruginosa strains from the chronic lung infections of cystic fibrosis (CF) patients are phenotypically and genotypically diverse. Using strain PAO1 whole genome DNA microarrays, we assessed the genomic variation in P. aeruginosa strains isolated from young children with CF (6 months to 8 years of age) as well as from the environment. Eighty-nine to 97% of the PAO1 open reading frames were detected in 20 strains by microarray analysis, while subsets of 38 gene islands were absent or divergent. No specific pattern of genome mosaicism defined strains associated with CF. Many mosaic regions were distinguished by their low G + C content; their inclusion of phage related or pyocin genes; or by their linkage to a vgr gene or a tRNA gene. Microarray and phenotypic analysis of sequential isolates from individual patients revealed two deletions of greater than 100 kbp formed during evolution in the lung. The gene loss in these sequential isolates raises the possibility that acquisition of pyomelanin production and loss of pyoverdin uptake each may be of adaptive significance. Further characterization of P. aeruginosa diversity within the airways of individual CF patients may reveal common adaptations, perhaps mediated by gene loss, that suggest new opportunities for therapy.     

Inter- and intraclonal diversity of the Pseudomonas aeruginosa proteome manifests within the secretome

The proteomes of cultured Pseudomonas aeruginosa isolates from chronically infected cystic fibrosis (CF) lungs were compared by using genetically divergent clones and isogenic morphotypes of one strain. Cellular extracts gave very similar protein patterns in two-dimensional gels, suggesting that the conserved species-specific core genome encodes proteins that are expressed under standard culture conditions in vitro. In contrast, the protein profiles of extracts of culture supernatants were dependent on the growth phase, and there were significant differences between clones. The profiles also varied within clonally related morphotypes from one CF patient, including a hyperpiliated small-colony variant. Mass spectrometry revealed that this variant overexpressed proteins secreted by the type I secretion system (including proteins involved in iron acquisition) and by the type III secretion system. Furthermore, the proteins in the supernatant extracts from the small-colony variant which were recognized by sera from different CF patients varied greatly. We concluded that the secretome expression is a sensitive measure of P. aeruginosa strain variation.     

The Genome of Pseudomonas fluorescens Strain R124 Demonstrates Phenotypic Adaptation to the Mineral Environment

Microbial adaptation to environmental conditions is a complex process, including acquisition of positive traits through horizontal gene transfer or the modification of existing genes through duplication and/or mutation. In this study, we examined the adaptation of a Pseudomonas fluorescens isolate (R124) from the nutrient-limited mineral environment of a silica cave in comparison with P. fluorescens isolates from surface soil and the rhizosphere. Examination of metal homeostasis gene pathways demonstrated a high degree of conservation, suggesting that such systems remain functionally similar across chemical environments. The examination of genomic islands unique to our strain revealed the presence of genes involved in carbohydrate metabolism, aromatic carbon metabolism, and carbon turnover, confirmed through phenotypic assays, suggesting the acquisition of potentially novel mechanisms for energy metabolism in this strain. We also identified a twitching motility phenotype active at low-nutrient concentrations that may allow alternative exploratory mechanisms for this organism in a geochemical environment. Two sets of candidate twitching motility genes are present within the genome, one on the chromosome and one on a plasmid; however, a plasmid knockout identified the functional gene as being present on the chromosome. This work highlights the plasticity of the Pseudomonas genome, allowing the acquisition of novel nutrient-scavenging pathways across diverse geochemical environments while maintaining a core of functional stress response genes.     

Cystic fibrosis-niche adaptation of Pseudomonas aeruginosa reduces virulence in multiple infection hosts

The opportunistic pathogen Pseudomonas aeruginosa is able to thrive in diverse ecological niches and to cause serious human infection. P. aeruginosa environmental strains are producing various virulence factors that are required for establishing acute infections in several host organisms; however, the P. aeruginosa phenotypic variants favour long-term persistence in the cystic fibrosis (CF) airways. Whether P. aeruginosa strains, which have adapted to the CF-niche, have lost their competitive fitness in the other environment remains to be investigated. In this paper, three P. aeruginosa clonal lineages, including early strains isolated at the onset of infection, and late strains, isolated after several years of chronic lung infection from patients with CF, were analysed in multi-host model systems of acute infection. P. aeruginosa early isolates caused lethality in the three non-mammalian hosts, namely Caenorhabditis elegans, Galleria mellonella, and Drosophila melanogaster, while late adapted clonal isolates were attenuated in acute virulence. When two different mouse genetic background strains, namely C57Bl/6NCrl and Balb/cAnNCrl, were used as acute infection models, early P. aeruginosa CF isolates were lethal, while late isolates exhibited reduced or abolished acute virulence. Severe histopathological lesions, including high leukocytes recruitment and bacterial load, were detected in the lungs of mice infected with P. aeruginosa CF early isolates, while late isolates were progressively cleared. In addition, systemic bacterial spread and invasion of epithelial cells, which were detected for P. aeruginosa CF early strains, were not observed with late strains. Our findings indicate that niche-specific selection in P. aeruginosa reduced its ability to cause acute infections across a broad range of hosts while maintaining the capacity for chronic infection in the CF host.     

A novel Pseudomonas aeruginosa Bacteriophage,Ab31, a Chimera Formed from Temperate Phage PAJU2 and P. putida Lytic Phage AF: Characteristics and Mechanism of Bacterial Resistance

A novel temperate bacteriophage of Pseudomonas aeruginosa, phage vB_PaeP_Tr60_Ab31 (alias Ab31) is described. Its genome is composed of structural genes related to those of lytic P. putida phage AF, and regulatory genes similar to those of temperate phage PAJU2. The virion structure resembles that of phage AF and other lytic Podoviridae (S. enterica Epsilon 15 and E. coli phiv10) with similar tail spikes. Ab31 was able to infect P. aeruginosa strain PA14 and two genetically related strains called Tr60 and Tr162, out of 35 diverse strains from cystic fibrosis patients. Analysis of resistant host variants revealed different phenotypes, including induction of pigment and alginate overproduction. Whole genome sequencing of resistant variants highlighted the existence of a large deletion of 234 kbp in two strains, encompassing a cluster of genes required for the production of CupA fimbriae. Stable lysogens formed by Ab31 in strain Tr60, permitted the identification of the insertion site. During colonization of the lung in cystic fibrosis patients, P. aeruginosa adapts by modifying its genome. We suggest that bacteriophages such as Ab31 may play an important role in this adaptation by selecting for bacterial characteristics that favor persistence of bacteria in the lung.     

Diversity of the abundant pKLC102/PAGI-2 family of genomic islands in Pseudomonas aeruginosa

The known genomic islands of Pseudomonas aeruginosa clone C strains are integrated into tRNA(Lys) (pKLC102) or tRNA(Gly) (PAGI-2 and PAGI-3) genes and differ from their core genomes by distinctive tetranucleotide usage patterns. pKLC102 and the related island PAPI-1 from P. aeruginosa PA14 were spontaneously mobilized from their host chromosomes at frequencies of 10% and 0.3%, making pKLC102 the most mobile genomic island known with a copy number of 30 episomal circular pKLC102 molecules per cell. The incidence of islands of the pKLC102/PAGI-2 type was investigated in 71 unrelated P. aeruginosa strains from diverse habitats and geographic origins. pKLC102- and PAGI-2-like islands were identified in 50 and 31 strains, respectively, and 15 and 10 subtypes were differentiated by hybridization on pKLC102 and PAGI-2 macroarrays. The diversity of PAGI-2-type islands was mainly caused by one large block of strain-specific genes, whereas the diversity of pKLC102-type islands was primarily generated by subtype-specific combination of gene cassettes. Chromosomal loss of PAGI-2 could be documented in sequential P. aeruginosa isolates from individuals with cystic fibrosis. PAGI-2 was present in most tested Cupriavidus metallidurans and Cupriavidus campinensis isolates from polluted environments, demonstrating the spread of PAGI-2 across habitats and species barriers. The pKLC102/PAGI-2 family is prevalent in numerous beta- and gammaproteobacteria and is characterized by high asymmetry of the cDNA strands. This evolutionarily ancient family of genomic islands retained its oligonucleotide signature during horizontal spread within and among taxa.     

Whole genome shotgun sequencing of one Colombian clinical isolate of Mycobacterium tuberculosis reveals DosR regulon gene deletions

Several genomes of different Mycobacterium tuberculosis isolates have been completely sequenced around the world. The genomic information obtained have shown higher diversity than originally thought and specific adaptations to different human populations. Within this work, we sequenced the genome of one Colombian M.?tuberculosis virulent isolate. Genomic comparison against the reference genome of H37Rv and other strains showed multiple deletion and insertions that ranged between a few bases to thousands. Excluding PPE and PG-PGRS genes, 430 proteins present changes in at least 1 amino acid. Also, novel positions of the IS6110 mobile element were identified. This isolate is also characterized by a large genomic deletion of 3.6?kb, leading to the loss and modification of the dosR regulon genes, Rv1996 and Rv1997. To our knowledge, this is the first report of the genome sequence of a Latin American M.?tuberculosis clinical isolate.     

Microevolution of the major common Pseudomonas aeruginosa clones C and PA14 in cystic fibrosis lungs

Clones C and PA14 are the worldwide most abundant clonal complexes in the Pseudomonas aeruginosa population. The microevolution of clones C and PA14 was investigated in serial cystic fibrosis (CF) airway isolates collected over 20 years since the onset of colonization. Intraclonal evolution in CF lungs was resolved by genome sequencing of first, intermediate and late isolates and subsequent multimarker SNP genotyping of the whole strain panel. Mapping of sequence reads onto the P. aeruginosa PA14 reference genome unravelled an intraclonal and interclonal sequence diversity of 0.0035% and 0.68% respectively. Clone PA14 diversified into three branches in the patient's lungs, and the PA14 population acquired 15 nucleotide substitutions and a large deletion during the observation period. The clone C genome remained invariant during the first 3 years in CF lungs; however, 15 years later 947 transitions and 12 transversions were detected in a clone C mutL mutant strain. Key mutations occurred in retS, RNA polymerase, multidrug transporter, virulence and denitrification genes. Late clone C and PA14 persistors in the CF lungs were compromised in growth and cytotoxicity, but their mutation frequency was normal even in mutL mutant clades.     

Genome sequence analyses of Pseudomonas savastanoi pv. glycinea and subtractive hybridization-based comparative genomics with nine pseudomonads

Bacterial blight, caused by Pseudomonas savastanoi pv. glycinea (Psg), is a common disease of soybean. In an effort to compare a current field isolate with one isolated in the early 1960s, the genomes of two Psg strains, race 4 and B076, were sequenced using 454 pyrosequencing. The genomes of both Psg strains share more than 4,900 highly conserved genes, indicating very low genetic diversity between Psg genomes. Though conserved, genome rearrangements and recombination events occur commonly within the two Psg genomes. When compared to each other, 437 and 163 specific genes were identified in B076 and race 4, respectively. Most specific genes are plasmid-borne, indicating that acquisition and maintenance of plasmids may represent a major mechanism to change the genetic composition of the genome and even acquire new virulence factors. Type three secretion gene clusters of Psg strains are near identical with that of P. savastanoi pv. phaseolicola (Pph) strain 1448A and they shared 20 common effector genes. Furthermore, the coronatine biosynthetic cluster is present on a large plasmid in strain B076, but not in race 4. In silico subtractive hybridization-based comparative genomic analyses with nine sequenced phytopathogenic pseudomonads identified dozens of specific islands (SIs), and revealed that the genomes of Psg strains are more similar to those belonging to the same genomospecies such as Pph 1448A than to other phytopathogenic pseudomonads. The number of highly conserved genes (core genome) among them decreased dramatically when more genomes were included in the subtraction, suggesting the diversification of pseudomonads, and further indicating the genome heterogeneity among pseudomonads. However, the number of specific genes did not change significantly, suggesting these genes are indeed specific in Psg genomes. These results reinforce the idea of a species complex of P. syringae and support the reclassification of P. syringae into different species.     

A major Pseudomonas aeruginosa clone common to patients and aquatic habitats.

The genomic relatedness of 573 Pseudomonas aeruginosa strains from environmental and clinical habitats was examined by digesting the genome with the rare-cutting enzyme SpeI. Thirty-nine strains were collected from environmental habitats mainly of aquatic origin, like rivers, lakes, or sanitary facilities. Four hundred fifty strains were collected from 76 patients with cystic fibrosis (CF) treated at four different centers, and 25 additional clinical isolates were collected from patients suffering from other diseases. Twenty-nine P. aeruginosa isolates were collected from the environment of one CF clinic. Thirty strains from culture collections were of environmental and clinic origin. A common macrorestriction fingerprint pattern was found in 13 of 46 CF patients, 5 of 29 environmental isolates from the same hospital, in a single ear infection isolate from another hospital, and 8 of 38 isolates from aquatic habitats about 300 km away from the CF clinic. The data indicate that closely related variants of one major clone (called clone C) persisted in various spatially and temporally separated habitats. Southern analysis of the clonal variants with six gene probes and two probes for genes coding for rRNA revealed almost the same hybridization patterns. With the exception of the phenotypically rapidly evolving CF isolates, the close relatedness of the strains of the clone was also shown by their identical responses in pyocin typing, phage typing, and serotyping. Besides clone C, three other P. aeruginosa clones were isolated from more than one clinical or environmental source.     

Emergence of Pseudomonas aeruginosa strains producing high levels of persister cells in patients with cystic fibrosis

The majority of cystic fibrosis (CF) patients succumb to a chronic infection of the airway with Pseudomonas aeruginosa. Paradoxically, pathogenic strains are often susceptible to antibiotics, but the infection cannot be eradicated with antimicrobial therapy. We find that in a majority of patients with airway infections, late isolates of P. aeruginosa produce increased levels of drug-tolerant persister cells. The genomes of a clonal pair of early/late isolates from a single patient have been previously sequenced, and the late isolate (obtained at age 96 months) showed a 100-fold increase in persister levels. The 96-month isolate carries a large number of mutations, including a mutation in mutS that confers a hypermutator phenotype. There is also a mutation in the mexZ repressor controlling the expression of the MexXY-OprM multidrug pump, which results in a moderate increase in the ofloxacin, carbenicillin, and tobramycin MICs. Knocking out the mexXY locus restored the resistance to that of the parent strain but did not affect the high levels of persisters formed by the 96-month isolate. This suggests that the late isolate is a high-persister (hip) mutant. Increased persister formation was observed in exponential phase, stationary phase, and biofilm populations of the 96-month isolate. Analysis of late isolates from 14 additional patients indicated that 10 of them are hip mutants. Most of these hip mutants did not have higher drug resistance. Increased persister formation appears to be their sole mechanism for surviving chemotherapy. Taken together, these findings suggest a link between persisters and recalcitrance of CF infection and identify an overlooked culprit-high-persister mutants producing elevated levels of drug-tolerant cells. Persisters may play a similarly critical role in the recalcitrance of other chronic infections.     

An intragenic deletion in pilQ leads to nonpiliation of a Pseudomonas aeruginosa strain isolated from cystic fibrosis lung

Deficient motility is one of the characteristic hallmarks observed in Pseudomonas aeruginosa strains that chronically colonize the lungs of cystic fibrosis (CF) patients. Pseudomonas aeruginosa TB is a nonpiliated CF isolate known to be defective in twitching motility. Complementation confirmed a direct link of this phenotype to an intragenic out-of-frame deletion in pilQ (PA5040). Sequence alignment of pilQ derived from TB vs. PAO1 suggests that close direct repeats framing the deletion site may have triggered this mutation. This type of mutation could play a role in the emergence of pathoadaptive mutations of P. aeruginosa in the CF lung habitat.     

Genomics of an endemic cystic fibrosis Burkholderia multivorans strain reveals low within-patient evolution but high between-patient diversity

Burkholderia multivorans is a member of the Burkholderia cepacia complex (Bcc), notorious for its pathogenicity in persons with cystic fibrosis. Epidemiological surveillance suggests that patients predominantly acquire B. multivorans from environmental sources, with rare cases of patient-to-patient transmission. Here we report on the genomic analysis of thirteen isolates from an endemic B. multivorans strain infecting four cystic fibrosis patients treated in different pediatric cystic fibrosis centers in Belgium, with no evidence of cross-infection. All isolates share an identical sequence type (ST-742) but whole genome analysis shows that they exhibit peculiar patterns of genomic diversity between patients. By combining short and long reads sequencing technologies, we highlight key differences in terms of small nucleotide polymorphisms indicative of low rates of adaptive evolution within patient, and well-defined, hundred kbps-long segments of high enrichment in mutations between patients. In addition, we observed large structural genomic variations amongst the isolates which revealed different plasmid contents, active roles for transposase IS3 and IS5 in the deactivation of genes, and mobile prophage elements. Our study shows limited within-patient B. multivorans evolution and high between-patient strain diversity, indicating that an environmental microdiverse reservoir must be present for this endemic strain, in which active diversification is taking place. Furthermore, our analysis also reveals a set of 30 parallel adaptations across multiple patients, indicating that the specific genomic background of a given strain may dictate the route of adaptation within the cystic fibrosis lung.     

Purification of outer membrane vesicles from Pseudomonas aeruginosa and their activation of an IL-8 response

Considerable lung injury results from the inflammatory response to Pseudomonas aeruginosa infections in patients with cystic fibrosis (CF). The P. aeruginosa laboratory strain PAO1, an environmental isolate, and isolates from CF patients were cultured in vitro and outer membrane vesicles from those cultures were quantitated, purified, and characterized. Vesicles were produced throughout the growth phases of the culture and vesicle yield was strain-independent. Strain-dependent differences in the protein composition of vesicles were quantitated and identified. The aminopeptidase PaAP (PA2939) was highly enriched in vesicles from CF isolates. Vesicles from all strains elicited IL-8 secretion by lung epithelial cells. These results suggest that P. aeruginosa colonizing the CF lung may produce vesicles with a particular composition and that the vesicles could contribute to inflammation.     

Impact of large chromosomal inversions on the adaptation and evolution of Pseudomonas aeruginosa chronically colonizing cystic fibrosis lungs

Pseudomonas aeruginosa chronically colonizing the lungs of cystic fibrosis (CF) patients undergoes fast evolution leading to clonal divergence. More than half of the genotypes of P. aeruginosa clone C isolates exclusively from CF lung infection exhibit large chromosomal inversions (LCIs). To analyse the impact of LCIs, as a novel mechanism of bacterial adaptation, the underlying molecular mechanism was examined. Analysis of inversion breakpoints suggested an IS6100-induced coupled insertion-inversion mechanism. A selective advantage was created by insertion of IS6100 into wbpM, pilB and mutS which leads to common CF phenotypes such as O-antigen and type IV pili deficiency and hypermutability. Speciation in bacteria is accompanied by LCIs. Therefore adaptation by LCIs that allows persistence of P. aeruginosa in the CF lung and species diversification in that new ecological niche can serve as a model for bacterial genome evolution.