Longitudinal Study of Lactococcus Phages in a Canadian Cheese Factory.

The presence of virulent phages is closely monitored during cheese manufacturing, as these bacterial viruses can significantly slow down the milk fermentation process and lead to low-quality cheeses. From 2001 to 2020, whey samples from cheddar cheese production in a Canadian factory were monitored for the presence of virulent phages capable of infecting proprietary strains of Lactococcus cremoris and Lactococcus lactis used in starter cultures. Phages were successfully isolated from 932 whey samples using standard plaque assays and several industrial Lactococcus strains as hosts. A multiplex PCR assay assigned 97% of these phage isolates to the Skunavirus genus, 2% to the P335 group, and 1% to the Ceduovirus genus. DNA restriction profiles and a multilocus sequence typing (MLST) scheme distinguished at least 241 unique lactococcal phages from these isolates. While most phages were isolated only once, 93 of them (out of 241, 39%) were isolated multiple times. Phage GL7 was isolated 132 times from 2006 to 2020, demonstrating that phages can persist in a cheese factory for long periods of time. Phylogenetic analysis of MLST sequences showed that phages could be clustered based on their bacterial hosts rather than their year of isolation. Host range analysis showed that Skunavirus phages exhibited a very narrow host range, whereas some Ceduovirus and P335 phages had a broader host range. Overall, the host range information was useful in improving the starter culture rotation by identifying phage-unrelated strains and helped mitigating the risk of fermentation failure due to virulent phages. IMPORTANCE Although lactococcal phages have been observed in cheese production settings for almost a century, few longitudinal studies have been performed. This 20-year study describes the close monitoring of dairy lactococcal phages in a cheddar cheese factory. Routine monitoring was conducted by factory staff, and when whey samples were found to inhibit industrial starter cultures under laboratory conditions, they were sent to an academic research laboratory for phage isolation and characterization. This led to a collection of at least 241 unique lactococcal phages, which were characterized through PCR typing and MLST profiling. Phages of the Skunavirus genus were by far the most dominant. Most phages lysed a small subset of the Lactococcus strains. These findings guided the industrial partner in adapting the starter culture schedule by using phage-unrelated strains in starter cultures and removing some strains from the starter rotation. This phage control strategy could be adapted for other large-scale bacterial fermentation processes.

DNA tandem repeats contribute to the genetic diversity of Brevibacterium aurantiacum phages.

This report presents the characterization of the first virulent phages infecting Brevibacterium aurantiacum, a bacterial species used during the manufacture of surface-ripened cheeses. These phages were also responsible for flavour and colour defects in surface-ripened cheeses. Sixteen phages (out of 62 isolates) were selected for genome sequencing and comparative analyses. These cos-type phages with a long non-contractile tail currently belong to the Siphoviridae family (Caudovirales order). Their genome sizes vary from 35,637 to 36,825 bp and, similar to their host, have a high GC content (~61%). Genes encoding for an immunity repressor, an excisionase and a truncated integrase were found, suggesting that these virulent phages may be derived from a prophage. Their genomic organization is highly conserved, with most of the diversity coming from the presence of long (198 bp) DNA tandem repeats (TRs) within an open reading frame coding for a protein of unknown function. We categorized these phages into seven genomic groups according to their number of TR, which ranged from two to eight. Moreover, we showed that TRs are widespread in phage genomes, found in more than 85% of the genomes available in public databases.

The Tape Measure Protein Is Involved in the Heat Stability of Lactococcus lactis Phages.

