Next Article in Journal
Sex-Dependent Differences in Predictive Value of the C2HEST Score in Subjects with COVID-19—A Secondary Analysis of the COLOS Study
Next Article in Special Issue
Identification and Characterization of vB_PreP_EPr2, a Lytic Bacteriophage of Pan-Drug Resistant Providencia rettgeri
Previous Article in Journal
SARS-CoV-2 and Influenza Virus Co-Infection Cases Identified through ILI/SARI Sentinel Surveillance: A Pan-India Report
Previous Article in Special Issue
Immune Response of an Oral Enterococcus faecalis Phage Cocktail in a Mouse Model of Ethanol-Induced Liver Disease
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 14
Issue 3
10.3390/v14030626
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Insights into Gene Transcriptional Regulation of Kayvirus Bacteriophages Obtained from Therapeutic Mixtures

by
Sara Arroyo-Moreno
1,
Colin Buttimer
2,
Francesca Bottacini
1,
Nina Chanishvili
3,
Paul Ross
2,
Colin Hill
2,4 and
Aidan Coffey
1,2,*
1
Department of Biological Sciences, Munster Technological University, T12 P928 Cork, Ireland
2
APC Microbiome Ireland, University College, T12 YT20 Cork, Ireland
3
George Eliava Institute of Bacteriophage, Microbiology & Virology, Tbilisi 0160, Georgia
4
School of Microbiology, University College Cork, T12 YN60 Cork, Ireland
*
Author to whom correspondence should be addressed.
Viruses 2022, 14(3), 626; https://doi.org/10.3390/v14030626
Submission received: 25 January 2022 / Revised: 13 March 2022 / Accepted: 15 March 2022 / Published: 17 March 2022
(This article belongs to the Special Issue Bacteriophages from Biology to Therapy for the 21st Century)
Browse Figures
Versions Notes

Abstract

:
Bacteriophages (phages) of the genus Kayvirus of Staphylococcus aureus are promising agents for therapeutic applications. In this study, we isolated Kayvirus phages, SAM1 and SAM2, from the Fersisi commercial phage cocktail (George Eliava Institute, Tbilisi, Georgia), which exhibits high sequence homology with phage K (≥94%, BLASTn). We found that phages SAM1 and SAM2 infected 95% and 86% of 21 MRSA of differing sequence types (MLST, SCCmec type) obtained from the Irish National MRSA collection, respectively. We conducted differential transcriptomic analysis by RNA-Seq on phage SAM1 during host infection, showing differential expression of its genes at different points during host infection. This analysis also allowed the identification of potentially adverse outcomes in the application of these phages to target MRSA as therapy. The interaction of phage SAM1 on the host caused the upregulation of prophage genes. Additionally, phage infection was found to cause the slight upregulation of host genes implicated in virulence factors relating to hemolysins, immune evasion, and adhesion, but also the downregulation of genes associated with enterotoxins. The findings of this study give further insights into the biology of kayviruses and their use as therapeutics.

1. Introduction

Phages are the most abundant biological entities in the biosphere, playing a critical role in bacterial biology, diversity, and evolution. Their global population has been estimated to be on the order of 1031 [1]. Lytic phages have been employed to treat bacterial infections since their discovery in the early 20th century. However, the commercialization of antibiotics from the 1940s caused phage therapy, the use of phages for therapeutic purposes, to decline in Western countries. Recent decades have seen increased research into the development of phage therapy as a viable alternative to antibiotics to eliminate multidrug-resistant pathogens.
The genus Staphylococcus includes several species that inhabit humans and animals, with a number of these species being opportunistic pathogens of humans. Staphylococcus aureus can cause infections, which manifest through various symptoms, ranging from relatively mild to life-threatening. The treatment of these conditions is becoming increasingly problematic due to the widespread dissemination of antibiotic-resistant strains, such as methicillin-resistant Staphylococcus aureus (MRSA) [2]. Naturally occurring virulent phages that infect and kill a wide range of S. aureus strains may become an alternative treatment for otherwise incurable infections caused by antibiotic-resistant staphylococci. However, when selecting phages for therapeutic applications, care must be taken to avoid those with virulence factors or genes implicated in lysogeny [3,4,5].
The Kayvirus genus (Twort-like) of phages of S. aureus are a promising agent for phage therapy due to their broad host range and high killing activity against clinical isolates of this bacterium [6,7]. These phages are often found in commercial therapeutic phage preparations [8,9,10]. They have been demonstrated to be helpful in the treatment of human infections, such as diabetic toe ulcers and musculoskeletal infections, where antibiotic therapy has failed or had limited success [11,12]. Furthermore, they have the ability to mutate and evolve to overcome bacterial defences, enabling the updating of phage therapeutics [6]. This is an advantage over traditional antibiotics, as resistance can result in their reduced antibiotic efficacy, where they can become ineffective as therapeutics. The genus Kayvirus belongs to the subfamily Twortvirinae within the family Herelleviridae [13]. Phages in this group are closely related and are strictly lytic [14]. Morphologically, they are characterized by an isometric head and a long contractile tail (170–220 nm). Their average genome size is about 130 kbp, featuring a pair of long direct repeats at their ends, up to several-thousand bases in length [15]. Most phages from this genus utilize the backbone of teichoic acid, a ubiquitous molecule associated with the cell wall of Staphylococcus, as their host cell receptor [16].
This study examines the Kayvirus Staphylococcus phages SAM1 and SAM2 obtained through Fersisi commercial phage preparation (George Eliava Institute, Tbilisi, Georgia). The George Eliava Institute developed the product approximately 15–20 years ago, based on Pyophage, albeit with fewer phage components. We examine the host range of two phages on a panel of MRSA representative sequence types obtained from an Irish MRSA collection from St James’s Hospital, Dublin. We also performed RNA-Seq analysis of SAM1 during the infection of its host S. aureus strain E1185(IV)ST12 to improve current understanding of how Kayvirus phages transcribe their genes during host infection but also investigate the response of the host to such infection. This analysis gives insight into potential safety concerns of using these phages on MRSA in a therapeutic context. This includes the implications of virulent phage infection on host-associated pathogenicity islands (SaPIs) or prophages, as both elements are important in the horizontal transfer of antibiotic resistance and virulence genes [17,18,19].

2. Results

2.1. Host Range of Staphylococcus Phages SAM1 and SAM2

