
==== Front
J Bacteriol
J Bacteriol
jb
Journal of Bacteriology
0021-9193
1098-5530
American Society for Microbiology 1752 N St., N.W., Washington, DC

39133004
jb00145-24
10.1128/jb.00145-24
jb.00145-24
Research Article
editors-pickEditor’s PickbacteriologyBacteriologyVibrio cholerae pathogenicity island 2 encodes two distinct types of restriction systems
https://orcid.org/0000-0002-2466-9679
Vizzarro Grazia 1
https://orcid.org/0000-0002-7997-5526
Lemopoulos Alexandre 1
https://orcid.org/0000-0002-9644-5607
Adams David William 1 david.adams@epfl.ch

https://orcid.org/0000-0002-7024-1489
Blokesch Melanie 1 melanie.blokesch@epfl.ch

1 Laboratory of Molecular Microbiology, Global Health Institute, School of Life Sciences, Ecole Polytechnique Fédérale de Lausanne (EPFL) , Lausanne, Switzerland
Editor Bondy-Denomy Joseph University of California San Francisco , San Francisco, California, USA

Address correspondence to David William Adams, david.adams@epfl.ch
Address correspondence to Melanie Blokesch, melanie.blokesch@epfl.ch
The authors declare no conflict of interest.

9 2024
12 8 2024
12 8 2024
206 9 e00145-2405 4 2024
15 7 2024
Copyright © 2024 Vizzarro et al.
2024
Vizzarro et al.
https://creativecommons.org/licenses/by/4.0/ This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

ABSTRACT

In response to predation by bacteriophages and invasion by other mobile genetic elements such as plasmids, bacteria have evolved specialized defense systems that are often clustered together on genomic islands. The O1 El Tor strains of Vibrio cholerae responsible for the ongoing seventh cholera pandemic (7PET) contain a characteristic set of genomic islands involved in host colonization and disease, many of which contain defense systems. Notably, Vibrio pathogenicity island 2 contains several characterized defense systems as well as a putative type I restriction-modification (T1RM) system, which, interestingly, is interrupted by two genes of unknown function. Here, we demonstrate that the T1RM system is active, methylates the host genomes of a representative set of 7PET strains, and identify a specific recognition sequence that targets non-methylated plasmids for restriction. We go on to show that the two genes embedded within the T1RM system encode a novel two-protein modification-dependent restriction system related to the GmrSD family of type IV restriction enzymes. Indeed, we show that this system has potent anti-phage activity against diverse members of the Tevenvirinae, a subfamily of bacteriophages with hypermodified genomes. Taken together, these results expand our understanding of how this highly conserved genomic island contributes to the defense of pandemic V. cholerae against foreign DNA.

IMPORTANCE

Defense systems are immunity systems that allow bacteria to counter the threat posed by bacteriophages and other mobile genetic elements. Although these systems are numerous and highly diverse, the most common types are restriction enzymes that can specifically recognize and degrade non-self DNA. Here, we show that the Vibrio pathogenicity island 2, present in the pathogen Vibrio cholerae, encodes two types of restriction systems that use distinct mechanisms to sense non-self DNA. The first system is a classical Type I restriction-modification system, and the second is a novel modification-dependent type IV restriction system that recognizes hypermodified cytosines. Interestingly, these systems are embedded within each other, suggesting that they are complementary to each other by targeting both modified and non-modified phages.

KEYWORDS

V. cholerae
pathogenicity island
restriction systems
DNA modification
Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung (SNF) 310030_185022 Blokesch Melanie Howard Hughes Medical Institute (HHMI) 55008726 Blokesch Melanie cover-dateSeptember 2024
==== Body
pmcINTRODUCTION

Mobile genetic elements (MGEs) such as plasmids, transposons, and integrative-conjugative elements can confer significant fitness advantages by facilitating the transfer of beneficial traits to the host bacterium, including key virulence factors or antibiotic-resistant genes (1). However, their maintenance and or activity can also impose a metabolic burden on the host cell, while elements that integrate on the chromosome have the potential to disrupt important genomic features (2). Furthermore, the replication of some MGEs such as lytic bacteriophages (phages) can result in the death of the host cell (3). Indeed, predation by phages, which are ubiquitous bacterial viruses, has imposed a strong evolutionary pressure to develop multiple lines of defense against these MGEs, including a vast array of specialized defense systems (4, 5).

Upon recognizing an infection, these systems can either respond directly by degrading the invading non-self DNA and thus provide individual level protection or alternatively can sacrifice the host cell prior to phage-induced lysis to protect the surrounding population (abortive infection) (6). The most common and best-studied defense systems are restriction-modification (RM) systems, which use restriction enzymes to directly degrade non-self DNA (5). Types I–III RM systems are modification-blocked enzymes that recognize specific DNA sequences and only cut DNA when it is unmodified, while the corresponding sequences in the host genome are protected by epigenetic modification with a cognate methylase (7–9). In contrast, type IV systems are modification-dependent enzymes that can recognize and degrade invading DNA with specific modifications, which are used by some phages to avoid restriction by modification-blocked systems (9, 10).

Diverse defense systems, including RM systems, tend to cluster together within genomic islands known as “defense islands” (11–13). This pattern also applies to the defense systems identified so far in Vibrio cholerae, the causative agent of cholera. This bacterium features specialized islands crucial to its pathogenic evolution. Indeed, only certain V. cholerae strains, referred to as toxigenic isolates, can cause cholera. This ability is due to the presence of two key virulence/colonization factors: the cholera toxin (CT) and toxin-coregulated pilus, encoded on the CTXΦ prophage and the Vibrio pathogenicity island 1, respectively (14–16). The ongoing seventh cholera pandemic is caused by the O1 El Tor V. cholerae lineage (7PET), which uniquely carries the Vibrio seventh pandemic islands I and II (VSP-I and VSP-II), characteristic of the 7PET strains (17, 18). These genomic islands are implicated in defense, as they encode for instance CBASS and AvcD systems (VSP-I) and the Lamassu system DdmABC on VSP-II (19–23). Additionally, toxigenic V. cholerae strains carry the Vibrio pathogenicity island 2 (VPI-2), which is believed to enhance pathogenicity by giving the pathogen a competitive advantage in using sialic acid as a carbon source during gut colonization (24–26). This capability is encoded within the island’s nan-nag genomic region (24–26). Moreover, the island houses several genes believed to protect against MGEs, including (i) a predicted Zorya type I system, a phage defense system identified across a wide range of bacterial genomes and experimentally studied primarily through Escherichia coli homologs (12, 27); (ii) the DNA defense module DdmDE that targets and degrades small multicopy plasmids (23); and (iii) a gene cluster/operon predicted to encode a type I restriction-modification (T1RM) system (24). The presence of both predicted and established defense systems encoded within VPI-2 suggests that it may serve as a genuine defense island.

In this study, we set out to characterize the predicted T1RM operon within VPI-2. We show that the T1RM system promotes methylation of the genomes of 7PET V. cholerae strains, and identify a specific recognition sequence that can target non-self-derived plasmids for restriction. Furthermore, we discovered two genes embedded within the T1RM operon that form a novel modification-dependent restriction system related to the GmrSD family of type IV restriction enzymes, which we term TgvAB. When produced in E. coli, this system has potent anti-phage activity against phages with hypermodified genomes. Collectively, these findings enhance our understanding of how this highly conserved genomic island contributes to the defense of pandemic V. cholerae against foreign DNA.

RESULTS AND DISCUSSION

In silico analysis of VPI-2 and the T1RM cluster

Although VPI-2 was discovered over 20 years ago (24), the genes it carries have not yet been fully characterized. To begin bridging this knowledge gap, we started by re-evaluating the gene content of the island. Consistent with earlier findings (24), this revealed that VPI-2 is highly conserved among a set of 7PET O1 strains isolated between 1975 and 2011. Our analysis confirmed the island’s modular structure as described by Jermyn and Boyd (24, 28), including a predicted Zorya system (12) encoded by genes VC1761–64 [as per reference strain N16961; (29)], a predicted T1RM system (VC1765–69) (24), the DdmDE defense module (VC1770–71) (23), the nan-nag sialic acid utilization cluster (VC1773–1784), and a region with phage-like properties (VC1791–1809) (24) (Fig. 1a). Notably, O139 serogroup strains such as MO10 carry a highly truncated version of VPI-2 that retains only the phage-like region (24, 30) (Fig. 1a).

