ViPTree: the viral proteomic tree server

Widespread distribution of prophages signaling the potential for adaptability and pathogenicity evolution of Ralstonia solanacearum species complex

Integrated bacteriophages (prophages) can impact host cells, affecting their lifestyle, genomic diversity, and fitness. However, many basic aspects of how these organisms affect the host cell remain poorly understood. Ralstonia solanacearum is a gram-negative plant pathogenic bacterium that encompasses a great diversity of ecotypes regarded as a species complex (R. solanacearum Species Complex - RSSC). RSSC genomes have a mosaic structure containing numerous elements, signaling the potential for its evolution through horizontal gene transfer. Here, we analyzed 120 Ralstonia spp. genomes from the public database to identify prophage sequences. In total, 379 prophage-like elements were found in the chromosome and megaplasmid of Ralstonia spp. These elements encode genes related to host fitness, virulence factors, antibiotic resistance, and niche adaptation, which might contribute to RSSC adaptability. Prophage-like elements are widespread into the complex in different species and geographic origins, suggesting that the RSSC phages are ancestrally acquired. Complete prophages belonging to the families Inoviridae, Myoviridae, and Siphoviridae were found, being the members of Inoviridae the most abundant. Analysis of CRISPR-Cas spacer sequences demonstrated the presence of prophages sequences that indicate successive infection events during bacterial evolution. Besides complete prophages, we also demonstrated 14 novel putative prophages integrated into Ralstonia spp. genomes. Altogether, our results provide insights into the diversity of prophages in RSSC genomes and suggest that these elements may deeply affect the virulence and host adaptation and shaping the genomes among the strains of this important pathogen.

Induction and Genome Analysis of HY01, a Newly Reported Prophage from an Emerging Shrimp Pathogen Vibrio campbellii

Vibrio campbellii is an emerging aquaculture pathogen that causes luminous vibriosis in farmed shrimp. Although prophages in various aquaculture pathogens have been widely reported, there is still limited knowledge regarding prophages in the genome of pathogenic V. campbellii. Here, we describe the full-genome sequence of a prophage named HY01, induced from the emerging shrimp pathogen V. campbellii HY01. The phage HY01 was induced by mitomycin C and was morphologically characterized as long tailed phage. V. campbellii phage HY01 is composed of 41,772 bp of dsDNA with a G+C content of 47.45%. A total of 60 open reading frames (ORFs) were identified, of which 31 could be predicted for their biological functions. Twenty seven out of 31 predicted protein coding regions were matched with several encoded proteins of various Enterobacteriaceae, Pseudomonadaceae, Vibrionaceae, and other phages of Gram-negative bacteria. Interestingly, the comparative genome analysis revealed that the phage HY01 was only distantly related to Vibrio phage Va_PF430-3_p42 of fish pathogen V. anguillarum but differed in genomic size and gene organization. The phylogenetic tree placed the phage together with Siphoviridae family. Additionally, a survey of Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) spacers revealed two matching sequences between phage HY01 genome and viral spacer sequence of Vibrio spp. The spacer results combined with the synteny results suggest that the evolution of V. campbellii phage HY01 is driven by the horizontal genetic exchange between bacterial families belonging to the class of Gammaproteobacteria.

Characterization and genome analysis of phage AL infecting Pseudoalteromonas marina

Although Pseudoalteromonas is an abundant, ubiquitous, marine algae-associated bacterial genus, there is still little information on their phages. In the present study, a marine phage AL, infecting Pseudoalteromonas marina, was isolated from the coastal waters off Qingdao. The AL phage is a siphovirus with an icosahedral head of 53 ± 1 nm and a non-contractile tail, length of 99 ± 1 nm. A one-step growth curve showed that the latent period was approximately 70 min, the rise period was 50 min, and the burst size was 227 pfu/cell. The genome sequence of this phage is a 33,582 bp double-stranded DNA molecule with a GC content of 40.1 %, encoding 52 open reading frames (ORFs). The order of the functional genes, especially those related to the structure module, is highly conserved and basically follows the common pattern used by siphovirus. The stable order has been formed during the long-term evolution of phages in the siphovirus group, which has helped the phages to maintain their normal morphology and function. Phylogenetic trees based on the major capsid protein (mcp) and genome-wide sequence have shown that the AL phage is closely related to four Pseudoalteromonas phages, including PHS21, PHS3, SL25 and Pq0. Further analysis using all-to-all BLASTP also confirmed that this phage shared high sequence homology with the same four Pseudoalteromonas phages, with amino acid sequence identities ranging from 44 % to 71 %. In particular, their similarity in virion structure module may imply that these phages share common assembly mechanism characteristics and infection pathways. Pseudoalteromonas phage AL not only provides basic information for the further study of the evolution of Pseudoalteromonas phages and interactions between marine phage and host but also helps to explain the unknown viral sequences in the metagenomic databases.

Long-read metagenomics using PromethION uncovers oral bacteriophages and their interaction with host bacteria

Bacteriophages (phages), or bacterial viruses, are very diverse and highly abundant worldwide, including as a part of the human microbiomes. Although a few metagenomic studies have focused on oral phages, they relied on short-read sequencing. Here, we conduct a long-read metagenomic study of human saliva using PromethION. Our analyses, which integrate both PromethION and HiSeq data of >30 Gb per sample with low human DNA contamination, identify hundreds of viral contigs; 0-43.8% and 12.5-56.3% of the confidently predicted phages and prophages, respectively, do not cluster with those reported previously. Our analyses demonstrate enhanced scaffolding, and the ability to place a prophage in its host genomic context and enable its taxonomic classification. Our analyses also identify a Streptococcus phage/prophage group and nine jumbo phages/prophages. 86% of the phage/prophage group and 67% of the jumbo phages/prophages contain remote homologs of antimicrobial resistance genes. Pan-genome analysis of the phages/prophages reveals remarkable diversity, identifying 0.3% and 86.4% of the genes as core and singletons, respectively. Furthermore, our study suggests that oral phages present in human saliva are under selective pressure to escape CRISPR immunity. Our study demonstrates the power of long-read metagenomics utilizing PromethION in uncovering bacteriophages and their interaction with host bacteria.

Bacterial membrane vesicles and phages in blood after consumption of lacticaseibacillus rhamnosus JB-1

Gut microbiota have myriad roles in host physiology, development, and immunity. Though confined to the intestinal lumen by the epithelia, microbes influence distal systems via poorly characterized mechanisms. Recent work has considered the role of extracellular vesicles in interspecies communication, but whether they are involved in systemic microbe-host interaction is unclear. Here, we show that distinctive nanoparticles can be isolated from mouse blood within 2.5 h of consuming Lacticaseibacillus rhamnosus JB-1. In contrast to blood nanoparticles from saline-fed mice, they reproduced lipoteichoic acid-mediated immune functions of the original bacteria, including activation of TLR2 and increased IL-10 expression by dendritic cells. Like the fed bacteria, they also reduced IL-8 induced by TNF in an intestinal epithelial cell line. Though enriched for host neuronal proteins, these isolated nanoparticles also contained proteins and viral (phage) DNA of fed bacterial origin. Our data strongly suggest that oral consumption of live bacteria rapidly leads to circulation of their membrane vesicles and phages and demonstrate a nanoparticulate pathway whereby beneficial bacteria and probiotics may systemically affect their hosts.

Biological and Genomic Characterization of a Novel Jumbo Bacteriophage, vB_VhaM_pir03 with Broad Host Lytic Activity against Vibrio harveyi

