Comparative Genomic and Phylogenomic Analyses Reveal a Conserved Core Genome Shared by Estuarine and Oceanic Cyanopodoviruses

Genome sequences of the first Autographiviridae phages infecting marine Roseobacter

The ubiquitous and abundant marine phages play critical roles in shaping the composition and function of bacterial communities, impacting biogeochemical cycling in marine ecosystems. Autographiviridae is among the most abundant and ubiquitous phage families in the ocean. However, studies on the diversity and ecology of Autographiviridae phages in marine environments are restricted to isolates that infect SAR11 bacteria and cyanobacteria. In this study, ten new roseophages that infect marine Roseobacter strains were isolated from coastal waters. These new roseophages have a genome size ranging from 38 917 to 42 634 bp and G+C content of 44.6-50 %. Comparative genomics showed that they are similar to known Autographiviridae phages regarding gene content and architecture, thus representing the first Autographiviridae roseophages. Phylogenomic analysis based on concatenated conserved genes showed that the ten roseophages form three distinct subgroups within the Autographiviridae, and sequence analysis revealed that they belong to eight new genera. Finally, viromic read-mapping showed that these new Autographiviridae phages are widely distributed in global oceans, mostly inhabiting polar and estuarine locations. This study has expanded the current understanding of the genomic diversity, evolution and ecology of Autographiviridae phages and roseophages. We suggest that Autographiviridae phages play important roles in the mortality and community structure of roseobacters, and have broad ecological applications.

Cryo-EM structure of cyanophage P-SCSP1u offers insights into DNA gating and evolution of T7-like viruses

Cyanophages, together with their host cyanobacteria, play important roles in marine biogeochemical cycles and control of marine food webs. The recently identified MPP-C (Marine Picocyanobacteria Podovirus clade C) cyanophages, belonging to the T7-like podoviruses, contain the smallest genomes among cyanopodoviruses and exhibit distinct infection kinetics. However, understanding of the MPP-C cyanophage infection process is hindered by the lack of high-resolution structural information. Here, we report the cryo-EM structure of the cyanophage P-SCSP1u, a representative member of the MPP-C phages, in its native form at near-atomic resolution, which reveals the assembly mechanism of the capsid and molecular interaction of the portal-tail complex. Structural comparison of the capsid proteins of P-SCSP1u and other podoviruses with known structures provides insights into the evolution of T7-like viruses. Furthermore, our study provides the near-atomic resolution structure of portal-tail complex for T7-like viruses. On the basis of previously reported structures of phage T7, we identify an additional valve and gate to explain the DNA gating mechanism for the T7-like viruses.

The Impact of Viral Infection on the Chemistries of the Earth's Most Abundant Photosynthesizes: Metabolically Talented Aquatic Cyanobacteria

Cyanobacteria are the most abundant photosynthesizers on earth, and as such, they play a central role in marine metabolite generation, ocean nutrient cycling, and the control of planetary oxygen generation. Cyanobacteriophage infection exerts control on all of these critical processes of the planet, with the phage-ported homologs of genes linked to photosynthesis, catabolism, and secondary metabolism (marine metabolite generation). Here, we analyze the 153 fully sequenced cyanophages from the National Center for Biotechnology Information (NCBI) database and the 45 auxiliary metabolic genes (AMGs) that they deliver into their hosts. Most of these AMGs are homologs of those found within cyanobacteria and play a key role in cyanobacterial metabolism-encoding proteins involved in photosynthesis, central carbon metabolism, phosphate metabolism, methylation, and cellular regulation. A greater understanding of cyanobacteriophage infection will pave the way to a better understanding of carbon fixation and nutrient cycling, as well as provide new tools for synthetic biology and alternative approaches for the use of cyanobacteria in biotechnology and sustainable manufacturing.

Multi-year molecular quantification and 'omics analysis of Planktothrix-specific cyanophage sequences from Sandusky Bay, Lake Erie

Introduction

Planktothrix agardhii is a microcystin-producing cyanobacterium found in Sandusky Bay, a shallow and turbid embayment of Lake Erie. Previous work in other systems has indicated that cyanophages are an important natural control factor of harmful algal blooms. Currently, there are few cyanophages that are known to infect P. agardhii, with the best-known being PaV-LD, a tail-less cyanophage isolated from Lake Donghu, China. Presented here is a molecular characterization of Planktothrix specific cyanophages in Sandusky Bay.

Methods and Results

Putative Planktothrix-specific viral sequences from metagenomic data from the bay in 2013, 2018, and 2019 were identified by two approaches: homology to known phage PaV-LD, or through matching CRISPR spacer sequences with Planktothrix host genomes. Several contigs were identified as having viral signatures, either related to PaV-LD or potentially novel sequences. Transcriptomic data from 2015, 2018, and 2019 were also employed for the further identification of cyanophages, as well as gene expression of select viral sequences. Finally, viral quantification was tested using qPCR in 2015-2019 for PaV-LD like cyanophages to identify the relationship between presence and gene expression of these cyanophages. Notably, while PaV-LD like cyanophages were in high abundance over the course of multiple years (qPCR), transcriptomic analysis revealed only low levels of viral gene expression.

Discussion

This work aims to provide a broader understanding of Planktothrix cyanophage diversity with the goals of teasing apart the role of cyanophages in the control and regulation of harmful algal blooms and designing monitoring methodology for potential toxin-releasing lysis events.

Genomic Analysis of a New Freshwater Cyanophage Lbo240-yong1 Suggests a New Taxonomic Family of Bacteriophages

