A Novel Wide-Range Freshwater Cyanophage MinS1 Infecting the Harmful Cyanobacterium Microcystis aeruginosa

Microcystis aeruginosa, as one of the major players in algal bloom, produces microcystins, which are strongly hepatotoxic, endangering human health and damaging the ecological environment. Biological control of the overgrowth of Microcystis with cyanophage has been proposed to be a promising solution for algal bloom. In this study, a novel strain of Microcystis cyanophage, MinS1, was isolated. MinS1 contains an icosahedral head approximately 54 nm in diameter and a 260 nm-long non-contractile tail. The phage genome consists of a linear, double-stranded 49,966 bp DNA molecule, which shares very low homology with known phages in the NCBI database (only 1% of the genome showed weak homology with known phages when analyzed by megablast). The phage contains 75 ORFs, of which 23 ORFs were predicted to code for proteins of known function, 39 ORFs were predicted to code for proteins of unknown function, and 13 ORFs showed no similarity to any protein sequences. Transmission electron microscopy and phylogenetic analysis showed that MinS1 belongs to the family Siphoviridae. Various experiments confirmed that the phage could infect several different orders of cyanobacteria, including Chroococcales, Nostocales, Oscillatoriales, Hormogonales, and Synechococcales, indicating that it has a very broad host range. In addition, MinS1 has no known antibiotic tolerance genes, virulence genes, and tRNAs, and it is tolerant to temperature, pH, UV, and salinity, suggesting that MinS1 has good potential for application as a biological control agent against cyanobacterial blooms. This study expands the diversity and knowledge of cyanophages, and it provides useful information for the development of novel prevention and control measures against cyanobacterial blooms.


Introduction
Water eutrophication is a global water pollution problem, which causes massive production of cyanobacteria [1,2]. The latter challenge can also lead to cyanobacteria blooms [3], which have become a common occurrence in water bodies worldwide. Microcystis is one of the most pervasive bloom-forming cyanobacteria in freshwater ecosystems [4][5][6], consuming large amounts of dissolved oxygen and causing water quality degradation. Its metabolites, algal toxins, are hepatotoxic, neurotoxic, reproductive toxic, genotoxic, and tumor promoting, causing great economic losses to aquaculture and posing risks to the safety of aquatic products for consumption [7,8].
Viruses infecting cyanobacteria are referred to as cyanophages, and they can play major roles in the dynamics, genetic diversity, and structure of cyanobacterial communities [9,10]. These cyanophages inhibit the production and release of algal toxins and are considered to have significant potential as biological control agents for harmful cyanobacterial blooms. Compared with the existing physical and chemical cyanobacteria control methods, which are prone to take the good with the bad and secondary pollution, they have the advantages of low cost, not harming eukaryotes and other beneficial microorganisms, and not causing secondary pollution. However, most reported cyanophages have a long lysis cycle, a narrow cyanobacteria-killing spectrum, and strong host specificity [11,12]. Compared to marine cyanophages, which have been widely studied [13], there are very limited studies available pertaining to freshwater cyanophages [14]. Current research aims to build the foundations of cyanophage-based water bloom treatment, much of which are learnt from the role that these viruses play in bloom decay, particularly concentrating on an environmentally friendly treatment therapy to combat water blooms [15,16].
In the present work, we report on the morphological and biological characteristics, as well as the genomic information, of MinS1, which was isolated from freshwater in Fujian, China. Not only is MinS1 a novel M. aeruginosa cyanophage, but it also has a broad range of hosts. This work provides basic data for further understanding a wide range of cyanophage hosts and, additionally, enriches the database of freshwater cyanophage genome information, which is of great research significance. This work also provides basic data underpinning freshwater cyanophage genomes, and it supports the growing interest in using freshwater cyanophages to control bloom-forming cyanobacterium.