Vibrio harveyi is a Gram-negative marine bacterium that causes major disease outbreaks and economic losses in aquaculture. Phage therapy has been considered as a potential alternative to antibiotics however, candidate bacteriophages require comprehensive characterization for a safe and practical phage therapy. In this work, a lytic novel jumbo bacteriophage, vB_VhaM_pir03 belonging to the Myoviridae family was isolated and characterized against V. harveyi type strain DSM19623. It had broad host lytic activity against 31 antibiotic-resistant strains of V. harveyi, V. alginolyticus, V. campbellii and V. owensii. Adsorption time of vB_VhaM_pir03 was determined at 6 min while the latent-phase was at 40 min and burst-size at 75 pfu/mL. vB_VhaM_pir03 was able to lyse several host strains at multiplicity-of-infections (MOI) 0.1 to 10. The genome of vB_VhaM_pir03 consists of 286,284 base pairs with 334 predicted open reading frames (ORFs). No virulence, antibiotic resistance, integrase encoding genes and transducing potential were detected. Phylogenetic and phylogenomic analysis showed that vB_VhaM_pir03 is a novel bacteriophage displaying the highest similarity to another jumbo phage, vB_BONAISHI infecting Vibrio coralliilyticus. Experimental phage therapy trial using brine shrimp, Artemia salina infected with V. harveyi demonstrated that vB_VhaM_pir03 was able to significantly reduce mortality 24 h post infection when administered at MOI 0.1 which suggests that it can be an excellent candidate for phage therapy.

Diversity and Distribution of a Novel Genus of Hyperthermophilic Aquificae Viruses Encoding a Proof-Reading Family-A DNA Polymerase

Despite the high abundance of Aquificae in many geothermal systems, these bacteria are difficult to culture and no viruses infecting members of this phylum have been isolated. Here, we describe the complete, circular dsDNA Uncultivated Virus Genome (UViG) of Thermocrinis Octopus Spring virus (TOSV), derived from metagenomic data, along with eight related UViGs representing three additional viral species. Despite low overall similarity among viruses from different hot springs, the genomes shared a high degree of synteny, and encoded numerous genes for nucleotide metabolism, including a PolA-type DNA polymerase polyprotein with likely accessory functions, a DNA Pol III sliding clamp, a thymidylate kinase, a DNA gyrase, a helicase, and a DNA methylase. Also present were conserved genes predicted to code for phage capsid, large and small subunits of terminase, portal protein, holin, and lytic transglycosylase, all consistent with a distant relatedness to cultivated Caudovirales. These viruses are predicted to infect Aquificae, as multiple CRISPR spacers matching the viral genomes were identified within the genomes and metagenomic contigs from these bacteria. Based on the predicted atypical bi-directional replication strategy, low sequence similarity to known viral genomes, and unique position in gene-sharing networks, we propose a new putative genus, "Pyrovirus," in the order Caudovirales.

Dissolved Microcystin Release Coincident with Lysis of a Bloom Dominated by Microcystis spp. in Western Lake Erie Attributed to a Novel Cyanophage

Western Lake Erie (Laurentian Great Lakes) is prone to annual cyanobacterial harmful algal blooms (cHABs) dominated by Microcystis spp. that often yield microcystin toxin concentrations exceeding the federal EPA recreational contact advisory of 8 μg liter-1 In August 2014, microcystin levels were detected in finished drinking water above the World Health Organization 1.0 μg liter-1 threshold for consumption, leading to a 2-day disruption in the supply of drinking water for >400,000 residents of Toledo, Ohio (USA). Subsequent metatranscriptomic analysis of the 2014 bloom event provided evidence that release of toxin into the water supply was likely caused by cyanophage lysis that transformed a portion of the intracellular microcystin pool into the dissolved fraction, rendering it more difficult to eliminate during treatment. In August 2019, a similar increase in dissolved microcystins at the Toledo water intake was coincident with a viral lytic event caused by a phage consortium different in composition from what was detected following the 2014 Toledo water crisis. The most abundant viral sequence in metagenomic data sets was a scaffold from a putative member of the Siphoviridae, distinct from the Ma-LMM01-like Myoviridae that are typically documented to occur in western Lake Erie. This study provides further evidence that viral activity in western Lake Erie plays a significant role in transformation of microcystins from the particulate to the dissolved fraction and therefore requires monitoring efforts from local water treatment plants. Additionally, identification of multiple lytic cyanophages will enable the development of a quantitative PCR toolbox to assess viral activity during cHABs.IMPORTANCE Viral attack on cHABs may contribute to changes in community composition during blooms, as well as bloom decline, yet loss of bloom biomass does not eliminate the threat of cHAB toxicity. Rather, it may increase risks to the public by delivering a pool of dissolved toxin directly into water treatment utilities when the dominating Microcystis spp. are capable of producing microcystins. Detecting, characterizing, and quantifying the major cyanophages involved in lytic events will assist water treatment plant operators in making rapid decisions regarding the pool of microcystins entering the plant and the corresponding best practices to neutralize the toxin.

Seasonal Regime Shift in the Viral Communities of a Permafrost Thaw Lake

Permafrost thaw lakes including thermokarst lakes and ponds are ubiquitous features of Subarctic and Arctic landscapes and are hotspots of microbial activity. Input of terrestrial organic matter into the planktonic microbial loop of these lakes may greatly amplify global greenhouse gas emissions. This microbial loop, dominated in the summer by aerobic microorganisms including phototrophs, is radically different in the winter, when metabolic processes shift to the anaerobic degradation of organic matter. Little is known about the viruses that infect these microbes, despite evidence that viruses can control microbial populations and influence biogeochemical cycling in other systems. Here, we present the results of a metagenomics-based study of viruses in the larger than 0.22 µm fraction across two seasons (summer and winter) in a permafrost thaw lake in Subarctic Canada. We uncovered 351 viral populations (vOTUs) in the surface waters of this lake, with diversity significantly greater during the summer. We also identified and characterized several phage genomes and prophages, which were mostly present in the summer. Finally, we compared the viral community of this waterbody to other habitats and found unexpected similarities with distant bog lakes in North America.

Characterization of Novel Lytic Bacteriophages of Achromobacter marplantensis Isolated from a Pneumonia Patient

Achromobacter spp. are becoming increasingly associated with lung infections in patients suffering from cystic fibrosis (CF). A. marplatensis, which is closely related to A. xylosoxidans, has been isolated from the lungs of CF patients and other human infections. This article describes the isolation, morphology and characterization of two lytic bacteriophages specific for an A. marplatensis strain isolated from a pneumonia patient. This host strain was the causal agent of hospital acquired pneumonia–the first clinical report of such an occurrence. Full genome sequencing revealed bacteriophage genomes ranging in size from 45901 to 46,328 bp. Transmission electron microscopy revealed that the two bacteriophages AMA1 and AMA2 belonged to the Siphoviridae family. Host range analysis showed that their host range did not extend to A. xylosoxidans. The possibility exists for future testing of such bacteriophages in the control of Achromobacter infections such as those seen in CF and other infections of the lungs. The incidence of antibiotic resistance in this genus highlights the importance of seeking adjuncts and alternatives in CF and other lung infections.

Viruses of freshwater bloom-forming cyanobacteria: genomic features, infection strategies and coexistence with the host