A worldwide ecological issue, cyanobacterial blooms in marine and freshwater have caused enormous losses in both the economy and the environment. Virulent cyanophages-specifically, infecting and lysing cyanobacteria-are key ecological factors involved in limiting the overall extent of the population development of cyanobacteria. Over the past three decades, reports have mainly focused on marine Prochlorococcus and Synechococcus cyanophages, while information on freshwater cyanophages remained largely unknown. In this study, a novel freshwater cyanophage, named Lbo240-yong1, was isolated via the double-layer agar plate method using Leptolyngbya boryana FACHB-240 as a host. Transmission electron microscopy observation illustrated the icosahedral head (50 ± 5 nm in diameter) and short tail (20 ± 5 nm in length) of Lbo240-yong1. Experimental infection against 37 cyanobacterial strains revealed that host-strain-specific Lbo240-yong1 could only lyse FACHB-240. The complete genome of Lbo240-yong1 is a double-stranded DNA of 39,740 bp with a G+C content of 51.99%, and it harbors 44 predicted open reading frames (ORFs). A Lbo240-yong1 ORF shared the highest identity with a gene of a filamentous cyanobacterium, hinting at a gene exchange between the cyanophage and cyanobacteria. A BLASTn search illustrated that Lbo240-yong1 had the highest sequence similarity with the Phormidium cyanophage Pf-WMP4 (89.67% identity, 84% query coverage). In the proteomic tree based on genome-wide sequence similarities, Lbo240-yong1, three Phormidium cyanophages (Pf-WMP4, Pf-WMP3, and PP), one Anabaena phage (A-4L), and one unclassified Arthronema cyanophage (Aa-TR020) formed a monophyletic group that was more deeply diverging than several other families. Pf-WMP4 is the only member of the independent genus Wumpquatrovirus that belongs to the Caudovircetes class. Pf-WMP3 and PP formed the independent genus Wumptrevirus. Anabaena phage A-4L is the only member of the independent Kozyakovvirus genus. The six cyanopodoviruses share similar gene arrangements. Eight core genes were found in them. We propose, here, to set up a new taxonomic family comprising the six freshwater cyanopodoviruses infecting filamentous cyanobacteria. This study enriched the field's knowledge of freshwater cyanophages.

Transcriptomics-Driven Characterization of LUZ100, a T7-like Pseudomonas Phage with Temperate Features

Autographiviridae is a diverse yet distinct family of bacterial viruses marked by a strictly lytic lifestyle and a generally conserved genome organization. Here, we characterized Pseudomonas aeruginosa phage LUZ100, a distant relative of type phage T7. LUZ100 is a podovirus with a limited host range which likely uses lipopolysaccharide (LPS) as a phage receptor. Interestingly, infection dynamics of LUZ100 indicated moderate adsorption rates and low virulence, hinting at temperate characteristics. This hypothesis was supported by genomic analysis, which showed that LUZ100 shares the conventional T7-like genome organization yet carries key genes associated with a temperate lifestyle. To unravel the peculiar characteristics of LUZ100, ONT-cappable-seq transcriptomics analysis was performed. These data provided a bird's-eye view of the LUZ100 transcriptome and enabled the discovery of key regulatory elements, antisense RNA, and transcriptional unit structures. The transcriptional map of LUZ100 also allowed us to identify new RNA polymerase (RNAP)-promoter pairs that can form the basis for biotechnological parts and tools for new synthetic transcription regulation circuitry. The ONT-cappable-seq data revealed that the LUZ100 integrase and a MarR-like regulator (proposed to be involved in the lytic/lysogeny decision) are actively cotranscribed in an operon. In addition, the presence of a phage-specific promoter transcribing the phage-encoded RNA polymerase raises questions on the regulation of this polymerase and suggests that it is interwoven with the MarR-based regulation. This transcriptomics-driven characterization of LUZ100 supports recent evidence that T7-like phages should not automatically be assumed to have a strictly lytic life cycle. IMPORTANCE Bacteriophage T7, considered the "model phage" of the Autographiviridae family, is marked by a strictly lytic life cycle and conserved genome organization. Recently, novel phages within this clade have emerged which display characteristics associated with a temperate life cycle. Screening for temperate behavior is of utmost importance in fields like phage therapy, where strictly lytic phages are generally required for therapeutic applications. In this study, we applied an omics-driven approach to characterize the T7-like Pseudomonas aeruginosa phage LUZ100. These results led to the identification of actively transcribed lysogeny-associated genes in the phage genome, pointing out that temperate T7-like phages are emerging more frequent than initially thought. In short, the combination of genomics and transcriptomics allowed us to obtain a better understanding of the biology of nonmodel Autographiviridae phages, which can be used to optimize the implementation of phages and their regulatory elements in phage therapy and biotechnological applications, respectively.

Abundant and cosmopolitan lineage of cyanopodoviruses lacking a DNA polymerase gene

Cyanopodoviruses affect the mortality and population dynamics of the unicellular picocyanobacteria Prochlorococcus and Synechococcus, the dominant primary producers in the oceans. Known cyanopodoviruses all contain the DNA polymerase gene (DNA pol) that is important for phage DNA replication and widely used in field quantification and diversity studies. However, we isolated 18 cyanopodoviruses without identifiable DNA pol. They form a new MPP-C clade that was separated from the existing MPP-A, MPP-B, and P-RSP2 clades. The MPP-C phages have the smallest genomes (37.3-37.9 kb) among sequenced cyanophages, and show longer latent periods than the MPP-B phages. Metagenomic reads of both clades are highly abundant in surface waters, but the MPP-C phages show higher relative abundance in surface waters than in deeper waters, while MPP-B phages have higher relative abundance in deeper waters. Our study reveals that cyanophages with distinct genomic contents and infection kinetics can exhibit different depth profiles in the oceans.

Widespread yet persistent low abundance of TIM5-like cyanophages in the oceans

Ocean ecosystems are inhabited by a diverse set of viruses that impact microbial mortality and evolution. However, the distribution and abundances of specific viral lineages, particularly those from the large bank of rare viruses, remains largely unknown. Here, we assessed the diversity and abundance of the TIM5-like cyanophages. The sequencing of three new TIM5-like cyanophage genomes and environmental amplicons of a signature gene from the Red Sea revealed highly conserved gene content and sequence similarity. We adapted the polony method, a solid-phase polymerase chain reaction assay, to quantify TIM5-like cyanophages during three 2000 km expeditions in the Pacific Ocean and four annual cycles in the Red Sea. TIM5-like cyanophages were widespread, detected at all latitudes and seasons surveyed throughout the photic zone. Yet they were generally rare, ranging between <100 and 4000 viruses·ml^-1 . Occasional peaks in abundance of 10- to 100-fold were observed, reaching 71,000 viruses·ml^-1 . These peaks were ephemeral and seasonally variable in the Red Sea. Infection levels, quantified during one such peak, were very low. These characteristics of low diversity and abundance, as well as variable outbreaks, distinguishes the TIM5-like lineage from other major cyanophage lineages and illuminates that rare virus lineages can be persistent and widespread in the oceans.

A Cyanophage MarR-Type Transcription Factor Regulates Host RNase E Expression during Infection

