Increasing prevalence of bacteriocin carriage in a six-year hospital cohort of E. faecium
=========================================================================================

* Andrea Garretto
* Suzanne Dawid
* Robert Woods

## Abstract

Vancomycin resistant enterococci (VRE) are important pathogens in hospitalized patients, however, the factors involved in VRE colonization of hospitalized patients are not well characterized. Bacteriocins provide a competitive advantage to enterococci in experimental models of colonization, but little is known about bacteriocin content in samples derived from humans and even less is known about their dynamics in the clinical setting. To identify bacteriocins which may be relevant in the transmission of VRE, we present a systematic analysis of bacteriocin content in the genomes of 2,428 patient derived *E. faecium* isolates collected over a six-year period from a single hospital system. We used computational methods to broadly search for bacteriocin structural genes and a functional assay to look for phenotypes consistent with bacteriocin expression. We identified homology to 15 different bacteriocins with two having high presence in this clinical cohort. Bacteriocin 43 (bac43) was found in a total of 58% of isolates, increasing from 8% to 91% presence over the six-year collection period. There was little genetic variation in the bac43 structural or immunity genes across isolates. The enterocin A structural gene was found in 98% of isolates but only 0.3% of isolates had an intact enterocin A gene cluster and displayed a bacteriocin producing phenotype. This study presents a wide survey of bacteriocins from hospital isolates and identified bac43 as highly conserved, increasing in prevalence, and phenotypically functional. This makes bac43 an interesting target for future investigation for a potential role in *E. faecium* transmission.

**Importance** While enterococci are a normal inhabitant of the human gut, vancomycin-resistant *E. faecalis* and *E. faecium* are urgent public health threats responsible for hospital associated infections. Bacteriocins are ribosomally synthesized antimicrobial proteins and are commonly used by bacteria to provide a competitive advantage in polymicrobial environments. Bacteriocins have the potential be used by *E. faecium* to invade and dominate the human gut leading to a greater propensity for transmission. In this work, we explore bacteriocin content in a defined clinically derived population of *E. faecium* using both genetic and phenotypic studies. We show that one highly active bacteriocin is increasing in prevalence over time and demonstrates great potential relevance to *E. faecium* transmission.

## Introduction

Vancomycin-resistant enterococci (VRE) are healthcare-associated pathogens that are transmitted on surfaces and between patients and healthcare workers, resulting in an estimated 54,500 infections and 5,400 deaths each year(1). Most clinical infections with VRE are caused by either *E. faecalis* of *E. faecium*, with vancomycin-resistant *E. faecium* being distinguished from other common nosocomial pathogens by its increased propensity to transmit in hospital settings(2). VRE is endemic in many health care settings resulting in opportunities for patients to acquire VRE during their hospital stay(3, 4). When VRE transmits to a previously uncolonized patient it can, in some patients, dominate in the gastrointestinal community before causing infections(5). The mechanisms that promote this colonization and gut domination by VRE, however, are unclear prompting a need to better understand the underlying biology.

Bacteriocins are ribosomally synthesized antimicrobial proteins and are commonly used by bacteria to provide a competitive advantage in polymicrobial environments(6, 7) . Bacteriocins are produced by many lactic acid bacteria, including enterococci(8) and their production has been suggested to play a role in transmission and intestinal domination, but this role is almost entirely unexplored in the clinical setting. Within the gut, enterococci compete extensively among themselves and with other species(9). In a mouse model of enterococcal gut colonization, it was demonstrated that *E. faecalis* harboring bacteriocin 21 (bac-21) was able to outcompete *E. faecalis* strains lacking this bacteriocin(10). The bacteriocin-producing *E. faecalis* were also able to maintain long-term gut colonization in contrast to the non-bacteriocin producing strain. Gut colonization is a necessary first step in gut dominance(5, 11), and patients colonized with high densities of VRE have a greater amount of shedding into the environment leading to more opportunities for transmission(12). Therefore, bacteriocin production has the potential to increase gut dominance and lead to greater transmission yet an understanding of the specific bacteriocins that could play this role is lacking.

Here, we present genomic analysis of a cohort of 2,428 *E. faecium* blood and perirectal isolates collected over a six-year period from a single hospital system. This sizeable cohort allows for an expansive search for a variety of potential bacteriocin loci, an analysis of how bacteriocin carriage changes over time, and the ability to shed light on how these bacteriocins could promote *E. faecium* transmission in this important setting. We identified many previously described bacteriocins in the patient derived isolates, most of which were found in limited numbers or did not have consistent phenotypic evidence of activity. However, bac43, a plasmid encoded bacteriocin, was found in 58% of isolates and became more prevalent over the study period. Isolates encoding bac43 had phenotypic evidence of bacterial inhibition and were found throughout the phylogenetic tree with evidence of both horizontal and vertical acquisition events.

## Methods

### Ethics

This study was approved by the University of Michigan Institutional Review Board (ID no. HUM00102282), which determined that informed consent was not required as all data utilized were collected for patient treatment purposes.

### Sample collection

*E. faecium* samples were collected from the Michigan Medicine hospital system between 2016 and 2021. Perirectal (PR) swabs were collected as part of the VRE surveillance screening protocol. Patients were screened upon admission to any of the hospital’s nine high risk units and then weekly during their stay in those units. Blood (BL) samples were collected from all patients with *E. faecium* bloodstream infection.

Per hospital protocol, perirectal (PR) swabs were streaked on VRE select plates (Bio-Rad) by the clinical microbiology laboratory. Each positive VRE select plate was transferred to our research laboratory where an isolated pink colony was selected at random, re-streaked twice on Brain Heart Infusion (BHI) agar (Becton, Dickinson and Company, 211065), and stored in BHI broth (Becton, Dickinson and Company, 237500) supplemented with 20% glycerol at -80°C.

Blood (BL) samples were obtained from the clinical microbiology laboratory having been streaked on either BBL TSA II 5% SB or BBL chocolate II agar plates. Random isolated colonies were selected and re-streaked twice on BHI plates with overnight growth at 35°C. If there were multiple clearly distinct morphologies, one of each was selected. Freezer stocks were created from a single colony stored in BHI broth supplemented with 20% glycerol at -80°C.

### Sequencing and assembly

