PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of frontmicrobioLink to Publisher's site
 
Front Microbiol. 2017; 8: 1754.
Published online 2017 September 12. doi:  10.3389/fmicb.2017.01754
PMCID: PMC5601072

Global Survey and Genome Exploration of Bacteriophages Infecting the Lactic Acid Bacterium Streptococcus thermophilus

Abstract

Despite the persistent and costly problem caused by (bacterio)phage predation of Streptococcus thermophilus in dairy plants, DNA sequence information relating to these phages remains limited. Genome sequencing is necessary to better understand the diversity and proliferative strategies of virulent phages. In this report, whole genome sequences of 40 distinct bacteriophages infecting S. thermophilus were analyzed for general characteristics, genomic structure and novel features. The bacteriophage genomes display a high degree of conservation within defined groupings, particularly across the structural modules. Supporting this observation, four novel members of a recently discovered third group of S. thermophilus phages (termed the 5093 group) were found to be conserved relative to both phage 5093 and to each other. Replication modules of S. thermophilus phages generally fall within two main groups, while such phage genomes typically encode one putative transcriptional regulator. Such features are indicative of widespread functional synteny across genetically distinct phage groups. Phage genomes also display nucleotide divergence between groups, and between individual phages of the same group (within replication modules and at the 3′ end of the lysis module)—through various insertions and/or deletions. A previously described multiplex PCR phage detection system was updated to reflect current knowledge on S. thermophilus phages. Furthermore, the structural protein complement as well as the antireceptor (responsible for the initial attachment of the phage to the host cell) of a representative of the 5093 group was defined. Our data more than triples the currently available genomic information on S. thermophilus phages, being of significant value to the dairy industry, where genetic knowledge of lytic phages is crucial for phage detection and monitoring purposes. In particular, the updated PCR detection methodology for S. thermophilus phages is highly useful in monitoring particular phage group(s) present in a given whey sample. Studies of this nature therefore not only provide information on the prevalence and associated threat of known S. thermophilus phages, but may also uncover newly emerging and genomically distinct phages infecting this dairy starter bacterium.

Keywords: panvirome, methyltransferase, transcriptional regulator, structural proteome, antireceptor

Introduction

The problem of phage predation of Streptococcus thermophilus in the dairy industry has been well described (Caldwell et al., 1996; Bruttin et al., 1997; Quiberoni et al., 2006; Garneau and Moineau, 2011), though the precise impact on the fermentation process can only be estimated due to logistical limitations and commercial sensitivities. A crucial first step in tackling this problem is the availability of comprehensive genetic information on both host and virus. The advent of the “genomics age,” as facilitated through advanced sequencing technologies (reviewed in the context of dairy starter selection by Kelleher et al., 2015), has enabled the accumulation of genetic data on many (groups of) phages, not least those infecting the most commonly employed dairy starter bacterium, Lactococcus lactis (Fortier et al., 2006; Dupuis and Moineau, 2010; Murphy et al., 2014; Mahony et al., 2015). Despite having access to full genomic data sets of both S. thermophilus hosts (Goh et al., 2011; Sun et al., 2011; Wu et al., 2014) and their infecting phages (Guglielmotti et al., 2009b; Mills et al., 2011; Ali et al., 2014), the amount of publicly released genetic data relating to lytic phages of S. thermophilus remains rather limited.

To date, the complete genomes of 20 phages infecting S. thermophilus have been published: O1205 (Stanley et al., 1997), Sfi19 and Sfi21 (Desiere et al., 1998), DT1 (Tremblay and Moineau, 1999), Sfi11 (Lucchini et al., 1999), 7201 (Stanley et al., 2000), 2972 (Levesque et al., 2005), 858 (Deveau et al., 2008), ALQ13.2 and Abc2 (Guglielmotti et al., 2009b), 5093 (Mills et al., 2011), TP-J34L and TP-778L (Ali et al., 2014), 9871, 9872, 9873, and 9874 (McDonnell et al., 2016), and (very recently) CHPC577, CHPC926, and CHPC1151 (Szymczak et al., 2016). Their availability revealed that S. thermophilus phage genomes possess a modular structure, while it also allowed an analysis of their evolution and relatedness (Lucchini et al., 1999; Proux et al., 2002), thereby providing insights into some unusual genetic lineages (Mills et al., 2011; McDonnell et al., 2016; Szymczak et al., 2016; discussed further below). It has, furthermore, improved our understanding on how phage-host interactions cause iterative genomic changes. For example, it is now known that the S. thermophilus CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats)-Cas phage-resistance system relies on the acquisition of short genomic regions from the infecting phage (Hols et al., 2005; Barrangou et al., 2007; Deveau et al., 2008), which in turn is counteracted by the accumulation of point mutations in the phage genome (Deveau et al., 2008), leading to iterative genomic alterations.

The emergence of the 5093 group (Mills et al., 2011) was a significant discovery in terms of expanding the biological diversity of S. thermophilus phages. While previously described phages (and their genomes) could be assigned to either the so-called cos-containing or pac-containing groups (Le Marrec et al., 1997), the genome of phage 5093 instead displays a striking sequence similarity to phages that infect non-dairy streptococcal species. Furthermore, the apparent absence of an antireceptor (the phage tail tip component involved in the initial adsorption to the host cell) suggests that this phage is fundamentally different in its adsorption mechanism compared to previously identified S. thermophilus phages (Duplessis and Moineau, 2001). The expanding genetic diversity of S. thermophilus phages is further highlighted by the recent discovery of the so-called 987 group phages (McDonnell et al., 2016). Members of this novel phage group appear to have evolved through a genetic exchange event between an S. thermophilus phage and an unknown member of the P335 group of L. lactis phages (Labrie et al., 2008). Recently, the genome sequences of further phages, exhibiting similarity to 5093, as well as two additional phages exhibiting similarity to P335 group L. lactis phages, have been determined (Szymczak et al., 2016).

It is clear that there is a need to expand the amount of genetic information available in this field, particularly in light of the emergence of genetically divergent phage groups described above. Here, we present whole genome sequences of 40 phages infecting S. thermophilus. This sequencing project triples the number of available full genomic sequences of S. thermophilus phages. The sequenced genomes were assessed for general structural characteristics as well as modular configurations. A multiplex PCR detection and classification system described previously was updated to include recently discovered phage groups. Novel features which may confer a selective advantage to the phage in terms of infection or persistence capability are indicated. In addition, the structural protein complement and receptor-binding protein (termed the “antireceptor”) of a representative 5093 group phage was defined, significantly updating the genomic annotation of this phage group.

Materials and methods

Phage propagation, enumeration, and storage

Streptococcus thermophilus strains were routinely grown from single colonies or from 10% Reconstituted Skimmed Milk (RSM) stocks overnight at 42°C in M17 Broth (Oxoid, Hampshire, U.K.) supplemented with 0.5% lactose (Sigma-Aldrich, St. Louis, MO, U.S.A.; LM17). L. lactis NZ9000 was maintained at 30°C in M17 broth containing 0.5% glucose (Sigma-Aldrich; GM17), while NZ9000 derivatives containing pNZ8048-based constructs were maintained in GM17 with the addition of 5 μg/ml chloramphenicol (Sigma-Aldrich; GM17 + Cm5). Phage enumeration was performed using a standard method (Lillehaug, 1997), for which LM17 broth was supplemented with 0.25% glycine (Sigma-Aldrich), 10 mM CaCl2 (Sigma-Aldrich) and either 10 g/L (solid agar base) or 4 g/L (semi-solid overlay) Technical Agar (Merck, Darmstadt, Germany). Whey samples from dairy plants producing fermented milk products were obtained and analyzed for the presence of phages against S. thermophilus using spot and plaque assay methods as described above. Single plaque isolates were then propagated as follows: 10 ml LM17 broth was inoculated with 100 μl of the appropriate host strain and grown at 42°C for 1.5–2.0 h. Using a sterile pipette tip, a single, well defined plaque was added to the growing culture, thoroughly mixed and incubated for a further 2–4 h or overnight. The lysed culture was centrifuged and the supernatant filtered (0.45 μm). The filtered supernatant was stored at 4°C and used as stock for subsequent assays. This process was repeated at least twice prior to DNA sample preparation and subsequent genome sequencing.

Phage purification

Individual phages (Table (Table1)1) were propagated as described above in a 2 L volume before the addition of (poly)ethylene glycol (Sigma-Aldrich) to a final volume of 10% (w/v) and NaCl (Sigma-Aldrich) to a final concentration of 0.5 M. The mixtures were incubated at 4°C for at least 6 h to encourage precipitation before centrifugation at 17,700 × g for 15 min (to concentrate) and resuspension in 5 ml TBT Buffer (100 mM NaCl, 100 mM Tris-HCl (pH 7), 10 mM MgCl2, 20 mM CaCl2; Sigma-Aldrich). The suspension was extracted at least twice using an equal volume of chloroform (Fisher Scientific, Waltham, MA, U.S.A.) and phages were purified by a discontinuous (3 M/5 M) cesium chloride (Sigma-Aldrich) gradient centrifugation at 76,000 × g for 2.5 h. Translucent blue bands visible at the interface of the gradient after centrifugation were carefully removed using a syringe and dialyzed against 50 ml TBT overnight at 4°C. Phage preparations were stored at 4°C until required for electron microscopy and DNA extraction.

