Next Article in Journal
Detection of Four Porcine Enteric Coronaviruses Using CRISPR-Cas12a Combined with Multiplex Reverse Transcriptase Loop-Mediated Isothermal Amplification Assay
Previous Article in Journal
Development of Dog Vaccination Strategies to Maintain Herd Immunity against Rabies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Two Novel Lytic Bacteriophages Infecting Enterococcus spp. Are Promising Candidates for Targeted Antibacterial Therapy

1
Scientific and Educational Center “Molecular Bases of Interaction of Microorganisms and Human” of the WCRC “Center for Personalized Medicine”, Institute of Experimental Medicine, 197022 Saint Petersburg, Russia
2
Smorodintsev Research Institute of Influenza, Ministry of Health of the Russian Federation, 197376 Saint Petersburg, Russia
3
Laboratory of Pathomorphology, Almazov National Research Centre, 197341 Saint Petersburg, Russia
*
Authors to whom correspondence should be addressed.
Viruses 2022, 14(4), 831; https://doi.org/10.3390/v14040831
Submission received: 15 March 2022 / Revised: 8 April 2022 / Accepted: 14 April 2022 / Published: 16 April 2022
(This article belongs to the Section Bacterial Viruses)

Abstract

:
The rapid emergence of antibiotic resistance is of major concern globally. Among the most worrying pathogenic bacteria are vancomycin-resistant enterococci. Phage therapy is a highly promising method for controlling enterococcal infections. In this study, we described two virulent tailed bacteriophages possessing lytic activity against Enterococcus faecalis and E. faecium isolates. The SSsP-1 bacteriophage belonged to the Saphexavirus genus of the Siphoviridae family, and the GVEsP-1 bacteriophage belonged to the Schiekvirus genus of Herelleviridae. The genomes of both viruses carried putative components of anti-CRISPR systems and did not contain known genes coding for antibiotic-resistance determinants and virulence factors. The conservative arrangement of protein-coding sequences in Saphexavirus and Schiekvirus genomes taken together with positive results of treating enterococcal peritonitis in an animal infection model imply the potential suitability of GVEsP-1 and SSsP-1 bacteriophages for clinical applications.

1. Introduction

The bacteria of the genus Enterococcus form a part of the human microbiome. In healthy individuals, they reside mainly in the gastrointestinal tract and participate in food digestion [1]. However, adapted to the hospital environment, vancomycin-resistant and multidrug-resistant Enterococcus strains are frequent cause of lethal system blood infections in comparison with antibiotic-sensitive strains [2,3,4]. In leukemia patients, vancomycin-resistant enterococcal infections have associated mortality rates estimated at 57% [5]. Infections caused by antibiotic-resistant strains of enterococci harm individuals and cause damage to society as a whole [6]. The two Enterococcus species encountered most often in multidrug-resistant nosocomial infections are E. faecium and E. faecalis [7,8]. Of these two species, drug-resistant E. faecium is on the ESKAPE list compiled by the Infectious Diseases Society of America, highlighting the utmost importance of developing alternative therapies against this pathogen [9]. Moreover, taking into account a number of factors, including mortality, health-care burden and others, the World Health Organization considered the development of effective drugs against vancomycin-resistant strains of E. faecium to be a high priority [10].
In recent decades, the decreasing rates of progress in the discovery of antibiotics have promoted research interest in bacteriophage therapy [11,12]. Though phage therapy has not become mainstream yet, there are a number of well-described successfully treated cases of surgical infections at present [13,14,15]. Enterococcal viruses are also promising agents for treating oral infections [16] and intestinal dysbiosis [17]. Among their potential uses is development for bio-sanitizing formulations [18]. Since the application of therapeutic bacteriophages is naturally limited by their narrow host ranges and host resistance [19], there is a continuous need for improvements in the spectrum of available bacteriophages.
Both the advantage and disadvantage of phage therapy as compared with antibiotics is its selectivity [20]. Bacteriophages usually target a limited number of strains within bacterial species. Cross-genus tropism of some phages is known, but the examples are quite rare [21]. There are several requirements for therapeutic bacteriophages. The list includes natural origin of a bacterial virus, known taxonomic identity, absence of harmful genetic determinants, absence of transducing activity and high in vitro efficacy [22]. Here we describe two novel virulent Enterococcus bacteriophages. We analyze their host range, morphology, whole-genome sequences and potential to eradicate a model infection in mice.

2. Materials and Methods

2.1. Bacteriophage Isolation and Preparation

Four E. faecalis strains were used for viral propagation in vitro and in vivo experiments (Table 1). Our bacteriophages originated from fecal samples and river water. The isolation of bacterial viruses was performed according to slightly modified classical protocols [23,24]. The protocol for phage isolation from feces included the initial step of resuspension; all other steps were similar. The fecal sample was resuspended by vortexing in sterile SM buffer. The SM buffer consisted of 10 mM NaCl, 10 mM MgSO4, 50 mM Tris-HCl pH 7.5 and 0.05% gelatin. In the next step, 15 mL of the diluted fecal sample or river water were centrifuged at 4500× g for 5 min and filtered through 0.45 µm pore syringe filters with PES membranes (Jet Biofil, YongHe Development Zone, Guangzhou, China). One milliliter of filtrate was supplied with 200 µL 5× Todd-Hewitt broth (HiMedia Laboratories, Mumbai, India). The sample of river water was also supplemented with CaCl2 and MgSO4 to obtain equal final concentrations at 10 mM. Then, the samples were inoculated with 20 µL of 18–24 h bacterial cultures and incubated at 30 °C overnight. The resulting lysate was centrifuged at 4500× g and filtered through 0.45 µm pore syringe filters. The presence of bacteriophages was revealed by spot test with a series of 10-fold dilutions. The pure viral cultures were obtained by triple propagation in sensitive bacterial strains.
Phage stocks used in all described here in vitro and in vivo experiments and DNA sequencing were prepared in the following way. The precipitation of viral particles was performed by supplementing 50 mL of viral cultures with PEG6000 up to 10% w/v and NaCl to a final concentration of 1 M and mixed, then kept at 4 °C overnight. This was followed by centrifugation at 5000× g for 60 min. The resulting pellet was resuspended in SM buffer. The solution was filtered with the use of 0.45 µm pore syringe filters.

2.2. Transmission Electron Microscopy

