Previous Article | Next Article ![]()
Journal of Bacteriology, June 2007, p. 4161-4167, Vol. 189, No. 11
0021-9193/07/$08.00+0 doi:10.1128/JB.00106-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
,
Skin Research Centre, Institute of Molecular & Cellular Biology, Garstang Building, Faculty of Biological Sciences, University of Leeds, Leeds LS2 9JT, United Kingdom,1 Sarum Biosciences Ltd., 8 Centre One, Lysander Way, Old Sarum Park, Salisbury SP4 6BU, United Kingdom2
Received 19 January 2007/ Accepted 22 March 2007
|
|
|---|
|
|
|---|
Bacteriophages that infect P. acnes can be readily isolated from human skin. The study of these bacteriophages has until now been limited to the development of phage typing systems to distinguish the different serotypes of P. acnes (19, 33). Bacteriophages are the most abundant organisms on Earth, and they play an important role in bacterial diversity and pathogenesis. Greater study of P. acnes bacteriophages should make a significant contribution to the study of propionibacterial genetics and to the wider field of phage biology and evolution.
With the sequencing of bacteriophages comes the potential to genetically manipulate the host bacterium. Although the physiology of P. acnes has been well studied in vitro (4, 13, 14), genetic studies have been limited by the lack of genetic transfer systems. The inability to manipulate P. acnes genetically has severely hindered progress toward understanding the role of this organism in acne and other diseases. Similar problems existed several years ago with the study of pathogenic mycobacteria, which are closely related to propionibacteria taxonomically. The development of transformation systems and plasmid vectors for mycobacteria has allowed greater study of the genetics of these organisms (18, 27). In addition, several mycobacteriophages have been characterized in detail, and the genomes of several such phages have been sequenced and analyzed (10, 11, 15, 24, 25). These phages have been proposed to be useful tools for the genetic manipulation of mycobacteria and for understanding mycobacterial gene expression. Studies of them have also made a significant contribution to the study of phage evolution and have demonstrated that there is a high degree of mosaicism among the mycobacteriophages.
Further studies of P. acnes bacteriophages may also lead to the development of phage therapy for acne. Such therapy may overcome the current problems surrounding long-term use of antibiotics and the emergence of resistance in P. acnes and other skin commensals. We present here the first genome sequence of a P. acnes bacteriophage along with an extensive sequence analysis and characterization of the phage morphology and host range.
|
|
|---|
Phage purification and DNA extraction. Particles of phage PA6 were prepared using a standard plate lysate procedure (28). An aliquot of a stock of PA6 sufficient to produce nearly confluent bacterial lysis was added to 100 µl of late-log-phase P. acnes AT1 resuspended to an optical density at 600 nm of 2.5 in SM buffer (20 mM Tris-HCl [pH 7.5], 100 mM NaCl, 60 mM MgSO4, 0.01% [wt/vol] gelatin). Phage were allowed to adsorb for 15 min at 34°C before the mixture was added to 3 ml molten top agarose (0.7% [wt/vol] agarose type VII; Sigma) and poured onto a 90-mm reinforced clostridial agar plate. Ten plates were prepared in this way and incubated anaerobically overnight at 34°C. Phage were harvested by overlaying each plate with 3 ml SM buffer and incubating it for 4 h at 4°C. Buffer was removed from the plates and centrifuged to remove bacteria (8,000 x g, 10 min, 4°C). Phage DNA was subsequently extracted using the Wizard Lambda Preps DNA purification system (Promega) according to the manufacturer's instructions.
DNA sequencing assembly and analysis. Sequencing of the phage genome was carried out using a shotgun strategy by Lark Technologies Inc. Fragmented phage DNA was cloned into the pUC19 vector, and clones were chosen at random for sequencing. Dye terminator sequencing reactions were performed with ABI3100 and ABI3700 automatic sequencers (Applied Biosystems). DNA sequencing data were assembled and edited using the Phred-Phrap software package (CodonCode Corp.). Single-stranded extensions not identified after shotgun sequencing were isolated by self-ligation of the PA6 genome, PCR over the ligated region, and sequencing of the PCR product. Open reading frames (ORFs) were identified using GeneMark (1) and FramePlot (17). Assignment of an ORF was based on the presence of a possible initiation codon (AUG, GUG, or UUG) preceded by a potential ribosome binding site and the third-letter G+C content of codons within the ORF. Codon usage was determined using the Countcodon program (http://www.kazusa.or.jp/codon/countcodon.html), and protein transmembrane prediction was carried out using the TMHMM server (http://www.cbs.dtu.dk/services/TMHMM-2.0/).
Electron microscopy. Purified phage from a plate lysate were placed onto carbon-coated copper grids and negatively stained with 1% (wt/vol) uranyl acetate. Particles were visualized by using a Philips CM10 transmission electron microscope.
Nucleotide sequence accession number. The GenBank accession number for the bacteriophage PA6 genome is DQ431235.
|
|
|---|
![]() View larger version (89K): [in a new window] |
FIG. 1. Transmission electron micrograph of negatively stained PA6 particles. Bar = 100 nm.
|
DNA sequence and organization of the PA6 genome. The linear, double-stranded DNA genome of phage PA6 was sequenced using a shotgun strategy. A single contig consisting of 29,739 bp was assembled from 681 sequences and was identified as the full-length phage genome, and >9-fold coverage was achieved. An additional 13-base single-stranded extension (CCTCGTACGGCTT) was identified at each end of the PA6 genome. The overall G+C content of the PA6 genome was 54.0%, which is slightly lower than the G+C content determined for the genome of P. acnes (60%) (3). A total of 48 ORFs were identified in the genome of PA6. Seventeen ORFs possessed potential GUG initiation codons, and the remaining ORFs had AUG initiation codons. All three of the standard termination codons were used, and 29 ORFs had a UGA termination codon. No tRNA genes were found in the genome. The coordinates of the protein-encoding genes identified are shown in Table 1, and a diagram of the PA6 genome is presented in Fig. 2.
|
View this table: [in a new window] |
TABLE 1. Coordinates of PA6 genes and molecular masses, putative functions, and homologues of predicted proteins
|
![]() View larger version (20K): [in a new window] |
FIG. 2. Genome organization of PA6. The 29,739-bp genome is represented by a horizontal bar with 1-kb intervals marked. Putative ORFs are indicated by boxes at three different heights representing the three frames on each strand. Rightward-transcribed genes are above the genome, and leftward-transcribed genes are below the genome. Predicted functions and homologues are indicated. ORFs with phage and nonphage homologues are indicated by gray and black boxes, respectively. ORFs unique to PA6 are indicated by white boxes.
|
Genes encoding structural and assembly proteins. The left arm of the PA6 genome contains genes encoding proteins of the viral particle. The organization of these genes is similar to that of many of the mycobacteriophages, and as expected, the morphology of the PA6 virion is similar to that observed for these phages (25). Many of the predicted gene products are related to those of mycobacteriophages. Predicted functions have been assigned based on homology to other phage proteins having known or predicted functions and on matches to the COG (http://www.ncbi.nlm.nih.gov/COG/) and Pfam (http://www.sanger.ac.uk/Software/Pfam/) databases.
PA6 gp2 is a putative terminase large subunit and shows a high degree of homology to other members of this protein family (43% similarity to TM4 gp4 and 47% similarity to the lactococcal phage r1t). The smaller gp1 protein shows weak homology to gp32 of phage VWB that is predicted to be the terminase small subunit (32). This homology, along with the size and position of gene 1, makes the gp1 protein the likely terminase small subunit of PA6. This family of proteins is involved in the packaging of phage DNA into the head of the virus particle (12). This process also involves the portal protein (gp3 of PA6) that lies at the head-tail interface and enables passage of genomic DNA into phage heads during packaging and ejection during infection. The gp5 and gp6 proteins are most likely components of the phage head. gp5 is likely to be the scaffold protein that directs assembly of the head proteins, and gp6 is likely to be the major head protein that shows homology to the gp7 head protein of mycobacteriophage Che9d (55% similarity over 312 residues). Three proteins were identified that are proposed components of the phage tail. gp11 exhibits homology with the major tail protein of lactococcal phage
LC3. gp14 is a putative tape measure protein that in bacteriophage
determines tail length (20, 21). The PA6 tape measure protein does not contain any of the three sequence motifs predicted to have a role in cell wall catabolism, one of which shows homology to the resuscitating-promoting factor family of proteins thought to be involved in stimulating dormant cell growth (22, 26). gp15 is proposed to be a minor tail protein based on homology to the proven tail protein (gp27) of phage L5 (15). One gene of note in this region is ORF16. The predicted product of this gene is a 42.8-kDa protein. The first 80 amino acids of the N terminus of this protein show a low level of homology to a family of protease enzymes that specifically bind to and degrade proteins with nonpolar C termini. It is possible that this protein is involved in some proteolytic process during phage assembly. Alternatively, this protein may be involved in entry of the phage into the host cell through degradation of peptide linkages in the cell wall to allow injection of viral DNA. The location of the gene encoding gp16 suggests that it may be a component of the PA6 tail, making this function not unlikely.
Lysis functions. ORF20, located near the end of the section of rightward-transcribed genes, encodes a putative N-acetylmuramoyl-L-alanine amidase. This enzyme is most likely involved in cell lysis following phage assembly rather than in phage entry as the gene immediately downstream of ORF20 (ORF21) encodes a putative holin. Following phage assembly, holin proteins assemble to form pores in the cellular membrane, allowing the amidase enzyme access to the surrounding peptidoglycan (35). PA6 gp21 appears to belong to the class II group of holins. Members of this group are usually 65 to 95 residues long and possess two transmembrane domains. Although gp21 is larger (113 residues), it is predicted to contain two transmembrane domains instead of the three transmembrane domains of the class I holins. The PA6 amidase exhibits homology with amidases from the PG1 and Che8 mycobacteriophages. However, the highest degrees of similarity are with an N-acetylmuramoyl-L-alanine amidase and likely an autolysin from its host, P. acnes (67% identity for amino acids 2 to 145). This is perhaps not unexpected considering their shared substrate.
DNA replication functions. The leftward-transcribed section of the genome contains genes predicted to encode proteins involved in phage DNA replication. PA6 gp31 and gp32 are both putative DNA primases. PA6 gp31 exhibits homology with one of the DNA primases of Bxb1, gp50 (35% identity over 154 residues). The organization of these genes is similar to that of the DNA primase-encoding genes of mycobacteriophages Bxb1, L5, and D29, with an extensive overlap of genes 31 and 32. It has been suggested that in these mycobacteriophages, the primase proteins are produced through a programmed translational frameshift resulting in a single functional protein (24). PA6 gp34 is proposed to be a DNA helicase based on its match to the COG DNA helicase group. This protein contains an ATP-binding domain and exhibits homology to Bxb1 gp57 and L5 gp65, which are also predicted DNA helicases. Other PA6 genes in this region encode proteins with no predicted functions but with homology to proteins of other phages (Table 1).
One significant aspect of the PA6 genome is the absence of any gene encoding an obvious integrase-type enzyme. Plaques of PA6 on lawns of P. acnes are clear with hazy centers and not turbid, as is characteristic of truly lysogenic phages. This plaque morphology has been observed for TM4, which is not temperate but may be capable of forming pseudolysogens in some strains of mycobacteria (11). Also absent from the PA6 genome is any gene encoding a repressor-like protein that would indicate a lysogenic function. Although the presence of such a protein cannot be completely ruled out, the lack of these two lysogeny-related proteins strongly suggests that PA6 is a lytic phage and does not integrate into the host genome. Other PA6 proteins possibly involved in DNA regulation/replication include gp37 and gp48. These proteins are predicted to be an exonuclease and an endonuclease, respectively, due to matches to the COG RecB exonuclease group and Pfam HNH endonuclease family, respectively.
Other genes and features of interest.
Two other predicted ORFs of PA6 are worthy of mention. ORF19 encodes a protein with considerable homology to the human collagen
1(III) chain precursor. In particular, gp19 contains many G-X-Y repeats that are characteristic of many collagens (Fig. 3). The gene encoding gp19 is located in the vicinity of the genes encoding the structural proteins that constitute the phage tail. The presence of collagen-like proteins in phage tails has been noted previously (29). The function of such proteins is unknown, but it has been suggested that they stabilize the tail structure. It is tempting to speculate that in PA6 this protein may have evolved a second function, binding to the host cell and, more specifically, to potential collagen-binding surface proteins of P. acnes. The P. acnes genome contains several genes encoding putative surface adhesins that are likely to be important for the colonization of human skin, one or more of which may bind collagen.
![]() View larger version (17K): [in a new window] |
FIG. 3. Comparison of PA6 gp19 with the human collagen 1(III) chain. Protein sequences are aligned, and identical amino acid residues are marked with asterisks. The glycine residues of the repetitive G-X-Y motifs are in bold type. The numbers indicate the positions of the amino acid sequences in the proteins.
|
Within the 1.4-kb noncoding region of the PA6 genome there is an approximately 500-bp section that contains numerous direct, inverted, and symmetric repeats (Fig. 4). This 500-bp section has a particularly low G+C content (35%). Such repeat sequences may be binding sites for transcription factors and/or may form secondary structures important for the regulation of gene transcription and DNA replication. The gene organization and direction of transcription indicate that this region is likely to contain promoters for phage gene transcription.
![]() View larger version (26K): [in a new window] |
FIG. 4. Repeat motifs located in the noncoding region of the PA6 genome. The noncoding region between ORF45 and ORF46 contains several direct (DR), inverted (IR), and symmetric (SR) repeats that are indicated by arrows. The numbers indicate the positions in the PA6 genome, and the start of ORF46 is in bold type.
|
|
|
|---|
The host range of PA6 is restricted to P. acnes. Although it infects and lyses all strains of P. acnes tested to date, PA6 cannot infect other members of the genus Propionibacterium. This restricted range is likely to be dictated by the phage receptor and/or other components of the phage tail involved in entry of the phage into the cell. The PA6 minor tail subunit gp15 may be involved in host specificity. This protein has homology to D29 gp27, L5 gp27, and Bxb1 gp23, which have been proposed to be host specificity determinants (24). Isolation of phages capable of infecting the closely related cutaneous propionibacteria P. granulosum and P. avidum would allow comparisons of these genes to be made. However, phages with the ability to infect these species have not been isolated to date.
The sequencing of PA6 makes a valuable contribution to the study of bacteriophage evolution. The close relationship between mycobacteria and propionibacteria is reflected in the relationships between their respective phages, and PA6 is more closely related to mycobacteriophages than to phages of any other bacterial group. More significantly, PA6 should contribute to the study of propionibacterial genetics. Finally, the sequencing of PA6 should greatly enhance the development of a potential bacteriophage therapy to treat acne and therefore overcome the significant problems associated with long-term antibiotic therapy and bacterial resistance.
Published ahead of print on 30 March 2007. ![]()
Supplemental material for this article may be found at http://jb.asm.org/. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»