Cyanophage Isolation and Purification
The 30 mL samples of surface cyanobacterial bloom water used for cyanophage separation were collected from the Mayang Stream (24 • 32 37.81" N 117 • 46 26.40" E), in Fujian, China, and centrifuged at 12,000 g for 20 min at 4 • C. The temperature, date, and time of the water collection were 20 • C, 23 November 2020, at 10:17 a.m., respectively. The supernatant was filtered through nitrocellulose membranes with pore sizes of 0.45 µm and 0.22 µm. The filtrate was then added to a ten-times-larger volume of algal solution of the logarithmic growth M. aeruginosa strain FACHB-905, which can produce microcystin. The host concentration used in the experiment was 3.4 × 10 6 CFU/mL. Then, the solution was mixed and incubated in a light incubator (Ningbo Jiangnan instrument Factory, model number: GXZ-280B) at 25 • C, 2000 Lux, with a 12 h:12 h light-dark cycle for about three days until yellowing was evident. The procedure was repeated three times. The resulting lysate was then serially diluted ten-fold using BG-11 medium (Qingdao Haibo Biotechnology Co., Ltd. product number: HB8793) and used in double-layer agar plate spread experiments [17]. Once the growth of phagocytic spots was observed, several individual spots were dug up and suspended in 5 mL of FACHB-905 solution during the logarithmic phase. The aforementioned phagocytic spot experiment was repeated after yellowing until phagocytic spots consisting of uniform shapes and sizes had formed on each plate.

Transmission Electron Microscopy
A total of 30 mL of phage lysate was centrifuged at 6000× g for 10 min, and the supernatant was filtered through a 0.22 µm microporous filter. Thereafter, the supernatant was added to a 50 mL ultracentrifuge tube, and a 30% (w:v) sucrose solution was slowly injected from the bottom of the sample using a pipette and then centrifuged at 4 • C at 35,000× g for 1 h. The phage precipitate was resuspended into 200 µL of 1 × TNE (10 mmol/L Tris-HCl, 1 mmol/L EDTA, 100 mmol/L NaCl) in order to obtain the ultrapurified phage. Next, 1 µL of the ultra-purified phage (remaining samples could be stored directly at 4 • C) was dropped onto a grid surface, 2% uranyl acetate negative staining was undertaken, and then the excess staining solution was immediately removed. Phage morphology was then examined under a 100 kV transmission electron microscope (Hitachi H-7650, Tokyo, Japan) [18].

Host Range
Thirty cyanobacterial strains (containing five orders and eight families) were cultured in our laboratory, identified using 16S rRNA genes, and cultured to a concentration of approximately 3.4 × 10 6 PFU/mL in logarithmic stage to serve as test hosts. All information about test hosts can be found on Freshwater Algae Culture Collection at the Institute of Hydrobiology, and the culture conditions were the same as above. Then, 250 µL of algal solution and 50 µL of phage lysate (experimental group) were added into 96-well plate wells (Thermo Scientific, Waltham, MA, USA), incubated in a light incubator at 25 • C, 2000 Lux, with a 12 h:12 h light-dark cycle. An equal volume of 250 µL of algal solution and 50 µL of BG11 medium was also added into the other wells (control group); these were incubated under the same conditions as indicated above. Three parallel experiments were set up for each group. After 7 days, the lytic nature of the phage was determined by measuring the OD680, as well as the number, density, integrity, transparency, and edge clarity of the cells when observed using light microscopy.
In addition, 30 cyanobacterial strains were incubated until their log phase was achieved, and then they were spread in double-layer plates. Then, 5 µL of each phage solution was spotted onto each plate, which was dried and placed upside down in a light incubator in order to observe the appearance of any spots [19].

One-Step Growth Curve
One-step growth curve experiments were performed according to existing methods, with some minor modifications. A total of 5 mL of purified phage supernatant was added to 25 mL of the host strain and stirred for 30 min at 25 • C, 100 rpm on a shaker, followed by centrifugation at 6000 g for 10 min. The supernatant was discarded, and the precipitate was washed twice using PBS and then suspended in 30 mL of BG11 medium; this was then incubated in a lighted incubator. Samples were collected at 0, 6,12,18,24,30,36,42,48,54,60,66, and 72 h. All of the 1 mL samples were centrifuged (6000 g, 10 min) and filtered (using 0.45 µm and 0.22 µm filters), and titration was undertaken using the double-layer agar plate method. The burst size was calculated as the ratio of the number of free phages at the end of the logarithmic period to the number of infected host cyanobacteria cells at the beginning of the latent period. The experiment was repeated in triplicate.

