Research article | Open | Published:
Complete genomic sequences of Propionibacterium freudenreichii phages from Swiss cheese reveal greater diversity than Cutibacterium (formerly Propionibacterium) acnes phages
BMC Microbiologyvolume 18, Article number: 19 (2018)
A remarkable exception to the large genetic diversity often observed for bacteriophages infecting a specific bacterial host was found for the Cutibacterium acnes (formerly Propionibacterium acnes) phages, which are highly homogeneous. Phages infecting the related species, which is also a member of the Propionibacteriaceae family, Propionibacterium freudenreichii, a bacterium used in production of Swiss-type cheeses, have also been described and are common contaminants of the cheese manufacturing process. However, little is known about their genetic composition and diversity.
We obtained seven independently isolated bacteriophages that infect P. freudenreichii from Swiss-type cheese samples, and determined their complete genome sequences. These data revealed that all seven phage isolates are of similar genomic length and GC% content, but their genomes are highly diverse, including genes encoding the capsid, tape measure, and tail proteins. In contrast to C. acnes phages, all P. freudenreichii phage genomes encode a putative integrase protein, suggesting they are capable of lysogenic growth. This is supported by the finding of related prophages in some P. freudenreichii strains. The seven phages could further be distinguished as belonging to two distinct genomic types, or ‘clusters’, based on nucleotide sequences, and host range analyses conducted on a collection of P. freudenreichii strains show a higher degree of host specificity than is observed for the C. acnes phages.
Overall, our data demonstrate P. freudenreichii bacteriophages are distinct from C. acnes phages, as evidenced by their higher genetic diversity, potential for lysogenic growth, and more restricted host ranges. This suggests substantial differences in the evolution of these related species from the Propionibacteriaceae family and their phages, which is potentially related to their distinct environmental niches.
Most bacteriophage populations display a wide range of genetic diversity, including those phages infecting Mycobacterium, Staphylococcus, and Pseudomonas spp. [1,2,3]. In contrast, Cutibacterium acnes bacteriophages, whose host was formerly a member of the genus Propionibacterium and was recently reclassified , have been found to have a surprisingly limited genetic diversity . Propionibacteria and Cutibacteria are Gram-positive, anaerobic to aerotolerant, rod-shaped bacteria of the Actinobacteria class, originally named for their unique metabolism, specifically the ability to synthesize propionic acid via the activity of transcarboxylase enzymes during fermentation . C. acnes, along with C. granulosum and C. avidum, were formerly referred to as the cutaneous propionibacteria and are all members of the human skin microbiome. The classical type of propionibacteria, including P. freudenreichii, P. acidifaciens, P. cyclohexanicum, and P. australiense, have been found in milk and cheese, as well as in various other biotopes, such as silage, soil, pasture, and the human mouth [7, 8]. C. acnes and P. freudenreichii are not highly related at the nucleotide level, but show a high synteny with regards to their protein coding sequences . C. acnes also has a lower average percent GC (60% vs. 67%) and encodes a larger proportion of genes involved in processes such as host association, tissue degradation, and iron acquisition [4, 9].
P. freudenreichii is specifically utilized in the manufacture of Swiss-type cheeses, such as Emmental, for ripening the culture. Through the Wood-Werkman cycle, this bacterium drives fermentation of lactate into acetate, propionate (which gives Swiss-type cheese its unique nutty and sweet flavor), and carbon dioxide (responsible for the holes or “eyes” in the cheese) (3 C3H6O3 → 2 C2H5CO2 + C2H3O2 + CO2) . An estimated one billion living cells of P. freudenreichii are present in one gram of Emmental cheese . In contrast to most lactic acid bacteria, this bacterium mainly breaks down lipids, forming free fatty acids. Recent research has focused on the possible probiotic benefits incurred from consuming P. freudenreichii. Notably, these bacteria have anti-inflammatory and probiotic properties for the gastrointestinal tract in animals  and have also been found to display cytotoxic activity against colon cancer cell lines in vitro [13, 14].
Bacteriophages that infect P. freudenreichii have been isolated from Swiss-type cheese [15, 16] and are common contaminants of the cheese manufacturing process. The genome of a single filamentous variant of P. freudenreichii phage has been described , the first filamentous phages of any Gram-positive bacterium to be reported. Most P. freudenreichii phages, however, are of the Siphoviridae type [15, 18], similar to the C. acnes phages. Recently, the complete genomes of two dsDNA tailed P. freudenreichii phages were reported . These are distinct from the C. acnes phages but highly similar to one another, raising the question as to whether this group of phages also displays restricted genomic diversity. Here we describe and compare the genome sequences of seven P. freudenreichii phages isolated from Swiss-type cheeses and a related prophage in P. freudenreichii CIRmBIA139.
Bacterial strains (Additional file 1: Table S1) were from either the American Type Culture Collection (ATCC) or from the Collection du Laboratoire de Recherche de Technologie Laitière, Institut National de la Recherche Agronomique, Rennes, France. Colonies were grown on Brucella agar with 5% sheep’s blood, hemin, and vitamin K (Thermo Fisher Scientific, Remel Products, Lenexa, KS, USA) at 30 °C for 5-7 days under anaerobic conditions in a sealed anaerobic box with oxygen-absorbing, carbon dioxide-generating AnaeroPack-Anaero® sachets (Mitsubishi Gas Chemical Co., Inc. [MGC], Tokyo, Japan). Liquid cultures inoculated from single colonies were grown anaerobically at 30 °C to mid-log phase (OD600 = 0.4-0.7) in yeast extract sodium lactate (YEL) media, containing (per Liter), 10 g pancreatic digest of casein, 10 g yeast extract, 10 g sodium lactate, 2.5 g KH2PO4, 5 mg MnSO4 .
Isolation of bacteriophage from cheese
The P. freudenreichii phages, B22, E1, E6, G4, and B3 were previously isolated from Swiss-type cheese as reported [15, 16]. Two new phages, Anatole and Doucette , were isolated at UCLA using a similar protocol. Ten grams of each cheese sample was shredded and dispersed with a Cuisinart Mini-Prep® Plus Processor in 90 ml of YEL media for 5 min. The suspension of cheese + YEL was enriched for 5 days at 30 °C in a sealed box with an AnaeroPack-Anaero® sachet (MGC). The suspension was then centrifuged at 3000Xg for 20 min and filtered through a Corning® 0.2 μm syringe filter (Corning Inc., Corning, NY, USA).
Mid-log phase P. freudenreichii TL110 (200 μl at OD600 = 0.4-0.7) was incubated with 100 μl filtered supernatant from cheese enrichments for 20 min at 30 °C. This was then plated as a top agar overly in 3 ml molten YEL top agar (0.7% agar) on YEL agar plates (1.5% agar). Plates were incubated for 5 days at 30 °C in a sealed box with an AnaeroPack-Anaero® (MGC). Plates were observed each subsequent day after plating, and plaques were observed on day 5.
Isolated plaques were picked from top agar lawns using a sterile micropipette tip into 100 μl YEL media and stored at 4 °C. These were serially diluted, and 10 μl of each dilution was incubated with 200 μl of TL110 strain at 30 °C for 20 min, and plated using the soft agar method, described above. In total, three rounds of plaque purification were performed for each phage. Lysates were then collected from plates showing a webbed lysis pattern indicative of multiple rounds of phage infection and filtered through a 0.2-μm syringe filter (Corning, Inc.).
Phage electron microscopy
High titer lysates at 109-1011 plaque-forming units (PFU)/ml were subjected to ultracentrifugation at 20,000Xg for 2 h, and the resulting pellets were resuspended in 1 ml phage buffer (10 mM Tris-HCl, pH 7.5; 10 mM MgSO4; 68.5 mM NaCl; 1 mM CaCl2). Purified phage lysates were spotted in 5 μl aliquots onto freshly glow-discharged 400-mesh carbon-Formvar-coated copper grids and allowed to sit for approximately 1 min. Grids were washed with distilled water and stained using 1% (w/v) uranyl acetate. Images were captured using an FEI Morgagni transmission electron microscope at 56,000X magnification.
Phage genome sequencing and analysis
High titer lysates (109-1011 PFU/ml) were used for phage DNA extraction. Briefly, DNase I (1 μg/ml final concentration) and RNase I (1 μg/ml final concentration) were added to the lysates for 30 min at room temperature to degrade any host nucleic acids in the lysate. Phage DNA was then isolated using the Promega Wizard DNA Clean Up Kit (Promega Corp., Madison, WI, USA), as described previously . Purified genomic DNA was subjected to whole-genome sequencing using the Illumina MiSeq platform by the Technology Center for Genomics and Bioinformatics (TCGB) at UCLA. The approximate coverage obtained for each phage is as follows: Anatole, 11,161; B3, 9371; B22, 9019; Doucette, 3409; E1, 1000; E6, 14,437; and G4, 13,283. Raw Sequence reads were assembled using Newbler Assembler and Consed or the ABySS software package with manual correction. Genome ends were deduced from large read-start buildups located on opposite strands, separated by a 12-base region of lower coverage. Finished sequences were analyzed and annotated using the genome editor DNAMaster (http://cobamide2.bio.pitt.edu). GenBank accession numbers are listed in Table 1. Genome comparisons were performed and represented using Phamerator and the database ‘Actinobacteriophage_685’ .
Isolation of bacteriophages
To investigate the diversity of bacteriophages that infect P. freudenreichii, we isolated P. freudenreichii phages from Swiss-type cheese by culturing with, or without, indicator strains of host bacteria (Table 1). A total of seven independent bacteriophages isolated from different cheese samples were included in our analyses (Table 1). Phages B22, E1, E6, and G4 were isolated previously from French Emmental-type cheese in 1992, whereas Phage B3 was isolated from French Comté-type cheese in 1991 [15, 16]. These phages, with the exception of G4, were regrown in cultures containing P. freudenreichii strain TL110 as the host; phage G4 was grown with TL19 as the host. Phages E6 and G4 were originally propagated on strains, LS2500 and Z435, respectively. The phages named ‘Anatole’ and ‘Doucette’ were isolated in 2015 by processing Emmental-type cheese of French factory origin in YEL media. Enriched cheese filtrates were amplified by infecting a sensitive P. freudenreichii strain (TL110, Table 1), and two rounds of phage purification were performed prior to generation of high titer lysates.
Morphology of P. freudenreichii plaques and phage virions
Isolated phages infecting P. freudenreichii strains (Table 1) were plated onto YEL agar plates, with plaques appearing after culturing for 5 days at 30 °C under anaerobic conditions. All plaques appear relatively clear and small, with similar morphologies and sizes. P. freudenreichii strain TL110 was used as the host to prepare high-titer phage lysates with titers 109-1011 PFU/ml. These were analyzed by electron microscopy, and all phages were found to have a siphoviral morphology with isometric heads of ~ 50 nm in diameter and long flexible tails (Fig. 1); the tails of Anatole and B22 are ~ 160 nm and ~ 190 nm in length, respectively (Fig. 1). The others are predicted to have similarly sized tails, based on the length of their tape measure genes, the product of which determines phage tail length . The predicted tape measure protein lengths are as follows: 1132 amino acids (aa) in Anatole, E1, and B3; 1326 aa in B22; and 1364 aa in Doucette, G4, and E6. Morphologically, the P. freudenreichii phages resemble C. acnes phages, which have a slightly shorter tail length of ~ 150 nm .
Host range analysis
To determine the host range specificities of these P. freudenreichii phage isolates, we infected different P. freudenreichii strains, including TL110, TL18, TL19, and TL29, with phages E6, B3, E1, G4, B22, Doucette, and Anatole. As shown in Table 2, TL18 is resistant to all seven phage isolates, whereas TL110 and TL29 are sensitive to all phages except G4. Of note, G4 is the only phage isolate that can specifically infect TL19, indicating that this phage has a host range that is more selective, and distinct from, than the other P. freudenreichii phage isolates. In addition, we also included seven P. freudenreichii strains from ATCC in the host range analysis, including ATCC 39393, 9616, 9615, 9614, and 13,273 (subsp. shermanii) and ATCC 9617 and 6207 (subsp. freudenreichii). Our results indicate that these strains are resistant to infection by all phages tested. Collectively, these host range data demonstrate that the P. freudenreichii phages show a higher degree of host specificity than the phages infecting C. acnes .
Genomic DNA characterization of P. freudenreichii phages
Total genomic DNA was prepared from high titer P. freudenreichii phage lysates, and the complete genome sequences were determined. We observed that the genome sizes among these sequenced phages are similar and vary from 35,209 bp to 38,555 bp, which is larger than those of the C. acnes phages that average 29,432 bp in length. However, the overall genome sizes for the P. freudenreichii phages and C. acnes phages are only about half the average size of the mycobacteriophages, whose hosts are related members of the high GC Gram-positive Actinobacteria . All P. freudenreichii viral genomes have defined ends with 12 bp 3′ single-stranded extensions, although the sequences vary among the phages: 5’-CAAACAAGTCAT in Doucette, 5’-CCAACAAGTCAT in B22, E6 and G4, and 5’-CTCCCGCTCGAA in Anatole, B3, and E1.
The percent GC content was also observed to vary little among the P. freudenreichii phages, ranging from 64.40% to 65.70% GC. The average percent GC content of P. freudenreichii phages is 64.94%, similar to that of their P. freudenreichii bacterial hosts (67.3%) and other Propionibacterium species, which range from 53%-68% . This is in contrast to the C. acnes phages, which show a more restricted range of GC% content (53.76%-54.37%) with an average of 54.04%, considerably lower than the ~ 60% GC content of its host C. acnes . The low variation in the GC% contents of P. freudenreichii phages suggests that they may have shared phylogenetically similar hosts in their recent evolutionary histories .
Nucleotide sequence identity among P. freudenreichii phages
Fully sequenced genomes for the seven P. freudenreichii phages sequenced in this study, as well as the previously sequenced P. freudenreichii phage, PFR1 , and the C. acnes phage, P100.1 , were compared according to their shared genome contents and average nucleotide sequence similarities (Fig. 2, Additional file 1: Table S2). The phages all share some nucleotide sequence similarity, although there is more extensive similarity among phages Doucette, B22, E6, and G4, compared to Anatole, E1, and B3, which also share extensive similarity with each other. For example, Doucette and B22 share 97% identity, with a coverage of 76%; whereas Doucette and Anatole share 91% identity, but with a coverage of only 24%. Conversely, Anatole and B3 share 98% identity with a coverage of 91%. On this basis we propose organizing the phages into two clusters based on previously described criteria, with span-length match being the primary factor that distinguishes the two clusters (as described further below) and cluster names deriving from the Actinophage Database [28, 29]: Cluster BW contains Doucette, B22, E6, and G4, and Cluster BV contains Anatole, E1, and B3 (Fig. 2). We further note that the previously characterized P. freudenreichii phage, PFR1, is more distantly related to the phages sequenced in this study and groups in a separate cluster designated BX. From the dotplot, we can observe discrete regions of nucleotide identity between PFR1 and the Cluster BW phages (Doucette, B22, E6, and G4) and more limited identity between PFR1 and the Cluster BV phages. No observable identity similarity can be detected between C. acnes phage P100.1 and the group of P. freudenreichii phages.
The members of Cluster BW share ~ 90% average nucleotide identity (Additional file 1: Table S2), but as illustrated by comparison of B22 and Doucette, there are regions of near-identity interspersed with regions that are more diverse (Fig. 2). The left end of E6 is very similar to that of Doucette, but the right ends are more diverse (Fig. 2). The Cluster BV genomes share > 96.4% average nucleotide identity to one another (Additional file 1: Table S2). Of note, phages Anatole and E1 have ~ 99% identity at the nucleotide level and differ by only two single nucleotides, despite the fact that E1 was isolated in France in 1995  and Anatole was isolated in the United States in 2014. B3 has 96.5% identity to E1 in overall genome sequence, although there are several gaps or regions of difference between the two genomes, particularly in the right genome arms, as can be seen in the dotplot in Fig. 2. B3 is also slightly longer than E1, with genome lengths of 35,948 bp and 35,209 bp, respectively. The nucleotide sequence identities between Clusters BW and Cluster BV range from 82% to 96%, but the matching regions only span 31-45% of the genome lengths (aligned using BLASTN), mostly in the 3′ halves of the genomes (Fig. 2). These parameters justify the separation of the phages into the separate BW and BV clusters. The annotated genomes of Anatole and Doucette are shown in Fig. 3, and genome annotations of B22, E6, G4, E1, and B3 are shown in Additional file 2: Figure S1, Additional file 3: Figure S2, Additional file 4: Figure S3, Additional file 5: Figure S4, and Additional file 6: Figure S5, respectively.
Genome annotations of Cluster BW P. freudenreichii phages
Annotation of the open reading frames in the Doucette genome, a 37,429 bp representative from Cluster BW, reveals 61 predicted protein-coding genes (Fig. 3, Table 1). The overall genome organizations of Doucette and other members of Cluster BW (B22, E6, and G4; Additional file 2: Figure S1, Additional file 3: Figure S2, Additional file 4: Figure S3) are similar to each other, although there are notable differences, particularly in the right arms of the genomes (Fig. 4). The virion structure and assembly genes are organized in the left arms of the genomes. We note that the portal and major capsid proteins are related to those in mycobacteriophage Gaia (62% and 37% aa identity, respectively), and the major tail subunits have C-terminal Ig-like domains, as described for other phage structural proteins . The Doucette and E6 tail proteins (gp21 and gp22, respectively) differ from B22 and E6 (gp21) in that they contain a glycosyl hydrolase family 43 domain that is typically an endolysin component. These are likely to be tail proteins, however, as the lysis cassette is readily identifiable and is located further to the right in the genomes. Although this glycosyl hydrolase domain is atypical of phage tail proteins, we note that several mycobacteriophage tail proteins contain putative D-ala-D-ala carboxypeptidase domains . We therefore predict that the cell wall target of these enzyme activities may be cleaved or remodeled during the process of phage adsorption and DNA injection. The lysis cassette also varies among the Cluster BV phages, with G4 having a different endolysin and holin gene than is found in Doucette, B22, and E6 (Fig. 3, Additional file 2: Figure S1, Additional file 3: Figure S2, Additional file 4: Figure S3).
The Cluster BW phages have several genes near the center of the genomes that are typically associated with temperate growth (e.g., putative integrase and repressor). Doucette gp31 is a tyrosine integrase, which is more distantly related to the integrases encoded by the other Cluster BW and BV phages than these proteins are to one another. That is, the predicted integrase proteins from B22 (gp29), E6 (gp28), and G4 (gp29), as well as those from Anatole (gp28), E1 (gp28), and B3 (gp28) share ~ 96-99% aa identity, whereas they have only ~ 54-55% sequence identity to Doucette gp31. Comparison of the sequences upstream of Doucette gene 31 with P. freudenreichii genomes identified a 45 bp region of sequence identity corresponding to the common core attP (Doucette coordinates 22,180 – 22,224) and an attB site (coordinates 934,527 – 934,571 in P. freudenreichii subsp. shermanii CIRM-BIA1) overlapping a tRNApro gene. We predict that Doucette integrates at this attB site in lysogenic establishment. The other Cluster BW genomes differ in their sequences upstream of the integrase gene, and a similar analysis identified a 49 bp common core corresponding to the attP common core in B22 (coordinates 20,750 – 20,798) and an attB in the P. freudenreichii genome (subsp. shermanii CIRM-BIA1, coordinates 376,562 – 376,610), such that the attB site overlaps a tRNAgly gene. The three Cluster BV phages have closely related attP sites and integrase genes and are predicted to integrate at the same attB site.
Three of the Cluster BW phages (Doucette, B22, and G4) encode putative repressors (genes 33, 31, and 31, respectively) that are identical to each other, and also identical to the repressor (gp30) of Cluster BV phage B3. We therefore predict that these form a homoimmune group of temperate phages. The E6 repressor (gp30) differs from these (40% aa identity), but is identical to the repressors of Cluster BV phages Anatole and E1, likely forming a separate homoimmune group. In each of the genomes, the repressor is transcribed leftwards and divergently from an operon of rightwards-transcribed genes, in which the first is a putative cro-like gene. There is considerable variation in the Cluster BW gene content in the Integrase/Repressor regions. G4 contains three leftwards-transcribed genes (26, 27, 28) of unknown function (Additional file 4: Figure S3) to the left of the integrase, and there are three small genes common to Doucette (28, 29, 30) and B22 (26, 27, 28) in this region that are absent from the other phages; the function of these is also unknown. E6 is unusual in that it contains a second set of divergently transcribed genes (34 and 35) organized as an additional repressor-cro pair of genes, suggesting that it may have complex immunity properties. It also has an insertion of a large gene (33) that contains a putative DNA binding motif but is otherwise of unknown function.
The Cluster BW right arm genes also vary among the genomes, although they all code for a RecET-like system, a WhiB-like regulator, a single-stranded DNA binding protein, and an HNH endonuclease (Fig. 3, Additional file 2: Figure S1, Additional file 3: Figure S2, Additional file 4: Figure S3); E6 contains a transposase gene (gp55) of the IS1341 group of insertion sequences  (Additional file 3: Figure S2). Doucette and B22 both encode DNA methylase genes, although they are not related to each other and are at different genomic loci. We note that, just downstream of the putative antirepressor genes (e.g. Doucette 35), Doucette, B22, and E6 contain two genes that are tandemly repeated (e.g. Doucette 38 and 39, B22 37 and 38, and E6 39 and 40); this gene is also duplicated in G4, but the two copies, genes 37 and 55, are separated by ~ 5.5 kb. However, the intergenomic similarities are closer than the intragenic similarities (e.g. Doucette gp38 and gp39 share 55% aa identity, and Doucette gp38 and B22 gp39 share 93% identity; Additional file 7: Figure S6).
Genome annotations of Cluster BV P. freudenreichii phages
Annotation of the Anatole genome, a 35,284 bp representative of Cluster BV, reveals 55 predicted protein-coding genes, and no tRNA genes (Fig. 3, Table 1). The genome organizations of Anatole (Fig. 3) and E1 (Additional file 5: Figure S4) are near-identical in their annotations, but differ from B3 (Additional file 6: Figure S5) at several locations. The virion structure and assembly genes are organized in the genome left arms (Anatole genes 3 – 21) and arranged syntenically with structural genes of other phages with siphoviral morphologies. Several of the structural proteins (e.g. portal, prohead protease, capsid, major tail subunit) have sequence similarity to phages of other actinobacterial hosts, including Cluster N mycobacteriophages, corynebacteriophage BFK20, and microbacteriophage Min1, although amino acid sequence similarities are modest (27-45% aa identity). The minor tail proteins implicated in host recognition are not related to other phages other than Anatole gp21, which is related to the tail proteins of Cluster BW phages Doucette and E6. The lysis cassette – including endolysin and holin genes (Anatole 24 and 25, respectively) – are located to the right of the virion structural genes (Fig. 3).
Similar to the Cluster BW phages, the centers of the genomes include several genes associated with temperate growth. There are integration cassettes containing a tyrosine-integrase, with an attP located upstream, as described for Cluster BW phages B22, E6, and G4. The Anatole 49 bp attP common core is located at coordinates 21,035 – 21,083. To the right of the integration genes are two leftwards-transcribed genes, including the likely phage repressor (Anatole 30), which is divergently transcribed from a cro-like gene (31); immediately downstream are two genes (32, 33) related to anti-repressors (Fig. 3). As noted above, the Anatole and E1 repressors are identical to those of BW phage E6, whereas the B3 repressor is identical to the Doucette, B22, and G4 repressors. Phage B3 differs from Anatole and E1 in that it contains only a single antirepressor gene (Additional file 6: Figure S5).
The genes in the rightwards-transcribed right arm include a RecET-like recombination system (Anatole gp37, gp38), a Holliday Junction resolvase (Anatole gp44), a WhiB-like regulator (Anatole gp40), a single stranded DNA binding protein (Anatole gp47), and a transposase (Anatole gp50) of the IS1341 group of insertion sequences . As noted above, this transposase is also found in E6 (gp55). It is currently unclear whether this transposon may be active, as we have not been able to identify flanking inverted repeats flanking the transposon in any phage in which it is located. Phage B3 differs from Anatole and E1 in the right arm in two aspects. First, it contains a different recT-like gene (B3 gene 37; Additional file 6: Figure S5) that is more similar to Doucette 41 (63% aa identity) than to the Anatole and E1 recT genes. Second, genes 42 and 43 of Anatole and E1 are replaced with B3 genes 41 – 45, all of which have closely-related homologues in Doucette and four genes (with the exception of B3 42) in B22 (76 – 100% aa identity). Most are of unknown function, although G3 gp42 is a probably a DNA N6-methyltransferase.
Comparison between P. freudenreichii phages and other phages
P. freudenreichii phages were also compared with other phages in terms of their genome structures and sequences, especially those that infect the related host, C. acnes. Although these two groups of phages are not closely related at the sequence level, the overall arrangement of genes involved in packaging, head assembly, tail assembly, lysis, regulation, and replication is similar between P. freudenreichii phages and C. acnes phages; the C. acnes phages, however, lack genes for lysogenic growth. While most of the P. freudenreichii phage genes are transcribed in the rightwards direction, the non-structural genes are transcribed leftwards in the C. acnes phages . We noted previously that C. acnes phages display only very limited diversity, which is in marked contrast to the P. freudenreichii phages that not only show a greater diversity, but also considerable genome modularity, consistent with the mosaic architectures observed in other phage populations . In particular, there has been considerable swapping of genes between the Cluster BW and BV phages, illustrated notably by the distribution of immunity and integrase functions. This difference in degrees of diversity is evident in a heat map showing pairwise shared gene content between the C. acnes and P. freudenreichii phages (Fig. 5a). The phylogenetic relationship between the phages (Fig. 5b) is illustrated by a network analysis using SplitsTree , where the Cluster BV and Cluster BW phages are on distinct branches, and both are distantly related to other phages of the Actinobacteria phylum (Fig. 5b).
Similarities to P. freudenreichii prophages
Comparison of the P. freudenreichii phage genome sequences reveals several close matches to P. freudenreichii genomes corresponding to integrated prophages. The P. freudenreichii CIRM-BIA 139 genome includes two regions (contigs 1 and 17) matching the Cluster BW phages that together span the entire prophage from the attL left integration site to the rightmost integration site, attR (Fig. 6). Interestingly, the prophage is most closely related to the Cluster BW phages, but is mosaic, and contains parts of the individual phages but in a new combination. The prophage segment from attL to cos is most similar to B22, and most of the predicted gene products have amino acid sequences with close to 100% identity, including the DNA methylase. The putative repressor is 100% identical to those predicted in G4 and B22, and thus, the prophage likely shares immunity with the BW phages. The region from cos to the tape measure protein gene is similar in all four Cluster BW phages, but the genes from the tape measure protein gene to attR have more complex relationships, with CIRM-BIA 139 contig 1 00035 and 00030 gene products being 100% identical to B22 gp22 and gp23, respectively (Fig. 6), whereas the prophage genes 00025 - 00005 are similar to genes 24-28 of phage G4 (89 – 97% amino acid identity). P. freudenreichii strain CIRM-BIA 119 appears to carry an identical prophage to the one in BIA 139, and also spans two contigs. Searching with the Anatole genome did not identify any additional prophages.
The Propionibacteriaceae family is comprised of a diverse group of organisms, containing species such as C. acnes, a human commensal proposed to be involved in the pathogenesis of acne vulgaris, and P. freudenreichii, which is used in the manufacture of Swiss-type cheeses, as well as in the fermentative production of vitamin B12 and propionic acid. In both cases, the roles of phages targeting these bacterial hosts have not been well defined. For example, it is unclear how the resident C. acnes phage population shapes that of its host to promote health or disease. P. freudenreichii phages are often present during cheese ripening and can be isolated from mature cheeses; the semi-solid state of the cheese paste likely acts to limit phage diffusion, allowing phages to co-exist within large populations of sensitive hosts . Recent phage genome sequencing efforts by us and others have revealed a strikingly limited diversity amongst C. acnes phages, and phenotypic analyses indicate that these phages generally display broad host ranges, with most C. acnes strains being highly phage sensitive [5, 35, 36]. In contrast, despite the report of a single genome sequence for a ssDNA filamentous phage active on P. freudenreichii  and the recently published genomes of two highly similar dsDNA P. freudenreichii phages , little is known about the genetic composition of the phages infecting this organism and whether they share features in common with C. acnes phages. Although these two bacterial species have recently been re-classified into distinct genera , their shared metabolic features (e.g. generation of propionic acid) and degree of evolutionary relatedness, support the potential utility of comparing their respective phages.
Here, we report the full-length genome sequences for seven dsDNA tailed P. freudenreichii phage isolates. Sequence analysis indicates that these phages fall into two major clusters with more divergent nucleotide and amino acid sequences, as compared to the C. acnes phages. Further, P. freudenreichii phages display more restricted host ranges, with most P. freudenreichii strains showing high levels of phage resistance. Notably, although we obtained over 80 cheese samples from across the United States, Switzerland, and France, which were originally produced in the United States, France, Switzerland, Spain, and Italy, we were only able to isolate phages from Swiss-type cheeses produced in French factories. The reason for this is unclear, although some of this difficulty is likely due to low levels of phage contamination within individual samples, as well as the lack of phage sensitivity displayed by most P. freudenreichii strains. Because the ‘eyes’, or ‘holes’ in Swiss-type cheese arise due to the CO2 byproduct produced from the added P. freudenreichii strains, cheeses with fewer to no holes present might indicate higher likelihood of phage contamination and/or the use of bacterial strains that are more phage sensitive. Few studies have investigated the microbial composition of cheese, although, a recent analysis of the cheese-rind microbiome from over 130 samples of diverse origin has revealed group of 24 widely distributed bacterial and fungal genera . This group does not include Propionibacterium, although members of this genus were detected at low abundance in numerous samples. We therefore searched our phage genome sequences against a subset of the cheese-rind datasets, focusing on those from French cheeses, those with a natural rind (as is present on Swiss-type cheeses), and samples containing representatives from the Propionibacterium genus. From these analyses, we detected a small number of 100 bp reads with > 40% identity to our P. freudenreichii phages in five of the samples, however, coverage ranged from 0.3-4.2%, which was not sufficient to map these reads onto any one phage genome.
P. freudenreichii phages share limited genetic similarity with phages infecting other members of the Actinobacteria phylum, including the Cluster N mycobacteriophages, the corynebacteriophage BFK20, and microbacteriophage Min1. When the C. acnes phages are used as a model for comparison, we observe little to no relationship between these two groups of phages, although they display similar genome architectures. Further, P. freudenreichii phages and C. acnes phages share virtually no genes in common, and while C. acnes phages demonstrate a strikingly limited diversity, P. freudenreichii phages contain more diverse genomes. Of particular interest is the observation that whereas all known C. acnes phages lack putative lysogeny functions [5, 35, 38, 39], the P. freudenreichii phages characterized here all encode putative tyrosine integrases and repressor proteins, strongly suggesting they are temperate.
Our host range analyses also revealed critical differences between these two groups of phages. Whereas C. acnes phages generally show broad host ranges, infecting most laboratory and clinical C. acnes isolates, the host ranges of P. freudenreichii phages are more restricted, and many bacterial strains show high levels of phage resistance. This may result from the presence of prophages in some P. freudenreichii strains, which has been suggested experimentally [15, 18] and is further supported by the presence of prophages closely related to the Cluster BW phages in the fully sequenced strains, CIRM-BIA 139 and CIRM-BIA 119. In contrast, and consistent with their broad host ranges, C. acnes prophages do not resemble any of the sequenced C. acnes phages [5, 35, 38, 39]. We therefore speculate that superinfection immunity due to the presence of prophages in P. freudenreichii strains could contribute to the observed host range restrictions, and although we did not detect spontaneous release of viable phages from our host strains (L. Marinelli, data not shown), further research will be needed to rule out this possibility.
Finally, it is also possible that the lack of phage sensitivity among strains of P. freudenreichii may also be due, in part, to other factors, such as increased phage receptor diversity on the bacterial cell surface and/or the presence of restriction/modification systems in these bacteria. In addition, a number of P. freudenreichii strains are predicted to contain clustered regularly interspaced short palindromic repeat (CRISPR) loci, as predicted by the CRISPRfinder online program (http://crispr.i2bc.paris-saclay.fr/Server/index.php, data not shown), and these may play a role in the defense against phage infection in these bacteria. Thus, future studies will focus on investigating this possibility and determining the mechanism(s) that contribute to phage resistance in P. freudenreichii.
Collectively, these findings suggest substantial differences in the evolution of C. acnes vs. P. freudenreichii and their phages, and we hypothesize that these have arisen due to the distinct ecological niches inhabited by these organisms. C. acnes reside within the pilosebaceous follicles of human skin, a lipid-rich, oxygen-limited environment that restricts growth of most other organisms. As such, this bacterium and its phages have likely evolved in isolation with few opportunities for genetic exchange and recombination, which may be reflected in the limited diversity and lack of lysogeny in this population. Conversely, P. freudenreichii and its phages are present both in soil environments and in industrial settings, and are also likely to transit through the human GI tract, and in all of these places, they are likely to encounter numerous other microbial species. Thus, we postulate that this dynamic environment has selected for increased diversity and the potential for lysogeny in this group of phages. Future investigations will be focused on determining the frequency with which these phages are able to lysogenize sensitive host strains and whether widespread lysogeny is a primary mediator of the restricted host ranges observed for these phages.
This report, by describing the complete genomes of seven dsDNA P. freudenreichii phages, has broadened our understanding of the phages specific for members of the Propionibacteriaceae family . We propose that these data, together with the existing P. freudenreichii phage genomes, may be utilized to further explore how phages influence the production of Swiss-type cheeses and potentially, to develop tools for accelerating Propionibacterium-mediated cheese ripening, which is an expensive step during cheese production. In addition, our host range data have identified bacterial strains that are highly resistant to different types P. freudenreichii phages. Unlike the C. acnes phages, all sequenced P. freudenreichii phages contain genes encoding a putative integrase, and our data suggest that lysogeny may be a strategy commonly employed by these phages. Therefore, this information could be utilized to design integration vectors, a system that has been widely employed for genetic manipulation in the related mycobacteria , and which may be used to improve the limited genetic tools available for both P. freudenreichii and C. acnes.
American type culture collection
Clustered regularly interspaced short palindromic repeat
Mitsubishi gas chemical co.
Technology center for genomics and bioinformatics
Yeast extract sodium lactate
Pope WH, Bowman CA, Russell DA, Jacobs-Sera D, Asai DJ, Cresawn SG, Jacobs WR, Hendrix RW, Lawrence JG, Hatfull GF. Whole genome comparison of a large collection of mycobacteriophages reveals a continuum of phage genetic diversity. Elife. 2015;4:e06416.
Kwan T, Liu J, Dubow M, Gros P, Pelletier J. The complete genomes and proteomes of 27 Staphylococcus aureus bacteriophages. Proc Natl Acad Sci U S A. 2005;102(14):5174–9.
Kwan T, Liu J, Dubow M, Gros P, Pelletier J. Comparative genomic analysis of 18 Pseudomonas aeruginosa bacteriophages. J Bacteriol. 2006;188(3):1184–7.
Scholz CF, Kilian M. The natural history of cutaneous propionibacteria, and reclassification of selected species within the genus Propionibacterium to the proposed novel genera Acidipropionibacterium gen. nov., Cutibacterium gen. nov. and Pseudopropionibacterium gen. nov. Int J Syst Evol Microbiol. 2016;66(11):4422–32.
Marinelli LJ, Fitz-Gibbon S, Hayes C, Bowman C, Inkeles M, Loncaric A, Russell DA, Jacobs-Sera D, Cokus S, Pellegrini M, et al. Propionibacterium acnes bacteriophages display limited genetic diversity and broad killing activity against bacterial skin isolates. MBio. 2012;3(5):e00279–12.
Wood HG. Metabolic cycles in the fermentation by propionic acid bacteria. Curr Top Cell Regul. 1981;18:255–87.
Downes J, Wade WG. Propionibacterium acidifaciens sp. nov., isolated from the human mouth. Int J Syst Evol Microbiol. 2009;59(Pt 11):2778–81.
de Freitas R, Chuat V, Madec MN, Nero LA, Thierry A, Valence F, de Carvalho AF. Biodiversity of dairy Propionibacterium isolated from dairy farms in Minas Gerais, Brazil. Int J Food Microbiol. 2015;203:70–7.
Falentin H, Deutsch SM, Jan G, Loux V, Thierry A, Parayre S, Maillard MB, Dherbecourt J, Cousin FJ, Jardin J, et al. The complete genome of Propionibacterium freudenreichii CIRM-BIA1, a hardy actinobacterium with food and probiotic applications. PLoS One. 2010;5(7):e11748.
Thierry A, Deutsch SM, Falentin H, Dalmasso M, Cousin FJ, Jan G. New insights into physiology and metabolism of Propionibacterium freudenreichii. Int J Food Microbiol. 2011;149(1):19–27.
Gagnaire V, Jardin J, Rabah H, Briard-Bion V, Jan G. Emmental Cheese Environment Enhances Propionibacterium freudenreichii Stress Tolerance. PLoS One. 2015;10(8):e0135780.
Cousin FJ, Foligne B, Deutsch SM, Massart S, Parayre S, Le Loir Y, Boudry G, Jan G. Assessment of the probiotic potential of a dairy product fermented by Propionibacterium freudenreichii in piglets. J Agric Food Chem. 2012;60(32):7917–27.
Cousin FJ, Jouan-Lanhouet S, Theret N, Brenner C, Jouan E. Le Moigne-Muller G, Dimanche-Boitrel MT, Jan G: The probiotic Propionibacterium freudenreichii as a new adjuvant for TRAIL-based therapy in colorectal cancer. Oncotarget. 2016;7(6):7161–78.
Jan G, Belzacq AS, Haouzi D, Rouault A, Metivier D, Kroemer G, Brenner C. Propionibacteria induce apoptosis of colorectal carcinoma cells via short-chain fatty acids acting on mitochondria. Cell Death Differ. 2002;9(2):179–88.
Gautier M, Rouault A, Lortal S, Leloir Y, Patureau D. Characterization of a phage infecting Propionibacterium freudenreichii. Lait. 1992;72:431–5.
Gautier M, Rouault A, Sommer P, Briandet R. Occurrence of Propionibacterium freudenreichii bacteriophages in swiss cheese. Appl Environ Microbiol. 1995;61(7):2572–6.
Chopin MC, Rouault A, Ehrlich SD, Gautier M. Filamentous phage active on the gram-positive bacterium Propionibacterium freudenreichii. J Bacteriol. 2002;184(7):2030–3.
Herve C, Coste A, Rouault A, Fraslin JM, Gautier M. First evidence of lysogeny in Propionibacterium freudenreichii subsp. shermanii. Appl Environ Microbiol. 2001;67(1):231–8.
Brown TL, Tucci J, Dyson ZA, Lock P, Adda CG, Petrovski S. Dynamic interactions between prophages induce lysis in Propionibacterium acnes. Res Microbiol. 2017;168(2):103–12.
Malik AC, Reinbold GW, Vedamuthu ER. An evaluation of the taxonomy of Propionibacterium. Can J Microbiol. 1968;14(11):1185–91.
Titus E, Galdone P. Anatole. New York: McGraw-Hill; 1956.
Pope WH, Jacobs-Sera D, Russell DA, Peebles CL, Al-Atrache Z, Alcoser TA, Alexander LM, Alfano MB, Alford ST, Amy NE, et al. Expanding the Diversity of Mycobacteriophages: Insights into Genome Architecture and Evolution. PLoS One. 2011;6(1):e16329.
Cresawn SG, Bogel M, Day N, Jacobs-Sera D, Hendrix RW, Hatfull GF. Phamerator: a bioinformatic tool for comparative bacteriophage genomics. BMC Bioinformatics. 2011;12:395.
Katsura I, Hendrix RW. Length determination in bacteriophage lambda tails. Cell. 1984;39(3 Pt 2):691–8.
Hatfull GF. Molecular Genetics of Mycobacteriophages. Microbiol Spectr. 2014;2(2):1–36.
Meurice G, Jacob D, Deborde C, Chaillou S, Rouault A, Leverrier P, Jan G, Thierry A, Maillard MB, Amet P, et al. Whole genome sequencing project of a dairy Propionibacterium freudenreichii subsp. shermanii genome: progress and first bioinformatic analysis. Lait. 2004;84:15–24.
Jacobs-Sera D, Marinelli LJ, Bowman C, Broussard GW, Guerrero Bustamante C, Boyle MM, Petrova ZO, Dedrick RM, Pope WH. Science Education Alliance Phage Hunters Advancing G et al: On the nature of mycobacteriophage diversity and host preference. Virology. 2012;434(2):187–201.
Hatfull GF, Jacobs-Sera D, Lawrence JG, Pope WH, Russell DA, Ko CC, Weber RJ, Patel MC, Germane KL, Edgar RH, et al. Comparative Genomic Analysis of 60 Mycobacteriophage Genomes: Genome Clustering, Gene Acquisition, and Gene Size. J Mol Biol. 2010;397(1):119–43.
Russell DA, Hatfull GF. PhagesDB: the actinobacteriophage database. Bioinformatics. 2017;33(5):784–6.
Fraser JS, Yu Z, Maxwell KL, Davidson AR. Ig-like domains on bacteriophages: a tale of promiscuity and deceit. J Mol Biol. 2006;359(2):496–507.
Hatfull GF. The secret lives of mycobacteriophages. Adv Virus Res. 2012;82:179–288.
Siguier P, Gourbeyre E, Varani A, Ton-Hoang B, Chandler M. Everyman's Guide to Bacterial Insertion Sequences. Microbiol Spectr. 2015;3(2):MDNA3-0030-2014.
Pedulla ML, Ford ME, Houtz JM, Karthikeyan T, Wadsworth C, Lewis JA, Jacobs-Sera D, Falbo J, Gross J, Pannunzio NR, et al. Origins of highly mosaic mycobacteriophage genomes. Cell. 2003;113(2):171–82.
Huson DH. SplitsTree: analyzing and visualizing evolutionary data. Bioinformatics. 1998;14(1):68–73.
Liu J, Yan R, Zhong Q, Ngo S, Bangayan NJ, Nguyen L, Lui T, Liu M, Erfe MC, Craft N, et al. The diversity and host interactions of Propionibacterium acnes bacteriophages on human skin. ISME J. 2015;9(9):2078–93.
Brown TL, Petrovski S, Dyson ZA, Seviour R, Tucci J. The Formulation of Bacteriophage in a Semi Solid Preparation for Control of Propionibacterium acnes Growth. PLoS One. 2016;11(3):e0151184.
Wolfe BE, Button JE, Santarelli M, Dutton RJ. Cheese rind communities provide tractable systems for in situ and in vitro studies of microbial diversity. Cell. 2014;158(2):422–33.
Farrar MD, Howson KM, Bojar RA, West D, Towler JC, Parry J, Pelton K, Holland KT. Genome sequence and analysis of a Propionibacterium acnes bacteriophage. J Bacteriol. 2007;189(11):4161–7.
Lood R, Collin M. Characterization and genome sequencing of two Propionibacterium acnes phages displaying pseudolysogeny. BMC Genomics. 2011;12:198.
Krumsiek J, Arnold R, Rattei T. Gepard: a rapid and sensitive tool for creating dotplots on genome scale. Bioinformatics. 2007;23(8):1026–8.
We would like to thank Yasinee Rotratsirikun and Rubeen Virani for excellent technical support, and Jessica Jeong for assisting with the host range assays and phage lysate preparation. We also thank the Technology Center for Genomics and Bioinformatics (TCGB) at UCLA for providing sequencing services and for technical support. This study was inspired by the book Anatole by Eve Titus.
This work was partially funded through grants (K01AR066571, F32AR060655) from the National Institutes of Arthritis and Musculoskeletal and Skin Diseases of the National Institutes of Health. The funding bodies played no role in the design of the study, in the collection, analysis, or interpretation of data, or in writing the manuscript.
Availability of data and materials
The bacteriophage sequence data that were generated in the current study are available on GenBank, with the following accession numbers: KX620750, KX620748, KX620752, KX620751, KX620753, KX620754, KX620749. Phylogeny data have been deposited in the TreeBASE repository and are available at the following Accession URL: http://purl.org/phylo/treebase/phylows/study/TB2:S21946
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. P. freudenreichii strains used in this study. Table S2. Average nucleotide identities of P. freudenreichii phages. (DOCX 71 kb)
Figure S1. Genome organization of phage B22. Predicted genes are shown as boxes either above or below the genome corresponding to rightwards- and leftwards-transcription, respectively. The gene numbers are shown within each colored box and its phamily number is shown above with the number of phamily members shown in parentheses; coloring reflects the phamily assignment. Grouping of genes into phamily of related sequences and map generation was performed using Phamerator (30) and the database ‘Actinobacteriophage_685’. Putative gene functions are listed above the genes. (PDF 164 kb)
Figure S2. Genome organization of phage E6. See Additional file 2: Figure S1 for details. (PDF 161 kb)
Figure S3. Genome organization of phage G4. See Additional file 2: Figure S1 for details. (PDF 185 kb)
Figure S4. Genome organization of phage E1. See Additional file 2: Figure S1 for details. (PDF 152 kb)
Figure S5. Genome organization of phage B3. See Additional file 2: Figure S1 for details. (PDF 159 kb)
Figure S6. Phylogenetic tree of tandemly repeated genes in P. freudenreichii phages. The four Cluster BW phages all encode two copies of related genes that are tandemly repeated in Doucette, B22, and E6, but which are separated by 17 genes in G4. The Cluster BV phage Anatole contains a single copy of this gene. Genes duplicated in one genome (e.g. Doucette 38 and 39) are more distantly related (51% aa identity) that with the corresponding gene in another genome (e.g. Doucette gp39 and B22 gp38 share 93% aa identity). Sequences were aligned using ClustalX and the tree drawn using NJPlot. Bootstrap values from 1000 iterations are shown. (PDF 22 kb)