Next Article in Journal
Spontaneous Calcium Bursts Organize the Apical Actin Cytoskeleton of Multiciliated Cells
Previous Article in Journal
The Role of IL-17A in Mediating Inflammatory Responses and Progression of Neurodegenerative Diseases
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Plasmid-Encoded Surface Polysaccharide Partly Blocks Ceduovirus Infection in Lactococci

by
Claudia Rendueles
1,
Javier Nicolás Garay-Novillo
2,
Martin Holm Rau
3,
Paula Gaspar
3,
José Ángel Ruiz-Masó
4,
Jennifer Mahony
5,6,
Ana Rodríguez
1,
José Luis Barra
2,
Gloria del Solar
4 and
Beatriz Martínez
1,*
1
Instituto de Productos Lacteos de Asturias (IPLA), CSIC, C/Francisco Pintado Fe, 26, 33011 Oviedo, Asturias, Spain
2
Departamento de Quimica Biologica Ranwel Caputto, CIQUIBIC-CONICET, Facultad de Ciencias Quimicas, Universidad Nacional de Cordoba, Cordoba X5000HUA, Argentina
3
R&D, Microbe and Culture Research, Novonesis A/S, 2970 Hørsholm, Denmark
4
Centro de Investigaciones Biológicas Margarita Salas (CIB), CSIC, c/ Ramiro de Maetzu, 9, 28040 Madrid, Spain
5
School of Microbiology, University College Cork, T12 K8AF Cork, Ireland
6
APC Microbiome Ireland, University College Cork, T12 YT20 Cork, Ireland
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2025, 26(6), 2508; https://doi.org/10.3390/ijms26062508
Submission received: 20 February 2025 / Revised: 5 March 2025 / Accepted: 9 March 2025 / Published: 11 March 2025
(This article belongs to the Special Issue Exploring Phage–Host Interactions: Novel Findings and Perspectives)

Abstract

:
Bacteriophages (or phages) remain the leading cause of failure in dairy fermentations. Thereby, phage-resistant Lactococcus lactis and Lactococcus cremoris dairy starters are in continuous demand. In this work, our goal was to identify phage defense mechanisms against ceduoviruses encoded by two wild isolates of dairy origin named L. lactis IPLA517 and IPLA1064. These strains were previously subjected to experimental evolution to select derivatives that are resistant to the bacteriocin Lcn972. It was observed that the Lcn972R derivatives became sensitive to phage infection; however, the underlying mechanism was not defined. The long-read sequencing technologies applied in this work reveal that all of the Lcn972R derivatives shared the loss of a 41 kb endogenous plasmid (p41) that harbors a putative exopolysaccharide (EPS) gene cluster with significant homology to one described in Lactococcus garvieae. Using a CRISPR-Cas9-based approach, p41 was selectively cured from L. lactis IPLA1064. Phage infection assays with three ceduoviruses demonstrated that curing p41 restored phage sensitivity at levels comparable to the Lcn972R-IPLA1064 derivatives. Phage adsorption to Δp41 cells was also increased, consistent with the hypothesis of EPS production hindering access to the phage receptor protein Pip. Our results reinforce the role of EPSs in protecting Lactococcus against phage infection, a phenomenon that is rarely reported for ceduoviruses. Moreover, the results also exemplify the likely horizontal gene transfer that can occur between L. lactis and L. garvieae in a dairy environment.

Graphical Abstract

1. Introduction

Lactococcus lactis and Lactococcus cremoris are key microorganisms in the production of fermented dairy products as starter cultures. Their significance extends beyond lactic acid production, encompassing the synthesis of metabolites that contribute to flavor development and biopreservation, thereby improving the overall quality and safety of the final products [1]. In this context, L. lactis must be able to tolerate adverse environmental conditions, including temperature fluctuations, elevated concentrations of salts, and, of particular interest in this study, the presence of bacteriophages [2].
As one of the main causes of failed fermentations or low-quality products, phages infecting L. lactis have been extensively studied. Lactococcal phages that are frequently associated with fermentation failures belong to Skunavirus (previously described as the 936 group), Ceduovirus (formerly known as the c2 group), and the P335 phages. [3,4]. The bacterial anti-phage arsenal, defined as genes or gene clusters that confer full or partial resistance to phage infection in the host bacterium, is of significant importance. These defense systems can be encoded in the chromosome or in mobile genetic elements and, in some cases, a distinction is made between active and passive resistance mechanisms [5]. Traditionally, numerous phage defense systems have been described in Lactococcus, with abortive infection (Abi) systems being the most ubiquitous. However, others such as restriction modification (R/M) and superinfection exclusion (Sie) systems are also frequently found [6]. Contrary to other bacterial genera, the CRISPR-Cas system does not appear to be a major component of the anti-phage arsenal in Lactococcus, as only a few CRISPR-Cas systems have been described in this genus [7]. Moreover, although its mode of action remains unclear in the majority of cases, a wide range of novel phage defense systems have emerged. In Lactococcus, systems such as PARIS, CBASS, or Lamassu have recently been described [8]. Additionally, passive mechanisms that hinder access to phage receptors (e.g., production of exopolysaccharides), thereby reducing the phage adsorption, may also exist [5,9]. As for phage receptors, cell wall polysaccharides (CWPSs) have been identified as the receptors for skunaviruses and P335-like phages [10]. An EPS has also been shown to be the phage receptor of some P335 phages [11]. In contrast, ceduoviruses bind reversibly to an unknown saccharide motif and irreversibly to the phage infection protein (Pip) or Pip-like proteins such as YjaE [12].
In previous adaptive evolution experiments carried out in our laboratory, in which L. lactis was exposed to the cell-wall-active bacteriocin lactococcin Lcn972 [13], it was observed that some of the evolved derivatives that developed resistance to Lcn972 (Lcn972R) became sensitive to phages CHPC1130 and CHPC1183, both classified as Ceduovirus members. Upon analysis of their draft genomes, different mutations were found but a common pattern could not be identified that explained their phage-sensitive phenotype. Nonetheless, a single contig sequence, harboring a putative EPS biosynthetic gene cluster and flanked by insertion sequences, was present in the parental strains and was missing in the Lcn972R phage-sensitive clones. It was hypothesized that plasmid-borne EPS production may be responsible for the phage-resistant phenotype observed in the parental strains based on non-targeted curing experiments [13]. However, the presence of additional mutations hindered the ability to directly correlate the presence of a particular plasmid with the resistant phenotype. Thus, our goal in the present work was to demonstrate the role of the putative plasmid in preventing phage infection. Long-read sequencing technologies led to the identification of p41, a 41 kb plasmid-encoding the EPS biosynthetic gene cluster. Selectively curing p41 restored sensitivity to ceduoviruses, confirming the role of p41 in phage resistance.

2. Results

2.1. Chromosome- and Plasmid-Encoded Phage-Resistance Mechanisms in L. lactis IPLA517 and L. lactis IPLA1064

Hybrid genome assemblies of the phage-resistant parental strains L. lactis IPLA517 and IPLA1064 (Table 1) were firstly conducted to determine their plasmid content and recognizable phage-resistance mechanisms. Genome sequence analysis confirmed the close relatedness of the two parental strains, with the main difference being their plasmid content. These results agree with their common origin and the plasmid extraction analysis as previously described [14].
After manual editing of the putative plasmid sequences up to seven plasmids were identified whose main features are summarized in Table 2. Both strains share the plasmids p6, p41, p45, and p66, whereas p20 was found exclusively in L. lactis IPLA517. Plasmid p17 from L. lactis IPLA1064 is identical to p15 from IPLA517 but with a 1.7 kb insertion, likely the result of a recombination event. It is noteworthy that, except for p66, all of them have highly related RepB proteins (Figure S1).
These plasmids encode important traits in the adaptation of Lactococcus to the dairy niches. Plasmid p66 harbors the lac and opp operons involved in the lactose utilization and oligopeptide transport, respectively, which are frequently present in the same plasmid. This plasmid also codes for the PepO and PepF peptidases. An additional peptidase (PepT) is encoded in the p15 and p17 plasmids and, as usual, the proteinase PrtP-encoding gene is separated from the opp operon and found in p45. Other features commonly found in lactococcal plasmids include a complete Type I restriction/modification system in p45, as well as a single subunit of another R/M system in p6, a cadmium resistance gene in p20, along with functions involved in cell surface modifications like a putative pilus-encoding gene in p20 and an EPS gene cluster in p41 (Table 2).
According to the search with PADLOC and DefenseFinder (Table 3), the anti-phage arsenal in L. lactis IPLA517 and IPLA1064 is basically composed of two R/M systems and one Abi system. With the exception of the Type I R/M system, which is plasmid-encoded (p45), the remaining systems Type IV R/M and AbiH are located in the chromosome. Some other genes in the chromosome were shown to have homology to other genes that may form part of resistance systems recently described, but complete systems were not found. According to PHASTEST, both strains share one incomplete and two intact prophages in their chromosomes (Table 4), as frequently observed in strains of dairy origin [17]. Lastly, sequence comparisons of the CWPS gene cluster of both strains revealed a high identity level to that of L. lactis UC06 (Figure S2). The CWPSs from this strain belong to the C5 subtype whose structure is formed by repeating phospho-pentasaccharide subunits that contain arabinitol residues, rather uncommon in the lactococcal pellicle [18].