Phage Stability under Environmental Stress
Cyanobacterial cleavage by cyanophages within natural environments is affected by water temperature, pH, salinity, ions, and UV. To examine the tolerance levels of the phage to these conditions, a published method was used with appropriate modifications [20,21]. The stability of the phage under different temperature treatments was analyzed by taking 1 mL of purified phage solution at different temperatures (0, 25, 40, 60, and 80 • C) at both 30 min and 60 min. Phage stability at different pH levels (2,3,4,5,6,7,8,9,10,11,12) was evaluated by adding the purified phage solution to different pH buffers for 1 h at 25 • C. To evaluate the effects of differing salt conditions on the phage, purified phage pellets were incubated at different NaCl concentrations (0%, 5%, 10%, 15% and 20%) at 25 • C for 1 h. Phage metal ion sensitivity was determined by incubating the phage in different metal ions (Ca 2+ , Mg 2+ , Fe 2+ ), at a final concentration of 10 mM, for 25 min at 25 • C. Phage sensitivity to UV was assessed by irradiating purified phage particles under a UV lamp (253.7 nm) for different amounts of time (0, 20, 40, 60, and 80 min).
Untreated phage solution was used as a control, and each of the treatments was repeated three times. The solutions were all titrated using the double-layer agar plate method. In order to corroborate the double-layer agar plate results, the treated phage from each group was added to the FACHB-905 host cyanobacteria to measure the OD680 values of the mixture after a 4-day incubation period. In these experiments, the original phage titer was 6.3 × 10 6 PFU/mL, and the host concentration was 3.9 × 10 7 PFU/mL. One-way ANOVA and SPSS 13.0 Duncan's new multiple range test were used to determine statistical significance. GraphPad Prism (8.0.2) was used for curve plotting and data statistical analysis.

DNA Extraction, Sequencing and Assembly
Phage DNA was extracted using a phenol-chloroform extraction method [22]. A 2 × 300 bp paired-end DNA library was constructed according to the manufacturer's instructions for the NEBNext Ultra™ II DNA Library PrepKit for Illumina, and the phage was genomically sequenced using the Illumina MiSeq (San Diego, CA, USA) sequencing platform. Low-quality (Q value < 20) reads and adaptors were filtered out using fastp. Clean reads were assembled using SPAdes 3.13.0 software (http://cab.spbu.ru/software/ spades/, 29 July 2021) [23]. Analysis of the phage ends was performed using an established method [24].

Isolation and Morphology of Cyanophage MinS1
M. aeruginosa cyanophage MinS1 (MinS1) [27] was isolated from the surface freshwater samples obtained from the Mayang Stream, located in Fujian, China, which had an outbreak of cyanobacterial blooms (more specific information is shown in Table S1). The plaques resulting from MinS1 lysis appeared clear and circular, with diameters of 3-4 mm, following a 4-day incubation (25 • C, 2000 Lux, 12 h:12 h light-dark cycle) on host algal plates at a multiplicity of infection (MOI) of 0.1 (Figure 1a). TEM of the purified phage particles showed that MinS1 had an isometric hexagon head measuring about 54 nm in diameter and a non-retractable long tail around 260 nm long ( Figure 1b). The head diameter of all Microcystis cyanophages ranged from 42 nm to 100 nm, with MinS1 having the second smallest head among these phagosomes. MinS1 is morphologically most similar to Mic1 [28] and has a non-contractible tail second only to Mic1 (400 nm) in length. Based on its morphology, and comparisons to the current International Committee on Taxonomy of Viruses (ICTV) classification system, MinS1 therefore belongs to the family Siphoviridae from the order Caudovirales. MinS1 is conserved at the China General Microbiological Culture Collection Center (CGMCC) under CGMCC No. 23089.