For transmission electron microscopy, the studied viral stocks were applied to copper grids (300 mesh, Sigma-Aldrich, St. Louis, MO, USA), coated with a collodion film substrate. After adsorption of particles from suspension to the supporting film for 1–2 min, the meshes were washed twice with distilled water. Further, negative contrasting of the samples was carried out for 1–2 min in sodium salt of phosphoric–tungstic acid, 2% aqueous solution (Sigma-Aldrich), pH 7.2. Afterwards, the meshes were dried and examined on a transmission electron microscope JEM 1011 (JEOL, Tokyo, Japan). Instrumental magnification ranged from 50,000× to 200,000×. Electron micrographs were obtained using a Morada high-resolution digital camera (Olympus, Tokyo, Japan). In each case, 20–30 fields of view were examined. Phage size measurements were calculated from TEM images (GVEsP-1, n = 10; SSsP-1, n = 4). The variability of phage dimensions was assessed by calculating the standard deviations of the mean values.

2.3. In Vitro Efficacy

The length of the latent period was assessed by one-step growth experiments, performed as described by Adams [23], with few modifications. To determine the optimal multiplicity of infection (MOI) of phages, E. faecalis strains Serg and CCUG 52538 were infected with SSsP-1 and GVEsP-1 phages, respectively, with different MOIs (0.001, 0.01, 0.1, 1, 10, 100). After six hours of incubation, MOIs of propagated phages were counted in supernatant using the double-layer agar method. For one step growth experiment, 0.1 mL of phage stock (3 × 108 PFU) was mixed with 0.9 mL of bacterial test culture (3 × 109 CFU) (MOI = 0.01). After 10 min of incubation at 37 °C, the mixture was centrifuged at 10,000× g for 5 min, the supernatant was discarded and the pellet was resuspended in a fresh 10 mL of Todd-Hewitt broth. Then, portions of 0.1 mL of broth were taken every 5 min and titrated against a host culture. The experiment was performed in triplicate. To calculate the burst size of each bacteriophage, we used the method described in reference [25].

2.4. Determining the Host Range of Viruses

Host range was determined with the use of an in-house collection of 82 bacterial strains (Table S1). The list of Enterococcus spp. strains included E. faecalis (n = 39), E. faecium (n = 23), E. casseliflavus (n = 1), E. durans (n = 1), E. hirae (n = 1) and E. gallinarum (n = 2). We also tested the phages with Staphylococcus aureus (n = 5), S. epidermidis (n = 1), Streptococcus agalactiae (n = 3), S. pyogenes (n = 4), Bifidobacterium longum (n = 1) and Escherichia coli (n = 1) isolates. Bacterial species identification was carried out using a BactoSCREEN MALDI-TOF MS system (Lytech, Moscow, Russia). All Enterococcus and Streptococcus strains were grown on Todd-Hewitt broth and Todd-Hewitt agar at 37 °C (HiMedia Laboratories, Mumbai, India). Staphylococcus spp. strains were grown on Mannitol salt broth and Mannitol salt agar (SRCAMB, Obolensk, Russia). The Escherichia coli strain was grown on LB agar and LB broth (BioFroxx, Einhausen, Germany). The Bifidobacterium longum strain was tested on Bifidum agar and Bifidum broth (SRCAMB) under anaerobic conditions. Host range was determined by spot test with eight serial 10-fold dilutions of phage stock on a bacterial lawn.

2.5. DNA Isolation and Whole-Genome Sequencing

To isolate viral DNA, the standard phenol/chloroform DNA extraction protocol was used [26]. The phage genome sequences were obtained with the use of the Illumina MiSeq platform (Illumina, San Diego, CA, USA). For GVEsP-1, library preparation was carried out with a Nextera XT DNA Library Preparation Kit (Illumina), resulting in paired-end 300 bp reads. The library for SSsP-1 genome was obtained with the use of a NEBNext Ultra II DNA Library Prep Kit (New England Biolabs, Ipswich, MA, USA) with an average read size of 350 bp. The raw reads were quality controlled using FastQC v0.11 (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/, accessed on 1 May 2020). The viral genomes were assembled de novo by SPAdes 3.13.0 [27]. The nucleotide sequences of the genomes were deposited in the NCBI Nucleotide database. The GenBank/ENA/DDBJ accession number for GVEsP-1 is MZ333462. The accession number for the SSsP-1 genome is MZ333457.

2.6. Bioinformatic Analysis

We visualized the whole-genome structure of the two studied viruses using CGView 1.7 Server [28]. The genetic structures of the genera Schiekvirus and Saphexavirus were visualized by phylogenetic network analysis. For each of the two studied viruses, the sequences of five protein-coding genes were used for BLASTN searches against the NCBI Nucleotide database. The marker genes were supplemented with genes having the widest possible taxonomic coverage. Selected GVEsP-1 sequences coded for the baseplate assembly protein, capsid and scaffold protein, DNA helicase, DNA polymerase and major capsid protein. The studied SSsP-1 sequences included the genes coding for the minor capsid protein, minor structural protein, portal protein, replicative DNA helicase and terminase large subunit. The concatenated alignments of these genes were used to calculate NeighborNet phylogenetic networks in SplitsTree 4.14.2 software [29].
To analyze protein-coding gene synteny in Schiekvirus spp. and Saphexavirus spp., we downloaded the available genomes from GenBank. The genomes were de novo annotated with the use of Prokka 1.14.6 [30]. For each genus, the resulting annotations were used to prepare five datasets for the Synima v 1.0 synteny imaging tool [31]. Four datasets were obtained by dividing all predicted protein-coding sequences from a particular genus into four equal parts according to their GC content, and the fifth dataset contained all CDSs. Synima was run with all five datasets independently and the resulting images were superimposed in a vector graphics editor. Data formatting was performed with the use of custom Python scripts, available at https://github.com/Ivan-Pchelin/scripts-for-synteny-visualization (accessed on 14 February 2022). Potential anti-CRISPR loci were predicted using AcrFinder [32]. The search for known antibiotic-resistance determinants and bacterial virulence factors was performed with the use of ABRicate v 0.8 [33] and its in-built databases ARG-ANNOT, CARD, NCBI AMRFinderPlus, Resfinder and VFDB [34,35,36,37,38].

2.7. Peritoneal Infection Model in Mice