Freshwater bloom-forming cyanobacteria densely grow in the aquatic environments, leading to an increase in the viral-contact rate. They possess numerous antiviral genes, as well as cell differentiation- and physiological performance-related genes, owing to genome expansion. Their genomic features and unique lifestyles suggest that they coexist with cyanoviruses in ways different from marine cyanobacteria. Furthermore, genome contents of isolated freshwater bloom-forming cyanobacterial viruses have little in common with those of marine cyanoviruses studied to date. They lack the marine cyanoviral hallmark genes that sustain photosynthetic activity and redirect host metabolism to viral reproduction; therefore, they are predicted to share metabolisms and precursor pools with host cyanobacteria to ensure efficient viral reproduction and avoid nutrient deficiencies and antiviral response. Additionally, cyanovirus-cyanobacteria coexistence strategies may change as bloom density increases. Diverse genotypic populations of cyanoviruses and hosts coexist and fluctuate under high viral-contact rate conditions, leading to their rapid coevolution through antiviral responses. The ancestral and newly evolved genotypes coexist, thereby expanding the diversity levels of host and viral populations. Bottleneck events occurring due to season-related decreases in bloom-forming species abundance provide each genotype within cyanobacterial population an equal chance to increase in prevalence during the next bloom and enhance further diversification.

From Orphan Phage to a Proposed New Family-the Diversity of N4-Like Viruses

Escherichia phage N4 was isolated in 1966 in Italy and has remained a genomic orphan for a long time. It encodes an extremely large virion-associated RNA polymerase unique for bacterial viruses that became characteristic for this group. In recent years, due to new and relatively inexpensive sequencing techniques the number of publicly available phage genome sequences expanded rapidly. This revealed new members of the N4-like phage group, from 33 members in 2015 to 115 N4-like viruses in 2020. Using new technologies and methods for classification, the Bacterial and Archaeal Viruses Subcommittee of the International Committee on Taxonomy of Viruses (ICTV) has moved the classification and taxonomy of bacterial viruses from mere morphological approaches to genomic and proteomic methods. The analysis of 115 N4-like genomes resulted in a huge reassessment of this group and the proposal of a new family "Schitoviridae", including eight subfamilies and numerous new genera.

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.

Isolation and Characterization of Shewanella Phage Thanatos Infecting and Lysing Shewanella oneidensis and Promoting Nascent Biofilm Formation

Species of the genus Shewanella are widespread in nature in various habitats, however, little is known about phages affecting Shewanella sp. Here, we report the isolation of phages from diverse freshwater environments that infect and lyse strains of Shewanella oneidensis and other Shewanella sp. Sequence analysis and microscopic imaging strongly indicate that these phages form a so far unclassified genus, now named Shewanella phage Thanatos, which can be positioned within the subfamily of Tevenvirinae (Duplodnaviria; Heunggongvirae; Uroviricota; Caudoviricetes; Caudovirales; Myoviridae; Tevenvirinae). We characterized one member of this group in more detail using S. oneidensis MR-1 as a host. Shewanella phage Thanatos-1 possesses a prolate icosahedral capsule of about 110 nm in height and 70 nm in width and a tail of about 95 nm in length. The dsDNA genome exhibits a GC content of about 34.5%, has a size of 160.6 kbp and encodes about 206 proteins (92 with an annotated putative function) and two tRNAs. Out of those 206, MS analyses identified about 155 phage proteins in PEG-precipitated samples of infected cells. Phage attachment likely requires the outer lipopolysaccharide of S. oneidensis, narrowing the phage's host range. Under the applied conditions, about 20 novel phage particles per cell were produced after a latent period of approximately 40 min, which are stable at a pH range from 4 to 12 and resist temperatures up to 55°C for at least 24 h. Addition of Thanatos to S. oneidensis results in partial dissolution of established biofilms, however, early exposure of planktonic cells to Thanatos significantly enhances biofilm formation. Taken together, we identified a novel genus of Myophages affecting S. oneidensis communities in different ways.

Large freshwater phages with the potential to augment aerobic methane oxidation

There is growing evidence that phages with unusually large genomes are common across various microbiomes, but little is known about their genetic inventories or potential ecosystem impacts. In the present study, we reconstructed large phage genomes from freshwater lakes known to contain bacteria that oxidize methane. Of manually curated genomes, 22 (18 are complete), ranging from 159 kilobase (kb) to 527 kb in length, were found to encode the pmoC gene, an enzymatically critical subunit of the particulate methane monooxygenase, the predominant methane oxidation catalyst in nature. The phage-associated PmoC sequences show high similarity to (>90%), and affiliate phylogenetically with, those of coexisting bacterial methanotrophs, including members of Methyloparacoccus, Methylocystis and Methylobacter spp. In addition, pmoC-phage abundance patterns correlate with those of the coexisting bacterial methanotrophs, supporting host-phage relationships. Future work is needed to determine whether phage-associated PmoC has similar functions to additional copies of PmoC encoded in bacterial genomes, thus contributing to growth on methane. Transcriptomics data from Lake Rotsee (Switzerland) showed that some phage-associated pmoC genes were highly expressed in situ and, of interest, that the most rapidly growing methanotroph was infected by three pmoC-phages. Thus, augmentation of bacterial methane oxidation by pmoC-phages during infection could modulate the efflux of this potent greenhouse gas into the environment.

Characterization of vB_StuS_MMDA13, a Newly Discovered Bacteriophage Infecting the Agar-Degrading Species Sphingomonas turrisvirgatae

Members of Sphingomonas genus have gained a notable interest for their use in a wide range of biotechnological applications, ranging from bioremediation to the production of valuable compounds of industrial interest. To date, knowledge on phages targeting Sphingomonas spp. are still scarce. Here, we describe and characterize a lytic bacteriophage, named vB_StuS_MMDA13, able to infect the Sphingomonas turrisvirgatae MCT13 type strain. Physiological characterization demonstrated that vB_StuS_MMDA13 has a narrow host range, a long latency period, a low burst size, and it is overall stable to both temperature and pH variations. The phage has a double-stranded DNA genome of 63,743 bp, with 89 open reading frames arranged in two opposite arms separated by a 1186 bp non-coding region and shows a very low global similarity to any other known phages. Interestingly, vB_StuS_MMDA13 is endowed with an original nucleotide modification biosynthetic gene cluster, which greatly differs from those of its most closely related phages of the Nipunavirus genus. vB_StuS_MMDA13 is the first characterized lytic bacteriophage of the Siphoviridae family infecting members of the Sphingomonas genus.

The Isolation and Genome Sequencing of Five Novel Bacteriophages From the Rumen Active Against Butyrivibrio fibrisolvens

Although the prokaryotic communities of the rumen microbiome are being uncovered through genome sequencing, little is known about the resident viral populations. Whilst temperate phages can be predicted as integrated prophages when analyzing bacterial and archaeal genomes, the genetics underpinning lytic phages remain poorly characterized. To the five genomes of bacteriophages isolated from rumen-associated samples sequenced and analyzed previously, this study adds a further five novel genomes and predictions gleaned from them to further the understanding of the rumen phage population. Lytic bacteriophages isolated from fresh ovine and bovine fecal and rumen fluid samples were active against the predominant fibrolytic ruminal bacterium Butyrivibrio fibrisolvens. The double stranded DNA genomes were sequenced and reconstructed into single circular complete contigs. Based on sequence similarity and genome distances, the five phages represent four species from three separate genera, consisting of: (1) Butyrivibrio phages Arian and Bo-Finn; (2) Butyrivibrio phages Idris and Arawn; and (3) Butyrivibrio phage Ceridwen. They were predicted to all belong to the Siphoviridae family, based on evidence in the genomes such as size, the presence of the tail morphogenesis module, genes that share similarity to those in other siphovirus isolates and phylogenetic analysis using phage proteomes. Yet, phylogenomic analysis and sequence similarity of the entire phage genomes revealed that these five phages are unique and novel. These phages have only been observed undergoing the lytic lifecycle, but there is evidence in the genomes of phages Arawn and Idris for the potential to be temperate. However, there is no evidence in the genome of the bacterial host Butyrivibrio fibrisolvens of prophage genes or genes that share similarity with the phage genomes.