Life Cycle
To understand the growth kinetics of MinS1 in the host M. aeruginosa strain FACHB-905, one-step growth curves of MinS1 were performed using a modified soft-agar overlay method [29]. The results show that the MinS1 latent period lasted 36-42 h and was followed by a plateau period after 60 h ( Figure 2). Burst size is often considered for phage usage as biocontrol because it indicates the lytic ability of the phage. The MinS1 burst size was around 34 PFU per cell, and the total time duration of the one-step growth curve experiment lasted 72 h. The latent period and burst size of Microcystis cyanophages appeared to be highly variable, ranging from 6 h to 108 h and from 28 to 127 PFU per cell. It is important to note that different methods were used to count viral abundance. To date, only five strains of the Microcystis phage with one-step growth curves have been reported in the literature ( Table 2), three of which were hosted by the same host as MinS1, M. aeruginosa. The latent period for MaMV-DC was 24-48 h, with an outbreak of 80 PFU per cell [30]. Ma-LBP had a latent period of 11.2 h, followed by a burst size of 80 PFU per cell [31]. Another Microcystis Viruses 2022, 14, 433 5 of 17 cyanophage, Ma-LMM01, revealed a latent period of 6-12 h with a burst size of 50-120 PFU per cell [32]. In general, burst size is considered to be influenced by many factors, including bacterial/viral size, metabolic activity of the host, and phage and host characteristics [33]. A correlation between outbreak size and the environmental trophic state has also been proposed [34,35]; however, this concept still requires further validation. Similar to the physiological and biological characteristics of phages, the features that influence the burst size of each phage should not be neglected.

Life Cycle
To understand the growth kinetics of MinS1 in the host M. aeruginosa strain FACHB-905, one-step growth curves of MinS1 were performed using a modified soft-agar overlay method [29]. The results show that the MinS1 latent period lasted 36-42 h and was followed by a plateau period after 60 h ( Figure 2). Burst size is often considered for phage usage as biocontrol because it indicates the lytic ability of the phage. The MinS1 burst size was around 34 PFU per cell, and the total time duration of the one-step growth curve experiment lasted 72 h. The latent period and burst size of Microcystis cyanophages appeared to be highly variable, ranging from 6 h to 108 h and from 28 to 127 PFU per cell. It is important to note that different methods were used to count viral abundance. To date, only five strains of the Microcystis phage with one-step growth curves have been reported in the literature ( Table 1), three of which were hosted by the same host as MinS1, M. aeruginosa. The latent period for MaMV-DC was 24-48 h, with an outbreak of 80 PFU per cell [30]. Ma-LBP had a latent period of 11.2 h, followed by a burst size of 80 PFU per cell [31]. Another Microcystis cyanophage, Ma-LMM01, revealed a latent period of 6-12 h with a burst size of 50-120 PFU per cell [32]. In general, burst size is considered to be influenced by many factors, including bacterial/viral size, metabolic activity of the host, and phage and host characteristics [33]. A correlation between outbreak size and the environmental trophic state has also been proposed [34,35]; however, this concept still requires further validation. Similar to the physiological and biological characteristics of phages, the features that influence the burst size of each phage should not be neglected.