The marine picocyanobacterium Prochlorococcus contributes significantly to global primary production, and its abundance and diversity is shaped in part by viral infection. Here, we identified a cyanophage-encoded MarR-type transcription factor that induces the gene expression of host Prochlorococcus MED4 endoribonuclease (RNase) E during phage infection. The increase in rne transcript levels relies on the phage (p)MarR-mediated activation of an alternative promoter that gives rise to a truncated yet enzymatically fully functional RNase E isoform. In this study, we demonstrate that pMarR binds to an atypical activator site downstream of the transcriptional start site and that binding is enhanced in the presence of Ca2+ ions. Furthermore, we show that dimeric pMarR interacts with the α subunit of RNA polymerase, and we identified amino acid residues S66, R67, and G106, which are important for Ca2+ binding, DNA binding, and dimerization of pMarR, respectively.

A Unique Set of Auxiliary Metabolic Genes Found in an Isolated Cyanophage Sheds New Light on Marine Phage-Host Interactions

Cyanophages, viruses that infect cyanobacteria, are abundant and widely distributed in aquatic ecosystems, playing important roles in regulating the abundance, activity, diversity, and evolution of cyanobacteria. A T4-like cyanophage, S-SCSM1, infecting Synechococcus and Prochlorococcus strains of different ecotypes, was isolated from the South China Sea in this study. For the first time, a mannose-6-phosphate isomerase (MPI) gene was identified in the cultured cyanophage. At least 11 phylogenetic clusters of cyanophage MPIs were retrieved and identified from the marine metagenomic data sets, indicating that cyanophage MPIs in the marine environment are extremely diverse. The existence of 24 genes encoding 2-oxoglutarate (2OG)-Fe(II) oxygenase superfamily proteins in the S-SCSM1 genome emphasizes their potential importance and diverse functions in reprogramming host metabolism during phage infection. Novel cell wall synthesis and modification genes found in the S-SCSM1 genome indicate that diverse phenotypic modifications imposed by phages on cyanobacterial hosts remain to be discovered. Two noncoding RNAs of cis-regulatory elements in the S-SCSM1 genome were predicted to be associated with host exopolysaccharide metabolism and photosynthesis. The isolation and genomic characterization of cyanophage S-SCSM1 provide more information on the genetic diversity of cyanophages and phage-host interactions in the marine environment. IMPORTANCE Cyanophages play important ecological roles in aquatic ecosystems. Genomic and proteomic characterizations of the T4-like cyanophage S-SCSM1 indicate that novel and diverse viral genes and phage-host interactions in the marine environment remain unexplored. The first identified mannose-6-phosphate isomerase (MPI) gene from a cultured cyanophage was found in the S-SCSM1 genome, although MPIs were previously found in viral metagenomes at high frequencies similar to those of the cyanophage photosynthetic gene psbA. The presence of 24 genes encoding 2-oxoglutarate (2OG)-Fe(II) oxygenase superfamily proteins, novel cell wall synthesis and modification genes, a nonbleaching protein A gene, and 2 noncoding RNAs of cis-regulatory elements in the S-SCSM1 genome as well as the presence of a virion-associated regulatory protein indicate the diverse functions that cyanophages have in reprogramming the metabolism and modifying the phenotypes of hosts during infection.

Cyanophages from a less virulent clade dominate over their sister clade in global oceans

Environmental virus communities are highly diverse. However, the infection physiology underlying the evolution of diverse phage lineages and their ecological consequences are largely unknown. T7-like cyanophages are abundant in nature and infect the marine unicellular cyanobacteria, Synechococcus and Prochlorococcus, important primary producers in the oceans. Viruses belonging to this genus are divided into two distinct phylogenetic clades: clade A and clade B. These viruses have narrow host-ranges with clade A phages primarily infecting Synechococcus genotypes, while clade B phages are more diverse and can infect either Synechococcus or Prochlorococcus genotypes. Here we investigated infection properties (life history traits) and environmental abundances of these two clades of T7-like cyanophages. We show that clade A cyanophages have more rapid infection dynamics, larger burst sizes and greater virulence than clade B cyanophages. However, clade B cyanophages were at least 10-fold more abundant in all seasons, and infected more cyanobacteria, than clade A cyanophages in the Red Sea. Models predicted that steady-state cyanophage abundances, infection frequency, and virus-induced mortality, peak at intermediate virulence values. Our findings indicate that differences in infection properties are reflected in virus phylogeny at the clade level. They further indicate that infection properties, together with differences in subclade diversity and host repertoire, have important ecological consequences with the less aggressive, more diverse virus clade having greater ecological impacts.

Genomic Characterization of a Novel Freshwater Cyanophage Reveals a New Lineage of Cyanopodovirus

Cyanobacteria are one of the dominant autotrophs in tropical freshwater communities, yet phages infecting them remain poorly characterized. Here we present the characterization of cyanophage S-SRP02, isolated from a tropical freshwater lake in Singapore, which infects Synechococcus sp. Strain SR-C1 isolated from the same lake. S-SRP02 represents a new evolutionary lineage of cyanophage. Out of 47 open reading frames (ORFs), only 20 ORFs share homology with genes encoding proteins of known function. There is lack of auxiliary metabolic genes which was commonly found as core genes in marine cyanopodoviruses. S-SRP02 also harbors unique structural genes highly divergent from other cultured phages. Phylogenetic analysis and viral proteomic tree further demonstrate the divergence of S-SRP02 from other sequenced phage isolates. Nonetheless, S-SRP02 shares synteny with phage genes of uncultured phages obtained from the Mediterranean Sea deep chlorophyll maximum fosmids, indicating the ecological importance of S-SRP02 and its related viruses. This is further supported by metagenomic mapping of environmental viral metagenomic reads onto the S-SRP02 genome.

Novel Freshwater Cyanophages Provide New Insights into Evolutionary Relationships between Freshwater and Marine Cyanophages

Cyanobacteria and cyanophages are present widely in both freshwater and marine environments. However, freshwater cyanophages remain unknown largely due to the small numbers of cyanophage isolates despite their ecological and environmental significance. In this study, we present the characterization of two novel lytic freshwater cyanophages isolated from a tropical inland lake in Singapore, namely, cyanopodovirus S-SRP01 and cyanomyovirus S-SRM01, infecting two different strains of Synechococcus spp. Functional annotation of S-SRP01 and S-SRM01 genomes revealed a high degree of homology with marine cyanophages. Phylogenetic trees of concatenated genes and whole-genome alignment provided further evidence that S-SRP01 is close evolutionarily to marine cyanopodoviruses, while S-SRM01 is evolutionarily close to marine cyanomyoviruses. Few genetic similarities between freshwater and marine cyanophages have been identified in previous studies. The isolation of S-SRP01 and S-SRM01 expand current knowledge on freshwater cyanophages infecting Synechococcus spp. Their high degree of gene sharing provides new insights into the evolutionary relationships between freshwater and marine cyanophages. This relatedness is further supported by the discovery of similar phenomenon from other freshwater viral metagenomes. IMPORTANCE This study expands the current knowledge on freshwater cyanophage isolates and cyanophage genetic diversity, indicating that freshwater and marine cyanophages infecting Synechococcus spp. may share close genetic similarity and evolutionary relationships.