Outbred male mice (n = 54) weighing 20–30 g were acquired from the laboratory breeding nursery of the Russian Academy of Sciences (Rappolovo, Leningrad Region, Russia). Two infection model experiments, one for each virus, were performed. In each case, the mice were divided into three equal groups. In the GVEsP-1 experiment, the groups numbered 8 animals. In the SSsP-1 experiment, each group numbered 10 animals. The Enterococcus faecalis strains CCUG 52538 (viral strain GVEsP-1) and Serg (viral strain SSsP-1) were used for modelling of the infection.
The animals were inoculated intraperitoneally with a lethal dosage of E. faecalis. The bacterial load was 2–4 × 109 CFU for the SSsP-1 experiment and 3–5 × 109 CFU for the GVEsP-1 experiment. In each experiment, the first group of mice did not receive treatment. The second group received 1 mL−1 3 × 109 phage stock per os 3 h after infection. The third group received an intra-abdominal injection of the same phage preparation 3 h after infection. The animal infection model experiment continued for 7 days. During the course of the experiment, lethal cases of infection were autopsied to check the circulation of bacterial pathogens and bacteriophages in blood and bloodstream organs (liver, heart, spleen). The organs were homogenized and cultured to detect viable enterococci. The obtained bacteria were tested for phage sensitivity. After seven days of the experiment, all surviving mice were executed and autopsied as described above. Statistical analysis of the experiment included estimation of the survival distributions by the Kaplan–Meier method and further comparison of the distributions by performing a logrank test in R 4.1.2 [39] with the use of the survival 3.3–1 package [40]. The survival curves were visualized in R using the survminer 0.4.9 package [41].

3. Results and Discussion

3.1. Isolation and Identification of Bacteriophages

The bacteriophage SSsP-1 was isolated from a sample of feces enriched with the E. faecalis strain Serg. The virus GVEsP-1 was isolated from a water sample from To Lich River (Hanoi, Vietnam) using the bacterial strain E. faecalis 5arctic. Taxonomic identification through MegaBLAST searches against viral sequences of the NCBI Nucleotide database placed GVEsP-1 in the Shiekvirus genus of the family Herelleviridae. The phage SSsP-1 clustered within the Saphexavirus genus of Siphoviridae. The identity of the studied genomes and their closest matches in the database was 92% for GVEsP-1 (compared to Enterococcus phage vB_EfaM_A2) and 77% for SSsP-1 (compared to Enterococcus phage vB_EfaS_IME198). Therefore, both bacteriophages belonged to undescribed viral species, given the currently applied 95% threshold [42].
From the early stages of research, the phages showed different behavior and different tropism to enterococcal strains. GVEsP-1 formed small colonies on a double agar layer. It was active against 61% of E. faecalis strains and 22% of E. faecium strains from our collection. On the double agar plates with sensitive cultures, the phage SSsP-1 initially formed smooth plaques. After a number of passages its plaques changed morphology, being surrounded by a pronounced halo zone. The phage SSsP-1 was active against 36% of our E. faecalis isolates. Neither phage infected other tested bacterial species, including E. gallinarum, E. casseliflavus, E. hirae, E. durans, Staphylococcus spp., Streptococcus spp., Bifidobacterium longum and Escherichia coli. In all cases of successful lysis of bacterial cultures, spot tests of serial dilutions revealed the presence of viable viral progeny. Known E. faecalis viruses were active against 7.6–70.5% of bacterial strains [43,44,45,46,47,48,49]. Therefore, the host range of GVEsP-1 can be considered broad, and the host range of SSsP-1 can be thought of as moderate.
Searches in the NCBI and the International Committee on Taxonomy of Viruses (ICTV) databases revealed the presence of one Streptococcus phage SP-QS1 with a Streptococcus pneumoniae host in the Saphexavirus genus [Almaghrabi et al., unpublished, GenBank accession NC_021868]. Formerly, enterococci were classified as group D streptococci by Lancefield [50]. In 1984, they were separated into their own genus after studies based on nucleic acid hybridization showed a more distant relationship to enterococci [51]. Therefore, we tested the activity of SSsP-1 and GVEsP-1 against streptococcal and staphylococcal strains. As mentioned earlier, the two phages did not infect any streptococcal or staphylococcal strains. Bifidobacterium longum and Escherichia coli are phylogenetically distant from Enterococcus and belong to Proteobacteria and Actinobacteria phyla, respectively. Given that ICTV considers host range as one of the criteria of taxonomic classification of viruses [52], the phages expectedly did not infect these strains.

3.2. Phage Life Cycle and Morphology

The optimal MOIs for both phages were estimated at 0.01. One-step growth experiments for both phages did not show any prominent differences between Saphexavirus and Shiekvirus bacteriophages. In GVEsP-1, the length of the latent period was 20 min (Figure 1A). In SSsP-1, this phase took 18 min (Figure 1B). The average burst size for the phages SSsP-1 and GVEsP-1 was 66 ± 6 and 94 ± 4 PFU/cell, respectively. The GVEsP-1 phage possessed an icosahedral capsid with dimensions along the main axis of the particles of 77.8 ± 10.4 nm and dimensions along the cross axis of 77.7 ± 4.6 nm. A long contractile non-flexible needle-like tail in a sheath and a baseplate receptor under the sheath measured 160.5 ± 15.4 nm in length. No whiskers or legs were observed (Figure 2A). The capsid of SSsP-1 had an oval shape and was 87.7 ± 4.1 nm long, its width 43.2 ± 1.9 nm. The virus had a long, flexible tail measuring 117.5 ± 5.4 nm. A putative phage receptor was located on the baseplate (Figure 2B). These data were in agreement with a series of studies on other Saphexavirus and Schiekvirus viruses [44,53,54,55,56].

3.3. Genome Structure and Conservation