Table 1
Characteristics of bacteriophage genomes sequenced during this study.

DNA preparation and restriction profile analyses

Phage DNA was prepared using existing methods (Sambrook et al., 1989; Moineau et al., 1994), and as described previously (McDonnell et al., 2016). At least 5 μg DNA was extracted using this method and quantified using a NanoDrop 2000 (Thermo Scientific). Total genomic DNA was qualitatively analyzed by agarose (1%; Sigma-Aldrich) gel electrophoresis, on which the samples were separated at 100 V for at least 20 min.

Unique phages were distinguished (prior to DNA sequencing) by restriction profile analysis using EcoRV or HaeIII FastDigest (Thermo Scientific) enzymes, whereby a reaction mixture containing 8 μl phage DNA, 1 μl enzyme, 2 μl 10X FastDigest Buffer and 9 μl sterile distilled water (sdH2O) was incubated at 37°C for 17 min prior to electrophoretic examination as described above.

Multiplex PCR detection of phages

PCR multiplex detection of S. thermophilus phage groups was performed using phage DNA as template and is a modification of a previously described method (Quiberoni et al., 2006). Two primer sets, each with a target in the 987 and 5093 group phages, were added to the PCR reaction in addition to the two previously described primer sets (targeting cos- and pac-containing phages; Table Table2).2). PCR reactions were conducted in a total volume of 25 μl using Taq polymerase (Qiagen, Hilden, Germany). Reaction conditions were as follows: initial denaturation of 95°C for 2 min, 30 cycles of 95°C for 15 s, 55°C for 30 s, and 72°C for 1 min, with a final extension of 72°C for 10 min.

Table 2
PCR primers applied during this study.

DNA sequencing and in silico analysis

Phage DNA samples exhibiting unique restriction pattern profiles were sequenced using one of two methods, as follows: (i) random shotgun sequencing using pyrosequencing technology was performed (Macrogen, Inc., Geumcheon-gu, Seoul, South Korea) using a 454 FLX instrument yielding at least 65-fold coverage for each genome, with the exception of phage 7632, which yielded 17-fold coverage. Individual sequence files were assembled using GSassembler (454 Lifesciences, Branford, CT, U.S.A.), generating a consensus sequence for each phage; (ii) The remaining genomes were determined by GenProbio srl (University of Parma, Italy) using a MiSeq sequencer (Illumina, USA). Genomic libraries were constructed employing the TruSeq Nano DNA LT Kit (Illumina) and using 200 ng of genomic DNA, which was fragmented with a Bioruptor NGS ultrasonicator (Diagenode, USA) followed by size evaluation using Tape Station 2200 (Agilent Technologies). Library samples were loaded into a Flow Cell V3 600 cycles (Illumina) according to the technical support guide. The generated paired-end reads (2 × 251 bp read lengths) were depleted of adapter sequences, quality filtered and assembled through the MIRA software program (version 4.0.2; Chevreux et al., 2004) and the MEGAnnotator pipeline (which was also used to perform the initial gene annotation; Lugli et al., 2016), where a coverage level of at least 158-fold was achieved. For all genomes, remaining gaps were closed and quality improvement of the genomes was carried out (particularly across homopolymeric tracts) by PCR and additional Sanger sequencing (performed by Eurofins MWG, Ebersberg, Germany). ORFs on each genome, representing putative protein products, were predicted using a Heuristic approach (Genemark; Besemer and Borodovsky, 1999) and annotated using the Basic Local Alignment Search Tool (Altschul et al., 1990), as well as Pfam (Finn et al., 2014), and HHpred (Soding et al., 2005). Amino acid identity between phage-encoded proteins was determined using BLASTP. Percent identity and divergence tables (Supplementary Figure S1) were constructed using Megalign [DNASTAR; version 7.1.0 (44)]. Phage VR2 region genetic relatedness trees were generated using the MEGA program (Tamura et al., 2013). Pan-genome analysis of a total of 57 S. thermophilus phage genomes was performed using PGAP v1.0 (Zhao et al., 2012), which employs the Heaps law pan-genome model (Tettelin et al., 2005). The pan-genome profile was built by firstly organizing the ORF content of each genome into functional gene clusters using the Gene Family method (Zhao et al., 2012).

Electron microscopic analysis

Cesium chloride purified phage samples were prepared as outlined above and subjected to electron microscopic analysis as described previously (Vegge et al., 2005; Casey et al., 2015).

Structural protein identification

Methanol-chloroform extraction of phage proteins, SDS-PAGE visualization and preparation of phage structural protein samples for mass spectrometry were performed as described previously (Casey et al., 2015). Electrospray ionization-tandem mass spectrometry (ESI-MS/MS) was performed as previously described (Ceyssens et al., 2008; Vanheel et al., 2012). For each peptide identified, coverage levels of two unique amino acid strings, or at least 5% of the total protein, were used as cut-off values when identifying gene products as components of the viral particle (Cornelissen et al., 2012).

Protein purification and adsorption inhibition assay

Proteins were purified using a method adapted from a previously described methodology (Collins et al., 2013). Firstly, ORF180095 (encoding RBP0095), and its attached 6xHis N-terminal-specifying purification tag and appropriate restriction enzyme recognition sequence was PCR amplified (using the primers described in Table Table2)2) employing Phusion polymerase (New England Biolabs, Ipswich, MA, U.S.A.) and cloned behind the Nisin-inducible promoter of lactococcal plasmid pNZ8048 (Kuipers et al., 1998). Plasmid DNA was dialysed against sdH2O for 10 min and transformed into electrocompetent L. lactis NZ9000 cells (Kuipers et al., 1998). Plasmid DNA was then extracted using the GeneJet Plasmid Miniprep Kit (Thermo Scientific) and subjected to Sanger sequencing (as above) to confirm the expected sequence of the recombinant plasmid, which was termed pNZ8048+RBP0095. NZ9000 containing plasmid pNZ8048+RBP0095 were grown to an OD600nm of 0.2 before nisin induction [10 ng/ml; using Nisaplin (Danisco, Copenhagen, Denmark)], cell lysis, and sonication as previously described (Collins et al., 2013). Target protein purification was then performed using a Ni-nitrilotriacetic acid (NTA) agarose (Qiagen) column (Bio-Rad, Hercules, CA, U.S.A.) as previously described (Collins et al., 2013). Protein fractions were eluted using varying concentrations of imidazole buffer (as per manufacturer's instructions) and separated on a 12.5% SDS-PAGE gel at 160 V for 80 min. Fractions containing bands of the correct size with minimal contamination were dialysed against 100 ml protein buffer (as above) three times for 40 min to remove remaining imidazole. Dialysed fractions were stored at 4°C for use in subsequent assays.

Adsorption inhibition assays were performed using a method adapted from two previously reported methods (Garvey et al., 1996; Collins et al., 2013). Briefly, the S. thermophilus strain to be tested was grown to at least OD600nm 0.5 (yet not above 0.54) and resuspended in 225 μl 14 strength Ringer's solution (Merck). 50 μl of varying concentrations of purified antireceptor protein (or protein buffer control) was added to cells (or 14 strength Ringer's solution control) and incubated for 1 h at 42°C. Phage lysates (225 μl), diluted in 14 strength Ringer's solution to a concentration of approximately 1 × 105 pfu/ml were added to the cells or control and the mixture was incubated at 42°C for a further 12 min before centrifugation at 20,000 × g for 3 min to pellet bound phages with the cells. Remaining phages were quantified by plaque assay (as described above). Adsorption to wild-type and antireceptor-incubated cells was calculated as follows: (control phage titer − sample phage titer)/(control phage titer) × 100. Adsorption inhibition, expressed as a percentage of phage adsorption to wild type cells, was calculated as follows: (% adsorption on WT − % adsorption on incubated cells)/(% adsorption on WT) × 100.

Results and discussion

Isolation of phages, genotype grouping, selection for sequencing and PCR detection

A total of 40 phage isolates (Table (Table1)1) were genetically characterized as part of this study. These isolates represent a subset of a larger phage collection of 120 phages isolated from whey samples which were obtained from dairy processing plants in various geographical locations from different continents (including Europe, North-America and Asia) and at various time points (being collected in the period 2006–2012, inclusive). Plaque assays were performed on suitable host strains, and single plaques were propagated as described in the Materials and Methods section. The resulting high titre, purified phage preparations were then employed in a phage-host survey, using a total of 91 potential host strains, in order to determine their host ranges. In the majority of cases, isolated phages infected their primary (i.e., the host on which they were originally propagated) host and one further strain—with a small number of phages infecting between three and six S. thermophilus strains—thus overall revealing a narrow host range (data not shown), which is consistent with previous observations for phages infecting this species (Binetti et al., 2005; Zinno et al., 2010). Following this analysis, and in order to assess the distribution of phage groups in the collection, phage isolates were subjected to multiplex polymerase chain reaction (PCR) grouping using the primer sets described by Quiberoni et al. (2006), which are designed to detect cos- and pac-site containing phage groups. While the vast majority of phages could be grouped using this method, eight isolates did not yield an amplicon using this PCR system employing either fresh lysate or extracted DNA as template material (discussed further below).

