Lactococcal 936-type phages and dairy fermentation problems: from detection to evolution and prevention

The so-called 936-type phages are the most frequently encountered lactococcal phage species in dairy fermentations, where they cause slow or even failed fermentations with concomitant economic losses. Several dairy phage population studies, performed in different geographical locations, have detailed their dominance in dairy phage populations, while various phage-resistance mechanisms have been assessed in a bid to protect against this virulent phage group. The impact of thermal and chemical treatments on 936 phages is an important aspect for dairy technologists and has been assessed in several studies, and has indicated that these phages have adapted to better resist such treatments. The abundance of 936 phage genome sequences has permitted a focused view on genomic content and regions of variation, and the role of such variable regions in the evolution of these phages. Here, we present an overview on detection and global prevalence of the 936 phages, together with their tolerance to industrial treatments and anti-phage strategies. Furthermore, we present a comprehensive review on the comparative genomic analyses of members of this fascinating phage species.


INTRODUCTION
Phages infecting Lactococcus lactis are well-studied and represent a rather diverse group of bacterial viruses that are classified into 10 species based on morphology and sequence relatedness . Eight of these ten species are Siphoviridae members, possessing long non-contractile tails, while the two remaining species belong to the Podoviridae family that are characterized by short tails . Representatives of each of these species have been isolated through phage population studies of (predominantly) dairy facilities, although the most frequently isolated phages belong to the so-called 936, P335, and c2 species. The most frequently isolated phages in lactococcal dairy fermentations belong to the 936-type species and these phages therefore pose the highest threat to this industry. Furthermore, the virulent nature of these phages combined with their increased tolerance to sanitation and thermal treatments highlight the requirement for constant monitoring and analysis of this phage species.