The GVEsP-1 genome sequence had a length of 149,913 bp with 194 predicted coding sequences, while the SSsP-1 genome measured 57,270 bp and contained 93 predicted CDSs. The GC content was 37% in GVEsP-1 and 40% in SSsP-1. The coding sequences in both genomes resided on the strands with positive GC skew (Figure 3). Additionally, both phages contained a cluster of genes coding for tRNAs. In the GVEsP-1 genome, there were 25 tRNA-coding genes, which is more than can be found in the genome of the closely related phage EFDG1. The SSsP-1 genome contained three tRNA genes. Neither phage carried integration-related genes and therefore they were obligately lytic.
To visualize the evolutionary relationships between the studied viruses within the genera, we calculated phylogenetic networks, using five protein-coding genes in both cases. The obtained alignments covered all the sequences available in GenBank for the Shiekvirus and Saphexavirus genomes. The Shiekvirus phylogenetic network had elements with tree-like structures, whereas in the Saphexavirus network, reticulation events prevailed (Figure 4). In both cases, the structure of the networks suggested either a significant impact of horizontal gene exchange in bacteriophage evolution [57] or the recent origin of the viruses from a common ancestor. The latter scenario may be supported by the high similarity of the phage SSsP-1 and the virus Entf1, isolated in the course of another Russian study (Rubalskii et al., unpublished, GenBank accession MK800154). However, recently diverged bacteriophage genomes are, at the same time, subjected to horizontal exchange of genetic information [58]. Still, a search for antibiotic-resistance determinants and virulence factors with the ABRicate tool in Schiekvirus and Saphexavirus, including the GVEsP-1 and SSsP-1 genomes, retrieved zero matches, which implies the safety of their potential clinical use.
The whole-genome synteny in Schiekvirus and Saphexavirus revealed conservative sets and orders of protein-coding sequences, implying high degrees of structural conservation. The number of rearrangement events was neglectable (Figure 5). Therefore, the studied bacteriophages likely share a comparable biology, replication cycle and ecology with their congeneric viruses.

3.4. Potential Anti-CRISPR Loci

CRISPR–Cas systems provide a defense against heterologous genetic material, such as viruses, plasmids and other mobile genetic elements, for bacteria and archaea. They were first described in 2013 in several bacteriophages of Pseudomonas aeruginosa [59]. On the viral side, anti-CRISPR systems protect phage DNA through six theoretical mechanisms, including prevention of the insertion of viral DNA into the genome of a host cell, disruption of the synthesis of Cas proteins, blockage of crRNA synthesis, inhibition of crRNA loading onto Cas proteins, inhibition of DNA binding by Cas proteins and inhibition of DNA cleavage. The two latter mechanisms have been studied in vitro by molecular biology techniques [60]. The best-described anti-CRISPR proteins are known from the Siphoviridae phage family, myoviridae phages and prophages [59,61]. SSsP-1 and GVEsP-1 belong to siphoviridae and myoviridae morphotypes of obligately lytic bacteriophages, respectively.
In the SSsP-1 bacteriophage genome, most proteins identified by AcrFinder as potential anti-CRISPR loci did not have functional annotations. Still, in putative anti-CRISPR loci, pAcrS1, pAcrS3 and pAcrS4, there were proteins annotated as DNA-binding proteins and HNH homing endonucleases (Table 2). DNA-binding protein can probably block DNA binding by Cas. HNH homing endonuclease can be a protein with a HNH domain that is bound to the Cas protein HNH domain as it was described by Harrington et al. [62]. The anti-CRISPR loci in the GVEsP-1 genome were predicted with less confidence, so it was unclear whether they indeed had anti-CRISPR functions.
To sum up, well-known anti-CRISPR mechanisms interrupt CRISPR–Cas immunity at stages of expression and interference. Enterococcus faecalis is a regular carrier of a type II CRISPR–Cas system [63]. Well-known anti-CRISPR type II proteins exploit mechanisms of DNA-binding inhibition through binding to the HNH domain of Cas nucleases [62]. In the studied bacteriophage genomes, we found putative components of an anti-CRISPR system involving inhibition of DNA binding by Cas proteins. Potentially, the presence of anti-CRISPR proteins can make therapeutic enterococcal bacteriophages more efficient from a practical point of view.

3.5. Mouse Infection Experiments

In our peritoneal infection model experiments, untreated animals infected with E. faecalis CCUG 52538 died by the fourth day after inoculation. Six out of ten mice infected with E. faecalis Serg died by the end of the first day. All treated mice survived until the end of the experiments (Figure 6). Therefore, phage administration route did not affect the mortality rate of mice.
The differences in survival distributions between the treated and untreated animal groups were statistically significant (p < 0.001 by the logrank test), implying the potential of the studied bacteriophages for therapeutic applications [64]. The positive results of our in vivo experiment corroborate an earlier conclusion that enterococcal viruses from the families Herelleviridae and Siphoviridae are good candidates for phage therapy [65]. Our culture tests with blood and bloodstream organs of autopsied animals did not reveal the presence of viable phage particles. This can be explained by the deactivation of virions by the animal immunity system [66].

4. Conclusions

The high activities of two novel tailed bacteriophages against Enterococcus spp. were demonstrated. Both phages were able to protect mice from lethal enterococcal infection. The phages had conservative genomic structures and contained putative components of anti-CRISPR systems. Given the relative genetic stability of viruses from the genera Saphexavirus and Schiekvirus, the results of our animal infection experiment may indicate the suitability of studied bacteriophages for treating septic enterococcal infections.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v14040831/s1, Table S1: Host range of GVEsP-1 and SSsP-1.

Author Contributions

Study design, A.E.G. and P.V.T.; investigation, P.V.T. and A.N.G.; formal analysis, P.V.T., I.M.P. and D.V.A.; writing—original draft preparation P.V.T. and I.M.P.; writing—review and editing, I.M.P., A.V.D., D.V.A. and A.E.G.; supervision and funding acquisition, D.V.A., A.E.G., O.V.S. and A.V.D. All authors have read and agreed to the published version of the manuscript.

Funding

The work was supported by the Ministry of Science and Higher Education of the Russian Federation, no. 075-15-2020-902, Federal State Budgetary Scientific Institution “Institute of Experimental Medicine” (FSBSI “IEM”) Scientific and Educational Center “Molecular bases of interaction of microorganisms and human” of the World-Class Research Center “Center for Personalized Medicine” at FSBSI “IEM”.

Institutional Review Board Statement

The experiments were performed in compliance with the principles of humanity (of the European Communities no. 86/609 EU) and were approved by the local ethics committee at FSBSI “IEM”, approval no. 3/21, issued on 27 October 2021.

Informed Consent Statement

Not applicable.

Data Availability Statement

The nucleotide sequences for the phages GVEsP-1 and SSsP-1 have been deposited in the NCBI Nucleotide database under the GenBank/ENA/DDBJ accession numbers MZ333462 and MZ333457, respectively.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study, in the collection, analyses, or interpretation of data, in the writing of the manuscript or in the decision to publish the results.