Phylogenetic systematics of Butyrivibrio and Pseudobutyrivibrio genomes illustrate vast taxonomic diversity, open genomes and an abundance of carbohydrate-active enzyme family isoforms

Butyrivibrio and Pseudobutyrivibrio dominate in anaerobic gastrointestinal microbiomes, particularly the rumen, where they play a key role in harvesting dietary energy. Within these genera, five rumen species have been classified (Butyrivibrio fibrisolvens, Butyrivibrio hungatei, Butyrivibrio proteoclasticus, Pseudobutyrivibrio ruminis and Pseudobutyrivibrio xylanivorans) and more recently an additional Butyrivibrio sp. group was added. Given the recent increase in available genomes, we re-investigated the phylogenetic systematics and evolution of Butyrivibrio and Pseudobutyrivibrio. Across 71 genomes, we show using 16S rDNA and 40 gene marker phylogenetic trees that the current six species designations (P. ruminis, P. xylanivorans, B. fibrisolvens, Butyrivibrio sp., B. hungatei and B. proteclasticus) are found. However, pangenome analysis showed vast genomic variation and a high abundance of accessory genes (91.50–99.34 %), compared with core genes (0.66–8.50 %), within these six taxonomic groups, suggesting incorrectly assigned taxonomy. Subsequent pangenome accessory genomes under varying core gene cut-offs (%) and average nucleotide identity (ANI) analysis suggest the existence of 42 species within 32 genera. Pangenome analysis of those that still group within B. fibrisolvens, B. hungatei and P. ruminis, based on revised ANI phylogeny, also showed possession of very open genomes, illustrating the diversity that exists even within these groups. All strains of both Butyrivibrio and Pseudobutyrivibrio also shared a broad range of clusters of orthologous genes (COGs) (870), indicating recent evolution from a common ancestor. We also demonstrate that the carbohydrate-active enzymes (CAZymes) predominantly belong to glycosyl hydrolase (GH)2, 3, 5, 13 and 43, with numerous within family isoforms apparent, likely facilitating metabolic plasticity and resilience under dietary perturbations. This study provides a major advancement in our functional and evolutionary understanding of these important anaerobic bacteria.

Pseudomonas Phage MD8: Genetic Mosaicism and Challenges of Taxonomic Classification of Lambdoid Bacteriophages

Pseudomonas phage MD8 is a temperate phage isolated from the freshwater lake Baikal. The organisation of the MD8 genome resembles the genomes of lambdoid bacteriophages. However, MD8 gene and protein sequences have little in common with classified representatives of lambda-like phages. Analysis of phage genomes revealed a group of other Pseudomonas phages related to phage MD8 and the genomic layout of MD8-like phages indicated extensive gene exchange involving even the most conservative proteins and leading to a high degree of genomic mosaicism. Multiple horizontal transfers and mosaicism of the genome of MD8, related phages and other λ-like phages raise questions about the principles of taxonomic classification of the representatives of this voluminous phage group. Comparison and analysis of various bioinformatic approaches applied to λ-like phage genomes demonstrated different efficiency and contradictory results in the estimation of genomic similarity and relatedness. However, we were able to make suggestions for the possible origin of the MD8 genome and the basic principles for the taxonomic classification of lambdoid phages. The group comprising 26 MD8-related phages was proposed to classify as two close genera belonging to a big family of λ-like phages.

Genomic diversity, life strategies and ecology of marine HTVC010P-type pelagiphages

SAR11 bacteria dominate ocean surface bacterioplankton communities, and play an important role in marine carbon and nutrient cycling. The biology and ecology of SAR11 are impacted by SAR11 phages (pelagiphages) that are highly diverse and abundant in the ocean. Among the currently known pelagiphages, HTVC010P represents an extremely abundant but under-studied phage group in the ocean. In this study, we have isolated seven new HTVC010P-type pelagiphages, and recovered 77 nearly full-length HTVC010P-type metagenomic viral genomes from marine metagenomes. Comparative genomic and phylogenomic analyses showed that HTVC010P-type pelagiphages display genome synteny and can be clustered into two major subgroups, with subgroup I consisting of strictly lytic phages and subgroup II mostly consisting of phages with potential lysogenic life cycles. All but one member of the subgroup II contain an integrase gene. Site-specific integration of subgroup II HTVC010P-type pelagiphage was either verified experimentally or identified by in silico genomic sequence analyses, which revealed that various SAR11 tRNA genes can serve as the integration sites of HTVC010P-type pelagiphages. Moreover, HTVC010P-type pelagiphage integration was confirmed by the detection of several Global Ocean Survey (GOS) fragments that contain hybrid phage–host integration sites. Metagenomic recruitment analysis revealed that these HTVC010P-type phages were globally distributed and most lytic subgroup I members exhibited higher relative abundance. Altogether, this study significantly expands our knowledge about the genetic diversity, life strategies and ecology of HTVC010P-type pelagiphages.

The genomic content and context of auxiliary metabolic genes in roseophages

Marine bacteriophages frequently possess auxiliary metabolic genes (AMGs) that accelerate host metabolism during phage infection. The significance of AMGs in phage infecting the ecologically important Roseobacter clade, found predominantly in marine environments, remains to be determined. Here, we analysed the distribution and genomic context of 180 AMGs, annotated into 20 types, across 50 roseophage genomes. Roseophages share seven high-frequency AMGs (trx, grx, RNR, thyX, DCD, phoH, and mazG)， most of them involved in the nucleotide biosynthesis pathway that represent conserved intra and inter operational taxonomic units (OTUs), and share ≥97% full-length DNA sequence similarity. Sporadic AMGs (dUTPase, lexA, degS, Que, NAPRT, AHL, pcnB, ctrA, RTX, RNR-nrdA, RNR-nrdE, wclP, and flgJ), present in only one or two OTUs, show high functional diversity. The roseophage AMG repertoire weakly correlates with environmental factors, while host range partially explains the sporadic AMG distribution. Locally co-linear blocks distribution index (LDI) analysis indicated that high-frequency roseopodovirus AMGs are restricted to particular genomic islands, possibly originating from limited historical acquisition events. Low-frequency roseopodovirus AMGs and all roseosiphovirus AMGs have high LDI values, implying multiple historical acquisition events. In summary, roseophages have acquired a range of AMGs through horizontal gene transfer, and the forces shaping the evolution of roseophages are described.