Genome Analysis of Two Novel Synechococcus Phages That Lack Common Auxiliary Metabolic Genes: Possible Reasons and Ecological Insights by Comparative Analysis of Cyanomyoviruses

The abundant and widespread unicellular cyanobacteria Synechococcus plays an important role in contributing to global phytoplankton primary production. In the present study, two novel cyanomyoviruses, S-N03 and S-H34 that infected Synechococcus MW02, were isolated from the coastal waters of the Yellow Sea. S-N03 contained a 167,069-bp genome comprising double-stranded DNA with a G + C content of 50.1%, 247 potential open reading frames and 1 tRNA; S-H34 contained a 167,040-bp genome with a G + C content of 50.1%, 246 potential open reading frames and 5 tRNAs. These two cyanophages contain fewer auxiliary metabolic genes (AMGs) than other previously isolated cyanophages. S-H34 in particular, is currently the only known cyanomyovirus that does not contain any AMGs related to photosynthesis. The absence of such common AMGs in S-N03 and S-H34, their distinct evolutionary history and ecological features imply that the energy for phage production might be obtained from other sources rather than being strictly dependent on the maintenance of photochemical ATP under high light. Phylogenetic analysis showed that the two isolated cyanophages clustered together and had a close relationship with two other cyanophages of low AMG content. Comparative genomic analysis, habitats and hosts across 81 representative cyanomyovirus showed that cyanomyovirus with less AMGs content all belonged to Synechococcus phages isolated from eutrophic waters. The relatively small genome size and high G + C content may also relate to the lower AMG content, as suggested by the significant correlation between the number of AMGs and G + C%. Therefore, the lower content of AMG in S-N03 and S-H34 might be a result of viral evolution that was likely shaped by habitat, host, and their genomic context. The genomic content of AMGs in cyanophages may have adaptive significance and provide clues to their evolution.

Application of Bacteriophages to Control Vibrio alginolyticus Contamination in Oyster (Saccostrea glomerata) Larvae

Mortalities of bivalve larvae and spat linked with Vibrio spp. infection have been described in hatcheries since 1959, causing potential development of resistant bacteria. A reliable and sustainable solution to this problem is yet to be developed. Potential treatment of bacterial infection with bacteriophages is gaining interest in aquaculture as a more sustainable option for managing Vibrio spp. infection. This study assessed the effectiveness of bacteriophages (Φ-5, Φ-6, and Φ-7) against pathogenic Vibrio isolates (USC-26004 and USC-26005). These phage isolates were found to belong to the Myoviridae viral family. A total of 212 ORFs of Φ-5 were identified and annotated. The genome of this phage contained putative thymidine kinase and lysin enzyme. During infections with phages, the OD values of the isolates USC-26005 and USC-26004 remained stable at a much lower reading compared to the control after 9 h of incubation. Mortality rate of oyster (Saccostrea glomerata) larvae was 28.2 ± 3.5% in the bacteriophage treatment group, compared to 77.9 ± 9.1% in the bacterial treatment group after 24 h incubation. Findings of this study indicate that lytic phages might be utilized as potential bio-control agents of luminescent bacterial disease in oyster hatcheries.

Challenges of Studying the Human Virome - Relevant Emerging Technologies

In this review we provide an overview of current challenges and advances in bacteriophage research within the growing field of viromics. In particular, we discuss, from a human virome study perspective, the current and emerging technologies available, their limitations in terms of de novo discoveries, and possible solutions to overcome present experimental and computational biases associated with low abundance of viral DNA or RNA. We summarize recent breakthroughs in metagenomics assembling tools and single-cell analysis, which have the potential to increase our understanding of phage biology, diversity, and interactions with both the microbial community and the human body. We expect that these recent and future advances in the field of viromics will have a strong impact on how we develop phage-based therapeutic approaches.

Computational approaches in viral ecology

Dynamic virus-host interactions play a critical role in regulating microbial community structure and function. Yet for decades prior to the genomics era, viruses were largely overlooked in microbial ecology research, as only low-throughput culture-based methods of discovering viruses were available. With the advent of metagenomics, culture-independent techniques have provided exciting opportunities to discover and study new viruses. Here, we review recently developed computational methods for identifying viral sequences, exploring viral diversity in environmental samples, and predicting hosts from metagenomic sequence data. Methods to analyze viruses in silico utilize unconventional approaches to tackle challenges unique to viruses, such as vast diversity, mosaic viral genomes, and the lack of universal marker genes. As the field of viral ecology expands exponentially, computational advances have become increasingly important to gain insight into the role viruses in diverse habitats.

Extreme Viral Partitioning in a Marine-Derived High Arctic Lake

The Arctic is warming at an accelerating pace, and the rise in temperature has increasing impacts on the Arctic biome. Lakes are integrators of their surroundings and thus excellent sentinels of environmental change. Despite their importance in the regulation of key microbial processes, viruses remain largely uncharacterized in Arctic lacustrine environments. We sampled a highly stratified meromictic lake near the northern limit of the Canadian High Arctic, a region in rapid transition due to climate change. We found that the different layers of the lake harbored viral communities that were strikingly dissimilar and highly divergent from known viruses. Viruses were more abundant in the deepest part of the lake containing ancient Arctic Ocean seawater that was trapped during glacial retreat and were genomically unlike any viruses previously described. This research demonstrates the complexity and novelty of viral communities in an environment that is vulnerable to ongoing perturbation.

High-latitude, perennially stratified (meromictic) lakes are likely to be especially vulnerable to climate warming because of the importance of ice in maintaining their water column structure and associated distribution of microbial communities. This study aimed to characterize viral abundance, diversity, and distribution in a meromictic lake of marine origin on the far northern coast of Ellesmere Island, in the Canadian High Arctic. We collected triplicate samples for double-stranded DNA (dsDNA) viromics from five depths that encompassed the major features of the lake, as determined by limnological profiling of the water column. Viral abundance and virus-to-prokaryote ratios were highest at greater depths, while bacterial and cyanobacterial counts were greatest in the surface waters. The viral communities from each zone of the lake defined by salinity, temperature, and dissolved oxygen concentrations were markedly distinct, suggesting that there was little exchange of viral types among lake strata. Ten viral assembled genomes were obtained from our libraries, and these also segregated with depth. This well-defined structure of viral communities was consistent with that of potential hosts. Viruses from the monimolimnion, a deep layer of ancient Arctic Ocean seawater, were more diverse and relatively abundant, with few similarities to available viral sequences. The Lake A viral communities also differed from published records from the Arctic Ocean and meromictic Ace Lake in Antarctica. This first characterization of viral diversity from this sentinel environment underscores the microbial richness and complexity of an ecosystem type that is increasingly exposed to major perturbations in the fast-changing Arctic.

