Core and accessory genome architecture in a group of Pseudomonas aeruginosa Mu-like phages

Background

Bacteriophages that infect the opportunistic pathogen Pseudomonas aeruginosa have been classified into several groups. One of them, which includes temperate phage particles with icosahedral heads and long flexible tails, bears genomes whose architecture and replication mechanism, but not their nucleotide sequences, are like those of coliphage Mu. By comparing the genomic sequences of this group of P. aeruginosa phages one could draw conclusions about their ontogeny and evolution.

Results

Two newly isolated Mu-like phages of P. aeruginosa are described and their genomes sequenced and compared with those available in the public data banks. The genome sequences of the two phages are similar to each other and to those of a group of P. aeruginosa transposable phages. Comparing twelve of these genomes revealed a common genomic architecture in the group. Each phage genome had numerous genes with homologues in all the other genomes and a set of variable genes specific for each genome. The first group, which comprised most of the genes with assigned functions, was named “core genome”, and the second group, containing mostly short ORFs without assigned functions was called “accessory genome”. Like in other phage groups, variable genes are confined to specific regions in the genome.

Conclusion

Based on the known and inferred functions for some of the variable genes of the phages analyzed here, they appear to confer selective advantages for the phage survival under particular host conditions. We speculate that phages have developed a mechanism for horizontally acquiring genes to incorporate them at specific loci in the genome that help phage adaptation to the selective pressures imposed by the host.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1146) contains supplementary material, which is available to authorized users.     

The Susceptibility of Pseudomonas aeruginosa Strains from Cystic Fibrosis Patients to Bacteriophages

Phage therapy may become a complement to antibiotics in the treatment of chronic Pseudomonas aeruginosa infection. To design efficient therapeutic cocktails, the genetic diversity of the species and the spectrum of susceptibility to bacteriophages must be investigated. Bacterial strains showing high levels of phage resistance need to be identified in order to decipher the underlying mechanisms. Here we have selected genetically diverse P. aeruginosa strains from cystic fibrosis patients and tested their susceptibility to a large collection of phages. Based on plaque morphology and restriction profiles, six different phages were purified from “pyophage”, a commercial cocktail directed against five different bacterial species, including P. aeruginosa. Characterization of these phages by electron microscopy and sequencing of genome fragments showed that they belong to 4 different genera. Among 47 P. aeruginosa strains, 13 were not lysed by any of the isolated phages individually or by pyophage. We isolated two new phages that could lyse some of these strains, and their genomes were sequenced. The presence/absence of a CRISPR-Cas system (Clustered Regularly Interspaced Short Palindromic Repeats and Crisper associated genes) was investigated to evaluate the role of the system in phage resistance. Altogether, the results show that some P. aeruginosa strains cannot support the growth of any of the tested phages belonging to 5 different genera, and suggest that the CRISPR-Cas system is not a major defence mechanism against these lytic phages.     

Genome landscapes and bacteriophage codon usage

Across all kingdoms of biological life, protein-coding genes exhibit unequal usage of synonymous codons. Although alternative theories abound, translational selection has been accepted as an important mechanism that shapes the patterns of codon usage in prokaryotes and simple eukaryotes. Here we analyze patterns of codon usage across 74 diverse bacteriophages that infect E. coli, P. aeruginosa, and L. lactis as their primary host. We use the concept of a “genome landscape,” which helps reveal non-trivial, long-range patterns in codon usage across a genome. We develop a series of randomization tests that allow us to interrogate the significance of one aspect of codon usage, such as GC content, while controlling for another aspect, such as adaptation to host-preferred codons. We find that 33 phage genomes exhibit highly non-random patterns in their GC3-content, use of host-preferred codons, or both. We show that the head and tail proteins of these phages exhibit significant bias towards host-preferred codons, relative to the non-structural phage proteins. Our results support the hypothesis of translational selection on viral genes for host-preferred codons, over a broad range of bacteriophages.     

The isolation and characterization of two Stenotrophomonas maltophilia bacteriophages capable of cross-taxonomic order infectivity

Background

A rapid worldwide increase in the number of human infections caused by the extremely antibiotic resistant bacterium Stenotrophomonas maltophilia is prompting alarm. One potential treatment solution to the current antibiotic resistance dilemma is “phage therapy”, the clinical application of bacteriophages to selectively kill bacteria.

Results