Novel Virus on Filamentous Arthronema africanum Cyanobacterium

Widely distributed in water environments and in soil, cyanobacteria are hosts of lysogenic or lytic bacterioviruses. A novel, probably lysogenic virus (phage) for which the name Arthronema africanum virus TR020 (Aa-TR020) is proposed, has been isolated from filamentous freshwater cyanobacterium Arthronema africanum. The virus formed turbid plaques on plate culture of A. africanum strain 1980/01 but not on other Arthronema strain and other bacterial species. The genome of Aa-TR020 is linear molecule of dsDNA, 44,805 bp in length with 216 bp long terminal repeats and with G + C content of 46%. Fifty-five genes organized on plus and minus strands were predicted there. The genome size, gene arrangement, and selected protein sequences showed relatedness to Phormidium virus Pf-WMP3 and other viruses known to infect cyanobacteria and classified in the family Podoviridae.

Frequency of mispackaging of Prochlorococcus DNA by cyanophage

Prochlorococcus cells are the numerically dominant phototrophs in the open ocean. Cyanophages that infect them are a notable fraction of the total viral population in the euphotic zone, and, as vehicles of horizontal gene transfer, appear to drive their evolution. Here we examine the propensity of three cyanophages-a podovirus, a siphovirus, and a myovirus-to mispackage host DNA in their capsids while infecting Prochlorococcus, the first step in phage-mediated horizontal gene transfer. We find the mispackaging frequencies are distinctly different among the three phages. Myoviruses mispackage host DNA at low and seemingly fixed frequencies, while podo- and siphoviruses vary in their mispackaging frequencies by orders of magnitude depending on growth light intensity. We link this difference to the concentration of intracellular reactive oxygen species and protein synthesis rates, both parameters increasing in response to higher light intensity. Based on our findings, we propose a model of mispackaging frequency determined by the imbalance between the production of capsids and the number of phage genome copies during infection: when protein synthesis rate increase to levels that the phage cannot regulate, they lead to an accumulation of empty capsids, in turn triggering more frequent host DNA mispackaging errors.

Temporal transcriptomes of a marine cyanopodovirus and its Synechococcus host during infection

Marine picocyanobacteria belonging to genera Synechococcus and Prochlorococcus are genetically diverged and distributed into distinct biogeographical patterns, and both are infected by genetically closely related cyanopodoviruses. Previous studies have not fully explored whether the two virus–host systems share similar gene expression patterns during infection. Whole‐genome expression dynamics of T7‐like cyanopodovirus P‐SSP7 and its host Prochlorococcus strain MED4 have already been reported. Here, we conducted genomic and transcriptomic analyses on T7‐like cyanopodovirus S‐SBP1 during its infection on Synechococcus strain WH7803. S‐SBP1 has a latent period of 8 h and phage DNA production of 30 copies per cell. In terms of whole‐genome phylogenetic relationships and average nucleotide identity, S‐SBP1 was most similar to cyanopodovirus S‐RIP2, which also infects Synechococcus WH7803. Three hypervariable genomic islands were identified when comparing the genomes of S‐SBP1 and S‐RIP2. Single nucleotide variants were also observed in three S‐SBP1 genes, which were located within the island regions. Based on RNA‐seq analysis, S‐SBP1 genes clustered into three temporal expression classes, whose gene content was similar to that of P‐SSP7. Thirty‐two host genes were upregulated during phage infection, including those involved in carbon metabolism, ribosome components, and stress response. These upregulated genes were similar to those upregulated by Prochlorococcus MED4 in response to infection by P‐SSP7. Our study demonstrates a programmed temporal expression pattern of cyanopodoviruses and hosts during infection.

Bioinformatics as a Tool for the Structural and Evolutionary Analysis of Proteins

This chapter deals with the topic of bioinformatics, computational, mathematics, and statistics tools applied to biology, essential for the analysis and characterization of biological molecules, in particular proteins, which play an important role in all cellular and evolutionary processes of the organisms. In recent decades, with the next generation sequencing technologies and bioinformatics, it has facilitated the collection and analysis of a large amount of genomic, transcriptomic, proteomic, and metabolomic data from different organisms that have allowed predictions on the regulation of expression, transcription, translation, structure, and mechanisms of action of proteins as well as homology, mutations, and evolutionary processes that generate structural and functional changes over time. Although the information in the databases is greater every day, all bioinformatics tools continue to be constantly modified to improve performance that leads to more accurate predictions regarding protein functionality, which is why bioinformatics research remains a great challenge.

Exploration into the origins and mobilization of di-hydrofolate reductase genes and the emergence of clinical resistance to trimethoprim

Trimethoprim is a synthetic antibacterial agent that targets folate biosynthesis by competitively binding to the di-hydrofolate reductase enzyme (DHFR). Trimethoprim is often administered synergistically with sulfonamide, another chemotherapeutic agent targeting the di-hydropteroate synthase (DHPS) enzyme in the same pathway. Clinical resistance to both drugs is widespread and mediated by enzyme variants capable of performing their biological function without binding to these drugs. These mutant enzymes were assumed to have arisen after the discovery of these synthetic drugs, but recent work has shown that genes conferring resistance to sulfonamide were present in the bacterial pangenome millions of years ago. Here, we apply phylogenetics and comparative genomics methods to study the largest family of mobile trimethoprim-resistance genes (dfrA). We show that most of the dfrA genes identified to date map to two large clades that likely arose from independent mobilization events. In contrast to sulfonamide resistance (sul) genes, we find evidence of recurrent mobilization in dfrA genes. Phylogenetic evidence allows us to identify novel dfrA genes in the emerging pathogen Acinetobacter baumannii, and we confirm their resistance phenotype in vitro. We also identify a cluster of dfrA homologues in cryptic plasmid and phage genomes, but we show that these enzymes do not confer resistance to trimethoprim. Our methods also allow us to pinpoint the chromosomal origin of previously reported dfrA genes, and we show that many of these ancient chromosomal genes also confer resistance to trimethoprim. Our work reveals that trimethoprim resistance predated the clinical use of this chemotherapeutic agent, but that novel mutations have likely also arisen and become mobilized following its widespread use within and outside the clinic. Hence, this work confirms that resistance to novel drugs may already be present in the bacterial pangenome, and stresses the importance of rapid mobilization as a fundamental element in the emergence and global spread of resistance determinants.