Fig 1 VPI-2 exhibits high conservation across 7PET V. cholerae strains. (a) Comparative genome alignment of the VPI-2 across a selection of 7PET O1 and O139 strains, isolated between 1975 and 2011. The genomes of these strains are displayed alongside their designated strain names (left) and their isolation dates and locations (right). Coding sequences within the genomes are represented by arrows, with gray bars connecting them to indicate amino acid identity percentages at or above a threshold of 0.93. Instances of lower identity are highlighted in black boxes. Gene locus tags are derived from the reference genome of strain N16961. Predicted or established functions are labeled above each cluster. (b) Close-up examination of the VC1765–69 gene cluster in strains A1552 and DRC193A reveals three genes responsible for the components of the putative T1RM system (hsdR, hsdS, hsdM). A comparative alignment highlights the disruption of the T1RM cluster in strain DRC193A, caused by an IS256 transposon insertion within hsdS.

Given the observed conservation of VPI-2 and the presence of established defense systems, we explored the possibility that the putative T1RM system was also actively involved in restricting foreign DNA. Interestingly, the previously annotated T1RM region sits within a five-gene operon, of which three genes encode homologs of the known T1RM components (Fig. 1b). These host-specificity determinant (hsd) genes encode the specificity subunit HsdS, which recognizes a specific DNA recognition sequence; the methylase subunit HsdM, which methylates (and therefore protects) the recognition sequences in the host genome; and HsdR, the restriction enzyme subunit, which upon encountering foreign DNA with an unmethylated recognition sequence translocates the flanking DNA and cleaves at variable distances from the recognition site (7, 31–33). These components function together as multi-subunit complexes capable of both methylating and restricting DNA. Importantly, restriction requires a pentameric complex of 2HsdR + 2HsdM + HsdS, and although HsdR is dispensable for methylation, HsdS is required for both activities (7, 31). Interestingly, two genes of unknown function are embedded within the T1RM cluster (VC1767–66; Fig. 1b), which we characterize in the subsequent sections below.

Deciphering the recognition motif of VPI-2’s T1RM system

If the T1RM system is active in V. cholerae, then we predicted that we should be able to detect a specific methylation signature that is absent in strains lacking this system. To test this hypothesis, we used SMRT PacBio whole-genome sequencing, which can detect the presence of various DNA modifications including methylation, to determine the methylomes of a selection of 7PET O1 serogroup strains (strains as in Fig. 1a), as well as those of control strains lacking the T1RM system (see Materials and Methods) (34–36). As shown in Fig. 2a, this analysis revealed a unique 13-nucleotide motif with methylation marks located on the second nucleotide within the sequence GATGNNNNNNCTT (m6A: GATGNNNNNNCTT:2). Upon further examination, we discovered that this DNA motif is present in over 600 copies throughout the genome of each strain and is modified in nearly 100% of cases in all O1 serogroup strains, except DRC193A (Fig. 2a). This phenotype is likely explained by the interruption of hsdS in this strain by an IS256-like transposase gene (37) (Fig. 1b). Finally, and as expected, the O139 serogroup strain MO10, which is missing the T1RM-encoding region of VPI-2 (Fig. 1a), and both a VPI-2 and a VC1765–69-deficient deletion strain (Table S1) all lacked this particular methylation mark (Fig. 2a).

Fig 2 Type I RM system’s role in chromosomal methylation and plasmid restriction. (a) SMRT sequencing uncovers a distinctive modified DNA motif across various 7PET V. cholerae strains. The gray bars show the number of the DNA motif (GATGNNNNNNCTT) in each genome, while the blue bars denote the percentage of this motif methylated in each strain’s genome (see secondary y-axis on the right). (b) Diagrams of plasmid constructs. Two plasmids were engineered to either contain (Pmotif+) or lack (Pmotif-) the T1RM recognition motif. The latter plasmid was generated through the introduction of silent mutations. (c, d) The T1RM system hinders plasmid acquisition. Transformation assays compare the uptake of the two plasmids (Pmotif+ and Pmotif-) in strains A1552, A1552∆ddmABC, and A1552∆ddmABC∆VC1765–69 (c) or A1552∆ddmABC and A1552∆ddmABC∆hsdR (d). Statistical differences were calculated on log-transformed data using a two-way analysis of variance corrected for multiple comparisons with Šidák’s method. *P < 0.05; **P < 0.01; ***P < 0.001; ns, not significant. <d.l., below detection limit.

The T1RM impairs plasmid acquisition

Having identified the methylated recognition motif, we next tested the ability of this motif to target plasmids for restriction by the VPI-2 T1RM system. Serendipitously, we realized that the recognition sequence is present within the widely used gentamicin resistance cassette aacC1. We therefore created plasmid derivatives carrying aacC1 either with the putative recognition sequence intact (Pmotif+) or with silent mutations that disrupt the nucleotide recognition sequence while preserving the protein coding sequence (Pmotif-) (Fig. 2b). We then purified these plasmids from E. coli and used them as substrates in an electroporation-based transformation assay to compare their transformation frequencies in various backgrounds. As shown in Fig. 2c, transformation with plasmid Pmotif+ was below the detection limit in the wild-type (WT) background (strain A1552) even though transformants could readily be obtained with plasmid Pmotif-. Furthermore, this disparity between the acquisition of the two plasmids became even stronger in the absence of the DdmABC system (23) (Fig. 2c), which is known to target derivatives of this high-copy number plasmid (38).

To determine if the plasmid restriction was mediated by the T1RM system, we removed either the entire five-gene restriction cluster or only the hsdR gene from the A1552∆ddmABC background and then assessed the plasmid transformability of the resulting strains. As shown in Fig. 2c and d, both deletions led to the recovery of Pmotif+ transformants. Moreover, the transformation difference between the Pmotif+ and Pmotif- plasmids was now no longer statistically significant. Consequently, we conclude that the T1RM system is active, that it methylates a specific recognition sequence, and that when this sequence is present on non-self DNA, the acquisition of this non-methylated DNA is restricted in an HsdR-dependent manner.

Genes embedded in the T1RM cluster protect against phages with modified genomes

Type I restriction-modification systems are recognized for their important role in defending the cell against phage infection (39). Therefore, we aimed to investigate the ability of the entire RM cluster, including the two embedded genes, to protect against viral infections. However, given that commonly used Vibrio phages, such as ICP1, ICP2, and ICP3, are typically isolated using VPI-2-carrying 7PET strains as the host [for example, strain E7946 and its derivatives (40)] and that they coevolved with these strains in cholera endemic areas (41, 42), it is unlikely that any defense system encoded on VPI-2 would provide protection against these phages, potentially due to the presence of phage-encoded anti-defense mechanism. Therefore, we engineered the E. coli strain MG1655 to carry an arabinose-inducible version of the entire five-gene RM cluster (VC1769–65), which was integrated into its chromosome. Utilizing this strain and a strain without the cluster as a control, we screened for protection against the BASEL collection, a recently established phage collection that represents the natural diversity of E. coli phages (43). As shown in Fig. 3, we noted a reduction in the efficiency of plaquing of at least 1,000-fold compared to the non-defense control upon infection with members of the Tevenvirinae subfamily. The Tevenvirinae subfamily is characterized by their unique cytosine modifications, which play a crucial role in their defense against RM systems like the T1RM (43). Specifically, Tequatrovirus group phages feature cytosines that are hydroxymethyl-glucosylated, while Mosigviruses possess cytosines that are hydroxymethyl-arabinosylated (10, 44).

Fig 3 Protection against Tevenvirinae by Tgv proteins encoded by the T1RM-embedded genes. (a) Observed defense activity against the BASEL phage collection. Protection levels (fold protection, as shown by the color code on the right) were determined by comparing plaque formation in strains with the system to those without, using 10-fold serial dilution assays. Data represent the average of two replicates. (b) Phage plaque assays on E. coli strains harboring an empty transposon (control, Ctrl) or the two T1RM-embedded genes (tgvAB), using a 10-fold serial dilution. (c) Defense against Tevenvirinae does not depend on HsdR’s restriction activity. Protection against members of the Tevenvirinae was assessed in E. coli strains expressing either the native five-gene cluster or only tgvAB, as shown in panel a. Additionally, strains containing the five-gene cluster with site-directed mutations hsdR[K172A] and tgvB[H571A] were evaluated to ascertain the individual contributions of the T1RM and TgvAB systems to defense capability. Data represent the average of three replicates. Details as in panel a. (d, e) Growth curves of E. coli cultures carrying an empty transposon (no system) or TntgvAB (+ tgvAB), without (NO phage) or with exposure to phages, initiated at time 0 with various multiplicities of infection (MOIs) (0.2, 2, 5, or 10). (d) Tequatroviruses and (e) Mosigviruses were used for infection. The presented data are the average of three independent experiments (±SD, illustrated with error bars).