References

  1. Rowland, I.; Gibson, G.; Heinken, A.; Scott, K.; Swann, J.; Thiele, I.; Tuohy, K. Gut microbiota functions: Metabolism of nutrients and other food components. Eur. J. Nutr. 2018, 57, 1–24. [Google Scholar] [CrossRef] [Green Version]
  2. De Oliveira, D.M.P.; Forde, B.M.; Kidd, T.J.; Harris, P.N.A.; Schembri, M.A.; Beatson, S.A.; Paterson, D.L.; Walker, M.J. Antimicrobial resistance in ESKAPE pathogens. Clin. Microbiol. Rev. 2020, 33, e00181-19. [Google Scholar] [CrossRef]
  3. Gao, W.; Howden, B.P.; Stinear, T.P. Evolution of virulence in Enterococcus faecium, a hospital-adapted opportunistic pathogen. Curr. Opin. Microbiol. 2018, 41, 76–82. [Google Scholar] [CrossRef]
  4. Fraser, J.L.; Mwatondo, A.; Alimi, Y.H.; Varma, J.K.; Vilas, V.J.D.R. Healthcare-associated outbreaks of bacterial infections in Africa, 2009-2018: A review. Int. J. Infect. Dis. 2021, 103, 469–477. [Google Scholar] [CrossRef]
  5. Kirkizlar, T.A.; Akalin, H.; Kirkizlar, O.; Ozkalemkas, F.; Ozkocaman, V.; Kazak, E.; Ozakin, C.; Bulbul, E.N.; Ozboz, E.S.; Ali, R. Vancomycin-resistant enterococci infection and predisposing factors for infection and mortality in patients with acute leukaemia and febrile neutropenia. Leuk. Res. 2020, 99, 106463. [Google Scholar] [CrossRef]
  6. Serra-Burriel, M.; Keys, M.; Campillo-Artero, C.; Agodi, A.; Barchitta, M.; Gikas, A.; Palos, C.; López-Casasnovas, G. Impact of multi-drug resistant bacteria on economic and clinical outcomes of healthcare-associated infections in adults: Systematic review and meta-analysis. PLoS ONE. 2020, 15, e0227139. [Google Scholar] [CrossRef]
  7. Gilmore, M.S.; Lebreton, F.; van Schaik, W. Genomic transition of enterococci from gut commensals to leading causes of multidrug-resistant hospital infection in the antibiotic era. Curr. Opin. Microbiol. 2013, 16, 10–16. [Google Scholar] [CrossRef] [Green Version]
  8. Zhen, X.; Lundborg, C.S.; Sun, X.; Hu, X.; Dong, H. Economic burden of antibiotic resistance in ESKAPE organisms: A systematic review. Antimicrob. Resist. Infect. Control. 2019, 8, 137. [Google Scholar] [CrossRef] [Green Version]
  9. Santajit, S.; Indrawattana, N. Mechanisms of antimicrobial resistance in ESKAPE pathogens. Biomed Res. Int. 2016, 2016, 2475067. [Google Scholar] [CrossRef] [Green Version]
  10. Tacconelli, E.; Carrara, E.; Savoldi, A.; Harbarth, S.; Mendelson, M.; Monnet, D.L.; Pulcini, C.; Kahlmeter, G.; Kluytmans, J.; Carmeli, Y.; et al. Discovery, research, and development of new antibiotics: The WHO priority list of antibiotic-resistant bacteria and tuberculosis. Lancet. Infect. Dis. 2018, 18, 318–327. [Google Scholar] [CrossRef]
  11. Baëtz, B.; Boudrioua, A.; Hartke, A.; Giraud, C. Alternatives to fight vancomycin-resistant staphylococci and enterococci. Antibiotics 2021, 10, 1116. [Google Scholar] [CrossRef]
  12. McCallin, S.; Brüssow, H. Clinical Trials of Bacteriophage Therapeutics. In Bacteriophages: Biology, Technology, Therapy; Springer: Berlin/Heidelberg, Germany, 2021; pp. 1099–1127. [Google Scholar] [CrossRef]
  13. Neuts, A.S.; Berkhout, H.J.; Hartog, A.; Goosen, J.H.M. Bacteriophage therapy cures a recurrent Enterococcus faecalis infected total hip arthroplasty? A case report. Acta Orthop. 2021, 20, 678–680. [Google Scholar] [CrossRef]
  14. Paul, K.; Merabishvili, M.; Hazan, R.; Christner, M.; Herden, U.; Gelman, D.; Khalifa, L.; Yerushalmy, O.; Coppenhagen-Glazer, S.; Harbauer, T.; et al. Bacteriophage rescue therapy of a vancomycin-resistant Enterococcus faecium infection in a one-year-old child following a third liver transplantation. Viruses 2021, 13, 1785. [Google Scholar] [CrossRef]
  15. Rubalskii, E.; Ruemke, S.; Salmoukas, C.; Boyle, E.C.; Warnecke, G.; Tudorache, I.; Shrestha, M.; Schmitto, J.D.; Martens, A.; Rojas, S.V.; et al. Bacteriophage therapy for critical infections related to cardiothoracic surgery. Antibiotics 2020, 9, 232. [Google Scholar] [CrossRef]
  16. Khalifa, L.; Shlezinger, M.; Beyth, S.; Houri-Haddad, Y.; Coppenhagen-Glazer, S.; Beyth, N.; Hazan, R. Phage therapy against Enterococcus faecalis in dental root canals. J. Oral Microbiol. 2016, 8, 32157. [Google Scholar] [CrossRef] [Green Version]
  17. Chatterjee, A.; Johnson, C.N.; Luong, P.; Hullahalli, K.; McBride, S.W.; Schubert, A.M.; Palmer, K.L.; Carlson, P.E., Jr.; Duerkop, B.A. Bacteriophage resistance alters antibiotic-mediated intestinal expansion of enterococci. Infect. Immun. 2019, 87, e00085-19. [Google Scholar] [CrossRef] [Green Version]
  18. Bhandare, S.; Goodridge, L. Bacteriophages as bio-sanitizers in food production and healthcare settings. In Bacteriophages: Biology, Technology, Therapy; Springer: Berlin/Heidelberg, Germany, 2021; pp. 769–788. [Google Scholar] [CrossRef]
  19. Taati Moghadam, M.; Amirmozafari, N.; Shariati, A.; Hallajzadeh, M.; Mirkalantari, S.; Khoshbayan, A.; Masjedian Jazi, F. How phages overcome the challenges of drug resistant bacteria in clinical infections. Infect. Drug Resist. 2020, 13, 45–61. [Google Scholar] [CrossRef] [Green Version]
  20. Gordillo Altamirano, F.L.; Barr, J.J. Phage therapy in the postantibiotic era. Clin. Microbiol. Rev. 2019, 32, e00066-18. [Google Scholar] [CrossRef] [Green Version]
  21. Olsen, R.H.; Siak, J.S.; Gray, R.H. Characteristics of PRD1, a plasmid-dependent broad host range DNA bacteriophage. J. Virol. 1974, 14, 689–699. [Google Scholar] [CrossRef] [Green Version]
  22. Pirnay, J.P.; Blasdel, B.G.; Bretaudeau, L.; Buckling, A.; Chanishvili, N.; Clark, J.R.; Corte-Real, S.; Debarbieux, L.; Dublanchet, A.; De Vos, D.; et al. Quality and safety requirements for sustainable phage therapy products. Pharm. Res. 2015, 32, 2173–2179. [Google Scholar] [CrossRef] [Green Version]
  23. Adams, M.H. Bacteriophages; Interscience Publishers: New York, NY, USA; London, UK, 1959. [Google Scholar]
  24. Kutter, E.; Sulakvelidze, A. Bacteriophages: Biology and Applications; CRC Press: Boca Raton, FL, USA, 2004. [Google Scholar]
  25. Sullivan, M. Calculating Burst Size. Available online: https://www.protocols.io/view/Calculating-burst-size-3byl47erlo5d/v1.dx.doi.org/10.17504/protocols.io.ddm245 (accessed on 8 April 2022).
  26. Maniatis, T.; Fritsch, E.F.; Sambrook, J. Molecular Cloning. A Laboratory Manual; Cold Spring Harbor Laboratory: Cold Spring Harbor, NY, USA, 1982. [Google Scholar]
  27. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef] [Green Version]
  28. Stothard, P.; Wishart, D.S. Circular genome visualization and exploration using CGView. Bioinformatics 2005, 21, 537–539. [Google Scholar] [CrossRef] [Green Version]
  29. Huson, D.H.; Bryant, D. Application of phylogenetic networks in evolutionary studies. Mol. Biol. Evol. 2006, 23, 254–267. [Google Scholar] [CrossRef]
  30. Seemann, T. Prokka: Rapid prokaryotic genome annotation. Bioinformatics 2014, 30, 2068–2069. [Google Scholar] [CrossRef]
  31. Farrer, R.A. Synima: A synteny imaging tool for annotated genome assemblies. BMC Bioinformatics 2017, 18, 507. [Google Scholar] [CrossRef]
  32. Yi, H.; Huang, L.; Yang, B.; Gomez, J.; Zhang, H.; Yin, Y. AcrFinder: Genome mining anti-CRISPR operons in prokaryotes and their viruses. Nucleic Acids Res. 2020, 48, W358–W365. [Google Scholar] [CrossRef]
  33. Seemann, T. ABRicate. Available online: https://github.com/tseemann/abricate (accessed on 27 December 2021).
  34. Gupta, S.K.; Padmanabhan, B.R.; Diene, S.M.; Lopez-Rojas, R.; Kempf, M.; Landraud, L.; Rolain, J.M. ARG-ANNOT, a new bioinformatic tool to discover antibiotic resistance genes in bacterial genomes. Antimicrob. Agents Chemother. 2014, 58, 212–220. [Google Scholar] [CrossRef] [Green Version]
  35. Jia, B.; Raphenya, A.R.; Alcock, B.; Waglechner, N.; Guo, P.; Tsang, K.K.; Lago, B.A.; Dave, B.M.; Pereira, S.; Sharma, A.N.; et al. CARD 2017: Expansion and model-centric curation of the comprehensive antibiotic resistance database. Nucleic Acids Res. 2017, 45, D566–D573. [Google Scholar] [CrossRef]
  36. Feldgarden, M.; Brover, V.; Haft, D.H.; Prasad, A.B.; Slotta, D.J.; Tolstoy, I.; Tyson, G.H.; Zhao, S.; Hsu, C.H.; McDermott, P.F. Validating the AMRFinder tool and resistance gene database by using antimicrobial resistance genotype-phenotype correlations in a collection of isolates. Antimicrob. Agents Chemother. 2019, 63, e00483-19. [Google Scholar] [CrossRef] [Green Version]
  37. Zankari, E.; Hasman, H.; Cosentino, S.; Vestergaard, M.; Rasmussen, S.; Lund, O.; Aarestrup, F.M.; Larsen, M.V. Identification of acquired antimicrobial resistance genes. J. Antimicrob. Chemother. 2012, 67, 2640–2644. [Google Scholar] [CrossRef]
  38. Chen, L.; Zheng, D.; Liu, B.; Yang, J.; Jin, Q. VFDB 2016: Hierarchical and refined dataset for big data analysis—10 years on. Nucleic Acids Res. 2016, 44, D694–D697. [Google Scholar] [CrossRef]
  39. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2021; Available online: https://www.R-project.org/ (accessed on 8 November 2021).
  40. Therneau, T. A Package for Survival Analysis in R. R Package Version 3.3-1. 2022. Available online: https://CRAN.R-project.org/package=survival (accessed on 7 April 2022).
  41. Kassambara, A.; Kosinski, M.; Biecek, P. Survminer: Drawing Survival Curves Using ‘ggplot2’. R Package Version 0.4.9. 2021. Available online: https://CRAN.R-project.org/package=survminer (accessed on 7 April 2022).
  42. Turner, D.; Kropinski, A.M.; Adriaenssens, E.M. A roadmap for genome-based phage taxonomy. Viruses 2021, 13, 506. [Google Scholar] [CrossRef]
  43. Rigvava, S.; Tchgkonia, I.; Jgenti, D.; Dvalidze, T.; Carpino, J.; Goderdzishvili, M. Comparative analysis of the biological and physical properties of Enterococcus faecalis bacteriophage vB_EfaS_GEC-EfS_3 and Streptococcus mitis bacteriophage vB_SmM_GEC-SmitisM_2. Can. J. Microbiol. 2013, 59, 18–21. [Google Scholar] [CrossRef]
  44. Zhang, W.; Mi, Z.; Yin, X.; Fan, H.; An, X.; Zhang, Z.; Chen, J.; Tong, Y. Characterization of Enterococcus faecalis phage IME-EF1 and its endolysin. PLoS ONE 2013, 8, e80435. [Google Scholar] [CrossRef]
  45. Eladero, V.; Egomez-Sordo, C.; Esanchez-Llana, E.; del Rio, B.; Redruello, B.; Efernandez, M.; Martin, M.C.; Alvarez, M.A. Q69 (an E. faecalis-infecting bacteriophage) as a biocontrol agent for reducing tyramine in dairy products. Front. Microbiol. 2016, 7, 445. [Google Scholar] [CrossRef]
  46. Cheng, M.; Liang, J.; Zhang, Y.; Hu, L.; Gong, P.; Cai, R.; Zhang, L.; Zhang, H.; Ge, J.; Ji, Y.; et al. The Bacteriophage EF-P29 efficiently protects against lethal vancomycin-resistant Enterococcus faecalis and alleviates gut microbiota imbalance in a murine bacteremia model. Front. Microbiol. 2017, 8, 837. [Google Scholar] [CrossRef]
  47. Rahmat Ullah, S.; Andleeb, S.; Raza, T.; Jamal, M.; Mehmood, K. Effectiveness of a lytic phage SRG1 against vancomycin-resistant Enterococcus faecalis in compost and soil. Biomed. Res. Int. 2017, 2017, 9351017. [Google Scholar] [CrossRef] [Green Version]
  48. Wang, R.; Xing, S.; Zhao, F.; Li, P.; Mi, Z.; Shi, T.; Liu, H.; Tong, Y. Characterization and genome analysis of novel phage vB_EfaP_IME195 infecting Enterococcus faecalis. Virus Genes 2018, 54, 804–811. [Google Scholar] [CrossRef]
  49. Lee, D.; Im, J.; Na, H.; Ryu, S.; Yun, C.H.; Han, S.H. The novel Enterococcus phage vB_EfaS_HEf13 has broad lytic activity against clinical isolates of Enterococcus faecalis. Front Microbiol. 2019, 10, 2877. [Google Scholar] [CrossRef]
  50. Lancefield, R.C. A serological differentiation of human and other groups of hemolytic streptococci. J. Exp. Med. 1933, 57, 571–595. [Google Scholar] [CrossRef]
  51. Fiore, E.; Van Tyne, D.; Gilmore, M.S. Pathogenicity of enterococci. Microbiol. Spectr. 2019, 7. [Google Scholar] [CrossRef]
  52. Barylski, J.; Kropinski, A.M.; Alikhan, N.F.; Adriaenssens, E.M.; Ictv Report Consortium. ICTV virus taxonomy profile: Herelleviridae. J. Gen. Virol. 2020, 101, 362–363. [Google Scholar] [CrossRef]
  53. Yang, D.; Chen, Y.; Sun, E.; Hua, L.; Peng, Z.; Wu, B. Characterization of a lytic bacteriophage vB_EfaS_PHB08 harboring endolysin Lys08 against Enterococcus faecalis biofilms. Microorganisms 2020, 8, 1332. [Google Scholar] [CrossRef]
  54. Wang, Q.; Liu, N. Complete genome sequence of bacteriophage EFC1 infecting Enterococcus faecalis from chicken. 2021. preprint. Available online: https://doi.org/10.21203/rs.3.rs-1118970/v1 (accessed on 22 December 2021).
  55. Xiang, Y.; Li, W.; Song, F.; Yang, X.; Zhou, J.; Yu, H.; Ji, X.; Wei, Y. Biological characteristics and whole-genome analysis of the Enterococcus faecalis phage PEf771. Can. J. Microbiol. 2020, 66, 505–520. [Google Scholar] [CrossRef]
  56. Khalifa, L.; Brosh, Y.; Gelman, D.; Coppenhagen-Glazer, S.; Beyth, S.; Poradosu-Cohen, R.; Que, Y.A.; Beyth, N.; Hazan, R. Targeting Enterococcus faecalis biofilms with phage therapy. Appl. Environ. Microbiol. 2015, 81, 2696–2705. [Google Scholar] [CrossRef] [Green Version]
  57. Lima-Mendez, G.; Van Helden, J.; Toussaint, A.; Leplae, R. Reticulate representation of evolutionary and functional relationships between phage genomes. Mol. Biol. Evol. 2008, 25, 762–777. [Google Scholar] [CrossRef] [Green Version]
  58. Casjens, S.R. Comparative genomics and evolution of the tailed-bacteriophages. Curr. Opin. Microbiol. 2005, 8, 451–458. [Google Scholar] [CrossRef]
  59. Maxwell, K.L. The anti-CRISPR story: A battle for survival. Mol. Cell. 2017, 68, 8–14. [Google Scholar] [CrossRef] [Green Version]
  60. Bondy-Denomy, J.; Pawluk, A.; Maxwell, K.L.; Davidson, A.R. Bacteriophage genes that inactivate the CRISPR/Cas bacterial immune system. Nature 2013, 493, 429–432. [Google Scholar] [CrossRef] [Green Version]
  61. Malone, L.M.; Warring, S.L.; Jackson, S.A.; Warnecke, C.; Gardner, P.P.; Gumy, L.F.; Fineran, P.C. A jumbo phage that forms a nucleus-like structure evades CRISPR–Cas DNA targeting but is vulnerable to type III RNA-based immunity. Nat. Microbiol. 2020, 5, 48–55. [Google Scholar] [CrossRef]
  62. Harrington, L.B.; Doxzen, K.W.; Ma, E.; Liu, J.J.; Knott, G.J.; Edraki, A.; Garcia, B.; Amrani, N.; Chen, J.S.; Cofsky, J.C.; et al. A broad-spectrum inhibitor of CRISPR-Cas9. Cell 2017, 170, 1224–1233.e15. [Google Scholar] [CrossRef] [Green Version]
  63. Lyons, C.; Raustad, N.; Bustos, M.A.; Shiaris, M. Incidence of Type II CRISPR1-Cas systems in Enterococcus is species-dependent. PLoS ONE 2015, 10, e0143544. [Google Scholar] [CrossRef] [Green Version]
  64. Goh, H.M.S.; Yong, M.H.A.; Chong, K.K.L.; Kline, K.A. Model systems for the study of enterococcal colonization and infection. Virulence. 2017, 8, 1525–1562. [Google Scholar] [CrossRef] [Green Version]
  65. Olsen, N.S.; Nielsen, K.J.; Plöger, M.; Kot, W.; Hansen, L.H. Enterococcus phage Nonaheksakonda infecting clinical isolates of Enterococcus faecalis represents a new lineage in the family Siphoviridae. Arch. Virol. 2021, 166, 593–599. [Google Scholar] [CrossRef]
  66. Dąbrowska, K. Phage therapy: What factors shape phage pharmacokinetics and bioavailability? Systematic and critical review. Med. Res. Rev. 2019, 39, 2000–2025. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The length of the latent period in the development of bacteriophages assessed by one-step growth experiments. (A) Enterococcus phage GVEsP-1. (B) Enterococcus phage SSsP-1. The shaded areas delineate standard deviation in three biological replicates.