Towards that end, phages DLP1 and DLP2 (vB_SmaS-DLP_1 and vB_SmaS-DLP_2, respectively) were isolated against S. maltophilia strain D1585. Host range analysis for each phage was conducted using 27 clinical S. maltophilia isolates and 11 Pseudomonas aeruginosa strains. Both phages exhibit unusually broad host ranges capable of infecting bacteria across taxonomic orders. Transmission electron microscopy of the phage DLP1 and DLP2 morphology reveals that they belong to the Siphoviridae family of bacteriophages. Restriction fragment length polymorphism analysis and complete genome sequencing and analysis indicates that phages DLP1 and DLP2 are closely related but different phages, sharing 96.7 % identity over 97.2 % of their genomes. These two phages are also related to P. aeruginosa phages vB_Pae-Kakheti_25 (PA25), PA73, and vB_PaeS_SCH_Ab26 (Ab26) and more distantly related to Burkholderia cepacia complex phage KL1, which together make up a taxonomic sub-family. Phages DLP1 and DLP2 exhibited significant differences in host ranges and growth kinetics.

Conclusions

The isolation and characterization of phages able to infect two completely different species of bacteria is an exciting discovery, as phages typically can only infect related bacterial species, and rarely infect bacteria across taxonomic families, let alone across taxonomic orders.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1848-y) contains supplementary material, which is available to authorized users.     

Isolation and Characterization of a “phiKMV-Like” Bacteriophage and Its Therapeutic Effect on Mink Hemorrhagic Pneumonia

The objective of this study was to investigate the potential of using phages as a therapy against hemorrhagic pneumonia in mink both in vitro and in vivo. Five Pseudomonas aeruginosa (P. aeruginosa) strains were isolated from lungs of mink with suspected hemorrhagic pneumonia and their identity was confirmed by morphological observation and 16S rDNA sequence analysis. Compared to P. aeruginosa strains isolated from mink with hemorrhagic pneumonia in 2002, these isolates were more resistant to antibiotics selected. A lytic phage vB_PaeP_PPA-ABTNL (PPA-ABTNL) of the Podoviridae family was isolated from hospital sewage using a P. aeruginosa isolate as host, showing broad host range against P. aeruginosa. A one-step growth curve analysis of PPA-ABTNL revealed eclipse and latent periods of 20 and 35 min, respectively, with a burst size of about 110 PFU per infected cell. Phage PPA-ABTNL significantly reduced the growth of P. aeruginosa isolates in vitro. The genome of PPA-ABTNL was 43,227 bp (62.4% G+C) containing 54 open reading frames and lacked regions encoding known virulence factors, integration-related proteins and antibiotic resistance determinants. Genome architecture analysis showed that PPA-ABTNL belonged to the “phiKMV-like Viruses” group. A repeated dose inhalational toxicity study using PPA-ABTNL crude preparation was conducted in mice and no significantly abnormal histological changes, morbidity or mortality were observed. There was no indication of any potential risk associated with using PPA-ABTNL as a therapeutic agent. The results of a curative treatment experiment demonstrated that atomization by ultrasonic treatment could efficiently deliver phage to the lungs of mink and a dose of 10 multiplicity of infection was optimal for treating mink hemorrhagic pneumonia. Our work demonstrated the potential for phage to fight P. aeruginosa involved in mink lung infections when administered by means of ultrasonic nebulization.     

Genome-driven elucidation of phage-host interplay and impact of phage resistance evolution on bacterial fitness

When considering the interactions between bacteriophages and their host, the issue of phage-resistance emergence is a key element in understanding the ecological impact of phages on the bacterial population. It is also an essential parameter for the implementation of phage therapy to combat antibiotic-resistant pathogens. This study investigates the phenotypic and genetic responses of five Pseudomonas aeruginosa strains (PAO1, A5803, AA43, CHA, and PAK) to the infection by seven phages with distinct evolutionary backgrounds and recognised receptors (LPS/T4P). Emerging phage-insensitivity was generally accompanied by self and cross-resistance mechanisms. Significant differences were observed between the reference PAO1 responses compared to other clinical representatives. LPS-dependent phage infections in clinical strains selected for mutations in the “global regulatory” and “other” genes, rather than in the LPS-synthesis clusters detected in PAO1 clones. Reduced fitness, as proxied by the growth rate, was correlated with large deletion (20–500 kbp) and phage carrier state. Multi-phage resistance was significantly correlated with a reduced growth rate but only in the PAO1 population. In addition, we observed that the presence of prophages decreased the lytic phage maintenance seemingly protecting the host against carrier state and occasional lytic phage propagation, thus preventing a significant reduction in bacterial growth rate.Subject terms: Bacteriophages, Biodiversity     