To determine which part of the RM cluster was responsible for this protection, we created E. coli strains that either carried the T1RM cluster or just the embedded two-gene cluster independently. Strikingly, this revealed that the two-gene cluster was responsible for this protection (Fig. 3a and b). Furthermore, the two genes did not provide protection when expressed individually, indicating a necessity for their combined action to achieve the observed anti-phage activity (Fig. 3a). For reasons explained below, we named these two genes as type I-embedded GmrSD-like system of VPI-2, tgvA (VC1767) and tgvB (VC1766).

Interestingly, the strain carrying only the tgvAB genes showed significantly less protection (between 10- and 100-fold reduction) against several tested phages compared to the strain harboring the entire gene cluster (Fig. 3a and c). To investigate if restriction by the T1RM system contributes to this increased protection, we engineered a variant of the five-gene cluster containing a site-directed mutant of hsdR designed to inactivate restriction. Specifically, we replaced the catalytic lysine in the PD-(D/E)XK nuclease motif of the encoded HsdR protein with alanine (HsdR[K172A]), a modification known to abolish restriction activity in other T1RM systems (45–49). Notably, this mutated construct displayed a similar protection pattern to the unmodified gene cluster (Fig. 3c), indicating that the T1RM’s restriction activity does not play a direct role in the anti-phage defense. A plausible explanation for the reduced effectiveness of the standalone tgvAB operon, compared to its performance within the entire gene cluster, could be an imbalance in the stoichiometry of TgvA and TgvB proteins. This imbalance might be caused by the use of the artificial PBAD promoter and a non-native ribosome-binding site upstream of tgvA in the two-gene operon only construct.

The reasons for the lack of protection against the BASEL phages by the T1RM system, despite the presence of the recognition motif in 95.6% of these phages (Table S2), might be attributable to multiple factors. RM systems are the most prevalent defense mechanisms in bacterial genomes (5). As a result, many phages have evolved counter-defense mechanisms. For instance, phage T7 produces an Ocr (overcomes classical restriction) protein that mimics DNA to sequester RM enzymes, thereby preventing the restriction of its own DNA (50–52). Furthermore, the initial characterization of the BASEL phage collection by Maffei et al. showed that the T1RM systems tested were largely ineffective against phages outside the Drexlerviridae family (43). Closer inspection of their findings showed that fewer than 10% of the phages in the BASEL collection were significantly restricted (i.e., by a factor of 10 or more) by the tested T1RM systems, despite approximately 80% of the phages’ genomes containing the cognate recognition sites. Thus, our findings align with these previous observations.

To dissect the underlying mechanism of anti-phage defense by TgvAB, we monitored the growth kinetics of E. coli strains infected with increasing multiplicities of infection (MOIs) for both Tequatrovirus (Fig. 3d) and Mosigvirus phages (Fig. 3e). As expected, cultures of the no system control strain grew and then lysed in an MOI-dependent manner (Fig. 3d and e). In contrast, TgvAB producing cultures infected with the Tequatrovirus Bas35 continued to grow at rates indistinguishable from those of the no phage control up to and including MOI 5, before being partially overcome at MOI 10 (Fig. 3d). This phenotype is consistent with TgvAB acting directly to target the invading phage. However, TgvAB producing cultures infected with either the Tequatrovirus Bas40 or the Mosigviruses Bas46 and 47 all showed more variable levels of protection (Fig. 3d and e). Indeed, while protection was robust at MOI 0.2, at higher MOIs, we observed growth inhibition and even partial lysis. Nevertheless, given that the cultures mostly continued to grow past the point at which they lysed in the no system control, together with the direct protection observed against Bas35 at all tested MOIs, we conclude that TgvAB likely also acts directly against these phages, but that they are better able to overwhelm the system at high MOI.

The TgvAB defense system is a member of the GmrSD family of type IV restriction enzymes

Bioinformatic analysis of the TgvAB system revealed that TgvA (VC1767) and TgvB (VC1766) both possess an N-terminal DUF262 domain, while TgvB additionally contains a C-terminal DUF1524 domain (Fig. 4a and b). Interestingly, previous work by Machnicka et al. found that GmrS and GmrD proteins contain the DUF262 and DUF1524 domains, respectively, typically coming together to form GmrSD fusion proteins (53). Notably, the TgvB homolog from classical biotype V. cholerae (VC0395_A1364) was also identified as a GmrSD homolog in this study (53). These double domain forms of GmrSD function as modification-dependent type IV restriction enzymes, and are known to specifically recognize and cleave DNA containing sugar-modified hydroxymethyl cytosines. However, they exhibit no activity against unmodified DNA (53–56). Given that such modifications are typical of the Tevenvirinae (10) and the specific protective effect we observed against them (Fig. 3), this suggests that TgvAB may function in a similar manner. Importantly, and in contrast to classical single protein GmrSD such as Eco94GmrSD (Fig. 4a) (54), our phage infection assay revealed that TgvA and TgvB cannot function independently, and that both proteins are required for anti-phage activity.

Fig 4 The two-protein TgvAB defense system is a member of the GmrSD family of type IV restriction enzymes. (a) Structural models of Eco94GmrSD of E. coli STEC_94C and TgvA (VC1767) and TgvB (VC1766) of V. cholerae 7PET strains. The models, produced via AlphaFold (ColabFold), portray the domains with corresponding colors, while also highlighting the residues characteristic to the DUF262 and DUF1524 domains. Images were generated using ChimeraX 1.7.1. (b) Schematics displaying conserved domains identified in the TgvA and TgvB proteins. (c, d) Zoomed view of the conserved (c) DGQQR motif found in the DUF262 region of TgvA and TgvB and (d) of the His-Me finger motif within the DUF1524 of TgvB, highlighting the catalytic histidine (H) situated at the terminus of the β1 strand, the Asparagine (N) residue positioned in the loop region and the final N residue within the α-helix. (e) Site-directed mutagenesis removed the antiviral effect. The level of protection was evaluated as described in Fig. 3. Mutagenesis aimed at disrupting NTPase or endonuclease functions exerted by DUF262 and DUF1524, respectively. The data are averages from three independent experiments (±SD, as shown by the error bars).

Machnicka et al. showed that the predominant form of GmrSD is as a single multi-domain protein containing an N-terminal DUF262(GmrS) domain and a C-terminal DUF1524(GmrD) domain, separated by an alpha helical linker region (53). This domain organization was subsequently confirmed by crystal structures of the related GmrSD family members BrxU, which also recognizes and degrades DNA containing modified cytosines, and the phosphorothioate modification sensing enzyme SspE (57–59). Furthermore, biochemical experiments with these enzymes have shown that the N-terminal DUF262 likely functions as DNA modification sensor, and uses nucleotide binding and hydrolysis to regulate the activity of the C-terminal DUF1524, which functions as a nuclease to degrade non-self DNA (57, 59). Strikingly, structural modeling of Eco94GmrSD and TgvAB using AlphaFold (60) revealed that TgvB is predicted to share a similar domain architecture, although in the case of TgvA, this similarity is limited to the N-terminal DUF262 domain (Fig. 4a and b). Moreover, the top hits in structural alignments of the TgvAB models were SspE and BrxU, reinforcing the idea that these proteins are related.