Figure 1. The length of the latent period in the development of bacteriophages assessed by one-step growth experiments. (A) Enterococcus phage GVEsP-1. (B) Enterococcus phage SSsP-1. The shaded areas delineate standard deviation in three biological replicates.
Viruses 14 00831 g001
Figure 2. Morphology of the bacteriophages as revealed by transmission electron microscopy. (A) Enterococcus phage GVEsP-1. (B) Enterococcus phage SSsP-1. Scale bars = 100 nm.
Figure 2. Morphology of the bacteriophages as revealed by transmission electron microscopy. (A) Enterococcus phage GVEsP-1. (B) Enterococcus phage SSsP-1. Scale bars = 100 nm.
Viruses 14 00831 g002
Figure 3. Genome organization of the bacteriophages. (A) Enterococcus phage GVEsP-1. (B) Enterococcus phage SSsP-1.
Figure 3. Genome organization of the bacteriophages. (A) Enterococcus phage GVEsP-1. (B) Enterococcus phage SSsP-1.
Viruses 14 00831 g003
Figure 4. Evolutionary relationships of studied bacteriophages within Schiekvirus and Saphexavirus visualized by NeighborNet phylogenetic networks. (A) The position of GVEsP-1 within Schiekvirus. (B) The position of SSsP-1 within Saphexavirus. The viruses described in the present work are highlighted in purple.
Figure 4. Evolutionary relationships of studied bacteriophages within Schiekvirus and Saphexavirus visualized by NeighborNet phylogenetic networks. (A) The position of GVEsP-1 within Schiekvirus. (B) The position of SSsP-1 within Saphexavirus. The viruses described in the present work are highlighted in purple.
Viruses 14 00831 g004
Figure 5. Protein-coding gene synteny in Schiekvirus (A) and Saphexavirus (B). The black axes correspond to the whole-genome sequences, whereas the colored connections indicate homologous coding sequences, determined by BLAST searches. All CDSs in the dataset were divided into quartiles according to their GC content. Color scale from purple to green indicates homology between the sequences within the same quartiles. The connections between CDSs from different quartiles are colored orange.
Figure 5. Protein-coding gene synteny in Schiekvirus (A) and Saphexavirus (B). The black axes correspond to the whole-genome sequences, whereas the colored connections indicate homologous coding sequences, determined by BLAST searches. All CDSs in the dataset were divided into quartiles according to their GC content. Color scale from purple to green indicates homology between the sequences within the same quartiles. The connections between CDSs from different quartiles are colored orange.
Viruses 14 00831 g005
Figure 6. Survival curves in the animal infection model experiment. (A) Mice infected with E. faecalis strain CCUG 52538 and treated with the phage GVEsP-1. (B) Mice infected with E. faecalis strain Serg and treated with the phage SSsP-1.
Figure 6. Survival curves in the animal infection model experiment. (A) Mice infected with E. faecalis strain CCUG 52538 and treated with the phage GVEsP-1. (B) Mice infected with E. faecalis strain Serg and treated with the phage SSsP-1.
Viruses 14 00831 g006
Table 1. Enterococcus faecalis strains used for experiments.
Table 1. Enterococcus faecalis strains used for experiments.
IsolateSource
SergUrine of patient with urinary tract infection, Saint Petersburg, Russia
5arcticOrnithogenic soil associated with Rissa tridactyla, Svalbard
ATCC 29212Swedish Institute for Infectious Disease Control (SMI)
CCUG 52538Swedish Institute for Infectious Disease Control (SMI)
Table 2. The putative anti-CRISPR loci detected by ArcFinder. The described version of the SSsP-1 genome is MZ333457.1; the version of the GVEsP-1 genome is MZ333462.1.
Table 2. The putative anti-CRISPR loci detected by ArcFinder. The described version of the SSsP-1 genome is MZ333457.1; the version of the GVEsP-1 genome is MZ333462.1.
PhagePutative Anti-CRISPR LocusStrandNumber of ORFsStartEndKnown ORF Annotations
SSsP-1pAcrS1-105293163DNA-binding protein
SSsP-1pAcrS2-9561111,590DNA-binding protein
SSsP-1pAcrS3-214,07114,808HNH endonuclease
SSsP-1pAcrS4-422,88024,454HNH homing endonuclease
SSsP-1pAcrS5+225,84226,812HNH homing endonuclease
SSsP-1pAcrS6-754,58556,973None
GVEsP-1pAcrG1+249475568DNA-binding protein
GVEsP-1pAcrG2+1122,50126,731Phosphoesterase
GVEsP-1pAcrG3+630,68832,712None
GVEsP-1pAcrG4+1137,91840,949None
GVEsP-1pAcrG5-352,26953,288None
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tkachev, P.V.; Pchelin, I.M.; Azarov, D.V.; Gorshkov, A.N.; Shamova, O.V.; Dmitriev, A.V.; Goncharov, A.E. Two Novel Lytic Bacteriophages Infecting Enterococcus spp. Are Promising Candidates for Targeted Antibacterial Therapy. Viruses 2022, 14, 831. https://doi.org/10.3390/v14040831

AMA Style

Tkachev PV, Pchelin IM, Azarov DV, Gorshkov AN, Shamova OV, Dmitriev AV, Goncharov AE. Two Novel Lytic Bacteriophages Infecting Enterococcus spp. Are Promising Candidates for Targeted Antibacterial Therapy. Viruses. 2022; 14(4):831. https://doi.org/10.3390/v14040831

Chicago/Turabian Style

Tkachev, Pavel V., Ivan M. Pchelin, Daniil V. Azarov, Andrey N. Gorshkov, Olga V. Shamova, Alexander V. Dmitriev, and Artemiy E. Goncharov. 2022. "Two Novel Lytic Bacteriophages Infecting Enterococcus spp. Are Promising Candidates for Targeted Antibacterial Therapy" Viruses 14, no. 4: 831. https://doi.org/10.3390/v14040831

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

Article Metrics

Back to TopTop