Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Novel Drexlerviridae bacteriophage KMI8 with specific lytic activity against Klebsiella michiganensis and its biofilms

  • Heng Ku,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Department of Pharmacy and Biomedical Science, La Trobe Institute for Molecular Science, La Trobe University, Victoria, Australia

  • Mwila Kabwe,

    Roles Formal analysis, Investigation, Methodology, Visualization, Writing – review & editing

    Affiliation Department of Pharmacy and Biomedical Science, La Trobe Institute for Molecular Science, La Trobe University, Victoria, Australia

  • Hiu Tat Chan,

    Roles Conceptualization, Writing – review & editing

    Affiliations Department of Physiology, Anatomy and Microbiology, La Trobe University, Victoria, Australia, Department of Microbiology, Royal Melbourne Hospital, Victoria, Australia

  • Cassandra Stanton,

    Roles Investigation

    Affiliation Department of Physiology, Anatomy and Microbiology, La Trobe University, Victoria, Australia

  • Steve Petrovski,

    Roles Investigation, Supervision, Validation, Writing – review & editing

    Affiliation Department of Physiology, Anatomy and Microbiology, La Trobe University, Victoria, Australia

  • Steven Batinovic,

    Roles Data curation, Formal analysis, Investigation, Methodology, Writing – review & editing

    Affiliation Department of Physiology, Anatomy and Microbiology, La Trobe University, Victoria, Australia

  • Joseph Tucci

    Roles Conceptualization, Formal analysis, Investigation, Methodology, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    J.Tucci@latrobe.edu.au

    Affiliation Department of Pharmacy and Biomedical Science, La Trobe Institute for Molecular Science, La Trobe University, Victoria, Australia

Abstract