DNA extraction from clones was performed using the Omega Bio-tek Mag-BIND Bacterial DNA 96 Kit (M2350). Illumina library preparation was performed using Collibri PCR-free ES DNA Library Prep Kit for Illumina - Invitrogen (CAT# A39123196) and sequencing was performed using Illumina Nova Seq or HiSeq platform. Reads were processed using a snakemake pipeline for read trimming using Trimmomatic (v0.39)(13), *de novo* assembly using Unicycler (v0.4.8)(14), and annotation using Prokka (v1.14.6)(15). Multilocus sequence types were identified using the software mlst (v2.22.0)(16).

### Hybrid assemblies

Select genomes were assembled using a combination of long and short read data. Long read sequencing was performed using a MinION with a Flongle flow cell from Oxford Nanopore Technologies. Libraries were prepared using the Rapid Barcoding Sequencing Kit (SQK-RBK004) and the Flongle flow cells were primed using the Custom Kit (CUST-KIT). Sequencing data were obtained using the MinKNOW software from Oxford Nanopore Technologies. Adapters were trimmed from long-read sequences using Porechop(v0.2.3)(17), *de novo* assemblies were generated using Unicycler(v0.4.8)(14), and assemblies were annotated using Prokka(v1.14.6)(15).

### Identification of bacteriocins in genome assemblies

A list of bacteriocin protein sequences was obtained from BAGEL4(18), including bacteriocins from classes 1-3. Sequences identical at the amino acid level were removed using the sRNAtoolbox(19). Bacteriocin sequences within our assembled genomes were identified by creating custom BLAST databases of each assembly(20). tblastn was then used to query bacteriocin sequences against each genome assembly database. Output files were filtered using R(v4.3.1). If multiple hits were made to overlapping regions of a genome, the hit with the greatest bitscore was kept. Hits with an e-value greater than 1e-08 and a percent identity less than 60% were removed. Hits from all assemblies were clustered using CD-HIT(21) based on 90% sequence identity and 75% minimal alignment coverage for the longer sequence.

Searches were also performed using nucleotide sequences of other genes involved in bacteriocin production and immunity for enterocin A, bac43, enterocin NKR-5-3B (Supplemental Table 1). These genes were used as the query against each of the custom BLAST databases for each genome assembly.

### Variant calling and phylogenetic analysis

Variants were identified using a custom snakemake pipeline. Trimmed reads were mapped to an internal reference strain *E. faecium* BL00198-1 and variants were called using snippy(22). A core genome was produced using snippy-core; samples with less than 250,000 base pairs aligned were excluded. The core alignment was then cleaned using snippy-clean_full_aln, and a phylogenetic tree was produced using Gubbins(23). Gubbins was run for six iterations where the first tree built using ‘rapidnj’ and ‘JC’ model. The remaining five iterations were built using ‘raxmlng’ and the ‘GTR’ model. Samples with less than 75% homology to the reference strain were excluded. The tree figure was rendered and annotated using ggtree in R(v4.3.1). MLSTs with less than 20 isolates associated with them are classified as “Other”. Bacteriocin sequences with less than five isolates associated with them are classified as “Other”.

### Variation analysis of enterocin A gene cluster

Variants in the enterocin A gene cluster and the bacteriocin 43 encoded plasmid were identified using an internal snakemake pipeline. For enterocin A, trimmed reads were mapped to internal reference with the entire enterocin A gene cluster (BL02040-1). Variants were called using FreeBayes(24) with the following parameters: minimum base quality greater than or equal to 10 and minimum alternate fraction less than or equal to 0.5. Variants were annotated using SnpEff(25). Insertion elements were identified and annotated using panISa(26) and ISFinder(27). Areas of low read coverage indicating gene deletions were identified using SAMtools(28) depth and filtering to positions with at least 10% of the average coverage across the whole gene locus. Variants, including insertion elements, identified in deleted regions were excluded using custom R scripts (v4.3.1). Read pileups were visualized using igv-reports(29). Variants and low coverage regions were manually inspected using these read maps; mapped regions were identified in isolate genome assemblies to confirm presence or absence of variants and low coverage regions.

### Variation analysis of bac43 encoding plasmid pDT1

Trimmed reads were mapped to an internal reference isolate, PR46485-1-C1, which encodes a plasmid identical to the previously published pDT1 sequence(30). Variant calling, annotation, identification of deletions, and insertion element identification was repeated for pDT1 with the same methodology as for enterocin A. Additional analyses for pDT1 included coverage calculations across the plasmid and bac43 specific genes using SAMtools(28) coverage. Isolates with less than 45% mean depth or less than 55% coverage across the plasmid were excluded from further analysis. Variant positions were accepted if the RO (reference observations) to DP (total reads mapped) ratio was less than 10%, and if the RO was less than 30 reads. Percent identity was then calculated as the number of non-variant positions divided by the total number of bases called. Isolates were finally filtered by those with at least 99.5% identity to the aligned portions of the pDT1 plasmid or if the BLAST analysis identified the full structural gene of bac43.

### Visualizing genetic associate among variants

Genetic associations between variants in the enterocin A gene cluster and pDT1 were visualized using PHYLOViZ(31) (v2.0) using the goeBURST Full MST algorithm (goeBURST distance). Subcluster annotations were made based on shared variants and/or variations within the same genomic regions. Gene homology maps were created using Easyfig(32) (v2.2.2) with a minimum identity threshold of 60%.

### Bacteriocin activity assay

Putative bacteriocin activity was assessed using a modified soft-agar assay for bacteriocin production and immunity(30). Isolates were streaked from freezer stock onto BHI agar plates and incubated at 35°C overnight. Single colonies were resuspended in 50ul saline and poked into a soft agar lawn containing a bacteriocin sensitive strain. If results were discordant with functional predictions, a higher inoculum was tested by using a toothpick to poke a single colony directly into the lawn. After incubation at 35°C overnight, putative killing effect was determined by the presence of a halo of bacterial clearing around the poke into the bacteriocin sensitive lawn. Immunity to enterocin A was determined poking an enterocin A producer (BL02040-1) into a soft agar bacterial lawn of the strain being tested and looking for the absence of a halo of bacterial clearing.

Immunity to bac43 was determined by the ability to grow in the presence of supernatant derived from a bac43 producing strain (PR46485-1-C1). Bac43 supernatant was obtained by incubation of PR46485-1-C1 in 45mL BHI at 35°C, 130rpm, for 24 hours. The culture was centrifuged for 15 minutes at 3220rcf and the supernatant was filter sterilized twice (0.22um filter). Colony resuspensions were plated on BHI agar supplemented with 10% bac43 supernatant. Plates were incubated at 35°C overnight. Immunity to bac43 was determined by bacterial growth on BHI plates supplemented with 10% bac43 lysate.

Isolates were chosen for the bacteriocin production and resistance assays as follows. When sequences were predicted to be nonfunctional (i.e. deletions or insertion elements present in bacteriocin genes) only one isolate with that mutation was tested. If the isolate with the presumed non-functional bacteriocin was found to be negative on the bacteriocin production and immunity assays, other isolates carrying the mutation were not tested. For sequence isolates from multiple MLSTs, one isolate was tested from each of the MLSTs present in at least 20 of all isolates.

### Data availability

All scripts and pipelines are available on GitHub ([https://github.com/woodslab/bacteriocin\_search](https://github.com/woodslab/bacteriocin\_search)). The whole genome sequencing project for our reference *E. faecium* isolates have been deposited at DDBJ/ENA/GenBank: [BL00198](http://medrxiv.org/lookup/external-ref?link\_type=GEN&access\_num=BL00198&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F18%2F2024.07.17.24310592.atom)- 1 (accession JBEFKN000000000), [BL02040](http://medrxiv.org/lookup/external-ref?link\_type=GEN&access_num=BL02040&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F18%2F2024.07.17.24310592.atom)- 1 (accession JBEFKM000000000), [PR46485](http://medrxiv.org/lookup/external-ref?link_type=GEN&access_num=PR46485&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F18%2F2024.07.17.24310592.atom)-1-C1 (JBEFKL000000000). The versions described in this paper are JBEFKN010000000, JBEFKM010000000, and JBEFKL010000000, respectively.

## Results

### Samples included

The cohort for this study included the first *E. faecium* isolate from each patient with a VRE *faecium* positive rectal surveillance swab between 2016 and 2021. This consisted of 2300 patients (Table 1). Additionally, 278 patients were found to have *E. faecium* growth in blood cultures. The first isolate from surveillance and blood culture per patient was sequenced. One-hundred and twenty-eight patients had both a perirectal (PR) swab and a blood (BL) sample taken during their stay in the hospital resulting in a total of 2428 *E. faecium* isolates.

View this table:
[Table 1:](http://medrxiv.org/content/early/2024/07/18/2024.07.17.24310592/T1)

Table 1: Sample summary.

### Bacteriocins identified in the clinical cohort

BLAST hits to bacteriocin structural genes in our cohort were clustered based on sequence identity. There were 17 bacteriocin clusters identified with homology to 15 different bacteriocins (Figure 1, supplemental table 3). Most bacteriocin sequence clusters were rare, being found in less than 3% of isolates. The majority of these low frequency clusters have previously been characterized from *E. faecium*(33–39), but enterocin SE-K4(40) and bac31(41) were characterized from *E. faecalis* and pneumolancidin(42) was characterized from *S. pneumoniae*. We chose the four most common cluster hits to investigate further as their higher prevalence suggested a potentially relevant role in hospital transmission. This includes the enterocin A cluster (98% of isolates), two cluster hits to bac43 (57.9% of isolates and 40.7% of isolates, respectively), and enterocin NKR-5-3B cluster (19.2% of isolates)(Figure 2A-D).

![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/18/2024.07.17.24310592/F1.medium.gif)

[Figure 1:](http://medrxiv.org/content/early/2024/07/18/2024.07.17.24310592/F1)

Figure 1: Frequency of bacteriocin gene sequence hits.
Bacteriocin sequences were classified into 17 different clusters with homology to 15 different bacteriocins. Four bacteriocin clusters were identified as present in over 19%: enterocin A cluster, bac43 full cluster, bac43 partial cluster, and enterocin NKR-5-3B cluster.

![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/18/2024.07.17.24310592/F2.medium.gif)

[Figure 2:](http://medrxiv.org/content/early/2024/07/18/2024.07.17.24310592/F2)

Figure 2: Gene maps and homology to previously identified bacteriocins.
Genes related to bacteriocin production are blue, hypothetical proteins are yellow, and insertion elements are red. A. Internal reference of the enterocin A gene cluster. B. Internal reference of the bacteriocin 43 encoding plasmid pDT1. C. Homology of the bacteriocin 43 genes to the bac43 partial cluster BLAST hit. D. Homology of the enterocin NKR-5-3B genes to the enterocin NKR-5-3B cluster BLAST hit.

The enterocin NKR-5-3B cluster shares 87.3% amino acid homology to the structural gene *enkB* of enterocin NKR-5-3B and was identified in 19.2% of isolates. This class IIc circular bacteriocin was previously characterized on a 9kb chromosomal region of *E. faecium*(43). The previously described functional operon includes five genes associated with enterocin NKR-5-3B: *enkB* the bacteriocin structural prepeptide, *enkB1* and *enkB3* for secretion and self-immunity, *enkB2* for maturation, and *enkB4* for self-immunity. None of the clinical isolates contained the entire previously described 9kb region. Rather we found homology to *enkB* and approximately 200 bp of *enkB1* (Figure 2D). Given the lack of predicted function, this cluster was not further investigated.

### Analysis of the enterocin A clusters

The enterocin A structural gene, *entA,* was identified in 98.6% of isolates(44, 45). Enterocin A is a class IIa bacteriocin previously characterized in *E. faecium* strain CTC492. In that strain, the enterocin A gene cluster is encoded on a 10.5kb region of the bacterial chromosome(45) (Figure 2A), and the locus includes genes for structural protein (*entA)*, self-immunity (*entI)*, the peptide pheromone (*entF)*, the two-component regulatory system (*entK* and *entR)*, and the peptide transporter (*entT and entD)*.

We identified variant positions, insertion elements, and deletions by comparing isolates to an internal reference. This reference sequence has a missense variant in the *entI* gene (position 246717 in BL02040-1, G to A) compared to the previously published sequence(45). We identified 87 unique sequences of the enterocin A gene cluster in our collection, 61 of which were found in only a single isolate. These sequences were assigned to five subclusters labeled A through E (Figure 3). Subcluster A is centered around sequence 1, which is found in 84% of all isolates. Sequence 1 isolates have a large deletion spanning *entFKRTD*. Of the 41 other sequences in subcluster A, 29 share this deletion. Four sequences in subcluster A (4, 6, 9, and 12) have all the genes required for enterocin A function intact. In subcluster B, four sequences (3, 49, 73, and 85) also have all the enterocin A genes. The other sequence in subcluster B (35) has an insertion element in histidine protein kinase *entK* likely interrupting function. Subcluster C sequences have an insertion element that is likely inactivating the regulator gene *entR.* Subcluster D shares the insertion element in *entR* with subcluster C plus an additional insertion element the immunity gene *entI* also likely interrupting function. Subcluster E has three sequences which all have the complete enterocin A gene cluster but with approximately 30 synonymous and missense mutations.

![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/18/2024.07.17.24310592/F3.medium.gif)

[Figure 3:](http://medrxiv.org/content/early/2024/07/18/2024.07.17.24310592/F3)

Figure 3: Enterocin A sequence clusters.
Isolates were assessed for SNPs, insertion elements, and gene deletions within the enterocin A gene locus and are clustered based on matching patterns of variations (sequences). Each node is a unique sequence of the enterocin A locus and branch numbers represent the number of variations between the sequences. Node size increases based on the number of isolates belonging to that sequence. The results of phenotypic assays on tested strains are noted on each sequence node. Inhibitory sequences (blue and red) have both the full enterocin A locus and can kill a bacteriocin sensitive strain, inhibitory/other sequences (red) also encode other bacteriocins. Not inhibitory/disrupted sequences (yellow) either lacked inhibition in phenotypic assays and/or had disruptions in the enterocin A locus that would predict nonfunction. Sequence 3 isolates are marked as non-inhibitory due to mixed phenotypic results.

### Analysis of the bac43 gene clusters

Bacteriocin 43 (bac43) is a class IIa *sec*-dependent bacteriocin previously isolated from clinical vancomycin resistant *E. faecium* isolates(30). The 6.1kb plasmid pDT1 encodes the bac43 structural gene *bacA*, the gene for self-immunity *bacB*, plasmid mobilization genes *mobC* and *mobA*, and plasmid replication genes *repA* and *repB*(30)(Figure 2B).

Bac43 was identified as the closest bacteriocin hit in two common clusters. One of these, referred to as ‘bac43 partial’ cluster, was found in 40.7% of isolates and aligned only to part of the *bacA* coding region, has no homology to the immunity gene, and the start codon found on *bacA* was absent. Thus, we do not expect it to be functional. The aligned bac43 partial sequence displayed only 87.5% amino acid homology to previously published *bacA*, and were not found on contigs with other pDT1 genes (Figure 2C).

The second cluster, ‘bac43 full,’ was identified in 57.9% of isolates and displayed 99.5–100% amino acid homology to the published *bacA*. Sequences falling into the bac43 full cluster also have homology to the immunity gene *bacB* and are often found on contigs consistent with the pDT1 plasmid.

Variants in the pDT1 plasmid were identified by comparing isolates to an internal reference genome carrying a plasmid identical to the one previously published(30). Four isolates discordant with phenotypic predictions were assayed using PCR for the presence of the *bacA* (Supplemental Table 2). PCR revealed an absence of this gene, and these isolates were removed from further analysis. We identified 106 unique sequences of the bac43 encoding plasmid pDT1, 77 of which were found in only one isolate. These sequences are organized into nine subclusters (Figure 4). Subcluster A is centered around sequence 3 which is found in 39.7% of isolates and is identical to the published pDT1 sequence(30). Most of the variation within this subcluster is outside of *bacA* and *bacB*, but four sequences had a mutation in one of these genes: a synonymous variant in *bacA* (sequence 82), an insertion element in *bacB* (sequence 14), and two with insertion elements in different positions in *bacA* (sequences 57 and 84). Additionally, sequence 65 (subcluster H) has a missense variant in *bacB*. Sequences 9 and 73 (subcluster D) have lost the bac43 genes but retained portions of the pDT1 plasmid. The remainder of the sequences are characterized by insertion elements and/or alterations to the plasmid backbone.

![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/18/2024.07.17.24310592/F4.medium.gif)

[Figure 4:](http://medrxiv.org/content/early/2024/07/18/2024.07.17.24310592/F4)

Figure 4: Clustered variations within the bac43 encoding plasmid pDT1.
Isolates were assessed for variants, insertion elements, and gene deletions within the bac43 encoding plasmid pDT1 and are clustered based on matching patterns of variations (sequences). Each node is a unique sequence of pDT1 and branch numbers represent the number of variations between the sequences. Node color represents phenotypic testing results. Node size is determined by number of isolates belonging to that sequence. Inhibitory sequences (blue) can kill a bacteriocin sensitive strain.

### Phenotypic evaluation of bacteriocin containing isolates

Phenotypic assays were performed to determine whether strains produce a phenotype of (1) killing of a strain with no bacteriocin sequences and (2) immunity to killing by strains that encode a known bacteriocin producer (for enterocin A) or a crude extract (for bac43).

For strains encoding the enterocin NKR-5-3B-like sequence, we tested seven isolates from different MLST backgrounds. As expected, none of the isolates demonstrated inhibition against a bacteriocin sensitive strain.

For enterocin A, eleven sequences were classified as potentially functional based on the absence of deletions or insertion elements into the gene cluster. Isolates encoding six of these sequences showed both consistent killing of the sensitive strain and resistance to killing by an enterocin A producer consistent with a fully functional enterocin A locus (sequences 6, 9, 12, 26, 73, and 85)(Figure 5A); four of these isolates (sequences 6, 9, 26, and 85), however, also encode other bacteriocins such as enterocin B, enterocin P, enterocin X𝛼, and enterocin SE-K4 which could not be ruled out as the source of inhibition. The five sequence 3 isolates that were tested gave mixed results. Four of the five tested isolates encoding sequence 3 demonstrated an inhibitory and resistant phenotype that could be consistent with a functional enterocin A locus, but the contribution of other bacteriocins could not be excluded. The only sequence 3 isolate that encodes no other bacteriocins had no inhibitory activity but did demonstrate resistance to killing by enterocin A. This pattern of being resistant but noninhibitory was also demonstrated by sequence 49. As neither of these isolates encode other bacteriocins the enterocin A resistance could stem from *entI* expression or chromosomal mutations. The remaining three sequences did not have an inhibitory phenotype. The sequence 4 isolate showed resistance without inhibition but unlike the preceding isolates, the tested isolate also encodes other bacteriocins. Sequence 15 and 19 lack both inhibitory activity and enterocin A resistance; both have numerous mutations that could have resulted in non-function (Figure 5B).

![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/18/2024.07.17.24310592/F5.medium.gif)

[Figure 5:](http://medrxiv.org/content/early/2024/07/18/2024.07.17.24310592/F5)

Figure 5: Enterocin A phenotypic analysis.
A. Overview of enterocin A isolates tested, functional expectations, and experimental results. B. Gene maps aligned with the enterocin A reference locus for sequences that were expected to be functional (i.e. no deletions or insertion elements within the bacteriocin genes).

Sixteen sequences were predicted to be non-functional due to deletions of or insertions into genes essential for enterocin A production. All of the tested isolates had a phenotype that was neither resistant nor inhibitory in isolates that did not encode other bacteriocins. Seven isolates did encode other bacteriocins and these had both inhibitory activity and grew in the presence of the enterocin A producer, suggesting the other bacteriocin loci, rather than the enterocin A locus, were the cause of the assay results.

Phenotypic testing was performed on 132 isolates representing 106 sequences with pDT1 variants (Figure 6A). Ninety-nine of these had 100% conserved sequences of *bacA* and *bacB* but other mutations in the pDT1 plasmid; all 124 isolates tested with these 99 sequences were able to kill the bac43 sensitive strain and were resistant to killing by a bac43 producer. Among the five sequences with mutations in the bac43 genes, three retained both an inhibitory phenotype and bac43 resistance: sequence 65 (missense variant in *bacB*), 82 (synonymous variant in *bacA*), and 14 (insertion element in *bacB*). However, no inhibition or resistance to bac43 were observed by isolates with an insertion element in *bacA* (sequences 84 and 57). Because the variant analysis aligned sequences to the full pDT1 plasmid, two sequences (9 and 73) were identified that lost the bac43 genes but retained portions of the plasmid (Figure 6B); these isolates did not demonstrate inhibition or resistance to bac43.

![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/18/2024.07.17.24310592/F6.medium.gif)

[Figure 6:](http://medrxiv.org/content/early/2024/07/18/2024.07.17.24310592/F6)

Figure 6: Bac43 phenotypic analysis.
A. Overview of bac43 isolates tested, functional expectations, and experimental results. B. Gene maps aligned with the bac43 encoding plasmid pDT1 for sequences with variants in the bac43 genes *bacA* and *bacB*.

### Phylogenetic analysis

A phylogenetic tree of the isolates was created to examine the evolutionary relationships between the bac43 and enterocin A sequences. The enterocin A sequences cluster together supporting a predominantly vertical mode of transmission (Figure 7A). Sequences are not restricted to MLSTs: for example, sequence 1 with a deletion of *entFKRTD* is found in multiple MLST backgrounds suggesting this wide-spread deletion may have occurred before MLST differentiation. Five of the seven phenotypically functional isolates encoding the full enterocin A locus (six sequence variations plus the enterocin A reference) passed the homology thresholds for inclusion on the tree. Four of these isolates cluster together while one is found in another cluster surrounded by nonfunctional sequence 1 isolates consistent with horizontal gene transfer of the intact gene locus.

![Figure 7:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/18/2024.07.17.24310592/F7.medium.gif)

[Figure 7:](http://medrxiv.org/content/early/2024/07/18/2024.07.17.24310592/F7)

Figure 7: Phylogenetic mapping of bacteriocin genotype and phenotype.
Core genome phylogeny with recombination removed by Gubbins. Outermost heatmap represents the bacteriocin sequence. Sequences with less than five isolates are classified as “Other.” Innermost heatmap is the isolate’s MLST. MLSTs with less than 20 isolates are classified as “Other.” A. Enterocin A sequences. Enterocin A isolates with the full enterocin A locus and an inhibitory phenotype are each indicated by an “*” B. Bac43 sequences.

For bac43, similar sequences tend to cluster together consistent with prominence of vertical transmission (Figure 7B). Sequence 3 is mostly found within MLST 117 and sequence 2 is associated with MLST 80. However, occasionally the same sequence can be found in distant clades supporting horizontal transmission events. Functional bac43 is found in 18 of the 59 MLSTs (seven of nine MLSTs found in more than 20 isolates).

### Temporal trends

The frequency of bac43 increased over the study period, rising from 8% in the first quarter of 2016 to 91% presence by the end of 2021 (Figure 8A). At the beginning of the collection period, no single MLST was present in more than 20% of isolates. However, MLST 117 subsequently rose to become the dominant sequence type with MLST 80 steadily increasing (Figure 8B). Isolates encoding bac43 are highly associated with MLSTs 117 and 80 accounting for 66.6% and 18.5% of bac43 isolates, respectively. By contrast, the six isolates that demonstrated functional enterocin A activity were dispersed over only the first five years of the study and had no clear association with MLSTs.

![Figure 8:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/18/2024.07.17.24310592/F8.medium.gif)

[Figure 8:](http://medrxiv.org/content/early/2024/07/18/2024.07.17.24310592/F8)

Figure 8: Temporal trends in bacteriocin variants and MLST.
Presence of (A) functional bac43 and enterocin A sequences by quarter year and (B) MLSTs over the six-year study period, over the study period by quarter year. MLSTs with less than 20 isolates are classified as “Other.” Data includes first isolate per patient only.

### Prevalence by isolate source

Phenotypically functional enterocin A was found in 0.09% of perirectal surveillance swab (PR) isolates and 1.9% of blood (BL) isolates with significantly more enterocin A isolates in the BL cohort (p-value 0.00028)(Table 2). Conversely, bac43 is present in significantly more PR isolates (p-value 0.0086). However, our collection protocol only identified vancomycin resistant PR isolates, while BL includes both vancomycin resistant and vancomycin sensitive isolates. When we restrict the comparison to isolates that were vancomycin resistant, there is no significant differences in bacteriocin presence between the PR and BL isolates (Table 3).

View this table:
[Table 2:](http://medrxiv.org/content/early/2024/07/18/2024.07.17.24310592/T2)

Table 2: Bacteriocin presence between perirectal (PR) and blood (BL) isolates.

View this table:
[Table 3:](http://medrxiv.org/content/early/2024/07/18/2024.07.17.24310592/T3)

Table 3: Bacteriocin presence between vancomycin resistant perirectal (VRE PR) and vancomycin resistant blood (VRE BL) isolates.

## Discussion

Enteroccocci occupy the intestinal tract of a broad array of host species and compete intensively with other inhabitants of this niche(5). Bacteriocins are commonly found in enterococci and believed to be effective weapons in this competition(6), but little is known about the role they play in clinical settings. Here we investigate the epidemiology and evolution of bacteriocins in *E. faecium* in a clinical setting by analyzing bacteriocin genetic content and phenotype of patient-derived isolates collected over a six-year period. Bacteriocin 43 (bac43) stands out as being the most common functional bacteriocin and shows a dramatic increase in frequency over time.

Bac43 was first isolated from a collection of VRE isolates from this same hospital in the late 1990s(30). From those isolates collected over five years, bac43 was identified in 3.3%. From our current study, we find the bac43 genes and its encoding plasmid to be mostly conserved with 40% of our isolates encoding a plasmid identical to the initially described plasmid(30). When considering both sets of data, we see that bac43 has been present within this hospital for over two decades with a recent rise in frequency. However, this observational genomic survey cannot identify a causal role for bac43 in this increasing frequency of over time. It is tempting to suggest bac43 could impact colonization and transmission given the strong phenotype it produces in the lab. Indeed, prior work using a mouse model of intestinal colonization showed that *E. faecalis* harboring a plasmid encoded bacteriocin bac-21 were able to eliminate bac-21 sensitive strains over the course of their experiment(10). That study also demonstrated bac-21 plasmid was actively conjugated within the GI tract(10). If pDT1 acts similarly in being advantageous to *E. faecium* for invasion and domination of the human gut, this could explain how bac43 has come to dominate our clinical cohort. However, there are alternative explanations. The strong association between bac43 and the MLSTs that have become dominant in our cohort (117 and 80), suggest these could be rising due to hitchhiking with some other features of these genomes. It is also possible that bac43 could be acting as a selfish plasmid, mostly serving to maintain its own existence with little benefit to the host bacterium. Finally, we have focused on the intraspecific competition, but it is also possible that bac43 could impact interspecific competition, especially among closely related species(30).

Prior studies have identified bac43 from clinical isolates of vancomycin resistant and sensitive enterococci (30, 46, 47) and demonstrated a correlation between bac43 and vancomycin resistant *E. faecium*(47). Consistently, we also find that bac43 is associated with vancomycin resistance. After correcting for this association, we see a similar presence of bac43 in strains from both colonization and infection.

This study also adds clarity to the contribution of enterocin A through genomic and phenotypic analyses. While the structural gene for enterocin A is found in 98% of *E. faecium* isolates within this cohort, the vast majority of these appear to be non-functional. Seven of the 11 sequences with a full enterocin A locus (0.3% of isolates) demonstrated an inhibitory phenotype. We also identified three isolates with the full enterocin A locus that did not have an inhibitory phenotype but did appear to be resistant to killing by enterocin A. While these isolates still retain the enterocin A immunity gene, this resistance could also reflect mutations in the bacteriocin target site, mannose phosphotransferase system(48).

Caution should be taken when interpreting the phenotypic results, as the inhibition and resistance assays we performed could be confounded by the presence of other bacteriocins. When present, other bacteriocins may result in killing of the sensitive strain used in the assay. Additionally, apparent resistance to the enterocin A producer could be due to other bacteriocin loci conferring resistance to enterocin A, or to suppression of the growth of the enterocin A producer. The presence of other bacteriocins may explain why several isolates that had only portions of the enterocin A locus and were not predicted to be functional had both inhibition and enterocin A immunity. Future work to elucidate the molecular function of the diverse bacteriocin loci described in this paper are warranted.

This study represents a systematic analyses of bacteriocins from a clinical cohort and identifies bacteriocin 43 (bac43) as being mostly conserved and increasing in prevalence in clinical *E. faecium* isolates over a six-year collection period. As bacteriocins have been shown to be effective weapons in bacterial competition, bac43 is an important target for future investigation of the role of bacteriocins in the transmission of vancomycin resistant *E. faecium*.

## Supplemental data

View this table:
[Supplemental Table 1:](http://medrxiv.org/content/early/2024/07/18/2024.07.17.24310592/T4)

Supplemental Table 1: Genes for enterocin A, bac43, and enterocin NKR-5-3B

View this table:
[Supplemental Table 2:](http://medrxiv.org/content/early/2024/07/18/2024.07.17.24310592/T5)

Supplemental Table 2: PCR primers used for confirming bacteriocin gene presence

View this table:
[Supplemental Table 3:](http://medrxiv.org/content/early/2024/07/18/2024.07.17.24310592/T6)

Supplemental Table 3: All bacteriocin cluster hits and the number of isolates they are found in.

## Acknowledgements

AG was supported by the Molecular Mechanisms in Microbial Pathogenesis Training Program (T32 AI007528) and the University of Michigan Integrated Training in Microbial Systems (grant from Burroughs Wellcome Fund: Institutional Program Unifying Population and Laboratory Based Sciences). This work was supported by the National Institutes of Health under Grant No. R01 AI143852 to R.J.W.

*   Received July 17, 2024.
*   Revision received July 17, 2024.
*   Accepted July 18, 2024.


*   © 2024, Posted by Cold Spring Harbor Laboratory

This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available for use under a CC0 license

## References

1.  1.2019. VRE in Healthcare Settings | HAI | CDC. [https://www.cdc.gov/hai/organisms/vre/vre.html](https://www.cdc.gov/hai/organisms/vre/vre.html). Retrieved 14 October 2021.
    
    

2.  2.Sundermann AJ, Chen J, Kumar P, Ayres AM, Cho ST, Ezeonwuka C, Griffith MP, Miller JK, Mustapha MM, Pasculle AW, Saul MI, Shutt KA, Srinivasa V, Waggle K, Snyder DJ, Cooper VS, Van Tyne D, Snyder GM, Marsh JW, Dubrawski A, Roberts MS, Harrison LH. 2021. Whole-Genome Sequencing Surveillance and Machine Learning of the Electronic Health Record for Enhanced Healthcare Outbreak Detection. Clin Infect Dis Off Publ Infect Dis Soc Am 75:476–482.
    
    

3.  3.Khader K, Thomas A, Huskins WC, Stevens V, Keegan LT, Visnovsky L, Samore MH, for the Centers for Disease Control and Prevention (CDC) Prevention Epicenter Program and for the CDC Modeling Infectious Diseases in Healthcare Program. 2021. Effectiveness of Contact Precautions to Prevent Transmission of Methicillin-Resistant Staphylococcus aureus and Vancomycin-Resistant Enterococci in Intensive Care Units. Clin Infect Dis 72:S42–S49.
    
    

4.  4.Huskins WC, Huckabee CM, O’Grady NP, Murray P, Kopetskie H, Zimmer L, Walker ME, Sinkowitz-Cochran RL, Jernigan JA, Samore M, Wallace D, Goldmann DA. 2011. Intervention to Reduce Transmission of Resistant Bacteria in Intensive Care. N Engl J Med 364:1407–1418.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa1000373&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21488763&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F18%2F2024.07.17.24310592.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000289467200005&link_type=ISI) 

5.  5.Banla LI, Salzman NH, Kristich CJ. 2019. Colonization of the mammalian intestinal tract by enterococci. Curr Opin Microbiol 47:26–31.
    
    

6.  6.Heilbronner S, Krismer B, Brötz-Oesterhelt H, Peschel A. 2021. The microbiome-shaping roles of bacteriocins. Nat Rev Microbiol 19:726–739.
    
    

7.  7.Kerr B. 2007. The Ecological and Evolutionary Dynamics of Model Bacteriocin Communities, p. . In Riley, MA, Chavan, MA (eds.), Bacteriocins: ecology and evolution. Springer, Berlin ; New York.
    
    

8.  8.Montealegre MC, Singh KV, Murray BE. 2016. Gastrointestinal Tract Colonization Dynamics by Different Enterococcus faecium Clades. J Infect Dis 213:1914–1922.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiv597&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26671890&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F18%2F2024.07.17.24310592.atom) 

9.  9.Jett BD, Huycke MM, Gilmore MS. 1994. Virulence of enterococci. Clin Microbiol Rev 7:462–478.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiY21yIjtzOjU6InJlc2lkIjtzOjc6IjcvNC80NjIiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wNy8xOC8yMDI0LjA3LjE3LjI0MzEwNTkyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

10. 10.Kommineni S, Bretl DJ, Lam V, Chakraborty R, Hayward M, Simpson P, Cao Y, Bousounis P, Kristich CJ, Salzman NH. 2015. Bacteriocin production augments niche competition by enterococci in the mammalian gastrointestinal tract. Nature 526:719–722.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature15524&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26479034&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F18%2F2024.07.17.24310592.atom) 

11. 11.Ubeda C, Taur Y, Jenq RR, Equinda MJ, Son T, Samstein M, Viale A, Socci ND, van den Brink MRM, Kamboj M, Pamer EG. 2010. Vancomycin-resistant Enterococcus domination of intestinal microbiota is enabled by antibiotic treatment in mice and precedes bloodstream invasion in humans. J Clin Invest 120:4332–4341.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1172/JCI43918&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21099116&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F18%2F2024.07.17.24310592.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000284971400019&link_type=ISI) 

12. 12.Chanderraj R, Millar JA, Patel TS, Read AF, Washer L, Kaye KS, Woods RJ. 2019. Vancomycin-Resistant Enterococcus Acquisition in a Tertiary Care Hospital: Testing the Roles of Antibiotic Use, Proton Pump Inhibitor Use, and Colonization Pressure. Open Forum Infect Dis 6.
    
    

13. 13.Bolger AM, Lohse M, Usadel B. 2014. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinforma Oxf Engl 30:2114–2120.
    
    

14. 14.Wick RR, Judd LM, Gorrie CL, Holt KE. 2017. Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads. PLOS Comput Biol 13:e1005595.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pcbi.1005595&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28594827&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F18%2F2024.07.17.24310592.atom) 

15. 15.Seemann T. 2014. Prokka: rapid prokaryotic genome annotation. Bioinformatics 30:2068–2069.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btu153&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24642063&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F18%2F2024.07.17.24310592.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000339814300017&link_type=ISI) 

16. 16.Seemann T. mlst.
    
    

17. 17.Wick RR, Judd LM, Gorrie CL, Holt KE. 2017. Completing bacterial genome assemblies with multiplex MinION sequencing. Microb Genomics 3:e000132.
    
    

18. 18.van Heel AJ, de Jong A, Song C, Viel JH, Kok J, Kuipers OP. 2018. BAGEL4: a user-friendly web server to thoroughly mine RiPPs and bacteriocins. Nucleic Acids Res 46:W278–W281.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gky383&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29788290&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F18%2F2024.07.17.24310592.atom) 

19. 19.Aparicio-Puerta E, Lebrón R, Rueda A, Gómez-Martín C, Giannoukakos S, Jaspez D, Medina JM, Zubkovic A, Jurak I, Fromm B, Marchal JA, Oliver J, Hackenberg M. 2019. sRNAbench and sRNAtoolbox 2019: intuitive fast small RNA profiling and differential expression. Nucleic Acids Res 47:W530–W535.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkz415&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31114926&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F18%2F2024.07.17.24310592.atom) 

20. 20.Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, Madden TL. 2009. BLAST+: architecture and applications. BMC Bioinformatics 10:421.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2105-10-421&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20003500&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F18%2F2024.07.17.24310592.atom) 

21. 21.Huang Y, Niu B, Gao Y, Fu L, Li W. 2010. CD-HIT Suite: a web server for clustering and comparing biological sequences. Bioinforma Oxf Engl 26:680–682.
    
    

22. 22.Seemann T. 2015. snippy: fast bacterial variant calling from NGS reads. Perl.
    
    

23. 23.Croucher NJ, Page AJ, Connor TR, Delaney AJ, Keane JA, Bentley SD, Parkhill J, Harris SR. 2015. Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins. Nucleic Acids Res 43:e15.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gku1196&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25414349&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F18%2F2024.07.17.24310592.atom) 

24. 24.Garrison E, Marth G. 2012. Haplotype-based variant detection from short-read sequencing. arXiv:1207.3907. arXiv.
    
    

25. 25.Cingolani P, Platts A, Wang LL, Coon M, Nguyen T, Wang L, Land SJ, Lu X, Ruden DM. 2012. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff. Fly (Austin) 6:80–92.
    
    

26. 26.Treepong P, Guyeux C, Meunier A, Couchoud C, Hocquet D, Valot B. 2018. panISa: ab initio detection of insertion sequences in bacterial genomes from short read sequence data. Bioinforma Oxf Engl 34:3795– 3800.
    
    

27. 27.Siguier P, Perochon J, Lestrade L, Mahillon J, Chandler M. 2006. ISfinder: the reference centre for bacterial insertion sequences. Nucleic Acids Res 34:D32–36.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkj014&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16381877&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F18%2F2024.07.17.24310592.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000239307700009&link_type=ISI) 

28. 28.Danecek P, Bonfield JK, Liddle J, Marshall J, Ohan V, Pollard MO, Whitwham A, Keane T, McCarthy SA, Davies RM, Li H. 2021. Twelve years of SAMtools and BCFtools. GigaScience 10:giab008.
    
    

29. 29.2024. igvteam/igv-reports. Python. igvteam.
    
    

30. 30.Todokoro D, Tomita H, Inoue T, Ike Y. 2006. Genetic Analysis of Bacteriocin 43 of Vancomycin-Resistant Enterococcus faecium. Appl Environ Microbiol 72:6955–6964.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWVtIjtzOjU6InJlc2lkIjtzOjEwOiI3Mi8xMS82OTU1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDcvMTgvMjAyNC4wNy4xNy4yNDMxMDU5Mi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

31. 31.Nascimento M, Sousa A, Ramirez M, Francisco AP, Carriço JA, Vaz C. 2017. PHYLOViZ 2.0: providing scalable data integration and visualization for multiple phylogenetic inference methods. Bioinformatics 33:128–129.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btw582&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27605102&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F18%2F2024.07.17.24310592.atom) 

32. 32.Sullivan MJ, Petty NK, Beatson SA. 2011. Easyfig: a genome comparison visualizer. Bioinformatics 27:1009–1010.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btr039&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21278367&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F18%2F2024.07.17.24310592.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000289162000018&link_type=ISI) 

33. 33.Cintas LM, Casaus P, Håvarstein LS, Hernández PE, Nes IF. 1997. Biochemical and genetic characterization of enterocin P, a novel sec-dependent bacteriocin from Enterococcus faecium P13 with a broad antimicrobial spectrum. Appl Environ Microbiol 63:4321–4330.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWVtIjtzOjU6InJlc2lkIjtzOjEwOiI2My8xMS80MzIxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDcvMTgvMjAyNC4wNy4xNy4yNDMxMDU5Mi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

34. 34.Casaus P, Nilsen T, Cintas LM, Nes IF, Hernández PE, Holo H. 1997. Enterocin B, a new bacteriocin from Enterococcus faecium T136 which can act synergistically with enterocin A. Microbiology 143:2287– 2294.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1099/00221287-143-7-2287&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9245817&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F18%2F2024.07.17.24310592.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1997XL97200021&link_type=ISI) 

35. 35.Inoue T, Tomita H, Ike Y. 2006. Bac 32, a novel bacteriocin widely disseminated among clinical isolates of Enterococcus faecium. Antimicrob Agents Chemother 50:1202–1212.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWFjIjtzOjU6InJlc2lkIjtzOjk6IjUwLzQvMTIwMiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA3LzE4LzIwMjQuMDcuMTcuMjQzMTA1OTIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

36. 36.Hu C-B, Malaphan W, Zendo T, Nakayama J, Sonomoto K. 2010. Enterocin X, a Novel Two-Peptide Bacteriocin from Enterococcus faecium KU-B5, Has an Antibacterial Spectrum Entirely Different from Those of Its Component Peptides. Appl Environ Microbiol 76:4542–4545.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWVtIjtzOjU6InJlc2lkIjtzOjEwOiI3Ni8xMy80NTQyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDcvMTgvMjAyNC4wNy4xNy4yNDMxMDU5Mi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

37. 37.Cintas LM, Casaus P, Holo H, Hernandez PE, Nes IF, Håvarstein LS. 1998. Enterocins L50A and L50B, Two Novel Bacteriocins from Enterococcus faecium L50, Are Related to Staphylococcal Hemolysins. J Bacteriol 180:1988–1994.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MjoiamIiO3M6NToicmVzaWQiO3M6MTA6IjE4MC84LzE5ODgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wNy8xOC8yMDI0LjA3LjE3LjI0MzEwNTkyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

38. 38.Park SH, Itoh K, Fujisawa T. 2003. Characteristics and identification of enterocins produced by Enterococcus faecium JCM 5804T. J Appl Microbiol 95:294–300.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1046/j.1365-2672.2003.01975.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12859761&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F18%2F2024.07.17.24310592.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000184110900012&link_type=ISI) 

39. 39.Cintas LM, Casaus P, Herranz C, Håvarstein LS, Holo H, Hernández PE, Nes IF. 2000. Biochemical and Genetic Evidence that Enterococcus faecium L50 Produces Enterocins L50A and L50B, the sec-Dependent Enterocin P, and a Novel Bacteriocin Secreted without an N-Terminal Extension Termed Enterocin Q. J Bacteriol 182:6806–6814.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MjoiamIiO3M6NToicmVzaWQiO3M6MTE6IjE4Mi8yMy82ODA2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDcvMTgvMjAyNC4wNy4xNy4yNDMxMDU5Mi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

40. 40. EGUCHI T, KAMINAKA K,  SHIMA J, KAWAMOTO S,  MORI K, CHOI S-H,  DOI K, OHMOMO S, Ogata S. 2001. Isolation and Characterization of Enterocin SE-K4 Produced by Thermophilic Enterococci, Enterococcus faecalis K-4. Biosci Biotechnol Biochem 65:247–253.
    
    

41. 41.Tomita H, Fujimoto S, Tanimoto K, Ike Y. 1996. Cloning and genetic organization of the bacteriocin 31 determinant encoded on the Enterococcus faecalis pheromone-responsive conjugative plasmid pYI17. J Bacteriol 178:3585–3593.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MjoiamIiO3M6NToicmVzaWQiO3M6MTE6IjE3OC8xMi8zNTg1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDcvMTgvMjAyNC4wNy4xNy4yNDMxMDU5Mi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

42. 42.Maricic N, Anderson ES, Opipari AE, Yu EA, Dawid S. 2016. Characterization of a Multipeptide Lantibiotic Locus in Streptococcus pneumoniae. mBio 7:e01656–15.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1128/mBio.01656-15&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26814178&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F18%2F2024.07.17.24310592.atom) 

43. 43.Perez RH, Ishibashi N, Inoue T, Himeno K, Masuda Y, Sawa N, Zendo T, Wilaipun P, Leelawatcharamas V, Nakayama J, Sonomoto K. 2015. Functional Analysis of Genes Involved in the Biosynthesis of Enterocin NKR-5-3B, a Novel Circular Bacteriocin. J Bacteriol 198:291–300.
    
    

44. 44.Aymerich T, Holo H, Håvarstein LS, Hugas M, Garriga M, Nes IF. 1996. Biochemical and genetic characterization of enterocin A from Enterococcus faecium, a new antilisterial bacteriocin in the pediocin family of bacteriocins. Appl Environ Microbiol 62:1676–1682.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWVtIjtzOjU6InJlc2lkIjtzOjk6IjYyLzUvMTY3NiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA3LzE4LzIwMjQuMDcuMTcuMjQzMTA1OTIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

45. 45.O’Keeffe T, Hill C, Ross RP. 1999. Characterization and Heterologous Expression of the Genes Encoding Enterocin A Production, Immunity, and Regulation inEnterococcus faecium DPC1146. Appl Environ Microbiol 65:1506–1515.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWVtIjtzOjU6InJlc2lkIjtzOjk6IjY1LzQvMTUwNiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA3LzE4LzIwMjQuMDcuMTcuMjQzMTA1OTIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

46. 46.Wagner TM, Pöntinen AK, Fenzel CK, Engi D, Janice J, Almeida-Santos AC, Tedim AP, Freitas AR, Peixe L, van Schaik W, Johannessen M, Hegstad K. 2024. Interactions between commensal Enterococcus faecium and Enterococcus lactis and clinical isolates of Enterococcus faecium. FEMS Microbes 5:xtae009.
    
    

47. 47.Tedim AP, Almeida-Santos AC, Lanza VF, Novais C, Coque TM, Freitas AR, Peixe L. 2023. Bacteriocin Distribution Patterns in Enterococcus faecium and Enterococcus lactis: Bioinformatic Analysis Using a Tailored Genomics Framework. bioRxiv doi:10.1101/2023.11.13.566347.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMy4xMS4xMy41NjYzNDd2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA3LzE4LzIwMjQuMDcuMTcuMjQzMTA1OTIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

48. 48.Geldart K, Kaznessis YN. 2017. Characterization of Class IIa Bacteriocin Resistance in Enterococcus faecium. Antimicrob Agents Chemother 61:e02033–16, e02033-16.