IMPORTANCE The Arctic is warming at an accelerating pace, and the rise in temperature has increasing impacts on the Arctic biome. Lakes are integrators of their surroundings and thus excellent sentinels of environmental change. Despite their importance in the regulation of key microbial processes, viruses remain largely uncharacterized in Arctic lacustrine environments. We sampled a highly stratified meromictic lake near the northern limit of the Canadian High Arctic, a region in rapid transition due to climate change. We found that the different layers of the lake harbored viral communities that were strikingly dissimilar and highly divergent from known viruses. Viruses were more abundant in the deepest part of the lake containing ancient Arctic Ocean seawater that was trapped during glacial retreat and were genomically unlike any viruses previously described. This research demonstrates the complexity and novelty of viral communities in an environment that is vulnerable to ongoing perturbation.

Tara Oceans: towards global ocean ecosystems biology

A planetary-scale understanding of the ocean ecosystem, particularly in light of climate change, is crucial. Here, we review the work of Tara Oceans, an international, multidisciplinary project to assess the complexity of ocean life across comprehensive taxonomic and spatial scales. Using a modified sailing boat, the team sampled plankton at 210 globally distributed sites at depths down to 1,000 m. We describe publicly available resources of molecular, morphological and environmental data, and discuss how an ecosystems biology approach has expanded our understanding of plankton diversity and ecology in the ocean as a planetary, interconnected ecosystem. These efforts illustrate how global-scale concepts and data can help to integrate biological complexity into models and serve as a baseline for assessing ecosystem changes and the future habitability of our planet in the Anthropocene epoch.

In silico Prediction of Virus-Host Interactions for Marine Bacteroidetes With the Use of Metagenome-Assembled Genomes

Bacteroidetes is one of the most abundant heterotrophic bacterial taxa in the ocean and play crucial roles in recycling phytoplankton-derived organic matter. Viruses of Bacteroidetes are also expected to have an important role in the regulation of host communities. However, knowledge on marine Bacteroidetes viruses is biased toward cultured viruses from a few species, mainly fish pathogens or Bacteroidetes not abundant in marine environments. In this study, we investigated the recently reported 1,811 marine viral genomes to identify putative Bacteroidetes viruses using various in silico host prediction techniques. Notably, we used microbial metagenome-assembled genomes (MAGs) to augment the marine Bacteroidetes reference genomic data. The examined viral genomes and MAGs were derived from simultaneously collected samples. Using nucleotide sequence similarity-based host prediction methods, we detected 31 putative Bacteroidetes viral genomes. The MAG-based method substantially enhanced the predictions (26 viruses) when compared with the method that is solely based on the reference genomes from NCBI RefSeq (7 viruses). Previously unrecognized genus-level groups of Bacteroidetes viruses were detected only by the MAG-based method. We also developed a host prediction method based on the proportion of Bacteroidetes homologs in viral genomes, which detected 321 putative Bacteroidetes virus genomes including 81 that were newly recognized as Bacteroidetes virus genomes. The majority of putative Bacteroidetes viruses were detected based on the proportion of Bacteroidetes homologs in both RefSeq and MAGs; however, some were detected in only one of the two datasets. Putative Bacteroidetes virus lineages included not only relatives of known viruses but also those phylogenetically distant from the cultured viruses, such as marine Far-T4 like viruses known to be widespread in aquatic environments. Our MAG and protein homology-based host prediction approaches enhanced the existing knowledge on the diversity of Bacteroidetes viruses and their potential interaction with their hosts in marine environments.

Taxonomy of prokaryotic viruses: 2018-2019 update from the ICTV Bacterial and Archaeal Viruses Subcommittee

This article is a summary of the activities of the ICTV's Bacterial and Archaeal Viruses Subcommittee for the years 2018 and 2019. Highlights include the creation of a new order, 10 families, 22 subfamilies, 424 genera and 964 species. Some of our concerns about the ICTV's ability to adjust to and incorporate new DNA- and protein-based taxonomic tools are discussed.

Isolation and complete genome sequence of a novel cyanophage, S-B05, infecting an estuarine Synechococcus strain: insights into environmental adaptation

A new cyanophage, S-B05, infecting a phycoerythrin-enriched (PE-type) Synechococcus strain was isolated by the liquid infection method, and its morphology and genetic features were examined. Phylogenetic analysis and morphological observation confirmed that S-B05 belongs to the family Myoviridae of the order Caudovirales. Its genome was fully sequenced, and found to be 208,857 bp in length with a G + C content of 39.9%. It contained 280 potential open reading frames and 123 conserved domains. Ninety-eight functional genes responsible for cyanophage structuring and packaging, DNA replication and regulation, and photosynthesis were identified, as well as genes encoding 172 hypothetical proteins. The genome of S-B05 is most similar to that of Prochlorococcus phage P-TIM68. Homologues of open reading frames of S-B05 can be found in various marine environments, as revealed by comparison of the S-B05 genome sequence to sequences in marine viral metagenomic databases. The presence of auxiliary metabolic genes (AMGs) related to photosynthesis, carbon metabolism, and phosphorus assimilation, as well as the phylogenetic relationships based on AMGs and the complete genome sequence, reflect the phage-host interaction mechanism or the specific adaptation strategy of the host to environmental conditions. The genome sequence information reported here will provide an important basis for further study of the adaptive evolution and ecological role of cyanophages and their hosts in the marine environment.

Genomic Analysis of Mic1 Reveals a Novel Freshwater Long-Tailed Cyanophage

Lake Chaohu, one of the five largest freshwater lakes in China, has been suffering from severe cyanobacterial blooms in the summer for many years. Cyanophages, the viruses that specifically infect cyanobacteria, play a key role in modulating cyanobacterial population, and thus regulate the emergence and decline of cyanobacterial blooms. Here we report a long-tailed cyanophage isolated from Lake Chaohu, termed Mic1, which specifically infects the cyanobacterium Microcystis aeruginosa. Mic1 has an icosahedral head of 88 nm in diameter and a long flexible tail of 400 nm. It possesses a circular genome of 92,627 bp, which contains 98 putative open reading frames. Genome sequence analysis enabled us to define a novel terminase large subunit that consists of two types of intein, indicating that the genome packaging of Mic1 is under fine control via posttranslational maturation of the terminase. Moreover, phylogenetic analysis suggested Mic1 and mitochondria share a common evolutionary origin of DNA polymerase γ gene. All together, these findings provided a start-point for investigating the co-evolution of cyanophages and its cyanobacterial hosts.

Genomic and protein structure modelling analysis depicts the origin and pathogenicity of 2019-nCoV, a new coronavirus which caused a pneumonia outbreak in Wuhan, China

Background: A pandemic outbreak caused by a novel coronavirus, 2019-nCoV, has originated from Wuhan, China and spread to many countries around the world. The outbreak has led to around 45 thousand cases and over one thousand death so far. Methods: Phylogenetic analysis and sequence alignment were used to align the whole genome sequence of 2019-nCoV with other over 200 sequences of coronaviruses to predict the origin of this novel virus. In addition, protein modeling and analysis were performed to access the potential binding of the spike protein of 2019-nCoV with human cell receptor, angiotensin-converting enzyme 2 (ACE2). Results: Detailed genomic and structure-based analysis of a new coronavirus, namely 2019-nCoV, showed that the new virus is a new type of bat coronavirus and is genetically fairly distant from the human SARS coronavirus. Structure analysis of the spike (S) protein of this new virus showed that its S protein only binds much weaker to the ACE2 receptor on human cells whereas the human SARS coronavirus exhibits strongly affinity to the ACE receptor. Conclusions: These findings suggest that the new virus should theoretically not be able to cause very serious human infection when compared to human SARS virus. However, the lower pathogenicity of this new virus may lead to longer incubation time and better adaption to human, which may favor its efficient transmission in human. These data are important to guide design of infection control policy and inform the public on the nature of threat imposed by 2019-nCov. Most importantly, using the analysis platform that we have developed, we should be able to predict whether the new mutations could lead to the increase of infectivity of the mutated virus in a very short time.