Host Range
Host infectivity tests showed that the phage MinS1 had polyvalent infectivity, as shown by the infection of 19 out of 30 cyanobacteria strains tested, containing five orders, Chroococcales, Nostocales, Oscillatoriales, Hormogonales, and Synechococcales (Table 1 and Supplementary Figure S1). This represents the broadest host range of the six reported Microcystis cyanophage strains (Ma-LMM01, MaMV-DC, Ma-LBP, vB_MelS-Me-ZS1, phiMa05 [36], and Mic1). This wide host range indicates that MinS1 has potentially significant applications, as water blooms are usually caused by multiple cyanobacteria [37]. Generally, Myoviruses have shown the broadest host range among the three families comprising the tailed viruses, whereas Podoviruses have the narrowest range [38]. Interestingly, a narrow host range was observed for Myoviruses Ma-LMM01, MaMV-DC, and phiMa05 and a Siphoviruses Mic1, whereas Siphoviruses vB_MelS-Me-ZS1 had a wide host range infecting 12 of the 15 host algal strains. Moreover, the host range of Ma-LBP has not previously been reported; it is a member of the Podoviridae family (Table 2). Host range may be influenced by the number of tested host algal strains. Moreover, the cyanophages' interaction in the environment, for example, their dynamics during cyanoHABs, may also affect their ability to infect algal strains. However, undeniably, this confirms the specificity and complexity of cyanophage-host interactions and the diversity within the Microcystis phages. Symbols: "+" suspective; "−" unsuspective; "*" means the strains can produce cyanotoxin according to references [39,40].
The interactions between cyanobacterial strains and environmental microorganisms and the co-evolution of phagosomes and cyanobacteria make it impossible for phagosomes to kill all of them, even if they are broad-spectrum phagosomes, and some of them will be selectively retained. In addition, the interaction between a phage with a variable host range and cyanobacteria can catalyze mutations and tolerant strains, promoting continuous turnover and, thus, reaching a dynamic equilibrium. Therefore, the application of broadspectrum phages does not destroy biodiversity. No reduction in species diversity has occurred during natural evolution due to the presence of phages. 3.4. Thermolability, pH, UV, Ions, Sensitivity, and Saline Stability As a biocontrol agent with potential applications, the stability of the cyanophage was determined using both double-layer agar and liquid infection OD680 measuring methods under several stress conditions, including differing salinity, pH, UV, temperature, and metal ion levels. MinS1 tolerated a wide range of pH values from 3-12. However, extreme acidic conditions, such as pH 2, resulted in a loss of cyanophage activity (Figure 3a). After both 30 min and 60 min incubation times, the cyanophage MinS1 showed stability from 0 • C to 40 • C. In particular, at 0 • C, the phage was more stable with better activity, and this result also suggests that phages are more stable when stored at 0 • C. The activity of this phage was, however, lost following exposure to temperatures of 60 • C and higher (Figure 3b). After 1 h of incubation, MinS1 withstood a wide variation of salinities ranging from 1-20%, and the highest titer of phage was observed when NaCl concentrations were at 15% (Figure 3c). MinS1 was isolated from freshwater, yet it exhibited a preference for high salinity. We reviewed the salinity tolerance of some phages isolated from freshwater and found that not only MinS1 but also other phages, such as vB_ValP_IME234, have a high salinity tolerance [41]. However, further studies are needed to confirm why phages prefer salt. As shown in Figure 3d, UV irradiation had an effect on the activity of the phage, gradually decreasing activity levels as exposure time increased. In addition, MinS1 activity was affected by Fe 2+ , but it was tolerant to some extent when exposed to Ca 2+ and Mg 2+ (Figure 3e). The yellowing status of the phage infected with host cyanobacteria under several groups of different treatments is shown in Supplementary Figure S2. Our study is the first report on the stability of multiple environmental stresses on cyanophages. In conclusion, MinS1 showed varying degrees of tolerance to temperature, pH, UV, ions, and salinity, and these tolerance characteristics indicate that MinS1 might provide more opportunities for survival in an aquatic environment, suggesting that it could be an alternative biological approach toward inhibiting cyanobacterial growth and reducing the accumulation of microcystins. Besides this, the development of multi-cyanophage cocktails may have antagonistic effects on water blooms formed by multiple cyanobacterial outbreaks [42].

Genome Features
The cyanophage genome was sequenced (a 2 × 300 paired-end sequence method) on an Illumina MiSeq sequencing platform. MinS1 was found to be a GC-rich double-stranded DNA cyanophage, 49,966 bp long, with very low homology to known sequences within the NCBI database (1% homology coverage and 74.06% identity to the Myoviridae sp. ct6sY2). MinS1 had no obvious termini and no antibiotic tolerance genes, virulence genes or tRNAs. linity, and these tolerance characteristics indicate that MinS1 might provide more opportunities for survival in an aquatic environment, suggesting that it could be an alternative biological approach toward inhibiting cyanobacterial growth and reducing the accumulation of microcystins. Besides this, the development of multi-cyanophage cocktails may have antagonistic effects on water blooms formed by multiple cyanobacterial outbreaks [42].  All assays were performed in triplicate, and the control group data for pH, temperature, and NaCl content were the same as those of the experimental group with pH of 7, temperature of 25 ℃, and NaCl content of 0%.