Epidemiological Fingerprinting of Pseudomonas aeruginosa by the Production of and Sensitivity to Pyocin and Bacteriophage

A new method has been devised to trace cross-infection by Pseudomonas aeruginosa. Unknown strains growing logarithmically in liquid media were treated with mitomycin C to induce the liberation of pyocin and phage. The lysates were then tested against 27 selected indicator strains, and the zones of clearing were differentiated as to killing by pyocin or lysis by phage. Twenty-four standard pyocin-phage lysates were then applied to each of the unknowns, and the sensitivity pattern was recorded. Thus, an “epidemiological fingerprint” consisting of 51 operational characteristics was established for each isolate. Organisms from the same source had identical or similar fingerprints, but organisms from different origins could easily be distinguished. Pyocin production, pyocin sensitivity, and phage production were found to be stable genetic characters; however, spontaneous mutations in phage sensitivity were frequently encountered. The epidemiological fingerprint has proven to be a sensitive tool in establishing the identity or dissimilarity of unknown strains. This method has been of great value in tracing the epidemiology of P. aeruginosa in the hospital environment. Each of the 157 P. aeruginosa strains tested has been typable by this method.     

Characterization of the core and accessory genomes of Pseudomonas aeruginosa using bioinformatic tools Spine and AGEnt

Background

Pseudomonas aeruginosa is an important opportunistic pathogen responsible for many infections in hospitalized and immunocompromised patients. Previous reports estimated that approximately 10% of its 6.6 Mbp genome varies from strain to strain and is therefore referred to as “accessory genome”. Elements within the accessory genome of P. aeruginosa have been associated with differences in virulence and antibiotic resistance. As whole genome sequencing of bacterial strains becomes more widespread and cost-effective, methods to quickly and reliably identify accessory genomic elements in newly sequenced P. aeruginosa genomes will be needed.

Results

We developed a bioinformatic method for identifying the accessory genome of P. aeruginosa. First, the core genome was determined based on sequence conserved among the completed genomes of twelve reference strains using Spine, a software program developed for this purpose. The core genome was 5.84 Mbp in size and contained 5,316 coding sequences. We then developed an in silico genome subtraction program named AGEnt to filter out core genomic sequences from P. aeruginosa whole genomes to identify accessory genomic sequences of these reference strains. This analysis determined that the accessory genome of P. aeruginosa ranged from 6.9-18.0% of the total genome, was enriched for genes associated with mobile elements, and was comprised of a majority of genes with unknown or unclear function. Using these genomes, we showed that AGEnt performed well compared to other publically available programs designed to detect accessory genomic elements. We then demonstrated the utility of the AGEnt program by applying it to the draft genomes of two previously unsequenced P. aeruginosa strains, PA99 and PA103.

Conclusions

The P. aeruginosa genome is rich in accessory genetic material. The AGEnt program accurately identified the accessory genomes of newly sequenced P. aeruginosa strains, even when draft genomes were used. As P. aeruginosa genomes become available at an increasingly rapid pace, this program will be useful in cataloging the expanding accessory genome of this bacterium and in discerning correlations between phenotype and accessory genome makeup. The combination of Spine and AGEnt should be useful in defining the accessory genomes of other bacterial species as well.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-737) contains supplementary material, which is available to authorized users.     

Phenotypic and Genotypic Characteristics of Neisseria meningitidis Disease-Causing Strains in Argentina, 2010

Phenotypic and genotypic characterization of 133 isolates of Neisseria meningitidis obtained from meningococcal disease cases in Argentina during 2010 were performed by the National Reference Laboratory as part of a project coordinated by the PAHO within the SIREVA II network. Serogroup, serotype, serosubtype and MLST characterization were performed. Minimum Inhibitory Concentration to penicillin, ampicillin, ceftriaxone, rifampin, chloramphenicol, tetracycline and ciprofloxacin were determined and interpreted according to CLSI guidelines. Almost 49% of isolates were W135, and two serotype:serosubtype combinations, W135∶2a:P1.5,2:ST-11 and W135∶2a:P1.2:ST-11 accounted for 78% of all W135 isolates. Serogroup B accounted for 42.1% of isolates, and was both phenotypically and genotypically diverse. Serogroup C isolates represented 5.3% of the dataset, and one isolate belonging to the ST-198 complex was non-groupable. Isolates belonged mainly to the ST-11 complex (48%) and to a lesser extent to the ST-865 (18%), ST-32 (9,8%) and the ST-35 complexes (9%). Intermediate resistance to penicillin and ampicillin was detected in 35.4% and 33.1% of isolates respectively. Two W135∶2a:P1.5,2:ST-11:ST-11 isolates presented resistance to ciprofloxacin associated with a mutation in the QRDR of gyrA gene Thr91-Ile. These data show serogroup W135 was the first cause of disease in Argentina in 2010, and was strongly associated with the W135∶2a:P1.5,2:ST-11 epidemic clone. Serogroup B was the second cause of disease and isolates belonging to this serogroup were phenotypically and genotypically diverse. The presence of isolates with intermediate resistance to penicillin and the presence of fluorquinolone-resistant isolates highlight the necessity and importance of maintaining and strengthening National Surveillance Programs.     

