Next Article in Journal
The Utility of Data Transformation for Alignment, De Novo Assembly and Classification of Short Read Virus Sequences
Next Article in Special Issue
Unprecedented Diversity of Lactococcal Group 936 Bacteriophages Revealed by Amplicon Sequencing of the Portal Protein Gene
Previous Article in Journal
In Depth Breadth Analyses of Human Blockade Responses to Norovirus and Response to Vaccination
Previous Article in Special Issue
A Cryptic Non-Inducible Prophage Confers Phage-Immunity on the Streptococcus thermophilus M17PTZA496
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Isolation and Characterization of Lactobacillus brevis Phages

1
School of Microbiology, University College of Cork, T12 YT20 Cork, Ireland
2
APC Microbiome Ireland, University College of Cork, T12 YT20 Cork, Ireland
3
Laboratory of Probiogenomics, Department of Chemistry, Life Sciences, and Environmental Sustainability, University of Parma, 43124, Parma, Italy
4
Department Microbiology and Biotechnology, Federal Research Centre of Nutrition and Food, Max Rubner-Institut, 24103, Kiel, Germany
5
Department Physiology Biochemistry and Immunology, Biomedical Research Institute, Hasselt University, B-3590 Diepenbeek, Belgium
6
HEINEKEN Global Innovation and Research, Heineken Supply Chain B.V, 2382 Zoeterwoude, The Netherlands
*
Author to whom correspondence should be addressed.
Viruses 2019, 11(5), 393; https://doi.org/10.3390/v11050393
Submission received: 26 March 2019 / Revised: 23 April 2019 / Accepted: 24 April 2019 / Published: 26 April 2019

Abstract

:
Lactobacillus brevis has been widely used in industry for fermentation purposes. However, it is also associated with the spoilage of foods and beverages, in particular, beer. There is an increasing demand for natural food preservation methods, and in this context, bacteriophages possess the potential to control such spoilage bacteria. Just a few studies on phages infecting Lactobacillus brevis have been performed to date and in the present study, we report the isolation and characterization of five virulent phages capable of infecting Lb. brevis strains. The analysis reveals a high diversity among the isolates, with members belonging to both, the Myoviridae and Siphoviridae families. One isolate, designated phage 3-521, possesses a genome of 140.8 kb, thus representing the largest Lb. brevis phage genome sequenced to date. While the isolated phages do not propagate on Lb. brevis beer-spoiling strains, phages showed activity against these strains, impairing the growth of some Lb. brevis strains. The results highlight the potential of bacteriophage-based treatments as an effective approach to prevent bacterial spoilage of beer.

1. Introduction

Lactobacillus brevis is a Gram-positive, heterofermentative lactic acid bacterium (LAB) that grows optimally at 30 °C and pH 4–6 [1]. Lb. brevis is used in the production of fermented foods [1,2]. Recently strains of this species have also been characterized as “probiotic” potentially promoting gut microbiota fitness and consumer health [2,3]. Conversely, Lb. brevis strains are also associated with food and beverage spoilage, particularly that of beer [4,5]. Beer is generally regarded as a harsh environment for microorganisms [5,6]. Indeed, the reduced availability of oxygen and nutrients coupled with the presence of an acidic environment, ethanol, carbon dioxide and hop compounds represent considerable challenges to microbial growth [5,6]. Despite the nature of beer and the array of antimicrobial compounds it contains, bacterial strains have emerged that can tolerate and grow in the presence of these hurdles [6,7]. This bacterial growth is mostly attributed to certain LAB, especially Lb. brevis and may result in the production of malodorous compounds, acidity and/or turbidity, thereby negatively impacting on the organoleptic properties of the final product [5,7,8]. Current approaches to increase the safety of beer include pasteurization, filtration, suitable materials and process packaging, strict cleaning and sanitation practices [7]. However, non-pasteurized beer products are in high demand, thus increasing the risk of microbial spoilage, for example by LAB, in particular when filtration cannot be applied [9]. The overuse of chemical sanitizers has led to an increase in biocidal-resistance of these food-spoilage bacteria [10]. Moreover, chemical sanitizers may be corrosive and/or toxic, thereby limiting the range of sanitizers that may be employed safely in industry. Various alternative strategies have been implemented to control bacterial spoilage using antimicrobials such as bacteriocins [11,12,13], and bioremediation using bacteriophages has re-appeared as a potential procedure for limiting spoilage bacteria in food and beverages [14,15,16,17]. Bacteriophages present an interesting bio-remediation approach, because they are naturally ubiquitous and specific to their bacterial host [18]. The impact of bacteriophages in preventing/limiting spoilage has been thoroughly explored in the case of food fermentation applications [16], although at the same time the prevention or limitation of bacterial spoilage of fermented beverages such as beer using bacteriophages is poorly studied [17]. To date, the genome sequences of approximately 50 Lb. brevis strains (and their associated prophages) are available on the NCBI database while only one lytic phage (SA-C12) [17] and one temperate phage (LBR48) have been described [19]. Here, we report the isolation and characterization of phages active against Lb. brevis strains in order to increase our understanding of the diversity and therapeutic potential of Lb. brevis phages.

2. Materials and Methods

2.1. Bacterial Strains and Cultivation Media

The Lactobacillus brevis strains used in this study are listed in Table 1. Bacterial stock cultures were stored in 20% glycerol at −80 °C. Bacteria and phages were cultured and/or propagated in MRS broth (Oxoid Ltd., Hampshire, UK) at 30 °C without agitation. MRS broth was supplemented with 10 mM CaCl2 (Sigma-Aldrich, St. Louis, MO, USA) where appropriate. Soft agar was prepared with MRS broth supplemented with 0.4% agar and 0.5% glycine [20].

2.2. Phage Isolation and Enrichment

Environmental samples were clarified by centrifugation at 4000× g for 10 min followed by filtration through a 0.45 µm filter (Sarsted, Nümbrecht, Germany) and stored at 4 °C until required. The filtrate was added to equal amounts of MRS broth supplemented with 10 mM CaCl2 and inoculated with an early log-phase host culture (Table 1). After incubation at 30 °C overnight, the culture was centrifuged at 4000× g for 10 min. This enrichment procedure was repeated twice. The filtered sample was then evaluated for the presence of phages active against a panel of Lb. brevis strains (Table 1). Each of the environmental samples were enriched and tested separately on each of the Lb. brevis strains listed in Table 1.