ISOLATION OF 936-TYPE PHAGES FROM GEOGRAPHICALLY DISTINCT LOCATIONS
The 936-type phage species were shown to clearly dominate lactococcal phage isolation studies as deduced from an extensive number of isolation studies across various geographical locations ( Table 1).
In 1993, 22 phages of the 936-type species with its classic morphology of a long non-contractile tail and isometric protein head were isolated from cheese whey (Casey et al., 1993). Similarly, over an 11 month period, twenty-two 936-type phages were isolated from a buttermilk fermentation facility, accounting for 80% of the total phage species isolated (Moineau et al., 1996). Interestingly, the latter study identified raw milk as the major source of phage contamination as buttermilk fermentation occurs in closed vats. It was also observed that 936-type phage isolates from distinct geographical regions had significant similarities in terms of host range and phage restriction profiles. In 1992, 30 lactococcal phages were isolated from Canadian (province of Quebec) cheese factories, of which eleven were found to belong to the 936-type species as determined by electron microscopy, DNA homology, and protein profiles (Moineau et al., 1992). Interestingly, phages belonging to the c2-type species were found to be the most frequently isolated in this particular study. In a more recent survey, phages belonging to the 936-type species have emerged to dominate phage populations in Canadian whey samples. Out of 71 lactococcal phages isolated against a mixed starter culture, the 936 species accounted for 74% and were found to exhibit a wide host spectrum (Bissonnette et al., 2000).
In stark contrast to the Canadian, Irish, and U.S. studies, analysis of phage prevalence of six Belarusian dairy regions showed that just four out of 23 isolated lactococcal phages were shown to belong to the 936-type species (Raiski and Belyasova, 2009). The 936-type phages BIM BV-27, E11, and E12 all had similar restriction profiles yet were isolated from central, western, and south-western regions in Belarus, further emphasizing the widespread proliferation of this phage species with limited changes in their genetic content. This is not precluding the possibility that these fermentation plants obtained their starter cultures from the same supplier(s), thus representing a limited number of strains with a consequently restricted level of phage diversity. Ecological analysis  (Szczepanska et al., 2007). Two particular isolates of the 936-type species were found in more than one location, bIBB5g 1 and bIBBEg 1 , which were isolated from three and seven of the 17 tested locations, respectively. Further affirming the dominance of the 936-type phages in a global context, members of this species were the shown to be the most frequently isolated in Norwegian and U.S. whey samples originating from Dutch type cheese fermentations (Kleppena et al., 2011). In total, twenty-five 936-type phages were detected in the latter study, which were shown to be distinct based on the sequence of a DNA fragment encoding structural phage proteins. This section of the genome contains a unique region that was identified as a mutation hotspot which defined the individual isolates, highlighting the diversification of these phages. Phages have also been isolated from various locations within the dairy facility other than the dairy product production line. Several studies have demonstrated that phages have the ability to become airborne within a dairy facility, more than likely due the use of open vats, processing conditions, and movement of equipment (Neve et al., 1994(Neve et al., , 1995(Neve et al., , 2003Verreault et al., 2011). In the Verreault et al. (2011) study, aerosol samples were taken at the filling section at the end of the production line within a cheese factory, where 936-type phages were detected in the air at levels of 10 4 phage genomes/m 3 , a level significant enough to affect cheese production if access to open vats was to be gained. Furthermore, phages were isolated from surfaces such as doors, cleaning materials, offices, and floors, with the highest concentration of 936-type phages at 9106 ± 789 genomes cm −2 . Such contamination levels create an environment whereby staff could spread phages across the plant, making efforts to control such populations very difficult. Not only are the 936-type phages prominent within the lactococcal dairy industry as evident from biodiversity studies, this particular species appears to possess an intrinsic ability to increase its virulence. Josephsen et al. (1999Josephsen et al. ( ) compared 936-type phages isolated between 1989Josephsen et al. ( and 1994Josephsen et al. ( , and 936-type isolates from 1982Josephsen et al. ( to 1986 showing that both phage sets shared homology to each other but that the more recent isolates exhibited a broader host range and a greater burst size, thus indicating that the phages found in this particular plant had evolved to become more infective. The 936-type phages have also been shown to persist in a specific dairy factory during an extended period of time, as demonstrated by the isolation and subsequent re-isolation of the 936-type phage CB13/GR7 (Rousseau and Moineau, 2009a). Moreover, this study suggested that newly isolated phages were derived in some part from the older isolates, highlighting that the 936 phage species continues to retain dominance due to its genetic elasticity and ability to evolve to a changing environment.

DETECTION OF 936-TYPE PHAGES
Two main methods for the detection and enumeration of phages from industrial dairy plants are used. The first method, which is predominantly used by the dairy industry for phage monitoring purposes, is the so-called Heap-Lawrence test, which is based on the acidification rate of milk and provides a reliable test for the presence of phages, while it also gives an indication of phage titer (Heap and Lawrence, 1976). The second one, used most intensely by research laboratories is the standard doublelayer plaque assay (Lillehaug, 1997), which following overnight incubation of the plates allows for the quantitative analysis of phage levels. In recent years, a number of alternative molecular methods have been developed to identify phages present in dairy samples. One such method is the use of polymerase chain reaction (PCR), which allows the reliable detection and distinction of the three most common phage species (936, c2, and P335) in a single multiplex primer-based PCR reaction . In terms of the 936-type phage, the major structural protein was shown to be highly conserved (approximately 86% identity at DNA level) among the phages examined (p2, Q7, and Q11 sk1, bIL170, and F4-1; Labrie and Moineau, 2000). Primers were designed to target this gene, generating an amplicon of 179 base pairs (bp) and allowing a detection limit of 10 4 pfu/ml (plaque forming units per ml of sample). A similar method was developed (del Rio et al., 2007) based on the msp gene of phages p2 and Q7, the orf11 gene of phage sk1, the orf13 gene of phage bIL170, and Frontiers in Microbiology | Food Microbiology the mcp gene of phage F4-1, all belonging to the 936-type species and lowering the detection limit to 10 3 pfu/ml. While both PCR methods are effective at detecting and discriminating phages in raw milk and whey, they do not permit the isolation of individual phages. The receptor binding protein (RBP) of the 936-type phage is one of the more divergent regions of the genome and is believed to be responsible for the phage's host range (Dupont et al., 2004). Using this information, Dupont et al. (2005) developed an optimized magnetic capture hybridization (MCH) PCR-based assay to detect 936-type phage and to further detect the different subgroups of this species with primers targeting the variable region of the RBP gene at a detection limit of 10 2 -10 3 pfu/ml (Dupont et al., 2005). However, recent studies have highlighted that emerging 936-type phages no longer fall into these originally identified RBP-based subgroups, thus advocating primer redesign for future multiplex-based phage detection and subgroup assignment (Castro-Nallar et al., 2012). Real-time (RT) PCR, another molecular technique that allows the rapid detection and quantification of microorganisms has also been used for the detection 936-type phages (Ly-Chatain et al., 2011;Verreault et al., 2011). As mentioned above, the Verreault et al. (2011) study reported the detection of airborne 936-type phages within a dairy environment using RT-PCR with a detection limit of 10 2 -10 4 pfu/ml. The Ly-Chatain et al. (2011) study applied RT-PCR to DNA isolated from milk and whey, and achieved a detection limit of 10 2 pfu/ml. Both studies would permit rapid sampling and identification of phages in given industrial plant to gage phage contamination levels, however, they do not permit the isolation of individual isolates within a given species and also would not detect the emergence of rare or new phages. One method to overcome such limitations employs the use of flow cytometry (Michelsen et al., 2007). Phage infection of a host cell can be monitored in real time by measuring loss of host cell mass using light scatter and DNA fluorescence staining, allowing a detection limit of 10 4 -10 5 pfu/ml. It is clear from the numerous biodiversity studies on lactococcal phages that a large quantity of diverse 936-type phages can be isolated from any given factory. Host range analyses and restriction profiling of phage genomes is time consuming and may require numerous attempts to achieve suitable patterns for comparative analysis between phages. A multilocus sequence typing (MLST) scheme was developed using primers that target the middle region of conserved genes found in 936-type phages (Moisan and Moineau, 2012). This approach utilizes a single PCR reaction using five primers sets targeting the genes encoding the major capsid protein (MCP), tape measure protein (TMP), lysin (lys), large terminase (terL), and the major tail protein (MTP). Of note, the PCR can be performed on phage lysates, eliminating the need to isolate phage DNA.
A number of methods are now available to both industry and academics alike to characterize the biodiversity of phages in a dairy factory, and to enumerate and detect phages that belong to the 936-type species. For industrial purposes, rapid detection of phages by PCR-based methods may be more applicable so a timely intervention can be introduced to prevent fermentation failure. For research, it is necessary to isolate individual 936-type phages present within a specific environment to allow continued study into why these phages are so prevalent and how they overcome technological hurdles. The standard plaque assay followed by molecular characterization may remain the best methodologies to perform such studies.

PHYSICAL APPROACHES TOWARD INACTIVATION OF 936-TYPE PHAGES
In order to prevent the negative consequences of phage infection in the dairy factory, i.e., the development of slow or dead vats, various counter measures are employed, such as thermal and biocidal treatments, strain rotation regimes and the application of bacteriophage-insensitive strains (Lavigne et al., 2009;Mills et al., 2010). A number of reports have been published on the thermal and biocidal elimination of phages for lactobacilli, streptococcal, and lactococcal phages (Quiberoni et al., 2003;Muller-Merbach et al., 2005a;Atamer and Hinrichs, 2009), which have indicated that the level of resistance/sensitivity varies between phage species, and that certain phages appear to be quite resilient to such treatments (Avsaroglu et al., 2007).

THERMAL INACTIVATION
In the dairy industry a number of temperature treatments are applied to the incoming raw milk depending on the product. These include both low temperature long time (63 • C for 30 min) and high temperature short time (72 • C for 15 s) for Cheddar cheese, while either a treatment of 80 • C for 30 min or 95 • C for 10 min is used for yogurt production (Soukoulis et al., 2007). These temperature treatments are applied to achieve adequate pasteurization of the raw milk, however, most of these are less harsh than the temperature exposure of 90 • C for 15 min to inactivate phage recommended by the International Dairy Federation to achieve complete phage inactivation (Svensson and Christiansso, 1991).
Raw milk is a constant source of phages (McIntyre et al., 1991) that threaten fermentation processes in a dairy plant due to the ability of a certain percentage of the phage population to remain viable following pasteurization, as was confirmed by Madera et al. (2004). This study further demonstrated that the 936-type phages were significantly more resistant to heat treatments as compared to the other two predominant lactococcal phages (i.e., those that belong to the P335 and c2-type species). The ability of 936-type phage to withstand pasteurization was also observed for phages P001 and P008, with active phage particles still detectable following a 10 min and 60 min exposure, respectively (Muller-Merbach et al., 2005b). The assumption of raw milk being an important source for phage contamination was further supported by the ability of phage P008 to remain detectable at 80 • C (Muller-Merbach et al., 2005a). Additionally, this study demonstrated that milk in fact protects against thermal inactivation of phages, an observation corroborated by several other studies (Quiberoni et al., 2003;. In addition to surviving pasteurization, several lactococcal 936-type phages have emerged with remarkably higher thermal resistance. Atamer and Hinrichs (2009) reported the isolation of twenty-two 936-type phages that exhibited the ability to survive a temperature exposure of 80 • C for 5 min. P335-and c2-type phages were also isolated in this study, however, only one of the P335-type phages and none of the c2-type phages were shown to be resistant to 80 • C. Active phage particles for four of the 936-type phages, www.frontiersin.org i.e., P680, P656, P1532, and P4565, were also observed after a 5 min exposure to 85 and 90 • C. Further characterization of P680 showed that a 6-log reduction in infectivity was achieved following exposure to 90 • C for 15 min (Atamer and Hinrichs, 2010). These studies highlight the emergence of highly thermo-tolerant phage within the 936-type phage species, and, along with the protective nature of milk, it may account for the dominance of this particular species in dairy industries globally.

BIOCIDAL INACTIVATION
A range of strategies are utilized within a dairy setting to control phage populations, one such approach being the use of biocidal agents. These include sodium hypochlorite, peracetic acid, and ethanol, and the impact of these agents on bacteriophage infectivity has been recently reviewed (Guglielmotti et al., 2011). Despite extensive characterization of the biocidal inactivation of phage infecting lactic acid bacteria (LAB) species other than L. lactis (Quiberoni et al., 1999;Quiberoni et al., 2003;Capra et al., 2004), only a limited number of investigations have been carried out to study inactivation of 936-type phages using such antimicrobial agents. One study showed that the 936-type phage QF12 was completely inactivated at a concentration of 200 parts per million (ppm) sodium hypochlorite for 45 min and a T 99 value (the time required to inactivate 99% of phage infectivity) of 1.9 min and 3.8 min in 100% ethanol and isopropanol, respectively (Suarez and Reinheimer, 2002). Peracetic acid, a biocidal agent with a low pH, is known to be effective in reducing infectivity of bacteriophages, including those belonging to the 936-type species (Suarez and Reinheimer, 2002). Other investigations into the effectiveness of biocidal inactivation of lactococcal phages have reported that the impact of such chemicals can be biocide-and phagedependent, however, the species identity of the tested phages was not stated (Parada and Fabrizio, 2001;Avsaroglu et al., 2007;Buzrula et al., 2007).
A more recently described approach with limited research carried out to date is the use of TiO 2 to induce a photolytic reaction generating highly oxidizing species when exposed to ultra-violet (UV) light (Briggiler Marcóa et al., 2011). Two 936-type phages, CHD and QF9, were shown to be completely eliminated within 120 and 60 min of exposure, respectively. While this would be considered a rather lengthy procedure the study also suggested that levels required have no associated negative impact on human health. Therefore, and with the right design, equipment can be incorporated in a factory setting to constantly expose the processing plant to photocatalytic activity, so as to reduce the overall number phages while also minimizing the spread of phages between surfaces, by staff and through the air.

PHAGES RESPONSE TO PHAGE-RESISTANCE SYSTEMS
In addition to physical and chemical barriers that phages encounter in fermentation facilities, starter cultures themselves present many biological hurdles to phage infection in the form of phage-resistance systems. During the past three decades, a very significant volume of data relating to the response of various lactococcal phages to phage-resistance systems has emerged. Many of these studies also highlight that while significant levels of protection against these phages are afforded by these systems, phage variants may arise that can bypass these systems and therefore have the potential to persist in dairy fermentation facilities (Gabs and Josephsen, 2003;Fortier et al., 2005;Haaber et al., 2008). Abortive infection (Abi) systems and restriction/modification (R/M) systems are widespread in the genomes and plasmids of lactococcal strains, and have been studied extensively with a view to improving dairy starter cultures (Gabs and Josephsen, 2003;O'Driscoll et al., 2004;Chopin et al., 2005). These systems dominated studies of lactococcal phage-resistance in the 1990s and the first decade of the current century (Garvey et al., 1995;Daly et al., 1996), and consequently there is a wealth of information relating to the response of 936-type phages to such resistance systems. For example, all Abi systems identified to date provide resistance against 936-type phages with the exception of AbiR and AbiZ, for which the effects have not been well defined (Twomey et al., 2000;Chopin et al., 2005;Durmaz and Klaenhammer, 2007). Several studies of R/M systems also detail the sensitivity of 936-type phages to such systems (Twomey et al., 2000;Boucher et al., 2001;Kong and Josephsen, 2002;O'Driscoll et al., 2004). Similarly, superinfection exclusion (Sie) systems encoded by the genomes of lactococcal prophages have been tested against a variety of 936 phages and are effective against members of this phage species (McGrath et al., 2002;Mahony et al., 2008). Interestingly, a prophage-encoded Sie system from a Streptococcus thermophilus-infecting phage is also effective against the 936-type phages, when applied in a lactococcal host background (Sun et al., 2006). This unexpected (as S. thermophilus is not known to be sensitive to 936-type phages) finding suggests that this gene was acquired by horizontal transfer from a lactococcal starter and may be reminiscent of the complex and symbiotic evolutionary relationships between these starter bacteria.
While the 936-type phages appear to be sensitive to the majority of lactococcal anti-phage strategies, it is apparent that variant phages, which have overcome the imposed resistance system, are readily isolated in many cases. Possibly the best described example of this is the mutation of the sak gene of the 936 phage, p2, to overcome the Abi system AbiK (Boucher et al., 2000;Bouchard and Moineau, 2004). Similarly, genes that are involved in sensitivity to the abortive infection AbiD1 have been identified in another 936type phage, bIL66 (Bidnenko et al., 1995). The deduced product of one of these genes, Orf3, displays similarity to a putative Holliday junction endonuclease of lactococcal phages and has homologs on the genomes of many 936 phages at the rightward end of the genome within the so-called "middle-expressed" module. The relative ease with which such mutations are acquired highlights the powerful adaptive responses of these phages to biological hurdles.

PHAGE GENOMICS: EVOLUTIONARY ASPECTS
To date, complete genomes of forty-one 936-type phages have been sequenced and are publicly available (Chandry et al., 1997;Crutz-Le Coq et al., 2002;Mahony et al., 2006;Hejnowicz et al., 2009;Rousseau and Moineau, 2009a;Castro-Nallar et al., 2012). The first sequence to be published was that of sk1 and in the 15 years ensuing this, a wealth of sequence data has become available permitting extensive comparative analyses within this species (Chandry et al., 1997). The genome sequence of bIL170 provided molecular proof for early suggestions of a bipartite sub-classification based on the Frontiers in Microbiology | Food Microbiology phage's ability to infect either L. lactis ssp. lactis or L. lactis ssp. cremoris (Crutz-Le Coq et al., 2002;Dupont et al., 2004). Several comparative genomic analyses have since reported that there are a number of genomic regions of divergence, most particularly within their presumed replication module (Mahony et al., 2006;Rousseau and Moineau, 2009a;Castro-Nallar et al., 2012).
The genomes of the 936-type phages are generally between 26 and 32 kb in length, and contain up to 65 open reading frames organized into three transcriptional units: the first in a rightward orientation encoding the morphogenesis module, which is head-to-head with the second transcriptional unit. The latter is involved in DNA replication and is followed by the third unit, for which a function is not established yet and which is in a back-to-back orientation with the replication module (Figure 1; Chandry et al., 1997). These three transcriptional units are conserved in all sequenced 936 phages with point mutations and insertions/deletions being the main source of divergence between such phages.
As mentioned above, the early expressed replication module is where the majority of divergence is observed between the sequenced phages, however, many of the predicted genes have no functional assignments with the exception of a putative DNA polymerase subunit and a single-stranded binding protein (Mahony et al., 2006;Rousseau and Moineau, 2009a). This makes the replication module essentially a black box, although some elements within this module have been subjected to detailed analysis such as the sak gene of p2 (Boucher et al., 2000;Bouchard and Moineau, 2004). This is in stark contrast to the well-studied P335 phages, for which many of the replication genes have been assigned functions (McGrath et al., 1999;Labrie et al., 2008).
The most intensely characterized of the three functional modules is the morphogenesis module and predicted functions have been assigned to most of the genes in this operon either through sequence relatedness to previously characterized phages or functional assessments. Within the structural module (the late-expressed transcriptional unit) genes involved in packaging (large and small terminase subunits), the capsid, portal protein, neck passage structure, tail tape measure protein and its associated chaperone, tail structural genes and the tail tip (RBP-encoding gene) and genes involved in lysis have been identified (Figure 1; Dupont et al., 2004;Mahony et al., 2006;Siponen et al., 2009a). It is the divergence of individual components of this module that has formed the basis for much of the more recent in-depth subclassification of the 936 phages (Mahony et al., 2006;Rousseau and Moineau, 2009a;Castro-Nallar et al., 2012). Of the abovementioned gene products, the presence of a so-called"neck passage structure" and the divergence of genes encoding the tail tape measure, the RBP and the lysin have been implicated as the major points of diversification within the structural module of these phages (Figure 1; Mahony et al., 2006).
The neck passage structure is encoded by gene l12 of bIL170, which is located between the genes specifying the head and tail structural elements, and has representative homologs on many of the sequenced 936 phages (Crutz-Le Coq et al., 2002). It bears similarity to a host recognition domain identified in the temperate P335-type phages bIL309 and BK5-T. It is a structural element that is not required for propagation or tail assembly but is proposed to be involved in host range extension or determination (Crutz-Le Coq et al., 2006). This gene, or clear homologs thereof, has been identified in the genomes of bIL170, P008, 712, bIBB29, CB13, CB14, CB19, and CB20, as well as in 16 out of 28 phages isolated and sequenced in a recently published phylogeographical study (Mahony et al., 2006;Rousseau and Moineau, 2009a;Castro-Nallar et al., 2012). While this gene product is not essential, it has been suggested that its presence may provide a competitive advantage to phages possessing it (Crutz-Le Coq et al., 2006;Mahony et al., 2006). The presence and evolution of such genes provide one of the primary indicators of the environmental pressure imposed on these phages to adapt to their surroundings; the bacterial hosts available to the phages, the seemingly FIGURE 1 | Schematic representation of the typical genome architecture of a typical 936 phage. Genes with a predicted or proven function are assigned colors: green arrows indicate genes encoding the terminases and portal protein involved in phage DNA packaging; red arrows indicate the capsid genes; Blue arrows indicate tail structural genes; the purple arrow represents the receptor binding-encoding gene; pink arrows represent the holin and lysin of the lysis module; the black arrow represents a (subunit of) DNA polymerase. The transcriptional units as shown for sk1 are highlighted above the genomic drawing as are the functional modules of the genome. Below the genomic drawing are yellow highlighters of the dominant points of divergence in the genomes of the 936 phages: the neck passage structure, the tail tape measure, the receptor binding protein, lysin, and the origin of replication, notwithstanding that a general divergence is observed in the replication module in general. The indicated scale at the bottom of the figure represents base numbering starting at the cos site.
www.frontiersin.org inexhaustible number of host-encoded phage-resistance mechanisms and harsh thermal and chemical treatments which phages have to cope with in modern dairy facilities. The presence or absence of this gene was one of the first points upon which phylogenetic sub-grouping of the 936 phages was based as those phages that were known at that time to infect L. lactis ssp. lactis strains possessed a homolog of this gene (bIL170, P008), while those known to predominantly infect L. lactis ssp. cremoris strains did not (jj50, sk1; Mahony et al., 2006). Additionally, it was proposed that there may a further sub-group that is intermediate in terms of host range and these (represented by 712) possessed a truncated derivative of the gene relative to bIL170 (Mahony et al., 2006).
More significantly, there has been a tremendous surge in data relating to the identification, structural analysis and phylogenetics of RBPs of the 936-type phages (Dupont et al., 2004;Ricagno et al., 2006;Spinelli et al., 2006;Tremblay et al., 2006;Sciara et al., 2010;Castro-Nallar et al., 2012). This increased interest in the phylogeny of the 936-type RBPs has permitted a more in-depth analysis of the relationship between host range and the phylogeny of the RBPs. While it was initially suggested (Mahony et al., 2006) that there may be a tripartite subgrouping of the phylogeny of the 936-type phages based on sequence divergence of particular genomic regions, including the gene encoding the RBP, it now appears that this may be much more complex as additional genome sequences become available (Figure 2; Mahony et al., 2006). The RBP sequences are conserved at the N-terminus of the encoded proteins, while the C-terminus is variable reflecting the varied interactions that the phages has with carbohydrates on the cell surface of the lactococcal host (Siponen et al., 2009b). While there are many 936-type phages that seemingly infect either L. lactis ssp. lactis or cremoris strains, it is also apparent that many of the phages within this species are capable of crossing the subspecies barrier and many more whose host range is as yet undefined but may potentially be predicted based on the RBP phylogeny (Figure 2). This may implicate two scenarios relating to the RBP sequence: (1) the RBP C-terminus is capable of interacting with multiple carbohydrates or (2) the RBP C-terminus is interacting with a carbohydrate element which is not specific to a particular subspecies. The latter scenario seems the more plausible although this is currently conjecture. It is also noteworthy that while there is significant interest in the phylogeny of the RBPs of these phages, it is of little benefit to understanding the relationship between host range and RBP sequence without information on the host range of the phages. For example, the RBP phylogenetic tree presented here (Figure 2) highlights sub-groups capable of infecting either L. lactis ssp. lactis or cremoris. These sub-groups are interjected with sub-groups for which host range data is unknown or poorly defined (Figure 2). Predictions on host range could be predicted based on the phylogeny of the RBPs but this is merely speculation until biological data is available.
Sequence comparisons of the tail tape measure, lysine, neck passage structure, and RBP proteins encoded by the various 936-type phages have shown that they seem to follow parallel evolutionary pathways in that observed divergence in one of these proteins in a given phage is normally accompanied by a similar level of divergence of any of the other mentioned proteins in that same phage (Mahony et al., 2006). A similar observation is made when examining the genomic termini (cohesive ends) and the origin of replication of the phage at the rightward end of the genome Phages that have been demonstrated to infect specifically L. lactis ssp. cremoris or ssp. lactis, are highlighted and P475 which infects both subspecies is also highlighted. *Indicates phages that have been propagated on a L. lactis ssp. cremoris strain but whose host range data is unknown. The remainder (ASCC collection) have a known host range but this data is unpublished. This phylogenetic tree highlights the divergent nature of the RBP sequences and while many of the sub-groups at the top of the tree are L. lactis ssp. cremoris-infecting phages and the lower half of the tree may be dominated by lactis-infecting phages, there are many of unknown host range. (Mahony et al., 2006). This indicates a common ancestral pathway for the 936 phages with diversification in certain elements to provide new derivatives with increased competitiveness. As more 936 genome sequences become available, a more confident phylogeny of these phages will be established and thus continued monitoring of dairy plants and environmental samples will define newer populations and will provide information on how the species is adapting to overcome new anti-phage strategies so as to infect starter strain mixes in the dairy setting.

FUTURE PERSPECTIVES
The 936-type phages continue to dominate as the single most problematic lactococcal dairy phage species. Although significant data exists relating to the 936 phages, it is evident through genomic and evolutionary studies and furthermore through physical characterization studies that this species is continually evolving to overcome biological and physical hurdles pitted against them. This evokes a warning that constant monitoring of this species at the genomic level is of great significance in order to understand the phage population dynamics within any plant. The constant evolution of newer isolates with increasing resistance to thermal, chemical, and biological barriers requires rapid responses from the dairy industry by developing resistant starter strains and identifying alternative phage-resistant starter strains. This fascinating species of lactococcal phages is an excellent model to study the evolutionary processes within a closely related species of phages and highlights the relative importance of these phages in the dairy industry but also to the broader study of phage biology and evolution.