2.2. Loss of an EPS-Encoding 41 kb Plasmid (p41) After Experimental Evolution Under Lcn972 Pressure

Variant analyses confirmed the previously described unique pattern of mutations in the Lcn972R mutants L. lactis IPLA517-B5, IPLA517-C6, IPLA1064-C11, and IPLA1064-E11 (Table 1), which were predominantly found in the two-component system TCS-G, dlt and pyrH genes, among others [13]. Nonetheless, the results also reveal that all of them had lost the plasmid p41, while the plasmid p20, exclusively present in L. lactis IPLA517, was absent in both of the Lcn972R mutants IPLA517-B5 and IPLA517-C6. No alteration of the prophage content or the phage-resistance systems found in the parental strains was identified, pointing to the loss of p41 as the most plausible explanation for their phage susceptible phenotype.
p41 is a 41.3 kb plasmid with two distinct domains (Figure 1A). One harbors the replication functions (repB and repA) and mobilization genes (mobC), while the other contains a 22.9 kb putative EPS biosynthetic gene cluster flanked by insertion sequences. BLASTN analyses of this putative cluster revealed a high degree of identity (99%) with a similar cluster present in a plasmid from Lactococcus garvieae M14 (GenBank accession no. CCXC01000006.1), isolated from Algerian fermented milk, and a certain degree of identity (40–90%) within specific genes from other Lactococcus eps clusters, as shown in Figure 1B. Among them, L. lactis and L. cremoris are the most represented species, followed by L. garvieae and L. petauri. Although some are chromosomally encoded, most of the eps clusters with identity to the one in p41 are located within relatively large plasmids, ranging in size from 20 to 63 kb. The eps clusters exhibit extensive variability in terms of their size and gene content with both the 5’ and 3’ ends highly conserved, while the region likely determining the EPS decoration is unique (Figure 1B). All CDSs have the same orientation, except for the last orfY in the 3’ end. The 5’ conserved region is composed of six CDSs with regulatory and chain-length determination functions, while the 3’ region has another two CDSs of unknown function. The variable region of the p41-eps cluster consists of eighteen CDSs, most of which encode putative glycosyltransferases, acetyltransferases, and enzymes likely involved in sugar-nucleotide biosynthesis. Genes predicted to encode elongation and transport functions suggest that this EPS would be synthesized via the Wzx/Wzy pathway (Figure 1).

2.3. Selective Curing of p41 Restores Phage Susceptibility

In order to demonstrate the role of the p41 plasmid in phage resistance, attempts to transfer the plasmid into phage-susceptible hosts were initially made. However, although the plasmid seemed to be easily mobilizable by conjugation, it was not stable in the plasmid-free L. cremoris MG1614 and L. lactis IL1403 used as receptors, and growing for a few generations resulted in the loss of the plasmid. Therefore, we opted to selectively cure p41 from L. lactis IPLA1064 using a CRISPR-Cas approach. This strain was chosen because of the lower number of plasmids compared to IPLA517 (Table 2).
A p41-specific 20bp-sgRNA plasmid (pILCp41LC9) was generated in E. coli and introduced into L. lactis IPLA1064 by electroporation. This plasmid targets the glycosyltransferase gene IPLA1064_002441 and contains the cas9 gene under the control of an acid-inducible weak promoter, which was induced by growing the transformants at pH 5.5 for 45 generations, until the plasmid p41 was eliminated, as judged by plasmid-specific PCR assays (Figure 2). Subsequently, the CRISPR-Cas9 pILCp41LC9 plasmid was removed by growing the transformants at pH 7 without antibiotic selection for 150 generations. In this way, a derivative named L. lactis IPLA1064-Δp41, which lacks the p41 plasmid, was obtained. Notably, additional clones that had lost plasmid p41 along with another endogenous plasmid were also detected during the curing process at pH 5.5.
The impact of losing p41 on phage infection was first screened using L. lactis IPLA1064 and its Lcn972R derivatives as hosts in double-layer agar assays. Apart from the phages CHPC1130 and CHPC1183 described previously [13], the phage c2 was also able to form clear lysis plaques on the Lcn972R derivatives. In terms of the efficiency of plaquing (EOP), phage c2 showed the largest difference in phage infectivity among the Lcn972R derivatives and the parental L. lactis IPLA1064, as the EOP of c2 on L. lactis IPLA1064 was 0.0003 ± 0.0004, while it was of 0.001 ± 0.0004 and 0.005 ± 0.002 for CHPC1130 and CHPC1183, respectively (Figure 3). Consistently, the Lcn972R derivatives lysed when infected in broth with c2 at different multiplicities of infection (MOIs) and the percentage of inhibition was always higher, while lysis did not occur in IPLA1064 (Figure S3).
When infectivity was tested on L. lactis IPLA1064-Δp41 (Figure 3A), the EOP values of phage c2 and CHPC1130 resembled those of the phage-susceptible L. lactis IPLA1064-E11. Similarly, the size and appearance of the lysis plaques were comparable, as illustrated in Figure 3B for the phage CHPC1130. Infectivity of CHPC1183 on L. lactis IPLA1064-Δp41 was not fully restored (Figure 3A), likely meaning that other mutations present in L. lactis IPLA1064-E11 might promote infection by this phage. Overall, these results establish a link between the presence of p41 and a hindered phage infection.

2.4. Phage Adsorption Is Impaired in p41-Bearing Cells

Considering that the main feature of p41 is an EPS biosynthetic gene cluster and that ceduoviruses make use of a protein receptor, we postulated that phage adsorption should be impaired in L. lactis IPLA1064. To test this hypothesis, phage adsorption experiments were carried out with phage c2. As shown in Figure 4, significant differences were noted within the tested lactococcal strains. The percentage of adsorbed phages to Lcn972R or L. lactis IPLA1064-Δp41 cells approximately doubled that of the parental L. lactis IPLA1064, reinforcing the role of the p41-encoded EPS in hindering access to the phage receptor.

3. Discussion