Next, to further investigate the relative contributions of the DUF262 and DUF1524 domains to TgvAB function, we used the structural modeling and alignments to identify key residues in each domain. For both TgvA and B, the three highly conserved motifs characteristically associated with the DUF262 domain [i.e., (i) QR, (ii) DGQQR, and (iii) FxxxN] were readily identifiable (Fig. 4a through d) (53). Notably, the DGQQR motif is thought to form part of a nucleotide-binding pocket and to be required for nucleotide hydrolysis. Indeed, site-directed mutants of either TgvA or TgvB encoding substitutions in this motif previously shown to disrupt NTPase activity (57–59) all resulted in a total loss of anti-phage activity (Fig. 4c through e). In contrast, the DUF1524 domain contains a highly conserved H…N…H/N motif, which belongs to the His-Me finger nuclease superfamily and that assumes a characteristic ββα fold (61, 62). Such a motif was readily apparent in the C-terminal domain of the predicted TgvB structure, and consistent with previous findings (54, 57–59), substitutions designed to disrupt either the catalytic histidine (TgvB[H571A]) or the metal-binding asparagine (TgvB[N602A]) were sufficient to abolish anti-phage activity (Fig. 3c, 4d and e).

Overall, our results suggest that the TgvAB system senses phages with hypermodified cytosines in a manner that requires the DUF262 domains of both TgvA and B, and that the His-Me nuclease domain of TgvB likely functions as the effector against phage DNA. The notion that TgvAB recognizes modified DNA aligns with findings from Gomez and Waters, who concurrently described the TgvAB system in their study (63). They demonstrated that T2/T4 phages lacking functional glucosyltransferase genes, essential for the glucosylation of hydroxymethylated cytosines, could evade the defense exerted by TgvAB (63). This further supports the specificity of TgvAB toward modified DNA. Nevertheless, why TgvB alone is not sufficient for phage protection remains unclear. One possibility is that TgvA is required to overcome a phage encoded inhibitor. For example, some GmrSD family enzymes such as Eco94GmrSD are inhibited by the protein IPI*, which is co-injected into the host cell with the T4 genome (54, 55, 64). However, TgvA could equally also play a regulatory or structural role, and further work will therefore be needed to clarify these possibilities.

Occurrence of the tgvAB system within and outside T1RM clusters

To investigate the prevalence of tgvAB homologs within the T1RM cluster, we examined the distribution of the specific five-gene operon within 41,142 bacterial genomes (see Materials and Methods for details). This in silico analysis revealed that the gene architecture found in VPI-2 of V. cholerae is also present in a variety of other bacterial genera (Fig. 5a) with 79 hits within this genome database, including Shewanella, Acinetobacter, and Pseudoalteromonas species (see Table S3 for species-level details). This wider distribution indicates the potential functional conservation of these gene arrangements across different gram-negative bacteria. However, the genus Vibrio was still most prominently featured in these findings with 52 hits (Fig. 5a). Precisely, apart from V. cholerae, species such as Vibrio vulnificus, Vibrio antiquarius, Vibrio nigripulchritudo, Vibrio parahaemolyticus, and the unclassified Vibrio strain B1ASS3 (Vibrio sp.) were identified to carry similar gene clusters (Fig. 5a). Despite the presence of these diverse Vibrio species, V. cholerae 7PET strains were the most commonly identified with 38 hits (Fig. 5a), likely reflecting their prominent representation in the NCBI database.

Fig 5 Phylogenetic distribution of the restriction systems. (a) The presence of the five-gene cluster (VC1765–69) was assessed across 41,142 bacterial genomes. The results revealed its distribution beyond V. cholerae, which represented 65.8% of all hits. *V. cholerae O37 serogroup strains are known to be closely related to classical O1 strains with highly similar chromosomal backbones. (b, c) Exploration of the (b) T1RM system (VC1769-68-65) and (c) TgvAB system (VC1767–66) across the bacterial genomes demonstrates their assessment at the order level of taxonomy. Orders represented in less than 1% of instances were consolidated into a singular category labeled “Others” for the visualization. For details at the species level, see Tables S3 to S5.

Subsequent analysis focused on the independent occurrences of the T1RM and TgvAB systems. As expected, the T1RM system was widespread (2,808 hits) across numerous bacterial orders (Fig. 5b and Table S4 for species-level details). Homologs of the tgvAB operon alone were slightly less common with 1,341 hits (Table S5 for species-level details), yet 17 times more prevalent than the instances of the five-gene operon described above. Indeed, as shown in Fig. 5c, the occurrence of TgvAB homologs spans a wide array of bacterial orders, with species found in the human gut, like Bacteroides fragilis (Bacteroidales), to organisms isolated from permafrost, such as Psychrobacter cryohalolentis (Moraxellales).

Conclusion

In this study, we aimed to characterize the predicted restriction gene cluster of VPI-2. We showed that the T1RM system actively methylates the genomes of 7PET V. cholerae strains, while restricting unmethylated foreign DNA. Additionally, we identified a novel two-protein modification-dependent restriction system, TgvAB, which is embedded within the T1RM cluster. Interestingly, Picton et al. demonstrated that the TgvB homolog BrxU, along with the bacteriophage exclusion (BREX) system (65), works in concert to offer complementary resistance against both modified and non-modified phages (53, 57). Therefore, it is tempting to speculate that the embedding of the tgvAB operon within the T1RM cluster serves a similar complementary role in V. cholerae. Supporting this notion, Machnicka et al. noted that GmrSD homologs are frequently encoded within type I RM loci. An example includes the gene encoding the DUF262 domain-containing protein RloF of Campylobacter jejuni, which is situated between hsdR and hsdS of a T1RM operon (66), similar to the positioning of tgvAB described in this study. That defense systems tend to cluster together within defense islands has been established over several years (11, 12, 67). However, this concept was recently extended by Payne and colleagues by identifying specific genes embedded within multi-gene defense clusters, highlighting the complex organization and integration of these systems within bacterial genomes (68). Notably, their research found GmrSD-like genes embedded within Hma (helicase, methylase, ATPase) defense gene clusters. However, unlike the HEC-05 (= BrxU) and HEC-06 GmrSD-like proteins identified in their work, which function independently (57, 68), our findings indicate that the TgvAB defense operates as a two-protein system, underscoring the diversity in bacterial defense strategies.

MATERIALS AND METHODS

Bacterial strains, plasmids, and culture conditions

The bacterial strains and the plasmids used in this study are listed in Table S1. pUC18-mini-Tn7T-Gm-lacZ was a gift from Herbert Schweizer via Addgene plasmid #63120 (69). The primary V. cholerae strain used, A1552, is a fully sequenced toxigenic O1 El Tor Inaba strain, representing the ongoing seventh cholera pandemic (70, 71). Unless stated otherwise, bacteria were aerobically cultured in lysogeny broth (LB; 1% tryptone, 0.5% yeast extract, 1% sodium chloride; Carl Roth, Switzerland) with shaking at 180 rpm, or on LB agar plates at either 30°C or 37°C. When required, antibiotic selection was applied using ampicillin (100 µg/mL), kanamycin (75 µg/mL), and gentamicin (25 or 50 µg/mL). For natural transformation, chitin powder (Alfa Aesar via Thermo Fisher, USA) was combined with half-concentrated Instant Ocean medium (Aquarium Systems) and sterilized by autoclaving prior to adding the bacterial cultures.

Conjugation with MFDpir (72) was used to introduce the mini-Tn7 transposon derivatives into E. coli strain MG1655 on agar plates containing 0.3 mM diaminopimelic acid (Sigma-Aldrich). To induce expression from the PBAD promoter, cultures were grown in media containing 0.2% L-arabinose. For bacteriophages experiments, LB medium was supplemented with 5 mM CaCl2 + 20 mM MgSO4. Double-layer LB plates were prepared by adding 0.5% agar for semi-solid agar and 1.5% agar for the solid base.

Genetic engineering of strains and plasmids

Standard molecular cloning techniques were utilized for the cloning process (73) using the following enzymes: Pwo polymerase (Roche), Q5 High-Fidelity Polymerase (New England Biolabs), GoTaq Polymerase (Promega), restriction enzymes (New England Biolabs), and T4 DNA ligase (New England Biolabs). Enzymes were used according to the manufacturer’s instructions. All constructs were verified through PCR and/or Sanger or Nanopore sequencing (performed by Microsynth AG, Switzerland) and analyzed using SnapGene version 4.3.11.

V. cholerae strains were created through natural transformation and Flp recombination (TransFLP) (74–76) or through allelic exchange using derivatives of the suicide plasmid pGP704-Sac28 (77) and SacB-based counter-selection on NaCl-free LB plates with 10% sucrose. Mini-Tn7 transposons, containing araC and the gene(s) of interest regulated by the arabinose-inducible promoter PBAD, were inserted in E. coli downstream of glmS via triparental mating, following established protocols (78). Site-directed mutations in these constructs were introduced by inverse PCR prior to their transposition into the E. coli chromosome.