In a previous study, we assessed the host range of Staphylococcus phage K on a panel of 21 Irish MRSA isolates obtained from the Irish MRSA collection at St. James’s Hospital Dublin, where phage K was found to form plaques on 43% (9/21) of the tested isolates [14]. Among the resistant strains, two were selected to isolate novel phages from Fersisi therapeutic phage preparation from the George Eliava Institute, namely strains E1185(IV)ST12 and E1139(IV)ST45. In the case of these two isolates, phage K showed limited infectivity on E1185(IV)ST12 (efficiency of plaquing, 1.16 × 10−6) and no ability to plaque on strain E1139(IV)ST45.
Phages SAM1 and SAM2 were isolated using the resistant strain E1139(IV)ST45 and E1185(IV)ST12, respectively. Host range studies for both phages were performed on 21 MRSA strains and the usual propagating host for phage K, strain DPC5246 (Table 1). All the MRSA isolates were typed based on MLST (based on seven loci) and the composition of their SCCmec element [20,21]. SAM1 was found to be capable of lysing 21 of the 22 S. aureus isolates. The efficiency of the plaquing (EOP) values for SAM1 were found not to be less than 1 × 10−3, with a maximum value of 1.2. For phage SAM2, it could form plaques on 19 of the 22 staphylococcal isolates, with EOP values ranging from 7.37 × 10−1 to 75.47. This phage had a low EOP on strain E1185(IV)ST12, the host for phage SAM1. Hence, there was complete coverage of the 22-strain collection with these two phages.

2.2. Genomic Features of Phages SAM1 and SAM2

The genome sequences obtained for phages SAM1 and SAM2 were 139,310 bp and 137,617 bp in length (excluding terminal repeats), with 212 and 211 CDSs, respectively. Four tRNA genes were found to be associated with each of their genomes, which possessed a GC content of 30%. These genomic features were similar to those reported for other phages of the Kayvirus group [10,14]. The CDSs were identified to encode proteins associated with transcription, DNA replication, virion structure, host lysis, and putative genes for homing endonucleases. No genes were predicted to be implicated in bacterial virulence, antibiotic resistance, or lysogeny (Supplementary Information S1, Tables S1 and S2).
A nucleotide pairwise sequence alignment based on BLASTn showed the genomes of SAM1 and SAM2 to possess 97% nucleotide sequence identity, while they shared 94% and 95% identity with phage K, respectively. A comparison of the proteomes of the phages SAM1, SAM2, and K showed that 19 ORFs were not present in all three phages. Twelve ORFs were found to be associated with phage SAM1 and SAM2, which was not found with phage K. The predicted size of their gene products ranged from 3.7 to 46.9 kDa. Apart from the ORF SAM1_142/SAM2_141, annotated as a uracil-DNA glycosylase, no function could be assigned to their gene products (Supplementary Information, Table S3).

2.3. Phage SAM1 Gene Transcription Analysis

The RNA sequencing of a single infection cycle of phage SAM1 on the S. aureus strain E1185(IV)ST12 at a multiplicity of infection (MOI) of 10 was performed. Under the conditions set out for this experiment, phage SAM1 possessed a latent period of 50 min. Therefore, we collected samples at 15, 35, and 45 min for analysis, with biological triplicates taken for each time-point. Unique reads that mapped to the phage or host over specified time-points ranged between 42–51% and 49–59% of the total reads, respectively (Supplementary Information S2, Figure S1). A principal component analysis (PCA) of the reads that mapped to SAM1 showed that the samples collected from the same time-point were more similar than those obtained from different times (Supplementary Information S2, Figure S2 and Table S4).
Gene expression was found to have occurred throughout the genome of phage SAM1 among the time-points selected for the RNA-Seq analysis. No region on the genome of this phage was found to have an absence of aligned reads (Figure 1). A previous transcriptional study of kayvirus phage vB_SauM-515A1 produced similar findings [15]. However, locations on the genome of SAM1 with low read depth were identified between genes SAM_75 and SAM_76, and between SAM1_213 and SAM1_214. The former region’s low GC content (20%) indicates that it may be the proximal locus for the origin of replication for the SAM1 genome [22]. However, a differential expression analysis conducted using DESeq2 showed significantly differentially expressed genes (p < 0.05) of the phage throughout infection (Figure 2). Differences in gene expression were typically observed in genes belonging to the same transcriptional unit, encoding proteins of related function and, possibly, co-transcribed. Notably, a sharp decrease in sequence coverage often matched the predicted boundaries of these transcriptional units, which were typically situated between the boundaries of the ORFs (Figure 3). These low-coverage regions likely correlated with the common start and endpoints of the mRNA transcripts.
A comparison of the transcription at 15 min with other time-points (35 and 45 min) of phage infection showed a preferential expression for transcriptional units with genes implicated in nucleotide metabolism and DNA replication, repair, or recombination (as observed for transcriptional units spanning from SAM1_58–71 to SAM_124–154 and SAM_158–162). This observation was also made for three of the four tRNA genes of the phage, which would make these tRNAs available throughout infection. A similar comparison showed that transcriptional units with genes implicated in virion structure (SAM_076–123 and SAM_155–157) and cell lysis, such as endolysin (SAM_072) and holin (SAM_073), were most upregulated at 35 min compared to 15 or 45 min.
The most notable differential gene expression was found in gene SAM1_045, which was highly upregulated at 45 min compared to its expression at 15 min. Its gene product is a putative cell-wall binding domain (LysM, HHpred: probability > 95%). This differential expression suggests that the gene product likely plays a role in host cell lysis. Examples of multimeric endolysins include PlyC of Streptococcus phage C1, where the cell-wall binding domain and peptidoglycan degradation are encoded on separate genes [23]. However, SAM_072 encodes an endolysin, where a cell-wall binding domain (SH3, IPR003646) is associated with a peptidoglycan domain (CHAP, IPR007921). Thus, the putative role of the gene product SAM1_045 at the late stage of infection is unclear.
A complete list of SAM1 genes that showed differential expression (p < 0.05, DESeq2) during host infection are presented in Supplementary Information S1, Table S5.
Kayvirus phages do not possess genes for RNA polymerase and therefore exploit those of the host bacterium. An examination of genes upregulated at 15 and 35 min showed the presence of putative promoters that were highly similar, resembling the canonical σ70 promoter with a consensus sequence -35 motif (TTGACA) and -10 motif (TATAAT). The genes upregulated at 45 min possessed similar promoters, but with a -10 motif that had a more divergent consensus sequence (Supplementary information S1, Table S5 and Figures S3–S5). These predicted promoters of SAM1 were similar to those previously reported among other kayviruses [15]. Additionally, one of the transcriptional units upregulated at 15 min terminated with a gene encoding an RNA polymerase sigma factor (SAM1_154, HHpred: probability > 99%), which was upregulated at this time-point but downregulated at subsequent time-points. It is expected that this factor may play a role in the temporal gene expression of this phage, as such elements cause the alteration of promoter sequences recognized by RNA polymerase. In the case of SAM1, possibly causing the host RNA polymerase to recognize promotors with the alternative -10 motif identified among the genes upregulated at the later stages of phage infection.