In this work, we present evidence to establish the basis for the phage-resistance phenotype observed in L. lactis dairy isolates, whose Lcn972R-evolved clones became susceptible to c2-like phages. To this end, their genomes were scrutinized for phage-resistance mechanisms, their plasmid complement was analyzed and, finally, using a CRISPR/Cas9 approach, the phage-resistant phenotype was linked to the presence of the p41 plasmid harboring an EPS biosynthetic cluster.
Genome analyses facilitated the identification of the plasmid content of L. lactis IPLA517 and IPLA1064. Consistent with other dairy lactococci that often possess a large number of plasmids (up to 12), these strains harbor 6 and 5 plasmids, respectively, that encode for non-essential but beneficial features. Co-existence of plasmids with highly related RepB proteins has been previously reported in L. lactis, suggesting that even minor alterations in the origins of replication or the replication proteins are sufficient to prevent the incompatibility of plasmids. [19,20,21].
Regarding their anti-phage arsenal, in addition to the plasmid-encoded R/M system in plasmid p45, we found an Abi system and another R/M system in the chromosome. Incomplete R/M systems such as the one encoded by the p6 plasmid are also frequently found and may act as a complement to other complete systems [8,22,23,24,25]. However, no changes in the anti-phage systems that explain their increased sensitivity to phages were found in the Lcn972R-evolved clones. It is also worth noting that, apart from c2, CHPC1130 and CHPC1183, no other phages in our collection were able to infect L. lactis IPLA517 and IPLA1064. This could be attributed to the combination of the diverse anti-phage systems and their CWPS type, as the C5 subtype is among the less common variants [18].
Genome analyses of the Lcn972R-evolved clones also revealed that, as suspected, some plasmids have been lost, a phenomenon that is frequently observed under stressful conditions when the plasmids in question do not encode essential functions [26]. For example, the plasmid p41-encoding the production of a putative novel EPS had been lost in all Lcn972R clones. EPS biosynthetic clusters, which are present either on plasmids or on the chromosome, are commonly found in Lactococcus [27]. It has been demonstrated that they may confer a selective advantage to the bacteria, providing protection against adverse conditions such as dehydration, antibiotics, toxic compounds or phages [28]. The modular organization of the EPS cluster in p41, with both 5’ and 3’ extremes conserved and all the CDSs in the same orientation except for orfY, is consistent with what has been previously reported [27,29]. Previous works support the idea that some EPSs may also impair phage adsorption by masking the phage receptors [9,28], pointing to the plasmid p41 behind phage resistance in the parental strain.
After the fruitless attempts at transferring the p41 plasmid to demonstrate its role in protecting against phage infection, we used a CRISPR-Cas9 system for lactic acid bacteria to selectively eliminate p41. This approach is known to be effective in the curation of lactococcal plasmids [30,31]. Elimination of p41 was observed after 45 generations at pH 5.5 in accordance with what was previously described when using the same system [16]. Curing of the CRISPR plasmid, on the contrary, was particularly challenging in this study and may be attributed to the plasmid content of this particular strain, which contains four additional plasmids. The fitness cost of maintaining an additional plasmid may not be as significant as in L. cremoris MG1363, a plasmid-free strain [16]. It should be noted that during this process other plasmids were also lost in addition to p41 in some of the tested clones. However, as this was only checked after the cas9 induction, it is not possible to determine whether this loss was due to a putative lack of specificity of the system or to the instability of these plasmids at pH 5.5. The latter seems to be the most likely scenario, as previous studies have demonstrated plasmid instability during prolonged culture of lactococcal cells due to acidification [26,32]. Additionally, although it has not been tested here, it is well documented that bacterial cells with a high plasmid content tend to exhibit significant variability in plasmid numbers and, consequently, in plasmid stability. This variability can result in the emergence of diverse phenotypes when working with single colonies [19].
The selective elimination of p41 allowed us to confirm its role in Ceduovirus resistance because after curing p41 in L. lactis IPLA1064, EOP values resembled those of the Lcn972R clones. Likewise, the observation that phage adsorption to p41-bearing cells is hampered is compatible with the notion that EPS production is blocking access to the phage receptor Pip, at least partially. Consistently, loss of the ability to synthesize the EPS may also explain the increase in surface hydrophobicity previously observed in the Lcn972R derivatives [13]. To the best of our knowledge only the plasmid pCI658 has been shown to be involved in resistance to Ceduovirus [9], while the majority of studies have focused on the role of EPS in resistance to skunaviruses infection [9,33,34]. Prior studies with non-classified phages have also supported the idea of EPSs blocking phage receptors [28,35]. In contrast, EPSs can sometimes be necessary for infection by P335-like phages. This highlights the diverse role of EPSs in phage infection, which is dependent on both the EPSs and the phage involved [11]. Similarly, EPSs may play a dual role either as phage receptor [36] or as a phage-resistance mechanism [37] in Streptococcus thermophilus. In other Gram-positive bacteria, such as Streptococcus pneumoniae [38], Staphylococcus aureus [39], and Streptococcus pyogenes [40], capsular polysaccharides have been shown to prevent phage infection.
The observation that the EPS gene cluster in p41 is also present in L. garvieae is not entirely unexpected. Although L. garvieae is most commonly recognized as a fish pathogen and is rarely associated with human infections, dairy strains have also been described [41,42,43]. This finding reinforces the notion of horizontal gene transfer between both species, particularly when both shared a common niche, as is the case for L. lactis IPLA1064 and L. garvieae M14, which were both isolated from dairy products. While EPS biosynthetic gene clusters are also commonly found in L. garvieae, they have not been as thoroughly characterized as in L. lactis [44,45,46]. Similarly, phage resistance in L. garvieae has not been extensively studied. To our knowledge, only one study has linked the presence of a polysaccharide capsule to a lower infection rate by some phages [47]. Considering the similarity found among certain L. garvieae and other lactococcal phages infecting dairy starters, as well as the documented occurrence of cross-infection between dairy and non-dairy phages [41,48,49], it seems plausible that similar phage-resistance strategies may be employed. It would be of interest to ascertain whether the p41-encoded EPS could confer resistance to L. garvieae phages.

4. Materials and Methods

4.1. Bacterial Strains and Bacteriophages

Bacterial strains and phages used in this work are listed in Table 1. L. lactis parental strains and their Lcn972R derivatives were routinely grown statically in M17 (Formedium, Hunstanton, United Kingdom) supplemented with 0.5% lactose (LM17) at 30 °C. L. cremoris NZ9000 was similarly grown but in M17 (Formedium) supplemented with 0.5% glucose (GM17).
Phages were propagated from previous stocks in GM17 with CaCl2 10 mM (GM17*) using L. cremoris NZ9000 as the host strain. Phage enumeration (pfu ml−1) was performed by the double-layer agar assay in LM17 or GM17 with CaCl2 10 mM and 0.5% glycine (LM17** or GM17**). The top agar (0.4% w/v agar) was inoculated at 4% v/v with an overnight culture of the appropriate host strain and lysate dilutions prepared in SM buffer (20 mM Tris HCl, pH 7.5, 100 mM NaCl, 10 mM Ca(NO3)2, and 10 mM MgSO4). The EOP was determined by dividing the titre (PFU mL−1) of the test strain by that of the L. lactis IPLA1064-E11, that was used as the reference.

4.2. Genome Analysis

