Sequence analysis and molecular characterization of the Lactococcus lactis temperate bacteriophage BK5-T

Exploring Structural Diversity among Adhesion Devices Encoded by Lactococcal P335 Phages with AlphaFold2

Bacteriophages, or phages, are the most abundant biological entities on Earth. They possess molecular nanodevices to package and store their genome, as well as to introduce it into the cytoplasm of their bacterial prey. Successful phage infection commences with specific recognition of, and adhesion to, a suitable host cell surface. Adhesion devices of siphophages infecting Gram-positive bacteria are very diverse and remain, for the majority, poorly understood. These assemblies often comprise long, flexible, and multi-domain proteins, which limit their structural analyses by experimental approaches. The protein structure prediction program AlphaFold2 is exquisitely adapted to unveil structural and functional details of such molecular machineries. Here, we present structure predictions of adhesion devices from siphophages belonging to the P335 group infecting Lactococcus spp., one of the most extensively applied lactic acid bacteria in dairy fermentations. The predictions of representative adhesion devices from types I-IV P335 phages illustrate their very diverse topology. Adhesion devices from types III and IV phages share a common topology with that of Skunavirus p2, with a receptor binding protein anchored to the virion by a distal tail protein loop. This suggests that they exhibit an activation mechanism similar to that of phage p2 prior to host binding.

Sequence and comparative analysis of Leuconostoc dairy bacteriophages

Bacteriophages attacking Leuconostoc species may significantly influence the quality of the final product. There is however limited knowledge of this group of phages in the literature. We have determined the complete genome sequences of nine Leuconostoc bacteriophages virulent to either Leuconostoc mesenteroides or Leuconostoc pseudomesenteroides strains. The phages have dsDNA genomes with sizes ranging from 25.7 to 28.4 kb. Comparative genomics analysis helped classify the 9 phages into two classes, which correlates with the host species. High percentage of similarity within the classes on both nucleotide and protein levels was observed. Genome comparison also revealed very high conservation of the overall genomic organization between the classes. The genes were organized in functional modules responsible for replication, packaging, head and tail morphogenesis, cell lysis and regulation and modification, respectively. No lysogeny modules were detected. To our knowledge this report provides the first comparative genomic work done on Leuconostoc dairy phages.

Identification of a new P335 subgroup through molecular analysis of lactococcal phages Q33 and BM13

Lactococcal dairy starter strains are under constant threat from phages in dairy fermentation facilities, especially by members of the so-called 936, P335, and c2 species. Among these three phage groups, members of the P335 species are the most genetically diverse. Here, we present the complete genome sequences of two P335-type phages, Q33 and BM13, isolated in North America and representing a novel lineage within this phage group. The Q33 and BM13 genomes exhibit homology, not only to P335-type, but also to elements of the 936-type phage sequences. The two phage genomes also have close relatedness to phages infecting Enterococcus and Clostridium, a heretofore unknown feature among lactococcal P335 phages. The Q33 and BM13 genomes are organized in functionally related clusters with genes encoding functions such as DNA replication and packaging, morphogenesis, and host cell lysis. Electron micrographic analysis of the two phages highlights the presence of a baseplate more reminiscent of the baseplate of 936 phages than that of the majority of members of the P335 group, with the exception of r1t and LC3.

Mobile CRISPR/Cas-mediated bacteriophage resistance in Lactococcus lactis