Three main criteria were used to select phages for sequencing: (i) phages which did not return a PCR product using the above mentioned phage grouping PCR and which were therefore considered novel, (ii) persistence and/or prevalence (as determined by similar phages being isolated persistently from the same dairy factories, or phages infecting two or more distinct S. thermophilus strains), and (iii) diversity (as determined by host range and restriction profiling of phage genomic DNA). For the purpose of this study, the results (and discussion thereof) below pertain to the subset of 40 phages whose genomes were sequenced, unless otherwise indicated.

As indicated in Table Table1,1, cos-containing phages were the most frequently detected phage group (both in the 120-phage collection as a whole, as well as within the selected 40 phages), a dominance also observed—though to varying degrees—in previously analyzed phage collections (Zago et al., 2003; Guglielmotti et al., 2009a; Zinno et al., 2010). The reason for the observed dominance of cos-containing phages is not known.

As noted above, a total of eight phage isolates did not yield a PCR product using the standard cos- and pac-containing phage detection PCR. Genome sequencing was performed on these eight isolates, which revealed that four of these represent a novel group, designated the 987 group (these findings have been published elsewhere; McDonnell et al., 2016), with the remainder belonging to the 5093 group (and described below). These genome sequences were utilized to design a multiplex PCR primer set (Table (Table2),2), which allows detection of representatives of all four S. thermophilus-infecting phage groups (Figure (Figure1,1, Lanes 5–12). The two newly designed primer sets were successfully tested on all four members of both the 5093 and 987 phage groups isolated in this study. Furthermore, considering that the primers used to detect the 987 (targeting Open Reading Frame 6 [ORF6] of phage 9871, which encodes a conserved scaffolding protein) and 5093 (targeting ORF3 of phage 0095, which encodes a portal protein) phage groups produce products of different sizes (707 base pairs [bp] and 983 bp, respectively), these newly designed primer sets are fully compatible with each other and with the cos and pac phage-detecting pairs (Quiberoni et al., 2006). This multiplex PCR method thus embodies a very useful updated tool to detect representatives of all currently known phage groups infecting this starter species.

Figure 1
Electrophoretic examination of PCR products amplified from phage DNA (as template) using the multiplex PCR phage detection primers described in Table Table2.2. Lane 1: 1 kb Full scale DNA Ladder (Fisher Scientific), Lane 2: 7201 (cos-containing ...

Morphological analysis

All S. thermophilus phages identified to date are classified as Siphoviridae, corresponding to type B as described by Bradley (1967), and thus possess icosahedral heads and non-contractile tails. Representative members of the phage collection employed in this study were selected for electron microscopic (EM) analysis, with members of the cos-containing, pac-containing and 5093 groups shown in Figures 2A–C, respectively), while EM analysis of members of the 987 group has been presented elsewhere (McDonnell et al., 2016). All analyzed phages display the expected morphology, particularly in terms of tail length, possessing tails that are longer than those of (the majority of) their lactococcal counterparts, i.e., 936, c2 and many P335 phages (Pedersen et al., 2000). Phages 7573 and 7951 exhibit an extended, thin tail fiber (Figures 2A,B; indicated by arrows), although genome analysis failed to identify candidate genes which may encode these interesting structures. Phage 0092 is a member of the 5093 group, whose morphology differs from the cos- and pac-containing group phages in that the tail tips are characterized by globular features that protrude from the base of the tail (Figure (Figure2C;2C; indicated by an arrow).

Figure 2
Electron microscopic analysis of representative cos-containing (phage 7573; A), pac-containing (7951; B) and 5093 group (phage 0092; C) phages visualized using uranyl acetate staining.

Phage genome features

The general characteristics of the sequenced phage genomes are given in Table Table1,1, with predicted ORF schematics of representative phages displayed in Figure Figure3.3. An overall comparative analysis of all 40 genomes (based on percent nucleotide identity as well as divergence score) is provided in Supplementary Figure S1. The complete nucleotide sequences of a total of 26 cos-containing phages, 10 pac-containing phage and four 5093 group phages were determined, with the genomes of the four 987 group phages presented elsewhere (McDonnell et al., 2016). Genomes ranged in size from 33.3 to 38.1 kilobasepairs (kb) (the longest and shortest genomes sequenced are highlighted in bold in Table Table1).1). The number of predicted ORFs in each genome ranged from 42 to 55 (increasing ORF number generally corresponding with increasing genome size), with GC content averaging 38.95%—consistent with the GC% content of previously published S. thermophilus genomes (Bolotin et al., 2004; Treu et al., 2014; Wu et al., 2014; Labrie et al., 2015). The 40 sequenced phage genomes are divided into three groups, based on amino acid identity of predicted ORF products to those of previously defined cos-containing, pac-containing and 5093 group phages (Le Marrec et al., 1997; Mills et al., 2011), which is in agreement with the results of the PCR-based typing system (discussed above).

Figure 3
Comparative analysis of the genetic organization and content of sequenced phage genomes. (A) Phages infecting S. thermophilus strain ST68757. (B) Phages infecting S. thermophilus strain ST67009. Predicted ORFs (indicated by arrows) and gene products (putative ...

The overall genomic organization, and amino acid conservation within defined groups, of the sequenced phages is comparable with that of previously published examples (Levesque et al., 2005; Guglielmotti et al., 2009b; Mills et al., 2011), indicating that this modular arrangement is highly conserved amongst these phages, even between distinct phage groups across a large number of samples. The structural protein-encoding genes of representatives of these phage groups have been defined previously (Levesque et al., 2005; Duplessis et al., 2006; Guglielmotti et al., 2009b; Mills et al., 2011; Szymczak et al., 2016), with the known structural protein-encoding genes in the 5093 group being expanded upon here. The functional annotation of the genes encoding tail components of S. thermophilus virions (specifically those encoding the distal tail (Dit), tail-associated lysin (TAL), antireceptor and baseplate components) was carefully performed in this study, in particular for members of the 5093 group, and also exhibit functional synteny between phage groups despite substantial nucleotide divergence (Figure (Figure3;3; McDonnell et al., 2016).

Figure Figure3A3A details the predicted ORFs of five phages infecting S. thermophilus ST67009 (one cos-containing phage and four 5093 group phages), confirming the genetic conservation of members of the 5093 lineage. The amino acid divergence between proteins encoded by the structural regions of this and other groups, previously observed by Mills et al. (2011), is evident. This is despite the conservation of the 7201-like/Group II replication module (discussed below) in all five phages.

Streptococcus thermophilus phage replication modules are defined here as the genomic regions immediately following the lysogeny replacement module (i.e., that region containing presumably remnant genes involved in the phage lysogenic cycle; Lucchini et al., 1999), and preceding the gene encoding the small subunit of the terminase. These regions contain genes encoding predicted primosome components, DnaC-like proteins (helicases), single-stranded DNA and DNA-binding proteins, endodeoxyribonucleases related to RusA, primases and replisome organizers, amongst various genes with hypothetical functions. Akin to the structural modules, S. thermophilus phage replication modules can be divided into (at least) two distinct groups, namely the Sfi21- and 7201-like, or Group I and II (Stanley et al., 2000; Brussow and Desiere, 2001); the former being detected most frequently in phages in the present study (present in 27 out of 40 examined phages), which is consistent with previous studies (Desiere et al., 1997; Stanley et al., 2000; Brussow and Desiere, 2001). The replication modules of two phages (7951 and 7955; Table Table1)1) do not appear to fall into either of the two above mentioned categories (based on BLAST searches), and therefore may be included in the “non-I/non-II” grouping as previously proposed by Stanley et al. (2000). In contrast with other features of these phages, such as structural protein content (Le Marrec et al., 1997), there appears to be no correlation between possession of the Sfi21- or 7201-like replication modules and mode of DNA packaging, e.g., in the cases of phages 7573 and 7574 (Figure (Figure3B).3B). This phenomenon may be due to modular rearrangements (Lucchini et al., 1999) producing a seemingly random distribution of discrete genomic segments in S. thermophilus phages. The “rightward” genomic region (defined here as the region at the 3′ end of the lysis module) is often characterized by insertions, deletions, and point mutations leading to a large degree of heterogeneity (Lucchini et al., 1999; Figure Figure3).3). Indeed, this was found to be the case in these segments of the phage genomes analyzed in the present study, which were (in some cases) observed to encode gene products of apparently non-streptococcal phage origin (e.g., RecT recombinases and transposases, indicated in the legends of Figures 3A,B).

Interestingly, several proteins predicted to be involved in (mostly host) DNA recombination are encoded by 12 phages in the present collection-either singly or in combination (Table (Table1).1). Among those identified are RecT and ERF family proteins (Hall and Kolodner, 1994; Noirot and Kolodner, 1998; Passy et al., 1999). These protein superfamilies have previously been detected in phages of low GC, Gram-positive bacteria, including S. thermophilus phage 7201 (Stanley et al., 2000) and L. lactis phage r1t (van Sinderen et al., 1996), and are generally associated with the presence of endodeoxyribonucleases of the RusA family (Sharples et al., 1994; Macmaster et al., 2006), as well as MTases and single-stranded DNA-binding proteins (Iyer et al., 2002), as is the case in a number of phages sequenced in this study (Table (Table1).1). In a similar fashion, recently described bacterial RecT-encoding genes have been associated with genes encoding putative exonuclease genes, and this “modular” combination was also observed on the genomes of four phages sequenced in the present study (Datta et al., 2008; phages 5641, 7574, 7951, and 7955; Table Table11).

The presence of HNH endonuclease-encoding genes in a 3′ “terminal” position (i.e., present downstream of predicted transcriptional regulator-encoding genes but preceding the deduced cos-sites) was found to correlate with mode of DNA packaging in each phage. Those phages utilizing the cos mode of DNA packaging possessed a terminal HNH endonuclease, whereas those predicted to utilize the pac mode did not (Table (Table1).1). This finding is consistent with previous studies examining the role of HNH endonucleases in the DNA packaging process, and specifically in cos-site cleavage (Kala et al., 2014).

Panvirome analysis

A pan-genome analysis of all phages (“pan-virome”) sequenced in this study, the 987 group phages as well as 13 previously sequenced (and publically available) S. thermophilus phage genomes was conducted as described in the Methods section, the results of which are shown in Figure Figure4.4. The exponential values of each analysis (Figures 4A–D, inset) indicate that, at least in the case of the cos-containing, pac-containing and 987 groups, sequencing of additional genomes of phages in these groups may not lead to a significant increase in known genetic diversity. However, due to the low number of complete genomes available for the 987 group phages (a total of 4 analyzed in this study), this conclusion cannot be drawn, despite the exponent value of <0.5 (Figure (Figure4D).4D). The inclusion of further genome sequences, such as those of the recently published P335-like S. thermophilus phages (Szymczak et al., 2016), in future studies, may be pertinent in order to overcome this limitation. An exponent value of >0.5 which was calculated for the 5093 group phages indicates an “open” pan genome and highlights the usefulness of whole genome sequencing of additional members of this novel group. In general, and as discussed above, conservation was observed within phage groups across the structural regions of the genomes, while the lysis and replication modules displayed an increased genetic divergence. In addition, no gene was identified with homologs across every individual member of the four phage groups (core genome = 0), indicating the high level of genetic diversity present in phages infecting S. thermophilus.

Figure 4
Panvirome analysis of the four currently known groups of phages infecting S. thermophilus. (A) cos-containing group (n = 31), (B) pac-containing group (n = 17), (C) 5093 group (n = 5), (D) 987 group (n = 4). Exponent and R2 values are given in the inset ...

Notable genomic features

Methyltransferases

The process of self-methylation of genomic DNA has been acknowledged as a protective strategy employed by phages against host-encoded restriction-modification (R/M) systems (reviewed by Murphy et al., 2013), which have been detected previously in S. thermophilus (Burrus et al., 2001; Goh et al., 2011; Labrie et al., 2015). In total, 12 predicted MTase-encoding genes were detected on 7 phage genomes sequenced as part of this study. In order to prove that phage-encoded MTases protect such S. thermophilus phages against DNA restriction, one deoxyadenosine (DAM) MTase-positive phage (indicated by the similarity of the protein product of this MTase to GATC-specific methytransferases, observed using the REBASE tool; Roberts et al., 2015), and one MTase-negative phage infecting the same strain (as control) were subjected to restriction profile analysis using DpnI and DpnII restriction endonucleases, as described previously (Murphy et al., 2014). It is known that DpnI only targets DNA which is DAM-methylated, while DAM-methylation protects DNA from restriction by DpnII (Lacks and Greenberg, 1977). Supplementary Figure S2 clearly shows that while the genome of phage 9901 (MTase-negative) is restricted by DpnII (and not by DpnI; Lanes 2 and 3), the genome of 9902 (MTase-positive) was restricted by DpnI (and not by DpnII; Lanes 4 and 5), thus validating the protective effect and functional activity of this MTase.

Predicted MTase-encoding genes were also present on the genomes of each of the four 5093 group phages {0092 to 0095 (and, interestingly, also on 5093; Mills et al., 2011)}, although these were of a different type than those discussed above, and predicted (using a BLAST search; Altschul et al., 1990) to represent cytosine-5 methyltransferases (C5 MTase; reviewed by Kumar et al., 1994). A striking feature of these particular MTases is their apparent total overlap, being produced as translation products coded by the same DNA in the same position on the phage genome, but in different reading frames (Figure (Figure3A,3A, indicated by “C”). One of these products bears amino acid similarity to a C5 MTase present in the temperate Streptococcus pneumoniae phage (Obregon et al., 2003), where a similar overlap was identified as previously also found in a conjugative streptococcal transposon (Sampath and Vijayakumar, 1998). The putative origin of these genes is consistent with the amino acid similarity of many gene products in the 5093 group phages to non-dairy streptococcal phages, a connection which has been explored previously in prophages of pathogenic streptococci (Desiere et al., 2001).

Introns

Introns which interrupt lysin-encoding genes are known to be widely distributed in S. thermophilus phages (Foley et al., 2000; Ali et al., 2014). A number of putative introns were indeed detected on the genomes of phages in the current study, though their function is unknown. In total, 19 group IA introns were identified in lysis modules, based on an interrupted lysin-encoding gene and/or the presence (with varying nucleotide identity) of a 14 bp consensus sequence correlated with the possession of an intron (Foley et al., 2000).

Transcriptional regulation

Late transcriptional regulators (Ltr's) in phages infecting Gram-positive bacteria have been implicated in DNA packaging and lysis (Quiles-Puchalt et al., 2013), and a role for Ltr encoded by certain S. thermophilus phages has previously been proposed by Ventura et al. (2002). Ali et al. (2014) observed an ArpU-like transcriptional regulator-encoding gene preceding the gene specifying the small subunit of the terminase in temperate phages TP-J34 and TP-778L, similar in organization to the Ltr-encoding gene which was recently identified as a member of a superfamily of Ltr's present in phages of Gram-positive bacteria (Quiles-Puchalt et al., 2013). The genomic position of these genes is similar to the gene encoding a late transcriptional regulator identified in L. lactis phage TP901-1 (Brondsted et al., 2001). These findings prompted us to conduct an in silico analysis of sequenced S. thermophilus phages to determine the prevalence of Ltr homologs and functional equivalents in phages of this species. All previously sequenced S. thermophilus phages (including the recently described 987 group phages), as well as those sequenced as part of this study were observed to encode an Ltr-like protein, which either belongs to the ArpU family or the Ltrb family (DUF1492; Quiles-Puchalt et al., 2013). With the exception of phage ALQ13.2 (Guglielmotti et al., 2009b) and 4761 (this study), the possession of a regulator of the Ltrb family correlates with phages utilizing the cos mode of DNA packaging, and possession of ArpU with the pac mode. This correlation is consistent with the role of the Ltr protein families in regulating the expression of the genes encoding DNA packaging machinery in other phages (Quiles-Puchalt et al., 2013).

Phage VR2 sequence clustering

The antireceptor-encoding genes of the phages sequenced in this study were generally conserved (with the exception of the 5093 group phages, whose antireceptor-encoding genes are discussed separately below), and found (using the NCBI CDD tool; Marchler-Bauer et al., 2015) to contain distinct domains, separated in some cases by collagen-like repeats as described previously (Duplessis and Moineau, 2001). S. thermophilus phage antireceptor proteins have previously been shown to harbor variable regions, one of which (VR2) is purported to be the main host range determinant (Duplessis and Moineau, 2001). An unrooted relatedness tree showing the genetic distance between the deduced amino acid sequences of these VR2 regions in the present phage collection is shown in Supplementary Figure S3 and in which VR2 regions are color-coded by host. In the majority of cases, the VR2 regions cluster according to host strain, but is not phage group-dependent, i.e., the VR2 regions of cos- and pac-containing phages are not aligned in the same cluster (as has similarly been observed; Binetti et al., 2005). The presence of outliers is not surprising, as it has been suggested that the VR2 region may not be the only host determinant operating in S. thermophilus phages (Duplessis et al., 2006). A comprehensive analysis of the multiple variable regions in the tail gene-encoding modules of S. thermophilus phages, as well as sequence and structural information on host-encoded receptors, may be required to reconcile the observed VR2 region anomaly.

5093 phage group structural protein identification

The structural protein complement of phage 0095, as a representative of the 5093 phage group, was determined by mass spectrometry (Figure (Figure5)5) as described in the Materials and Methods section. Three protein products hypothesized to form the phage tail tip and so-called “initiator complex” (Dit, TAL and antireceptor) were all detected as structural components of the 0095 virion. As expected, the presumed major capsid protein and major tail protein were also detected (Figure (Figure5),5), of which homologs had previously been identified for phage 5093 (Mills et al., 2011) and CHPC1151 (Szymczak et al., 2016). Furthermore, the predicted portal protein, the putative tape measure protein (TMP), and several minor structural proteins were identified (Figure (Figure5),5), significantly expanding the experimentally-determined identification of structural proteins for this phage group. The confirmation of the structural nature of these protein products, as well as the genomic position (and, indeed, order) of the corresponding genes confirms an overall functional synteny in this region across all four S. thermophilus-infecting phage groups (namely the cos-containing, pac-containing, 5093 and 987 groups).

Figure 5
Structural proteome analysis of phage 0095. (A) SDS-PAGE gel (12%) showing the structural protein profile of phage 0095. Lane 1: Broad range protein ladder (New England Biolabs); Lane 2: phage 0095 protein extraction. (B) Deduced structural proteins (and ...

Identification of the 5093 phage group antireceptor

The 5093 genome appears to lack an assigned antireceptor-encoding gene (Mills et al., 2011). Phages are expected to encode an antireceptor as such proteins are responsible for the initial binding of the phage to the bacterial cell surface thereby determining host specificity (Duplessis and Moineau, 2001). Assuming that the conserved modular genomic structure of S. thermophilus phages (discussed above) also applies to phage 5093, ORF295093 is predicted to encode the TMP, known to be responsible for phage tail length (Pedersen et al., 2000). This is borne out by its large size (the encoded protein is approximately 1500 amino acids in length) and the presence of extensive helical structures, as predicted by the PSIPRED tool (Jones, 1999). Downstream of ORF295093, ORF305093 is presumed to specify the distal tail protein (Dit), based on a HHPred-mediated match for this product to the Dit protein of Bacillus phage SPP1 (P = 6.2 × 10−31). The product of the next gene (ORF315093) shows similarity to phage endopeptidases using BLAST and Pfam searches, and its position on the genome is consistent with a common lactococcal phage tail structure, namely the tail-associated lysin or TAL (Kenny et al., 2004; Mc Grath et al., 2006; Stockdale et al., 2013).

ORF325093, ORF335093 and ORF345093 are predicted to encode two hypothetical proteins and a hydrolase, respectively. However, in the four 5093 group phages sequenced in the current study, these three deduced ORFs form instead a single ORF (for example ORF180095 in the case of phage 0095; Figure Figure3A),3A), the product of which is 821 amino acids in length. We propose that this single ORF encodes the antireceptor (henceforth termed RBP, including a subscript number to indicate the particular phage it is derived from, e.g., RBP0095) in this group of phages (Figure (Figure3A).3A). This hypothesis is based on: (i) the gene position at the 3′ end of the tail morphogenesis module, but preceding the lysis module, consistent with typical RBP-encoding gene positions (Duplessis and Moineau, 2001; Mahony et al., 2013; Casey et al., 2014) and (ii) the presence, at the C-terminal end of the protein, of a GDSL-family esterase domain, detected using Pfam, CDD, and HHPred. GDSL lipases are members of the diverse SGNH hydrolase superfamily, and examples include the Axe2 acetylxylan esterase, responsible for the cleavage of acetyl side chains from xylooligosaccharides, such as xylan (Lansky et al., 2014). These similarities suggest a carbohydrate binding function for this protein product. Taken together, and further corroborated below using an adsorption inhibition assay, these findings lead us to assign the antireceptor function to this gene product. It is not known if phage 5093 possesses several truncated forms of this particular ORF or if the observed differences are due to sequencing errors.

To verify its role as the antireceptor, RBP0095 was heterologously expressed and purified (Figure (Figure6A),6A), resulting in a protein with a molecular weight of approximately 90 kilodaltons (kDa), consistent with the calculated molecular weight (of 92.1 kDa). Figure Figure6B6B clearly shows that this purified antireceptor, when incubated with ST67009 cells (host for phage 0095) at varying concentrations, inhibits adsorption of phage 0095 to the host in a dose-dependent manner. Phage 0091 (a cos-containing phage also infecting ST67009) was included as a negative control to indicate the specificity of RBP0095 in inhibiting the adsorption of phage 0095 (but not that of phage 0091; Figure Figure6B).6B). The concentrations observed to exert blocking of phage adsorption were comparable to those previously observed using the 987 group antireceptor (McDonnell et al., 2016) and that of lactococcal phages Tuc2009 and TP901-1 (Collins et al., 2013). In the present study, specific concentrations of RBP0095 (Figure (Figure6B)6B) were selected to illustrate an even distribution of adsorption blocking levels. This dose-dependent reaction indicates that RBP0095 is indeed responsible for host interaction of this phage, which may be extended to the other ST67009-infecting 5093 group phages based on the high level of amino acid identity (99%) between the gene products (Figure (Figure3A).3A). Considering the lytic nature of these phages, and the fact that they appear to be emerging in commercial fermentations (having been either absent or undetected until 2011), this information may be industrially valuable, and the above data represent a significant step forward in the annotation and characterisation of these recently discovered (and potentially important) phages.

Figure 6
Phage 0095 adsorption inhibition analysis using varying concentrations of purified RBP0095 on S. thermophilus ST67009 by blocking assay. (A) SDS-PAGE gel (12%) showing purified antireceptor of phage 0095. Lane 1: Blue prestained protein standard, Broad ...

Conclusions

The genomes of phages infecting S. thermophilus, which are a recurring problem in the dairy industry, are conserved in terms of modular genomic arrangements, but divergent between defined cos-containing, pac-containing, 987 and 5093 groupings. Complete genome sequencing has enabled the development of a multiplex detection PCR capable of rapidly identifying individual members of these four phage groups, and shed light on their mechanisms of proliferation, their adaptation to their environment and to host-encoded defenses. These include - but are not limited to - genome modular rearrangements, point mutations and phage-encoded methyltransferases. Structural protein identification through mass spectrometry and antireceptor definition by means of an adsorption inhibition assay were used here to characterize a representative member of the recently discovered 5093 group, substantially improving the annotation of members of this group.

To our knowledge this study represents the largest single (published) undertaking of S. thermophilus phage whole genome sequencing. The increased database of phage genomes provides a valuable knowledge resource to the dairy industry which relies on up-to-date phage detection and phage-host interaction information in order to successfully implement rotational schemes which minimize the economic impact of phage fermentation contamination.

Declaration

The datasets generated and/or analyzed during the current study are available in the GenBank repository, https://www.ncbi.nlm.nih.gov/genbank/, under accession numbers KY705251-KY705290, inclusive.

Author contributions

BM was involved in the experimental design and work and prepared the manuscript. JM was involved in the experimental design and work and data analysis interpretation. LH performed the initial processing of whey samples and isolation of phages. HN performed the electron microscopy. JN performed the mass spectrometry. GL and MV performed the Illumina sequencing of phage genomes. TK was involved in the experimental design and manuscript editing. DS was involved in the experimental design, data interpretation and manuscript editing and preparation. All authors read and approved the manuscript.

Conflict of interest statement

A patent application describing elements of this work has been submitted. BM is funded by DSM Food Specialties, and LH and TK are employees of DSM Food Specialties. The other authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors gratefully acknowledge the technical assistance of Erick Royackers in deducing the structural protein complement of phage 0095. We also gratefully acknowledge the technical assistance of Philip Kelleher in performing the panvirome analysis and of Angela Back in performing the electron microscopic analysis. The authors gratefully acknowledge the financial support of DSM Food Specialties. JM is in receipt of a Starting Investigator Research Grant (SIRG) (Ref. No. 15/SIRG/3430) funded by Science Foundation Ireland (SFI). DvS is supported by a Principal Investigator award (Ref. No. 13/IA/1953) through Science Foundation Ireland (SFI).

Glossary

Abbreviations

CRISPR
clustered regularly interspaced short palindromic repeats
DAM
deoxyadenosine
MTase
methyltransferase-encoding gene
PCR
polymerase chain reaction
BIM
bacteriophage-insensitive mutant
bp
base pair
ORF
open reading frame
Dit
distal tail
TAL
tail-associated lysin
R/M
restriction-modification
Ltr
late transcriptional regulator
TMP
tape measure protein
RBP
receptor binding protein
kDa
kilodaltons
RSM
reconstituted skimmed milk.

Supplementary material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb.2017.01754/full#supplementary-material

References

  • Ali Y., Koberg S., Hessner S., Sun X., Rabe B., Back A., et al. . (2014). Temperate Streptococcus thermophilus phages expressing superinfection exclusion proteins of the Ltp type. Front. Microbiol. 5:98. 10.3389/fmicb.2014.00098 [PMC free article] [PubMed] [Cross Ref]
  • Altschul S. F., Gish W., Miller W., Myers E. W., Lipman D. J. (1990). Basic local alignment search tool. J. Mol. Biol. 215, 403–410. 10.1016/S0022-2836(05)80360-2 [PubMed] [Cross Ref]
  • Barrangou R., Fremaux C., Deveau H., Richards M., Boyaval P., Moineau S., et al. . (2007). CRISPR provides acquired resistance against viruses in prokaryotes. Science 315, 1709–1712. 10.1126/science.1138140 [PubMed] [Cross Ref]
  • Besemer J., Borodovsky M. (1999). Heuristic approach to deriving models for gene finding. Nucleic Acids Res. 27, 3911–3920. 10.1093/nar/27.19.3911 [PMC free article] [PubMed] [Cross Ref]
  • Binetti A. G., Del Rio B., Martin M. C., Alvarez M. A. (2005). Detection and characterization of Streptococcus thermophilus bacteriophages by use of the antireceptor gene sequence. Appl. Environ. Microbiol. 71, 6096–6103. 10.1128/AEM.71.10.6096-6103.2005 [PMC free article] [PubMed] [Cross Ref]
  • Bolotin A., Quinquis B., Renault P., Sorokin A., Ehrlich S. D., Kulakauskas S., et al. . (2004). Complete sequence and comparative genome analysis of the dairy bacterium Streptococcus thermophilus. Nat. Biotechnol. 22, 1554–1558. 10.1038/nbt1034 [PubMed] [Cross Ref]
  • Bradley D. E. (1967). Ultrastructure of bacteriophage and bacteriocins. Bacteriol. Rev. 31, 230–314. [PMC free article] [PubMed]
  • Brondsted L., Pedersen M., Hammer K. (2001). An activator of transcription regulates phage TP901-1 late gene expression. Appl. Environ. Microbiol. 67, 5626–5633. 10.1128/AEM.67.12.5626-5633.2001 [PMC free article] [PubMed] [Cross Ref]
  • Brussow H., Desiere F. (2001). Comparative phage genomics and the evolution of Siphoviridae: insights from dairy phages. Mol. Microbiol. 39, 213–222. 10.1046/j.1365-2958.2001.02228.x [PubMed] [Cross Ref]
  • Bruttin A., Desiere F., d'Amico N., Guerin J. P., Sidoti J., Huni B., et al. . (1997). Molecular ecology of Streptococcus thermophilus bacteriophage infections in a cheese factory. Appl. Environ. Microbiol. 63, 3144–3150. [PMC free article] [PubMed]
  • Burrus V., Bontemps C., Decaris B., Guedon G. (2001). Characterization of a novel type II restriction-modification system, Sth368I, encoded by the integrative element ICESt1 of Streptococcus thermophilus CNRZ368. Appl. Environ. Microbiol. 67, 1522–1528. 10.1128/AEM.67.4.1522-1528.2001 [PMC free article] [PubMed] [Cross Ref]
  • Caldwell S. L., McMahon D. J., Oberg C. J., Broadbent J. R. (1996). Development and characterization of lactose-positive Pediococcus species for milk fermentation. Appl. Environ. Microbiol. 62, 936–941. [PMC free article] [PubMed]
  • Casey E., Mahony J., Neve H., Noben J.-P., Dal Bello F., van Sinderen D. (2015). Novel phage group infecting Lactobacillus delbrueckii subsp. lactis, as revealed by genomic and proteomic analysis of bacteriophage Ldl1. Appl. Environ. Microbiol. 81, 1319–1326. 10.1128/AEM.03413-14 [PMC free article] [PubMed] [Cross Ref]
  • Casey E., Mahony J., O'Connell-Motherway M., Bottacini F., Cornelissen A., Neve H., et al. . (2014). Molecular characterization of three Lactobacillus delbrueckii subsp. bulgaricus phages. Appl. Environ. Microbiol. 80, 5623–5635. 10.1128/AEM.01268-14 [PMC free article] [PubMed] [Cross Ref]
  • Ceyssens P.-J., Mesyanzhinov V., Sykilinda N., Briers Y., Roucourt B., Lavigne R., et al. . (2008). The genome and structural proteome of YuA, a new Pseudomonas aeruginosa phage resembling M6. J. Bacteriol. 190, 1429–1435. 10.1128/JB.01441-07 [PMC free article] [PubMed] [Cross Ref]
  • Chevreux B., Pfisterer T., Drescher B., Driesel A. J., Muller W. E., Wetter T., et al. . (2004). Using the miraEST assembler for reliable and automated mRNA transcript assembly and SNP detection in sequenced ESTs. Genome Res. 14, 1147–1159. 10.1101/gr.1917404 [PubMed] [Cross Ref]
  • Collins B., Bebeacua C., Mahony J., Blangy S., Douillard F. P., Veesler D., et al. . (2013). Structure and functional analysis of the host recognition device of lactococcal phage Tuc2009. J. Virol. 87, 8429–8440. 10.1128/JVI.00907-13 [PMC free article] [PubMed] [Cross Ref]
  • Cornelissen A., Ceyssens P. J., Krylov V. N., Noben J. P., Volckaert G., Lavigne R. (2012). Identification of EPS-degrading activity within the tail spikes of the novel Pseudomonas putida phage AF. Virology 434, 251–256. 10.1016/j.virol.2012.09.030 [PubMed] [Cross Ref]
  • Datta S., Costantino N., Zhou X., Court D. L. (2008). Identification and analysis of recombineering functions from Gram-negative and Gram-positive bacteria and their phages. Proc. Natl. Acad. Sci. U.S.A. 105, 1626–1631. 10.1073/pnas.0709089105 [PubMed] [Cross Ref]
  • Desiere F., Lucchini S., Brussow H. (1998). Evolution of Streptococcus thermophilus bacteriophage genomes by modular exchanges followed by point mutations and small deletions and insertions. Virology 241, 345–356. 10.1006/viro.1997.8959 [PubMed] [Cross Ref]
  • Desiere F., Lucchini S., Bruttin A., Zwahlen M. C., Brussow H. (1997). A highly conserved DNA replication module from Streptococcus thermophilus phages is similar in sequence and topology to a module from Lactococcus lactis phages. Virology 234, 372–382. 10.1006/viro.1997.8643 [PubMed] [Cross Ref]
  • Desiere F., McShan W. M., van Sinderen D., Ferretti J. J., Brussow H. (2001). Comparative genomics reveals close genetic relationships between phages from dairy bacteria and pathogenic streptococci: evolutionary implications for prophage-host interactions. Virology 288, 325–341. 10.1006/viro.2001.1085 [PubMed] [Cross Ref]
  • Deveau H., Barrangou R., Garneau J. E., Labont,é J., Fremaux C., Boyaval P., et al. . (2008). Phage response to CRISPR-encoded resistance in Streptococcus thermophilus. J. Bacteriol. 190, 1390–1400. 10.1128/JB.01412-07 [PMC free article] [PubMed] [Cross Ref]
  • Duplessis M., Levesque C. M., Moineau S. (2006). Characterization of Streptococcus thermophilus host range phage mutants. Appl. Environ. Microbiol. 72, 3036–3041. 10.1128/AEM.72.4.3036-3041.2006 [PMC free article] [PubMed] [Cross Ref]
  • Duplessis M., Moineau S. (2001). Identification of a genetic determinant responsible for host specificity in Streptococcus thermophilus bacteriophages. Mol. Microbiol. 41, 325–336. 10.1046/j.1365-2958.2001.02521.x [PubMed] [Cross Ref]
  • Dupuis M. E., Moineau S. (2010). Genome organization and characterization of the virulent lactococcal phage 1358 and its similarities to Listeria phages. Appl. Environ. Microbiol. 76, 1623–1632. 10.1128/AEM.02173-09 [PMC free article] [PubMed] [Cross Ref]
  • Finn R. D., Bateman A., Clements J., Coggill P., Eberhardt R. Y., Eddy S. R., et al. . (2014). Pfam: the protein families database. Nucleic Acids Res. 42, D222–D230. 10.1093/nar/gkt1223 [PMC free article] [PubMed] [Cross Ref]
  • Foley S., Bruttin A., Brussow H. (2000). Widespread distribution of a group I intron and its three deletion derivatives in the lysin gene of Streptococcus thermophilus bacteriophages. J. Virol. 74, 611–618. 10.1128/JVI.74.2.611-618.2000 [PMC free article] [PubMed] [Cross Ref]
  • Fortier L. C., Bransi A., Moineau S. (2006). Genome sequence and global gene expression of Q54, a new phage species linking the 936 and c2 phage species of Lactococcus lactis. J. Bacteriol. 188, 6101–6114. 10.1128/JB.00581-06 [PMC free article] [PubMed] [Cross Ref]
  • Garneau J. E., Moineau S. (2011). Bacteriophages of lactic acid bacteria and their impact on milk fermentations. Microb. Cell Fact. 10(Suppl. 1):S20. 10.1186/1475-2859-10-S1-S20 [PMC free article] [PubMed] [Cross Ref]
  • Garvey P., Hill C., Fitzgerald G. (1996). The lactococcal plasmid pNP40 encodes a third bacteriophage resistance mechanism, one which affects phage DNA penetration. Appl. Environ. Microbiol. 62, 676–679. [PMC free article] [PubMed]
  • Goh Y. J., Goin C., O'Flaherty S., Altermann E., Hutkins R. (2011). Specialized adaptation of a lactic acid bacterium to the milk environment: the comparative genomics of Streptococcus thermophilus LMD-9. Microb Cell Fact. 10(Suppl. 1):S22. 10.1186/1475-2859-10-S1-S22 [PMC free article] [PubMed] [Cross Ref]
  • Guglielmotti D. M., Binetti A., Reinheimer J., Quiberoni A. (2009a). Streptococcus thermophilus phage monitoring in a cheese factory: phage characteristics and starter sensitivity. Int. Dairy J. 19, 476–480. 10.1016/j.idairyj.2009.02.009 [Cross Ref]
  • Guglielmotti D. M., Deveau H., Binetti A. G., Reinheimer J. A., Moineau S., Quiberoni A. (2009b). Genome analysis of two virulent Streptococcus thermophilus phages isolated in Argentina. Int. J. Food Microbiol. 136, 101–109. 10.1016/j.ijfoodmicro.2009.09.005 [PubMed] [Cross Ref]
  • Hall S. D., Kolodner R. D. (1994). Homologous pairing and strand exchange promoted by the Escherichia coli RecT protein. Proc. Natl. Acad. Sci. U.S.A. 91, 3205–3209. 10.1073/pnas.91.8.3205 [PubMed] [Cross Ref]
  • Hols P., Hancy F., Fontaine L., Grossiord B., Prozzi D., Leblond-Bourget N., et al. . (2005). New insights in the molecular biology and physiology of Streptococcus thermophilus revealed by comparative genomics. FEMS Microbiol. Rev. 29, 435–463. 10.1016/j.femsre.2005.04.008 [PubMed] [Cross Ref]
  • Iyer L. M., Koonin E. V., Aravind L. (2002). Classification and evolutionary history of the single-strand annealing proteins, RecT, Redbeta, ERF and RAD52. BMC Genomics 3:8. 10.1186/1471-2164-3-8 [PMC free article] [PubMed] [Cross Ref]
  • Jones D. T. (1999). Protein secondary structure prediction based on position-specific scoring matrices. J. Mol. Biol. 292, 195–202. 10.1006/jmbi.1999.3091 [PubMed] [Cross Ref]
  • Kala S., Cumby N., Sadowski P. D., Hyder B. Z., Kanelis V., Davidson A. R., et al. . (2014). HNH proteins are a widespread component of phage DNA packaging machines. Proc. Natl. Acad. Sci. U.S.A. 111, 6022–6027. 10.1073/pnas.1320952111 [PubMed] [Cross Ref]
  • Kelleher P., Murphy J., Mahony J., van Sinderen D. (2015). Next-generation sequencing as an approach to dairy starter selection. Dairy Sci. Technol. 95, 545–568. 10.1007/s13594-015-0227-4 [PMC free article] [PubMed] [Cross Ref]
  • Kenny J. G., McGrath S., Fitzgerald G. F., van Sinderen D. (2004). Bacteriophage Tuc2009 encodes a tail-associated cell wall-degrading activity. J. Bacteriol. 186, 3480–3491. 10.1128/JB.186.11.3480-3491.2004 [PMC free article] [PubMed] [Cross Ref]
  • Kuipers O. P., de Ruyter P. G., Kleerebezem M., de Vos W. M. (1998). Quorum sensing-controlled gene expression in lactic acid bacteria. J. Biotechnol. 64, 15–21. 10.1016/S0168-1656(98)00100-X [Cross Ref]
  • Kumar S., Cheng X., Klimasauskas S., Mi S., Posfai J., Roberts R. J., et al. . (1994). The DNA (cytosine-5) methyltransferases. Nucleic Acids Res. 22, 1. 10.1093/nar/22.1.1 [PMC free article] [PubMed] [Cross Ref]
  • Labrie S. J., Josephsen J., Neve H., Vogensen F. K., Moineau S. (2008). Morphology, genome sequence, and structural proteome of type phage P335 from Lactococcus lactis. Appl. Environ. Microbiol. 74, 4636–4644. 10.1128/AEM.00118-08 [PMC free article] [PubMed] [Cross Ref]
  • Labrie S. J., Tremblay D. M., Plante P. L., Wasserscheid J., Dewar K., Corbeil J., et al. . (2015). Complete genome sequence of Streptococcus thermophilus SMQ-301, a model strain for phage-host interactions. Genome Announc. 3:e00480–15. 10.1128/genomeA.00480-15 [PMC free article] [PubMed] [Cross Ref]
  • Lacks S., Greenberg B. (1977). Complementary specificity of restriction endonucleases of Diplococcus pneumoniae with respect to DNA methylation. J. Mol. Biol. 114, 153–168. 10.1016/0022-2836(77)90289-3 [PubMed] [Cross Ref]
  • Lansky S., Alalouf O., Solomon H. V., Alhassid A., Govada L., Chayen N. E., et al. . (2014). A unique octameric structure of Axe2, an intracellular acetyl-xylooligosaccharide esterase from Geobacillus stearothermophilus. Acta Crystallogr. D Biol. Crystallogr. 70(Pt 2), 261–278. 10.1107/S139900471302840X [PubMed] [Cross Ref]
  • Le Marrec C., van Sinderen D., Walsh L., Stanley E., Vlegels E., Moineau S., et al. . (1997). Two groups of bacteriophages infecting Streptococcus thermophilus can be distinguished on the basis of mode of packaging and genetic determinants for major structural proteins. Appl. Environ. Microbiol. 63, 3246–3253. [PMC free article] [PubMed]
  • Levesque C., Duplessis M., Labonte J., Labrie S., Fremaux C., Tremblay D., et al. . (2005). Genomic organization and molecular analysis of virulent bacteriophage 2972 infecting an exopolysaccharide-producing Streptococcus thermophilus strain. Appl. Environ. Microbiol. 71, 4057–4068. 10.1128/AEM.71.7.4057-4068.2005 [PMC free article] [PubMed] [Cross Ref]
  • Lillehaug D. (1997). An improved plaque assay for poor plaque-producing temperate lactococcal bacteriophages. J. Appl. Microbiol. 83, 85–90. 10.1046/j.1365-2672.1997.00193.x [PubMed] [Cross Ref]
  • Lucchini S., Desiere F., Brussow H. (1999). Comparative genomics of Streptococcus thermophilus phage species supports a modular evolution theory. J. Virol. 73, 8647–8656. [PMC free article] [PubMed]
  • Lugli G. A., Milani C., Mancabelli L., van Sinderen D., Ventura M. (2016). MEGAnnotator: a user-friendly pipeline for microbial genomes assembly and annotation. FEMS Microbiol. Lett. 363:fnw049. 10.1093/femsle/fnw049 [PubMed] [Cross Ref]
  • Macmaster R., Sedelnikova S., Baker P. J., Bolt E. L., Lloyd R. G., Rafferty J. B. (2006). RusA Holliday junction resolvase: DNA complex structure–insights into selectivity and specificity. Nucleic Acids Res. 34, 5577–5584. 10.1093/nar/gkl447 [PMC free article] [PubMed] [Cross Ref]
  • Mahony J., Martel B., Tremblay D. M., Neve H., Heller K. J., Moineau S., et al. . (2013). Identification of a new P335 subgroup through molecular analysis of lactococcal phages Q33 and BM13. Appl. Environ. Microbiol. 79, 4401–4409. 10.1128/AEM.00832-13 [PMC free article] [PubMed] [Cross Ref]
  • Mahony J., Randazzo W., Neve H., Settanni L., van Sinderen D. (2015). Lactococcal 949 group phages recognize a carbohydrate receptor on the host cell surface. Appl. Environ. Microbiol. 81, 3299–3305. 10.1128/AEM.00143-15 [PMC free article] [PubMed] [Cross Ref]
  • Marchler-Bauer A., Derbyshire M. K., Gonzales N. R., Lu S., Chitsaz F., Geer L. Y., et al. . (2015). CDD: NCBI's conserved domain database. Nucleic Acids Res. 43, D222–D226. 10.1093/nar/gku1221 [PMC free article] [PubMed] [Cross Ref]
  • McDonnell B., Mahony J., Neve H., Hanemaaijer L., Noben J. P., Kouwen T., et al. . (2016). Identification and analysis of a novel group of bacteriophages infecting the lactic acid bacterium Streptococcus thermophilus. Appl. Environ. Microbiol. 82, 5153–5165. 10.1128/AEM.00835-16 [PMC free article] [PubMed] [Cross Ref]
  • Mc Grath S., Neve H., Seegers J. F., Eijlander R., Vegge C. S., Brondsted L., et al. . (2006). Anatomy of a lactococcal phage tail. J. Bacteriol. 188, 3972–3982. 10.1128/JB.00024-06 [PMC free article] [PubMed] [Cross Ref]
  • Mills S., Griffin C., O'Sullivan O., Coffey A., McAuliffe O., Meijer W., et al. (2011). A new phage on the ‘Mozzarella’ block: bacteriophage 5093 shares a low level of homology with other Streptococcus thermophilus phages. Int. Dairy J. 21, 963–969. 10.1016/j.idairyj.2011.06.003 [Cross Ref]
  • Moineau S., Pandian S., Klaenhammer T. R. (1994). Evolution of a lytic bacteriophage via DNA acquisition from the Lactococcus lactis chromosome. Appl. Environ. Microbiol. 60, 1832–1841. [PMC free article] [PubMed]
  • Murphy J., Klumpp J., Mahony J., O'Connell-Motherway M., Nauta A., van Sinderen D. (2014). Methyltransferases acquired by lactococcal 936-type phage provide protection against restriction endonuclease activity. BMC Genomics 15:831. 10.1186/1471-2164-15-831 [PMC free article] [PubMed] [Cross Ref]
  • Murphy J., Mahony J., Ainsworth S., Nauta A., van Sinderen D. (2013). Bacteriophage orphan DNA methyltransferases: insights from their bacterial origin, function, and occurrence. Appl. Environ. Microbiol. 79, 7547–7555. 10.1128/AEM.02229-13 [PMC free article] [PubMed] [Cross Ref]
  • Noirot P., Kolodner R. D. (1998). DNA strand invasion promoted by Escherichia coli RecT protein. J. Biol. Chem. 273, 12274–12280. 10.1074/jbc.273.20.12274 [PubMed] [Cross Ref]
  • Obregon V., Garcia J. L., Garcia E., Lopez R., Garcia P. (2003). Genome organization and molecular analysis of the temperate bacteriophage MM1 of Streptococcus pneumoniae. J. Bacteriol. 185, 2362–2368. 10.1128/JB.185.7.2362-2368.2003 [PMC free article] [PubMed] [Cross Ref]
  • Passy S. I., Yu X., Li Z., Radding C. M., Egelman E. H. (1999). Rings and filaments of beta protein from bacteriophage lambda suggest a superfamily of recombination proteins. Proc. Natl. Acad. Sci. U.S.A. 96, 4279–4284. 10.1073/pnas.96.8.4279 [PubMed] [Cross Ref]
  • Pedersen M., Ostergaard S., Bresciani J., Vogensen F. K. (2000). Mutational analysis of two structural genes of the temperate lactococcal bacteriophage TP901-1 involved in tail length determination and baseplate assembly. Virology 276, 315–328. 10.1006/viro.2000.0497 [PubMed] [Cross Ref]
  • Proux C., van Sinderen D., Suarez J., Garcia P., Ladero V., Fitzgerald G. F., et al. . (2002). The dilemma of phage taxonomy illustrated by comparative genomics of Sfi21-like Siphoviridae in lactic acid bacteria. J. Bacteriol. 184, 6026–6036. 10.1128/JB.184.21.6026-6036.2002 [PMC free article] [PubMed] [Cross Ref]
  • Quiberoni A., Tremblay D., Ackermann H. W., Moineau S., Reinheimer J. A. (2006). Diversity of Streptococcus thermophilus phages in a large-production cheese factory in Argentina. J. Dairy Sci. 89, 3791–3799. 10.3168/jds.S0022-0302(06)72420-1 [PubMed] [Cross Ref]
  • Quiles-Puchalt N., Tormo-Mas M. A., Campoy S., Toledo-Arana A., Monedero V., Lasa I., et al. . (2013). A super-family of transcriptional activators regulates bacteriophage packaging and lysis in Gram-positive bacteria. Nucleic Acids Res. 41, 7260–7275. 10.1093/nar/gkt508 [PMC free article] [PubMed] [Cross Ref]
  • Roberts R. J., Vincze T., Posfai J., Macelis D. (2015). REBASE–a database for DNA restriction and modification: enzymes, genes and genomes. Nucleic Acids Res. 43, D298–D299. 10.1093/nar/gku1046 [PMC free article] [PubMed] [Cross Ref]
  • Sambrook J., Fritsch E. F., Maniatis T. (1989). Molecular Cloning. New York, NY: Cold Spring Harbour Laboratory Press.
  • Sampath J., Vijayakumar M. N. (1998). Identification of a DNA cytosine methyltransferase gene in conjugative transposon Tn5252. Plasmid 39, 63–76. 10.1006/plas.1997.1316 [PubMed] [Cross Ref]
  • Sharples G. J., Chan S. N., Mahdi A. A., Whitby M. C., Lloyd R. G. (1994). Processing of intermediates in recombination and DNA repair: identification of a new endonuclease that specifically cleaves Holliday junctions. EMBO J. 13, 6133–6142. [PubMed]
  • Soding J., Biegert A., Lupas A. N. (2005). The HHpred interactive server for protein homology detection and structure prediction. Nucleic Acids Res. 33, W244–W248. 10.1093/nar/gki408 [PMC free article] [PubMed] [Cross Ref]
  • Stanley E., Fitzgerald G. F., Le Marrec C., Fayard B., van Sinderen D. (1997). Sequence analysis and characterization of phi O1205, a temperate bacteriophage infecting Streptococcus thermophilus CNRZ1205. Microbiology 143(Pt 11), 3417–3429. 10.1099/00221287-143-11-3417 [PubMed] [Cross Ref]
  • Stanley E., Walsh L., van der Zwet A., Fitzgerald G. F., van Sinderen D. (2000). Identification of four loci isolated from two Streptococcus thermophilus phage genomes responsible for mediating bacteriophage resistance. FEMS Microbiol. Lett. 182, 271–277. 10.1111/j.1574-6968.2000.tb08907.x [PubMed] [Cross Ref]
  • Stockdale S. R., Mahony J., Courtin P., Chapot-Chartier M. P., van Pijkeren J. P., Britton R. A., et al. . (2013). The lactococcal phages Tuc2009 and TP901-1 incorporate two alternate forms of their tail fiber into their virions for infection specialization. J. Biol. Chem. 288, 5581–5590. 10.1074/jbc.M112.444901 [PMC free article] [PubMed] [Cross Ref]
  • Sun Z., Chen X., Wang J., Zhao W., Shao Y., Wu L., et al. . (2011). Complete genome sequence of Streptococcus thermophilus strain ND03. J. Bacteriol. 193, 793–794. 10.1128/JB.01374-10 [PMC free article] [PubMed] [Cross Ref]
  • Szymczak P., Janzen T., Neves A. R., Kot W., Hansen L. H., Lametsch R., et al. (2016). Novel variants of Streptococcus thermophilus bacteriophages indicate genetic recombination across phages from different bacterial species. Appl, Environ, Microbiol. 83:e02748–16. 10.1128/AEM.02748-16 [PMC free article] [PubMed] [Cross Ref]
  • Tamura K., Stecher G., Peterson D., Filipski A., Kumar S. (2013). MEGA6: molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 30, 2725–2729. 10.1093/molbev/mst197 [PMC free article] [PubMed] [Cross Ref]
  • Tettelin H., Masignani V., Cieslewicz M. J., Donati C., Medini D., Ward N. L., et al. . (2005). Genome analysis of multiple pathogenic isolates of Streptococcus agalactiae: implications for the microbial “pan-genome.” Proc. Natl. Acad. Sci. U.S.A. 102, 13950–13955. 10.1073/pnas.0506758102 [PubMed] [Cross Ref]
  • Tremblay D. M., Moineau S. (1999). Complete genomic sequence of the lytic bacteriophage DT1 of Streptococcus thermophilus. Virology 255, 63–76. 10.1006/viro.1998.9525 [PubMed] [Cross Ref]
  • Treu L., Vendramin V., Bovo B., Campanaro S., Corich V., Giacomini A. (2014). Genome sequences of Streptococcus thermophilus strains MTH17CL396 and M17PTZA496 from fontina, an Italian PDO cheese. Genome Announc. 2:e00067–1. 10.1128/genomeA.00067-14 [PMC free article] [PubMed] [Cross Ref]
  • Vanheel A., Daniels R., Plaisance S., Baeten K., Hendriks J., Leprince P., et al. . (2012). Identification of protein networks involved in the disease course of experimental autoimmune encephalomyelitis, an animal model of multiple sclerosis. PLoS ONE 7:e35544. 10.1371/journal.pone.0035544 [PMC free article] [PubMed] [Cross Ref]
  • van Sinderen D., Karsens H., Kok J., Terpstra P., Ruiters M. H., Venema G., et al. . (1996). Sequence analysis and molecular characterization of the temperate lactococcal bacteriophage r1t. Mol. Microbiol. 19, 1343–1355. 10.1111/j.1365-2958.1996.tb02478.x [PubMed] [Cross Ref]
  • Vegge C. S., Brøndsted L., Neve H., Mc Grath S., van Sinderen D., Vogensen F. K. (2005). Structural characterization and assembly of the distal tail structure of the temperate lactococcal bacteriophage TP901-1. J. Bacteriol. 187, 4187–4197. 10.1128/JB.187.12.4187-4197.2005 [PMC free article] [PubMed] [Cross Ref]
  • Ventura M., Foley S., Bruttin A., Chennoufi S. C., Canchaya C., Brussow H. (2002). Transcription mapping as a tool in phage genomics: the case of the temperate Streptococcus thermophilus phage Sfi21. Virology 296, 62–76. 10.1006/viro.2001.1331 [PubMed] [Cross Ref]
  • Wu Q., Tun H. M., Leung F. C., Shah N. P. (2014). Genomic insights into high exopolysaccharide-producing dairy starter bacterium Streptococcus thermophilus ASCC 1275. Sci. Rep. 4:4974. 10.1038/srep04974 [PMC free article] [PubMed] [Cross Ref]
  • Zago M., Carminati D., Giraffa G. (2003). Characterization of Streptococcus thermophilus phages from cheese. Ann. Microbiol. 53, 171–178.
  • Zhao Y., Wu J., Yang J., Sun S., Xiao J., Yu J. (2012). PGAP: pan-genomes analysis pipeline. Bioinformatics 28, 416–418. 10.1093/bioinformatics/btr655 [PMC free article] [PubMed] [Cross Ref]
  • Zinno P., Janzen T., Bennedsen M., Ercolini D., Mauriello G. (2010). Characterization of Streptococcus thermophilus lytic bacteriophages from mozzarella cheese plants. Int. J. Food Microbiol. 138, 137–144. 10.1016/j.ijfoodmicro.2009.12.008 [PubMed] [Cross Ref]

Articles from Frontiers in Microbiology are provided here courtesy of Frontiers Media SA