A hybrid genome assembly approach was undertaken based on both the available Illumina short reads (Bioproject PRJNA492214) and new long-reads, created using Oxford Nanopore Technology (Oxford, UK). For long-read sequencing, genomic DNA from L. lactis IPLA517 and L. lactis IPLA1064 was extracted by genXone, S.A. (Złotniki, Poland). A total of 200 mg of the pelleted material was homogenized in a BeadBlaster™ 24 (Benchmark Scientific, Sayreville, NJ, USA). Samples were further lysed in buffer G2 from the Genomic DNA Buffer Set (QIAGEN, Hilden, Germany) with Proteinase K and RNase A, and DNA was isolated by using QIAGEN Genomic-tip 20/G kit following the Qiagen Genomic DNA Handbook. Sequencing was conducted using a GridION instrument (ONT, Oxford, UK). Illumina short reads were trimmed using AdapterRemoval v2.2.4 [50] with the following non-default parameters: “--minquality 20 --minlength 30 --trimqualities --trimns --trim5p 15”. Bacteriophage phiX DNA (sequencing control) was removed by only retrieving reads that did not map to the phiX genome, using bwa-mem2 (default parameters). No trimming or adapter removal were performed for the long reads. Hybrid genome assembly was initiated using Flye v2.9 with default parameters. (https://github.com/mikolmogorov/Flye). This was followed by one round of long-read polishing using Medaka v1.5.0 (https://github.com/nanoporetech/medaka) and two rounds of short-read polishing. First Polca [51] and masurca v4.0.8 [52] were used with default parameters, and second, polypolish v0.5.0 [53] and bwa-mem2 [54] with default parameters, in addition to the polypolish_insert_filter parameter.
Variant analysis was conducted with the breseq bacterial variant calling pipeline [55] using the default parameters. For a variant to be called, at least 80% of reads mapping to a position was required to harbor the variant. The PGAP [56] annotated hybrid assembly was applied as a reference while the corresponding short reads of each strain were used as input.
The presence of known phage-resistance systems was analyzed by using the available bioinformatical tools PADLOC and DefenseFinder [57,58]. Additionally, prophages encoded within the genome were identified using PHASTEST [59].

4.3. Plasmid Assemblies and Confirmatory PCRs Reactions

Putative plasmid sequences were manually edited and confirmed by PCR. Primers (Table S1) were designed considering the sequences at the beginning and the end of the putative plasmid contigs and PCR reactions were performed on total DNA in order to close and confirm their sequences. PCR reactions were performed using the 2× Taq master Mix (AMPLIQON, Odense, Denmark) following the manufacturer’s instructions. The amplified DNA fragments were purified with the PCR DNA and Gel Band Purification kit (Cytiva, Buckinghamshire, UK), and their nucleotide sequences determined by Sanger sequencing (Macrogen, Madrid, Spain). The PCR conditions were as follows for all the plasmids, with the exception of p20 where the annealing step was performed at 50 °C: an initial denaturation cycle at 98 °C for 10 s, 35 cycles of a denaturation step at 98 °C for 30 s, an annealing step at 55 °C for 30 s, an extension step at 72 °C for 1 min, and a final extension cycle at 72 °C for 10 min.

4.4. Selective Plasmid Curing

A CRISPR-Cas9 plasmid harboring a specific single-guide RNA (sgRNA) with a spacer for the elimination of the p41 (sgp41) plasmid was constructed. This 20 bp spacer was designed using “CCTop—CRISPR/Cas9 online predictor” (https://cctop.cos.uni-heidelberg.de/) [60] with no possible mismatches to avoid off-target occurrence in the L. lactis IPLA1064 chromosome or plasmids. The 5′ phosphorylated primers P41sgRNAFor and RNDcrePR1Rev (Table S1) were used to amplify and generate the vector with sgp41 sequence by an inverse PCR approach using DNA of pILCsgLC9 as template. [16]. The PCR product corresponding to the linear vector with the p41sgRNA was then subjected to intramolecular ligation, heated to 65 °C to inactivate T4 DNA ligase and digested with DpnI to degrade methylated DNA from the template. Finally, the ligation mixture was used to transform E. coli NZY5α (NZYTech, Lisbon, Portugal) competent cells. The pILCp41LC9 construct obtained in E. coli was then introduced into L. lactis IPLA1064-WT by electro-transformation, as previously described, using chloramphenicol (Cm) 5 µg/mL as the selective marker [61]. The expression of the cas9 gene, located under the control of the acid-inducible weak P170-578 promoter [16], was induced by cultivation in LM17,f adjusted to pH 5.5 by adding HCl 37% (LM17-5.5). To achieve the selective elimination of p41, a single transformant colony was inoculated into liquid LM17-5.5 medium supplemented with Cm and incubated for 45 generations of serial cultivation. Five colonies obtained by plating the previous culture on LM17-5.5 agar plates containing Cm were tested by colony PCR using multiple primer pairs for parallel amplification of a target sequence from each endogenous plasmid. A colony retaining all the plasmids except p41 was selected. Next, pILCp41LC9 was cured by growing the selected clone for 150 generations in liquid LM17 medium in the absence of Cm. Three Cm-sensitive colonies were selected that did not show the pILCp41LC9 (tested with primers_Cas9ForEnd and sgp41R, Supplemental Material) but contained all the endogenous plasmids except p41 (L. lactis IPLA1064-Δp41). Detection of p41 and pILCp41LC9, as well as the rest of the plasmids, was performed by PCR with the appropriate primers (Table S1). In all cases, PCRs were performed using Phusion DNA Polymerase (New England Biolabs, Ipswich, MA, US). The correct nucleotide sequence of all novel constructs was confirmed by automated DNA sequencing.

4.5. Phage Assays

Progression of infection by the phage c2 in broth was followed by measuring the optical density (OD) at 600 nm each 15 min for 24 h at 30 °C in a microtiter plate reader (Tecan Trading AG). Each L. lactis strain was inoculated at 2% v/v, grown until an OD of 0.2 (approx. 5 × 107 cfu mL−1) and subsequently diluted 1:10 or 1:100. The three cultures were mixed with the phage c2 in a 3:1 ratio in a ninety-six well-plate to achieve a multiplicity of infection (MOI) of approximately 10, 1, and 0.1. CaCl2 was added at a final concentration of 10 mM. The area under the curve (AUC) after 24 h of incubation and the percentage of inhibition (% PI) were determined as previously described [62].
Phage adsorption assays were performed as previously described with some modifications [63]. Each L. lactis culture was grown until an OD of 0.5 and mixed with the phage c2 at a MOI of 0.01 in the presence of 10 mM CaCl2. Samples were incubated for 10 min at 30 °C and, after centrifugation, the non-adsorbed phages (residual) were quantified. A sample without cells was equally treated to determine the initial phage titre (control). Experiments were carried out with at least three biological replicates and the percentage of adsorption was calculated as (1−residual/control) × 100.

4.6. Statistical Analyses

The EOP and adsorption data were subjected to statistical analysis using a one-way ANOVA, followed by a Tukey’s multiple comparisons test, employing GraphPad Prism 6 software (Boston, MA, USA). Differences were considered significant when p < 0.05.

5. Conclusions

In this work, we demonstrated that the p41 plasmid containing an EPS biosynthetic gene cluster was responsible for the Ceduovirus resistance in L. lactis IPLA1064. This EPS may block access to the phage receptor, impeding phage adsorption and reducing the efficiency of the infection. Further work could be conducted in two distinct areas in relation to this study. First, the chemical and technological characterizations of the EPS and the EPS-producing L. lactis strains could facilitate their potential application in the dairy industry, e.g., as a texturing agent. Conversely, it would be beneficial to ascertain if it is also implicated in phage resistance in L. garvieae and whether or not the plasmid p41 can be transferred to clinical or pathogenic strains, which would make the potential use of phage therapy in this field more challenging.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms26062508/s1. References [64,65,66] are cited in the Supplementary Material.

Author Contributions

Conceptualization, A.R., G.d.S. and B.M.; Formal analysis, C.R., M.H.R., A.R., G.d.S. and B.M.; Funding acquisition, A.R. and B.M.; Investigation, C.R., P.G. and J.Á.R.-M.; Methodology, J.N.G.-N., J.M. and J.L.B.; Writing—original draft, C.R. and B.M.; Writing—review and editing, C.R., J.N.G.-N., M.H.R., P.G., J.Á.R.-M., J.M., A.R., J.L.B., G.d.S. and B.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by grant no. PID2020-119697RB-I00, funded by MCIN/AEI/10.13039/501100011033, awarded to B.M. J.M. was funded by Research Ireland under grant no. 20/FFP-P/8664. C.R. is a fellow of the program “Ayudas Severo Ochoa” of the Principality of Asturias (BP20 006) and was awarded mobility grants from the Spanish Society for Microbiology (Ayuda César Nombela CN-SEM-2024-005) and the Federation of European Microbiological Societies (FEMS Research and Training, grant no. 1771).

Data Availability Statement

The original contributions presented in this study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author. Manually curated plasmid sequences have been deposited at the NCBI, with the following accession nos: PV037728-IPLA1064_p17; PV037729-IPLA1064_p41; PV037730-IPLA1064_p45; PV037731-IPLA1064_p6; PV037732-IPLA1064_p66; PV037733-IPLA517_p15; and PV037734-IPLA517_p20. DNA sequences have been deposited in the SRA archive under the project PRJNA492214.

Acknowledgments

We thank T. Janzen (Novonesis, Denmark) for sharing phages CHPC1130 and CHPC1183 and K. Al-Nakeeb (Novonesis, Denmark) for his help with the hybrid genome assemblies.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
EPSExopolysaccharide
AbiAbortive infection
R/MRestriction modification
SieSuperinfection exclusion
CWPSsCell wall polysaccharides
PipPhage infection protein
Lcn972RLactococcin 972-resistant derivatives
EOPEfficiency of plaquing
MOIMultiplicity of infection
AUCArea under the curve
PIPercentage of inhibition
ODOptical density

References

  1. Bourdichon, F.; Arias, E.; Babuchowski, A.; Bückle, A.; Dal Bello, F.; Dubois, A.; Fontana, A.; Fritz, D.; Kemperman, R.; Laulund, S.; et al. The Forgotten Role of Food Cultures. FEMS Microbiol. Lett. 2021, 368, fnab085. [Google Scholar] [CrossRef]
  2. Papadimitriou, K.; Alegría, Á.; Bron, P.A.; de Angelis, M.; Gobbetti, M.; Kleerebezem, M.; Lemos, J.A.; Linares, D.M.; Ross, P.; Stanton, C.; et al. Stress Physiology of Lactic Acid Bacteria. Microbiol. Mol. Biol. Rev. 2016, 80, 837–890. [Google Scholar] [CrossRef] [PubMed]
  3. Pujato, S.A.; Quiberoni, A.; Mercanti, D.J. Bacteriophages on Dairy Foods. J. Appl. Microbiol. 2019, 126, 14–30. [Google Scholar] [CrossRef] [PubMed]
  4. Romero, D.A.; Magill, D.; Millen, A.; Horvath, P.; Fremaux, C. Dairy Lactococcal and Streptococcal Phage–Host Interactions: An Industrial Perspective in an Evolving Phage Landscape. FEMS Microbiol. Rev. 2020, 44, 909–932. [Google Scholar] [CrossRef] [PubMed]
  5. Georjon, H.; Bernheim, A. The Highly Diverse Antiphage Defence Systems of Bacteria. Nat. Rev. Microbiol. 2023, 21, 686–700. [Google Scholar] [CrossRef]
  6. Philippe, C.; Cornuault, J.K.; de Melo, A.G.; Morin-Pelchat, R.; Jolicoeur, A.P.; Moineau, S. The never-ending battle between lactic acid bacteria and their phages. FEMS Microbiol. Rev. 2023, 47, fuad035. [Google Scholar] [CrossRef]
  7. Millen, A.M.; Horvath, P.; Boyaval, P.; Romero, D.A. Mobile CRISPR/Cas-Mediated Bacteriophage Resistance in Lactococcus lactis. PLoS ONE 2012, 7, e51663. [Google Scholar] [CrossRef]
  8. Grafakou, A.; Mosterd, C.; Beck, M.H.; Kelleher, P.; McDonnell, B.; de Waal, P.P.; van Rijswijck, I.M.H.; van Peij, N.N.M.E.; Cambillau, C.; Mahony, J.; et al. Discovery of Antiphage Systems in the Lactococcal Plasmidome. Nucleic Acids Res. 2024, 52, 9760–9776. [Google Scholar] [CrossRef]
  9. Forde, A.; Fitzgerald, G.F. Molecular Organization of Exopolysaccharide (EPS) Encoding Genes on the Lactococcal Bacteriophage Adsorption Blocking Plasmid, PCI658. Plasmid 2003, 49, 130–142. [Google Scholar] [CrossRef]
  10. Mahony, J. Cell Surface Polysaccharides Represent a Common Strategy for Adsorption among Phages Infecting Lactic Acid Bacteria: Lessons from Dairy Lactococci and Streptococci. mSystems 2021, 6, 10–1128. [Google Scholar] [CrossRef]
  11. Millen, A.M.; Romero, D.A.; Horvath, P.; Magill, D.; Simdon, L. Host-Encoded, Cell Surface-Associated Exopolysaccharide Required for Adsorption and Infection by Lactococcal P335 Phage Subtypes. Front. Microbiol. 2022, 13, 971166. [Google Scholar] [CrossRef]
  12. Millen, A.M.; Romero, D.A. Genetic Determinants of Lactococcal C2viruses for Host Infection and Their Role in Phage Evolution. J. Gen. Virol. 2016, 97, 1998–2007. [Google Scholar] [CrossRef]
  13. López-González, M.J.; Escobedo, S.; Rodríguez, A.; Rute Neves, A.; Janzen, T.; Martínez, B. Adaptive Evolution of Industrial Lactococcus lactis under Cell Envelope Stress Provides Phenotypic Diversity. Front. Microbiol. 2018, 9, 2654. [Google Scholar] [CrossRef]
  14. Martinez, B.; Suarez, J.E.; Rodriguez, A. Antimicrobials Produced by Wild Lactococcal Strains Isolated from Homemade Cheeses. J. Food Prot. 1995, 58, 1118–1123. [Google Scholar] [CrossRef]
  15. Kuipers, O.P.; De Ruyter, P.G.G.A.; Kleerebezem, M.; De Vos, W.M. Quorum Sensing-Controlled Gene Expression in Lactic Acid Bacteria. J. Biotechnol. 1998, 64, 15–21. [Google Scholar] [CrossRef]
  16. Garay-Novillo, J.N.; Ruiz-Masó, J.A.; del Solar, G.; Barra, J.L. Easy-Curing and PH-Regulated CRISPR-Cas9 Plasmids for Gene Editing and Plasmid Curing in Lactococcus cremoris. Microb. Biotechnol. 2024, 17, e70060. [Google Scholar] [CrossRef] [PubMed]
  17. Kelleher, P.; Mahony, J.; Schweinlin, K.; Neve, H.; Franz, C.M.; van Sinderen, D. Assessing the Functionality and Genetic Diversity of Lactococcal Prophages. Int. J. Food Microbiol. 2018, 272, 29–40. [Google Scholar] [CrossRef]
  18. Mahony, J.; Frantzen, C.; Vinogradov, E.; Sadovskaya, I.; Theodorou, I.; Kelleher, P.; Chapot-Chartier, M.P.; Cambillau, C.; Holo, H.; van Sinderen, D. The CWPS Rubik’s Cube: Linking Diversity of Cell Wall Polysaccharide Structures with the Encoded Biosynthetic Machinery of Selected Lactococcus lactis Strains. Mol. Microbiol. 2020, 114, 582–596. [Google Scholar] [CrossRef]
  19. Siezen, R.J.; Renckens, B.; Van Swam, I.; Peters, S.; Van Kranenburg, R.; Kleerebezem, M.; De Vos, W.M. Complete Sequences of Four Plasmids of Lactococcus lactis subsp. cremoris SK11 Reveal Extensive Adaptation to the Dairy Environment. Appl. Environ. Microbiol. 2005, 71, 8371–8382. [Google Scholar]
  20. Górecki, R.K.; Koryszewska-Bagińska, A.; Gołebiewski, M.; Zylińska, J.; Grynberg, M.; Bardowski, J.K. Adaptative Potential of the Lactococcus lactis IL594 Strain Encoded in Its 7 Plasmids. PLoS ONE 2011, 6, e22238. [Google Scholar] [CrossRef]
  21. Strahinic, I.; Kojic, M.; Tolinacki, M.; Fira, D.; Topisirovic, L. Molecular Characterization of Plasmids PS7a and PS7b from Lactococcus lactis subsp. lactis biovar. diacetylactis S50 as a Base for the Construction of Mobilizable Cloning Vectors. J. Appl. Microbiol. 2009, 106, 78–88. [Google Scholar]
  22. Ainsworth, S.; Mahony, J.; van Sinderen, D. The Plasmid Complement of Lactococcus lactis UC509.9 Encodes Multiple Bacteriophage Resistance Systems. Appl. Environ. Microbiol. 2014, 80, 4341–4349. [Google Scholar] [CrossRef]
  23. Ainsworth, S.; Stockdale, S.; Bottacini, F.; Mahony, J.; van Sinderen, D. The Lactococcus lactis Plasmidome: Much Learnt, yet Still Lots to Discover. FEMS Microbiol. Rev. 2014, 38, 1066–1088. [Google Scholar] [CrossRef]
  24. Cui, Y.; Hu, T.; Qu, X.; Zhang, L.; Ding, Z.; Dong, A. Plasmids from Food Lactic Acid Bacteria: Diversity, Similarity, and New Developments. Int. J. Mol. Sci. 2015, 16, 13172–13202. [Google Scholar] [CrossRef]
  25. Kelleher, P.; Mahony, J.; Bottacini, F.; Lugli, G.A.; Ventura, M.; Van Sinderen, D. The Lactococcus lactis Pan-Plasmidome. Front. Microbiol. 2019, 10, 707. [Google Scholar] [CrossRef]
  26. van Mastrigt, O.; Di Stefano, E.; Hartono, S.; Abee, T.; Smid, E.J. Large Plasmidome of Dairy Lactococcus lactis subsp. lactis biovar diacetylactis FM03P Encodes Technological Functions and Appears Highly Unstable. BMC Genom. 2018, 19, 620. [Google Scholar]
  27. Zeidan, A.A.; Poulsen, V.K.; Janzen, T.; Buldo, P.; Derkx, P.M.F.; Øregaard, G.; Neves, A.R. Polysaccharide Production by Lactic Acid Bacteria: From Genes to Industrial Applications. FEMS Microbiol. Rev. 2017, 41, S168–S200. [Google Scholar] [CrossRef]
  28. Looijesteijn, P.J.; Trapet, L.; De Vries, E.; Abee, T.; Hugenholtz, J. Physiological Function of Exopolysaccharides Produced by Lactococcus lactis. Int. J. Food Microbiol. 2001, 64, 71–80. [Google Scholar] [CrossRef]
  29. Dabour, N.; LaPointe, G. Identification and Molecular Characterization of the Chromosomal Exopolysaccharide Biosynthesis Gene Cluster from Lactococcus lactis subsp. cremoris SMQ-461. Appl. Environ. Microbiol. 2005, 71, 7414–7425. [Google Scholar] [CrossRef]
  30. Berlec, A.; Škrlec, K.; Kocjan, J.; Olenic, M.; Štrukelj, B. Single Plasmid Systems for Inducible Dual Protein Expression and for CRISPR-Cas9/CRISPRi Gene Regulation in Lactic Acid Bacterium Lactococcus lactis. Sci. Rep. 2018, 8, 1009. [Google Scholar] [CrossRef] [PubMed]
  31. van der Els, S.; James, J.K.; Kleerebezem, M.; Bron, P.A. Versatile Cas9-Driven Subpopulation Selection Toolbox for Lactococcus lactis. Appl. Environ. Microbiol. 2018, 84, e02752-17. [Google Scholar] [CrossRef]
  32. Sinha, R.P. A New Simple Method of Curing Plasmids in Lactic Streptococci (Streptococcus cremoris; Streptococcus lactis, Plasmids). FEMS Microbiol. Lett. 1989, 57, 349–352. [Google Scholar] [CrossRef] [PubMed]
  33. Deveau, H.; Van Calsteren, M.R.; Moineau, S. Effect of Exopolysaccharides on Phage-Host Interactions in Lactococcus lactis. Appl. Environ. Microbiol. 2002, 68, 4364–4369. [Google Scholar] [CrossRef]
  34. Millen, A.M.; Magill, D.; Romero, D.; Simdon, L. Evolved Distal Tail Protein of Skunaviruses Facilitates Adsorption to Exopolysaccharide-Encoding Lactococci. Microbiome Res. Rep. 2023, 2, 26. [Google Scholar] [CrossRef]
  35. Akçelik, M.; Şanlibaba, P. Characterisation of an Exopolysaccharide Preventing Phage Adsorption in Lactococcus lactis subsp. cremoris MA39. Turk. J. Vet. Anim. Sci. 2002, 26, 1151–1156. [Google Scholar]
  36. Rodríguez, C.; Van der Meulen, R.; Vaningelgem, F.; Font de Valdez, G.; Raya, R.; De Vuyst, L.; Mozzi, F. Sensitivity of capsular-producing Streptococcus thermophilus strains to bacteriophage adsorption. Lett. Appl. Microbiol. 2008, 46, 462–468. [Google Scholar] [CrossRef]
  37. Janzen, T.; Christiansen, D.E. Lactic Bacterium for Texturizing Food Products Selected on Basis of Phage Resistance. WO2011092300 A1, 4 August 2011. [Google Scholar]
  38. López, R.; García, E. Recent trends on the molecular biology of pneumococcal capsules, lytic enzymes, and bacteriophage. FEMS Microbiol. Rev. 2004, 28, 553–580. [Google Scholar] [PubMed]
  39. Lehmann, E.; van Dalen, R.; Gritsch, L.; Slavetinsky, C.; Korn, N.; Rohmer, C.; Krause, D.; Peschel, A.; Weidenmaier, C.; Wolz, C. The capsular polysaccharide obstructs wall teichoic acid functions in Staphylococcus aureus. J. Infect. Dis. 2024, 230, 1253–1261. [Google Scholar] [PubMed]
  40. Beerens, D.; Franch-Arroyo, S.; Sullivan, T.J.; Goosmann, C.; Brinkmann, V.; Charpentier, E. Survival strategies of Streptococcus pyogenes in response to phage infection. Viruses 2021, 13, 612. [Google Scholar] [CrossRef]
  41. Flórez, A.B.; Mayo, B. The Plasmid Complement of the Cheese Isolate Lactococcus garvieae IPLA 31405 Revealed Adaptation to the Dairy Environment. PLoS ONE 2015, 10, e0126101. [Google Scholar]
  42. Moumene, M.; Drissi, F.; Croce, O.; Djebbari, B.; Robert, C.; Angelakis, E.; Benouareth, D.E.; Raoult, D.; Merhej, V. Complete Genome Sequence and Description of Lactococcus garvieae M14 Isolated from Algerian Fermented Milk. New Microbes New Infect. 2016, 10, 122–131. [Google Scholar] [CrossRef]
  43. Ferrario, C.; Ricci, G.; Milani, C.; Lugli, G.A.; Ventura, M.; Eraclio, G.; Borgo, F.; Fortina, M.G. Lactococcus garvieae: Where Is It from? A First Approach to Explore the Evolutionary History of This Emerging Pathogen. PLoS ONE 2013, 8, e84796. [Google Scholar] [CrossRef]
  44. Ayyash, M.; Abu-Jdayil, B.; Itsaranuwat, P.; Almazrouei, N.; Galiwango, E.; Esposito, G.; Hunashal, Y.; Hamed, F.; Najjar, Z. Exopolysaccharide Produced by the Potential Probiotic Lactococcus garvieae C47: Structural Characteristics, Rheological Properties, Bioactivities and Impact on Fermented Camel Milk. Food Chem. 2020, 333, 127418. [Google Scholar] [CrossRef]
  45. Morita, H.; Toh, H.; Oshima, K.; Yoshizaki, M.; Kawanishi, M.; Nakaya, K.; Suzuki, T.; Miyauchi, E.; Ishii, Y.; Tanabe, S.; et al. Complete Genome Sequence and Comparative Analysis of the Fish Pathogen Lactococcus garvieae. PLoS ONE 2011, 6, e23184. [Google Scholar] [CrossRef]
  46. Niederle, M.V.; Bosch, J.; Ale, C.E.; Nader-Macías, M.E.; Aristimuño Ficoseco, C.; Toledo, L.F.; Valenzuela-Sánchez, A.; Soto-Azat, C.; Pasteris, S.E. Skin-Associated Lactic Acid Bacteria from North American Bullfrogs as Potential Control Agents of Batrachochytrium dendrobatidis. PLoS ONE 2019, 14, e0223020. [Google Scholar] [CrossRef]
  47. Hoai, T.D.; Nishiki, I.; Yoshida, T.; Nakai, T. Host Range and Influence of a Cell Capsule on the Phage Efficacy of Three Lactococcus garvieae Lytic Phages. Dis. Aquat. Organ. 2018, 128, 81–86. [Google Scholar] [CrossRef]
  48. Eraclio, G.; Fortina, M.G.; Labrie, S.J.; Tremblay, D.M.; Moineau, S. Characterization of Prophages of Lactococcus garvieae. Sci. Rep. 2017, 7, 1856. [Google Scholar] [CrossRef]
  49. Eraclio, G.; Tremblay, D.M.; Lacelle-Côté, A.; Labrie, S.J.; Fortina, M.G.; Moineau, S. A Virulent Phage Infecting Lactococcus garvieae, with Homology to Lactococcus lactis Phages. Appl. Environ. Microbiol. 2015, 81, 8358–8365. [Google Scholar] [CrossRef]
  50. Schubert, M.; Lindgreen, S.; Orlando, L. AdapterRemoval v2: Rapid Adapter Trimming, Identification, and Read Merging. BMC Res. Notes 2016, 9, 88. [Google Scholar] [CrossRef] [PubMed]
  51. Zimin, A.V.; Salzberg, S.L. The Genome Polishing Tool POLCA Makes Fast and Accurate Corrections in Genome Assemblies. PLoS Comput. Biol. 2020, 16, e1007981. [Google Scholar] [CrossRef] [PubMed]
  52. Zimin, A.V.; Marçais, G.; Puiu, D.; Roberts, M.; Salzberg, S.L.; Yorke, J.A. The MaSuRCA Genome Assembler. Bioinformatics 2013, 29, 2669–2677. [Google Scholar] [CrossRef] [PubMed]
  53. Wick, R.R.; Holt, K.E. Polypolish: Short-Read Polishing of Long-Read Bacterial Genome Assemblies. PLoS Comput. Biol. 2022, 18, e1009802. [Google Scholar] [CrossRef]
  54. Md, V.; Misra, S.; Li, H.; Aluru, S. Efficient Architecture-Aware Acceleration of BWA-MEM for Multicore Systems. In Proceedings of the 2019 IEEE International Parallel and Distributed Processing Symposium (IPDPS), Rio de Janeiro, Brazil, 20–24 May 2019; pp. 314–324. [Google Scholar]
  55. Deatherage, D.E.; Barrick, J.E. Identification of Mutations in Laboratory-Evolved Microbes from Next-Generation Sequencing Data Using Breseq. In Engineering and Analyzing Multicellular Systems: Methods and Protocols; Humana Press: New York, NY, USA, 2014; pp. 165–188. [Google Scholar]
  56. Tatusova, T.; DiCuccio, M.; Badretdin, A.; Chetvernin, V.; Ciufo, S.; Li, W. Prokaryotic Genome Annotation Pipeline; National Center for Biotechnology Information: Bethesda, MD, USA, 2013. [Google Scholar]
  57. Tesson, F.; Hervé, A.; Mordret, E.; Touchon, M.; d’Humières, C.; Cury, J.; Bernheim, A. Systematic and Quantitative View of the Antiviral Arsenal of Prokaryotes. Nat. Commun. 2022, 13, 2561. [Google Scholar] [CrossRef]
  58. Payne, L.J.; Meaden, S.; Mestre, M.R.; Palmer, C.; Toro, N.; Fineran, P.C.; Jackson, S.A. PADLOC: A Web Server for the Identification of Antiviral Defence Systems in Microbial Genomes. Nucleic Acids Res. 2022, 50, W541–W550. [Google Scholar] [CrossRef]
  59. Wishart, D.S.; Han, S.; Saha, S.; Oler, E.; Peters, H.; Grant, J.R.; Stothard, P.; Gautam, V. PHASTEST: Faster than PHASTER, Better than PHAST. Nucleic Acids Res. 2023, 51, W443–W450. [Google Scholar] [CrossRef]
  60. Stemmer, M.; Thumberger, T.; Del Sol Keyer, M.; Wittbrodt, J.; Mateo, J.L. CCTop: An Intuitive, Flexible and Reliable CRISPR/Cas9 Target Prediction Tool. PLoS ONE 2015, 10, e0124633. [Google Scholar] [CrossRef]
  61. Holo, H.; Nes, I.F. High-Frequency Transformation, by Electroporation, of Lactococcus lactis subsp. cremoris Grown with Glycine in Osmotically Stabilized Media. Appl. Environ. Microbiol. 1989, 55, 3119–3123. [Google Scholar]
  62. Xie, Y.; Wahab, L.; Gill, J.J. Development and Validation of a Microtiter Plate-Based Assay for Determination of Bacteriophage Host Range and Virulence. Viruses 2018, 10, 189. [Google Scholar] [CrossRef]
  63. McDonnell, B.; Mahony, J.; Neve, H.; Hanemaaijer, L.; Noben, J.P.; Kouwen, T.; van Sinderen, D. Identification and Analysis of a Novel Group of Bacteriophages Infecting the Lactic Acid Bacterium Streptococcus thermophilus. Appl. Environ. Microbiol. 2016, 82, 5153–5165. [Google Scholar] [CrossRef]
  64. Tamura, K.; Stecher, G.; Kumar, S. MEGA 11: Molecular Evolutionary Genetics Analysis Version 11. Mol. Biol. Evol. 2021, 38, 3022–3027. [Google Scholar] [CrossRef] [PubMed]
  65. Felsenstein, J. Confidence limits on phylogenies: An approach using the bootstrap. Evolution 1985, 39, 783–791. [Google Scholar] [CrossRef] [PubMed]
  66. Zuckerkandl, E.; Pauling, L. Evolutionary divergence and convergence in proteins. In Evolving Genes and Proteins; Bryson, V., Vogel, H.J., Eds.; Academic Press: New York, NY, USA, 1965; pp. 97–166. [Google Scholar]
Figure 1. EPS biosynthetic cluster carried by plasmid p41. (A) p41 Plasmid map. Color code for genes (arrows): mobilization (light blue); replication (red); insertion sequences/recombinases (black); hypothetical/unknown function (white); miscellaneous (grey); EPS assembly (green); EPS modulation (yellow); EPS glycosyl transferases/sugar modification (orange); and nucleotide sugar biosynthesis (pink). The figure was created with SnapGene (https://www.snapgene.com). (B) Sequence relatedness of lactococcal EPS biosynthetic gene clusters. Strain codes or plasmid names are shown to the left (GenBank accession number), from top to bottom: L. cremorisi, B40; L. lactis, KLDS 4.0325; L. petauri, R21-74-AL-1; L. lactis, IPLA1064; L. garviae, M14; L. cremoris, HO2; L. cremoris, DGCC6073. The EPS cluster from KLDS 4.0325 is chromosome-borne. Figure was created with Clinker (https://cagecat.bioinformatics.nl/tools/clinker).
Figure 1. EPS biosynthetic cluster carried by plasmid p41. (A) p41 Plasmid map. Color code for genes (arrows): mobilization (light blue); replication (red); insertion sequences/recombinases (black); hypothetical/unknown function (white); miscellaneous (grey); EPS assembly (green); EPS modulation (yellow); EPS glycosyl transferases/sugar modification (orange); and nucleotide sugar biosynthesis (pink). The figure was created with SnapGene (https://www.snapgene.com). (B) Sequence relatedness of lactococcal EPS biosynthetic gene clusters. Strain codes or plasmid names are shown to the left (GenBank accession number), from top to bottom: L. cremorisi, B40; L. lactis, KLDS 4.0325; L. petauri, R21-74-AL-1; L. lactis, IPLA1064; L. garviae, M14; L. cremoris, HO2; L. cremoris, DGCC6073. The EPS cluster from KLDS 4.0325 is chromosome-borne. Figure was created with Clinker (https://cagecat.bioinformatics.nl/tools/clinker).
Ijms 26 02508 g001
Figure 2. Plasmid-specific PCR assays conducted during the curing of p41. L. lactis IPLA1064 (1), L. lactis IPLA1064-∆p41/pILCp41LC9 (2), and L. lactis IPLA1064-∆p41 (3). Lane A: amplicon of pILCp41LC9. Lane B (from top to bottom) amplicons of plasmids p17, p41, and p6. Lane C: amplicons of plasmids p66 (top) and p45 (bottom). FastGene 1 kb DNA Marker Plus (Nippon Genetics Europe) was used as the DNA ladder.
Figure 2. Plasmid-specific PCR assays conducted during the curing of p41. L. lactis IPLA1064 (1), L. lactis IPLA1064-∆p41/pILCp41LC9 (2), and L. lactis IPLA1064-∆p41 (3). Lane A: amplicon of pILCp41LC9. Lane B (from top to bottom) amplicons of plasmids p17, p41, and p6. Lane C: amplicons of plasmids p66 (top) and p45 (bottom). FastGene 1 kb DNA Marker Plus (Nippon Genetics Europe) was used as the DNA ladder.
Ijms 26 02508 g002
Figure 3. Infection of L. lactis IPLA1064 and derivatives by phages c2, CHPC1130, and CHPC1183. (A) Efficiency of plaquing (EOP) with L. lactis IPLA1064-E11 as a reference. The results are the mean and standard deviation of at least three independent experiments. EOPs on L. lactis IPLA1064 are close to 0 (0.0003 ± 0.0004 for c2; 0.001 ± 0.0004 for CHPC1130; 0.005 ± 0.002 for CHPC1183). The results with different letters (a, b, c) represent significant differences (p < 0.05). (B) Double-layer agar spot test (phage lysate tenfold dilutions) showing CHPC1130 lysis plaques on L. lactis IPLA1064 (A), L. lactis IPLA1064-E11 (B), and two clones of L. lactis IPLA1064-∆p41 (C,D). Due to their small size, lysis plaques on L. lactis IPLA1064 are hardly visible (arrow).
Figure 3. Infection of L. lactis IPLA1064 and derivatives by phages c2, CHPC1130, and CHPC1183. (A) Efficiency of plaquing (EOP) with L. lactis IPLA1064-E11 as a reference. The results are the mean and standard deviation of at least three independent experiments. EOPs on L. lactis IPLA1064 are close to 0 (0.0003 ± 0.0004 for c2; 0.001 ± 0.0004 for CHPC1130; 0.005 ± 0.002 for CHPC1183). The results with different letters (a, b, c) represent significant differences (p < 0.05). (B) Double-layer agar spot test (phage lysate tenfold dilutions) showing CHPC1130 lysis plaques on L. lactis IPLA1064 (A), L. lactis IPLA1064-E11 (B), and two clones of L. lactis IPLA1064-∆p41 (C,D). Due to their small size, lysis plaques on L. lactis IPLA1064 are hardly visible (arrow).
Ijms 26 02508 g003
Figure 4. Phage c2 adsorption (%) to Lactococcus cells. The results are the mean and standard deviation of, at least, three independent experiments. Significant differences (p < 0.05) are represented by letters (a, b).
Figure 4. Phage c2 adsorption (%) to Lactococcus cells. The results are the mean and standard deviation of, at least, three independent experiments. Significant differences (p < 0.05) are represented by letters (a, b).
Ijms 26 02508 g004
Table 1. Bacterial strains, phages, and plasmids used in this study.
Table 1. Bacterial strains, phages, and plasmids used in this study.
StrainDescriptionReference
Lactococcus cremoris
NZ9000L. cremoris MG1363 pepN::nisRK. Host for phage propagation[15]
Lactococcus lactis
IPLA517Phage-resistant; 6 endogenous plasmids[14]
IPLA517-B5IPLA517 Lcn972R-evolved clone; phage-sensitive; 4 endogenous plasmids[13]
IPLA517-C6IPLA517 Lcn972R-evolved clone; phage-sensitive; 4 endogenous plasmids[13]
IPLA1064Phage-resistant; 5 endogenous plasmids[14]
IPLA1064-C11IPLA1064 Lcn972R-evolved clone; phage-sensitive; 4 endogenous plasmids[13]
IPLA1064-E11IPLA1064 Lcn972R-evolved clone; phage-sensitive; 4 endogenous plasmids[13]
IPLA1064-Δp41IPLA1064 cured of plasmid p41This work
Bacteriophages
c2Ceduovirus (GenBank NC_001706)Lab collection
CHPC1130CeduovirusNovonesis collection
CHPC1183Ceduovirus (GenBank MN689511)Novonesis collection
Plasmids
pILCsgLC9Derived from pIL253 with a CRISPR-Cas9 cassette for specific spacers cloning[16]
pILCp41LC9Derived from pILCsgLC9 with a specific spacer for p41 curationThis work
Table 2. Size and main features of endogenous plasmids found in L. lactis IPLA517 and L. lactis IPLA1064.
Table 2. Size and main features of endogenous plasmids found in L. lactis IPLA517 and L. lactis IPLA1064.
Plasmid (Accession Number)Size (bp)Present inIdentity IPLA1064 vs. IPLA517Main Features
p6
(PV037731)
5860Both100%RepB; R/M subunit
p15
(PV037733)
15315IPLA517NA 1RepB and RepX; MobC;
Pyrimidine synthesis genes (tmk, pyrF, cdd)
p17
(PV037728)
16890IPLA106499% to IPLA 517_p15 plus a 1678-bp insertionRepB and RepX; MobC;
Pyrimidine synthesis genes (tmk, pyrF, cdd)
p20
(PV037734)
20376IPLA517NA 1RepB
Pili genes; cadA
p41
(PV037729)
41314Both100%RepB and RepA; MobC;
EPS cluster
p45
(PV037730)
45496Both99%RepB; MobC;
Type I R/M system; PrtP and Clp (protease)
p66
(PV037732)
66477Both100%RepB; MobC
lac operon; opp operon; pepF and pepO
1 NA: not applicable.
Table 3. Phage-resistance systems present in both L. lactis IPLA517 and L. lactis IPLA1064 as identified through DefenseFinder and PADLOC.
Table 3. Phage-resistance systems present in both L. lactis IPLA517 and L. lactis IPLA1064 as identified through DefenseFinder and PADLOC.
Resistance SystemLocation (L. lactis IPLA1064 Coordinates)Hit GeneBlastp
(% Identity/Query Cover)
Identified by
Restriction-modification system type IVChromosome
(87858–89420)
RM—RM_Type_IVWP_081213772. McrB family protein (99.81%/100%)DefenseFinder
AbiHChromosome
(788281–789141)
AbiH—AbiHWP_153004554. AbiH family protein (99.65%/100%)DefenseFinder
Restriction-modification system type IPlasmid (p45)
(7069–12972)
RM__Type_I_REases; RM__Type_I_MTases; RM__Type_I_SWP_124137260.1. Type I restriction endonuclease subunit R (100%/100%);
WP_124137261.1. Type I restriction–modification system subunit M (100%/100%);
WP_124137264.1. Restriction endonuclease subunit S (100%/100%)
PADLOC and DefenseFinder
Table 4. Prophage content of both L. lactis IPLA517 and L. lactis IPLA1064 based on PHASTEST output.
Table 4. Prophage content of both L. lactis IPLA517 and L. lactis IPLA1064 based on PHASTEST output.
ProphageRegion LengthPosition in L. lactis IPLA1064Closest Phage (%Identity/Query Cover)Observations
#140 kb178022–218024ul36 (NC_004066)
(93.27%/50%)
Intact
#220.4 kb464894–485327bIL3140 (NC_002669)
(86.74%/2%)
Incomplete
#346.1 kb724573–770760PLgT-1 (NC_031016)
(82.58%/33%)
Intact
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Rendueles, C.; Garay-Novillo, J.N.; Rau, M.H.; Gaspar, P.; Ruiz-Masó, J.Á.; Mahony, J.; Rodríguez, A.; Barra, J.L.; del Solar, G.; Martínez, B. A Plasmid-Encoded Surface Polysaccharide Partly Blocks Ceduovirus Infection in Lactococci. Int. J. Mol. Sci. 2025, 26, 2508. https://doi.org/10.3390/ijms26062508

AMA Style

Rendueles C, Garay-Novillo JN, Rau MH, Gaspar P, Ruiz-Masó JÁ, Mahony J, Rodríguez A, Barra JL, del Solar G, Martínez B. A Plasmid-Encoded Surface Polysaccharide Partly Blocks Ceduovirus Infection in Lactococci. International Journal of Molecular Sciences. 2025; 26(6):2508. https://doi.org/10.3390/ijms26062508

Chicago/Turabian Style

Rendueles, Claudia, Javier Nicolás Garay-Novillo, Martin Holm Rau, Paula Gaspar, José Ángel Ruiz-Masó, Jennifer Mahony, Ana Rodríguez, José Luis Barra, Gloria del Solar, and Beatriz Martínez. 2025. "A Plasmid-Encoded Surface Polysaccharide Partly Blocks Ceduovirus Infection in Lactococci" International Journal of Molecular Sciences 26, no. 6: 2508. https://doi.org/10.3390/ijms26062508

APA Style

Rendueles, C., Garay-Novillo, J. N., Rau, M. H., Gaspar, P., Ruiz-Masó, J. Á., Mahony, J., Rodríguez, A., Barra, J. L., del Solar, G., & Martínez, B. (2025). A Plasmid-Encoded Surface Polysaccharide Partly Blocks Ceduovirus Infection in Lactococci. International Journal of Molecular Sciences, 26(6), 2508. https://doi.org/10.3390/ijms26062508

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