Iroki: automatic customization and visualization of phylogenetic trees

Phylogenetic trees are an important analytical tool for evaluating community diversity and evolutionary history. In the case of microorganisms, the decreasing cost of sequencing has enabled researchers to generate ever-larger sequence datasets, which in turn have begun to fill gaps in the evolutionary history of microbial groups. However, phylogenetic analyses of these types of datasets create complex trees that can be challenging to interpret. Scientific inferences made by visual inspection of phylogenetic trees can be simplified and enhanced by customizing various parts of the tree. Yet, manual customization is time-consuming and error prone, and programs designed to assist in batch tree customization often require programming experience or complicated file formats for annotation. Iroki, a user-friendly web interface for tree visualization, addresses these issues by providing automatic customization of large trees based on metadata contained in tab-separated text files. Iroki’s utility for exploring biological and ecological trends in sequencing data was demonstrated through a variety of microbial ecology applications in which trees with hundreds to thousands of leaf nodes were customized according to extensive collections of metadata. The Iroki web application and documentation are available at https://www.iroki.net or through the VIROME portal http://virome.dbi.udel.edu. Iroki’s source code is released under the MIT license and is available at https://github.com/mooreryan/iroki.

Bacteriophages of Klebsiella spp., their diversity and potential therapeutic uses

Klebsiella spp. are commensals of the human microbiota, and a leading cause of opportunistic nosocomial infections. The incidence of multidrug resistant (MDR) strains of Klebsiella pneumoniae causing serious infections is increasing, and Klebsiella oxytoca is an emerging pathogen. Alternative strategies to tackle infections caused by these bacteria are required as strains become resistant to last-resort antibiotics such as colistin. Bacteriophages (phages) are viruses that can infect and kill bacteria. They and their gene products are now being considered as alternatives or adjuncts to antimicrobial therapies. Several in vitro and in vivo studies have shown the potential for lytic phages to combat MDR K. pneumoniae infections. Ready access to cheap sequencing technologies has led to a large increase in the number of genomes available for Klebsiella-infecting phages, with these phages being heterogeneous at the whole-genome level. This review summarizes our current knowledge on phages of Klebsiella spp. and highlights technological and biological issues relevant to the development of phage-based therapies targeting these bacteria.

Analysis of Spounaviruses as a Case Study for the Overdue Reclassification of Tailed Phages

Tailed bacteriophages are the most abundant and diverse viruses in the world, with genome sizes ranging from 10 kbp to over 500 kbp. Yet, due to historical reasons, all this diversity is confined to a single virus order-Caudovirales, composed of just four families: Myoviridae, Siphoviridae, Podoviridae, and the newly created Ackermannviridae family. In recent years, this morphology-based classification scheme has started to crumble under the constant flood of phage sequences, revealing that tailed phages are even more genetically diverse than once thought. This prompted us, the Bacterial and Archaeal Viruses Subcommittee of the International Committee on Taxonomy of Viruses (ICTV), to consider overall reorganization of phage taxonomy. In this study, we used a wide range of complementary methods-including comparative genomics, core genome analysis, and marker gene phylogenetics-to show that the group of Bacillus phage SPO1-related viruses previously classified into the Spounavirinae subfamily, is clearly distinct from other members of the family Myoviridae and its diversity deserves the rank of an autonomous family. Thus, we removed this group from the Myoviridae family and created the family Herelleviridae-a new taxon of the same rank. In the process of the taxon evaluation, we explored the feasibility of different demarcation criteria and critically evaluated the usefulness of our methods for phage classification. The convergence of results, drawing a consistent and comprehensive picture of a new family with associated subfamilies, regardless of method, demonstrates that the tools applied here are particularly useful in phage taxonomy. We are convinced that creation of this novel family is a crucial milestone toward much-needed reclassification in the Caudovirales order.

Charting the diversity of uncultured viruses of Archaea and Bacteria

BACKGROUND

Viruses of Archaea and Bacteria are among the most abundant and diverse biological entities on Earth. Unraveling their biodiversity has been challenging due to methodological limitations. Recent advances in culture-independent techniques, such as metagenomics, shed light on the unknown viral diversity, revealing thousands of new viral nucleotide sequences at an unprecedented scale. However, these novel sequences have not been properly classified and the evolutionary associations between them were not resolved.

RESULTS

Here, we performed phylogenomic analysis of nearly 200,000 viral nucleotide sequences to establish GL-UVAB: Genomic Lineages of Uncultured Viruses of Archaea and Bacteria. The pan-genome content of the identified lineages shed light on some of their infection strategies, potential to modulate host physiology, and mechanisms to escape host resistance systems. Furthermore, using GL-UVAB as a reference database for annotating metagenomes revealed elusive habitat distribution patterns of viral lineages and environmental drivers of community composition.

CONCLUSIONS

These findings provide insights about the genomic diversity and ecology of viruses of prokaryotes. The source code used in these analyses is freely available at https://sourceforge.net/projects/gluvab/.

Charting the diversity of uncultured viruses of Archaea and Bacteria

Background

Viruses of Archaea and Bacteria are among the most abundant and diverse biological entities on Earth. Unraveling their biodiversity has been challenging due to methodological limitations. Recent advances in culture-independent techniques, such as metagenomics, shed light on the unknown viral diversity, revealing thousands of new viral nucleotide sequences at an unprecedented scale. However, these novel sequences have not been properly classified and the evolutionary associations between them were not resolved.

Results

Here, we performed phylogenomic analysis of nearly 200,000 viral nucleotide sequences to establish GL-UVAB: Genomic Lineages of Uncultured Viruses of Archaea and Bacteria. The pan-genome content of the identified lineages shed light on some of their infection strategies, potential to modulate host physiology, and mechanisms to escape host resistance systems. Furthermore, using GL-UVAB as a reference database for annotating metagenomes revealed elusive habitat distribution patterns of viral lineages and environmental drivers of community composition.

Conclusions

These findings provide insights about the genomic diversity and ecology of viruses of prokaryotes. The source code used in these analyses is freely available at https://sourceforge.net/projects/gluvab/.

A Novel Benthic Phage Infecting Shewanella with Strong Replication Ability

The coastal sediments were considered to contain diverse phages playing important roles in driving biogeochemical cycles based on genetic analysis. However, till now, benthic phages in coastal sediments were very rarely isolated, which largely limits our understanding of their biological characteristics. Here, we describe a novel lytic phage (named Shewanella phage S0112) isolated from the coastal sediments of the Yellow Sea infecting a sediment bacterium of the genus Shewanella. The phage has a very high replication capability, with the burst size of ca. 1170 phage particles per infected cell, which is 5–10 times higher than that of most phages isolated before. Meanwhile, the latent period of this phage is relatively longer, which might ensure adequate time for phage replication. The phage has a double-stranded DNA genome comprising 62,286 bp with 102 ORFs, ca. 60% of which are functionally unknown. The expression products of 16 ORF genes, mainly structural proteins, were identified by LC-MS/MS analysis. Besides the general DNA metabolism and structure assembly genes in the phage genome, there is a cluster of auxiliary metabolic genes that may be involved in 7-cyano-7-deazaguanine (preQ₀) biosynthesis. Meanwhile, a pyrophosphohydrolase (MazG) gene being considered as a regulator of programmed cell death or involving in host stringer responses is inserted in this gene cluster. Comparative genomic and phylogenetic analysis both revealed a great novelty of phage S0112. This study represents the first report of a benthic phage infecting Shewanella, which also sheds light on the phage–host interactions in coastal sediments.