Origin and Evolution of Studiervirinae Bacteriophages Infecting Pectobacterium: Horizontal Transfer Assists Adaptation to New Niches

Black leg and soft rot are devastating diseases causing up to 50% loss of potential potato yield. The search for, and characterization of, bacterial viruses (bacteriophages) suitable for the control of these diseases is currently a sought-after task for agricultural microbiology. Isolated lytic Pectobacterium bacteriophages Q19, PP47 and PP81 possess a similar broad host range but differ in their genomic properties. The genomic features of characterized phages have been described and compared to other Studiervirinae bacteriophages. Thorough phylogenetic analysis has clarified the taxonomy of the phages and their positioning relative to other genera of the Autographiviridae family. Pectobacterium phage Q19 seems to represent a new genus not described previously. The genomes of the phages are generally similar to the genome of phage T7 of the Teseptimavirus genus but possess a number of specific features. Examination of the structure of the genes and proteins of the phages, including the tail spike protein, underlines the important role of horizontal gene exchange in the evolution of these phages, assisting their adaptation to Pectobacterium hosts. The results provide the basis for the development of bacteriophage-based biocontrol of potato soft rot as an alternative to the use of antibiotics.

Genomic analysis of Synechococcus phage S-B43 and its adaption to the coastal environment

Synechococcus dominate picocyanobacterial communities in coastal environments. However, only a few Synechococcus phages have been described from the coastal seas of the Northwest Pacific Ocean. Here a new Synechococcus phage, S-B43 was isolated from the Bohai Sea, a semi-closed coastal sea of the Northwest Pacific Ocean. S-B43 is a member of Myoviridae, containing 275 predicted open reading frames. Fourteen auxiliary metabolic genes (AMG) were identified from the genome of S-B43, including five photosynthetic associated genes and several AMGs related to its adaption to the high turbidity and eutrophic coastal environment with a low ratio of phosphorus to nitrogen (HNLP). The occurrences of 31 AMGs among 34 cyanophage genomes indicates that AMGs zwf, gnd, speD, petF and those coding for FECH and thioredoxin were more common in coastal areas than in the open ocean and AMGs pebS and ho1 were more prevalent in the open ocean. The occurrence of cyanophage AMGs in different environments might be a reflection of the environmental adaption of their hosts. This study contributes to our understanding of the interactions between cyanobacteria and cyanophages and their environmental adaption to the coastal environment.

Comparative analysis of the total and active bacterial communities in the surface sediment of Lake Taihu

Bacterial communities play crucial roles in the biogeochemical cycle of the surface sediments of freshwater lakes, but previous studies on bacterial community changes in this habitat have mostly been based on the total bacterial community (DNA level), while an exploration of the active microbiota at the RNA level has been lacking. Herein, we analysed the bacterial communities in the surface sediments of Lake Taihu at the DNA and RNA levels. Using MiSeq sequencing and real-time quantification, we found that the sequencing and quantitative results obtained at the RNA level compared with the DNA level were more accurate in responding to the spatiotemporal dynamic changes of the bacterial community. Although both sequencing methods indicated that Proteobacteria, Chloroflexi, Acidobacteria, Nitrospirae, Bacteroidetes and Actinobacteria were the dominant phyla, the co-occurrence network at the RNA level could better reflect the close relationship between microorganisms in the surface sediment. Additionally, further analysis showed that Prochlorococcus and Microcystis were the most relevant and dominant genera of Cyanobacteria in the total and active bacterial communities, respectively; our results also demonstrated that the analysis of Cyanobacteria-related groups at the RNA level was more 'informative'.

Diverse, Abundant, and Novel Viruses Infecting the Marine Roseobacter RCA Lineage

Many major marine bacterial lineages such as SAR11, Prochlorococcus, SAR116, and several Roseobacter lineages have members that are abundant, relatively slow-growing, and genome streamlined. The isolation of phages that infect SAR11 and SAR116 have demonstrated the dominance of these phages in the marine virosphere. However, no phages have been isolated from bacteria in the Roseobacter RCA lineage, another abundant group of marine bacteria. In this study, seven RCA phages that infect three different RCA strains were isolated and characterized. All seven RCA phages belong to the Podoviridae family and have genome sizes ranging from 39.6 to 58.1 kb. Interestingly, three RCA phages (CRP-1, CRP-2, and CRP-3) show similar genomic content and architecture as SAR116 phage HMO-2011, which represents one of the most abundant known viral groups in the ocean. The high degree of homology among CRP-1, CRP-2, CRP-3, and HMO-2011 resulted in the contribution of RCA phages to the dominance of the HMO-2011-type group. CRP-4 and CRP-5 are similar to the Cobavirus group roseophages in terms of gene content and organization. The remaining two RCA phages, CRP-6 and CRP-7, show limited genomic similarity with known phages and represent two new phage groups. Metagenomic fragment recruitment analyses reveal that these RCA phage groups are much more abundant in the ocean than most existing marine roseophage groups. The characterization of these RCA phages has greatly expanded our understanding of the genomic diversity and evolution of marine roseophages and suggests the critical need for isolating phages from the abundant but "unculturable" bacteria.IMPORTANCE The RCA lineage of the marine Roseobacter group represents one of the slow-growing but dominant components of marine microbial communities. Although dozens of roseophages have been characterized, no phages infecting RCA strains have been reported. In this study, we reported on the first RCA phage genomes and investigated their distribution pattern and relative abundance in comparison with other important marine phage groups. Two of the four RCA phage groups were found closely related to previously reported SAR116 phage HMO-2011 and Cobavirus group roseophages, respectively. The remaining two groups are novel in the genome contents. Our study also revealed that RCA phages are widely distributed and exhibit high abundance in marine viromic data sets. Altogether, our findings have greatly broadened our understanding of RCA phages and emphasize the ecological and evolutionary importance of RCA phages in the marine virosphere.

Prevalence of viral photosynthesis genes along a freshwater to saltwater transect in Southeast USA