Lactococcus lactis is a biotechnological workhorse for food fermentations and potentially therapeutic products and is therefore widely consumed by humans. It is predominantly used as a starter microbe for fermented dairy products, and specialized strains have adapted from a plant environment through reductive evolution and horizontal gene transfer as evidenced by the association of adventitious traits with mobile elements. Specifically, L. lactis has armed itself with a myriad of plasmid-encoded bacteriophage defensive systems to protect against viral predation. This known arsenal had not included CRISPR/Cas (clustered regularly interspaced short palindromic repeats/CRISPR-associated proteins), which forms a remarkable microbial immunity system against invading DNA. Although CRISPR/Cas systems are common in the genomes of closely related lactic acid bacteria (LAB), none was identified within the eight published lactococcal genomes. Furthermore, a PCR-based search of the common LAB CRISPR/Cas systems (Types I and II) in 383 industrial L. lactis strains proved unsuccessful. Here we describe a novel, Type III, self-transmissible, plasmid-encoded, phage-interfering CRISPR/Cas discovered in L. lactis. The native CRISPR spacers confer resistance based on sequence identity to corresponding lactococcal phage. The interference is directed at phages problematic to the dairy industry, indicative of a responsive system. Moreover, targeting could be modified by engineering the spacer content. The 62.8-kb plasmid was shown to be conjugally transferrable to various strains. Its mobility should facilitate dissemination within microbial communities and provide a readily applicable system to naturally introduce CRISPR/Cas to industrially relevant strains for enhanced phage resistance and prevention against acquisition of undesirable genes.

Sequence analysis of Leuconostoc mesenteroides bacteriophage Phi1-A4 isolated from an industrial vegetable fermentation

Vegetable fermentations rely on the proper succession of a variety of lactic acid bacteria (LAB). Leuconostoc mesenteroides initiates fermentation. As fermentation proceeds, L. mesenteroides dies off and other LAB complete the fermentation. Phages infecting L. mesenteroides may significantly influence the die-off of L. mesenteroides. However, no L. mesenteroides phages have been previously genetically characterized. Knowledge of more phage genome sequences may provide new insights into phage genomics, phage evolution, and phage-host interactions. We have determined the complete genome sequence of L. mesenteroides phage Phi1-A4, isolated from an industrial sauerkraut fermentation. The phage possesses a linear, double-stranded DNA genome consisting of 29,508 bp with a G+C content of 36%. Fifty open reading frames (ORFs) were predicted. Putative functions were assigned to 26 ORFs (52%), including 5 ORFs of structural proteins. The phage genome was modularly organized, containing DNA replication, DNA-packaging, head and tail morphogenesis, cell lysis, and DNA regulation/modification modules. In silico analyses showed that Phi1-A4 is a unique lytic phage with a large-scale genome inversion ( approximately 30% of the genome). The genome inversion encompassed the lysis module, part of the structural protein module, and a cos site. The endolysin gene was flanked by two holin genes. The tail morphogenesis module was interspersed with cell lysis genes and other genes with unknown functions. The predicted amino acid sequences of the phage proteins showed little similarity to other phages, but functional analyses showed that Phi1-A4 clusters with several Lactococcus phages. To our knowledge, Phi1-A4 is the first genetically characterized L. mesenteroides phage.

Morphology, genome sequence, and structural proteome of type phage P335 from Lactococcus lactis

Lactococcus lactis phage P335 is a virulent type phage for the species that bears its name and belongs to the Siphoviridae family. Morphologically, P335 resembled the L. lactis phages TP901-1 and Tuc2009, except for a shorter tail and a different collar/whisker structure. Its 33,613-bp double-stranded DNA genome had 50 open reading frames. Putative functions were assigned to 29 of them. Unlike other sequenced genomes from lactococcal phages belonging to this species, P335 did not have a lysogeny module. However, it did carry a dUTPase gene, the most conserved gene among this phage species. Comparative genomic analyses revealed a high level of identity between the morphogenesis modules of the phages P335, ul36, TP901-1, and Tuc2009 and two putative prophages of L. lactis SK11. Differences were noted in genes coding for receptor-binding proteins, in agreement with their distinct host ranges. Sixteen structural proteins of phage P335 were identified by liquid chromatography-tandem mass spectrometry. A 2.8-kb insertion was recognized between the putative genes coding for the activator of late transcription (Alt) and the small terminase subunit (TerS). Four genes within this region were autonomously late transcribed and possibly under the control of Alt. Three of the four deduced proteins had similarities with proteins from Streptococcus pyogenes prophages, suggesting that P335 acquired this module from another phage genome. The genetic diversity of the P335 species indicates that they are exceptional models for studying the modular theory of phage evolution.