2.4. Host Gene Transcription Analysis during Phage SAM1 Infection

After sequencing with a combination of long and short reads, the genome of S. aureus strain E1185(IV)ST12 returned as a complete circular chromosome of 2,782,853 bp containing 2704 coding sequences (CDSs) and 59 tRNA genes. The isolate also possesses a plasmid 26,732 bp in size with 27 CDSs. Its genome was observed to contain the SCCmec_type_IVc(2B) element with the strain associated with MLST type ST12. RNAseq analysis showed alternation of the transcriptional regulation of genes of the host during phage SAM1 infection over selected time-points. A complete list of genes that showed differential expression (p < 0.05, DESeq2) during host infection are presented in Supplementary Information S1, Table S6.
A comparison of host gene expression at 15 and 45 min during phage infection (Table 2) shows slight upregulation of the virulence factor genes implicated in the formation of hemolysins (hyl, hlgC). This was also detected for gene-encoding proteins with a role in subverting the human immune system, such as staphylococcal protein A (sbi) and complement inhibitor (scn). Similar observations were made for the products implicated in adherence of the bacterium to surfaces, such as extracellular adherence protein (eap) and fibronectin-binding protein A (fnbA). The virulence-factor genes in which slight downregulation was observed were those for staphylococcal enterotoxins (sel26, selZ). Transcriptional upregulation was also detected in a prophage (spanning ORFs LUU_550–797) located on the genome of the S. aureus strain E1185(IV)ST12 belonging to the genus Phietavirus. Similar studies with Pseudomonas aeruginosa have also shown the transcriptional upregulation of prophage elements during virulent phage infection [24].

3. Discussion

Numerous studies report the inefficacy of phage K against several MRSA strains [9,14]. In the present study, two phages belonging to the same genus as phage K were isolated from the George Eliava Institute’s Fersisi commercial phage mix, namely SAM1 and SAM2. These two phages present a much broader host range against MRSA strains from the National Reference Collection in comparison with phage K, which effectively only lysed nine strains out of twenty-one [14]. This strain collection was used because the differing sensitivity/resistance profiles that can be observed with Kayvirus phages can inform us as to which phages have the best potential for use in therapeutic mixes that might be adapted to eliminate specific sequence types from specific MRSA infections in Ireland. The strains selected for the isolation and propagation of phages SAM1 and SAM2, namely S. aureus E1185(IV)ST12 and S. aureus E1139(IV)ST45, were previously shown to have poor sensitivity to phage K and, indeed, two other phages from the genus Kayvirus [14]. As phage K can only infect S. aureus E1185(IV)ST12 at an extremely low EOP and cannot produce plaques on S. aureus strain E1139(IV)ST45, it is interesting to observe how minor differences between the SAM1/SAM2 genomes and that of phage K might play a role in phage virulence. The SAM1 and SAM2 phage genomes shared 94% and 95% nucleotide sequence identity (BLASTn), respectively, with phage K. More significant steps now need to be taken to understand the factors these differences contribute to the Kayvirus host range. This knowledge would enable enhanced phage therapeutics based on these phages. However, this study indicates the potential of key individual component phages from the Fersisi commercial phage mix for the elimination of Irish MRSA strains.
A transcriptional analysis of Kayvirus phage SAM1 showed the temporal regulation of gene expression as observed for other myoviruses, such as phage PAK_P3 of the genus Nankokuvirus, infecting Pseudomonas [25]. The genes of this phage can be primarily categorized as those upregulated during the early, middle, or late stage of infection. Early genes encode hypothetical proteins of unknown function, middle genes focus mainly on nucleic acid metabolism, and late genes are for virion structure. Under the conditions set out by the one-step growth curve of Staphylococcus phage SAM2 on host S. aureus strain E1185(IV)ST12 in this study, a latent period of 50 min was determined with three-time-points being selected for RNA-Seq analysis. In the case of Staphylococcus phage SAM1, the genes upregulated in early infection (15mins) were primarily those involved in nucleotide metabolism and functions related to DNA replication and repair. The genes upregulated toward the middle (35 min)/late (45 min) stage of infection were those implicated in virion structure and host cell lysis. These changes in gene expression levels possibly related to the gene SAM1_154 encoding an RNA polymerase sigma factor, upregulated at the 15 min time-point and homologous with that of Bacillus phage SPO1 gene gp34, implicated in the late gene expression of that phage [26].
During the infection of S. aureus E1185(IV)ST12 by phage SAM1, genes of a prophage element were shown to be upregulated at 45 min, indicating that this prophage was exiting the host chromosome and entering the lytic cycle. This observation draws attention to the concern that such elements are known to be agents of horizontal gene transfer [17,18]. Additionally, they can harbor genes that influence host fitness (i.e., pathogenicity, prophage immunity, biofilm formation, and stress response), potentially promoting their spread among sensitive strains of S. aureus [27]. Further investigation is merited to gauge the impact of Kayvirus lytic infection on indigenous prophage induction in S. aureus. It is worth mentioning that prophage induction can also occur during the application of antibiotics [28]. The analysis in this study also showed the slight upregulation of host genes implicated in virulence factors relating to hemolysins, immune evasion, and adhesion during phage infection. However, it must be stated that the application of these phages in the treatment of severe staphylococcus infections has been shown to be safe: no adverse reactions were reported among thirteen tested patients [29]. Thus, the observations in this study should not be a barrier to the use of kayviruses as therapeutics.

4. Materials and Methods

4.1. Bacterial Strains, Phage, and Growth Requirements

Phages SAM1 and SAM2 were isolated from the Fersisi commercial phage cocktail from the George Eliava Institute of Bacteriophage, Microbiology and Virology, Tbilisi, Georgia. The MRSA strains utilized for host range studies were predominantly the Irish National MRSA Reference collection, except for DPC5246, the routine propagation of phage K in our laboratory [14,30,31]. The MRSA strains used for the isolation and routine propagation of phages SAM1 and SAM2 were S. aureus E1185(IV)ST12 and S. aureus E1139(IV)ST45, respectively.
S. aureus isolates were routinely cultured in brain–heart infusion broth (BHI; Sigma-Aldrich, St. Louis, MO, USA) at 37 °C or on BHI plates containing 1.5% (w/v) bacteriological agar (Sigma-Aldrich). All strains were stocked in BHI containing 40% (v/v) glycerol and stored at −80 °C. Phage propagation was performed using plaque assay plating technique, propagation with BHI broth, and plate propagation with SM buffer [32].

4.2. Phage Host Range