Bacteriophages encode host-acquired functional genes known as auxiliary metabolic genes (AMGs). Photosynthesis AMGs are commonly found in marine cyanobacteria-infecting Myoviridae and Podoviridae cyanophages, but their ecology remains understudied in freshwater environments. To advance knowledge of this issue, we analysed viral metagenomes collected in the summertime for four years from five lakes and two estuarine locations interconnected by the Chattahoochee River, Southeast USA. Sequences representing ten different AMGs were recovered and found to be prevalent in all sites. Most freshwater AMGs were 10-fold less abundant than estuarine and marine AMGs and were encoded by novel Myoviridae and Podoviridae cyanophage genera. Notably, several of the corresponding viral genomes showed endemism to a specific province along the river. This translated into psbA gene phylogenetic clustering patterns that matched a marine vs. freshwater origin indicating that psbA may serve as a robust classification and source-tracking biomarker. Genomes classified in a novel viral lineage represented by isolate S-EIVl contained psbA, which is unprecedented for this lineage. Collectively, our findings indicated that the acquisition of photosynthesis AMGs is a widespread strategy used by cyanophages in aquatic ecosystems, and further indicated the existence of viral provinces in which certain viral species and/or genotypes are locally abundant.

Reannotation of the Ribonucleotide Reductase in a Cyanophage Reveals Life History Strategies Within the Virioplankton

Ribonucleotide reductases (RNRs) are ancient enzymes that catalyze the reduction of ribonucleotides to deoxyribonucleotides. They are required for virtually all cellular life and are prominent within viral genomes. RNRs share a common ancestor and must generate a protein radical for direct ribonucleotide reduction. The mechanisms by which RNRs produce radicals are diverse and divide RNRs into three major classes and several subclasses. The diversity of radical generation methods means that cellular organisms and viruses typically contain the RNR best-suited to the environmental conditions surrounding DNA replication. However, such diversity has also fostered high rates of RNR misannotation within subject sequence databases. These misannotations have resulted in incorrect translative presumptions of RNR biochemistry and have diminished the utility of this marker gene for ecological studies of viruses. We discovered a misannotation of the RNR gene within the Prochlorococcus phage P-SSP7 genome, which caused a chain of misannotations within commonly observed RNR genes from marine virioplankton communities. These RNRs are found in marine cyanopodo- and cyanosiphoviruses and are currently misannotated as Class II RNRs, which are O2-independent and require cofactor B12. In fact, these cyanoviral RNRs are Class I enzymes that are O2-dependent and may require a di-metal cofactor made of Fe, Mn, or a combination of the two metals. The discovery of an overlooked Class I β subunit in the P-SSP7 genome, together with phylogenetic analysis of the α and β subunits confirms that the RNR from P-SSP7 is a Class I RNR. Phylogenetic and conserved residue analyses also suggest that the P-SSP7 RNR may constitute a novel Class I subclass. The reannotation of the RNR clade represented by P-SSP7 means that most lytic cyanophage contain Class I RNRs, while their hosts, B12-producing Synechococcus and Prochlorococcus, contain Class II RNRs. By using a Class I RNR, cyanophage avoid a dependence on host-produced B12, a more effective strategy for a lytic virus. The discovery of a novel RNR β subunit within cyanopodoviruses also implies that some unknown viral genes may be familiar cellular genes that are too divergent for homology-based annotation methods to identify.

Cylindrospermopsis raciborskii Virus and host: genomic characterization and ecological relevance

Cylindrospermopsis (Raphidiopsis) raciborskii is an invasive, filamentous, nitrogen-fixing cyanobacterium that forms frequent blooms in freshwater habitats. While viruses play key roles in regulating the abundance, production and diversity of their hosts in aquatic ecosystems, the role(s) of viruses in the ecology of C. raciborskii is almost unexplored. Progress in this field has been hindered by the absence of a characterized virus-host system in C. raciborskii. To bridge this gap, we sequenced the genome of CrV-01T, a previously isolated cyanosiphovirus, and its host, C. raciborskii strain Cr2010. Analyses suggest that CrV-01T represents a distinct clade of siphoviruses infecting, and perhaps lysogenizing, filamentous cyanobacteria. Its genome contains unique features that include an intact CRISPR array and a 12 kb inverted duplication. Evidence suggests CrV-01T recently gained the ability to infect Cr2010 and recently lost the ability to form lysogens. The cyanobacterial host contains a CRISPR-Cas system with CRISPR spacers matching protospacers within the inverted duplication of the CrV-01T genome. Examination of metagenomes demonstrates that viruses with high genetic identity to CrV-01T, but lacking the inverted duplication, are present in C. raciborskii blooms in Australia. The unique genomic features of the CrV/Cr2010 system offers opportunities to investigate in more detail virus-host interactions in an ecologically important bloom-forming cyanobacterium.

A New Freshwater Cyanosiphovirus Harboring Integrase

Pelagic cyanobacteria are key players in the functioning of aquatic ecosystems, and their viruses (cyanophages) potentially affect the abundance and composition of cyanobacterial communities. Yet, there are few well-described freshwater cyanophages relative to their marine counterparts, and in general, few cyanosiphoviruses (family Siphoviridae) have been characterized, limiting our understanding of the biology and the ecology of this prominent group of viruses. Here, we characterize S-LBS1, a freshwater siphovirus lytic to a phycoerythrin-rich Synechococcus isolate (Strain TCC793). S-LBS1 has a narrow host range, a burst size of ∼400 and a relatively long infecting step before cell lysis occurs. It has a dsDNA 34,641 bp genome with putative genes for structure, DNA packing, lysis, replication, host interactions, DNA repair and metabolism. S-LBS1 is similar in genome size, genome architecture, and gene content, to previously described marine siphoviruses also infecting PE-rich Synechococcus, e.g., S-CBS1 and S-CBS3. However, unlike other Synechococcus phages, S-LBS1 encodes an integrase, suggesting its ability to establish lysogenic relationships with its host. Sequence recruitment from viral metagenomic data showed that S-LBS1-like viruses are diversely present in a wide range of aquatic environments, emphasizing their potential importance in controlling and structuring Synechococcus populations. A comparative analysis with 16 available sequenced cyanosiphoviruses reveals the absence of core genes within the genomes, suggesting high degree of genetic variability in siphoviruses infecting cyanobacteria. It is likely that cyanosiphoviruses have evolved as distinct evolutionary lineages and that adaptive co-evolution occurred between these viruses and their hosts (i.e., Synechococcus, Prochlorococcus, Nodularia, and Acaryochloris), constituting an important driving force for such phage diversification.

Genetic Diversity and Cooccurrence Patterns of Marine Cyanopodoviruses and Picocyanobacteria