Characterization of 1706, a virulent phage from Lactococcus lactis with similarities to prophages from other Firmicutes

The virulent lactococcal phage 1706, isolated in 1995 from a failed cheese production in France, represents a new lactococcal phage species of the Siphoviridae family. This phage has a burst size of 160 and a latent period of 85 min. Its linear double-stranded DNA genome was composed of 55,597 bp with a 33.7% G+C content. Its deduced proteome (76 ORFs) shared limited similarities to other known phage proteins. SDS-PAGE coupled with LC-MS/MS analyses led to the identification of 15 structural proteins. The most striking feature of the 1706 proteome was that 22 ORFs shared similarities with proteins deduced from the genome of either Ruminococcus torques and/or Clostridium leptum. Both are Firmicutes bacteria found in the gut flora of humans. We also identified a four-gene module in phage 1706, most likely involved in host recognition that shared similarities with lactococcal prophages. We propose that the virulent phage 1706 infected another bacterial genus before picking up a lactococcal host recognition module.

The genome of epsilon15, a serotype-converting, Group E1 Salmonella enterica-specific bacteriophage

The genome sequence of the Salmonella enterica serovar Anatum-specific, serotype-converting bacteriophage epsilon15 has been completed. The nonredundant genome contains 39,671 bp and 51 putative genes. It most closely resembles the genome of phiV10, an Escherichia coli O157:H7-specific temperate phage, with which it shares 36 related genes. More distant relatives include the Burkholderia cepacia-specific phage, BcepC6B (8 similar genes), the Bordetella bronchiseptica-specific phage, BPP-1 (8 similar genes) and the Photobacterium profundum prophage, P Pphipr1 (6 similar genes). epsilon15 gene identifications based on homologies with known gene families include the terminase small and large subunits, integrase, endolysin, two holins, two DNA methylase enzymes (one adenine-specific and one cytosine-specific) and a RecT-like enzyme. Genes identified experimentally include those coding for the serotype conversion proteins, the tail fiber, the major capsid protein and the major repressor. epsilon15's attP site and the Salmonella attB site with which it interacts during lysogenization have also been determined.

Comparative analyses of prophage-like elements present in two Lactococcus lactis strains

In this study, we describe the genetic organizations of six and five apparent prophage-like elements present in the genomes of the Lactococcus lactis subsp. cremoris strains MG1363 and SK11, respectively. Phylogenetic investigation as well bioinformatic analyses indicates that all 11 prophages belong to subdivisions of the lactococcal P335 group of temperate bacteriophages.

Curated list of prokaryote viruses with fully sequenced genomes

Genome sequencing is of enormous importance for classification of prokaryote viruses and for understanding the evolution of these viruses. This survey covers 284 sequenced viruses for which a full description has been published and for which the morphology is known. This corresponds to 219 (4%) of tailed and 75 (36%) of tailless viruses of prokaryotes. The number of sequenced tailless viruses almost doubles if viruses of unknown morphology are counted. The sequences are from representatives of 15 virus families and three groups without family status, including eight taxa of archaeal viruses. Tailed phages, especially those with large genomes and hosts other than enterobacteria or lactococci, mycobacteria and pseudomonads, are vastly under investigated.

Genome sequence and global gene expression of Q54, a new phage species linking the 936 and c2 phage species of Lactococcus lactis