Virulent lactococcal phages are still a major risk for milk fermentation processes as they may lead to slowdowns and low-quality fermented dairy products, particularly cheeses. Some of the phage control strategies used by the industry rely on heat treatments. Recently, a few Lactococcus lactis phages were found to be highly thermo-resistant. To identify the genetic determinant(s) responsible for the thermal resistance of lactococcal phages, we used the virulent phage CB14 (of the Lactococcus lactis 936 [now Sk1virus] phage group) to select for phage mutants with increased heat stability. By treating phage CB14 to successive low and high temperatures, we were able to select two CB14 derivatives with increased heat stability. Sequencing of their genome revealed the same nucleotide sequences as the wild-type phage CB14, except for a same-sized deletion (120 bp) in the gene coding for the tape measure protein (TMP) of each phage mutant, but at a different position. The TMP protein sequences of these mutant phages were compared with their homologues in other wild-type L. lactis phages with a wide diversity in heat stability. Comparative analysis showed that the same nucleotide deletion appears to have also occurred in the gene coding for the TMP of highly thermo-resistant lactococcal phages P1532 and P680. We propose that the TMP is, in part, responsible for the heat stability of the highly predominant lactococcal phages of the Sk1virus group.IMPORTANCE Virulent lactococcal phages still represent a major risk for milk fermentation as they may lead to slowdowns and low-quality fermented dairy products. Heat treatment is one of the most commonly used methods to control these virulent phages in cheese by-products. Recently, a few Lactococcus lactis phages, members of the Sk1virus group, have emerged with high thermal stability. To our knowledge, the genetic determinant(s) responsible for this thermal resistance in lactococcal phages is unknown. A better understanding of the thermal stability of these emerging virulent lactococcal phages is needed to improve industrial control strategies. In this work, we report the identification of a phage structural protein that is involved in the heat stability of a virulent Sk1virus phage. Identifying such a genetic determinant for heat stability is a first step in understanding the emergence of this group of thermostable phages.

Genomic Diversity of Phages Infecting Probiotic Strains of Lactobacillus paracasei.

Strains of the Lactobacillus casei group have been extensively studied because some are used as probiotics in foods. Conversely, their phages have received much less attention. We analyzed the complete genome sequences of five L. paracasei temperate phages: CL1, CL2, iLp84, iLp1308, and iA2. Only phage iA2 could not replicate in an indicator strain. The genome lengths ranged from 34,155 bp (iA2) to 39,474 bp (CL1). Phages iA2 and iLp1308 (34,176 bp) possess the smallest genomes reported, thus far, for phages of the L. casei group. The GC contents of the five phage genomes ranged from 44.8 to 45.6%. As observed with many other phages, their genomes were organized as follows: genes coding for DNA packaging, morphogenesis, lysis, lysogeny, and replication. Phages CL1, CL2, and iLp1308 are highly related to each other. Phage iLp84 was also related to these three phages, but the similarities were limited to gene products involved in DNA packaging and structural proteins. Genomic fragments of phages CL1, CL2, iLp1308, and iLp84 were found in several genomes of L. casei strains. Prophage iA2 is unrelated to these four phages, but almost all of its genome was found in at least four L. casei strains. Overall, these phages are distinct from previously characterized Lactobacillus phages. Our results highlight the diversity of L. casei phages and indicate frequent DNA exchanges between phages and their hosts.

A virulent phage infecting Lactococcus garvieae, with homology to Lactococcus lactis phages.

A new virulent phage belonging to the Siphoviridae family and able to infect Lactococcus garvieae strains was isolated from compost soil. Phage GE1 has a prolate capsid (56 by 38 nm) and a long noncontractile tail (123 nm). It had a burst size of 139 and a latent period of 31 min. Its host range was limited to only two L. garvieae strains out of 73 tested. Phage GE1 has a double-stranded DNA genome of 24,847 bp containing 48 predicted open reading frames (ORFs). Putative functions could be assigned to only 14 ORFs, and significant matches in public databases were found for only 17 ORFs, indicating that GE1 is a novel phage and its genome contains several new viral genes and encodes several new viral proteins. Of these 17 ORFs, 16 were homologous to deduced proteins of virulent phages infecting the dairy bacterium Lactococcus lactis, including previously characterized prolate-headed phages. Comparative genome analysis confirmed the relatedness of L. garvieae phage GE1 to L. lactis phages c2 (22,172 bp) and Q54 (26,537 bp), although its genome organization was closer to that of phage c2. Phage GE1 did not infect any of the 58 L. lactis strains tested. This study suggests that phages infecting different lactococcal species may have a common ancestor.

A new Microviridae phage isolated from a failed biotechnological process driven by Escherichia coli.

Bacteriophages are present in every environment that supports bacterial growth, including man made ecological niches. Virulent phages may even slow or, in more severe cases, interrupt bioprocesses driven by bacteria. Escherichia coli is one of the most widely used bacteria for large-scale bioprocesses; however, literature describing phage-host interactions in this industrial context is sparse. Here, we describe phage MED1 isolated from a failed industrial process. Phage MED1 (Microviridae family, with a single-stranded DNA [ssDNA] genome) is highly similar to the archetypal phage phiX174, sharing >95% identity between their genomic sequences. Whole-genome phylogenetic analysis of 52 microvirus genomes from public databases revealed three genotypes (alpha3, G4, and phiX174). Phage MED1 belongs to the phiX174 group. We analyzed the distribution of single nucleotide variants in MED1 and 18 other phiX174-like genomes and found that there are more missense mutations in genes G, B, and E than in the other genes of these genomes. Gene G encodes the spike protein, involved in host attachment. The evolution of this protein likely results from the selective pressure on phages to rapidly adapt to the molecular diversity found at the surface of their hosts.