Picocyanobacteria Prochlorococcus and Synechococcus are abundant in the global oceans and subject to active viral infection. In this study, the genetic diversity of picocyanobacteria and the genetic diversity of cyanopodoviruses were synchronously investigated along water columns in the equatorial Indian Ocean and over a seasonal time course in the coastal Sanya Bay, South China Sea. Using the 16S-23S rRNA internal transcribed spacer (ITS)-based clone library and quantitative PCR (qPCR) analyses, the picocyanobacterial community composition and abundance were determined. Sanya Bay was dominated by clade II Synechococcus during all the seasons, and a typical population shift from high-light-adapted Prochlorococcus to low-light-adapted Prochlorococcus was found along the vertical profiles. Strikingly, the DNA polymerase gene sequences of cyanopodoviruses revealed a much greater genetic diversity than we expected. Nearly one-third of the phylogenetic groups were newly described here. No apparent seasonal pattern was observed for the Sanya Bay picocyanobacterial or cyanopodoviral communities. Different dominant cyanopodovirus lineages were identified for the coastal area, upper euphotic zone, and middle-to-lower euphotic zone of the open ocean. Diversity indices of both picocyanobacteria and cyanopodoviruses were highest in the middle euphotic zone and both were lower in the upper euphotic zone, reflecting a host-virus interaction. Cyanopodoviral communities differed significantly between the upper euphotic zone and the middle-to-lower euphotic zone, showing a vertical pattern similar to that of picocyanobacteria. However, in the surface waters of the open ocean, cyanopodoviruses exhibited no apparent biogeographic pattern, differing from picocyanobacteria. This study demonstrates correlated distribution patterns of picocyanobacteria and cyanopodoviruses, as well as the complex biogeography of cyanopodoviruses.IMPORTANCE Picocyanobacteria are highly diverse and abundant in the ocean and display remarkable global biogeography and a vertical distribution pattern. However, how the diversity and distribution of picocyanobacteria affect those of the viruses that infect them remains largely unknown. Here we synchronously analyzed the community structures of cyanopodoviruses and picocyanobacteria at spatial and temporal scales. Both spatial and temporal variations of cyanopodoviral communities can be linked to those of picocyanobacteria. The coastal area, upper euphotic zone, and middle-to-lower euphotic zone of the open ocean have distinct cyanopodoviral communities, showing horizontal and vertical variation patterns closely related to those of picocyanobacteria. These findings emphasize the driving force of host community in shaping the biogeographic structure of viruses. Our work provides important information for future assessments of the ecological roles of viruses and hosts for each other.

Co-existence of freshwater and marine T4-like myoviruses in a typical subtropical estuary

Viruses are the most abundant biological entities on Earth and play an important role in microbial community dynamics and biogeochemical cycling, yet their ecological characteristics in estuarine ecosystems are unclear. Here, virioplankton communities in a typical subtropical estuary, the Jiulong River estuary (JRE) in China, were investigated. The abundance of virioplankton ranged from 1.01 ± 0.05 × 107 to 1.62 ± 0.09 × 107 particles mL-1 in JRE, and the population size of viruses was correlated with temperature and nutrient levels. Three tailed viral morphotypes (myovirus, siphovirus and podovirus) were observed. Phylogenetic analysis showed that most of the g23 sequences in the JRE fell into three previously established groups (Marine, Paddy and Lake Groups) and two potential Estuary Groups. This demonstrates the co-existence of typical freshwater and marine T4-like myoviruses in the estuarine ecosystem, suggesting the movement of viruses and their hosts among biomes. Additionally, the spatial variation of g23 sequences suggests a geographic distribution pattern of T4-like myoviruses in the JRE, which might be shaped by the environmental gradient and/or their host distribution. These results provide valuable insights into the abundance, diversity and distribution patterns of virioplankton, as well as the factors influencing them, in subtropical estuarine ecosystems.

Pan-Genome Storage and Analysis Techniques

Computational pan-genome analysis has emerged from the rapid increase of available genome sequencing data. Starting from a microbial pan-genome, the concept has spread to a variety of species, such as plants or viruses. Characterizing a pan-genome provides insights into intra-species evolution, functions, and diversity. However, researchers face challenges such as processing and maintaining large datasets while providing accurate and efficient analysis approaches. Comparative genomics methods are required for detecting conserved and unique regions between a set of genomes. This chapter gives an overview of tools available for indexing pan-genomes, identifying the sub-regions of a pan-genome and offering a variety of downstream analysis methods. These tools are categorized into two groups, gene-based and sequence-based, according to the pan-genome identification method. We highlight the differences, advantages, and disadvantages between the tools, and provide information about the general workflow, methodology of pan-genome identification, covered functionalities, usability and availability of the tools.

Metabolic Genes within Cyanophage Genomes: Implications for Diversity and Evolution

Cyanophages, a group of viruses specifically infecting cyanobacteria, are genetically diverse and extensively abundant in water environments. As a result of selective pressure, cyanophages often acquire a range of metabolic genes from host genomes. The host-derived genes make a significant contribution to the ecological success of cyanophages. In this review, we summarize the host-derived metabolic genes, as well as their origin and roles in cyanophage evolution and important host metabolic pathways, such as the light-dependent reactions of photosynthesis, the pentose phosphate pathway, nutrient acquisition and nucleotide biosynthesis. We also discuss the suitability of the host-derived metabolic genes as potential diagnostic markers for the detection of genetic diversity of cyanophages in natural environments.

Characterization and Complete Genome Sequences of Three N4-Like Roseobacter Phages Isolated from the South China Sea

Three bacteriophages (RD-1410W1-01, RD-1410Ws-07, and DS-1410Ws-06) were isolated from the surface water of Sanya Bay, northern South China Sea, on two marine bacteria type strains of the Roseobacter lineage. These phages have an isometric head and a short tail, morphologically belonging to the Podoviridae family. Two of these phages can infect four of seven marine roseobacter strains tested and the other one can infect three of them, showing relatively broader host ranges compared to known N4-like roseophages. One-step growth curves showed that these phages have similar short latent periods (1-2 h) but highly variable burst sizes (27-341 pfu cell(-1)). Their complete genomes show high level of similarities to known N4-like roseophages in terms of genome size, G + C content, gene content, and arrangement. The morphological and genomic features of these phages indicate that they belong to the N4likevirus genus. Moreover, comparative genomic analysis based on 43 N4-like phages (10 roseobacter phages and 33 phages infecting other lineages of bacteria) revealed a core genome of 18 genes shared by all the 43 phages and 38 genes shared by all the ten roseophages. The 38 core genes of N4-like roseophages nearly make up 70 % of each genome in length. Phylogenetic analysis based on the concatenated core gene products showed that our phage isolates represent two new phyletic branches, suggesting the broad genetic diversity of marine N4-like roseophages remains.