Genome-resolved viral and cellular metagenomes revealed potential key virus-host interactions in a deep freshwater lake

Metagenomics has dramatically expanded the known virosphere, but freshwater viral diversity and their ecological interaction with hosts remain poorly understood. Here, we conducted a metagenomic exploration of planktonic dsDNA prokaryotic viruses by sequencing both virion (<0.22 μm) and cellular (0.22-5.0 μm) fractions collected spatiotemporally from a deep freshwater lake (Lake Biwa, Japan). This simultaneously reconstructed 183 complete (i.e., circular) viral genomes and 57 bacterioplankton metagenome-assembled genomes. Analysis of metagenomic read coverage revealed vertical partitioning of the viral community analogous to the vertically stratified bacterioplankton community. The hypolimnetic community was generally stable during stratification, but occasionally shifted abruptly, presumably due to lysogenic induction. Genes involved in assimilatory sulfate reduction were encoded in 20 (10.9%) viral genomes, including those of dominant viruses, and may aid viral propagation in sulfur-limited freshwater systems. Hosts were predicted for 40 (21.9%) viral genomes, encompassing 10 phyla (or classes of Proteobacteria) including ubiquitous freshwater bacterioplankton lineages (e.g., Ca. Fonsibacter and Ca. Nitrosoarchaeum). Comparison with viral genomes derived from published metagenomes revealed viral phylogeographic connectivity in geographically isolated habitats. Notably, analogous to their hosts, actinobacterial viruses were among the most diverse, ubiquitous and abundant viral groups in freshwater systems, with potential high lytic activity in surface waters.

Cooccurrence of Broad- and Narrow-Host-Range Viruses Infecting the Bloom-Forming Toxic Cyanobacterium Microcystis aeruginosa

Viruses play important roles in regulating the abundance and composition of bacterial populations in aquatic ecosystems. The bloom-forming toxic cyanobacterium Microcystis aeruginosa is predicted to interact with diverse cyanoviruses, resulting in Microcystis population diversification. However, current knowledge of the genomes from these viruses and their infection programs is limited to those of Microcystis virus Ma-LMM01. Here, we performed a time series sampling at a small pond in Japan during a Microcystis bloom and then investigated the genomic information and transcriptional dynamics of Microcystis-interacting viruses using metagenomic and metatranscriptomic approaches. We identified 15 viral genomic fragments classified into three groups, groups I (including Ma-LMM01), II (high abundance and transcriptional activity), and III (new lineages). According to the phylogenetic distribution of Microcystis strains possessing spacers against each viral group, the group II-original viruses interacted with all three phylogenetically distinct Microcystis population types (phylotypes), whereas the groups I and III-original viruses interacted with only one or two phylotypes, indicating the cooccurrence of broad- (group II) and narrow (groups I and III)-host-range viruses in the bloom. These viral fragments showed the highest transcriptional levels during daytime regardless of their genomic differences. Interestingly, M. aeruginosa expressed antiviral defense genes in the environment, unlike what was seen with an Ma-LMM01 infection in a previous culture experiment. Given that broad-host-range viruses often induce antiviral responses within alternative hosts, our findings suggest that such antiviral responses might inhibit viral multiplication, mainly that of broad-host-range viruses like those in group II.IMPORTANCE The bloom-forming toxic cyanobacterium Microcystis aeruginosa is thought to have diversified its population through the interactions between host and viruses in antiviral defense systems. However, current knowledge of viral genomes and infection programs is limited to those of Microcystis virus Ma-LMM01, which was a narrow host range in which it can escape from the highly abundant host defense systems. Our metagenomic approaches unveiled the cooccurrence of narrow- and broad-host-range Microcystis viruses, which included fifteen viral genomic fragments from Microcystis blooms that were classified into three groups. Interestingly, Microcystis antiviral defense genes were expressed against viral infection in the environment, unlike what was seen in a culture experiment with Ma-LMM01. Given that viruses with a broad host range often induce antiviral responses within alternative hosts, our findings suggest that antiviral responses inhibit viral reproduction, especially that of broad-range viruses like those in group II. This paper augments our understanding of the interactions between M. aeruginosa and its viruses and fills an important knowledge gap.

Viruses as key reservoirs of antibiotic resistance genes in the environment

Antibiotic resistance is a rapidly growing health care problem globally and causes many illnesses and deaths. Bacteria can acquire antibiotic resistance genes (ARGs) by horizontal transfer mediated by mobile genetic elements, where the role of phages in their dissemination in natural environments has not yet been clearly resolved. From metagenomic studies, we showed that the mean proportion of predicted ARGs found in prophages (0-0.0028%) was lower than those present in the free viruses (0.001-0.1%). Beta-lactamase, from viruses in the swine gut, represented 0.10 % of the predicted genes. Overall, in the environment, the ARG distribution associated with viruses was strongly linked to human activity, and the low dN/dS ratio observed advocated for a negative selection of the ARGs harbored by the viruses. Our network approach showed that viruses were linked to putative pathogens (Enterobacterales and vibrionaceae) and were considered key vehicles in ARG transfer, similar to plasmids. Therefore, these ARGs could then be disseminated at larger temporal and spatial scales than those included in the bacterial genomes, allowing for time-delayed genetic exchanges.

Transcriptome analysis of a thermophilic and hydrogenogenic carboxydotroph Carboxydothermus pertinax

A thermophilic and hydrogenogenic carboxydotroph, Carboxydothermus pertinax, performs hydrogenogenic CO metabolism in which CODH-II couples with distally encoded ECH. To enhance our knowledge of its hydrogenogenic CO metabolism, we performed whole transcriptome analysis of C. pertinax grown under 100% CO or 100% N2 using RNA sequencing. Of the 2577 genes, 36 and 64 genes were differentially expressed genes (DEGs) with false discovery rate adjusted P value < 0.05 when grown under 100% CO or 100% N2, respectively. Most of the DEGs were components of 23 gene clusters, suggesting switch between metabolisms via intensive expression changes in a relatively low number of gene clusters. Of the 9 significantly expressed gene clusters under 100% CO, CODH-II and ECH gene clusters were found. Only the ECH gene cluster was regulated by the CO-responsive transcriptional factor CooA, suggesting that others were separately regulated in the same transcriptional cascade as the ECH gene cluster. Of the 14 significantly expressed gene clusters under 100% N2, ferrous iron transport gene cluster involved in anaerobic respiration and prophage region were found. Considering that the expression of the temperate phage was strictly repressed under 100% CO, hydrogenogenic CO metabolism might be stable for C. pertinax.

Taxonomic assignment of uncultivated prokaryotic virus genomes is enabled by gene-sharing networks

Microbiomes from every environment contain a myriad of uncultivated archaeal and bacterial viruses, but studying these viruses is hampered by the lack of a universal, scalable taxonomic framework. We present vConTACT v.2.0, a network-based application utilizing whole genome gene-sharing profiles for virus taxonomy that integrates distance-based hierarchical clustering and confidence scores for all taxonomic predictions. We report near-identical (96%) replication of existing genus-level viral taxonomy assignments from the International Committee on Taxonomy of Viruses for National Center for Biotechnology Information virus RefSeq. Application of vConTACT v.2.0 to 1,364 previously unclassified viruses deposited in virus RefSeq as reference genomes produced automatic, high-confidence genus assignments for 820 of the 1,364. We applied vConTACT v.2.0 to analyze 15,280 Global Ocean Virome genome fragments and were able to provide taxonomic assignments for 31% of these data, which shows that our algorithm is scalable to very large metagenomic datasets. Our taxonomy tool can be automated and applied to metagenomes from any environment for virus classification.