The lytic lactococcal phage Q54 was previously isolated from a failed sour cream production. Its complete genomic sequence (26,537 bp) is reported here, and the analysis indicated that it represents a new Lactococcus lactis phage species. A striking feature of phage Q54 is the low level of similarity of its proteome (47 open reading frames) with proteins in databases. A global gene expression study confirmed the presence of two early gene modules in Q54. The unusual configuration of these modules, combined with results of comparative analysis with other lactococcal phage genomes, suggests that one of these modules was acquired through recombination events between c2- and 936-like phages. Proteolytic cleavage and cross-linking of the major capsid protein were demonstrated through structural protein analyses. A programmed translational frameshift between the major tail protein (MTP) and the receptor-binding protein (RBP) was also discovered. A "shifty stop" signal followed by putative secondary structures is likely involved in frameshifting. To our knowledge, this is only the second report of translational frameshifting (+1) in double-stranded DNA bacteriophages and the first case of translational coupling between an MTP and an RBP. Thus, phage Q54 represents a fascinating member of a new species with unusual characteristics that brings new insights into lactococcal phage evolution.

Complete nucleotide sequence and genome analysis of bacteriophage BFK20 — A lytic phage of the industrial producer Brevibacterium flavum

The entire double-stranded DNA genome of bacteriophage BFK20, a lytic phage of the Brevibacterium flavum CCM 251--industrial producer of L-lysine--was sequenced and analyzed. It consists of 42,968 base pairs with an overall molar G + C content of 56.2%. Fifty-five potential open reading frames were identified and annotated using various bioinformatics tools. Clusters of functionally related putative genes were defined (structural, lytic, replication and regulatory). To verify the annotation of structural proteins, they were resolved by 2D gel electrophoresis and were submitted to N-terminal amino acid sequencing. Structural proteins identified included the portal and major and minor tail proteins. Based on the overall genome sequence comparison, similarities with other known bacteriophage genomes include primarily bacteriophages from Mycobacterium spp. and some regions of Corynebacterium spp. genomes--possible prophages. Our results support the theory that phage genomes are mosaics with respect to each other.

Genome analysis of the obligately lytic bacteriophage 4268 of Lactococcus lactis provides insight into its adaptable nature

Analysis of the complete nucleotide sequence of the lactococcal phage 4268, which is lytic for the cheese starter Lactococcus lactis DPC4268, is presented. Phage 4268 has a linear genome of 36,596 bp, which is modularly organised and encompasses 49 open reading frames. Putative functions were assigned to approximately 45% of the predicted products of these open reading frames based on sequence similarity with known proteins, N-terminal sequence analysis and identification of conserved domains. Significantly, a segment of the genome has homology to the recently sequenced lysogenic module in lactococcal phage phi31 that contains a lytic switch but no phage integrase or attachment site. This suggests that it is derived from a prophage. A phage 4268-encoded and a host-encoded methylase were found to be highly similar, having only two nucleotide mismatches, suggesting that the phage acquired the methylase gene to protect it from a host endonuclease. Comparative genomic analysis revealed significant homology between phage 4268 and the lactococcal phage BK5-T. The comparative analysis also supported the classification of phage 4268 and other BK5-T-related phage as separate from the proposed P335 species of lactococcal phage.

Structural characterization and assembly of the distal tail structure of the temperate lactococcal bacteriophage TP901-1

The tail structures of bacteriophages infecting gram-positive bacteria are largely unexplored, although the phage tail mediates the initial interaction with the host cell. The temperate Lactococcus lactis phage TP901-1 of the Siphoviridae family has a long noncontractile tail with a distal baseplate. In the present study, we investigated the distal tail structures and tail assembly of phage TP901-1 by introducing nonsense mutations into the late transcribed genes dit (orf46), tal(TP901-1) (orf47), bppU (orf48), bppL (orf49), and orf50. Transmission electron microscopy examination of mutant and wild-type TP901-1 phages showed that the baseplate consisted of two different disks and that a central tail fiber is protruding below the baseplate. Evaluation of the mutant tail morphologies with protein profiles and Western blots revealed that the upper and lower baseplate disks consist of the proteins BppU and BppL, respectively. Likewise, Dit and Tal(TP901-1) were shown to be structural tail proteins essential for tail formation, and Tal(TP901-1) was furthermore identified as the tail fiber protein by immunogold labeling experiments. Determination of infection efficiencies of the mutant phages showed that the baseplate is fundamental for host infection and the lower disk protein, BppL, is suggested to interact with the host receptor. In contrast, ORF50 was found to be nonessential for tail assembly and host infection. A model for TP901-1 tail assembly, in which the function of eight specific proteins is considered, is presented.