Three Prochlorococcus cyanophage genomes: signature features and ecological interpretations

The oceanic cyanobacteria Prochlorococcus are globally important, ecologically diverse primary producers. It is thought that their viruses (phages) mediate population sizes and affect the evolutionary trajectories of their hosts. Here we present an analysis of genomes from three Prochlorococcus phages: a podovirus and two myoviruses. The morphology, overall genome features, and gene content of these phages suggest that they are quite similar to T7-like (P-SSP7) and T4-like (P-SSM2 and P-SSM4) phages. Using the existing phage taxonomic framework as a guideline, we examined genome sequences to establish “core” genes for each phage group. We found the podovirus contained 15 of 26 core T7-like genes and the two myoviruses contained 43 and 42 of 75 core T4-like genes. In addition to these core genes, each genome contains a significant number of “cyanobacterial” genes, i.e., genes with significant best BLAST hits to genes found in cyanobacteria. Some of these, we speculate, represent “signature” cyanophage genes. For example, all three phage genomes contain photosynthetic genes (psbA, hliP) that are thought to help maintain host photosynthetic activity during infection, as well as an aldolase family gene (talC) that could facilitate alternative routes of carbon metabolism during infection. The podovirus genome also contains an integrase gene (int) and other features that suggest it is capable of integrating into its host. If indeed it is, this would be unprecedented among cultured T7-like phages or marine cyanophages and would have significant evolutionary and ecological implications for phage and host. Further, both myoviruses contain phosphate-inducible genes (phoH and pstS) that are likely to be important for phage and host responses to phosphate stress, a commonly limiting nutrient in marine systems. Thus, these marine cyanophages appear to be variations of two well-known phages—T7 and T4—but contain genes that, if functional, reflect adaptations for infection of photosynthetic hosts in low-nutrient oceanic environments.     

Romulus and Remus,Two Phage Isolates Representing a Distinct Clade within the Twortlikevirus Genus,Display Suitable Properties for Phage Therapy Applications

The renewed interest in controlling Staphylococcus aureus infections using their natural enemies, bacteriophages, has led to the isolation of a limited number of virulent phages so far. These phages are all members of the Twortlikevirus, displaying little variance. We present two novel closely related (95.9% DNA homology) lytic myoviruses, Romulus and Remus, with double-stranded DNA (dsDNA) genomes of 131,333 bp and 134,643 bp, respectively. Despite their relatedness to Staphylococcus phages K, G1, ISP, and Twort and Listeria phages A511 and P100, Romulus and Remus can be proposed as isolates of a new species within the Twortlikevirus genus. A distinguishing feature for these phage genomes is the unique distribution of group I introns compared to that in other staphylococcal myoviruses. In addition, a hedgehog/intein domain was found within their DNA polymerase genes, and an insertion sequence-encoded transposase exhibits splicing behavior and produces a functional portal protein. From a phage therapy application perspective, Romulus and Remus infected approximately 70% of the tested S. aureus isolates and displayed promising lytic activity against these isolates. Furthermore, both phages showed a rapid initial adsorption and demonstrated biofilm-degrading capacity in a proof-of-concept experiment.     

Study of the diversity in a group of phages of <Emphasis Type="Italic">Pseudomonas aeruginosa</Emphasis> species PB1 (Myoviridae) and their behavior in adsorbtion-resistant bacterial mutants