The bacterial genus Klebsiella includes the closely related species K. michiganensis, K. oxytoca and K. pneumoniae, which are capable of causing severe disease in humans. In this report we describe the isolation, genomic and functional characterisation of the lytic bacteriophage KMI8 specific for K. michiganensis. KMI8 belongs to the family Drexlerviridae, and has a novel genome which shares very little homology (71.89% identity over a query cover of only 8%) with that of its closest related bacteriophages (Klebsiella bacteriophage LF20 (MW417503.1); Klebsiella bacteriophage 066039 (MW042802.1). KMI8, which possess a putative endosialidase (depolymerase) enzyme, was shown to be capable of degrading mono-biofilms of a strain of K. michiganensis that carried the polysaccharide capsule KL70 locus. This is the first report of a lytic bacteriophage for K. michiganensis, which is capable of breaking down a biofilm of this species.

Introduction

The closely related Klebsiella bacteria K. pneumoniae and K. oxytoca form part of the normal human flora, but are also known to cause serious infections. K. pneumoniae, a member of the group of ESKAPE pathogens (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa and Enterobacter spp), can cause pneumonia, urinary tract infections, and wound infections [13]. K. oxytoca is implicated in up to half of the cases of antibiotic associated haemorrhagic colitis and antibiotic associated diarrhoea [4, 5]. K. oxytoca has also been shown to contribute to the severe weight loss in cancer anorexia/cachexia, as its prevalence is increased, and it may alter gut barrier function, in this condition [6]. In 2013, a novel species of the genus Klebsiella, K. michiganensis, was first identified and characterised [7]. Its closest phylogenetic relative is K. oxytoca, with which it shares 99% nucleotide sequence identity in the 16S rRNA gene sequence [7]. Although there have only been limited studies on K. michiganensis to date, it is being recognised as an emerging pathogen [8]. It has been associated with sepsis in immunocompromised patients [9], and nosocomial infections [10], including in neonatal units [11] where it can cause necrotizing enterocolitis, a potentially fatal disease affecting the gut of premature infants [8]. In this connection, it has been suggested to be more clinically relevant to a subset of preterm infants than K. oxytoca [8]. Further, as K. michiganensis is often misidentified as K. oxytoca, despite clinical laboratories deploying the latest technology available to them [11], some of the conditions attributed to K. oxytoca may be caused by K. michiganensis instead. Similar to K. oxytoca and K. pneumoniae, K. michiganensis has been shown to carry a range of resistance mechanisms, including genomic and plasmid mediated extended spectrum β-lactamases [8, 11], carbapenemases [9, 10, 12], and multi drug resistance efflux pumps [10]. Some of these resistance factors may have been acquired through horizontal gene transfer with other enteric bacterial species [10].

Bacterial resistance to antibiotics, which was warned about by Alexander Fleming in his Nobel Prize speech in 1945 [13], is now one of the biggest problems facing medicine. Another factor which increases bacterial virulence is the capacity to form biofilms. While both K. pneumoniae and K. oxytoca are able to form microbial biofilms [1419], the capacity of K. michiganensis to do so has not been studied. The treatment of bacteria in these communities is hindered by the fact that antibiotics do not readily penetrate biofilms [20].

Bacteriophages may offer the hope of adjunct therapy for infectious diseases, and provide options in such treatment. These bacterial viruses are specific in their actions. They offer a useful tool in the selective manipulation of the microbiome, with minimal disruption to the non-targeted microflora, in a manner that antibiotics cannot achieve. They are also able to disrupt bacterial biofilms and kill antibiotic resistant bacteria [2124]. In the past decade or so, medical research in Western countries has begun to focus on the potential for application of bacteriophages, with successful administration in clinical scenarios having been reported [2527]. In some instances, not only was the clinical response noted, but also total eradication of the pathogen, a scenario which is not seen with conventional antibiotic treatment [28]. That bacteriophages administered intravenously in the treatment of severe bacterial infections are safe and well tolerated, has also been reported [29].

We and others have previously reported the isolation and characterisation of bacteriophages lytic against K. oxytoca and K. pneumoniae [3032]. While the host range of most bacteriophages is limited to strains within a single species, there have been reports of bacteriophages which are able to kill both K. oxytoca and K. pneumoniae [30, 31, 33]. This extended host range killing across species is most likely due to the very close homology between K. oxytoca and K. pneumoniae. K. pneumoniae associated capsules have been found among K. oxytoca isolates [34, 35], which may help to explain why some bacteriophages are able to lyse both strains. The initial report of the characterisation of K. michiganensis also described it as containing a capsule [7].

We report here the isolation, genomic and functional characterisation of KMI8, a novel lytic bacteriophage which is specific for K. michiganensis. KMI8 shares very little genetic similarity with other Klebsiella bacteriophages, and is shown here to be capable of degrading K. michiganensis biofilms in-vitro. This is the first report of a bacteriophage specific for K. michiganensis and which is capable of degrading its biofilms. The potential exists for further testing of bacteriophages such as KMI8 in the treatment of biofilms and infections of this bacteria.

Methods

Ethics approval

All methods were performed in accordance with the La Trobe University Ethics, Biosafety and Integrity guidelines and regulations. Clinical isolates of K. michiganensis, K. oxytoca and K. pneumoniae were obtained from specimen cultures as part of routine care. All samples were obtained with oral consent, de-identified and handled according to the La Trobe University Ethics committee. No ethical concerns were raised as human material was not the focus of the study and bacteria obtained were not traced back to the individual. The study protocols were approved by the La Trobe University Ethics Committee, reference number: S17-111.

Bacterial cultures and identification

The 11 Klebsiella strains used here (Table 1) were grown using tryptone soy broth (TSB) and 1.2% agar (Oxoid, Adelaide, Australia), incubated at 37ºC aerobically for two days. Preliminary identification of the strains was performed by matrix assisted laser desorption ionisation-time of flight mass spectrometry (MALDI-TOF, Bruker Daltonik, Germany), and 16S rRNA gene sequencing. For extraction of bacterial DNA, a loopful of colonies was resuspended in 50 μL of nuclease free water in a 1.5 mL microcentrifuge tube and briefly vortexed. To this was added 2 μL of 0.5 M ethylenediaminetetraacetic acid (EDTA) (pH 8.0), 2.5 μL of 10% (w/v) sodium dodecyl sulphate (SDS) and 2.5 μL of proteinase K (20 mg/mL). The contents were gently mixed by inverting and incubated at 55 °C for one h. After cooling to room temperature, 150 μL of nuclease free water and 200 μL of phenol: chloroform: isoamyl alcohol (29:28:1) were added and mixed by inverting 10 times and vortexed until cloudy. The cloudy suspension was centrifuged at 12000 x g for 10 min. The clear upper aqueous layer was pipetted into a fresh 1.5 mL microcentrifuge tube containing an equal volume of isopropanol and incubated overnight at -20 °C. The tube was then centrifuged for 10 min at 12000 x g. The supernatant was discarded and 200 μL of 70% ethanol added to the DNA pellet. This was centrifuged for 5 min at 12000 x g, then the ethanol discarded. The DNA pellet was air dried before resuspending in 50 μL of nuclease free water. PCR amplification of the 16S rRNA gene of each bacterial strain was performed using universal primers 27F (5’-AGAGTTTGATCMTGGCTCAG-3’), 519R (5’-GWATTACCGCGGCKGCTG-3’), 895F (5’-RCCTGGGGAGTRCRG-3) and 1492R (5’-GGTTACCTTGTTACGACTT-3’) [32, 3638]. Go Taq® Long PCR Master Mix (Promega, Sydney, Australia) was used for the PCR amplification. The samples were subjected to an initial denaturation step of 2 min at 95°C, followed by 30 amplification cycles (94°C, 20 s; 59°C, 20 s; 72°C, 30 s) and a final elongation step of 10 min at 72°C. PCR products were analysed by electrophoresis in 1.2% agarose gels, stained with ethidium bromide and visualised under ultraviolet light [39]. Products were cleaned with an Ultra Clean® DNA purification kit (MO BIO San Diego), before submission to the Australian Genome Research Facility for Sanger sequencing.

thumbnail
Table 1. Klebsiella spp. used for KMI8 isolation and host range studies.

https://doi.org/10.1371/journal.pone.0257102.t001

The 11 Klebsiella strains were then definitively identified by whole genome sequencing (WGS). DNA extracted from the Klebsiella bacterial strains was sequenced using an Illumina MiSeq® next generation sequencing platform. The DNA library was prepared by using the NEBNext® Ultra II DNA Library Prep Kit according to the manufacturer’s instructions and sequenced using a MiSeq® V2 300-cycle reagent or V3 600-cycle kit to generate 150 bp or 300 bp paired end reads. Short-read data were filtered using Trim Galore v0.6.4 with the default settings (Q scores of ≥20, with automatic adapter detection). Reads were assembled de novo using Unicycler v0.4.8. Open reading frames of these sequences were predicted using Prokka 1.14.5. Bacterial taxonomy was assessed based on average nucleotide identity (ANI) using JSpeciesW [40] and the GTDBK-Tk software toolkit for assigning objective taxonomic classifications to prokaryotic genomes. These results were confirmed by ANI analysis results from NCBI. The program Staramr V.0.7.2 (Galaxy Australia) and the databases ResFinder, PlasmidFinder, and PointFinder were used to identify antibiotic resistance genes in the genomes of KLEB011, KLEB014 & KLEB015. We employed the program CRISPRcasFinder to identify any potential CRISPR regions in the genomes of the Klebsiella strains used here. The K. michiganensis KLEB011 genome was analysed using the Kaptive Web online tool for the rapid typing of Klebsiella K (polysaccharide capsule) and O (lipopolysaccharide) loci [41].

Isolation and purification of bacteriophages

Wastewater from Victoria, Australia, was used to screen for bacteriophages. Briefly, bacterial lawns of K. michiganensis, K. pneumoniae and K. oxytoca on Tryptone soy agar (TSA) (Oxoid, Australia) were prepared by using a cotton swab to spread onto the 1.2% agar plates a sample from a log phase culture. Bacteria and wastewater samples were incubated for 2 d in TSB (Oxoid, Adelaide, Australia) at 37 ºC, then centrifuged (12,000 × g for 5 min) and filtered through 0.2 μM pore cellulose acetate filters (Advantec, Melbourne, Australia). The filtrate (10 μL) was placed onto the bacterial lawn and plates were incubated for 24 h. Any potential bacteriophage clearing was excised (along with a portion of agar) and resuspended in 500 μL of TSB in a 1.5 mL microcentrifuge tube, before centrifugation (12,000× g for 5 min) and a 10-fold serial dilution was completed. A total of 10 μL of each dilution was placed on a bacterial lawn to observe plaques. This serial dilution purification was repeated five times to ensure single virion infection. After isolation and purification of bacteriophage, an enrichment method was used to prepare high concentration phage stocks. 300μL of bacteriophage solution was incubated with 500μL of bacteria (log phase) in a final volume of 10mL TSB at 37°C. This process was repeated several times using the previous concentrated bacteriophage stocks until the bacteriophage concentration of >107 PFU/mL was reached.

Extraction of bacteriophage DNA

Five mL of bacteriophage stock (>109 PFU/mL) was used to extract bacteriophage DNA. All chemicals were obtained from Sigma-Aldrich (Sydney, Australia) unless stated otherwise. The bacteriophage stock was digested at room temperature for 30 min in 5 mmol/L MgCl2, 10 μg/mL DNase I and RNase A to remove any bacterial nucleic acids present. The bacteriophages were recovered by overnight polyethylene glycol (PEG) precipitation using 10% (w/v) PEG 8000 and 1 mol/L NaCl at 4°C. Precipitated bacteriophages were centrifuged (12,000× g) for 15 min and resuspended in nuclease free water (Promega, Sydney, Australia). Bacteriophage proteins were digested for one h at 55°C by adding 50 μg/mL Proteinase K, 20 mmol/L EDTA and 0.5% (v/v) SDS. Following this, an equal volume of phenol–chloroform–isoamyl alcohol (29:28:1) was added and centrifuged (12,000× g) for 5 min. The aqueous phase was removed and an equal volume of isopropanol was added to it. The DNA was precipitated overnight at −20°C, then collected by centrifugation (12,000× g) for 10 min. The DNA pellet was washed with 70% ethanol, air dried and resuspended in 25 μL of nuclease free water (Promega).

Host range analysis

To assess for lytic capacity, bacteriophage stocks (>107 PFU/mL) were diluted 10-fold serially and tested on K. michiganensis, K. oxytoca and K. pneumoniae strains. A bacterial lawn of the strain was prepared and 10 μL of each dilution of the stock was placed on the lawn. If individual plaques were observed, it was noted that the bacteriophage could target the strain. The efficiency of plating was calculated as the PFU/ml of the bacteriophage on each strain divided by, and relative to, the PFU/ml obtained on the strain on which it most efficiently lysed, multiplied by 100.

Transmission electron microscopy

Transmission electron microscopy (TEM) was performed using 400-mesh formvar and carbon copper grids (ProSciTech, Townsville, Australia). Bacteriophage stocks (> 107 PFU/ mL) were allowed to adsorb to the grid for one min before excess solution was removed using filter paper. The grids were then negatively stained three times with 2% (w/v) uranyl acetate for 20 s. Excess stain was removed by filter paper and grids were air-dried for 20 min before examination under a JEOL JEM-2100 transmission electron microscope. This was operated at an accelerating voltage of 200 kV and high-resolution digital images were recorded on a Gatan Orius SC200D 1 wide angle camera with Gatan Microscopy Suite and Digital Micrograph (Version 2.32.888.0) imaging software. Virions were measured using ImageJ software [42] (Version 1.8.0_112).

Bacteriophage genomic characterisation and phylogeny

DNA extracted from the bacteriophage was sequenced using an Illumina MiSeq® next generation sequencing platform. The DNA library was prepared by using a Nextera® XT DNA sample preparation kit according to the manufacturer’s instructions and sequenced using a MiSeq® V2 reagent kit (300 cycles), as 150 bp paired end reads. Short-read data were filtered using Trim Galore v0.6.4 with the default settings (Q scores of ≥20, with automatic adapter detection). Reads were assembled de novo within Unicycler v0.4.8. Open reading frames of these sequences were predicted using Glimmer, Prodigal V.2.6.2 and manually checked by Geneious prime V.2019.0.3. The bacteriophage genome was annotated using Blastp (Version 2.9.0) nr database with a cut off e-value of (1e-5), NCBI Web CD-search Tool (Pfam database) with a cut off e-value of 0.01 and HHpred (default database:PDB_mmCIF70_23_Jul; Pfam-A_V33.1; NCBI_Conserved Domain (CD)_v3.18). The proteomic tree was constructed with ViPTreeGen (Version 1.1.2) using sequence similarity distance based on tBLASTx results and construct (bio)nj tree. EvolView was used to annotate the proteomic tree. tRNAs were predicted by ARAGORN and tRNAscan-SE-2.0.5.

One step growth analysis

The one step growth analysis was performed on three replicates as described previously [43]. K. michiganensis strain KLEB011 was grown to exponential phase in TSB, and harvested by centrifugation at 12,000× g for 5 min. 1x106 bacterial cells were resuspended in fresh 900 μL TSB to which 100 μL of KMI8 bacteriophage stock was added at a multiplicity of infection (MOI) of 1, and incubated at 4°C for 5 min to allow for adsorption. Adsorbed bacteriophages were collected by centrifugation at 12,000× g for 5 mins and resuspended in 50 mL of fresh TSB. The mixture was incubated aerobically at 37°C. One mL aliquots were taken every 10 min and bacteriophage concentration (PFU/mL) assayed by centrifugation at 12,000× g for 2 min before 10-fold serial dilution and placing 10 μL of these dilutions on a fresh lawn of K. michiganensis strain KLEB11. These were incubated aerobically overnight to observe bacteriophage plaque formation.

Biofilm growth and quantification

The biofilm experiment was conducted aerobically using TSA broth (Oxoid, Scoresby, Australia) supplemented with 0.5% glucose in 96 well polystyrene plates (Greiner bio-one, Frickenhausen, Germany) coated with 0.5% gelatin (Sigma-Aldrich, Sydney, Australia). The K. michiganensis strain KLEB11 was cultured overnight, and 100 μL of 108 CFU/mL of this bacteria in exponential growth phase was added to each well with an equivalent volume of fresh TSA broth. The 96 well plates were incubated shaking at 37°C under aerobic conditions for 3 d with sterile broth replenishment every 24 h. After the biofilm was formed, bacteria in suspension were removed before 100 μL of bacteriophage KMI8 at 108 PFU/mL was added to each well and incubated for a further 2.5 h, before quantification assays completed as described previously [44]. Briefly, unattached K. michiganensis cells were washed off by submerging in Milli-Q® deionised water (Merck, Bayswater, Australia) for 5 min and air dried before staining with 200 μL of 0.1% (v/v) crystal violet for 10 min. Excess crystal violet stain was washed off twice by resuspension in Milli-Q® deionised water (Merck, Bayswater, Australia) for another 5 min before leaching out the cell bound crystal violet in 70% (v/v) ethanol. Absorbance in each well was determined using a FlexStation 3 plate reader (Molecular Devices, USA) at a wavelength of 600 nm.

Biofilm viability analysis

Biofilm viability was determined as previously described [44]. Briefly, the K. michiganensis biofilm was grown on gelatin-coated microscope coverslips and treated with bacteriophage KMI8 in the same manner as the 96 well plates described above. SYBR® gold and propidium iodide (PI) were used to stain nucleic acids of all cells (including live membrane intact and dead, membrane compromised cells). 2 μL of a 1 mg/mL concentration of PI in dimethyl sulfoxide (Sigma-Aldrich, Australia) was added to 0.2 μL of a 1 mg/mL concentration of SYBR® gold and resuspended in 100 μL of nuclease free water (Promega, Australia). The mixture was applied to the biofilm on slides and incubated for 15 min in the dark. Excess stain was then washed off in nuclease free water (Promega, Sydney, Australia) and the coverslips mounted on microscope slides using vectorshield® (Burlingame, CA, USA). Slides were immediately examined under the Fluoview Fv10i-Confocal laser-scanning microscope (Olympus Life Science, Notting Hill, Australia) at excitation wavelength of 485 nm and emissions at 535 nm and 635 nm for green and red fluorescence, respectively. SYBR® gold stains all cells green while PI stains all membrane compromised cells red. When these stains were combined and cells imaged as in these experiments, cells which were stained green indicated live cells, and those not green, (red or orange) indicated membrane compromised dead cells. A quantitative analysis of the viable bacterial cells, with and without bacteriophage treatment, was performed. For these analyses, we scanned and counted cells in 14 different areas of the coverslips for bacteriophage treated and untreated biofilms, across triplicate experiments.

Statistical analysis

The absorbance values quantifying the biofilms were analysed for normality using the Shapiro Wilk test and since all data were normal, their means compared by a paired T-test. For the biofilm viability examination using confocal microscopy, a Shapiro-Wilk test showed that some data were normally distributed while some were not normally distributed. As such, a non-parametric test (2-independent samples test), the Mann-Whitney U test, was used to analyse the data. A p-value of less than 0.05 was considered statistically significant. All statistical analysis was performed using the Statistical Package for Social Sciences (SPSS version 25).

Codon usage analysis

Geneious prime V.2020.0.3 was employed to extract nucleotide sequences from genome annotations of the K. michiganensis strain KLEB011 and KMI8. Biopython script was written and used to remove the start codons by slicing nucleotide sequence [45, 46]. K. michiganensis KLEB011 and KMI8 codon usage was calculated using python script and findings confirmed with the Sequence Manipulation Suite of JavaScript programs for analyzing and formatting protein and DNA sequences [47].

Codon usage analyses to investigate the adaptive role of bacteriophage encoded tRNAs have been described previously [48]. Briefly, for the 64 codons, a relative codon frequency, fi, was calculated using the following formula: fi = Number of codoni/ Number of all codons. The ratio of relative codon frequencies, ri, between bacteriophage and host was then assessed as: ri = fi bacteriophage/fihost. If ri was ≥ 1.1, the bacteriophage was considered to exhibit a higher relative codon frequency for codon i than the host; If ri was between 0.9 and 1.1, the bacteriophage and host were considered to have similar relative codon frequencies; if ri was < 0.9, the host was considered to exhibit a higher relative codon frequency [48].

Results

Definitive identification of Klebsiella strains for use in KMI8 isolation

Klebsiella bacterial strains used in this study were identified and differentiated using a combination of MALDI-TOF, 16S rRNA gene sequencing, and WGS (Table 1). The use of MALDI-TOF and 16S rRNA gene sequencing were successful in differentiating K. pneumoniae from the other species used in this study (Table 1). However, they were unsuccessful in differentiating K. oxytoca from K. michiganensis: these techniques incorrectly identified each of KLEB11-KLEB16 as K. oxytoca. In contrast, WGS analysis identified only KLEB12 as K. oxytoca, while KLEB11 and KLEB13-KLEB16 were identified as K. michiganensis (Table 1). Using a program employed by NCBI, JSpeciesWs, we identified the % genomic similarity of the strains used here. The program uses Average Nucleotide Identity (ANI) calculations based on BLAST+ to calculate the % genomic similarity between these strains. S1 Table shows the average nucleotide identity between the Klebsiella strains used in this study.

Identification of antibiotic resistance genes in the genomes of KLEB011, KLEB014 & KLEB015, the K. michiganensis strains lysed by KMI8

The program Staramr V.0.7.2 (Galaxy Australia) and databases ResFinder, PlasmidFinder, and PointFinder were used to identify antibiotic resistance genes in the genomes of KLEB011, KLEB014 & KLEB015 (Table 2).

thumbnail
Table 2. Antibiotic resistance genes in K. michiganensis strains KLEB011, KLEB014 & KLEB015.

https://doi.org/10.1371/journal.pone.0257102.t002

Capsule locus identification for K. michiganensis strain KLEB011 used in one step growth curve and biofilm assays

Use of the Kaptive Web tool for the rapid typing of Klebsiella K (polysaccharide capsule) and O (lipopolysaccharide) loci [41] revealed that the KLEB011 genome carried the capsule KL70 locus.

Isolation, host range analysis, and phenotypic characterization of novel K. michiganensis bacteriophage KMI8

KMI8, a novel bacteriophage isolated from samples of wastewater from Victoria, Australia, was capable of lysing three K. michiganensis strains (KLEB011, KLEB014 & KLEB015) but did not kill any of the K. oxytoca or K. pneumoniae strains tested. We employed the program CRISPRcasFinder to identify any potential CRISPR regions in the genomes of the Klebsiella strains used here. These genomes appear to be devoid of any functional anti-KMI8 CRISPR or any evidence of sufficient Cas genes for a functional CRISPR-Cas system. The efficiency of plating for KMI8 on its host K. michiganensis strains was KLEB011: 100%; KLEB014: 3.2%; KLEB015: 0.48% (calculated as the PFU/ml of KMI8 on each strain divided by, and relative to, the PFU/ml obtained on strain KLEB011, on which it most efficiently lysed, multiplied by 100).

Plaques formed by KMI8 on 1.2% agar were approximately 1 mm in diameter. TEM revealed that KMI8 was of the siphovirus morphotype (Fig 1). KMI8 had a capsid length of approximately 71 nm ± 2 nm; tail length and width approximately 190 nm ± 5 nm and 10 ± 0.1 nm, respectively. In a study of the growth kinetics of KMI8 grown on KLEB11, one step growth analysis showed latent and burst periods of 20 min and 30 min respectively (Fig 2), and a burst size of 12 ± 2 PFU/infected bacteria.

thumbnail
Fig 1. TEM characterisation of bacteriophage KMI8: Siphovirus with capsid length of 71nm ± 2nm; tail length and width approximately 190 nm ± 5 nm and 10 ± 0.1 nm, respectively.

https://doi.org/10.1371/journal.pone.0257102.g001

thumbnail
Fig 2. KMI8 one step growth curve.

Error bars represent standard error of mean, calculated from three independent experiments.

https://doi.org/10.1371/journal.pone.0257102.g002

Genome analysis of K. michiganensis bacteriophage KMI8

The genome sequence data was assembled de novo. The KMI8 bacteriophage genome (Accession number: MN101222) had a genome size of 52904 bp (with an average coverage of 474x) and was composed of 78 predicted open reading frames (ORFs) (Table 3 and Fig 3) of which 26.9% (21/78) were related to genes encoding putative functional proteins (Table 4). The GC content of KMI8 was 49.3% (Table 3) (the GC content of one of its hosts, the K. michiganensis strain KLEB11, was 55.7%) and the closest related organism is Klebsiella bacteriophage 066039 (MW042802) with 71.89% identity over a query cover of only 8%. No putative integrase genes, toxin genes, tRNA or Transfer-messenger RNA (tmRNA) genes were detected in the genome of KMI8. We assessed the bacteria it lyses, and did not detect the KMI8 genome in the genome of its host bacteria, nor in extracts of bacterial cells.

thumbnail
Fig 3. KMI8 (MN101222) genome.

Arrows indicate putative ORFs (numbered) and direction of transcription. Grey arrows indicate putative hypothetical genes; green arrows indicate putative terminase genes; blue arrow indicates putative capsid gene; yellow arrows indicate putative tail fibre genes; purple arrows indicate putative DNA manipulation and metabolism genes; red arrows indicate putative lysis genes. Scale bar = 1000bp.

https://doi.org/10.1371/journal.pone.0257102.g003

thumbnail
Table 4. KMI8 ORFs where predicted protein function is known.

Predicted function based on Pfam database of conserved functional domains.

https://doi.org/10.1371/journal.pone.0257102.t004

Disruption of K. michiganensis biofilm by bacteriophage KMI8

Bacteriophage KMI8 significantly lowered the K. michiganensis KLEB11 biofilm mass after 2.5 h exposure (p<0.001) (Fig 4). Confocal microscope imaging was used to visualise biofilm mass on glass slides. Control K. michiganensis biofilm mass (Fig 5B) and biofilm exposed to KMI8 for 2.5 h (Fig 5A) were stained with SYBR® Gold and PI to indicate membrane intact live cells (green colour) and membrane compromised dead cells (red/orange colour). Biofilm mass showed many live cells attached to the glass slide in the untreated control (Fig 5B) while there was a sparse population of membrane compromised/dead bacteria in the KMI8 treated biofilm (Fig 5A). Counting of the viable cells in confocal images of K. michiganensis biofilms before and after KMI8 treatment revealed that there was a significant difference between untreated and treated biofilms (P < 0.001) (Mann-Whitney U test) (Fig 6).

thumbnail
Fig 4. K. michiganensis biofilm absorbance measurements, showing significant reduction of K. michiganensis biofilm treated with KMI8 bacteriophage (p < 0.001).

https://doi.org/10.1371/journal.pone.0257102.g004

thumbnail
Fig 5. Confocal images of SYBR® gold and PI staining following (A) bacteriophage treated and (B) untreated K. michiganensis biofilm.

Green color indicates live cells, red/orange color indicates dead cells on the slides.

https://doi.org/10.1371/journal.pone.0257102.g005

thumbnail
Fig 6. Viable cells in K. michiganensis biofilms before and after KMI8 treatment.

There was a significant difference between untreated and treated biofilms (P < 0.001) (Mann-Whitney U test).

https://doi.org/10.1371/journal.pone.0257102.g006

Bacteriophage phylogeny

ViPTree online was used to generate a proteomic tree based on genome wide sequence similarities, computed by tBLASTx, for KMI8 and 59 related bacteriophages. KMI8, which belongs to the family Drexlerviridae, is part of a monophyletic clade with other Klebsiella bacteriophages of the family Drexlerviridae, genus Weberviridae (Fig 7). It is unclear, however, to which genus KMI8 can be classified at this stage.

thumbnail
Fig 7. Proteomic tree showing the genome-wide proteomic diversity of KMI8 and related bacteriophages.

Black box indicates KMI8.

https://doi.org/10.1371/journal.pone.0257102.g007

Codon usage analysis in KMI8 and K. michiganensis strain KLEB11

Analysis of the relative codon frequencies for KMI8 and K. michiganensis strain KLEB11 revealed that the frequency of usage of five codons was at least two and a half times higher in KMI8, compared to the host K. michiganensis strain KLEB11 (Fig 8). The ratio of usage of ACT (threonine) by KMI8, compared to the host bacteria, was over three times greater.

thumbnail
Fig 8. The ratio of relative codon frequencies between bacteriophage KMI8 and K. michiganensis strain KLEB011.

The figure shows codons for which KMI8 exhibited the 10 highest (blue) and 10 lowest (orange) relative codon frequencies compared to the host.

https://doi.org/10.1371/journal.pone.0257102.g008

Discussion

We describe here a novel lytic bacteriophage, KMI8, which has very little genetic similarity to other reported bacteriophages. Compared to its closest relative, the Klebsiella bacteriophage 066039 (accession MW042802), its genome has 71.89% identity over a query coverage of only 8%. KMI8 appears to be specific for K. michiganensis. It is capable of lysing several strains of this bacteria, which carry resistance genes for ampicillin and kanamycin, but its host range does not extend to the K. pneumoniae or K. oxytoca strains tested here. In connection to this, our study showed that techniques such as MALDI-TOF and 16S rRNA gene sequencing were not sufficient to differentiate between K. oxytoca and K. michiganensis, and as shown here and by others, the precise definition of the phylogeny of these Klebsiella species requires more extensive applications such as WGS [8, 10, 11]. This is the first report of a bacteriophage lytic for K. michiganensis. However, it is possible that if WGS were used to refine the taxonomic analysis of closely related hosts of previously reported Klebsiella bacteriophages, some of these hosts would be reclassified as K. michiganensis. Further, given the genetic similarity between K. michiganensis and K. oxytoca, it is possible that known K. oxytoca bacteriophages may also target K. michiganensis.

KMI8 has a siphovirus morphology, with a long flexible non-contractile tail. It is a member of the Drexlerviridae family of viruses, and forms a monophyletic group with other bacteriophages lytic for Klebsiella, and which are all Drexlerviridae viruses of the genus Webervirus. It is unclear, however, to which genus KMI8 can be classified at this stage. Analysis of the growth kinetics of KMI8 showed a long latent period (20 min), and a rise period of up to 30 min during which the virus was bursting from the host cell, followed by a plateau. The burst size was 12 ± 2 PFU/infected bacteria. This compares to latency of 30 min and burst of 113 PFU/infected bacteria for the most closely related bacteriophage whose growth kinetics have been reported, the Drexlerviridae Webervirus TSK1 (accession MH688453) [49], which lyses K. pneumoniae. The GC content of the KMI8 genome was 49.3%, compared to that of one of its hosts, the K. michiganensis strain KLEB11, which was 55.7%. The putative function of most of the genome of KMI8 is unknown. We were only able to ascribe functionality for 21 of the 78 open reading frames it contains, based on the Pfam database of conserved functional domains. Its genes whose functions were predicted included putative terminase genes, putative capsid genes, putative tail fibre genes, putative DNA manipulation and metabolism genes, and putative lysis genes.

Analysis of the relative codon frequencies for KMI8 and the K. michiganensis strain KLEB11 revealed that the frequency of usage of a range of codons was higher in the bacteriophage, compared to the host. Previous suggestions have been for an adaptive role of bacteriophage encoded tRNAs, in that they are selected in bacteriophages to compensate for differences in codon usage between the bacteriophage and the host [50]. However, even though the ratio of usage of codon ACT (threonine) by KMI8 compared to the host bacteria was over three times greater, KMI8 possessed no tRNA or tmRNA genes for the anticodon corresponding to this or other codons. The reasons for this are unclear.

Klebsiella bacteria are known to form biofilms, which increases their capacity for virulence. KMI8 was shown to be capable of degrading K. michiganensis mono biofilms, and this is the first report of bacteriophages degrading biofilms of this host. It is known that bacteria, including Klebsiella spp, can develop bacteriophage resistance via structural barriers such as capsular polysaccharides [51], and that capsule degrading depolymerase enzymes are important in allowing bacteriophages to break down K. pneumoniae biofilms [52]. The K. michiganensis strain KLEB011 which was used for the biofilm studies in our work, possessed the polysaccharide capsule KL70 locus, which is also found in K. pneumoniae [35, 41]. The genome of KMI8 possessed a putative endosialidase gene element (according to the Pfam database), embedded in an ORF which bears similarity to tail fibre elements (ORF 25). Endosialidases act as depolymerases, and are known to degrade the bacterial capsule polysaccharide polysialic acid [5355]. The capacity to degrade K. michiganensis biofilms is potentially of clinical importance. Gut microbial dysbiosis is known to contribute to carcinogenesis, and while organisms such as Fusobacterium nucleatum have been shown to promote cancer progression by forming biofilms around tumours [5658], K. michiganensis has not been directly implicated. However, extensive studies on the pathology caused by K. michiganensis are limited. Those that have been performed show that it is important clinically, as a nosocomial pathogen able to cause blood borne infection in immunocompromised patients, and severe, possibly fatal gut damage in preterm infants [811]. Apart from their incapacity to penetrate bacterial biofilms [20], antibiotics may increase colonisation and biofilm formation by Klebsiella spp [18]. Recent reports of the purification of enzymes capable of interfering with quorum sensing (involved in cell-cell communication in bacterial biofilms) offer new prospects for disruption/prevention of biofilm formation [59]. The use of bacteriophages similar to KMI8 which are specific for K. michiganensis also offer potential, especially in specificity in rectifying dysbiosis in the gut and invasion of other tissues. To this end, bacteriophage cocktails have been shown to be more effective than single bacteriophages in disrupting biofilms in vitro as well as in prevention and elimination of biofilm related infections [21]. As such, further analysis of bacteriophages which lyse K. michiganensis, and their capacity to be applied in cocktails to disrupt biofilms, is warranted.

Finally, as mentioned above, extensive studies on the pathology caused by K. michiganensis are limited, and the differentiation between this species and K. oxytoca requires genomic analyses such as WGS. This has important implications for Klebsiella speciation in diagnostic laboratories, as such molecular techniques are not readily available in routine clinical laboratories. Therefore, the occurrence and significance of K. michiganensis in clinical disease is likely to have been under-estimated, and this newly identified species may contribute to diseases previously reported to be caused by K. oxytoca and other members of the genus. This also highlights the importance of WGS of Klebsiella spp. when considering establishment of bacteriophage banks for integration in clinical practice [60]. Given the tendency to misidentify K. michiganensis as K. oxytoca in clinical settings, bacteriophages against K. michiganensis, such as KMI8, may be considered for inclusion in bacteriophage cocktails for future treatment of clinically diagnosed K. oxytoca infections.

References

  1. 1. Podschun R, Ullmann U. Klebsiella spp. as Nosocomial Pathogens: Epidemiology, Taxonomy, Typing Methods, and Pathogenicity Factors. 1998;11(4):589–603.
  2. 2. Paczosa MK, Mecsas J. Klebsiella pneumoniae: Going on the Offense with a Strong Defense 2016;80(3):629–61. pmid:27307579
  3. 3. Martin RM, Bachman MA. Colonization, Infection, and the Accessory Genome of Klebsiella pneumoniae. Frontiers in cellular and infection microbiology. 2018;8:4. pmid:29404282
  4. 4. Högenauer C, Langner C, Beubler E, Lippe IT, Schicho R, Gorkiewicz G, et al. Klebsiella oxytoca as a causative organism of antibiotic-associated hemorrhagic colitis. The New England journal of medicine. 2006;355(23):2418–26. pmid:17151365
  5. 5. Yilmaz M, Bilir YA, Aygün G, Erzin Y, Ozturk R, Celik AF. Prospective observational study on antibiotic-associated bloody diarrhea: report of 21 cases with a long-term follow-up from Turkey. European journal of gastroenterology & hepatology. 2012;24(6):688–94. pmid:22433794
  6. 6. Pötgens SA, Brossel H, Sboarina M, Catry E, Cani PD, Neyrinck AM, et al. Klebsiella oxytoca expands in cancer cachexia and acts as a gut pathobiont contributing to intestinal dysfunction. Scientific Reports. 2018;8(1):12321. pmid:30120320
  7. 7. Saha R, Farrance CE, Verghese B, Hong S, Donofrio RS. Klebsiella michiganensis sp. nov., a new bacterium isolated from a tooth brush holder. Current microbiology. 2013;66(1):72–8. pmid:23053492
  8. 8. Chen Y, Brook TC, Soe CZ, apos, Neill I, Alcon-Giner C, et al. Preterm infants harbour diverse Klebsiella populations, including atypical species that encode and produce an array of antimicrobial resistance- and virulence-associated factors. 2020;6(6).
  9. 9. Seiffert SN, Wüthrich D, Gerth Y, Egli A, Kohler P, Nolte O. First clinical case of KPC-3-producing Klebsiella michiganensis in Europe. New microbes and new infections. 2019;29:100516. pmid:30949345
  10. 10. Founou RC, Founou LL, Allam M, Ismail A, Essack SY. Genomic characterisation of Klebsiella michiganensis co-producing OXA-181 and NDM-1 carbapenemases isolated from a cancer patient in uMgungundlovu District, KwaZulu-Natal Province, South Africa. South African medical journal = Suid-Afrikaanse tydskrif vir geneeskunde. 2018;109(1):7–8. pmid:30606295
  11. 11. Chapman P, Forde BM, Roberts LW, Bergh H, Vesey D, Jennison AV, et al. Genomic Investigation Reveals Contaminated Detergent as the Source of an Extended-Spectrum-β-Lactamase-Producing Klebsiella michiganensis Outbreak in a Neonatal Unit. 2020;58(5):e01980–19. pmid:32102855
  12. 12. Pedersen T, Sekyere JO, Govinden U, Moodley K, Sivertsen A, Samuelsen Ø, et al. Spread of Plasmid-Encoded NDM-1 and GES-5 Carbapenemases among Extensively Drug-Resistant and Pandrug-Resistant Clinical Enterobacteriaceae in Durban, South Africa. Antimicrobial agents and chemotherapy. 2018;62(5). pmid:29507063
  13. 13. Fleming A. Penicillin. 1945.
  14. 14. Burgos-Garay M, Ganim C, de Man TJB, Davy T, Mathers AJ, Kotay S, et al. Colonization of carbapenem-resistant Klebsiella pneumoniae in a sink-drain model biofilm system. Infection control and hospital epidemiology. 2020:1–9.
  15. 15. Tan S, Gao J, Li Q, Guo T, Dong X, Bai X, et al. Synergistic effect of chlorogenic acid and levofloxacin against Klebsiella pneumonia infection in vitro and in vivo. Scientific Reports. 2020;10(1):20013. pmid:33203903
  16. 16. Yang K, Liu S, Li H, Du N, Yao J, He Q, et al. NDM-1-Positive K. pneumoniae at a Teaching Hospital in Southwestern China: Clinical Characteristics, Antimicrobial Resistance, Molecular Characterization, Biofilm Assay, and Virulence. The Canadian journal of infectious diseases & medical microbiology = Journal canadien des maladies infectieuses et de la microbiologie medicale. 2020;2020:9091360.
  17. 17. Amaretti A, Righini L, Candeliere F, Musmeci E, Bonvicini F, Gentilomi GA, et al. Antibiotic Resistance, Virulence Factors, Phenotyping, and Genotyping of Non-Escherichia coli Enterobacterales from the Gut Microbiota of Healthy Subjects. International journal of molecular sciences. 2020;21(5).
  18. 18. Alkhudhairy MK, Alshadeedi SMJ, Mahmood SS, Al-Bustan SA, Ghasemian A. Comparison of adhesin genes expression among Klebsiella oxytoca ESBL-non-producers in planktonic and biofilm mode of growth, and imipenem sublethal exposure. Microbial pathogenesis. 2019;134:103558. pmid:31136790
  19. 19. Ghasemian A, Mohabati Mobarez A, Najar Peerayeh S, Talebi Bezmin Abadi A, Khodaparast S, Mahmood SS. Expression of adhesin genes and biofilm formation among Klebsiella oxytoca clinical isolates from patients with antibiotic-associated haemorrhagic colitis. Journal of medical microbiology. 2019;68(7):978–85. pmid:31136296
  20. 20. Høiby N, Bjarnsholt T, Givskov M, Molin S, Ciofu O. Antibiotic resistance of bacterial biofilms. International journal of antimicrobial agents. 2010;35(4):322–32. pmid:20149602
  21. 21. Łusiak-Szelachowska M, Weber-Dąbrowska B, Górski A. Bacteriophages and Lysins in Biofilm Control. Virologica Sinica. 2020;35(2):125–33. pmid:32125643
  22. 22. Hansen MF, Svenningsen SL, Røder HL, Middelboe M, Burmølle M. Big Impact of the Tiny: Bacteriophage-Bacteria Interactions in Biofilms. Trends in microbiology. 2019;27(9):739–52. pmid:31128928
  23. 23. Shlezinger M, Khalifa L, Houri-Haddad Y, Coppenhagen-Glazer S, Resch G, Que YA, et al. Phage Therapy: A New Horizon in the Antibacterial Treatment of Oral Pathogens. Current topics in medicinal chemistry. 2017;17(10):1199–211. pmid:27770768
  24. 24. Jończyk-Matysiak E, Weber-Dąbrowska B, Owczarek B, Międzybrodzki R, Łusiak-Szelachowska M, Łodej N, et al. Phage-Phagocyte Interactions and Their Implications for Phage Application as Therapeutics. Viruses. 2017;9(6). pmid:28613272
  25. 25. Dedrick RM, Guerrero-Bustamante CA, Garlena RA, Russell DA, Ford K, Harris K, et al. Engineered bacteriophages for treatment of a patient with a disseminated drug-resistant Mycobacterium abscessus. Nature Medicine. 2019;25(5):730–3. pmid:31068712
  26. 26. Schooley RT, Biswas B, Gill JJ, Hernandez-Morales A, Lancaster J, Lessor L, et al. Development and Use of Personalized Bacteriophage-Based Therapeutic Cocktails To Treat a Patient with a Disseminated Resistant Acinetobacter baumannii Infection. 2017;61(10):e00954–17.
  27. 27. Maddocks S, Fabijan AP, Ho J, Lin RCY, Ben Zakour NL, Dugan C, et al. Bacteriophage Therapy of Ventilator-associated Pneumonia and Empyema Caused by Pseudomonas aeruginosa. American journal of respiratory and critical care medicine. 2019;200(9):1179–81. pmid:31437402
  28. 28. Wunderink RG. Turning the Phage on Treatment of Antimicrobial-Resistant Pneumonia. American journal of respiratory and critical care medicine. 2019;200(9):1081–2. pmid:31453719
  29. 29. Petrovic Fabijan A, Lin RCY, Ho J, Maddocks S, Ben Zakour NL, Iredell JR, et al. Safety of bacteriophage therapy in severe Staphylococcus aureus infection. Nature Microbiology. 2020;5(3):465–72. pmid:32066959
  30. 30. Kęsik-Szeloch A, Drulis-Kawa Z, Weber-Dąbrowska B, Kassner J, Majkowska-Skrobek G, Augustyniak D, et al. Characterising the biology of novel lytic bacteriophages infecting multidrug resistant Klebsiella pneumoniae. Virology Journal. 2013;10(1):100. pmid:23537199
  31. 31. Karumidze N, Kusradze I, Rigvava S, Goderdzishvili M, Rajakumar K, Alavidze Z. Isolation and characterisation of lytic bacteriophages of Klebsiella pneumoniae and Klebsiella oxytoca. Current microbiology. 2013;66(3):251–8. pmid:23143289
  32. 32. Brown TL, Petrovski S, Hoyle D, Chan HT, Lock P, Tucci J. Characterization and formulation into solid dosage forms of a novel bacteriophage lytic against Klebsiella oxytoca. PLOS ONE. 2017;12(8):e0183510. pmid:28817689
  33. 33. Ku H BT, Kabwe M, Petrovski S, Dyson ZA, Chan HT, Batinovic S & Tucci J Isolation and characterisation of novel Klebsiella bacteriophages: can comparative genomics shed light on the molecular mechanisms driving their host range specificity?. 2020.
  34. 34. Ishihara Y, Yagi T, Mochizuki M, Ohta M. Capsular types, virulence factors and DNA types of Klebsiella oxytoca strains isolated from blood and bile. Kansenshogaku zasshi The Journal of the Japanese Association for Infectious Diseases. 2012;86(2):121–6. pmid:22606879
  35. 35. Wyres KL, Wick RR, Gorrie C, Jenney A, Follador R, Thomson NR, et al. Identification of Klebsiella capsule synthesis loci from whole genome data. Microbial genomics. 2016;2(12):e000102. pmid:28348840
  36. 36. Hodkinson BP, Lutzoni F. A microbiotic survey of lichen-associated bacteria reveals a new lineage from the Rhizobiales. Symbiosis. 2009;49(3):163–80.
  37. 37. Turner S, Pryer KM, Miao VP, Palmer JD. Investigating deep phylogenetic relationships among cyanobacteria and plastids by small subunit rRNA sequence analysis. The Journal of eukaryotic microbiology. 1999;46(4):327–38. pmid:10461381
  38. 38. Schumann P. E. Stackebrandt and M. Goodfellow (Editors), Nucleic Acid Techniques in Bacterial Systematics (Modern Microbiological Methods). XXIX + 329 S., 46 Abb., 28 Tab. Chichester—New York—Brisbane—Toronto—Singapore 1991. John Wiley & Sons. $ 55.00. ISBN: 0–471–92906–9. 1991;31(6):479–80.
  39. 39. Kovtunovych G, Lytvynenko T, Negrutska V, Lar O, Brisse S, Kozyrovska N. Identification of Klebsiella oxytoca using a specific PCR assay targeting the polygalacturonase pehX gene. Research in Microbiology. 2003;154(8):587–92. pmid:14527660
  40. 40. Richter M, Rosselló-Móra R, Oliver Glöckner F, Peplies J. JSpeciesWS: a web server for prokaryotic species circumscription based on pairwise genome comparison. Bioinformatics (Oxford, England). 2016;32(6):929–31. pmid:26576653
  41. 41. Wick RR, Heinz E, Holt KE, Wyres KL. Kaptive Web: User-Friendly Capsule and Lipopolysaccharide Serotype Prediction for Klebsiella Genomes. 2018;56(6):e00197–18. pmid:29618504
  42. 42. Schneider CA, Rasband WS, Eliceiri KW. NIH Image to ImageJ: 25 years of image analysis. Nature methods. 2012;9(7):671–5. pmid:22930834
  43. 43. Chan HT, Ku H, Low YP, Brown T, Batinovic S, Kabwe M, et al. Characterization of Novel Lytic Bacteriophages of Achromobacter marplantensis Isolated from a Pneumonia Patient. Viruses. 2020;12(10). pmid:33049935
  44. 44. Kabwe M, Brown T, Speirs L, Ku H, Leach M, Chan HT, et al. Novel Bacteriophages Capable of Disrupting Biofilms From Clinical Strains of Aeromonas hydrophila. Frontiers in microbiology. 2020;11:194. pmid:32117183
  45. 45. Cock PJA, Antao T, Chang JT, Chapman BA, Cox CJ, Dalke A, et al. Biopython: freely available Python tools for computational molecular biology and bioinformatics. Bioinformatics (Oxford, England). 2009;25(11):1422–3. pmid:19304878
  46. 46. Chan HT, Ku H, Low YP, Batinovic S, Kabwe M, Petrovski S, et al. Characterization of Novel Lytic Bacteriophages of Achromobacter marplantensis Isolated from a Pneumonia Patient. 2020;12(10):1138.
  47. 47. Stothard P. The sequence manipulation suite: JavaScript programs for analyzing and formatting protein and DNA sequences. BioTechniques. 2000;28(6):1102, 4. pmid:10868275
  48. 48. Yoshikawa G, Askora A, Blanc-Mathieu R, Kawasaki T, Li Y, Nakano M, et al. Xanthomonas citri jumbo phage XacN1 exhibits a wide host range and high complement of tRNA genes. Scientific Reports. 2018;8(1):4486. pmid:29540765
  49. 49. Tabassum R, Shafique M, Khawaja KA, Alvi IA, Rehman Y, Sheik CS, et al. Complete genome analysis of a Siphoviridae phage TSK1 showing biofilm removal potential against Klebsiella pneumoniae. Scientific Reports. 2018;8(1):17904. pmid:30559386
  50. 50. Bailly-Bechet M, Vergassola M, Rocha E. Causes for the intriguing presence of tRNAs in phages. Genome Res. 2007;17(10):1486–95. pmid:17785533
  51. 51. Latka A, Maciejewska B, Majkowska-Skrobek G, Briers Y, Drulis-Kawa Z. Bacteriophage-encoded virion-associated enzymes to overcome the carbohydrate barriers during the infection process. Appl Microbiol Biotechnol. 2017;101(8):3103–19. pmid:28337580
  52. 52. Latka A, Drulis-Kawa Z. Advantages and limitations of microtiter biofilm assays in the model of antibiofilm activity of Klebsiella phage KP34 and its depolymerase. Scientific Reports. 2020;10(1):20338. pmid:33230270
  53. 53. Morley TJ, Willis LM, Whitfield C, Wakarchuk WW, Withers SG. A new sialidase mechanism: bacteriophage K1F endo-sialidase is an inverting glycosidase. The Journal of biological chemistry. 2009;284(26):17404–10. pmid:19411257
  54. 54. Jakobsson E, Jokilammi A, Aalto J, Ollikka P, Lehtonen JV, Hirvonen H, et al. Identification of amino acid residues at the active site of endosialidase that dissociate the polysialic acid binding and cleaving activities in Escherichia coli K1 bacteriophages. The Biochemical journal. 2007;405(3):465–72. pmid:17394421
  55. 55. Machida Y, Miyake K, Hattori K, Yamamoto S, Kawase M, Iijima S. Structure and function of a novel coliphage-associated sialidase. FEMS Microbiology Letters. 2000;182(2):333–7. pmid:10620688
  56. 56. Gur C, Maalouf N, Shhadeh A, Berhani O, Singer BB, Bachrach G, et al. Fusobacterium nucleatum supresses anti-tumor immunity by activating CEACAM1. Oncoimmunology. 2019;8(6):e1581531–e. pmid:31069151
  57. 57. Brewer ML, Dymock D, Brady RL, Singer BB, Virji M, Hill DJ. Fusobacterium spp. target human CEACAM1 via the trimeric autotransporter adhesin CbpF. J Oral Microbiol. 2019;11(1):1565043-. pmid:30719234
  58. 58. Kaplan CW, Lux R, Haake SK, Shi W. The Fusobacterium nucleatum outer membrane protein RadD is an arginine-inhibitable adhesin required for inter-species adherence and the structured architecture of multispecies biofilm. Molecular microbiology. 2009;71(1):35–47. pmid:19007407
  59. 59. Weiland-Bräuer N, Kisch MJ, Pinnow N, Liese A, Schmitz RA. Highly Effective Inhibition of Biofilm Formation by the First Metagenome-Derived AI-2 Quenching Enzyme. Frontiers in microbiology. 2016;7:1098-. pmid:27468282
  60. 60. Yerushalmy O, Khalifa L, Gold N, Rakov C, Alkalay-Oren S, Adler K, et al. The Israeli Phage Bank (IPB). Antibiotics (Basel). 2020;9(5):269.