Medusavirus, a Novel Large DNA Virus Discovered from Hot Spring Water

Recent discoveries of new large DNA viruses reveal high diversity in their morphologies, genetic repertoires, and replication strategies. Here, we report the novel features of medusavirus, a large DNA virus newly isolated from hot spring water in Japan. Medusavirus, with a diameter of 260 nm, shows a T=277 icosahedral capsid with unique spherical-headed spikes on its surface. It has a 381-kb genome encoding 461 putative proteins, 86 of which have their closest homologs in Acanthamoeba, whereas 279 (61%) are orphan genes. The virus lacks the genes encoding DNA topoisomerase II and RNA polymerase, showing that DNA replication takes place in the host nucleus, whereas the progeny virions are assembled in the cytoplasm. Furthermore, the medusavirus genome harbored genes for all five types of histones (H1, H2A, H2B, H3, and H4) and one DNA polymerase, which are phylogenetically placed at the root of the eukaryotic clades. In contrast, the host amoeba encoded many medusavirus homologs, including the major capsid protein. These facts strongly suggested that amoebae are indeed the most promising natural hosts of medusavirus, and that lateral gene transfers have taken place repeatedly and bidirectionally between the virus and its host since the early stage of their coevolution. Medusavirus reflects the traces of direct evolutionary interactions between the virus and eukaryotic hosts, which may be caused by sharing the DNA replication compartment and by evolutionarily long lasting virus-host relationships. Based on its unique morphological characteristics and phylogenomic relationships with other known large DNA viruses, we propose that medusavirus represents a new family, Medusaviridae IMPORTANCE We have isolated a new nucleocytoplasmic large DNA virus (NCLDV) from hot spring water in Japan, named medusavirus. This new NCLDV is phylogenetically placed at the root of the eukaryotic clades based on the phylogenies of several key genes, including that encoding DNA polymerase, and its genome surprisingly encodes the full set of histone homologs. Furthermore, its laboratory host, Acanthamoeba castellanii, encodes many medusavirus homologs in its genome, including the major capsid protein, suggesting that the amoeba is the genuine natural host from ancient times of this newly described virus and that lateral gene transfers have repeatedly occurred between the virus and amoeba. These results suggest that medusavirus is a unique NCLDV preserving ancient footprints of evolutionary interactions with its hosts, thus providing clues to elucidate the evolution of NCLDVs, eukaryotes, and virus-host interaction. Based on the dissimilarities with other known NCLDVs, we propose that medusavirus represents a new viral family, Medusaviridae.

Minimum Information about an Uncultivated Virus Genome (MIUViG)

We present an extension of the Minimum Information about any (x) Sequence (MIxS) standard for reporting sequences of uncultivated virus genomes. Minimum Information about an Uncultivated Virus Genome (MIUViG) standards were developed within the Genomic Standards Consortium framework and include virus origin, genome quality, genome annotation, taxonomic classification, biogeographic distribution and in silico host prediction. Community-wide adoption of MIUViG standards, which complement the Minimum Information about a Single Amplified Genome (MISAG) and Metagenome-Assembled Genome (MIMAG) standards for uncultivated bacteria and archaea, will improve the reporting of uncultivated virus genomes in public databases. In turn, this should enable more robust comparative studies and a systematic exploration of the global virosphere.

Potential Application of Bacteriophages in Enrichment Culture for Improved Prenatal Streptococcus agalactiae Screening

Vertical transmission of Streptococcus agalactiae can cause neonatal infections. A culture test in the late stage of pregnancy is used to screen for the presence of maternal S. agalactiae for intrapartum antibiotic prophylaxis. For the test, a vaginal–rectal sample is recommended to be enriched, followed by bacterial identification. In some cases, Enterococcus faecalis overgrows in the enrichment culture. Consequently, the identification test yields false-negative results. Bacteriophages (phages) can be used as antimicrobial materials. Here, we explored the feasibility of using phages to minimize false-negative results in an experimental setting. Phage mixture was prepared using three phages that specifically infect E. faecalis: phiEF24C, phiEF17H, and phiM1EF22. The mixture inhibited the growth of 86.7% (26/30) of vaginal E. faecalis strains. The simple coculture of E. faecalis and S. agalactiae was used as an experimental enrichment model. Phage mixture treatment led to suppression of E. faecalis growth and facilitation of S. agalactiae growth. In addition, testing several sets of S. agalactiae and E. faecalis strains, the treatment with phage mixture in the enrichment improved S. agalactiae detection on chromogenic agar. Our results suggest that the phage mixture can be usefully employed in the S. agalactiae culture test to increase test accuracy.

Biogeography and taxonomic overview of terrestrial hot spring thermophilic phages

Bacterial viruses ("phages") play important roles in the regulation and evolution of microbial communities in most ecosystems. Terrestrial hot springs typically contain thermophilic bacterial communities, but the diversity and impacts of its associated viruses ("thermophilic phages") are largely unexplored. Here, we provide a taxonomic overview of phages that have been isolated strictly from terrestrial hot springs around the world. In addition, we placed 17 thermophilic phage genomes in a global phylogenomic context to detect evolutionary patterns. Thermophilic phages have diverse morphologies (e.g., tailed, filamentous), unique virion structures (e.g., extremely long tailed siphoviruses), and span five taxonomic families encompassing strictly thermophilic phage genera. Within the phage proteomic tree, six thermophilic phage-related clades were identified, with evident genomic relatedness between thermophilic phages and archaeal viruses. Moreover, whole proteome analyses showed clustering between phages that infect distinct host phyla, such as Firmicutes and Deinococcus-Thermus. The potential for discovery of novel phage-host systems in terrestrial hot springs remain mostly untapped, thus additional emphasis on thermophilic phages in ecological prospecting is encouraged to gain insights into the microbial population dynamics of these environments.

ViCTree: an automated framework for taxonomic classification from protein sequences

Motivation

The increasing rate of submission of genetic sequences into public databases is providing a growing resource for classifying the organisms that these sequences represent. To aid viral classification, we have developed ViCTree, which automatically integrates the relevant sets of sequences in NCBI GenBank and transforms them into an interactive maximum likelihood phylogenetic tree that can be updated automatically. ViCTree incorporates ViCTreeView, which is a JavaScript-based visualization tool that enables the tree to be explored interactively in the context of pairwise distance data.

Results

To demonstrate utility, ViCTree was applied to subfamily Densovirinae of family Parvoviridae. This led to the identification of six new species of insect virus.

Availability and implementation

ViCTree is open-source and can be run on any Linux- or Unix-based computer or cluster. A tutorial, the documentation and the source code are available under a GPL3 license, and can be accessed at http://bioinformatics.cvr.ac.uk/victree_web/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Genome Sequences of 12 Mycobacteriophages Recovered from Archival Stocks in Japan

Using Mycobacterium smegmatis mc2155, 12 siphoviruses were recovered from long-term archival stocks stored in Japan. Their genome sequences were 46.0 to 61.3 kbp with 63 to 68% G+C contents, which allowed them to be categorized within cluster W and subclusters A1, A2, B3, A7, I1, and K4.