A group of 12 Pseudomonas aeruginosa virulent bacteriophages of different origin scored with regard to the plaque phenotype are assigned to PB1-like species based on the similarity in respect to morphology of particles and high DNA homology. Phages differ in restriction profile and the set of capsid major proteins. For the purpose of studying adsorption properties of these phages, 20 random spontaneous mutants of P. aeruginosa PAO1 with the disturbed adsorption placed in two groups were isolated. Mutants of the first group completely lost the ability to adsorb all phages of this species. It is assumed that their adsorption receptors are functionally inactive or lost at all, because the attempt to isolate phage mutants or detect natural phages of PB1 species capable of overcoming resistance of these bacteria failed. The second group includes five bacterial mutants resistant to the majority of phages belonging to species PB1. These mutants maintain the vigorous growth of phage SN and poor growth of phage 9/3, which forms turbid plaques with low efficiency of plating. In the background of weak growth, phage 9/3 yields plaques that grew well. The examination of the progeny of phage 9/3, which can grow on these bacteria, showed that its DNA differed from DNA of the original phage 9/3 by restriction profile and is identical to DNA of phage PB1 with regard to this trait. Data supported a suggestion that this phage variant resulted from recombination of phage 9/3 DNA with the locus of P. aeruginosa PAO1 genome encoding the bacteriocinogenic factor R. However, this variant of phage 9/3 did not manifest the ability to grow on phage-resistant mutants of the first group. Possible reasons for the difference between phages 9/3 or SN and the remaining phages of PB1 species are discussed. A preliminary formal scheme of the modular structure for adsorption receptors on the surface of P. aeruginosa PAO1 bacteria was constructed based on the analysis of growth of some other phage species on adsorption mutants of the first type.     

Whole-genome comparison of two Campylobacter jejuni isolates of the same sequence type reveals multiple loci of different ancestral lineage

Campylobacter jejuni ST-474 is the most important human enteric pathogen in New Zealand, and yet this genotype is rarely found elsewhere in the world. Insight into the evolution of this organism was gained by a whole genome comparison of two ST-474, flaA SVR-14 isolates and other available C. jejuni isolates and genomes. The two isolates were collected from different sources, human (H22082) and retail poultry (P110b), at the same time and from the same geographical location. Solexa sequencing of each isolate resulted in ~1.659 Mb (H22082) and ~1.656 Mb (P110b) of assembled sequences within 28 (H22082) and 29 (P110b) contigs. We analysed 1502 genes for which we had sequences within both ST-474 isolates and within at least one of 11 C. jejuni reference genomes. Although 94.5% of genes were identical between the two ST-474 isolates, we identified 83 genes that differed by at least one nucleotide, including 55 genes with non-synonymous substitutions. These covered 101 kb and contained 672 point differences. We inferred that 22 (3.3%) of these differences were due to mutation and 650 (96.7%) were imported via recombination. Our analysis estimated 38 recombinant breakpoints within these 83 genes, which correspond to recombination events affecting at least 19 loci regions and gives a tract length estimate of ~2 kb. This includes a ~12 kb region displaying non-homologous recombination in one of the ST-474 genomes, with the insertion of two genes, including ykgC, a putative oxidoreductase, and a conserved hypothetical protein of unknown function. Furthermore, our analysis indicates that the source of this recombined DNA is more likely to have come from C. jejuni strains that are more closely related to ST-474. This suggests that the rates of recombination and mutation are similar in order of magnitude, but that recombination has been much more important for generating divergence between the two ST-474 isolates.     

Genetic and Phenotypic Variability among Salmonella enterica Serovar Typhimurium Isolates from California Dairy Cattle and Humans

Fifty-six human and 24 adult dairy cattle isolates of Salmonella enterica serovar Typhimurium from a single county in California were compared using ribotyping, insertion sequence typing (IS200), pulsed-field gel electrophoresis, plasmid typing, phage typing, and antimicrobial resistance testing. The majority of the isolates fell into one of two groups which were phage types DT104 and DT193. Combining the information from all typing methods, a total of 45 different “clusters” were defined, with 35 of those including only a single isolate. The library of isolates had a high degree of variability, but antibiotic resistance and plasmid typing each defined single clusters in which human or bovine isolates predominated (χ², P < 0.05).     

Effects of bacteriophages on biofilm formation by strains of <Emphasis Type="Italic">Pseudomonas aeruginosa</Emphasis>

The effects of two Pseudomonas aeruginosa bacteriophages, vB-Pa 4 and vB-Pa 5, on the formation and development of biofilms of six polyresistant hospital strains of P. aeruginosa have been investigated. Pretreatment of bacteriophages prevented the formation or almost completely prevented the growth of adequate biofilms. The biofilms that had already formed were partially or completely destroyed after phage treatment. The results demonstrate the prospects of using isolated bacteriophages of P. aeruginosa to destroy biofilms and prevent their formation.     