Host range assay was performed for phages SAM1 and SAM2 using the plaque assay plating technique. This was performed in triplicate for three independent experiments. The efficiency of plaquing (EOP) was determined by dividing the phage titre on each test strain by the phage titre of the reference strain (S. aureus E1185(IV)ST12, in the case of phage SAM1, and S. aureus E1139(IV)ST45 for SAM2) [33].

4.3. DNA Extraction, Sequencing, and Genome Assembly

Phage DNA extraction was performed on high-titre phage lysates, as described by Ajuebor et al. in 2018 [14]. Briefly, phage samples were initially treated with MgCl2, followed by pre-treatment with DNase and RNase for 60 min at 37 °C. Next, subsequent treatments with SDS, EDTA, and proteinase K with further incubation for 60min at 55 °C were performed. DNA extractions were then performed with phenol/chloroform/isoamyl alcohol (25:24:1 v/v/v) and chloroform/isoamyl alcohol (24:1 v/v). DNA precipitation was achieved using sodium acetate and 95% (v/v) ethanol. DNA quality and quantity were estimated using a Nanodrop (ND-1000) and visualized in a 1% agarose gel electrophoresis.
Genomic sequencing was outsourced to genomic DNA MicrobesNG, Birmingham University, England, United Kingdom. For Illumina DNA sequencing, genomic DNA was isolated using the QIAamp DNA Mini Kit (Qiagen, Hilden, Germany). DNA was quantified in triplicates with a Quantit dsDNA HS assay in an Eppendorf AF2200 plate-reader. Genomic DNA libraries were prepared using Nextera XT Library Prep Kit (Illumina, San Diego, CA, USA) following the manufacturer’s protocol with the following modifications: two nanograms of DNA instead of one were used as input, and PCR elongation times was increased to 1 min from 30 s. DNA quantification and library preparation were carried out on a Hamilton Microlab STAR automated liquid handling system. Pooled libraries were quantified using the Kapa Biosystems Library Quantification Kit (Roche, Basel, Switzerland) for Illumina on a Roche Lightcycler 96 qPCR machine. Libraries were sequenced on an Illumina instrument using 250 bp paired-end protocol.
For Nanopore sequencing, broth cultures of each isolate were then pelleted and resuspended in the tube with cryopreservative (Microbank, Pro-Lab Diagnostics UK, Wirral, UK). Approximately 2 × 109 cells were used for high-molecular-weight (HMW) DNA extraction using Nanobind CCB Big DNA Kit (Circulomics, Baltimore, MA, USA). DNA was quantified with the Qubit dsDNA HS assay in an Invitrogen Qubit 3.0 (Thermo Fisher Scientific, Waltham, MA, USA). Long-read genomic DNA libraries were prepared with Oxford Nanopore SQK-RBK004 kit with Native Barcoding EXP-NBD104/114 (Oxford Nanopore Technologies, Oxford, UK) using 500 ng of HMW DNA. Barcoded samples were pooled together into a single sequencing library and loaded in a FLO-MIN106 (R.9.4) flow cell in a GridION (Oxford Nanopore Technologies).
Illumina reads were adapter-trimmed using Trimmomatic (v0.30) with a sliding window quality cut-off of Q15 [34]. For phage genomes (SAM1 and SAM2), assembly was performed on samples using SPAdes (v3.7) [35]. For bacterial host (S. aureus strain E1185(IV)ST12), genome assembly was conducted with Unicycler (v0.4.8) using short and long reads obtained from Illumina and Nanopore sequencing, respectively [36].

4.4. Genome Annotation and Comparison

RAST was employed to predict ORFs of sequenced phages, while ARAGORN was used to predict tRNA genes [37,38]. Further annotation was derived from BLASTP using the non-redundant nr/nt database [39], along with HHpred [40] and InterProscan [41]. Lipoproteins were identified using LipoP [42]. TMHMM was used to identify transmembrane domains [43]. Identification of Kayvirus virions structural proteins was determined by comparison with those found for Staphlococus phage ISP [44]. The proteomes of phages K (GenBank accession number: NC_005880), SAM1 and SAM2 were compared with Proteinortho using default parameters [45]. Potential promoters were identified by submitting up-stream sequences (250 bp) of predicted genes to multiple em for motif elicitation (MEME) [46]. The genome of Staphylococcus aureus strain E1185(IV)ST12 was annotated with PGAP [47]. Bacterial genomes (only E1185(IV)ST12 was made available for this manuscript) were analyzed with SCCmecFinder to find SCCmec element type [20] and MLST using PubMLST [48]. Virulence genes of the host were identified using the Virulence Factor Database (January 2022) [49].

4.5. RNA-seq Analysis of Phage SAM2

A one-step growth curve was created. S. aureus strain E1185(IV)ST12 was grown to an OD600 of 0.6 (corresponding to approximately 4.8 × 108 CFU/mL), followed by centrifugation of 2 mL volume in a microfuge to pellet bacteria. The bacterial pellets were resuspended in 1 mL of phage lysate to yield an approximate multiplicity of infection MOI of 10, following incubation at 37 °C for 1 min and subsequent centrifugation to eliminate unbound phages. Bacterial pellets with bound phages were resuspended in a total volume of 10 mL of BHI broth incubated aerobically in a water bath at 37 °C with agitation at 60 rpm. At time-points of 15, 30 and 45 min, 1.8 mL samples were taken into 1/10 volume of stop solution (10% phenol in 100% ethanol) and immediately placed on ice. These tubes were then centrifuged for 1 min at 16,000× g at room temperature. Pellets were resuspended in 800 μL of Tri Reagent (Sigma-Aldrich) and transferred to fresh tubes containing acid-washed beads. These tubes are placed in a ribolyser (MagNA Lyser, Roche), with beads subsequently removed by centrifugation. Supernatants were subject to total RNA extraction using a Direct-zol RNA miniprep kit (Zymo Research, Irvine, CA, USA, as per manufacturer instructions). RNA samples were subject to rRNA depletion before library creation. Paired-end sequencing was performed using an IlluminaHiSeq sequencer (Illumina Inc., San Diego, CA, USA).
RNA-seq reads were mapped to genomes (genomic features model GTF file5) of Staphylococcus phage SAM1 and that of host S. aureus strain E1185(IV)ST12 through Bowtie2 (v2.3.4.1) [50]. The alignment SAM files were further processed using Samtools (v1.7) to obtain BAM files necessary to obtain matrices with read counts per gene (normalized by gene length) [51]. Differential gene expression (DGE) analysis was performed using the R environment and the DESeq2 package available as part of the Bioconductor release [52]. As a pre-processing step, rows with zero counts (unmapped genes) were discarded from the count matrices. Differential expression analysis was performed on the count matrices using the DESeq function in DESeq2.