Sequence analysis of the Lactobacillus plantarum bacteriophage PhiJL-1

The complete genomic sequence of a Lactobacillus plantarum virulent phage PhiJL-1 was determined. The phage possesses a linear, double-stranded, DNA genome consisting of 36,677 bp with a G+C content of 39.36%. A total of 52 possible open reading frames (ORFs) were identified. According to N-terminal amino acid sequencing and bioinformatic analyses, proven or putative functions were assigned to 21 ORFs (41%), including 5 structural protein genes. The PhiJL-1 genome shows functionally related genes clustered together in a genome structure composed of modules for DNA replication, DNA packaging, head and tail morphogenesis, and lysis. This type of modular genomic organization was similar to several other phages infecting lactic acid bacteria. The structural gene maps revealed that the order of the head and tail genes is highly conserved among the genomes of several Siphoviridae phages, allowing the assignment of probable functions to certain uncharacterized ORFs from phage PhiJL-1 and other Siphoviridae phages.

Identification of the receptor-binding protein in 936-species lactococcal bacteriophages

The aim of this work was to identify genes responsible for host recognition in the lactococcal phages sk1 and bIL170 belonging to species 936. These phages have a high level of DNA identity but different host ranges. Bioinformatic analysis indicated that homologous genes, orf18 in sk1 and orf20 in bIL170, could be the receptor-binding protein (RBP) genes, since the resulting proteins were unrelated in the C-terminal part and showed homology to different groups of proteins hypothetically involved in host recognition. Consequently, chimeric bIL170 phages carrying orf18 from sk1 were generated. The recombinant phages were able to form plaques on the sk1 host Lactococcus lactis MG1614, and recombination was verified by PCR analysis directly with the plaques. A polyclonal antiserum raised against the C-terminal part of phage sk1 ORF18 was used in immunogold electron microscopy to demonstrate that ORF18 is located at the tip of the tail. Sequence analysis of corresponding proteins from other lactococcal phages belonging to species 936 showed that the N-terminal parts of the RBPs were very similar, while the C-terminal parts varied, suggesting that the C-terminal part plays a role in receptor binding. The phages investigated could be grouped into sk1-like phages (p2, fd13, jj50, and phi 7) and bIL170-like phages (P008, P113G, P272, and bIL66) on the basis of the homology of their RBPs to the C-terminal part of ORF18 in sk1 and ORF20 in bIL170, respectively. Interestingly, sk1-like phages bind to and infect a defined group of L. lactis subsp. cremoris strains, while bIL170-like phages bind to and infect a defined group of L. lactis subsp. lactis strains.

Essentiality of the early transcript in the replication origin of the lactococcal prolate phage c2

The genome of the prolate-headed lytic lactococcal bacteriophage c2 is organized into two divergently oriented blocks consisting of the early genes and the late genes. These blocks are separated by the noncoding origin of DNA replication. We examined the functional role of transcription of the origin in a plasmid model system. Deletion of the early promoter P(E)1 abolished origin function. Introduction of mutations into P(E)1 which did not eliminate promoter activity or replacement of P(E)1 with an unrelated but functional promoter did not abolish replication. The A-T-rich region upstream of P(E)1, which is conserved in prolate phages, was not required for plasmid replication. Replacement of the P(E)1 transcript template sequence with an unrelated sequence with a similar G+C content abolished replication, showing that the sequence encoding the transcript is essential for origin function. Truncated transcript and internal deletion constructs did not support replication except when the deletion was at the very 3' end of the DNA sequence coding for the transcript. The P(E)1 transcript could be detected for all replication-proficient constructs. Recloning in a plasmid vector allowed detection of P(E)1 transcripts from some fragments that did not support replication, indicating that stability of the transcript alone was not sufficient for replication. The data suggest that production of a transcript of a specific length and with a specific sequence or structure is essential for the function of the phage c2 origin in this model system.