PacBio (SMRT) sequencing

Genomic DNA was purified from overnight cultures using Qiagen’s Genomic-tip procedure combined with the Genomic DNA buffer set (Qiagen, Switzerland), following the manufacturer’s instructions. Sample processing, PacBio Single Molecule, Real-Time (SMRT) sequencing, and de novo genome assembly were performed at the University of Lausanne’s Genomic Technology Facility, as previously described (34). Note that the assembled genomes of strains A1552, C6706, C6709, P27459, E7946, DRC193A, and MO10 have been previously reported without analysis of their epigenetic modifications (34, 35, 71).

Electroporation-mediated transformation of V. cholerae using plasmids

To explore the T1RM system’s efficiency in restricting DNA with specific recognition sequences, we compared the uptake frequency of a plasmid harboring the putative recognition motif (Pmotif+) to that of a variant plasmid with silent mutations in aacC1 (Pmotif-) altering its sequence while maintaining the encoded aminoglycoside-3-O-acetyltransferase-I protein. Transformation frequencies were assessed through electroporation. V. cholerae competent cells were prepared by standard protocols (73), involving 1:100 dilution of overnight cultures, growth for 2 h and 30 min at 37°C (until an optical density at 600 nm [OD600] of ~1.0), and washing steps with cold 2 mM CaCl2 and 10% glycerol before shock-freezing. After 2 h at −80°C, electroporation with 300 ng plasmid was performed at 1.6 kV followed by recovery in 2xYT-rich medium at 30°C for 2  h. Cells were plated on LB agar with and without kanamycin and incubated at 37°C overnight. Transformation frequencies were calculated as the ratio of kanamycin-resistant transformants to the total number of bacteria.

Bacteriophage handling and culturing

The E. coli BASEL phage collection (43) was used in this study. To generate phage stocks, an E. coli MG1655ΔaraCBAD (79) overnight culture was diluted and grown to the exponential phase in LB medium supplemented with 5 mM CaCl2 and 20 mM MgSO4. Subsequently, the culture was 1:10 diluted in prewarmed medium, infected with 104 plaque-forming units/mL, and incubated under shaking conditions at 37°C for 5 h. Following incubation, centrifugation and filtration were used to clear the lysate, which was then treated with 1% chloroform and stored at +4°C. Phage titers were determined using plaque assays on the propagation strain.

Bacteriophage plaque assays

For plaque assays, E. coli MG1655ΔaraCBAD, either with the candidate defense system or the empty mini-Tn7 transposon control, was grown in LB medium. Overnight cultures were diluted 1:100 in LB medium supplemented with 0.2% arabinose, 5 mM CaCl2, and 20 mM MgSO4 and grown at 37°C with shaking for 2 h. Once reaching the exponential phase, the cultures were diluted 1:40 in 0.5% LB agar containing 5 mM CaCl2, 20 mM MgSO4, and 0.2% arabinose, then overlaid on 1.5% LB agar. Phage samples were serial diluted in LB medium with 5 mM CaCl2 and 20 mM MgSO4 and spotted onto the bacterial overlays. After overnight incubation at 37°C, plaques were counted to assess the defense system’s effectiveness compared to the mini-Tn7-carrying control strain (= fold protection).

Infection kinetics

The infection kinetics assay of Tequatroviruses (Bas35, Bas40) and Mosigviruses (Bas46, Bas47) was conducted as follows: overnight cultures of E. coli strains were diluted 1:100 in LB medium supplemented with 5 mM CaCl2, 20 mM MgSO4, and 0.2% arabinose. Bacterial cultures were then incubated at 37°C with shaking for 2 h. Subsequently, 20 µL of phage per well at MOIs of 0, 0.2, 5, or 10 was added in technical triplicate to a 96-well plate. The cultures were further diluted 1:10 in the same LB condition, and 180 µL of each diluted culture was then added to the wells. The SpectraMax i3x plate reader from Molecular Devices was utilized to assess bacterial growth at 37°C, with measurements taken at 6-minute intervals over a total of 49 cycles. To calculate the MOI, cultures of strains MG1655∆araCBAD-TnAraC (no system control) and MG1655∆araCBAD-TnTgvAB were cultured following the protocol outlined in the Bacteriophage Plaque Assays section. Colony-forming units (per milliliter) were quantified by spotting serially diluted cultures onto LB plates. The calculated values represent the average of three technical replicates.

Bioinformatics analyses

The VPI-2 genomic region of 7PET O1 strains and one O139 serogroup strain (MO10) was compared and visualized using Clinker software (v.0.0.25, default parameters) (80) after reannotation of the genome sequence of strain A1552 using the Prokaryotic Genome Annotation pipeline version 2023-10-03.build7061 (81) to unify the annotation method. For sequence similarity, NCBI’s blastp was utilized (default parameters, non-redundant protein database; accession August 2023), while structural modeling was conducted with ColabFold (1.5.2) (82) based on AlphaFold2 (60) using default settings. The DALI server was employed for structural similarity predictions against the Protein Data Bank (83, 84).

The distribution of the specific five-gene operon (VC1769–65) across bacterial species was examined with MacSyFinder (v.2.1) (85), using a comprehensive database of sequenced and fully assembled bacterial genomes (taxid:2) from the NCBI database (accession date 26 January 2024; GenBank database with all complete and chromosome level assemblies using the data sets utility from the NCBI command line tool). This analysis therefore covered a data set comprising 41,142 bacterial genomes, which altogether contained over 154 million protein sequences.

To build hidden Markov model (HMM) profiles for each target coding sequence within the VC1769–65 operon, homologous protein sequences were identified via PSI-BLAST searches in the NCBI database (three iterations), using the non-redundant protein sequence database (accessed in February 2024) with a cutoff e-value of 1e − 10.

After identifying homologous sequences for each CDS through PSI-BLAST, the sequences were aligned using MAFFT (v.7.508, --maxiterate 1000 –localpair parameters for higher accuracy alignments) (86). From these multiple alignments, HMM profiles were generated with HMMER (v.3.3.2, using hmmbuild with default parameters) (87), forming the basis for constructing different models in MacSyFinder. These models were used to search for the occurrence of the CDS in various combinations encompassing the T1RM and/or TgvAB system genes. The constructed models were then applied in a search across the bacterial protein database mentioned above.

Statistics and reproducibility

Results are derived from biologically independent experiments, as specified in the figure legends. Statistical analyses were conducted using Prism software (v.10.2.1; GraphPad).

ACKNOWLEDGMENTS

We thank members of the Blokesch lab and especially Nicolas Flaugnatti and Alexis Proutière for fruitful discussions and Sandrine Stutzmann, Laurie Righi, and Loriane Bader for technical assistance. We acknowledge the staff of the Lausanne Genomic Technologies Facility and Christian Iseli & Nicolas Guex from the EPFL/UNIL Bioinformatics Competence Center for sample processing, sequencing, genome assembly, and SMRT sequencing analysis. We are grateful to Alexander Harms for sharing the BASEL phage collection and for valuable discussions.

This work was supported by the Swiss National Science Foundation (310030_185022) and an International Research Scholarship by the Howard Hughes Medical Institute (HHMI) (55008726) awarded to M.B.

M.B. secured funding; M.B. and D.W.A. conceived and supervised the study; G.V., D.W.A., and M.B. designed the experiments and analyzed the data. G.V., D.W.A., and M.B. designed and constructed strains/plasmids. G.V. performed the experiments. M.B. initiated the SMRT sequencing. A.L. performed comparative genomic and conservation analyses. G.V., D.W.A., and M.B. wrote and revised the manuscript with input from A.L. All authors approved the final version of the paper.

DATA AVAILABILITY

All SMRT sequencing raw data have been made available on Zenodo (three data sets: 10.5281/zenodo.10838595; 10.5281/zenodo.10839511; 10.5281/zenodo.10839547).

SUPPLEMENTAL MATERIAL

The following material is available online at https://doi.org/10.1128/jb.00145-24.

10.1128/jb.00145-24.SuF1 Table S1 jb.00145-24-s0001.docx

Bacterial strains and plasmids used in this study.

10.1128/jb.00145-24.SuF2 Table S2 jb.00145-24-s0002.xlsx