4.6. Accession Numbers for Sequence Data

Genome sequences of Stapholococus phages SAM1 and SAM2 were deposited under GenBank accession numbers MT338525 and MT226657, respectively. The genome and plasmid of S. aureus strain E1185(IV)ST12 were deposited under Genbank accession numbers CP089586 and CP089587, respectively. Sequence reads from the transcriptional analysis of phage SAM1 infection of S. aureus strain E1185(IV)ST12 were submitted to GEO under accession number GSE192733.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v14030626/s1. Table S1. Annotation of the Staphylococcal phage vB_SauM_SAM1 genome. Table S2. Annotation of the Staphylococcal phage vB_SauM_SAM1 genome. Table S3. Genes of Staphylococcus phages SAM1, SAM2 and K fall into the same orthologous groups. Table S4. Filtered differential gene expression results of Staphylococcus phage SAM1 infection of Staphylococcus aureus strain E1185(IV)ST12. Table S5. Predicted promoter sequences of Staphylococcus phage SAM1. Table S6. Filtered differential gene expression results of Staphylococcus aureus strain E1185(IV)ST12 during infection by Staphylococcus phage SAM1. Figure S1. Percentage of unique reads mapping to the genomes of Staphylococcus phage SAM1 and S. aureus strain E1185(IV)ST12. Figure S2. Principal component analysis graph presenting the correlation between all the samples used in this study concerning reads mapped to the genome of Staphylococcus phage SAM1. Figure S3. Consensus sequence logos for a putative promotor of Staphylococcus phage SAM1 among genes up regulated at 15 min. Figure S4. Consensus sequence logos for a putative promotor of Staphylococcus phage SAM1 among genes up regulated at 35 min. Figure S5. Consensus sequence logos for a putative promotor of Staphylococcus phage SAM1 among genes up regulated at 45 min.

Author Contributions

Conceptualization, S.A.-M., C.B. and A.C.; methodology; S.A.-M., C.B. and F.B.; investigation, S.A.-M. and C.B.; resources, A.C.; writing—original draft preparation, S.A.-M. and C.B.; writing—review and editing, F.B., N.C., P.R., C.H. and A.C.; visualization, C.B. and F.B.; supervision, A.C.; project administration, C.B. and A.C.; funding acquisition, A.C. All authors have read and agreed to the published version of the manuscript.

Funding

His publication had the financial support of Science Foundation Ireland Project Ref 06/RFP/BIM004 and 12/R1/2335. During the work of this manuscript, Colin Buttimer was a holder of the Government of Ireland Postdoctoral Fellowship 2019 (project ID GOIPD/2019/1097).

Acknowledgments