2.3. Phage Detection, Purification and Host Range Analysis

The spot test method was applied in first instance to detect the presence of phages [20]. Soft agar (4 mL) was seeded with 200 µL of fresh overnight culture and poured onto an MRS agar plate supplemented with 10 mM CaCl2 and 0.5% glycine. On the lawn of the series of Lb. brevis strains, 10 µL of the enriched samples was spotted and incubated at 30 °C overnight. The presence of phages was demonstrated by the presence of a clear zone on the plate. Presumptive positive samples were confirmed by plaque assay using the double-layer agar plate method [20]. A 10 µL volume of the appropriate phage dilution and 200 µL of Lb. brevis culture were added to 4 mL of soft agar supplemented with 10 mM CaCl2, mixed and poured onto an MRS agar plate supplemented with 10 mM CaCl2 and 0.5% glycine. The plate was incubated at 30 °C overnight and resulting plaques were enumerated. Phages were purified by single-plaque isolation using an appropriate Lb. brevis host strain. A single plaque was picked from the bacterial lawn, transferred into a tube containing 10 mL MRS broth, 10 mM CaCl2 and 1% inoculum of the propagating Lb. brevis culture. The tube was incubated at 30 °C overnight. The phage lysate was centrifuged at 4000× g for 10 min at 4 °C. The supernatant was filtered (0.45 µm) and stored at 4 °C until required. Host range studies were performed using the spot and plaque assay techniques as described above where phage lysates were tested against available Lb. brevis strains (Table 1). The presence or absence of plaque formation was recorded indicating the susceptibility of Lb. brevis strains to isolated phages. Plaques were enumerated and phage titre determined as plaque-forming units (PFU/mL).

2.4. Phage Concentration and Purification

A 2 L phage lysate was centrifuged at 5000× g for 10 min, 0.5 M NaCl was added to the supernatant and incubated for 1 h at 4 °C. The preparation was centrifuged at 5000× g for 10 min and phages were precipitated by adding 10% (w/v) polyethylene glycol 8000 (Sigma-Aldrich) and incubated overnight at 4 °C. Phages were harvested by centrifugation (as described above) and resuspended in 4 mL SM buffer (50 mM Tris-HCl pH 7.5, 100 mM NaCl, 10 mM MgSO4). Phages were extracted with chloroform (1:1 phage suspension:chloroform) applying multiple extraction steps where necessary (typically two or three times). The phage lysate was purified on a discontinuous CsCl (Sigma-Aldrich) gradient [22] and dialyzed against phage buffer (50 mM Tris-HCl, 100 mM NaCl, 8 mM MgSO4) overnight at 4 °C.

2.5. Transmission Electron Microscopy

Purified bacteriophage lysates were analysed by electron microscopy, as previously described [23]. Negative staining was performed using 2% (w/v) uranyl acetate on freshly prepared ultrathin carbon films. Grids were analysed in a Tecnai 10 transmission electron microscope (FEI Thermo Fisher Scientific, Eindhoven, The Netherlands) at an acceleration voltage of 80 kV. Micrographs were taken with a MegaView G2 charge-coupled device camera (Emsis, Muenster, Germany).

2.6. Phage DNA Extraction and Sequencing