The genome and proteome of coliphage T1

The genome of enterobacterial phage T1 has been sequenced, revealing that its 50.7-kb terminally redundant, circularly permuted sequence contains 48,836 bp of nonredundant nucleotides. Seventy-seven open reading frames (ORFs) were identified, with a high percentage of small genes located at the termini of the genomes displaying no homology to existing phage or prophage proteins. Of the genes showing homologs (47%), we identified those involved in host DNA degradation (three endonucleases) and T1 replication (DNA helicase, primase, and single-stranded DNA-binding proteins) and recombination (RecE and Erf homologs). While the tail genes showed homology to those from temperate coliphage N15, the capsid biosynthetic genes were unique. Phage proteins were resolved by 2D gel electrophoresis, and mass spectrometry was used to identify several of the spots including the major head, portal, and tail proteins, thus verifying the annotation.

Distribution and composition of the lysis cassette ofLactococcus lactisphages and functional analysis of bacteriophage ul36 holin

The bacteriophage lysis cassette, which comprises a lysin and a holin gene, was analyzed in 18 Lactococcus lactis phages. A muramidase motif was found in the lysins of c2-like phages, while an amidase motif was observed in the lysins of 936-like phages. Both amidase and muramidase types were detected among the P335 phages. The P335 lysins were separated into three groups based on amino acid sequence identity. A class I holin was recognized in 936-like and c2-like phages, whereas P335-like phages possess class II holins. The P335 holins were further divided into four groups based on sequence identity. Only the holins of 936-like phages contained putative dual-start motifs. The unusual lysis cassette of the highly virulent P335-like phage ul36 contains a unique holin (orf74B) upstream of a lysin which is present in several other P335-like phages. Using the lambdadelta Sthf system, we demonstrated that gpORF74B induces cell lysis at the same time as lambdadelta Sthf::S105, the effector of lambda lysis. Transcriptional analysis of ul36 lysis cassette showed that first transcripts are detected 35 min after infection of L. lactis cells. The lysis clock of phage ul36 appears to be controlled by the late expression of the holin and lysin genes.

Transcriptional analysis of the genetic elements involved in the lysogeny/lysis switch in the temperate lactococcal bacteriophage phiLC3, and identification of the Cro-like protein ORF76

A transcriptional analysis of the lysogeny-related genes of the temperate bacteriophage Lactococcus lactis phiLC3 was performed using Northern blot hybridization during lysogeny and lytic infection by the phage. The lysogeny-related gene cluster was found to contain four promoters (P(1), P(2), Pint and P(173)), while the P(87) promoter directed transcription of orf80 and the putative gene orf87, which are located between the integrase gene and the cell lysis genes. The start sites of the transcripts were determined by primer extension. The divergently oriented lysogenic P(1) and lytic P(2) promoters located in the genetic switch region are responsible for transcription of orf286 which encodes the phage repressor, and the genes orf63 - orf76 - orf236 - orf110 - orf82 - orf57, respectively, while orf173 is transcribed from P(173). orf76 was identified as the gene encoding the Cro-like protein of phiLC3, and it was shown that ORF76 is able to bind specifically to the genetic switch region, albeit with lower affinity than does the phage repressor ORF286. ORF76 also competed with ORF286 for binding to this region. The functionality of P(1) and P(2), and their regulation by ORF286 and ORF76, was investigated using a reporter gene. In general, P(2) was a stronger promoter than P(1), but expression from both promoters, especially P(2), was regulated and modulated by flanking sequences and the presence of orf286 and orf76. ORF286 and ORF76 were both able to repress transcription from P(1) and P(2), while ORF286 was able to stimulate its own synthesis by tenfold. This work reveals the complex interplay between the regulatory elements that control the genetic switch between lysis and lysogeny in phiLC3 and other temperate phages of Lactococcus.