We want to acknowledge Padraig McKeown for his analysis of the genomes of SAM1 and SAM2, as well as Lorraine Draper and Colette Crowley for their guidance and assistance with RNA sample delivery.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kwan, T.; Liu, J.; DuBow, M.; Gros, P.; Pelletier, J. The complete genomes and proteomes of 27 Staphylococcus aureus bacteriophages. Proc. Natl. Acad. Sci. USA 2005, 102, 5174–5179. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Lee, A.S.; De Lencastre, H.; Garau, J.; Kluytmans, J.; Malhotra-Kumar, S.; Peschel, A.; Harbarth, S. Methicillin-resistant Staphylococcus aureus. Nat. Rev. Dis. Prim. 2018, 4, 1–23. [Google Scholar] [CrossRef] [PubMed]
  3. Xia, G.; Wolz, C. Phages of Staphylococcus aureus and their impact on host evolution. Infect. Genet. Evol. 2014, 21, 593–601. [Google Scholar] [CrossRef] [PubMed]
  4. Bondy-Denomy, J.; Qian, J.; Westra, E.R.; Buckling, A.; Guttman, D.S.; Davidson, A.R.; Maxwell, K.L. Prophages mediate defense against phage infection through diverse mechanisms. ISME J. 2016, 10, 2854–2866. [Google Scholar] [CrossRef] [Green Version]
  5. Łobocka, M.; Hejnowicz, M.S.; Dąbrowski, K.; Gozdek, A.; Kosakowski, J.; Witkowska, M.; Ulatowska, M.I.; Weber-Dąbrowska, B.; Kwiatek, M.; Parasion, S.; et al. Genomics of staphylococcal Twort-like phages--potential therapeutics of the post-antibiotic era. Adv. Virus Res. 2012, 83, 143–216. [Google Scholar] [PubMed]
  6. Botka, T.; Pantůček, R.; Mašlaňová, I.; Benešík, M.; Petráš, P.; Růžičková, V.; Havlíčková, P.; Varga, M.; Žemličková, H.; Koláčková, I.; et al. Lytic and genomic properties of spontaneous host-range Kayvirus mutants prove their suitability for upgrading phage therapeutics against staphylococci. Sci. Rep. 2019, 9, 5475. [Google Scholar] [CrossRef]
  7. Melo, L.D.R.; Brandão, A.; Akturk, E.; Santos, S.B.; Azeredo, J. Characterization of a new Staphylococcus aureus Kayvirus harboring a lysin active against biofilms. Viruses 2018, 10, 182. [Google Scholar] [CrossRef] [Green Version]
  8. Abdelkader, K.; Gerstmans, H.; Saafan, A.; Dishisha, T.; Briers, Y. The preclinical and clinical progress of bacteriophages and their lytic enzymes: The parts are easier than the whole. Viruses 2019, 11, 96. [Google Scholar] [CrossRef] [Green Version]
  9. El Haddad, L.; Ben Abdallah, N.; Plante, P.-L.; Dumaresq, J.; Katsarava, R.; Labrie, S.; Corbeil, J.; St-Gelais, D.; Moineau, S. Improving the Safety of Staphylococcus aureus Polyvalent Phages by Their Production on a Staphylococcus xylosus Strain. PLoS ONE 2014, 9, e102600. [Google Scholar] [CrossRef] [Green Version]
  10. McCallin, S.; Sarker, S.A.; Sultana, S.; Oechslin, F.; Brüssow, H. Metagenome analysis of Russian and Georgian Pyophage cocktails and a placebo-controlled safety trial of single phage versus phage cocktail in healthy Staphylococcus aureus carriers. Environ. Microbiol. 2018, 20, 3278–3293. [Google Scholar] [CrossRef]
  11. Fish, R.; Kutter, E.; Wheat, G.; Blasdel, B.; Kutateladze, M.; Kuhl, S. Bacteriophage treatment of intransigent Diabetic toe ulcers: A case series. J. Wound Care 2016, 25, S27–S33. [Google Scholar] [CrossRef]
  12. Onsea, J.; Soentjens, P.; Djebara, S.; Merabishvili, M.; Depypere, M.; Spriet, I.; De Munter, P.; Debaveye, Y.; Nijs, S.; Vanderschot, P.; et al. Bacteriophage Application for Difficult-To-Treat Musculoskeletal Infections: Development of a Standardized Multidisciplinary Treatment Protocol. Viruses 2019, 11, 891. [Google Scholar] [CrossRef] [Green Version]
  13. Adriaenssens, E.M.; Brister, J.R. How to Name and Classify Your Phage: An Informal Guide. Viruses 2017, 9, 70. [Google Scholar] [CrossRef] [Green Version]
  14. Ajuebor, J.; Buttimer, C.; Arroyo-Moreno, S.; Chanishvili, N.; Gabriel, E.M.; O’Mahony, J.; McAuliffe, O.; Neve, H.; Franz, C.; Coffey, A.; et al. Comparison of Staphylococcus phage K with close phage relatives commonly employed in phage therapeutics. Antibiotics 2018, 7, 37. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Kornienko, M.; Fisunov, G.; Bespiatykh, D.; Kuptsov, N.; Gorodnichev, R.; Klimina, K.; Kulikov, E.; Ilina, E.; Letarov, A.; Shitikov, E. Transcriptional landscape of Staphylococcus aureus Kayvirus Bacteriophage vB_SauM-515A1. Viruses 2020, 12, 1320. [Google Scholar] [CrossRef]
  16. Xia, G.; Corrigan, R.M.; Winstel, V.; Goerke, C.; Gründling, A.; Peschel, A. Wall teichoic acid-dependent adsorption of staphylococcal siphovirus and myovirus. J. Bacteriol. 2011, 193, 4006–4009. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Canchaya, C.; Fournous, G.; Chibani-Chennoufi, S.; Dillmann, M.L.; Brüssow, H. Phage as agents of lateral gene transfer. Curr. Opin. Microbiol. 2003, 6, 417–424. [Google Scholar] [CrossRef]
  18. McCarthy, A.J.; Witney, A.A.; Lindsay, J.A. Staphylococcus aureus temperate bacteriophage: Carriage and horizontal gene transfer is lineage associated. Front. Cell. Infect. Microbiol. 2012, 2, 6. [Google Scholar] [CrossRef] [Green Version]
  19. Chen, J.; Ram, G.; Penadés, J.R.; Brown, S.; Novick, R.P. Pathogenicity island-directed transfer of unlinked chromosomal virulence genes. Mol. Cell 2015, 57, 138–149. [Google Scholar] [CrossRef] [Green Version]
  20. Kaya, H.; Hasman, H.; Larsen, J.; Stegger, M.; Johannesen, T.B.; Allesøe, R.L.; Lemvigh, C.K.; Aarestrup, F.M.; Lund, O.; Larsen, A.R. SCC mec Finder, a Web-Based Tool for Typing of Staphylococcal Cassette Chromosome mec in Staphylococcus aureus Using Whole-Genome Sequence Data. Msphere 2018, 3, e00612-17. [Google Scholar] [CrossRef] [Green Version]
  21. Enright, M.C.; Day, N.P.J.; Davies, C.E.; Peacock, S.J.; Spratt, B.G. Multilocus sequence typing for characterization of methicillin-resistant and methicillin-susceptible clones of Staphylococcus aureus. J. Clin. Microbiol. 2000, 38, 1008–1015. [Google Scholar] [CrossRef] [Green Version]
  22. Rajewska, M.; Wegrzyn, K.; Konieczny, I. AT-rich region and repeated sequences—The essential elements of replication origins of bacterial replicons. FEMS Microbiol. Rev. 2012, 36, 408–434. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Nelson, D.; Schuch, R.; Chahales, P.; Zhu, S.; Fischetti, V.A. PlyC: A multimeric bacteriophage lysin. Proc. Natl. Acad. Sci. USA 2006, 103, 10765–10770. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Blasdel, B.G.; Chevallereau, A.; Monot, M.; Lavigne, R.; Debarbieux, L. Comparative transcriptomics analyses reveal the conservation of an ancestral infectious strategy in two bacteriophage genera. ISME J. 2017, 11, 1988–1996. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Chevallereau, A.; Blasdel, B.G.; De Smet, J.; Monot, M.; Zimmermann, M.; Kogadeeva, M.; Sauer, U.; Jorth, P.; Whiteley, M.; Debarbieux, L.; et al. Next-Generation “-omics” Approaches Reveal a Massive Alteration of Host RNA Metabolism during Bacteriophage Infection of Pseudomonas aeruginosa. PLoS Genet. 2016, 12, e1006134. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Stewart, C.R.; Casjens, S.R.; Cresawn, S.G.; Houtz, J.M.; Smith, A.L.; Ford, M.E.; Peebles, C.L.; Hatfull, G.F.; Hendrix, R.W.; Huang, W.M.; et al. The genome of Bacillus subtilis bacteriophage SPO1. J. Mol. Biol. 2009, 388, 48. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Bondy-Denomy, J.; Davidson, A.R. When a Virus is not a Parasite: The Beneficial Effects of Prophages on Bacterial Fitness. J. Microbiol. 2014, 52, 235–242. [Google Scholar] [CrossRef]
  28. Goerke, C.; Koller, J.; Wolz, C. Ciprofloxacin and Trimethoprim Cause Phage Induction and Virulence Modulation in Staphylococcus aureus. Antimicrob. Agents Chemother. 2006, 50, 171. [Google Scholar] [CrossRef] [Green Version]
  29. Petrovic Fabijan, A.; Lin, R.C.Y.; Ho, J.; Maddocks, S.; Ben Zakour, N.L.; Iredell, J.R.; Khalid, A.; Venturini, C.; Chard, R.; Morales, S.; et al. Safety of bacteriophage therapy in severe Staphylococcus aureus infection. Nat. Microbiol. 2020, 5, 465–472. [Google Scholar] [CrossRef]
  30. O’Flaherty, S.; Ross, R.P.; Meaney, W.; Fitzgerald, G.F.; Elbreki, M.F.; Coffey, A. Potential of the polyvalent anti-Staphylococcus bacteriophage K for control of antibiotic-resistant staphylococci from hospitals. Appl. Environ. Microbiol. 2005, 71, 1836–1842. [Google Scholar] [CrossRef] [Green Version]
  31. O’Flaherty, S.; Coffey, A.; Edwards, R.; Meaney, W.; Fitzgerald, G.F.; Ross, R.P. Genome of staphylococcal phage K: A new lineage of Myoviridae infecting gram-positive bacteria with a low G+C content. J. Bacteriol. 2004, 186, 2862–2871. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Bonilla, N.; Rojas, M.I.; Cruz, G.N.F.; Hung, S.H.; Rohwer, F.; Barr, J.J. Phage on tap-a quick and efficient protocol for the preparation of bacteriophage laboratory stocks. PeerJ 2016, 4, e2261. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Gutiérrez, D.; Vandenheuvel, D.; Martínez, B.; Rodríguez, A.; Lavigne, R.; García, P. Two Phages, phiIPLA-RODI and phiIPLA-C1C, Lyse Mono- and Dual-Species Staphylococcal Biofilms. Appl. Environ. Microbiol. 2015, 81, 3336–3348. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Bolger, A.M.M.; Lohse, M.; Usadel, B. Genome analysis Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.A.; Dvorkin, M.; Kulikov, A.S.S.; Lesin, V.M.M.; Nikolenko, S.I.I.; Pham, S.; Prjibelski, A.D.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]
  36. Wick, R.R.; Judd, L.M.; Gorrie, C.L.; Holt, K.E. Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput. Biol. 2017, 13, e1005595. [Google Scholar] [CrossRef] [Green Version]
  37. Laslett, D.; Canback, B. ARAGORN, a program to detect tRNA genes and tmRNA genes in nucleotide sequences. Nucleic Acids Res. 2004, 32, 11–16. [Google Scholar] [CrossRef]
  38. Brettin, T.; Davis, J.J.; Disz, T.; Edwards, R.A.; Gerdes, S.; Olsen, G.J.; Olson, R.; Overbeek, R.; Parrello, B.; Pusch, G.D.; et al. RASTtk: A modular and extensible implementation of the RAST algorithm for building custom annotation pipelines and annotating batches of genomes. Sci. Rep. 2015, 5, 8365. [Google Scholar] [CrossRef] [Green Version]
  39. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  40. Söding, J.; Biegert, A.; Lupas, A.N.; Soding, J.; Biegert, A.; Lupas, A.N. The HHpred interactive server for protein homology detection and structure prediction. Nucleic Acids Res. 2005, 33, W244–W248. [Google Scholar] [CrossRef] [Green Version]
  41. Jones, P.; Binns, D.; Chang, H.-Y.; Fraser, M.; Li, W.; McAnulla, C.; McWilliam, H.; Maslen, J.; Mitchell, A.; Nuka, G.; et al. InterProScan 5: Genome-scale protein function classification. Bioinformatics 2014, 30, 1236–1240. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Juncker, A.S.; Willenbrock, H.; von Heijne, G.; Brunak, S.; Nielsen, H.; Krogh, A. Prediction of lipoprotein signal peptides in Gram-negative bacteria. Protein Sci. 2003, 12, 1652–1662. [Google Scholar] [CrossRef] [PubMed]
  43. Käll, L.; Krogh, A.; Sonnhammer, E.L.L. A combined transmembrane topology and signal peptide prediction method. J. Mol. Biol. 2004, 338, 1027–1036. [Google Scholar] [CrossRef] [PubMed]
  44. Vandersteegen, K.; Mattheus, W.; Ceyssens, P.J.; Bilocq, F.; de Vos, D.; Pirnay, J.P.; Noben, J.P.; Merabishvili, M.; Lipinska, U.; Hermans, K.; et al. Microbiological and Molecular Assessment of Bacteriophage ISP for the Control of Staphylococcus aureus. PLoS ONE 2011, 6, e24418. [Google Scholar] [CrossRef] [Green Version]
  45. Lechner, M.; Findeiß, S.; Steiner, L.; Marz, M.; Stadler, P.F.; Prohaska, S.J. Proteinortho: Detection of (Co-)orthologs in large-scale analysis. BMC Bioinform. 2011, 12, 124. [Google Scholar] [CrossRef] [Green Version]
  46. Bailey, T.L.; Boden, M.; Buske, F.A.; Frith, M.; Grant, C.E.; Clementi, L.; Ren, J.; Li, W.W.; Noble, W.S. MEME SUITE: Tools for motif discovery and searching. Nucleic Acids Res. 2009, 37, W202–W208. [Google Scholar] [CrossRef] [PubMed]
  47. Tatusova, T.; Dicuccio, M.; Badretdin, A.; Chetvernin, V.; Nawrocki, E.P.; Zaslavsky, L.; Lomsadze, A.; Pruitt, K.D.; Borodovsky, M.; Ostell, J. NCBI prokaryotic genome annotation pipeline. Nucleic Acids Res. 2016, 44, 6614. [Google Scholar] [CrossRef]
  48. Jolley, K.A.; Bray, J.E.; Maiden, M.C.J. Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications. Wellcome Open Res. 2018, 3, 124. [Google Scholar] [CrossRef]
  49. Liu, B.; Zheng, D.; Jin, Q.; Chen, L.; Yang, J. VFDB 2019: A comparative pathogenomic platform with an interactive web interface. Nucleic Acids Res. 2019, 47, D687–D692. [Google Scholar] [CrossRef]
  50. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [Green Version]
  51. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R. The Sequence Alignment/Map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 1–21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Viruses 14 00626 g001 550