Genome Features
The cyanophage genome was sequenced (a 2 × 300 paired-end sequence method) on an Illumina MiSeq sequencing platform. MinS1 was found to be a GC-rich doublestranded DNA cyanophage, 49,966 bp long, with very low homology to known sequences within the NCBI database (1% homology coverage and 74.06% identity to the Myoviridae sp. ct6sY2). MinS1 had no obvious termini and no antibiotic tolerance genes, virulence genes or tRNAs.
A total of 75 ORFs of MinS1 were predicted using RAST [18], averaging 203 amino acid (aa) long, the smallest of which had only 24 aa and the largest as many as 1155 aa (Table 3). Contrary to evidence from cyanophage SH-68 (107 positive strands and 10 negative strands), MinS1 had just 5 positive strands but 70 negative strands, but MinS1 had the same coding density (91.8%) as that of M. elabens vB_MelS-Me-ZS1. Homology align- All assays were performed in triplicate, and the control group data for pH, temperature, and NaCl content were the same as those of the experimental group with pH of 7, temperature of 25 • C, and NaCl content of 0%.
A total of 75 ORFs of MinS1 were predicted using RAST [18], averaging 203 amino acid (aa) long, the smallest of which had only 24 aa and the largest as many as 1155 aa (Table 3). Contrary to evidence from cyanophage SH-68 (107 positive strands and 10 negative strands), MinS1 had just 5 positive strands but 70 negative strands, but MinS1 had the same coding density (91.8%) as that of M. elabens vB_MelS-Me-ZS1. Homology alignments of the protein sequences (nr) database were applied to the 75 ORFs' functional annotation; Viruses 2022, 14, 433 9 of 17 only 23 ORFs showed homology to genes with a known function, 39 were assigned to hypothetical proteins, and 13 had no homology to sequences within the database. This high rate of 69% of unknowns may be caused by the poor information known about the genomes of the cyanophages. Six functional modules (structure, replication, metabolism, regulation, packaging, and lysis) were divided among the MinS1 functionally annotated ORFs (Figure 4). The MinS1 genome therefore has similar functional modules to the other long-tailed phages and is capable of forming a complete phage independently.

Phylogenetic Analysis
To define the evolutionary status of newly isolated phages, sequences of TerL, cap proteins, and other aspects are often used for phylogenetic analysis. However, the div sity of phage genomes increases with the number of sequenced sequences; hence, it is n accurate to simply classify them based on morphological features or single-gene evo tionary lineages alone. Whole-genome-based proteome phylogeny analysis is increasin favored by researchers. The inter-genomic similarity calculated by VIRIDIC showed lit similarity between MinS1 and all of the other known Microcystis cyanophages (Figu 5b). Thirty representative phage strains from the nine families belonging to Caudovira The outermost circle represents 75 ORFs encoded within the genome, with different colors representing different functions (clockwise arrow indicates the forward reading frame, and counterclockwise arrow indicates the reverse reading frame); the gray circles in the middle represent the GC content (outwards indicates greater than the average GC content compared with the whole genome, and inwards indicates the opposite); the innermost circle represents the GC skew (G − C/G + C. Outwards indicates >0, and inwards indicates <0).
Lysis and packaging genes. ORF38 and ORF65 are classified into lysis and packaging clusters, encoding N-acetylmuramoyl-L-alanine amidase and the terminase large subunit (TerL), respectively. In general, lysozyme homologs (e.g., lysin, holin, endolysin, and others) are commonly found in phages and are thought to be responsible for cell lysis [43][44][45]. However, the above lytic enzymes were not found, but N-acetylmuramoyl-L-alanine amidase was annotated within MinS1. In enzymology, N-acetylmuramoyl-L-alanine amidase (which belongs to the hydrolase family) is able to cleave the amide bond between N-acetylmuramoyl and L-amino acids in bacterial cell walls. Autolysins and some phage lysins are examples of N-acetylmuramoyl-L-alanine amidases [46]. Thus, ORF38, which encodes N-acetylmuramoyl and L-amino acid protein, may be the key gene responsible for lysis. This suggests that MinS1 adopts a different lysis strategy to others previously discovered, which may also be related to its broad-spectrum capabilities. TerL is commonly found in phages and is responsible for phage assembly. Due to its conserved nature, it is often used to aid evolutionary analysis.
Structure and regulation genes. The structure collection consists of ORF17, ORF39, ORF46, and ORF53, which encode gp121, collagen-like protein, tail tape measure protein (TMP), and gp10, respectively. Among these proteins, TMP determines the tail length and allows phage genes to transfer into the host [47]; however, this protein is not unique to tailed phages and has also been observed in tailless phages [48]. The regulation gene cassette consists of four ORFs, including AAA ATPase (ORF4), kinase domain protein (ORF19), site-specific integrase (ORF29), and tyrosine-type recombinase/integrase (ORF31). The AAA family proteins often perform chaperone-like functions, assisting in the assembly, operation, or disassembly of protein complexes [49].
Replication and metabolism genes. From an evolutionary perspective, phages have undergone multiple gene exchange events in response to selection pressure from their hosts, which has, in turn, driven their diversity [50]. Here, two ORFs associated with gene transfer were found in the MinS1 genome, HNH endonuclease (ORF62) and the GIY-YIG nuclease family protein (ORF66) [51][52][53]. They play a key role in relation to phages acquiring new genes when competing with specific bacteria, which facilitates survival adaptation [54]. The remaining ORFs in this block focus on nuclease activities, such as the endonuclease/exonuclease/phosphatase family protein (ORF37), which is involved in cellular signaling [55]; 3 -5 exonuclease (ORF16), which controls the 3 -5 exonuclease activity of DNA polymerase I and other enzymes and catalyzes mismatched nucleotide hydrolysis; the Lsr2 family protein (ORF7) and helix-turn-helix domain-containing protein (ORF24, ORF25, and ORF34), which are capable of multivariate regulation of gene expression and metabolism; and the RusA family crossover junction endodeoxyribonuclease (ORF9), which can correct DNA repair defects.