TL,the new bacteriophage of <Emphasis Type="Italic">Pseudomonas aeruginosa</Emphasis> and its application for the search of halo-producing bacteriophages

The properties of new virulent bacteriophage TL of Pseudomonas aeruginosa belonging to the family Podoviridae (genome size of 46 kb) were investigated. This bacteriophage is capable of lysing the bacterial lawn in halo zones around negative colonies (NC) of other bacteriophages. TL forms large NC, that are hardly distinguishable on the lawn of P. aeruginisa PAO1. At the same time, on the lawns of some phage-resistant PAO1 mutants, as well as on those produced by a number of clinical isolates, TL forms more transparent NC. It is suggested that more effective growth of the bacteriophage TL NC is associated with the differences in outer lipopolysaccharide (LPS) layer of the cell walls of different bacterial strains, as well as of the bacteria inside and outside of the halos. This TL property was used to optimize selection of bacteriophages producing halos around NC on the lawn of P. aeruginosa PAO1. As a result, a group of bacteriophages differing in the patterns of interaction between their halos and TL bacteriophage, as well as in some characters was identified. Taking into consideration the importance of cell-surfaced structures of P. aeruginosa in manifestation of virulence and pathogenicity, possible utilization of specific phage enzymes, polysacchadide depolymerases, for more effective treatment of P. aeruginosa infections is discussed.     

Comparison of new giant bacteriophages OBP and Lu11 of soil pseudomonads with bacteriophages of the ϕKZ-supergroup of Pseudomonas aeruginosa

Study of two recently isolated giant bacteriophages Lu11 and OBP that are active on Pseudomonas putida var. Manila and Pseudomonas fluorescens, respectively, demonstrated their similarity in morphotype, genome size, and size of phage particles, with giant bacteriophages of Pseudomonas aeruginosa assigned to the supergroup of ?KZ-like phages of the family Myoviridae. This supergroup was designated in this manner according to the best studied phage ?KZ that belongs to the species of this group widely distributed in nature. Comparison of major polypeptide sizes of mature particles suggests similarity of certain proteins in the phages examined. In OBP particles visualized with an electron microscope, an “inner body” was detected, which points to specific DNA package intrinsic to phages of ?KZ group. In the meantime, phages Lu11 and OBP do not exhibit resemblance among themselves or with any of earlier described ?KZ-like phages in respect to detectable DNA homology. Note that phage Lu11 of P. putida var. Manila exhibits very slight homology with phage Lin68 of the family of P. aeruginosa ?KZ-like phages detected only in blot hybridization. This suggests the possible involvement of these phages in interspecies recombination (“gene shuffling”) between phages of various bacterial species. Results of partial sequencing of phage genomes confirmed the phylogenetic relatedness of phage OBP to phages of the ?KZ supergroup, whereas phage Lu11 most probably belongs to a novel species that is not a member of supergroup ?KZ composition. The results of the study are discussed in terms of the evolution of these phages.     

Diverse Phage-Encoded Toxins in a Protective Insect Endosymbiont

The lysogenic bacteriophage APSE infects “Candidatus Hamiltonella defensa,” a facultative endosymbiont of aphids and other sap-feeding insects. This endosymbiont has established a beneficial association with aphids, increasing survivorship following attack by parasitoid wasps. Although APSE and “Ca. Hamiltonella defensa” are effectively maternally transmitted between aphid generations, they can also be horizontally transferred among insect hosts, which results in genetically distinct “Ca. Hamiltonella defensa” strains infecting the same aphid species and sporadic distributions of both APSE and “Ca. Hamiltonella defensa” among hosts. Aphids infected only with “Ca. Hamiltonella defensa” have significantly less protection than those infected with both “Ca. Hamiltonella defensa” and APSE. This protection has been proposed to be connected to eukaryote-targeted toxins previously discovered in the genomes of two characterized APSE strains. In this study, we have sequenced partial genomes from seven additional APSE strains to address the evolution and extent of toxin variation in this phage. The APSE lysis region has been a hot spot for nonhomologous recombination of novel virulence cassettes. We identified four new toxins from three protein families, Shiga-like toxin, cytolethal distending toxin, and YD-repeat toxins. These recombination events have also resulted in reassortment of the downstream lysozyme and holin genes. Analysis of the conserved APSE genes flanking the variable toxin cassettes reveals a close phylogenetic association with phage sequences from two other facultative endosymbionts of insects. Thus, phage may act as a conduit for ongoing gene exchange among heritable endosymbionts.     

Pseudomonas aeruginosa Population Structure Revisited

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.