Figure 1. Transcription of Staphylococcus phage SAM1 genome during infection of Staphylococcus aureus at 15 min, 35 min, and 45 min using standard RNA-Seq.
Figure 1. Transcription of Staphylococcus phage SAM1 genome during infection of Staphylococcus aureus at 15 min, 35 min, and 45 min using standard RNA-Seq.
Viruses 14 00626 g001
Viruses 14 00626 g002 550
Figure 2. Time course transcriptional response of the genome of Staphylococcus phage SAM1. Values are reported as log2 fold expression change of transcripts from three independent experimental runs taken at time-points 15 min, 35 min and 45 min (15 min was set as a reference level for differential gene expression analysis); only genes with a significantly different expression between any of the time-points are shown (p < 0.05 in DESeq2, n = 169 (non-redundant genes)).
Figure 2. Time course transcriptional response of the genome of Staphylococcus phage SAM1. Values are reported as log2 fold expression change of transcripts from three independent experimental runs taken at time-points 15 min, 35 min and 45 min (15 min was set as a reference level for differential gene expression analysis); only genes with a significantly different expression between any of the time-points are shown (p < 0.05 in DESeq2, n = 169 (non-redundant genes)).
Viruses 14 00626 g002
Viruses 14 00626 g003 550
Figure 3. Comparison of RNA-Seq sequence depth with gene map (region spanning bases 45,000 to 60,000) of Staphylococcus phage SAM1. The gene map comprises of arrows representing genes, their length and location, with the direction of the arrow indicating direction transcription. It can be observed that areas of low sequence coverage often correlate with the gene cluster boundary (red lines) that are up regulated at similar time points during phage infection of the host.
Figure 3. Comparison of RNA-Seq sequence depth with gene map (region spanning bases 45,000 to 60,000) of Staphylococcus phage SAM1. The gene map comprises of arrows representing genes, their length and location, with the direction of the arrow indicating direction transcription. It can be observed that areas of low sequence coverage often correlate with the gene cluster boundary (red lines) that are up regulated at similar time points during phage infection of the host.
Viruses 14 00626 g003
Table
Table 1. Host ranges of staphylococcal phages SAM1 and SAM2 on 22 Staphylococcus aureus strains, comprising of 21 MRSAs of differing sequence types obtained from the Irish National MRSA collection and DPC5246, the propagating host of phage K. The table documents the efficiency of plaquing (EOP) values.
Table 1. Host ranges of staphylococcal phages SAM1 and SAM2 on 22 Staphylococcus aureus strains, comprising of 21 MRSAs of differing sequence types obtained from the Irish National MRSA collection and DPC5246, the propagating host of phage K. The table documents the efficiency of plaquing (EOP) values.
S. aureus StrainSCCmec TypeMLSTBacteriophages
SAM1SAM2
EOP ± SDEOP ± SD
DPC5246Not detected711.2 ± 0.43.72 ± 0.48
0.0066(IIv)ST239II(2A) **365.1 × 10−1 ± 0.14.37 ± 1.8
0.1206(IV)ST250IVc(2B)125.74 × 10−1 ± 0.4no plaques
0.1239(III)ST239III(3A)2391.4 × 10−1 ± 0.056.6 ± 2.46
0.1345(II)ST5VIII(4A)82.23 × 10−1 ±0.11.96 ± 0.51
0073(III)ST239III(3A)2391.21 × 10−1 ± 0.1no plaques
0104(III)ST239III(3A)2399.65 × 10−2 ± 0.013.5 ± 1.19
0220(II)ST5IV(2B) **82.22 × 10−2 ± 0.00314.68 ± 4.31
0242(IV)ST30II(2A)4961.95 × 10−1 ± 0.179.35 ± 3.15
0308(IA)ST247I(1B)2473.41 × 10−1 ± 0.2611.85 ± 3.24
3045(IIv)ST8II(2A)81.94 × 10−1 ± 0.144.54 ± 2.13
3144(IIv)ST8II(2A)82,17 × 10−1 ± 0.2210.46 ± 5.18
3488(vv)ST8IV(2B) **86.72 × 10−2 ± 0.0111.11 ± 4.35
3581(IA)ST247VIII(4A) **81 × 10−1 ± 0.015.51 ± 1.25
3594(II)ST36III(3A)2395.77 × 10−1 ± 0.575.47 ± 2.32
3596(IIv)ST8VIII(4A) **82.29 × 10−1 ± 0.0528.89 ± 16.48
E1038(IIv)ST8II(2A)89.7 × 10−1 ± 0.9713.47 ± 3.47
E1139(IV)ST45IVa(2B)45no plaques1 *
E1174(IV)ST22IV(2B)221.4 × 10−1 ± 0.017.37 × 10−1
E1185(IV)ST12IVc121 *3.5 × 10−4 ± 3.3 × 10−4
E1202(II)ST496VIII(4A) **83.59 × 10−1 ± 0.220.83 ± 3.24
M03/0073(III)ST239III(3A)2391.32 ± 0.716.54 ± 0.52
0104(III)ST239III(3A)2391.2 ± 0.43.72 ± 0.48
0220(II)ST5IV(2B) **85.1 × 10−1 ± 0.14.37 ± 1.8
0242(IV)ST30II(2A)4965.74 × 10−1 ± 0.4no plaques
* Host strain of phage. ** SCCmec type: two possible types suggested by SCCmecFinder with a hit of greatest coverage selected.
Table
Table 2. List of virulence-factor genes of S. aureus strain E1185(IV)ST12 that experience up/downregulation during infection of Staphlococcus phage SAM1 (p < 0.05, DESeq2).
Table 2. List of virulence-factor genes of S. aureus strain E1185(IV)ST12 that experience up/downregulation during infection of Staphlococcus phage SAM1 (p < 0.05, DESeq2).
Virulence Factors (Locus_tag, Product, Gene)log2 Fold Change
15 min35 min45 min
LUU82_11680—bi-component gamma-hemolysin HlgCB subunit C—hlgC −0.11173 −0.0293 0.141026
LUU82_04645—alpha-hemolysin—hyl −0.15039 −0.04266 0.193054
LUU82_09700—staphylococcal protein A—spa −0.26288 −0.0223 0.285182
LUU82_09135—complement inhibitor SCIN—scn −0.18095 −0.16701 0.347955
LUU82_04665—complement convertase inhibitor—efb −0.44123 −0.06633 −0.507566
LUU82_00160—extracellular adherence protein Eap/Map—eap −0.18624 −0.03037 0.216609
LUU82_11260—fibronectin-binding protein FnbA—fnbA −0.43413 −0.09529 0.529417
LUU82_02360—staphylococcal enterotoxin type—sel26 0.415231 −0.08836 −0.32687
LUU82_09885—staphylococcal enterotoxin type Z—selZ 0.413472 −0.08253 −0.33094
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Arroyo-Moreno, S.; Buttimer, C.; Bottacini, F.; Chanishvili, N.; Ross, P.; Hill, C.; Coffey, A. Insights into Gene Transcriptional Regulation of Kayvirus Bacteriophages Obtained from Therapeutic Mixtures. Viruses 2022, 14, 626. https://doi.org/10.3390/v14030626

AMA Style

Arroyo-Moreno S, Buttimer C, Bottacini F, Chanishvili N, Ross P, Hill C, Coffey A. Insights into Gene Transcriptional Regulation of Kayvirus Bacteriophages Obtained from Therapeutic Mixtures. Viruses. 2022; 14(3):626. https://doi.org/10.3390/v14030626

Chicago/Turabian Style

Arroyo-Moreno, Sara, Colin Buttimer, Francesca Bottacini, Nina Chanishvili, Paul Ross, Colin Hill, and Aidan Coffey. 2022. "Insights into Gene Transcriptional Regulation of Kayvirus Bacteriophages Obtained from Therapeutic Mixtures" Viruses 14, no. 3: 626. https://doi.org/10.3390/v14030626

APA Style

Arroyo-Moreno, S., Buttimer, C., Bottacini, F., Chanishvili, N., Ross, P., Hill, C., & Coffey, A. (2022). Insights into Gene Transcriptional Regulation of Kayvirus Bacteriophages Obtained from Therapeutic Mixtures. Viruses, 14(3), 626. https://doi.org/10.3390/v14030626

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