Prophage genomics

The majority of the bacterial genome sequences deposited in the National Center for Biotechnology Information database contain prophage sequences. Analysis of the prophages suggested that after being integrated into bacterial genomes, they undergo a complex decay process consisting of inactivating point mutations, genome rearrangements, modular exchanges, invasion by further mobile DNA elements, and massive DNA deletion. We review the technical difficulties in defining such altered prophage sequences in bacterial genomes and discuss theoretical frameworks for the phage-bacterium interaction at the genomic level. The published genome sequences from three groups of eubacteria (low- and high-G+C gram-positive bacteria and gamma-proteobacteria) were screened for prophage sequences. The prophages from Streptococcus pyogenes served as test case for theoretical predictions of the role of prophages in the evolution of pathogenic bacteria. The genomes from further human, animal, and plant pathogens, as well as commensal and free-living bacteria, were included in the analysis to see whether the same principles of prophage genomics apply for bacteria living in different ecological niches and coming from distinct phylogenetical affinities. The effect of selection pressure on the host bacterium is apparently an important force shaping the prophage genomes in low-G+C gram-positive bacteria and gamma-proteobacteria.

A comparative study and phage typing of silage-making Lactobacillus bacteriophages

To investigate basic characteristics of 10 virulent phages active on silage-making lactobacilli, morphological properties, host ranges, protein composition and genome characterization were separated into five groups based on host ranges and basic properties. The seven phages of groups I, II and V were active on Lactobacillus plantarum and Lactobacillus pentosus. Phage phiPY4 (group III) infected both L. casei and Lactobacillus rhamnosus. Phage phiPY5 (group IV) specifically infected Lactobacillus casei. Morphologically, three phages of groups I belonged to the Myoviridae family, while seven other phages of groups II, III and V belonged to the Siphoviridae family. SDS-PAGE profiles, restriction analysis, G + C contents of DNA and Dot blot hybridization revealed a high degree of homology in each group. Clustering derived from host range analysis was closely related to results of DNA and protein analyses. These phages may be applicable to phage typing for silage-making lactobacilli.

The dilemma of phage taxonomy illustrated by comparative genomics of Sfi21-like Siphoviridae in lactic acid bacteria

The complete genome sequences of two dairy phages, Streptococcus thermophilus phage 7201 and Lactobacillus casei phage A2, are reported. Comparative genomics reveals that both phages are members of the recently proposed Sfi21-like genus of Siphoviridae, a widely distributed phage type in low-GC-content gram-positive bacteria. Graded relatedness, the hallmark of evolving biological systems, was observed when different Sfi21-like phages were compared. Across the structural module, the graded relatedness was represented by a high level of DNA sequence similarity or protein sequence similarity, or a shared gene map in the absence of sequence relatedness. This varying range of relatedness was found within Sfi21-like phages from a single species as demonstrated by the different prophages harbored by Lactococcus lactis strain IL1403. A systematic dot plot analysis with 11 complete L. lactis phage genome sequences revealed a clear separation of all temperate phages from two classes of virulent phages. The temperate lactococcal phages share DNA sequence homology in a patchwise fashion over the nonstructural gene cluster. With respect to structural genes, four DNA homology groups could be defined within temperate L. lactis phages. Closely related structural modules for all four DNA homology groups were detected in phages from Streptococcus or Listeria, suggesting that they represent distinct evolutionary lineages that have not uniquely evolved in L. lactis. It seems reasonable to base phage taxonomy on data from comparative genomics. However, the peculiar modular nature of phage evolution creates ambiguities in the definition of phage taxa by comparative genomics. For example, depending on the module on which the classification is based, temperate lactococcal phages can be classified as a single phage species, as four distinct phage species, or as two if not three different phage genera. We propose to base phage taxonomy on comparative genomics of a single structural gene module (head or tail genes). This partially phylogeny-based taxonomical system still mirrors some aspects of the current International Committee on Taxonomy in Virology classification system. In this system the currently sequenced lactococcal phages would be grouped into five genera: c2-, sk1, Sfi11-, r1t-, and Sfi21-like phages.