Counts of the T1RM recognition motif in phage genomes from the BASEL collection.

10.1128/jb.00145-24.SuF3 Table S3 jb.00145-24-s0003.xlsx

Summary of matching hits for the T1RM + TgvAB (VC1769-65) model, detected by MacSyFinder v.2.

10.1128/jb.00145-24.SuF4 Table S4 jb.00145-24-s0004.xlsx

Summary of matching hits for the T1RM (VC1769-68-65) model, detected by MacSyFinder v.2.

10.1128/jb.00145-24.SuF5 Table S5 jb.00145-24-s0005.xlsx

Summary of matching hits for the TgvAB (VC1767-66) model, detected by MacSyFinder v.2.

ASM does not own the copyrights to Supplemental Material that may be linked to, or accessed through, an article. The authors have granted ASM a non-exclusive, world-wide license to publish the Supplemental Material files. Please contact the corresponding author directly for reuse.
==== Refs
REFERENCES

1 Rocha EPC, Bikard D. 2022. Microbial defenses against mobile genetic elements and viruses: who defends whom from what? PLoS Biol 20 :e3001514. doi:10.1371/journal.pbio.3001514 35025885
2 Baltrus DA. 2013. Exploring the costs of horizontal gene transfer. Trends Ecol Evol 28 :489–495. doi:10.1016/j.tree.2013.04.002 23706556
3 Stern A, Sorek R. 2011. The phage-host arms race: shaping the evolution of microbes. Bioessays 33 :43–51. doi:10.1002/bies.201000071 20979102
4 Dy RL, Richter C, Salmond GPC, Fineran PC. 2014. Remarkable mechanisms in microbes to resist phage infections. Annu Rev Virol 1 :307–331. doi:10.1146/annurev-virology-031413-085500 26958724
5 Georjon H, Bernheim A. 2023. The highly diverse antiphage defence systems of bacteria. Nat Rev Microbiol 21 :686–700. doi:10.1038/s41579-023-00934-x 37460672
6 Lopatina A, Tal N, Sorek R. 2020. Abortive infection: bacterial suicide as an antiviral immune strategy. Annu Rev Virol 7 :371–384. doi:10.1146/annurev-virology-011620-040628 32559405
7 Murray NE. 2000. Type I restriction systems: sophisticated molecular machines (a legacy of Bertani and Weigle). Microbiol Mol Biol Rev 64 :412–434. doi:10.1128/MMBR.64.2.412-434.2000 10839821
8 Tock MR, Dryden DTF. 2005. The biology of restriction and anti-restriction. Curr Opin Microbiol 8 :466–472. doi:10.1016/j.mib.2005.06.003 15979932
9 Loenen WAM, Raleigh EA. 2014. The other face of restriction: modification-dependent enzymes. Nucleic Acids Res 42 :56–69. doi:10.1093/nar/gkt747 23990325
10 Weigele P, Raleigh EA. 2016. Biosynthesis and function of modified bases in bacteria and their viruses. Chem Rev 116 :12655–12687. doi:10.1021/acs.chemrev.6b00114 27319741
11 Makarova KS, Wolf YI, Snir S, Koonin EV. 2011. Defense islands in bacterial and archaeal genomes and prediction of novel defense systems. J Bacteriol 193 :6039–6056. doi:10.1128/JB.05535-11 21908672
12 Doron S, Melamed S, Ofir G, Leavitt A, Lopatina A, Keren M, Amitai G, Sorek R. 2018. Systematic discovery of antiphage defense systems in the microbial pangenome. Science 359 :eaar4120. doi:10.1126/science.aar4120 29371424
13 Hussain FA, Dubert J, Elsherbini J, Murphy M, VanInsberghe D, Arevalo P, Kauffman K, Rodino-Janeiro BK, Gavin H, Gomez A, Lopatina A, Le Roux F, Polz MF. 2021. Rapid evolutionary turnover of mobile genetic elements drives bacterial resistance to phages. Science 374 :488–492. doi:10.1126/science.abb1083 34672730
14 Waldor MK, Mekalanos JJ. 1996. Lysogenic conversion by a filamentous phage encoding cholera toxin. Science 272 :1910–1914. doi:10.1126/science.272.5270.1910 8658163
15 Taylor RK, Miller VL, Furlong DB, Mekalanos JJ. 1987. Use of phoA gene fusions to identify a pilus colonization factor coordinately regulated with cholera toxin. Proc Natl Acad Sci U S A 84 :2833–2837. doi:10.1073/pnas.84.9.2833 2883655
16 Karaolis DKR, Johnson JA, Bailey CC, Boedeker EC, Kaper JB, Reeves PR. 1998. A Vibrio cholerae pathogenicity island associated with epidemic and pandemic strains. Proc Natl Acad Sci U S A 95 :3134–3139. doi:10.1073/pnas.95.6.3134 9501228
17 Balasubramanian D, López-Pérez M, Almagro-Moreno S. 2023. Cholera dynamics and the emergence of pandemic Vibrio cholerae, p 127–147. In Almagro-Moreno S, Pukatzki S (ed), Vibrio spp. infections. Springer International Publishing, Cham.
18 Dziejman M, Balon E, Boyd D, Fraser CM, Heidelberg JF, Mekalanos JJ. 2002. Comparative genomic analysis of Vibrio cholerae: genes that correlate with cholera endemic and pandemic disease. Proc Natl Acad Sci U S A 99 :1556–1561. doi:10.1073/pnas.042667999 11818571
19 Davies BW, Bogard RW, Young TS, Mekalanos JJ. 2012. Coordinated regulation of accessory genetic elements produces cyclic di-nucleotides for V. cholerae virulence. Cell 149 :358–370. doi:10.1016/j.cell.2012.01.053 22500802
20 Severin GB, Ramliden MS, Hawver LA, Wang K, Pell ME, Kieninger A-K, Khataokar A, O’Hara BJ, Behrmann LV, Neiditch MB, Benning C, Waters CM, Ng W-L. 2018. Direct activation of a phospholipase by cyclic GMP-AMP in El Tor Vibrio cholerae. Proc Natl Acad Sci U S A 115 :E6048–E6055. doi:10.1073/pnas.1801233115 29891656
21 Cohen D, Melamed S, Millman A, Shulman G, Oppenheimer-Shaanan Y, Kacen A, Doron S, Amitai G, Sorek R. 2019. Cyclic GMP–AMP signalling protects bacteria against viral infection. Nature 574 :691–695. doi:10.1038/s41586-019-1605-5 31533127
22 Hsueh BY, Severin GB, Elg CA, Waldron EJ, Kant A, Wessel AJ, Dover JA, Rhoades CR, Ridenhour BJ, Parent KN, Neiditch MB, Ravi J, Top EM, Waters CM. 2022. Phage defence by deaminase-mediated depletion of deoxynucleotides in bacteria. Nat Microbiol 7 :1210–1220. doi:10.1038/s41564-022-01162-4 35817890
23 Jaskólska M, Adams DW, Blokesch M. 2022. Two defence systems eliminate plasmids from seventh pandemic Vibrio cholerae. Nature 604 :323–329. doi:10.1038/s41586-022-04546-y 35388218
24 Jermyn WS, Boyd EF. 2002. Characterization of a novel Vibrio pathogenicity island (VPI-2) encoding neuraminidase (nanH) among toxigenic Vibrio cholerae isolates. Microbiology (Reading) 148 :3681–3693. doi:10.1099/00221287-148-11-3681 12427958
25 Almagro-Moreno S, Boyd EF. 2009. Sialic acid catabolism confers a competitive advantage to pathogenic Vibrio cholerae in the mouse intestine. Infect Immun 77 :3807–3816. doi:10.1128/IAI.00279-09 19564383
26 Yan J, Liu Q, Xue X, Li J, Li Y, Su Y, Cao B. 2023. The response regulator VC1795 of Vibrio pathogenicity island-2 contributes to intestinal colonization by Vibrio cholerae. Int J Mol Sci 24 :13523. doi:10.3390/ijms241713523 37686329
27 Hu H, Hughes TCD, Popp PF, Roa-Eguiara A, Martin FJO, Rutbeek NR, Hendriks IA, Payne LJ, Yan Y, Sousa VK de, Wang Y, Nielsen ML, Berry RM, Erhardt M, Jackson SA, Taylor NMI. 2023. Structure and mechanism of Zorya anti-phage defense system. bioRxiv. doi:10.1101/2023.12.18.572097
28 Jermyn WS, Boyd EF. 2005. Molecular evolution of Vibrio pathogenicity island-2 (VPI-2): mosaic structure among Vibrio cholerae and Vibrio mimicus natural isolates. Microbiology (Reading) 151 :311–322. doi:10.1099/mic.0.27621-0 15632448
29 Heidelberg JF, Eisen JA, Nelson WC, Clayton RA, Gwinn ML, Dodson RJ, Haft DH, Hickey EK, Peterson JD, Umayam L, et al. . 2000. DNA sequence of both chromosomes of the cholera pathogen Vibrio cholerae. Nature 406 :477–483. doi:10.1038/35020000 10952301
30 Dorman MJ, Domman D, Uddin MI, Sharmin S, Afrad MH, Begum YA, Qadri F, Thomson NR. 2019. High quality reference genomes for toxigenic and non-toxigenic Vibrio cholerae serogroup O139. Sci Rep 9 :5865. doi:10.1038/s41598-019-41883-x 30971707
31 Loenen WAM, Dryden DTF, Raleigh EA, Wilson GG. 2014. Type I restriction enzymes and their relatives. Nucleic Acids Res 42 :20–44. doi:10.1093/nar/gkt847 24068554
32 Dussoix D, Arber W. 1962. Host specificity of DNA produced by Escherichia coli. J Mol Biol 5 :37–49. doi:10.1016/S0022-2836(62)80059-X 13888713
33 Arber W. 1965. Host-controlled modification of bacteriophage. Annu Rev Microbiol 19 :365–378. doi:10.1146/annurev.mi.19.100165.002053 5318444
34 Stutzmann S, Blokesch M. 2020. Comparison of chitin-induced natural transformation in pandemic Vibrio cholerae O1 El Tor strains. Environ Microbiol 22 :4149–4166. doi:10.1111/1462-2920.15214 32860313
35 Lemopoulos A, Miwanda B, Drebes Dörr NC, Stutzmann S, Bompangue D, Muyembe-Tamfum J-J, Blokesch M. 2024. Genome sequences of Vibrio cholerae strains isolated in the DRC between 2009 and 2012. Microbiol Resour Announc 13 :e0082723. doi:10.1128/mra.00827-23 38345380
36 Beaulaurier J, Schadt EE, Fang G. 2019. Deciphering bacterial epigenomes using modern sequencing technologies. Nat Rev Genet 20 :157–172. doi:10.1038/s41576-018-0081-3 30546107
37 Hennig S, Ziebuhr W. 2010. Characterization of the transposase encoded by IS256, the prototype of a major family of bacterial insertion sequence elements. J Bacteriol 192 :4153–4163. doi:10.1128/JB.00226-10 20543074
38 Robins WP, Meader BT, Toska J, Mekalanos JJ. 2022. Cell density-dependent death triggered by viral palindromic DNA sequences. bioRxiv. doi:10.1101/2022.11.18.517080
39 Bickle TA, Krüger DH. 1993. Biology of DNA restriction. Microbiol Rev 57 :434–450. doi:10.1128/mr.57.2.434-450.1993 8336674
40 Seed KD, Bodi KL, Kropinski AM, Ackermann H-W, Calderwood SB, Qadri F, Camilli A. 2011. Evidence of a dominant lineage of Vibrio cholerae-specific lytic bacteriophages shed by cholera patients over a 10-year period in Dhaka, Bangladesh. mBio 2 :e00334-10. doi:10.1128/mBio.00334-10 21304168
41 Boyd CM, Angermeyer A, Hays SG, Barth ZK, Patel KM, Seed KD. 2021. Bacteriophage ICP1: a persistent predator of Vibrio cholerae. Annu Rev Virol 8 :285–304. doi:10.1146/annurev-virology-091919-072020 34314595
42 LeGault KN, Hays SG, Angermeyer A, McKitterick AC, Johura F-T, Sultana M, Ahmed T, Alam M, Seed KD. 2021. Temporal shifts in antibiotic resistance elements govern phage-pathogen conflicts. Science 373 :eabg2166. doi:10.1126/science.abg2166 34326207
43 Maffei E, Shaidullina A, Burkolter M, Heyer Y, Estermann F, Druelle V, Sauer P, Willi L, Michaelis S, Hilbi H, Thaler DS, Harms A. 2021. Systematic exploration of Escherichia coli phage-host interactions with the BASEL phage collection. PLoS Biol 19 :e3001424. doi:10.1371/journal.pbio.3001424 34784345
44 Thomas JA, Orwenyo J, Wang L-X, Black LW. 2018. The odd “RB” phage—identification of arabinosylation as a new epigenetic modification of DNA in T4-like phage RB69. Viruses 10 :313. doi:10.3390/v10060313 29890699
45 Steczkiewicz K, Muszewska A, Knizewski L, Rychlewski L, Ginalski K. 2012. Sequence, structure and functional diversity of PD-(D/E)XK phosphodiesterase superfamily. Nucleic Acids Res 40 :7016–7045. doi:10.1093/nar/gks382 22638584
46 Obarska-Kosinska A, Taylor JEN, Callow P, Orlowski J, Bujnicki JM, Kneale GG. 2008. HsdR subunit of the type I restriction-modification enzyme EcoR124I: biophysical characterisation and structural modelling. J Mol Biol 376 :438–452. doi:10.1016/j.jmb.2007.11.024 18164032
47 Uyen NT, Park S-Y, Choi J-W, Lee H-J, Nishi K, Kim J-S. 2009. The fragment structure of a putative HsdR subunit of a type I restriction enzyme from Vibrio vulnificus YJ016: implications for DNA restriction and translocation activity. Nucleic Acids Res 37 :6960–6969. doi:10.1093/nar/gkp603 19625490
48 Grinkevich P, Sinha D, Iermak I, Guzanova A, Weiserova M, Ludwig J, Mesters JR, Ettrich RH. 2018. Crystal structure of a novel domain of the motor subunit of the Type I restriction enzyme EcoR124 involved in complex assembly and DNA binding. J Biol Chem 293 :15043–15054. doi:10.1074/jbc.RA118.003978 30054276
49 Davies GP, Martin I, Sturrock SS, Cronshaw A, Murray NE, Dryden DTF. 1999. On the structure and operation of type I DNA restriction enzymes. J Mol Biol 290 :565–579. doi:10.1006/jmbi.1999.2908 10390354
50 Studier FW. 1975. Gene 0.3 of bacteriophage T7 acts to overcome the DNA restriction system of the host. J Mol Biol 94 :283–295. doi:10.1016/0022-2836(75)90083-2 1095770
51 Bandyopadhyay PK, Studier FW, Hamilton DL, Yuan R. 1985. Inhibition of the type I restriction-modification enzymes EcoB and EcoK by the gene 0.3 protein of bacteriophage T7. J Mol Biol 182 :567–578. doi:10.1016/0022-2836(85)90242-6 2989534
52 Atanasiu C, Su T-J, Sturrock SS, Dryden DTF. 2002. Interaction of the ocr gene 0.3 protein of bacteriophage T7 with EcoKI restriction/modification enzyme. Nucleic Acids Res 30 :3936–3944. doi:10.1093/nar/gkf518 12235377
53 Machnicka MA, Kaminska KH, Dunin-Horkawicz S, Bujnicki JM. 2015. Phylogenomics and sequence-structure-function relationships in the GmrSD family of Type IV restriction enzymes. BMC Bioinformatics 16 :336. doi:10.1186/s12859-015-0773-z 26493560
54 He X, Hull V, Thomas JA, Fu X, Gidwani S, Gupta YK, Black LW, Xu S. 2015. Expression and purification of a single-chain Type IV restriction enzyme Eco94GmrSD and determination of its substrate preference. Sci Rep 5 :9747. doi:10.1038/srep09747 25988532
55 Bair CL, Rifat D, Black LW. 2007. Exclusion of glucosyl-hydroxymethylcytosine DNA containing bacteriophages is overcome by the injected protein inhibitor IPI*. J Mol Biol 366 :779–789. doi:10.1016/j.jmb.2006.11.049 17188711
56 Bair CL, Black LW. 2007. A type IV modification dependent restriction nuclease that targets glucosylated hydroxymethyl cytosine modified DNAs. J Mol Biol 366 :768–778. doi:10.1016/j.jmb.2006.11.051 17188297
57 Picton DM, Luyten YA, Morgan RD, Nelson A, Smith DL, Dryden DTF, Hinton JCD, Blower TR. 2021. The phage defence island of a multidrug resistant plasmid uses both BREX and type IV restriction for complementary protection from viruses. Nucleic Acids Res 49 :11257–11273. doi:10.1093/nar/gkab906 34657954
58 Xiong X, Wu G, Wei Y, Liu L, Zhang Y, Su R, Jiang X, Li M, Gao H, Tian X, Zhang Y, Hu L, Chen S, Tang Y, Jiang S, Huang R, Li Z, Wang Y, Deng Z, Wang J, Dedon PC, Chen S, Wang L. 2020. SspABCD–SspE is a phosphorothioation-sensing bacterial defence system with broad anti-phage activities. Nat Microbiol 5 :917–928. doi:10.1038/s41564-020-0700-6 32251370
59 Gao H, Gong X, Zhou J, Zhang Y, Duan J, Wei Y, Chen L, Deng Z, Wang J, Chen S, Wu G, Wang L. 2022. Nicking mechanism underlying the DNA phosphorothioate-sensing antiphage defense by SspE. Nat Commun 13 :6773. doi:10.1038/s41467-022-34505-0 36351933
60 Jumper J, Evans R, Pritzel A, Green T, Figurnov M, Ronneberger O, Tunyasuvunakool K, Bates R, Žídek A, Potapenko A, et al. . 2021. Highly accurate protein structure prediction with AlphaFold. Nature 596 :583–589. doi:10.1038/s41586-021-03819-2 34265844
61 Jablonska J, Matelska D, Steczkiewicz K, Ginalski K. 2017. Systematic classification of the His-Me finger superfamily. Nucleic Acids Res 45 :11479–11494. doi:10.1093/nar/gkx924 29040665
62 Wu C-C, Lin JLJ, Yuan HS. 2020. Structures, mechanisms, and functions of His-Me finger nucleases. Trends Biochem Sci 45 :935–946. doi:10.1016/j.tibs.2020.07.002 32807610
63 Gomez JB, Waters CM. 2024. A Vibrio cholerae Type IV restriction system targets glucosylated 5-hydroxyl methyl cytosine to protect against phage infection. bioRxiv:2024.04.05.588314. doi:10.1101/2024.04.05.588314
64 Rifat D, Wright NT, Varney KM, Weber DJ, Black LW. 2008. Restriction endonuclease inhibitor IPI* of bacteriophage T4: a novel structure for a dedicated target. J Mol Biol 375 :720–734. doi:10.1016/j.jmb.2007.10.064 18037438
65 Goldfarb T, Sberro H, Weinstock E, Cohen O, Doron S, Charpak-Amikam Y, Afik S, Ofir G, Sorek R. 2015. BREX is a novel phage resistance system widespread in microbial genomes. EMBO J 34 :169–183. doi:10.15252/embj.201489455 25452498
66 Miller WG, Pearson BM, Wells JM, Parker CT, Kapitonov VV, Mandrell RE. 2005. Diversity within the Campylobacter jejuni type I restriction-modification loci. Microbiology (Reading) 151 :337–351. doi:10.1099/mic.0.27327-0 15699185
67 Hochhauser D, Millman A, Sorek R. 2023. The defense island repertoire of the Escherichia coli pan-genome. PLoS Genet 19 :e1010694. doi:10.1371/journal.pgen.1010694 37023146
68 Payne LJ, Hughes TCD, Fineran PC, Jackson SA. 2024. New antiviral defences are genetically embedded within prokaryotic immune systems. bioRxiv. doi:10.1101/2024.01.29.577857
69 Choi K-H, Gaynor JB, White KG, Lopez C, Bosio CM, Karkhoff-Schweizer RR, Schweizer HP. 2005. A Tn7-based broad-range bacterial cloning and expression system. Nat Methods 2 :443–448. doi:10.1038/nmeth765 15908923
70 Yildiz FH, Schoolnik GK. 1998. Role of rpoS in stress survival and virulence of Vibrio cholerae. J Bacteriol 180 :773–784. doi:10.1128/JB.180.4.773-784.1998 9473029
71 Matthey N, Drebes Dörr NC, Blokesch M. 2018. Long-read-based genome sequences of pandemic and environmental Vibrio cholerae strains. Microbiol Resour Announc 7 :e01574-18. doi:10.1128/MRA.01574-18 30574591
72 Ferrières L, Hémery G, Nham T, Guérout A-M, Mazel D, Beloin C, Ghigo J-M. 2010. Silent mischief: bacteriophage Mu insertions contaminate products of Escherichia coli random mutagenesis performed using suicidal transposon delivery plasmids mobilized by broad-host-range RP4 conjugative machinery. J Bacteriol 192 :6418–6427. doi:10.1128/JB.00621-10 20935093
73 Green MR, Sambrook J. 2012. Molecular cloning: a laboratory manual. 4th ed. Cold Spring Harbor Laboratory Press, Cold Spring Harbor, N.Y.
74 Marvig RL, Blokesch M. 2010. Natural transformation of Vibrio cholerae as a tool - Optimizing the procedure. BMC Microbiol 10 :155. doi:10.1186/1471-2180-10-155 20509862
75 De Souza Silva O, Blokesch M. 2010. Genetic manipulation of Vibrio cholerae by combining natural transformation with FLP recombination. Plasmid 64 :186–195. doi:10.1016/j.plasmid.2010.08.001 20709100
76 Blokesch M. 2012. TransFLP — A method to genetically modify Vibrio cholerae based on natural transformation and FLP-recombination. J Vis Exp 68 :3761. doi:10.3791/3761
77 Meibom KL, Li XB, Nielsen AT, Wu C-Y, Roseman S, Schoolnik GK. 2004. The Vibrio cholerae chitin utilization program. Proc Natl Acad Sci U S A 101 :2524–2529. doi:10.1073/pnas.0308707101 14983042
78 Bao Y, Lies DP, Fu H, Roberts GP. 1991. An improved Tn7-based system for the single-copy insertion of cloned genes into chromosomes of gram-negative bacteria. Gene 109 :167–168. doi:10.1016/0378-1119(91)90604-a 1661697
79 Aoki SK, Lillacci G, Gupta A, Baumschlager A, Schweingruber D, Khammash M. 2019. A universal biomolecular integral feedback controller for robust perfect adaptation. Nature 570 :533–537. doi:10.1038/s41586-019-1321-1 31217585
80 Gilchrist CLM, Chooi Y-H. 2021. clinker & clustermap.js: automatic generation of gene cluster comparison figures. Bioinformatics 37 :2473–2475. doi:10.1093/bioinformatics/btab007 33459763
81 Tatusova T, DiCuccio M, Badretdin A, Chetvernin V, Nawrocki EP, Zaslavsky L, Lomsadze A, Pruitt KD, Borodovsky M, Ostell J. 2016. NCBI prokaryotic genome annotation pipeline. Nucleic Acids Res 44 :6614–6624. doi:10.1093/nar/gkw569 27342282
82 Mirdita M, Schütze K, Moriwaki Y, Heo L, Ovchinnikov S, Steinegger M. 2022. ColabFold: making protein folding accessible to all. Nat Methods 19 :679–682. doi:10.1038/s41592-022-01488-1 35637307
83 Holm L. 2020. Using dali for protein structure comparison, p 29–42. In Gáspári Z (ed), Structural bioinformatics. Springer US, New York, NY.
84 Holm L. 2022. Dali server: structural unification of protein families. Nucleic Acids Res 50 :W210–W215. doi:10.1093/nar/gkac387 35610055
85 Abby SS, Néron B, Ménager H, Touchon M, Rocha EPC. 2014. MacSyFinder: a program to mine genomes for molecular systems with an application to CRISPR-Cas systems. PLoS One 9 :e110726. doi:10.1371/journal.pone.0110726 25330359
86 Katoh K, Standley DM. 2013. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol 30 :772–780. doi:10.1093/molbev/mst010 23329690
87 Eddy SR. 2011. Accelerated profile HMM searches. PLoS Comput Biol 7 :e1002195. doi:10.1371/journal.pcbi.1002195 22039361