Diversity of Salmonella enterica phages isolated from chicken farms in Kenya.

IMPORTANCE: Non-typhoidal Salmonella enterica infections are one of the leading causes of diarrhoeal diseases that spread to humans from animal sources such as poultry. Hence, keeping poultry farms free of Salmonella is essential for consumer safety and for a better yield of animal products. However, the emergence of antibiotic resistance due to over usage has sped up the search for alternative biocontrol methods such as the use of bacteriophages. Isolation and characterization of novel bacteriophages are key to adapt phage-based biocontrol applications. Here, we isolated and characterized Salmonella phages from samples collected at chicken farms and slaughterhouses in Kenya. The genomic characterization of these phage isolates revealed that they belong to four ICTV (International Committee on Taxonomy of Viruses) phage genera. All these phages are lytic and possibly suitable for biocontrol applications because no lysogenic genes or virulence factors were found in their genomes. Hence, we recommend further studies on these phages for their applications in Salmonella biocontrol.

Complete genome of Escherichia phages Carena and JoYop.

Escherichia phages Carena and JoYop were isolated from water samples in Abidjan (Cote d'Ivoire). Their genomes comprise 39,283 and 169,193 bp, encoding 44 and 246 predicted genes, respectively. Carena shares 93.4% nucleotide identity with Escherichia podophage CarlSpitteler (Berlinvirus), and JoYop shows 95.6% identity with Escherichia myophage ADUt (Tequatrovirus).

Involvement of the major capsid protein and two early-expressed phage genes in the activity of the lactococcal abortive infection mechanism AbiT.

The dairy industry uses the mesophilic, Gram-positive, lactic acid bacterium (LAB) Lactococcus lactis to produce an array of fermented milk products. Milk fermentation processes are susceptible to contamination by virulent phages, but a plethora of phage control strategies are available. One of the most efficient is to use LAB strains carrying phage resistance systems such as abortive infection (Abi) mechanisms. Yet, the mode of action of most Abi systems remains poorly documented. Here, we shed further light on the antiviral activity of the lactococcal AbiT system. Twenty-eight AbiT-resistant phage mutants derived from the wild-type AbiT-sensitive lactococcal phages p2, bIL170, and P008 were isolated and characterized. Comparative genomic analyses identified three different genes that were mutated in these virulent AbiT-insensitive phage derivatives: e14 (bIL170 [e14(bIL170)]), orf41 (P008 [orf41(P008)]), and orf6 (p2 [orf6(p2)] and P008 [orf6(P008)]). The genes e14(bIL170) and orf41(P008) are part of the early-expressed genomic region, but bioinformatic analyses did not identify their putative function. orf6 is found in the phage morphogenesis module. Antibodies were raised against purified recombinant ORF6, and immunoelectron microscopy revealed that it is the major capsid protein (MCP). Coexpression in L. lactis of ORF6(p2) and ORF5(p2), a protease, led to the formation of procapsids. To our knowledge, AbiT is the first Abi system involving distinct phage genes.

Complete Genome Sequences of 10 Lactococcal Skunavirus Phages Isolated from Cheddar Cheese Whey Samples in Canada.

We report the complete genome sequences of 10 virulent phages of the Skunavirus genus (Siphoviridae) that infect Lactococcus lactis strains used for cheddar cheese production in Canada. Their linear genomes range from 28,969 bp to 31,042 bp with GC contents of 34.1 to 35.1% and 55 to 60 predicted open reading frames (ORFs).

Lactococcal abortive infection protein AbiV interacts directly with the phage protein SaV and prevents translation of phage proteins.

AbiV is an abortive infection protein that inhibits the lytic cycle of several virulent phages infecting Lactococcus lactis, while a mutation in the phage gene sav confers insensitivity to AbiV. In this study, we have further characterized the effects of the bacterial AbiV and its interaction with the phage p2 protein SaV. First, we showed that during phage infection of lactococcal AbiV(+) cells, AbiV rapidly inhibited protein synthesis. Among early phage transcripts, sav gene transcription was slightly inhibited while the SaV protein could not be detected. Analyses of other phage p2 mRNAs and proteins suggested that AbiV blocks the activation of late gene transcription, probably by a general inhibition of translation. Using size exclusion chromatography coupled with on-line static light scattering and refractometry, as well as fluorescence quenching experiments, we also demonstrated that both AbiV and SaV formed homodimers and that they strongly and specifically interact with each other to form a stable protein complex.

Complete Genome Sequences for Two Myoviridae Strains Infecting Cyanobacteria in a Subarctic Lake.

We isolated two closely related strains that belong to the Myoviridae family and infect cyanobacteria in a shallow subarctic rock basin lake. Their host was identified as a member of the Synechococcus-Cyanobium complex. Sequenced genomes of the two phages were 244,930 bp and 243,633 bp. We describe their annotation and highlight some noteworthy features.

Author Correction: A mutation in the methionine aminopeptidase gene provides phage resistance in Streptococcus thermophilus.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Mobilome of Brevibacterium aurantiacum Sheds Light on Its Genetic Diversity and Its Adaptation to Smear-Ripened Cheeses.

Brevibacterium aurantiacum is an actinobacterium that confers key organoleptic properties to washed-rind cheeses during the ripening process. Although this industrially relevant species has been gaining an increasing attention in the past years, its genome plasticity is still understudied due to the unavailability of complete genomic sequences. To add insights on the mobilome of this group, we sequenced the complete genomes of five dairy Brevibacterium strains and one non-dairy strain using PacBio RSII. We performed phylogenetic and pan-genome analyses, including comparisons with other publicly available Brevibacterium genomic sequences. Our phylogenetic analysis revealed that these five dairy strains, previously identified as Brevibacterium linens, belong instead to the B. aurantiacum species. A high number of transposases and integrases were observed in the Brevibacterium spp. strains. In addition, we identified 14 and 12 new insertion sequences (IS) in B. aurantiacum and B. linens genomes, respectively. Several stretches of homologous DNA sequences were also found between B. aurantiacum and other cheese rind actinobacteria, suggesting horizontal gene transfer (HGT). A HGT region from an iRon Uptake/Siderophore Transport Island (RUSTI) and an iron uptake composite transposon were found in five B. aurantiacum genomes. These findings suggest that low iron availability in milk is a driving force in the adaptation of this bacterial species to this niche. Moreover, the exchange of iron uptake systems suggests cooperative evolution between cheese rind actinobacteria. We also demonstrated that the integrative and conjugative element BreLI (Brevibacterium Lanthipeptide Island) can excise from B. aurantiacum SMQ-1417 chromosome. Our comparative genomic analysis suggests that mobile genetic elements played an important role into the adaptation of B. aurantiacum to cheese ecosystems.

A mutation in the methionine aminopeptidase gene provides phage resistance in Streptococcus thermophilus.

Streptococcus thermophilus is a lactic acid bacterium widely used by the dairy industry for the manufacture of yogurt and specialty cheeses. It is also a Gram-positive bacterial model to study phage-host interactions. CRISPR-Cas systems are one of the most prevalent phage resistance mechanisms in S. thermophilus. Little information is available about other host factors involved in phage replication in this food-grade streptococcal species. We used the model strain S. thermophilus SMQ-301 and its virulent phage DT1, harboring the anti-CRISPR protein AcrIIA6, to show that a host gene coding for a methionine aminopeptidase (metAP) is necessary for phage DT1 to complete its lytic cycle. A single mutation in metAP provides S. thermophilus SMQ-301 with strong resistance against phage DT1. The mutation impedes a late step of the lytic cycle since phage adsorption, DNA replication, and protein expression were not affected. When the mutated strain was complemented with the wild-type version of the gene, the phage sensitivity phenotype was restored. When this mutation was introduced into other S. thermophilus strains it provided resistance against cos-type (Sfi21dt1virus genus) phages but replication of pac-type (Sfi11virus genus) phages was not affected. The mutation in the gene coding for the MetAP induces amino acid change in a catalytic domain conserved across many bacterial species. Introducing the same mutation in Streptococcus mutans also provided a phage resistance phenotype, suggesting the wide-ranging importance of the host methionine aminopeptidase in phage replication.

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

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

CRISPRStudio: A User-Friendly Software for Rapid CRISPR Array Visualization.

The CRISPR-Cas system biologically serves as an adaptive defense mechanism against phages. However, there is growing interest in exploiting the hypervariable nature of the CRISPR locus, often of viral origin, for microbial typing and tracking. Moreover, the spacer content of any given strain provides a phage resistance profile. Large-scale CRISPR typing studies require an efficient method for showcasing CRISPR array similarities across multiple isolates. Historically, CRISPR arrays found in microbes have been represented by colored shapes based on nucleotide sequence identity and, while this approach is now routinely used, only scarce computational resources are available to automate the process, making it very time-consuming for large datasets. To alleviate this tedious task, we introduce CRISPRStudio, a command-line tool developed to accelerate CRISPR analysis and standardize the preparation of CRISPR array figures. It first compares nucleotide spacer sequences present in a dataset and then clusters them based on sequence similarity to assign a meaningful representative color. CRISPRStudio offers versatility to suit different biological contexts by including options such as automatic sorting of CRISPR loci and highlighting of shared spacers, while remaining fast and user-friendly.

Characterization of the Escherichia coli Virulent Myophage ST32.

The virulent phage ST32 that infects the Escherichiacoli strain ST130 was isolated from a wastewater sample in China and analyzed. Morphological observations showed that phage ST32 belongs to the Myoviridae family, as it has an icosahedral capsid and long contractile tail. Host range analysis showed that it exhibits a broad range of hosts including non-pathogenic and pathogenic E. coli strains. Interestingly, phage ST32 had a much larger burst size when amplified at 20 °C as compared to 30 °C or 37 °C. Its double-stranded DNA genome was sequenced and found to contain 53,092 bp with a GC content of 44.14%. Seventy-nine open reading frames (ORFs) were identified and annotated as well as a tRNA-Arg. Only nineteen ORFs were assigned putative functions. A phylogenetic tree using the large terminase subunit revealed a close relatedness with four unclassified Myoviridae phages. A comparative genomic analysis of these phages showed that the Enterobacteria phage phiEcoM-GJ1 is the closest relative to ST32 and shares the same new branch in the phylogenetic tree. Still, these two phages share only 47 of 79 ORFs with more than 90% identity. Phage ST32 has unique characteristics that make it a potential biological control agent under specific conditions.

Characterization of prophages of Lactococcus garvieae.

This report describes the morphological characterization and genome analysis of an induced prophage (PLg-TB25) from a dairy strain of Lactococcus garvieae. The phage belongs to the Siphoviridae family and its morphology is typical of other lactococcal phages. A general analysis of its genome did not reveal similarities with other lactococcal phage genomes, confirming its novelty. However, similarities were found between genes of its morphogenesis cluster and genes of Gram-positive bacteria, suggesting that this phage genome resulted from recombination events that took place in a heterogeneous microbial environment. An in silico search for other prophages in 16 L. garvieae genomes available in public databases, uncovered eight seemingly complete prophages in strains isolated from dairy and fish niches. Genome analyses of these prophages revealed three novel L. garvieae phages. The remaining prophages had homology to phages of Lactococcus lactis (P335 group) suggesting a close relationship between these lactococcal species. The similarity in GC content of L. garvieae prophages to the genomes of L. lactis phages further supports the hypothesis that these phages likely originated from the same ancestor.

Characterization of two polyvalent phages infecting Enterobacteriaceae.

Bacteriophages display remarkable genetic diversity and host specificity. In this study, we explore phages infecting bacterial strains of the Enterobacteriaceae family because of their ability to infect related but distinct hosts. We isolated and characterized two novel virulent phages, SH6 and SH7, using a strain of Shigella flexneri as host bacterium. Morphological and genomic analyses revealed that phage SH6 belongs to the T1virus genus of the Siphoviridae family. Conversely, phage SH7 was classified in the T4virus genus of the Myoviridae family. Phage SH6 had a short latent period of 16 min and a burst size of 103 ± 16 PFU/infected cell while the phage SH7 latent period was 23 min with a much lower burst size of 26 ± 5 PFU/infected cell. Moreover, phage SH6 was sensitive to acidic conditions (pH < 5) while phage SH7 was stable from pH 3 to 11 for 1 hour. Of the 35 bacterial strains tested, SH6 infected its S. flexneri host strain and 8 strains of E. coli. Phage SH7 lysed additionally strains of E. coli O157:H7, Salmonella Paratyphi, and Shigella dysenteriae. The broader host ranges of these two phages as well as their microbiological properties suggest that they may be useful for controlling bacterial populations.