Characterization of the putative replisome organizer of the lactococcal bacteriophage r1t

Analysis of the nucleotide sequence of the genome of the lactococcal bacteriophage r1t showed that it may encode at least two proteins involved in DNA replication. On the basis of its similarity with the G38P protein encoded by the Bacillus subtilis phage SPP1, the product of orf11 (Pro11) is thought to be involved in the initiation of phage DNA replication. This protein was overexpressed in Lactococcus lactis and partially purified. Gel retardation analysis using various r1t DNA fragments indicates that Pro11 specifically binds to a sequence located within its cognate gene. DNase I footprinting showed that Pro11 protects a stretch of DNA of 47 bp. This region spans four 6-bp short direct repeats, which suggests that the region contains four binding sites for Pro11. 1,10-Phenanthroline-copper footprinting confirmed the protection of the hexamers. An asymmetric protection pattern of each strand was observed, suggesting that Pro11 contacts each DNA strand separately at contiguous hexamers. We propose a model for the binding of Pro11 to its target sites that may account for the torsion strain required for strand opening at the origin of replication.

Complete genomic sequence of bacteriophage ul36: demonstration of phage heterogeneity within the P335 quasi-species of lactococcal phages

The complete genomic sequence of the Lactococcus lactis virulent phage ul36 belonging to P335 lactococcal phage species was determined and analyzed. The genomic sequence of this lactococcal phage contained 36,798 bp with an overall G+C content of 35.8 mol %. Fifty-nine open reading frames (ORFs) of more than 40 codons were found. N-terminal sequencing of phage structural proteins as well as bioinformatic analysis led to the attribution of a function to 24 ORFs (41%). A lysogeny module was found within the genome of this virulent phage. The putative integrase gene seems to be the product of a horizontal transfer because it is more closely related to Streptococcus pyogenes phages than it is to L. lactis phages. Comparative genome analysis with six complete genomes of temperate P335-like phages confirmed the heterogeneity among phages of P335 species. A dUTPase gene is the only conserved gene among all P335 phages analyzed as well as the phage BK5-T. A genetic relationship between P335 phages and the phage-type of the BK5-T species was established. Thus, we proposed that phage BK5-T be included within the P335 species and thereby reducing the number of lactococcal phage species to 11.

Sequence analysis of the lactococcal bacteriophage bIL170: insights into structural proteins and HNH endonucleases in dairy phages

The complete 31754 bp genome of bIL170, a virulent bacteriophage of Lactococcus lactis belonging to the 936 group, was analysed. Sixty-four ORFs were predicted and the function of 16 of them was assigned by significant homology to proteins in databases. Three putative homing endonucleases of the HNH family were found in the early region. An HNH endonuclease with zinc-binding motif was identified in the late cluster, potentially being part of the same functional module as terminase. Three putative structural proteins were analysed in detail and show interesting features among dairy phages. Notably, gpl12 (putative fibre) and gpl20 (putative baseplate protein) of bIL170 are related by at least one of their domains to a number of multi-domain proteins encoded by lactococcal or streptococcal phages. A 110- to 150-aa-long hypervariable domain flanked by two conserved motifs of about 20 aa was identified. The analysis presented here supports the participation of some of these proteins in host-range determination and suggests that specific adsorption to the host may involve a complex multi-component system. Divergences in the genome of phages of the 936 group, that may have important biological properties, were noted. Insertions/deletions of units of one or two ORFs were the main source of divergence in the early clusters of the two entirely sequenced phages, bIL170 and sk1. An exchange of fragments probably affected the regions containing the putative origin of replication. It led to the absence in bIL170 of the direct repeats recognized in sk1 and to the presence of different ORFs in the ori region. Shuffling of protein domains affected the endolysin (putative cell-wall binding part), as well as gpl12 and gpl20.