Phylogenetic Analysis
To define the evolutionary status of newly isolated phages, sequences of TerL, capsid proteins, and other aspects are often used for phylogenetic analysis. However, the diversity of phage genomes increases with the number of sequenced sequences; hence, it is not accurate to simply classify them based on morphological features or single-gene evolutionary lineages alone. Whole-genome-based proteome phylogeny analysis is increasingly favored by researchers. The inter-genomic similarity calculated by VIRIDIC showed little similarity between MinS1 and all of the other known Microcystis cyanophages (Figure 5b). Thirty representative phage strains from the nine families belonging to Caudovirales in ICTV were selected, and a proteomic phylogenetic tree was constructed based on the MinS1 whole-genome sequence similarity. The proteomic tree (Figure 5a) showed that MinS1 was assigned into the cluster of Siphoviridae phages and was in the same branch with the Brevibacterium phage Lucky Barnes.

Conclusions
The present study described the characteristics and genome of MinS1, a novel strain of Siphoviridae freshwater cyanophage, with most of the predicted protein-coding genes showing no significant similarity to sequences within published databases. Its possible broad spectrum of genetic factors was explored through genomic analysis. In addition, MinS1 exhibited a broader host range than that of other known cyanophages and was tolerant to temperature, pH, UV, and salinity, suggesting that MinS1 has good potential for application as a biological control agent against cyanobacterial blooms. This research highlights our understanding of cyanophage biological characteristics, and it indicates that it has good potential for developing applications against water blooms caused by multiple cyanobacterial blooms.

Conclusions
The present study described the characteristics and genome of MinS1, a novel strain of Siphoviridae freshwater cyanophage, with most of the predicted protein-coding genes showing no significant similarity to sequences within published databases. Its possible broad spectrum of genetic factors was explored through genomic analysis. In addition, MinS1 exhibited a broader host range than that of other known cyanophages and was tolerant to temperature, pH, UV, and salinity, suggesting that MinS1 has good potential for application as a biological control agent against cyanobacterial blooms. This research highlights our understanding of cyanophage biological characteristics, and it indicates that it has good potential for developing applications against water blooms caused by multiple cyanobacterial blooms.

Supplementary Materials:
The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/v14020433/s1, Figure S1: (a) Phenotypic photographs of cyanobacterial cultures within the host range experiments. (b) Phage spots formed by the cyanophage MinS1 on susceptible cyanobacteria strains; Figure S2: Status of host cyanobacteria infected by MinS1 treated under different conditions. Table S1: MinS1 classification, general features, and genome sequencing information.