Isolation of Staphylococcus bacterial strains
Among the wound samples isolated from the patients at RIMS, Kadapa, many bacterial spp have been isolated and identified. From the list, we selected only the potential pathogenic i.e., Staphylococcus sps based on morphological, different culture methods (mannitol salt agar, blood agar) physiological and molecular characterization. Before collecting the samples, the wounds were cleaned with phosphate-buffered saline (PBS), and pus samples were obtained of any gender and age with sterile swabs and brought to the laboratory within 1 hour to avoid the wound swabs drying up. The isolates were identified as Staphylococcus sps by the application of culture and standard biochemical tests for identification of Staphylococcus sps as follows: Gram-staining, catalase, Indole, oxidase, mannitol, sucrose and lactose fermentation. The MDR clinical isolates were selected for further analysis. The Staphylococcus spp were preserved at -80 C in Luria Broth medium containing 30% glycerol until further processing.
Isolation of SVV09-A phages
Sewage samples for phage isolation were collected from the sources of public ponds in and around hospitals in Kadapa. The bacterial strain was grown overnight in Luria Bertani broth (LB) at 37°C with albumin. LB soft agar overlays were utilized for phage experiments (including isolation and plaque counting). A double-layer agar method was employed for phage isolation and propagation. Sewage water was centrifuged to remove debris, followed by inoculation with Staphylococcus ureilyticus, incubation, and filtration(20-24). A plaque assay was used for detection.
Phage purification
Phage propagation involves incubating isolated phages with particular host bacteria, followed by the centrifugation and filtration of the upper agar layer. To isolate phages, we choose one plaque and placed it into salt magnesium buffer (SM) combined with Manganese chloride (MnCl2) followed by centrifugation to settle debris, then filtered and the double layer agar technique was employed. Following 18-24 hours, the earlier procedures were repeated to isolate phages. The clear phages were kept in suitable SM buffer with minor adjustments (including the addition of 0.002% w/v gelatin and albumin) at 4◦C and with 30% glycerol at −80◦C
Host Range Analysis
The isolated phage host range was tested on several pathogenic bacterial strains. The bacterial strains that were tested ( Staphylococcus ureilyticus strains and Staphylococcus aureus strains) were clinical pathogens and wound samples collected at the Rajiv Gandhi Institute of Medical Sciences in Kadapa. The susceptibility of the phage was evaluated using the spot assay method. Plates were inverted and incubated overnight, then examined for plaque presence with negative control. In summary, 100 µl of overnight bacterial host cultures (108 –109 CFU/ml) were combined with 2.5 ml of 0.7 % soft agar at 45 °C. The mixture was subsequently introduced to a 1.5% solid agar plate with the addition of MnCl2 metal ion to increase the adsorption on double layer agar plates. Following solidification, 10 µl aliquots of phage suspension (1.0 × 108 PFU) were applied to the lawn of host bacteria. The plates were dried and incubated at 37°C for a duration of 18-24 hours(25-26). The area of clearance observed at the site of phage inoculation indicated that the host was susceptible to the corresponding phage.
Temperature, pH, and the impact of metal ions on the phage adsorption rate were among the other external variables that were examined.
Effect of temperature on the stability and viability of the phages
Thermal stability assessments were conducted to evaluate the effect of environmental factors on phage development, as stability is vital for the preservation of lytic phages(27). In this experiment, phage filtrates at a concentration of 1 × 10^9 PFU/ml were prepared in microcentrifuge tubes and exposed to varying temperatures: 30°C, 37°C (control), and 45°C. The samples were incubated for durations of 10, 30, 60, and 90 minutes. After incubation, the double layer agar technique was employed to assess the lytic efficacy of the phages under different temperature conditions, with the results being compared to those of the control at 37°C (24,28-29).
Stability of phages at different pH
The effects of varying acidic and alkaline pH levels on phages were systematically examined. Purified phages, maintained at known concentrations, were prepared in SM buffer across a range of pH values, specifically pH 2, pH 4, pH 6, pH 7, and pH 8. This incubation occurred for one hour at a controlled temperature of 37 °C, as referenced in (24,30). Subsequently, the phage lysate was serially diluted (up to a 10-fold dilution) with SM buffer, combining 1 mL of the diluted solution with 0.5 mL of host culture (S. ureilyticus) at an optical density (O.D.) of 0.6. This mixture was then incubated for 30 minutes under the initial conditions previously mentioned. Phages cultivated in the pH 7 solution were designated as the control group for this experiment.
Effect of metal ions on phage adsorption
Research was conducted to investigate the influence of calcium, magnesium, zinc, and manganese ions on the adsorption of phages. The primary objective was to assess how divalent metal ions affect the rate of phage adsorption, utilizing solutions of CaCl2, MgSO4, ZnCl2, and MnCl2 (31,32). An overnight culture of S. ureilyticus was prepared, achieving an O.D. of 0.6. This culture was then distributed into four autoclaved flasks, with 25 ml designated for each flask. One set of flasks was inoculated with 500 μL of phage containing 1x10^9 PFU, serving as the control group. The remaining set received 500 μL of phage in conjunction with 250 μL of 10 mM solutions of CaCl2, MgSO4, ZnCl2, and MnCl2. The flasks were incubated under continuous shaking at 120 rpm and 37°C. Samples were collected from both groups at specified intervals: 10, 20, 30, 40, 50, and 60 minutes.
Transmission Electron Microscopy
We conducted transmission electron microscopy (TEM) analysis to enhance our understanding of the structural determinants of the phage regarding its infectivity and host range, following the protocol (33). Transmission electron microscopy (TEM) was utilized for the identification and classification of bacteriophages at Jamia Hamdard University, New Delhi. Bacteriophage isolate was applied to the grids (carbon film copper grids) and negatively stained using 2% uranyl acetate. The stained samples were dried with filter paper and examined with a transmission electron microscope. Phage classification and identification were performed in accordance with the International Committee on Taxonomy of Viruses recommendations
Phage DNA extraction
Genomic DNA was isolated utilizing the Phenol-Chloroform extraction technique in conjunction with CTAB (35). To effectively degrade the bacterial genome, the sample treated with DNase and RNase were subjected to heating at 37 °C for a duration of one hour. Following this incubation, the enzymes were inactivated by further heating at 75 °C for ten minutes. A total of 0.5 mL of a 1% SDS solution was added to the remaining 500 µL of phage lysate. To denature the phage protein capsids, a 5 µL aliquot of proteinase K (at a concentration of 20 mg/mL) was added, and the mixture was incubated overnight at 56 °C. After the digestion, 1 mL of a phenol–chloroform–isoamyl alcohol solution was placed, and the mixture was centrifuged for five minutes at 6000 rpm. A quarter volume of sodium acetate was then incorporated, along with an equal volume of isopropanol and followed by centrifugation for fifteen minutes at 13,000 rpm, the supernatant was collected. These steps were followed by two sequential washing by the addition of 1 mL of 70% ethanol, and further centrifugate the content for two minutes at 13,000 rpm. Finally, the DNA pellet was air dried (36,37).
Phage DNA purification and Quantification
The phage genomic DNA was extracted and quantified using both the phenol-chloroform extraction method and the CTAB method for column purification. The evaluation was conducted with a Nanodrop™ Lite Spectrophotometer (Thermo Fisher Scientific Limited) to determine the DNA concentration, which was measured at 260 nm and yielded a total reading of 46.7 ng/µL. Furthermore, the purity, quality, and size of the DNA were analyzed through agarose gel electrophoresis.
Whole genome Illumina sequencing of SVV09-A phage:
Phage DNA was sequenced utilizing Illumina technology on the HiSeq 2500 platform (Illumina, USA) at Eurofins Pvt. Ltd. in Bangalore. This sequencing achieved a genome coverage of 30X, generating a total of 3.5 GB of data for the sample.
Library preparation and Data analysis
Paired end sequencing was prepared from the DNA sample using NEB Next® UltraTM II FS DNA Library Preparation kit for Illumina. 100ng of DNA sample was processed for enzymatic fragmentation using FS enzyme mix supplied in the kit to generate a mean fragment distribution of 200-300bp. The fragmentated DNA samples were then subjected to end-repair and adapter ligation as per the kit recommendation. The adapter ligated products were purified using AMPureXP beads and processed for PCR amplification with the index primers to facilitate the hybridization onto a flow cell. The purified PCR-enriched libraries were evaluated on the 4200 Tape Station System (Agilent Technologies) utilizing high sensitivity D1000 screen tape according to the manufacturer's guidelines. The PE (Pair End) Illumina libraries were introduced to the NovaSeq X Plus for cluster creation and were utilized to sequence the entire genome of the phages. The high-quality PE reads of the sample were assembled using metaviral SPAdes assembler (v3.15.5).
Insilco analysis of SVV09-A phage
Genome validation, Gene Prediction and Phylogeny
In order to evaluate the completeness of phage genome, the genome was validated using PhageScope. The assessment of completeness was conducted utilizing Average Amino acid Identity method (AAI). Gene prediction was performed using GeneMarkS (phage). The National Center for Biotechnology Information (NCBI) Basic Local Alignment Search Tool (BLAST) used to manually annotate the predicted proteins’ specific functions. Circular plot was generated using Proksee which helps in the identification of gene position, GC content, GC Skew. The Mega 11 tool was additionally utilized to generate a diagram illustrating the comparative analysis of phage genomes alongside their related homologs.
Open reading frames (ORFs) Identification
The ORF Finder tool from NCBI was utilized to identify possible open reading frames (ORFs) in the phage DNA sequence. The software offers the range of each ORF along with the translation of its associated protein(38, 39).
Endolysin gene structural identification and Phylogenetic analysis
To analyse the structural features and evolutionary relationships of the endolysin gene, we performed a BLASTP search for the endolysin protein sequence to identify homologous sequences across different species. NCBI tools, were used for pairwise alignments of similar sequences from various organisms (40). We retrieved FASTA format sequences closely related to the target endolysin protein and filtered them based on alignment and similarity scores to ensure their relevance to the endolysin protein family(38, 41). Additionally, the identification of transmembrane helices within the protein was conducted using DeepTMHMM version 1.0, ensuring a comprehensive understanding of the protein’s structural and functional properties(38).
The endolysin protein sequence was meticulously analyzed to determine its secondary structure using the SOPMA algorithm, a method described (42) and further utilized (43). To assess the protein's physical and chemical properties, and other parameters were calculated through the Expasy-Protparam tool (44,45). These parameters included the theoretical isoelectric point, which indicates the pH at which the protein has no net charge, and the molecular weight, reflecting the size of the protein in Daltons. Additionally, we determined the total counts of positive and negative amino acid residues, the extinction coefficient, which is important for spectroscopic analyses, the protein's half-life in vivo, the instability index suggesting the stability of the protein in a cellular environment, the aliphatic index reflecting the proportion of hydrophobic residues, and the grand average hydropathy (GRAVY) score, which provides insights into the hydrophilic or hydrophobic nature of the protein. To predict the protein's folding state, we employed the FoldIndex program (46), which utilizes the amino acid composition to estimate the likelihood of a given conformation. A three-dimensional representation of the endolysin protein was constructed in the Swiss-Model environment, a platform known for its ability to generate homology models based on template sequences. Following the model generation, we conducted a thorough analysis to identify the model with the highest sequence identity to the chosen template. This selection was based on comparative metrics, including the GMQE (Global Model Quality Estimate) score and the QMEAN score, both of which provide quantitative measures of model reliability. Elevated GMQE and QMEAN (Qualitative Model Energy Analysis) scores are indicative of a high-quality model, suggesting that it closely resembles the true structure of the protein (47). Finally, to assess the conformational viability of the predicted model, we generated a Ramachandran plot using the RAMPAGE server (48). This plot provides visual insights into the distribution of dihedral angles in the protein structure, allowing us to validate the sterically allowed conformations of amino acid residues. Through these comprehensive analyses, we gained valuable understanding of the structural and functional aspects of the endolysin protein.