Phage DNA was extracted using the Norgen Biotek Corp phage DNA isolation kit as per the manufacturer’s instructions (Norgen Biotek Corp., Thorold, Ontario, Canada). Phage genome sequencing was performed by GenProbio at the University of Parma, Italy. Genomes were sequenced with Illumina MiSeq Sequencing System and assembled with MIRA v4.0.2. De novo sequence assemblies and automated gene calling was performed using the MEGAnnotator pipeline [24] and assessed for predicted transfer RNA genes via tRNAscan-SE v1.2.1 [25]. Predicted open reading frames (ORFs) were determined via Prodigal v2.6 [26]. A BLASTP [27] analysis was performed to assign functional annotations to the predicted ORFs (https://blast.ncbi.nlm.nih.gov/Blast.cgi). The proposed functional annotations were further investigated by performing structural homology searches via HHpred [28] and querying the NCBI Conserved Domain Database (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi). The annotated genomes were manually inspected, edited and finalized using the Artemis visualization tool [29].

2.7. Phage Structural Proteome and Mass-Spectrometry

An aliquot (30 μL) of CsCl-purified phage sample was mixed with 10 μL of SDS loading buffer containing 50 mM β-mercaptoethanol. The structural protein profile was generated by standard Tris-glycine sodium dodecyl sulfate (SDS)–12% polyacrylamide gel electrophoresis (PAGE). Gel slices were then excised, trypsinized, and analysed using electrospray ionisation tandem mass spectrometry (ESI-MS/MS), as previously described [30].

2.8. Proteomic Tree

To study the relationship between Lactobacillus phages a proteomic tree was constructed. The genomes of the five Lb. brevis phages isolated as part of this study as well as all Lactobacillus phage genomes available on the NCBI database were downloaded. All predicted protein-encoding sequences were extracted and concatenated beginning with the ORF encoding for the small terminase subunit (TerS) [23]. The concatenated sequences were aligned using ClustalW [31]. The phylogenetic tree was constructed using the neighbour-joining method and bootstrapped employing 1000 replicates. The final tree was visualized using MEGA7 [32].

2.9. Phage Activity against Lb. brevis Beer-Spoiling Strains

To assess if the isolated phages could affect Lb. brevis beer-spoiling strains’ ability to grow, the strains were grown in MRS broth until an OD600 nm of 0.2 was reached, at which point phages were added at a MOI (Multiplicity Of Infection) of 1, along with 10 mM CaCl2. The optical density at 600 nm (OD600 nm) was recorded at 30-min intervals for 48 h to monitor the impact of the addition of phages on the growth of Lb. brevis beer-spoiling strains. A control culture was also employed where the strain was grown in the absence of phage lysate but treated identically in all other aspects, i.e., filtered MRS broth was added in place of phage lysate. Addition of calcium chloride, incubation time and temperature were identical for both scenarios.
Adsorption assays were adapted from a previously outlined protocol [33]. Briefly, strains were grown to mid-late exponential phase (OD600nm ~ 0.5), at which point they were harvested by centrifugation at 4000× g for 10 min and resuspended in 1/4-strength Ringer’s solution. Phages were added to the cells at a final titre of 106 PFU/mL followed by incubation at 30 °C for 15 min. The supernatant was retained after centrifugation and tested for the residual phage concentration by plaque assay as described above. Adsorption efficiency was calculated using the formula:
((Ci – Cr) / Ci) × 100
where Ci represents the total phage concentration used in the adsorption assay and Cr represents the residual phage concentration after the adsorption step.
The ability of phages to propagate and multiply within the host cell was also tested. Lb. brevis strains were grown to mid-late exponential phase (OD600nm ~ 0.5), at which point phages were added to the culture (T0), the mix was further incubated at 30 °C overnight (T1). The phage titre was enumerated at T0 and T1 to assess phage propagation efficiency.

2.10. Nucleotide Sequence Accession Numbers

The genome sequences of the phages isolated in this study were deposited in the GenBank database under accession numbers: 3-521: MK504444; 521B: MK504443; 3-SAC12: MK504442; SAC12B: MK504446; ATCCB: MK504445. The GenBank accession numbers of phage genome sequences applied in the proteomic tree preparation are as follows: Lactobacillus plantarum phage ATCC8014-B1: JX486087; Lactobacillus plantarum phage ATCC8014-B2: JX486088; Lactobacillus casei prophage A2: AJ251789; Lactobacillus helveticus phage AQ113: HE956704; Lactobacillus delbrueckii phage c5: EU340421; Lactobacillus casei phage J-1: KC171646; Lactobacillus delbrueckii phage JCL1032: EU409559; Lactobacillus gasseri phage kc5a: DQ320509; Lactobacillus paracasei phage Lb3381: FJ822135; Lactobacillus brevis phage LBR48: GU967410; Lactobacillus rhamnosus phage Lc-Nu: AY131267; Lactobacillus delbrueckii phage Ld3: KJ564038; Lactobacillus delbrueckii phage Ld17: KJ654037; Lactobacillus delbrueckii phage Ld25A: KJ654036; Lactobacillus delbrueckii phage Ldl1: KM514685; Lactobacillus fermentum phage LF1: HQ141410; Lactobacillus delbrueckii phage LL-H: EF455602; Lactobacillus delbrueckii phage LL-Ku: AY739900; Lactobacillus johnsonii phage Lj965: AY459535; Lactobacillus johnsonii phage Lj928: AY459533; Lactobacillus plantarum phage LP65: AY682195; Lactobacillus rhamnosus phage Lrm1: EU246945; Lactobacillus jensenii phage Lv1: EU871039; Lactobacillus gasseri phage phiadh: AJ131519; Lactobacillus casei phage phiAT3: AY605066; Lactobacillus plantarum phage phig1e: X98106; Lactobacillus delbrueckii phage phiJB: KF188409; Lactobacillus delbrueckii phage phiLdb: KF188410; Lactobacillus fermentum phage phiPYB5: GU323708; Lactobacillus casei phage PL1: KC171647; Lactobacillus brevis phage SA-C12: KU052488 and Lactobacillus plantarum phage Sha1: HQ141411.

3. Results

3.1. Phage Isolation and Host Range Profile

Lactobacillus brevis is a persistent problem in the brewing industry due to its ability to grow in, and spoil, beer. Therefore, a screen for phages capable of infecting Lb. brevis strains with potential industrial relevance was undertaken. In excess of 200 environmental samples were screened for the presence of phages active against Lb. brevis. These environmental samples included silage, fermented foods and wastewater samples, collected at different locations (Ireland, Belgium and The Netherlands) over a period of three years. Five distinct virulent phages capable of infecting one or more Lb. brevis strain(s) within our collection were isolated from two different Irish wastewater samples (collected in 2017 and 2018), purified and characterised. Two phages were isolated that infected each of two strains, namely Lb. brevis UCCLB521 and Lb. brevis SA-C12, while a further isolate was identified that targeted Lb. brevis ATCC 367. These isolates were propagated to a titre of 109 PFU/ml (except for ATCCB, where only a titre of 107 PFU/ml could be reached) and applied to a host range analysis (Table 2) against the collection of Lb. brevis strains available. This analysis highlighted the narrow host range of the isolated phages, while also highlighting the relative sensitivity of two strains. The Lb. brevis strains UCCLB521 and SA-C12 exhibited sensitivity to three and two phages, respectively. On Lb. brevis strains UCCLBBS449, UCCLB95 and RIBM 2-56, a clearing zone was observed on bacterial lawns used in the spot assay technique. However, propagation of the phages using these Lb. brevis strains as hosts was not possible (see results below).

3.2. Phage Morphology

The morphological diversity of the phage isolates was assessed by transmission electron microscopy (Figure 1). Lb. brevis phages 3-521, SAC12B and 521B possessed relatively short yet wide contractile tails and a large icosahedral head with a large complex baseplate structure at the distal end of the tail (Figure 1). These structural features are consistent with the typical attributes of Myoviridae phages [34] and revealed morphological similarity to the only virulent Lb. brevis phage identified to date, SA-C12 [17]. Lb. brevis phage 3-SAC12 possessed an icosahedral head, a defined baseplate structure and a long decorated contractile tail (Figure 1) and, therefore, also belongs to the Myoviridae family [34] and resembles the Lb. brevis temperate phage LBR48 [19]. Lb. brevis phage ATCCB was classified as a Siphoviridae phage due to the presence of a long non-contractile tail, a large icosahedral head and a discrete baseplate at the tip of the tail (Figure 1).

3.3. Lb. brevis Phages Comparative Analysis and Grouping

In order to evaluate the diversity of Lb. brevis phages and their phylogenetic links to phages of other lactobacilli, a proteomic tree was created gathering the five Lb. brevis phages characterized in this study as well as all previously sequenced Lactobacillus phages (Figure 2). The phylogenetic tree shows an interesting organization based seemingly on morphology rather than phage infecting-species. The right side of the tree displays exclusively phages belonging to the Siphoviridae family, while the left side predominantly gathered phages belonging to the Myoviridae family. It is noteworthy that Lb. brevis phages are quite diverse as they do not form a single cluster and are, in fact, spread across the phylogenetic tree with the exception of phages SAC12B and 521B, which form a clade next to the Lb. helveticus phage AQ113, a Myoviridae phage which shows similarity to phages of human gut-inhabiting species [35]. Lb. brevis phages 3-521, 3-SAC12, SA-C12, ATCCB and LBR48 all gathered closely on the tree in between Lb. plantarum phages 8014-B1, 8014-B2 and Lb. delbrueckii phage JCL1032, highlighting once again the interrelationships of the Lactobacillus phages. The relationship between these phages infecting similar host species might be explained by evolution over time from a common ancestor [23,36].

3.4. Genome Analysis

Genomic DNA of the five lytic phages was isolated and sequenced revealing significant genetic disparity between these phages. General genome characteristics of the phage isolates are summarized in Table 3. The Siphoviridae phage ATCCB possesses a genome of 80.5 kb while the Myoviridae phage genomes vary in size from ~ 41–141 kb (Table 3). The largest phage genome among the isolates is that of 3-521 with a genome of 140.8 kb, which now represents the largest known Lb. brevis phage genome sequenced to date. Interestingly, the Lb. brevis myophages SAC12B, 521B and 3-521 are more closely related to myophages of other Lactobacillus spp. harbouring a large genome size, such as Lb. casei Lb338-1 (142 kb) [37] and Lb. plantarum LP65 (131 kb) [38], than the previously characterized Lb. brevis phages (LBR48 and SA-C12) (Figure 2).
The genome of the phages investigated here display limited/no similarity to each other or to the genomes of other Lb. brevis phages, with the exception of phages 521B and SAC12B. These two phages share 97% nucleotide sequence identity (88% coverage) and their close relationship may be the result of their cohabitation within the same environment, as they were both isolated from the same wastewater sample in 2018 (Table 3). The absence of similarity with previously described Lb. brevis phages highlights the limited knowledge, and the apparent genetic diversity of these phages. The GC content of the phages is relatively low (Table 3) compared to that of the host (~ 46%), implying that they may have evolved recently to infect Lb. brevis strains. The genomes of phages 3-521, 521B and SAC12B appear complex due to their size and their high number of predicted ORFs but present a similar genome organization. The genomes were organised into four functional modules: DNA packaging, morphogenesis, DNA replication and lysis modules (Figure 3A). For phages 3-SAC12 and ATCCB, the genome organization is similar but with a lysis module preceding the replication module (Figure 3B,C). While these phage isolates are predicted to be obligatorily virulent, there are traces of temperate ancestry in some of their genomes. For example, phage 3-SAC12 possesses a predicted antirepressor-encoding gene (typically associated with lytic/lysogenic switch genomic regions) while ATCCB possesses a predicted recombinase/integrase-encoding gene [23].

3.5. Morphogenesis Module

The majority of the morphogenesis modules of Lb. brevis phages 3-521, 521B and SAC12B exhibited a high degree of synteny in the region encoding the portal protein through to the putative adsorption protein (Figure 3A). Phages 521B and SAC12B share more than 90% amino acid (aa) sequence identity, while 3-521 shares less than 50% aa sequence similarity with 521B and SAC12B across the morphogenesis module (Figure 3A). The most notable difference is the apparent insertion of an additional capsid-encoding protein in 3-521. The encoded predicted capsid protein, ORF573-521, is divided in two different proteins in 521B (ORF80521B-ORF81521B, 90% and 45% aa similarity, respectively) and in SAC12B (ORF79SAC12B-ORF80SAC12B, 26% and 83% aa similarity). In 521B and SAC12B, these protein-encoding genes are located upstream of the DNA packaging module in a divergently oriented cluster of genes of unknown function. Interestingly, it suggests the fusion of two ancestral phages into these two unique phages: 521B and SAC12B (Figure 3A).
The morphogenesis module of the Myoviridae Lb. brevis phage 3-SAC12 harbours genes encoding the phage capsid and tail structural components including a portal protein, two capsid proteins and a head-tail adaptor protein, the tail sheath protein, a tail tape measure protein (TMP), the major tail protein (MTP) and a “puncturing device” protein. This puncturing device comprises the tip of the central spike and is proposed to facilitate DNA ejection into the host cell [39]. Furthermore, at the distal tail region, there is a large organelle described as a baseplate complex that comprises three structural proteins (ORF213-SAC12, ORF223-SAC12 and ORF233-SAC12) and a protein that harbours a predicted carbohydrate binding domain (ORF243-SAC12) that we predict to bind to the host cell acting as the receptor binding protein (RBP) (Figure 3B).
The Siphoviridae phage ATCCB appears less complex in its morphogenesis module compared to the Myoviridae phages and genes encoding a portal protein, a prohead protease, a major capsid protein, a head-tail joining protein, four predicted tail proteins, a distal tail (Dit) protein, a tail fibre protein, a baseplate protein and a predicted attachment protein, assumed to be involved in host recognition and binding, were identified (Figure 3C).

3.6. Structural Proteome

The lytic Lb. brevis phages were analysed by mass spectrometry to identify their structural proteomes (Table 4). Most of the predicted proteins encoded within the morphogenesis module of the genomes of 521B, SAC12B and 3-521 were confirmed as structural proteins with the predicted portal protein, prohead protease, major capsid precursor protein, capsid protein, tail sheath protein, tail proteins, tape measure protein, tail lysin, tail fibre, baseplate proteins, putative receptor binding protein and adsorption protein; all identified as structural proteins using this approach (Table 4 and Figure 3A). The majority of the predicted structural proteins forming the capsid and the tail components were identified in 3-SAC12 and ATCCB (Table 4, Figure 3B,C). Some (presumed) structural proteins were not identified in the experimentally determined proteome, which was likely due to their small size or their low relative abundance.

3.7. Phage Activity against Lb. brevis Beer-Spoiling Strains

Phage adsorption experiments were performed in order to test the ability of the phages to recognize and bind to Lb. brevis beer-spoiling strains. Here, an adsorption efficiency higher than 50% was considered as significantly effective adsorption of the phage to the strain. Phages were tested against all Lb. brevis beer-spoiling strains and efficient phage adsorption was only observed in the cases described below (Figure 4D–F). Adsorption of the lytic phages 521B and 3-521 to their Lb. brevis host UCCLB521 showed more than 90% adsorption efficiency. Phages 521B and 3-521 were capable of high adsorption efficiencies to the Lb. brevis strains UCCLBBS449 and UCCLB95 (86.6 ± 4.7% and 98.9 ± 0.5%, respectively). Similarly, SAC12B adsorbed to its host strain SA-C12 and the beer-spoiling strain RIBM 2-56, with similar efficiencies (90.9 ± 0.9% versus 87.7 ± 0.0%). Lb. brevis strain UCCLBBS124 was not adsorbed efficiently by phages 3-521 and SAC12B (Figure 4E,F) and even if 521B showed an adsorption efficiency of 66.6 ± 7.0% on UCCLBBS124 (Figure 4D), no infection or effect of the phage was observed against this strain (Figure S1). Since most of the phages were capable of adsorbing to Lb. brevis beer-spoiling strains, experiments were performed to study their ability to affect growth of Lb. brevis strains in nutritive media (MRS broth). Lb. brevis beer-spoiling strains were grown in nutritive media until they reached an OD600nm of 0.2, at which point the relevant lytic phage showing adsorption capability (Figure 4D–F) were added to reach an MOI of 1. In some cases, the addition of the bacteriophage to the culture had a negative effect on growth of the Lb. brevis beer-spoiling strain, as the strains were not able to grow after addition of the phages even after 32 h of exposure (Figure 4A–C). Lytic phages 521B, 3-521 and SAC12B were shown to affect growth of Lb. brevis strains UCCLBBS449, UCCLB95 and RIBM 2-56, respectively (Figure 4A–C). No negative impact was observed on the growth of Lb. brevis beer-spoiling strain UCCLBBS124 even after addition of phages and was used as a control (Figure S1). Growth curves of the host strains Lb. brevis UCCLB521 and SA-C12, with and without phage treatment, are also described in Figure S1.
Lytic phages isolated as part of this study adsorb onto Lb. brevis beer-spoiling strains and negatively affect their growth. However, they failed to form visible plaques on the beer-spoiling strains, thus we aimed to evaluate the potential of these phages to propagate within the host cell. Plaque assays after the enrichment did not reveal phage propagation and multiplication within the host (data not shown) as the phage titre did not increase after the incubation period. However, while they did not infect the beer-spoiling strains, they did appear to affect the growth rate of Lb. brevis beer-spoiling strains. It did not seem that the phages propagate lytically on these strains, but the negative impact of phages on certain Lb. brevis beer-spoiling strains might be caused by a high-multiplicity phage adsorption and/or by the action of exogenous phage-encoded lysin on the bacteria [40].
The negative impact of phages on Lb. brevis strains growth presents potential for the application of such entities to control bacterial spoilage of beer. In the experiment presented above, the number of cells to which phages were added is high and most certainly exceeds levels encountered during the beer fermentation process.

4. Conclusions

In this study, the isolation and characterization of five Lb. brevis-infecting phages considerably increases knowledge of the genetic and morphological diversity of Lb. brevis phages, as only one lytic Lb. brevis phage had been isolated to date. Despite their shared host species, they show a high level of genetic diversity. Their morphology and genome size vary considerably with the largest phage isolated against Lb. brevis being that of 3-521 with a genome size of 141 kb. Some of the phages isolated as part of this study showed activity against Lb. brevis beer-spoiling strains preventing them from growing optimally, thus providing new approaches to control bacterial spoilage of beer. Indeed, such phages may be used in the future during beer fermentation to control and restrain growth of spoilage bacteria by bioremediation.
Interestingly, Lb. brevis bacteria are widely present in fermented foods, silage or microbiota; however, phages against this microorganism were not ubiquitously and easily isolated. Out of 200 environmental samples screened, only five lytic Lb. brevis phages were retrieved and only from Irish wastewater samples, indicating the rarity and the hurdle of isolating such entities. The study of Lb. brevis phages is in its infancy and many questions remain to be answered regarding their mode of action and their evolutionary strategies. For this reason, screening of phage populations from different sources is necessary to provide sufficient knowledge for their potential use in bioremediation applications.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4915/11/5/393/s1, Figure S1: Growth of (A) Lb. brevis beer-spoiling strain UCCLBBS124, (B) Lb. brevis strain UCCLB521 and (C) Lb. brevis strain SA-C12 when challenged with lytic phages (MOI = 1 when the culture reached an OD600nm of 0.2). A culture of the bacterial strain where no phage was added was used as a control.

Author Contributions

M.F. performed experiments, annotation, genomic and proteomic analysis. G.A.L. and M.V. carried out sequence assembly and annotation. H.N. and C.M.A.P.F. carried out electron microscopy analysis. J.-P.N. performed mass-spectrometry analysis. D.v.S., J.M., T.O.S. provided materials and strains. M.F., J.M., T.O.S. and D.v.S. were involved in project design and wrote the manuscript. All authors read and approved the final manuscript.

Funding

Marine Feyereisen is the recipient of an Irish Research Council Enterprise Partnership Scheme postgraduate scholarship (Ref. No. EPSPG/2015/7). Douwe van Sinderen is supported by a Principal Investigator award (Ref. No. 450 13/IA/1953) through Science Foundation Ireland (SFI). Jennifer Mahony is in receipt of the Starting Investigator Research Grant (SIRG) (Ref. No. 15/SIRG/3430) funded by Science Foundation Ireland (SFI).

Acknowledgments

The authors are grateful to Dagmar Matoulková of the Research Institute of Brewing and Malting, Prague, who kindly provided Lb. brevis RIBM 2-56. The authors also acknowledge the technical help of Angela Back (Max Rubner-Institut, Kiel) in the electron microscopic analysis.

Conflicts of Interest

The authors declare that T.O.S. is an employee of Heineken.

References

  1. Vos, P.; Garrity, G.; Jones, D.; Krieg, N.R.; Ludwig, W.; Rainey, F.A.; Schleifer, K.-H.; Whitman, W. Bergey’s Manual of Systematic Bacteriology Volume 3: The Firmicutes; Springer Science and Business Media: Athens, GA, USA, 2011; pp. 464–511. [Google Scholar]
  2. Salvetti, E.; Torriani, S.; Felis, G.E. The genus Lactobacillus: A taxonomic update. Probiotics Antimicrob. Proteins 2012, 4, 217–226. [Google Scholar] [CrossRef]
  3. Felis, G.E.; Dellaglio, F. Taxonomy of Lactobacilli and Bifidobacteria. Curr. Issues Intest. Microbiol. 2007, 8, 44–61. [Google Scholar]
  4. Bergsveinson, J.; Pittet, V.; Ewen, E.; Baecker, N.; Ziola, B. Genome sequence of rapid beer-spoiling isolate Lactobacillus brevis BSO 464. Genome Announce. 2015, 3, e01411–e01415. [Google Scholar] [CrossRef]
  5. Sakamoto, K.; Konings, W.N. Beer spoilage bacteria and hop resistance. Int. J. Food Microbiol. 2003, 89, 105–124. [Google Scholar] [CrossRef]
  6. Suzuki, K.; Iijima, K.; Sakamoto, K.; Sami, M.; Yamashita, H. A review of hop resistance in beer spoilage lactic acid bacteria. J. Institute Brewing 2006, 112, 173–191. [Google Scholar] [CrossRef]
  7. Vaughan, A.; O’Sullivan, T.; van Sinderen, D. Enhancing the microbiological stability of malt and beer—A review. J. Institute Brewing 2005, 111, 355–371. [Google Scholar] [CrossRef]
  8. Bokulich, N.A.; Bamforth, C.W. The Microbiology of malting and brewing. Microbiol. Mol. Biol. Rev. 2013, 77, 157–172. [Google Scholar] [CrossRef]
  9. Back, W.; Leibhard, M.; Bohak, I. Flash pasteurisation-membrane filtration. Biological security in comparison. Brauwelt Int. 1990, 12, 42–49. [Google Scholar]
  10. Davidson, P.M.; Harrison, M.A. Resistance and adaptation to food antimicrobials, sanitizers, and other process controls. Food Technol. 2002, 56, 69–78. [Google Scholar]
  11. Basanta, A.; Sánchez, J.; Gómez-Sala, B.; Herranz, C.; Hernández, P.E.; Cintas, L.M. Antimicrobial activity of Enterococcus faecium L50, a strain producing enterocins L50 (L50A and L50B), P and Q, against beer-spoilage lactic acid bacteria in broth, wort (hopped and unhopped), and alcoholic and non-alcoholic lager beers. Int. J. Food Microbiol. 2008, 125, 293–307. [Google Scholar] [CrossRef]
  12. Gil, G.; del Monaco, S.; Cerrutti, P.; Galvagno, M. Selective antimicrobial activity of chitosan on beer spoilage bacteria and brewing yeasts. Biotechnol. Lett. 2004, 26, 569–574. [Google Scholar] [CrossRef]
  13. Vaughan, A.; Rouse, S.; van Sinderen, D. Investigating the antimicrobial efficacy of a Lactococcal bacteriocin for the development of microbiologically stable beer. J. Institute Brewing 2004, 110, 181–188. [Google Scholar] [CrossRef]
  14. Mahony, J.; McAuliffe, O.; Ross, R.P.; van Sinderen, D. Bacteriophages as biocontrol agents of food pathogens. Curr. Opin. Biotechnol. 2011, 22, 157–163. [Google Scholar] [CrossRef]
  15. Coffey, B.; Mills, S.; Coffey, A.; McAuliffe, O.; Ross, R.P. Phage and their lysins as biocontrol agents for food safety applications. Annual Rev. Food Sci. Technol. 2010, 1, 449–468. [Google Scholar] [CrossRef]
  16. Endersen, L.; O’Mahony, J.; Hill, C.; Ross, R.P.; McAuliffe, O.; Coffey, A. phage therapy in the food industry. Annual Rev. Food Sci. Technol. 2014, 5, 327–349. [Google Scholar] [CrossRef]
  17. Deasy, T.; Mahony, J.; Neve, H.; Heller, K.J.; van Sinderen, D. Isolation of a virulent Lactobacillus brevis phage and its application in the control of beer spoilage. J. Food Protect. 2011, 74, 2157–2161. [Google Scholar] [CrossRef]
  18. Brüssow, H.; Desiere, F. Comparative phage genomics and the evolution of Siphoviridae: Insights from dairy phages. Mol. Microbiol. 2001, 39, 213–222. [Google Scholar] [CrossRef]
  19. Jang, S.H.; Yoon, B.H.; Chang, H.I. Complete nucleotide sequence of the temperate bacteriophage LBR48, a new member of the family Myoviridae. Archi. Virol. 2011, 156, 319–322. [Google Scholar] [CrossRef]
  20. Svensson, U.; Christiansson, A. Methods for phage monitoring. Bulletin Inter. Dairy Federation 1991, 263, 29–39. [Google Scholar]
  21. Makarova, K.; Slesarev, A.; Wolf, Y.; Sorokin, A.; Mirkin, B.; Koonin, E.; Pavlov, A.; Pavlova, N.; Karamychev, V.; Polouchine, N.; et al. Comparative genomics of the lactic acid bacteria. Proc. Nat. Acad. Sci. USA 2006, 103, 15611–15616. [Google Scholar] [CrossRef]
  22. Sambrook, J.; Russell, D.W. Purification of bacteriophage λ particles by isopycnic centrifugation through CsCl gradients. Cold Spring Harbor Protocols 2006, 1, pdb-prot3968. [Google Scholar] [CrossRef]
  23. Casey, E.; Mahony, J.; Neve, H.; Noben, J.-P.; Dal Bello, F.; van Sinderen, D. Novel phage group infecting Lactobacillus delbrueckii subsp. lactis, as revealed by genomic and proteomic analysis of bacteriophage Ldl1. Appl. Environ. Microbiol. 2015, 81, 1319–1326. [Google Scholar] [CrossRef]
  24. Lugli, G.A.; Milani, C.; Mancabelli, L.; van Sinderen, D.; Ventura, M. MEGAnnotator: A user-friendly pipeline for microbial genomes assembly and annotation. FEMS Microbiol. Lett. 2016, 363. [Google Scholar] [CrossRef]
  25. Lowe, T.M.; Eddy, S.R. tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25, 955–964. [Google Scholar] [CrossRef] [PubMed]
  26. Hyatt, D.; Chen, G.-L.; LoCascio, P.F.; Land, M.L.; Larimer, F.W.; Hauser, L.J. Prodigal: Prokaryotic gene recognition and translation initiation site identification. BMC Bioinformat. 2010, 11, 119. [Google Scholar] [CrossRef]
  27. Altschul, S.F.; Madden, T.L.; Schaffer, A.A.; Zhang, J.; Zhang, Z.; Miller, W.; Lipman, D.J. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res. 1997, 25, 3389–3402. [Google Scholar] [CrossRef]
  28. Soding, J.; Biegert, A.; Lupas, A.N. The HHpred interactive server for protein homology detection and structure prediction. Nucleic Acids Res. 2005, 33, W244–W248. [Google Scholar] [CrossRef]
  29. Rutherford, K.; Parkhill, J.; Crook, J.; Horsnell, T.; Rice, P.; Rajandream, M.A.; Barrell, B. Artemis: sequence visualization and annotation. Bioinformatics 2000, 16, 944–945. [Google Scholar] [CrossRef]
  30. Holtappels, M.; Vrancken, K.; Schoofs, H.; Deckers, T.; Remans, T.; Noben, J.-P.; Valcke, R. A comparative proteome analysis reveals flagellin, chemotaxis regulated proteins and amylovoran to be involved in virulence differences between Erwinia amylovora strains. J. Proteomics 2015, 123, 54–69. [Google Scholar] [CrossRef]
  31. Thompson, J.D.; Gibson, T.J.; Higgins, D.G. Multiple sequence alignment using ClustalW and ClustalX. Curr. Protocols Bioinformat. 2003, 1, 2–3. [Google Scholar] [CrossRef] [PubMed]
  32. Kumar, S.; Nei, M.; Dudley, J.; Tamura, K. MEGA: A biologist-centric software for evolutionary analysis of DNA and protein sequences. Briefings Bioinformat. 2008, 9, 299–306. [Google Scholar] [CrossRef] [PubMed]
  33. Collins, B.; Bebeacua, C.; Mahony, J.; Blangy, S.; Douillard, F.P.; Veesler, D.; Cambillau, C.; van Sinderen, D. Structure and functional analysis of the host recognition device of lactococcal phage Tuc2009. J. Virol. 2013, 87, 8429–8440. [Google Scholar] [CrossRef] [PubMed]
  34. Matthews, R.E.F.; Maurin, J. Classification and nomenclature of viruses. Ann. Microbiol. 1979, 1, 133–136. [Google Scholar] [CrossRef]
  35. Zago, M.; Scaltriti, E.; Rossetti, L.; Guffanti, A.; Armiento, A.; Fornasari, M.E.; Grolli, S.; Carminati, D.; Brini, E.; Pavan, P.; et al. Characterization of the genome of the dairy Lactobacillus helveticus bacteriophage ΦAQ113. Appl. Environment. Microbiol. 2013, 79, 4712–4718. [Google Scholar] [CrossRef]
  36. Tuohimaa, A.; Riipinen, K.-A.; Brandt, K.; Alatossava, T. The genome of the virulent phage Lc-Nu of probiotic Lactobacillus rhamnosus, and comparative genomics with Lactobacillus casei phages. Archi. Virol. 2006, 151, 947. [Google Scholar] [CrossRef] [PubMed]
  37. Alemayehu, D.; Ross, R.P.; O’Sullivan, O.; Coffey, A.; Stanton, C.; Fitzgerald, G.F.; McAuliffe, O. Genome of a virulent bacteriophage Lb338–1 that lyses the probiotic Lactobacillus paracasei cheese strain. Gene 2009, 448, 29–39. [Google Scholar] [CrossRef]
  38. Chibani-Chennoufi, S.; Dillmann, M.-L.; Marvin-Guy, L.; Rami-Shojaei, S.; Brüssow, H. Lactobacillus plantarum bacteriophage LP65: A new member of the SPO1-like genus of the family Myoviridae. J. Bacteriol. 2004, 186, 7069–7083. [Google Scholar] [CrossRef]
  39. Rossmann, M.G.; Mesyanzhinov, V.V.; Arisaka, F.; Leiman, P.G. The bacteriophage T4 DNA injection machine. Curr. Opin. Struct. Biol. 2004, 14, 171–180. [Google Scholar] [CrossRef]
  40. Abedon, S.T. Lysis from without. Bacteriophage 2011, 1, 46–49. [Google Scholar] [CrossRef]
Figure 1. Electron micrographs of lytic Lb. brevis phages 3-521 (A), 521B (B), 3-SAC12 (C), SAC12B (D) and ATCCB (E). Head diameter, tail length and width are also indicated, where “n” represents the number of phage particles measured. For phage 521B and SAC212B, only few particles were detected with original extended tail sheaths (i.e., 2–4 particles). Tail lengths of phages 3-521, 521-B and SAC12B are also including the complex baseplate structures.
Figure 1. Electron micrographs of lytic Lb. brevis phages 3-521 (A), 521B (B), 3-SAC12 (C), SAC12B (D) and ATCCB (E). Head diameter, tail length and width are also indicated, where “n” represents the number of phage particles measured. For phage 521B and SAC212B, only few particles were detected with original extended tail sheaths (i.e., 2–4 particles). Tail lengths of phages 3-521, 521-B and SAC12B are also including the complex baseplate structures.
Viruses 11 00393 g001
Figure 2. Proteomic tree of all Lactobacillus phages sequenced to date. Colour coding indicates the host species for each phage. Black circles indicate Myoviridae phages, while white diamonds indicate Siphoviridae phages.
Figure 2. Proteomic tree of all Lactobacillus phages sequenced to date. Colour coding indicates the host species for each phage. Black circles indicate Myoviridae phages, while white diamonds indicate Siphoviridae phages.
Viruses 11 00393 g002
Figure 3. Genomic organisation of lytic Lb. brevis phages 3-521, 521B and SAC12B (A), 3-SAC12 (B) and ATCCB (C). The scale at the bottom of genomes is in base pairs. Each arrow represents an ORF, with the colour representing the putative function of the encoded protein. Confirmed structural protein-encoding genes from mass spectrometry analysis are also highlighted (bold outline). TerS. Small terminase subunit, TerL. Large terminase subunit.
Figure 3. Genomic organisation of lytic Lb. brevis phages 3-521, 521B and SAC12B (A), 3-SAC12 (B) and ATCCB (C). The scale at the bottom of genomes is in base pairs. Each arrow represents an ORF, with the colour representing the putative function of the encoded protein. Confirmed structural protein-encoding genes from mass spectrometry analysis are also highlighted (bold outline). TerS. Small terminase subunit, TerL. Large terminase subunit.
Viruses 11 00393 g003
Figure 4. (AC) Growth of Lb. brevis beer-spoiling strains when challenged with lytic phages (MOI = 1 when the culture reached an OD600 nm of 0.2). A culture of the bacterial strain where no phage was added was used as a control. (DF) Adsorption assays of lytic phages 521B, 3-521 and SAC12B onto Lb. brevis strains. Respective Lb. brevis host strains were used as positive control while Lb. brevis UCCLBBS124 was used as a negative control.
Figure 4. (AC) Growth of Lb. brevis beer-spoiling strains when challenged with lytic phages (MOI = 1 when the culture reached an OD600 nm of 0.2). A culture of the bacterial strain where no phage was added was used as a control. (DF) Adsorption assays of lytic phages 521B, 3-521 and SAC12B onto Lb. brevis strains. Respective Lb. brevis host strains were used as positive control while Lb. brevis UCCLBBS124 was used as a negative control.
Viruses 11 00393 g004
Table 1. Lactobacillus brevis strains used for phage isolation and characterisation.
Table 1. Lactobacillus brevis strains used for phage isolation and characterisation.
Lactobacillus brevis StrainsIsolation Source
ATCC367 [21]Silage
UCCLBBS124Beer
UCCLB521Brewery
UCCLB556Brewery
SA-C12Silage
UCCLBBS449Beer
UCCLB94Beer
UCCLB95Beer
RIBM 2-56Beer
Table 2. Lb. brevis phage host range analysis.
Table 2. Lb. brevis phage host range analysis.
Lb. brevis Non-Beer Spoiling StrainsLb. brevis Beer-Spoiling Strains
ATCC367UCCLB521UCCLB556SA-C12UCCLBBS124UCCLBBS449UCCLB94UCCLB95RIBM 2-56
Phage3-521-+ *-----~-
521B-+ *---~---
3-SAC12---+ *-----
SAC12B-+-+ *----~
ATCCB+ *--------
+: strain susceptible to phage infection; -: strain resistant to phage infection; ~: clearing zone was observed; *: host strain.
Table 3. General characteristics of Lb. brevis phages.
Table 3. General characteristics of Lb. brevis phages.
PhageSample (Date)Isolation SourceGenome Size (bp)ORFsGC Content (%)% nt Identity
(% coverage)
Myoviridae3-521S1 (2017)Wastewater (Ireland)140,81615536.93
521BS2 (2018)Wastewater (Ireland)136,44218832.2797 (88) with SAC12B
SAC12BS2 (2018)Wastewater (Ireland)136,60819132.4197 (88) with 521B
3-SAC12S1 (2017)Wastewater (Ireland)41,2926140.01
SiphoviridaeATCCBS2 (2018)Wastewater (Ireland)80,5389630.80
Table 4. Structural proteins extracted from purified phage particles by ESI-MS/MS. A minimum of two independent unique peptides or 5% coverage were used as threshold values.
Table 4. Structural proteins extracted from purified phage particles by ESI-MS/MS. A minimum of two independent unique peptides or 5% coverage were used as threshold values.
PhageORFPutative FunctionNo. of PeptidesSequence Coverage (%)
521B80Probable capsid protein829.4
81Probable capsid protein928.4
86Structural protein316.1
88Lipoprotein550.8
106Structural protein437.9
121Portal protein1228.3
122Structural protein217.4
123Caudovirus prohead protease420.8
125Major capsid protein precursor1959.7
128Capsid protein316.8
130Gp91835
132Major tail sheath protein1640.4
133Tail protein559.9
136Tape measure protein2831.2
137Tail lysin1215.1
138Structural component of the tail fibre810
140Structural protein215.1
141Structural protein542.2
142Baseplate protein323.1
143Baseplate J-like protein615.1
144Baseplate protein79.4
146Tail protein1534.1
147Putative adsorption protein922.2
156DNA starvation/stationary phase protein648
185Structural protein338.5
3-52110dUTP diphosphatase 29.9
19Zn-dependent protease523.7
52Portal protein37.9
53Prohead protease18.3
55Major capsid protein1951.1
57Phage capsid and scaffold1918.6
60Structural protein420.1
65Tail sheath protein1425
66Putative tail protein653.3
69Tape measure protein1417.4
70Tail lysin protein1621.1
71gp67321.7
72Structural protein315.3
76Baseplate protein43
77Structural protein218
78Tail protein2224.5
79Tail associated protein2320.1
98Nucleoside 2-deoxyribosyltransferase212.4
101Structural protein16.1
106Structural protein762
108Tail protein118
117Adenyltransferase616.9
119ADP-ribose pyrophosphatase 15.3
124Structural protein615.7
126AAA superfamily ATPase926.2
128Phosphatase 45.4
3-SAC121Terminase small subunit16.6
3Portal protein1331.4
5Major capsid protein426.3
6Capsid protein926.7
10Putative head-tail adaptor19.8
11Structural protein19
12Sheath protein311.8
13Structural protein222.4
17Structural protein413.6
ATCCB70Baseplate protein69.2
79Major tail protein650.2
86Major capsid protein825.9
87Prohead protease27.5

Share and Cite

MDPI and ACS Style

Feyereisen, M.; Mahony, J.; Lugli, G.A.; Ventura, M.; Neve, H.; Franz, C.M.A.P.; Noben, J.-P.; O’Sullivan, T.; van Sinderen, D. Isolation and Characterization of Lactobacillus brevis Phages. Viruses 2019, 11, 393. https://doi.org/10.3390/v11050393

AMA Style

Feyereisen M, Mahony J, Lugli GA, Ventura M, Neve H, Franz CMAP, Noben J-P, O’Sullivan T, van Sinderen D. Isolation and Characterization of Lactobacillus brevis Phages. Viruses. 2019; 11(5):393. https://doi.org/10.3390/v11050393

Chicago/Turabian Style

Feyereisen, Marine, Jennifer Mahony, Gabriele A. Lugli, Marco Ventura, Horst Neve, Charles M. A. P. Franz, Jean-Paul Noben, Tadhg O’Sullivan, and Douwe van Sinderen. 2019. "Isolation and Characterization of Lactobacillus brevis Phages" Viruses 11, no. 5: 393. https://doi.org/10.3390/v11050393

APA Style

Feyereisen, M., Mahony, J., Lugli, G. A., Ventura, M., Neve, H., Franz, C. M. A. P., Noben, J. -P., O’Sullivan, T., & van Sinderen, D. (2019). Isolation and Characterization of Lactobacillus brevis Phages. Viruses, 11(5), 393. https://doi.org/10.3390/v11050393

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop