Assembly-free single-molecule sequencing recovers complete virus genomes from natural microbial communities

A review of neural networks for metagenomic binning

One of the main goals of metagenomic studies is to describe the taxonomic diversity of microbial communities. A crucial step in metagenomic analysis is metagenomic binning, which involves the (supervised) classification or (unsupervised) clustering of metagenomic sequences. Various machine learning models have been applied to address this task. In this review, the contributions of artificial neural networks (ANN) in the context of metagenomic binning are detailed, addressing both supervised, unsupervised, and semi-supervised approaches. 34 ANN-based binning tools are systematically compared, detailing their architectures, input features, datasets, advantages, disadvantages, and other relevant aspects. The findings reveal that deep learning approaches, such as convolutional neural networks and autoencoders, achieve higher accuracy and scalability than traditional methods. Gaps in benchmarking practices are highlighted, and future directions are proposed, including standardized datasets and optimization of architectures, for third-generation sequencing. This review provides support to researchers in identifying trends and selecting suitable tools for the metagenomic binning problem.

Distinct horizontal gene transfer potential of extracellular vesicles versus viral-like particles in marine habitats

Horizontal gene transfer (HGT) is enabled in part through the movement of DNA within two broad groups of small (<0.2 µm), diffusible nanoparticles: extracellular vesicles (EVs) and virus-like particles (VLPs; including viruses, gene transfer agents, and phage satellites). The information enclosed within these structures represents a substantial portion of the HGT potential available in planktonic ecosystems, but whether some genes might be preferentially transported through one type of nanoparticle versus another is unknown. Here we use long-read sequencing to compare the genetic content of EVs and VLPs from the oligotrophic North Pacific. Fractionated EV-enriched and VLP-enriched subpopulations contain diverse DNA from the surrounding microbial community, but differ in their capacity and encoded functions. The sequences carried by both particle types are enriched in mobile genetic elements (MGEs) as compared with other cellular chromosomal regions, and we highlight how this property enables novel MGE discovery. Examining the Pelagibacter mobilome reveals >7200 distinct chromosomal fragments and MGEs, many differentially partitioned between EVs and VLPs. Together these results suggest that distinctions in nanoparticle contents contribute to the mode and trajectory of microbial HGT networks and evolutionary dynamics in natural habitats.

Towards a unifying phylogenomic framework for tailed phages

Classifying viruses systematically has remained a key challenge of virology due to the absence of universal genes and vast genetic diversity of viruses. In particular, the most dominant and diverse group of viruses, the tailed double-stranded DNA viruses of prokaryotes belonging to the class Caudoviricetes, lack sufficient similarity in the genetic machinery that unifies them to reconstruct an inclusive, stable phylogeny of these genes. While previous approaches to organize tailed phage diversity have managed to distinguish various taxonomic levels, these methods are limited in scalability, reproducibility, and the inclusion of modes of evolution, like gene gains and losses, remain key challenges. Here, we present a novel, comprehensive, and reproducible framework for examining evolutionary relationships of tailed phages. In this framework, we compare phage genomes based on the presence and absence of a fixed set of gene families which are used as binary trait data that is input into maximum likelihood models. Our resulting phylogeny stably recovers known taxonomic families of tailed phages, with and without the inclusion of metagenome-derived phages. We also quantify the mosaicism of replication and structural genes among known families, and our results suggest that these exchanges likely underpin the emergence of new families. Additionally, we apply this framework to large phages (>100 kilobases) to map emergences of traits associated with genome expansion. Taken together, this evolutionary framework for charting and organizing tailed phage diversity improves the systemization of phage taxonomy, which can unify phage studies and advance our understanding of their evolution.

Virseqimprover: an integrated pipeline for viral contig error correction, extension, and annotation

Despite the recent surge of viral metagenomic studies, it remains a significant challenge to recover complete virus genomes from metagenomic data. The majority of viral contigs generated from de novo assembly programs are highly fragmented, presenting significant challenges to downstream analysis and inference. To address this issue, we have developed Virseqimprover, a computational pipeline that can extend assembled contigs to complete or nearly complete genomes while maintaining extension quality. Virseqimprover first examines whether there is any chimeric sequence based on read coverage, breaks the sequence into segments if there is, then extends the longest segment with uniform depth of coverage, and repeats these procedures until the sequence cannot be extended. Finally, Virseqimprover annotates the gene content of the resulting sequence. Results show that Virseqimprover has good performances on correcting and extending viral contigs to their full lengths, hence can be a useful tool to improve the completeness and minimize the assembly errors of viral contigs. Both a web server and a conda package for Virseqimprover are provided to the research community free of charge.

Complementary insights into gut viral genomes: a comparative benchmark of short- and long-read metagenomes using diverse assemblers and binners

BACKGROUND

Metagenome-assembled viral genomes have significantly advanced the discovery and characterization of the human gut virome. However, we lack a comparative assessment of assembly tools on the efficacy of viral genome identification, particularly across next-generation sequencing (NGS) and third-generation sequencing (TGS) data.

RESULTS

We evaluated the efficiency of NGS, TGS, and hybrid assemblers for viral genome discovery using 95 viral-like particle (VLP)-enriched fecal samples sequenced on both Illumina and PacBio platforms. MEGAHIT, metaFlye, and hybridSPAdes emerged as the optimal choices for NGS, TGS, and hybrid datasets, respectively. Notably, these assemblers recovered distinct viral genomes, demonstrating a remarkable degree of complementarity. By combining individual assembler results, we expanded the total number of nonredundant high-quality viral genomes by 4.83 ~ 21.7-fold compared to individual assemblers. Among them, viral genomes from NGS and TGS data have the least overlap, indicating the impact of data type on viral genome recovery. We also evaluated four binning methods, finding that CONCOCT incorporated more unrelated contigs into the same bins, while MetaBAT2, AVAMB, and vRhyme balanced inclusiveness and taxonomic consistency within bins.

CONCLUSIONS

Our findings highlight the challenges in metagenome-driven viral discovery, underscoring tool limitations. We advocate for combined use of multiple assemblers and sequencing technologies when feasible and highlight the urgent need for specialized tools tailored to gut virome assembly. This study contributes essential insights for advancing viral genome research in the context of gut metagenomics. Video Abstract.

Method evaluation for viruses in activated sludge: Concentration, sequencing, and identification

Activated sludge (AS) in wastewater treatment plants is one of the largest artificial microbial ecosystems on earth and it makes enormous contributions to human societies. Viruses are an important component in AS with a high abundance. However, their communities and functionalities have not been as widely explored as those of other microorganisms, such as bacteria. This gap is mainly due to technical challenges in effective viral concentration, extraction, and sequencing. In this study, we compared four kinds of concentration methods, two sequencing approaches, and four identification bioinformatic tools to evaluate the whole analysis workflow for viruses in AS. Results showed flocculation, filtration, and resuspension (FFR) could get the longest DNA lengths and ultracentrifugation obtained the highest DNA yields for viruses in AS. Based on the results of present study, FFR and tangential flow filtration with the membrane pore size of 100 kDa were most recommended to concentrate viruses in AS samples with huge volumes. Besides, different concentration methods could get different viral catalogs and thus multiple methods should be combined to get the whole picture of viruses in the system. In addition, geNomad was the most recommended identification tool for viruses in the present study and the long-read sequencing could improve the assembly statistics of viruses when compared with the short-read sequencing. For the 8192 viral operational taxonomic units in this study, 95.1 % of them were phages and belonged to the same lineage at the order level of Caudovirales. Virulent phages dominated the AS system and Pseudomonadota were the main host. Taken together, this study provides new insights into methods selection for virus research of AS.

Comparative Evaluation of Open-Source Bioinformatics Pipelines for Full-Length Viral Genome Assembly

The increasingly widespread application of next-generation sequencing (NGS) in clinical diagnostics and epidemiological research has generated a demand for robust, fast, automated, and user-friendly bioinformatics workflows. To guide the choice of tools for the assembly of full-length viral genomes from NGS datasets, we assessed the performance and applicability of four open-source bioinformatics pipelines (shiver-for which we created a user-friendly Dockerized version, referred to as dshiver; SmaltAlign; viral-ngs; and V-pipe) using both simulated and real-world HIV-1 paired-end short-read datasets and default settings. All four pipelines produced consensus genome assemblies with high quality metrics (genome fraction recovery, mismatch and indel rates, variant calling F1 scores) when the reference sequence used for assembly had high similarity to the analyzed sample. The shiver and SmaltAlign pipelines (but not viral-ngs and V-Pipe) also showed robust performance with more divergent samples (non-matching subtypes). With empirical datasets, SmaltAlign and viral-ngs exhibited an order of magnitude shorter runtime compared to V-Pipe and shiver. In terms of applicability, V-Pipe provides the broadest functionalities, SmaltAlign and dshiver combine user-friendliness with robustness, while the use of viral-ngs requires less computational resources compared to other pipelines. In conclusion, if a closely matched reference sequence is available, all pipelines can reliably reconstruct viral consensus genomes; therefore, differences in user-friendliness and runtime may guide the choice of the pipeline in a particular setting. If a matched reference sequence cannot be selected, we recommend shiver or SmaltAlign for robust performance. The new Dockerized version of shiver offers ease of use in addition to the accuracy and robustness of the original pipeline.

Strainy: phasing and assembly of strain haplotypes from long-read metagenome sequencing

Bacterial species in microbial communities are often represented by mixtures of strains, distinguished by small variations in their genomes. Short-read approaches can be used to detect small-scale variation between strains but fail to phase these variants into contiguous haplotypes. Long-read metagenome assemblers can generate contiguous bacterial chromosomes but often suppress strain-level variation in favor of species-level consensus. Here we present Strainy, an algorithm for strain-level metagenome assembly and phasing from Nanopore and PacBio reads. Strainy takes a de novo metagenomic assembly as input and identifies strain variants, which are then phased and assembled into contiguous haplotypes. Using simulated and mock Nanopore and PacBio metagenome data, we show that Strainy assembles accurate and complete strain haplotypes, outperforming current Nanopore-based methods and comparable with PacBio-based algorithms in completeness and accuracy. We then use Strainy to assemble strain haplotypes of a complex environmental metagenome, revealing distinct strain distribution and mutational patterns in bacterial species.

The coral Oculina patagonica holobiont and its response to confinement, temperature, and Vibrio infections

BACKGROUND

Extensive research on the diversity and functional roles of the microorganisms associated with reef-building corals has been promoted as a consequence of the rapid global decline of coral reefs attributed to climate change. Several studies have highlighted the importance of coral-associated algae (Symbiodinium) and bacteria and their potential roles in promoting coral host fitness and survival. However, the complex coral holobiont extends beyond these components to encompass other entities such as protists, fungi, and viruses. While each constituent has been individually investigated in corals, a comprehensive understanding of their collective roles is imperative for a holistic comprehension of coral health and resilience.

RESULTS

The metagenomic analysis of the microbiome of the coral Oculina patagonica has revealed that fungi of the genera Aspergillus, Fusarium, and Rhizofagus together with the prokaryotic genera Streptomyces, Pseudomonas, and Bacillus were abundant members of the coral holobiont. This study also assessed changes in microeukaryotic, prokaryotic, and viral communities under three stress conditions: aquaria confinement, heat stress, and Vibrio infections. In general, stress conditions led to an increase in Rhodobacteraceae, Flavobacteraceae, and Vibrionaceae families, accompanied by a decrease in Streptomycetaceae. Concurrently, there was a significant decline in both the abundance and richness of microeukaryotic species and a reduction in genes associated with antimicrobial compound production by the coral itself, as well as by Symbiodinium and fungi.

CONCLUSION

Our findings suggest that the interplay between microeukaryotic and prokaryotic components of the coral holobiont may be disrupted by stress conditions, such as confinement, increase of seawater temperature, or Vibrio infection, leading to a dysbiosis in the global microbial community that may increase coral susceptibility to diseases. Further, microeukaryotic community seems to exert influence on the prokaryotic community dynamics, possibly through predation or the production of secondary metabolites with anti-bacterial activity. Video Abstract.

Decomposing a San Francisco estuary microbiome using long-read metagenomics reveals species- and strain-level dominance from picoeukaryotes to viruses

Although long-read sequencing has enabled obtaining high-quality and complete genomes from metagenomes, many challenges still remain to completely decompose a metagenome into its constituent prokaryotic and viral genomes. This study focuses on decomposing an estuarine metagenome to obtain a more accurate estimate of microbial diversity. To achieve this, we developed a new bead-based DNA extraction method, a novel bin refinement method, and obtained 150 Gbp of Nanopore sequencing. We estimate that there are ~500 bacterial and archaeal species in our sample and obtained 68 high-quality bins (>90% complete, <5% contamination, ≤5 contigs, contig length of >100 kbp, and all ribosomal and tRNA genes). We also obtained many contigs of picoeukaryotes, environmental DNA of larger eukaryotes such as mammals, and complete mitochondrial and chloroplast genomes and detected ~40,000 viral populations. Our analysis indicates that there are only a few strains that comprise most of the species abundances.

IMPORTANCE

Ocean and estuarine microbiomes play critical roles in global element cycling and ecosystem function. Despite the importance of these microbial communities, many species still have not been cultured in the lab. Environmental sequencing is the primary way the function and population dynamics of these communities can be studied. Long-read sequencing provides an avenue to overcome limitations of short-read technologies to obtain complete microbial genomes but comes with its own technical challenges, such as needed sequencing depth and obtaining high-quality DNA. We present here new sampling and bioinformatics methods to attempt decomposing an estuarine microbiome into its constituent genomes. Our results suggest there are only a few strains that comprise most of the species abundances from viruses to picoeukaryotes, and to fully decompose a metagenome of this diversity requires 1 Tbp of long-read sequencing. We anticipate that as long-read sequencing technologies continue to improve, less sequencing will be needed.

Targeted whole-genome recovery of single viral species in a complex environmental sample

Characterizing unknown viruses is essential for understanding viral ecology and preparing against viral outbreaks. Recovering complete genome sequences from environmental samples remains computationally challenging using metagenomics, especially for low-abundance species with uneven coverage. We present an experimental method for reliably recovering complete viral genomes from complex environmental samples. Individual genomes are encapsulated into droplets and amplified using multiple displacement amplification. A unique gene detection assay, which employs an RNA-based probe and an exonuclease, selectively identifies droplets containing the target viral genome. Labeled droplets are sorted using a microfluidic sorter, and genomes are extracted for sequencing. We demonstrate this method's efficacy by spiking two known viral genomes, Simian virus 40 (SV40, 5,243 bp) and Human Adenovirus 5 (HAd5, 35,938 bp), into a sewage sample with a final abundance in the droplets of around 0.1% and 0.015%, respectively. We achieve 100% recovery of the complete sequence of the spiked-in SV40 genome with uniform coverage distribution. For the larger HAd5 genome, we cover approximately 99.4% of its sequence. Notably, genome recovery is achieved with as few as one sorted droplet, which enables the recovery of any desired genomes in complex environmental samples, regardless of their abundance. This method enables single-genome whole-genome amplification and targeting characterizations of rare viral species and will facilitate our ability to access the mutational profile in single-virus genomes and contribute to an improved understanding of viral ecology.

A Review of Probe-Based Enrichment Methods to Inform Plant Virus Diagnostics

Modern diagnostic techniques based on DNA sequence similarity are currently the gold standard for the detection of existing and emerging pathogens. Whilst individual assays are inexpensive to use, assay development is costly and carries risks of not being sensitive or specific enough to capture an increasingly diverse range of targets. Sequencing can provide the entire nucleic acid content of a sample and may be used to identify all pathogens present in the sample when the depth of coverage is sufficient. Targeted enrichment techniques have been used to increase sequence coverage and improve the sensitivity of detection within virus samples, specifically, to capture sequences for a range of different viruses or increase the number of reads from low-titre virus infections. Vertebrate viruses have been well characterised using in-solution hybridisation capture to target diverse virus families. The use of probes for genotyping and strain identification has been limited in plants, and uncertainty around sensitivity is an impediment to the development of a large-scale virus panel to use within regulatory settings and diagnostic pipelines. This review aims to compare significant studies that have used targeted enrichment of viruses to identify approaches to probe design and potential for use in plant virus detection and characterisation.

Assessing RNA integrity by digital RT-PCR: Influence of extraction, storage, and matrices

The development of high-throughput sequencing has greatly improved our knowledge of microbial diversity in aquatic environments and its evolution in highly diverse ecosystems. Relevant microbial diversity description based on high-throughput sequencing relies on the good quality of the nucleic acid recovered. Indeed, long genetic fragments are more informative for identifying mutation combinations that characterize variants or species in complex samples. This study describes a new analytical method based on digital Polymerase Chain Reaction (PCR) partitioning technology for assessing the fragmentation of nucleic acid and more specifically viral RNA. This method allows us to overcome limits associated with hydrolysis probe-based assay by focusing on the distance between different amplicons, and not, as usual, on the size of amplicons. RNA integrity can thus be determined as a new fragmentation index, the so-called Fragment size 50. The application of this method has provided information on issues that are inherent in environmental analyses, such as the storage impact of raw samples or extracted RNA, extraction methods, and the nature of the sample on the integrity of viral RNA. Finally, the estimation of fragment size by digital PCR (dPCR) showed a very strong similarity with the fragment size sequenced using Oxford Nanopore Technology. In addition to enabling objective improvements in analytical methods, this approach could become a systematic quality control prior to any long-read sequencing, avoiding insufficiently productive sequencing runs or biases in the representativeness of sequenced fragments.

Complete genome sequence of a Pseudomonas fluorescens bacteriophage UNO-G1W1 isolated from freshwater ice in Nebraska

We provide the complete genome sequence for a novel Pseudomonas fluorescens bacteriophage named UNO-G1W1. This phage was isolated from a single ice cover sampling. The genome was sequenced on the Nanopore MinION, generated with the direct terminal repeat-phage-pipeline and polished with Illumina short reads. Sequence identity classifies the phage as an otagovirus.

Sensitive and error-tolerant annotation of protein-coding DNA with BATH

Summary

We present BATH, a tool for highly sensitive annotation of protein-coding DNA based on direct alignment of that DNA to a database of protein sequences or profile hidden Markov models (pHMMs). BATH is built on top of the HMMER3 code base, and simplifies the annotation workflow for pHMM-based translated sequence annotation by providing a straightforward input interface and easy-to-interpret output. BATH also introduces novel frameshift-aware algorithms to detect frameshift-inducing nucleotide insertions and deletions (indels). BATH matches the accuracy of HMMER3 for annotation of sequences containing no errors, and produces superior accuracy to all tested tools for annotation of sequences containing nucleotide indels. These results suggest that BATH should be used when high annotation sensitivity is required, particularly when frameshift errors are expected to interrupt protein-coding regions, as is true with long-read sequencing data and in the context of pseudogenes.

Availability and implementation

The software is available at https://github.com/TravisWheelerLab/BATH.

Nanopore and Illumina sequencing reveal different viral populations from human gut samples

The advent of viral metagenomics, or viromics, has improved our knowledge and understanding of global viral diversity. High-throughput sequencing technologies enable explorations of the ecological roles, contributions to host metabolism, and the influence of viruses in various environments, including the human intestinal microbiome. However, bacterial metagenomic studies frequently have the advantage. The adoption of advanced technologies like long-read sequencing has the potential to be transformative in refining viromics and metagenomics. Here, we examined the effectiveness of long-read and hybrid sequencing by comparing Illumina short-read and Oxford Nanopore Technology (ONT) long-read sequencing technologies and different assembly strategies on recovering viral genomes from human faecal samples. Our findings showed that if a single sequencing technology is to be chosen for virome analysis, Illumina is preferable due to its superior ability to recover fully resolved viral genomes and minimise erroneous genomes. While ONT assemblies were effective in recovering viral diversity, the challenges related to input requirements and the necessity for amplification made it less ideal as a standalone solution. However, using a combined, hybrid approach enabled a more authentic representation of viral diversity to be obtained within samples.

Characterization and genomic analysis of the Lyme disease spirochete bacteriophage ϕBB-1

Lyme disease is a tick-borne infection caused by the spirochete Borrelia (Borreliella) burgdorferi. Borrelia species have highly fragmented genomes composed of a linear chromosome and a constellation of linear and circular plasmids some of which are required throughout the enzootic cycle. Included in this plasmid repertoire by almost all Lyme disease spirochetes are the 32-kb circular plasmid cp32 prophages that are capable of lytic replication to produce infectious virions called ϕBB-1. While the B. burgdorferi genome contains evidence of horizontal transfer, the mechanisms of gene transfer between strains remain unclear. While we know that ϕBB-1 transduces cp32 and shuttle vector DNA during in vitro cultivation, the extent of ϕBB-1 DNA transfer is not clear. Herein, we use proteomics and long-read sequencing to further characterize ϕBB-1 virions. Our studies identified the cp32 pac region and revealed that ϕBB-1 packages linear cp32s via a headful mechanism with preferential packaging of plasmids containing the cp32 pac region. Additionally, we find ϕBB-1 packages fragments of the linear chromosome and full-length plasmids including lp54, cp26, and others. Furthermore, sequencing of ϕBB-1 packaged DNA allowed us to resolve the covalently closed hairpin telomeres for the linear B. burgdorferi chromosome and most linear plasmids in strain CA-11.2A. Collectively, our results shed light on the biology of the ubiquitous ϕBB-1 phage and further implicates ϕBB-1 in the generalized transduction of diverse genes and the maintenance of genetic diversity in Lyme disease spirochetes.

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.

COBRA improves the completeness and contiguity of viral genomes assembled from metagenomes

Viruses are often studied using metagenome-assembled sequences, but genome incompleteness hampers comprehensive and accurate analyses. Contig Overlap Based Re-Assembly (COBRA) resolves assembly breakpoints based on the de Bruijn graph and joins contigs. Here we benchmarked COBRA using ocean and soil viral datasets. COBRA accurately joined the assembled sequences and achieved notably higher genome accuracy than binning tools. From 231 published freshwater metagenomes, we obtained 7,334 bacteriophage clusters, ~83% of which represent new phage species. Notably, ~70% of these were circular, compared with 34% before COBRA analyses. We expanded sampling of huge phages (≥200 kbp), the largest of which was curated to completion (717 kbp). Improved phage genomes from Rotsee Lake provided context for metatranscriptomic data and indicated the in situ activity of huge phages, whiB-encoding phages and cysC- and cysH-encoding phages. COBRA improves viral genome assembly contiguity and completeness, thus the accuracy and reliability of analyses of gene content, diversity and evolution.

The long and short of it: benchmarking viromics using Illumina, Nanopore and PacBio sequencing technologies

Viral metagenomics has fuelled a rapid change in our understanding of global viral diversity and ecology. Long-read sequencing and hybrid assembly approaches that combine long- and short-read technologies are now being widely implemented in bacterial genomics and metagenomics. However, the use of long-read sequencing to investigate viral communities is still in its infancy. While Nanopore and PacBio technologies have been applied to viral metagenomics, it is not known to what extent different technologies will impact the reconstruction of the viral community. Thus, we constructed a mock bacteriophage community of previously sequenced phage genomes and sequenced them using Illumina, Nanopore and PacBio sequencing technologies and tested a number of different assembly approaches. When using a single sequencing technology, Illumina assemblies were the best at recovering phage genomes. Nanopore- and PacBio-only assemblies performed poorly in comparison to Illumina in both genome recovery and error rates, which both varied with the assembler used. The best Nanopore assembly had errors that manifested as SNPs and INDELs at frequencies 41 and 157 % higher than found in Illumina only assemblies, respectively. While the best PacBio assemblies had SNPs at frequencies 12 and 78 % higher than found in Illumina-only assemblies, respectively. Despite high-read coverage, long-read-only assemblies recovered a maximum of one complete genome from any assembly, unless reads were down-sampled prior to assembly. Overall the best approach was assembly by a combination of Illumina and Nanopore reads, which reduced error rates to levels comparable with short-read-only assemblies. When using a single technology, Illumina only was the best approach. The differences in genome recovery and error rates between technology and assembler had downstream impacts on gene prediction, viral prediction, and subsequent estimates of diversity within a sample. These findings will provide a starting point for others in the choice of reads and assembly algorithms for the analysis of viromes.

Characterization and genomic analysis of the Lyme disease spirochete bacteriophage ϕBB-1

Lyme disease is a tick-borne infection caused by the spirochete Borrelia (Borreliella) burgdorferi. Borrelia species have highly fragmented genomes composed of a linear chromosome and a constellation of linear and circular plasmids some of which are required throughout the enzootic cycle. Included in this plasmid repertoire by almost all Lyme disease spirochetes are the 32-kb circular plasmid cp32 prophages that are capable of lytic replication to produce infectious virions called ϕBB-1. While the B. burgdorferi genome contains evidence of horizontal transfer, the mechanisms of gene transfer between strains remain unclear. While we know that ϕBB-1 transduces cp32 and shuttle vector DNA during in vitro cultivation, the extent of ϕBB-1 DNA transfer is not clear. Herein, we use proteomics and long-read sequencing to further characterize ϕBB-1 virions. Our studies identified the cp32 pac region and revealed that ϕBB-1 packages linear cp32s via a headful mechanism with preferentially packaging of plasmids containing the cp32 pac region. Additionally, we find ϕBB-1 packages fragments of the linear chromosome and full-length plasmids including lp54, cp26, and others. Furthermore, sequencing of ϕBB-1 packaged DNA allowed us to resolve the covalently closed hairpin telomeres for the linear B. burgdorferi chromosome and most linear plasmids in strain CA-11.2A. Collectively, our results shed light on the biology of the ubiquitous ϕBB-1 phage and further implicates ϕBB-1 in the generalized transduction of diverse genes and the maintenance of genetic diversity in Lyme disease spirochetes.

Sensitive and error-tolerant annotation of protein-coding DNA with BATH

We present BATH, a tool for highly sensitive annotation of protein-coding DNA based on direct alignment of that DNA to a database of protein sequences or profile hidden Markov models (pHMMs). BATH is built on top of the HMMER3 code base, and simplifies the annotation workflow for pHMM-based annotation by providing a straightforward input interface and easy-to-interpret output. BATH also introduces novel frameshift-aware algorithms to detect frameshift-inducing nucleotide insertions and deletions (indels). BATH matches the accuracy of HMMER3 for annotation of sequences containing no errors, and produces superior accuracy to all tested tools for annotation of sequences containing nucleotide indels. These results suggest that BATH should be used when high annotation sensitivity is required, particularly when frameshift errors are expected to interrupt protein-coding regions, as is true with long read sequencing data and in the context of pseudogenes.

VIGA: a one-stop tool for eukaryotic virus identification and genome assembly from next-generation-sequencing data

Identification of viruses and further assembly of viral genomes from the next-generation-sequencing data are essential steps in virome studies. This study presented a one-stop tool named VIGA (available at https://github.com/viralInformatics/VIGA) for eukaryotic virus identification and genome assembly from NGS data. It was composed of four modules, namely, identification, taxonomic annotation, assembly and novel virus discovery, which integrated several third-party tools such as BLAST, Trinity, MetaCompass and RagTag. Evaluation on multiple simulated and real virome datasets showed that VIGA assembled more complete virus genomes than its competitors on both the metatranscriptomic and metagenomic data and performed well in assembling virus genomes at the strain level. Finally, VIGA was used to investigate the virome in metatranscriptomic data from the Human Microbiome Project and revealed different composition and positive rate of viromes in diseases of prediabetes, Crohn's disease and ulcerative colitis. Overall, VIGA would help much in identification and characterization of viromes, especially the known viruses, in future studies.

Viral Diversity in Benthic Abyssal Ecosystems: Ecological and Methodological Considerations

Viruses are the most abundant 'biological entities' in the world's oceans. However, technical and methodological constraints limit our understanding of their diversity, particularly in benthic abyssal ecosystems (>4000 m depth). To verify advantages and limitations of analyzing virome DNA subjected either to random amplification or unamplified, we applied shotgun sequencing-by-synthesis to two sample pairs obtained from benthic abyssal sites located in the North-eastern Atlantic Ocean at ca. 4700 m depth. One amplified DNA sample was also subjected to single-molecule long-read sequencing for comparative purposes. Overall, we identified 24,828 viral Operational Taxonomic Units (vOTUs), belonging to 22 viral families. Viral reads were more abundant in the amplified DNA samples (38.5-49.9%) compared to the unamplified ones (4.4-5.8%), with the latter showing a greater viral diversity and 11-16% of dsDNA viruses almost undetectable in the amplified samples. From a procedural point of view, the viromes obtained by direct sequencing (without amplification step) provided a broader overview of both ss and dsDNA viral diversity. Nevertheless, our results suggest that the contextual use of random amplification of the same sample and long-read technology can improve the assessment of viral assemblages by reducing off-target reads.

Effects of microplastics on functional genes related to CH4 and N2O metabolism in bacteriophages during manure composting and its planting applications

Microplastics (MPs), as a new type of pollutant, widely exist in livestock and poultry breeding and agricultural soils. However, research on MPs pollution on greenhouse gas emissions in combined planting and breeding systems is lacking, especially from the perspective of phage horizontal gene transfer. Therefore, this paper explores the effects of MPs on functional genes related to CH4 and N2O metabolism in bacteriophages during manure composting and its planting applications. The results of the study indicated that the addition of MPs had an impact on both the physicochemical properties and microbial community structure of manure during the composting process and on the compost-applied rhizosphere soil of lactuca (Lactuca sativa). Specifically, on day 7 of composting, mcrA/pmoA and (nirS+nirK) levels in bacteria in the MP group significantly increased. Additionally, it was observed that the MP group had higher average temperatures during the high-temperature period of composting, which led to a rapid reduction in phages. However, the phage levels quickly recovered during the cooling period. Furthermore, the addition of MPs to the rhizosphere soil resulted in higher levels of nirK. These changes may affect greenhouse gas emissions.

Genomic diversity and biogeographic distributions of a novel lineage of bacteriophages that infect marine OM43 bacteria

The marine methylotrophic OM43 clade is considered an important bacterial group in coastal microbial communities. OM43 bacteria, which are closely related to phytoplankton blooms, have small cell sizes and streamlined genomes. Bacteriophages profoundly shape the evolutionary trajectories, population dynamics, and physiology of microbes. The prevalence and diversity of several phages that infect OM43 bacteria have been reported. In this study, we isolated and sequenced two novel OM43 phages, MEP401 and MEP402. These phages share 90% of their open reading frames (ORFs) and are distinct from other known phage isolates. Furthermore, a total of 99 metagenomic viral genomes (MVGs) closely related to MEP401 and MEP402 were identified. Phylogenomic analyses suggest that MEP401, MEP402, and these identified MVGs belong to a novel subfamily in the family Zobellviridae and that they can be separated into two groups. Group I MVGs show conserved whole-genome synteny with MEP401, while group II MVGs possess the MEP401-type DNA replication module and a distinct type of morphogenesis and packaging module, suggesting that genomic recombination occurred between phages. Most members in these two groups were predicted to infect OM43 bacteria. Metagenomic read-mapping analysis revealed that the phages in these two groups are globally ubiquitous and display distinct biogeographic distributions, with some phages being predominant in cold regions, some exclusively detected in estuarine stations, and others displaying wider distributions. This study expands our knowledge of the diversity and ecology of a novel phage lineage that infects OM43 bacteria by describing their genomic diversity and global distribution patterns. IMPORTANCE OM43 phages that infect marine OM43 bacteria are important for host mortality, community structure, and physiological functions. In this study, two OM43 phages were isolated and characterized. Metagenomic viral genome (MVG) retrieval using these two OM43 phages as baits led to the identification of two phage groups of a new subfamily in the family Zobellviridae. We found that group I MVGs share similar genomic content and arrangement with MEP401 and MEP402, whereas group II MVGs only possess the MEP401-type DNA replication module. Metagenomic mapping analysis suggests that members in these two groups are globally ubiquitous with distinct distribution patterns. This study provides important insights into the genomic diversity and biogeography of the OM43 phages in the global ocean.

Virus impacted community adaptation in oligotrophic groundwater environment revealed by Hi-C coupled metagenomic and viromic study

Viruses play a crucial role in microbial mortality, diversity and biogeochemical cycles. Groundwater is the largest global freshwater and one of the most oligotrophic aquatic systems on Earth, but how microbial and viral communities are shaped in this special habitat is largely unexplored. In this study, we collected groundwater samples from 23 to 60 m aquifers at Yinchuan Plain, China. In total, 1920 non-reductant viral contigs were retrieved from metagenomes and viromes constructed by Illumina and Nanopore hybrid sequencing. Only 3% of them could be clustered with known viruses, most of which were Caudoviricetes. Coupling 1.2 Tb Hi-C sequencing with CRISPR matching and homology search, we connected 469 viruses with their hosts while some viral clusters presented a broad-host-range trait. Meanwhile, a large proportion of biosynthesis related auxiliary metabolism genes were identified. Those characteristics might benefit viruses for a better survival in this special oligotrophic environment. Additionally, the groundwater virome showed genomic features distinct from those of the open ocean and wastewater treatment facilities in GC distribution and unannotated gene compositions. This paper expands the current knowledge of the global viromic records and serves as a foundation for a more thorough understanding of viruses in groundwater.

Environmental DNA/RNA for pathogen and parasite detection, surveillance, and ecology

Detection of pathogens, parasites, and other symbionts in environmental samples via eDNA/eRNA (collectively eNA) is an increasingly important source of information about their occurrence and activity. There is great potential for using such detections as a proxy for infection of host organisms in connected habitats, for pathogen monitoring and surveillance, and for early warning systems for disease. However, many factors require consideration, and appropriate methods developed and verified, in order that eNA detections can be reliably interpreted and adopted for surveillance and assessment of disease risk, and potentially inclusion in international standards, such as the World Organisation for Animal Health guidelines. Disease manifestation results from host-symbiont-environment interactions between hosts, demanding a multifactorial approach to interpretation of eNA signals.

Integrated molecular approaches for fermented food microbiome research

Molecular technologies, including high-throughput sequencing, have expanded our perception of the microbial world. Unprecedented insights into the composition and function of microbial communities have generated large interest, with numerous landmark studies published in recent years relating the important roles of microbiomes and the environment-especially diet and nutrition-in human, animal, and global health. As such, food microbiomes represent an important cross-over between the environment and host. This is especially true of fermented food microbiomes, which actively introduce microbial metabolites and, to a lesser extent, live microbes into the human gut. Here, we discuss the history of fermented foods, and examine how molecular approaches have advanced research of these fermented foods over the past decade. We highlight how various molecular approaches have helped us to understand the ways in which microbes shape the qualities of these products, and we summarize the impacts of consuming fermented foods on the gut. Finally, we explore how advances in bioinformatics could be leveraged to enhance our understanding of fermented foods. This review highlights how integrated molecular approaches are changing our understanding of the microbial communities associated with food fermentation, the creation of unique food products, and their influences on the human microbiome and health.

FastViromeExplorer-Novel: Recovering Draft Genomes of Novel Viruses and Phages in Metagenomic Data

Despite the recent surge of viral metagenomic studies, recovering complete virus/phage genomes from metagenomic data is still extremely difficult and most viral contigs generated from de novo assembly programs are highly fragmented, posing serious challenges to downstream analysis and inference. In this study, we develop FastViromeExplorer (FVE)-novel, a computational pipeline for reconstructing complete or near-complete viral draft genomes from metagenomic data. The FVE-novel deploys FVE to efficiently map metagenomic reads to viral reference genomes, performs de novo assembly of the mapped reads to generate contigs, and extends the contigs through iterative assembly to produce final viral scaffolds. We applied FVE-novel to an ocean metagenomic sample and obtained 268 viral scaffolds that potentially come from novel viruses. Through manual examination and validation of the 10 longest scaffolds, we successfully recovered 4 complete viral genomes, 2 are novel as they cannot be found in the existing databases and the other 2 are related to known phages. This hybrid reference-based and de novo assembly approach used by FVE-novel represents a powerful new approach for uncovering near-complete viral genomes in metagenomic data.

Novel integrative elements and genomic plasticity in ocean ecosystems

Horizontal gene transfer accelerates microbial evolution. The marine picocyanobacterium Prochlorococcus exhibits high genomic plasticity, yet the underlying mechanisms are elusive. Here, we report a novel family of DNA transposons-"tycheposons"-some of which are viral satellites while others carry cargo, such as nutrient-acquisition genes, which shape the genetic variability in this globally abundant genus. Tycheposons share distinctive mobile-lifecycle-linked hallmark genes, including a deep-branching site-specific tyrosine recombinase. Their excision and integration at tRNA genes appear to drive the remodeling of genomic islands-key reservoirs for flexible genes in bacteria. In a selection experiment, tycheposons harboring a nitrate assimilation cassette were dynamically gained and lost, thereby promoting chromosomal rearrangements and host adaptation. Vesicles and phage particles harvested from seawater are enriched in tycheposons, providing a means for their dispersal in the wild. Similar elements are found in microbes co-occurring with Prochlorococcus, suggesting a common mechanism for microbial diversification in the vast oligotrophic oceans.

Genomic analysis and characterization of phages infecting the marine Roseobacter CHAB-I-5 lineage reveal a globally distributed and abundant phage genus

Marine phages play an important role in marine biogeochemical cycles by regulating the death, physiological metabolism, and evolutionary trajectory of bacteria. The Roseobacter group is an abundant and important heterotrophic bacterial group in the ocean, and plays an important role in carbon, nitrogen, sulfur and phosphorus cycling. The CHAB-I-5 lineage is one of the most dominant Roseobacter lineages, but remains largely uncultured. Phages infecting CHAB-I-5 bacteria have not yet been investigated due to the lack of culturable CHAB-I-5 strains. In this study, we isolated and sequenced two new phages (CRP-901 and CRP-902) infecting the CHAB-I-5 strain FZCC0083. We applied metagenomic data mining, comparative genomics, phylogenetic analysis, and metagenomic read-mapping to investigate the diversity, evolution, taxonomy, and biogeography of the phage group represented by the two phages. The two phages are highly similar, with an average nucleotide identity of 89.17%, and sharing 77% of their open reading frames. We identified several genes involved in DNA replication and metabolism, virion structure, DNA packing, and host lysis from their genomes. Metagenomic mining identified 24 metagenomic viral genomes closely related to CRP-901 and CRP-902. Genomic comparison and phylogenetic analysis demonstrated that these phages are distinct from other known viruses, representing a novel genus-level phage group (CRP-901-type). The CRP-901-type phages do not contain DNA primase and DNA polymerase genes, but possess a novel bifunctional DNA primase-polymerase gene with both primase and polymerase activities. Read-mapping analysis showed that the CRP-901-type phages are widespread across the world's oceans and are most abundant in estuarine and polar waters. Their abundance is generally higher than other known roseophages and even higher than most pelagiphages in the polar region. In summary, this study has greatly expanded our understanding of the genetic diversity, evolution, and distribution of roseophages. Our analysis suggests that the CRP-901-type phage is an important and novel marine phage group that plays important roles in the physiology and ecology of roseobacters.

What the Phage: a scalable workflow for the identification and analysis of phage sequences

Phages are among the most abundant and diverse biological entities on earth. Phage prediction from sequence data is a crucial first step to understanding their impact on the environment. A variety of bacteriophage prediction tools have been developed over the years. They differ in algorithmic approach, results, and ease of use. We, therefore, developed "What the Phage" (WtP), an easy-to-use and parallel multitool approach for phage prediction combined with an annotation and classification downstream strategy, thus supporting the user's decision-making process by summarizing the results of the different prediction tools in charts and tables. WtP is reproducible and scales to thousands of datasets through a workflow manager (Nextflow). WtP is freely available under a GPL-3.0 license (https://github.com/replikation/What_the_Phage).

Evaluation of tangential flow filtration coupled to long-read sequencing for ostreid herpesvirus type 1 genome assembly

Whole-genome sequencing is widely used to better understand the transmission dynamics, the evolution and the emergence of new variants of viral pathogens. This can bring crucial information to stakeholders for disease management. Unfortunately, aquatic virus genomes are usually difficult to characterize because most of these viruses cannot be easily propagated in vitro. Developing methodologies for routine genome sequencing of aquatic viruses is timely given the ongoing threat of disease emergence. This is particularly true for pathogenic viruses infecting species of commercial interest that are widely exchanged between production basins or countries. For example, the ostreid herpesvirus type 1 (OsHV-1) is a Herpesvirus widely associated with mass mortality events of juvenile Pacific oyster Crassostrea gigas. Genomes of Herpesviruses are large and complex with long direct and inverted terminal repeats. In addition, OsHV-1 is unculturable. It therefore accumulates several features that make its genome sequencing and assembly challenging. To overcome these difficulties, we developed a tangential flow filtration (TFF) method to enrich OsHV-1 infective particles from infected host tissues. This virus purification allowed us to extract high molecular weight and high-quality viral DNA that was subjected to Illumina short-read and Nanopore long-read sequencing. Dedicated bioinformatic pipelines were developed to assemble complete OsHV-1 genomes with reads from both sequencing technologies. Nanopore sequencing allowed characterization of new structural variations and major viral isomers while having 99,98 % of nucleotide identity with the Illumina assembled genome. Our study shows that TFF-based purification method, coupled with Nanopore sequencing, is a promising approach to enable in field sequencing of unculturable aquatic DNA virus.

Marine viral particles reveal an expansive repertoire of phage-parasitizing mobile elements

Phage satellites are mobile genetic elements that propagate by parasitizing bacteriophage replication. We report here the discovery of abundant and diverse phage satellites that were packaged as concatemeric repeats within naturally occurring bacteriophage particles in seawater. These same phage-parasitizing mobile elements were found integrated in the genomes of dominant co-occurring bacterioplankton species. Like known phage satellites, many marine phage satellites encoded genes for integration, DNA replication, phage interference, and capsid assembly. Many also contained distinctive gene suites indicative of unique virus hijacking, phage immunity, and mobilization mechanisms. Marine phage satellite sequences were widespread in local and global oceanic virioplankton populations, reflecting their ubiquity, abundance, and temporal persistence in marine planktonic communities worldwide. Their gene content and putative life cycles suggest they may impact host-cell phage immunity and defense, lateral gene transfer, bacteriophage-induced cell mortality and cellular host and virus productivity. Given that marine phage satellites cannot be distinguished from bona fide viral particles via commonly used microscopic techniques, their predicted numbers (∼3.2 × 1026 in the ocean) may influence current estimates of virus densities, production, and virus-induced mortality. In total, the data suggest that marine phage satellites have potential to significantly impact the ecology and evolution of bacteria and their viruses throughout the oceans. We predict that any habitat that harbors bacteriophage will also harbor similar phage satellites, making them a ubiquitous feature of most microbiomes on Earth.

Marine viruses and climate change: Virioplankton, the carbon cycle, and our future ocean

Interactions between marine viruses and microbes are a critical part of the oceanic carbon cycle. The impacts of virus-host interactions range from short-term disruptions in the mobility of microbial biomass carbon to higher trophic levels through cell lysis (i.e., the viral shunt) to long-term reallocation of microbial biomass carbon to the deep sea through accelerating the biological pump (i.e., the viral shuttle). The biogeochemical backdrop of the ocean-the physical, chemical, and biological landscape-influences the likelihood of both virus-host interactions and particle formation, and the fate and flow of carbon. As climate change reshapes the oceanic landscape through large-scale shifts in temperature, circulation, stratification, and acidification, virus-mediated carbon flux is likely to shift in response. Dynamics in the directionality and magnitude of changes in how, where, and when viruses mediate the recycling or storage of microbial biomass carbon is largely unknown. Integrating viral infection dynamics data obtained from experimental models and field systems, with particle motion microphysics and global observations of oceanic biogeochemistry, into improved ecosystem models will enable viral oceanographers to better predict the role of viruses in marine carbon cycling in the future ocean.

Long-Read Metagenomics Improves the Recovery of Viral Diversity from Complex Natural Marine Samples

The recovery of DNA from viromes is a major obstacle in the use of long-read sequencing to study their genomes. For this reason, the use of cellular metagenomes (>0.2-μm size range) emerges as an interesting complementary tool, since they contain large amounts of naturally amplified viral genomes from prelytic replication. We have applied second-generation (Illumina NextSeq; short reads) and third-generation (PacBio Sequel II; long reads) sequencing to compare the diversity and features of the viral community in a marine sample obtained from offshore waters of the western Mediterranean. We found that a major wedge of the expected marine viral diversity was directly recovered by the raw PacBio circular consensus sequencing (CCS) reads. More than 30,000 sequences were detected only in this data set, with no homologues in the long- and short-read assembly, and ca. 26,000 had no homologues in the large data set of the Global Ocean Virome 2 (GOV2), highlighting the information gap created by the assembly bias. At the level of complete viral genomes, the performance was similar in both approaches. However, the hybrid long- and short-read assembly provided the longest average length of the sequences and improved the host assignment. Although no novel major clades of viruses were found, there was an increase in the intraclade genomic diversity recovered by long reads that produced an enriched assessment of the real diversity and allowed the discovery of novel genes with biotechnological potential (e.g., endolysin genes).

IMPORTANCE We explored the vast genetic diversity of environmental viruses by using a combination of cellular metagenome (as opposed to virome) sequencing using high-fidelity long-read sequences (in this case, PacBio CCS). This approach resulted in the recovery of a representative sample of the viral population, and it performed better (more phage contigs, larger average contig size) than Illumina sequencing applied to the same sample. By this approach, the many biases of assembly are avoided, as the CCS reads recovers (typically around 5 kb) complete genes and even operons, resulting in a better discovery of the viral gene diversity based on viral marker proteins. Thus, biotechnologically promising genes, such as endolysin genes, can be very efficiently searched with this approach. In addition, hybrid assembly produces more complete and longer contigs, which is particularly important for studying little-known viral groups such as the nucleocytoplasmic large DNA viruses (NCLDV).

Diversity and origins of bacterial and archaeal viruses on sinking particles reaching the abyssal ocean

Sinking particles and particle-associated microbes influence global biogeochemistry through particulate matter export from the surface to the deep ocean. Despite ongoing studies of particle-associated microbes, viruses in these habitats remain largely unexplored. Whether, where, and which viruses might contribute to particle production and export remain open to investigation. In this study, we analyzed 857 virus population genomes associated with sinking particles collected over three years in sediment traps moored at 4000 m in the North Pacific Subtropical Gyre. Particle-associated viruses here were linked to cellular hosts through matches to bacterial and archaeal metagenome-assembled genome (MAG)-encoded prophages or CRISPR spacers, identifying novel viruses infecting presumptive deep-sea bacteria such as Colwellia, Moritella, and Shewanella. We also identified lytic viruses whose abundances correlated with particulate carbon flux and/or were exported from the photic to abyssal ocean, including cyanophages. Our data are consistent with some of the predicted outcomes of the viral shuttle hypothesis, and further suggest that viral lysis of both autotrophic and heterotrophic prokaryotes may play a role in carbon export. Our analyses revealed the diversity and origins of prevalent viruses found on deep-sea sinking particles and identified prospective viral groups for future investigation into processes that govern particle export in the open ocean.

Evolving Perspective on the Origin and Diversification of Cellular Life and the Virosphere

The tree of life (TOL) is a powerful framework to depict the evolutionary history of cellular organisms through time, from our microbial origins to the diversification of multicellular eukaryotes that shape the visible biosphere today. During the past decades, our perception of the TOL has fundamentally changed, in part, due to profound methodological advances, which allowed a more objective approach to studying organismal and viral diversity and led to the discovery of major new branches in the TOL as well as viral lineages. Phylogenetic and comparative genomics analyses of these data have, among others, revolutionized our understanding of the deep roots and diversity of microbial life, the origin of the eukaryotic cell, eukaryotic diversity, as well as the origin, and diversification of viruses. In this review, we provide an overview of some of the recent discoveries on the evolutionary history of cellular organisms and their viruses and discuss a variety of complementary techniques that we consider crucial for making further progress in our understanding of the TOL and its interconnection with the virosphere.

Newly identified HMO-2011-type phages reveal genomic diversity and biogeographic distributions of this marine viral group

Viruses play critical roles in influencing biogeochemical cycles and adjusting host mortality, population structure, physiology, and evolution in the ocean. Marine viral communities are composed of numerous genetically distinct subfamily/genus-level viral groups. Among currently identified viral groups, the HMO-2011-type group is known to be dominant and broadly distributed. However, only four HMO-2011-type cultivated representatives that infect marine SAR116 and Roseobacter strains have been reported to date, and the genetic diversity, potential hosts, and ecology of this group remain poorly elucidated. Here, we present the genomes of seven HMO-2011-type phages that were isolated using four Roseobacter strains and one SAR11 strain, as well as additional 207 HMO-2011-type metagenomic viral genomes (MVGs) identified from various marine viromes. Phylogenomic and shared-gene analyses revealed that the HMO-2011-type group is a subfamily-level group comprising at least 10 discernible genus-level subgroups. Moreover, >2000 HMO-2011-type DNA polymerase sequences were identified, and the DNA polymerase phylogeny also revealed that the HMO-2011-type group contains diverse subgroups and is globally distributed. Metagenomic read-mapping results further showed that most HMO-2011-type phages are prevalent in global oceans and display distinct geographic distributions, with the distribution of most HMO-2011-type phages being associated with temperature. Lastly, we found that members in subgroup IX, represented by pelagiphage HTVC033P, were among the most abundant HMO-2011-type phages, which implies that SAR11 bacteria are crucial hosts for this viral group. In summary, our findings substantially expand current knowledge regarding the phylogenetic diversity, evolution, and distribution of HMO-2011-type phages, highlighting HMO-2011-type phages as major ecological agents that can infect certain key bacterial groups.

Whole-genome sequencing and genetic characteristics of representative porcine reproductive and respiratory syndrome virus (PRRSV) isolates in Korea

Background

Porcine reproductive and respiratory syndrome virus (PRRSV) is a macrophage-tropic arterivirus with extremely high genetic and pathogenic heterogeneity that causes significant economic losses in the swine industry worldwide. PRRSV can be divided into two species [PRRSV1 (European) and PRRSV2 (North American)] and is usually diagnosed and genetically differentiated into several lineages based on the ORF5 gene, which constitutes only 5% of the whole genome. This study was conducted to achieve nonselective amplification and whole-genome sequencing (WGS) based on a simplified sequence-independent, single-primer amplification (SISPA) technique with next-generation sequencing (NGS), and to genetically characterize Korean PRRSV field isolates at the whole genome level.

Methods

The SISPA-NGS method coupled with a bioinformatics pipeline was utilized to retrieve full length PRRSV genomes of 19 representative Korean PRRSV strains by de novo assembly. Phylogenetic analysis, analysis of the insertion and deletion (INDEL) pattern of nonstructural protein 2 (NSP2), and recombination analysis were conducted.

Results

Nineteen complete PRRSV genomes were obtained with a high depth of coverage by the SISPA-NGS method. Korean PRRSV1 belonged to the Korean-specific subtype 1A and vaccine-related subtype 1C lineages, showing no evidence of recombination and divergent genetic heterogeneity with conserved NSP2 deletion patterns. Among Korean PRRSV2 isolates, modified live vaccine (MLV)-related lineage 5 viruses, lineage 1 viruses, and nation-specific Korean lineages (KOR A, B and C) could be identified. The NSP2 deletion pattern of the Korean lineages was consistent with that of the MN-184 strain (lineage 1), which indicates the common ancestor and independent evolution of Korean lineages. Multiple recombination signals were detected from Korean-lineage strains isolated in the 2010s, suggesting natural interlineage recombination between circulating KOR C and MLV strains. Interestingly, the Korean strain GGYC45 was identified as a recombinant KOR C and MLV strain harboring the KOR B ORF5 gene and might be the ancestor of currently circulating KOR B strains. Additionally, two novel lineage 1 recombinants of NADC30-like and NADC34-like viruses were detected.

Conclusion

Genome-wide analysis of Korean PRRSV isolates retrieved by the SISPA-NGS method and de novo assembly, revealed complex evolution and recombination in the field. Therefore, continuous surveillance of PRRSV at the whole genome level should be conducted, and new vaccine strategies for more efficient control of the virus are needed.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12985-022-01790-6.

Unexpected myriad of co-occurring viral strains and species in one of the most abundant and microdiverse viruses on Earth

Viral genetic microdiversity drives adaptation, pathogenicity, and speciation and has critical consequences for the viral-host arms race occurring at the strain and species levels, which ultimately impact microbial community structure and biogeochemical cycles. Despite the fact that most efforts have focused on viral macrodiversity, little is known about the microdiversity of ecologically important viruses on Earth. Recently, single-virus genomics discovered the putatively most abundant ocean virus in temperate and tropical waters: the uncultured dsDNA virus vSAG 37-F6 infecting Pelagibacter, the most abundant marine bacteria. In this study, we report the cooccurrence of up to ≈1,500 different viral strains (>95% nucleotide identity) and ≈30 related species (80-95% nucleotide identity) in a single oceanic sample. Viral microdiversity was maintained over space and time, and most alleles were the result of synonymous mutations without any apparent adaptive benefits to cope with host translation codon bias and efficiency. Gene flow analysis used to delimitate species according to the biological species concept (BSC) revealed the impact of recombination in shaping vSAG 37-F6 virus and Pelagibacter speciation. Data demonstrated that this large viral microdiversity somehow mirrors the host species diversity since ≈50% of the 926 analyzed Pelagibacter genomes were found to belong to independent BSC species that do not significantly engage in gene flow with one another. The host range of this evolutionarily successful virus revealed that a single viral species can infect multiple Pelagibacter BSC species, indicating that this virus crosses not only formal BSC barriers but also biomes since viral ancestors are found in freshwater.

RNA Viruses in Aquatic Ecosystems through the Lens of Ecological Genomics and Transcriptomics

Massive amounts of data from nucleic acid sequencing have changed our perspective about diversity and dynamics of marine viral communities. Here, we summarize recent metatranscriptomic and metaviromic studies targeting predominantly RNA viral communities. The analysis of RNA viromes reaffirms the abundance of lytic (+) ssRNA viruses of the order Picornavirales, but also reveals other (+) ssRNA viruses, including RNA bacteriophages, as important constituents of extracellular RNA viral communities. Sequencing of dsRNA suggests unknown diversity of dsRNA viruses. Environmental metatranscriptomes capture the dynamics of ssDNA, dsDNA, ssRNA, and dsRNA viruses simultaneously, unravelling the full complexity of viral dynamics in the marine environment. RNA viruses are prevalent in large size fractions of environmental metatranscriptomes, actively infect marine unicellular eukaryotes larger than 3 µm, and can outnumber bacteriophages during phytoplankton blooms. DNA and RNA viruses change abundance on hourly timescales, implying viral control on a daily temporal basis. Metatranscriptomes of cultured protists host a diverse community of ssRNA and dsRNA viruses, often with multipartite genomes and possibly persistent intracellular lifestyles. We posit that RNA viral communities might be more diverse and complex than formerly anticipated and that the influence they exert on community composition and global carbon flows in aquatic ecosystems may be underestimated.

viralFlye: assembling viruses and identifying their hosts from long-read metagenomics data

Although the use of long-read sequencing improves the contiguity of assembled viral genomes compared to short-read methods, assembling complex viral communities remains an open problem. We describe the viralFlye tool for identification and analysis of metagenome-assembled viruses in long-read assemblies. We show it significantly improves viral assemblies and demonstrate that long-reads result in a much larger array of predicted virus-host associations as compared to short-read assemblies. We demonstrate that the identification of novel CRISPR arrays in bacterial genomes from a newly assembled metagenomic sample provides information for predicting novel hosts for novel viruses.

Whole Genome Assembly of Human Papillomavirus by Nanopore Long-Read Sequencing

Human papillomavirus (HPV) is a causal agent for most cervical cancers. The physical status of the HPV genome in these cancers could be episomal, integrated, or both. HPV integration could serve as a biomarker for clinical diagnosis, treatment, and prognosis. Although whole-genome sequencing by next-generation sequencing (NGS) technologies, such as the Illumina sequencing platform, have been used for detecting integrated HPV genome in cervical cancer, it faces challenges of analyzing long repeats and translocated sequences. In contrast, Oxford nanopore sequencing technology can generate ultra-long reads, which could be a very useful tool for determining HPV genome sequence and its physical status in cervical cancer. As a proof of concept, in this study, we completed whole genome sequencing from a cervical cancer tissue and a CaSki cell line with Oxford Nanopore Technologies. From the cervical cancer tissue, a 7,894 bp-long HPV35 genomic sequence was assembled from 678 reads at 97-fold coverage of HPV genome, sharing 99.96% identity with the HPV sequence obtained by Sanger sequencing. A 7904 bp-long HPV16 genomic sequence was assembled from data generated from the CaSki cell line at 3857-fold coverage, sharing 99.99% identity with the reference genome (NCBI: U89348). Intriguingly, long reads generated by nanopore sequencing directly revealed chimeric cellular-viral sequences and concatemeric genomic sequences, leading to the discovery of 448 unique integration breakpoints in the CaSki cell line and 60 breakpoints in the cervical cancer sample. Taken together, nanopore sequencing is a unique tool to identify HPV sequences and would shed light on the physical status of HPV genome in its associated cancers.

Direct RNA nanopore sequencing of Pseudomonas aeruginosa clone C transcriptomes

The transcriptomes of Pseudomonas aeruginosa clone C isolates NN2 and SG17M during the mid-exponential and early stationary phase of planktonic growth were evaluated by direct RNA sequencing on the nanopore platform and compared with established short-read cDNA sequencing on the Illumina platform. Fifty to ninety percent of the sense RNAs turned out to be rRNA molecules followed by similar proportions of mRNA transcripts and non-coding RNAs. Both platforms detected similar proportions of uncharged tRNAs and 29 yet undescribed antisense tRNAs. For example, the rarest arginine codon was paired with the most abundant tRNA^Arg, and the tRNA^Arg gene is missing for the most frequent arginine codon. More than 90% of the antisense RNA molecules were complementary to a coding sequence. The antisense RNAs were evenly distributed in the genomes. Direct RNA sequencing identified more than 4,000 distinct non-overlapping antisense RNAs during exponential and stationary growth. Besides highly expressed small antisense RNAs less than 200 bases in size, a population of longer antisense RNAs was sequenced that covered a broad range of a few hundred to thousands of bases and could be complementary to a contig of several genes. In summary, direct RNA sequencing identified yet undescribed RNA molecules and an unexpected composition of the pools of tRNAs, sense and antisense RNAs. IMPORTANCE Genome-wide gene expression of bacteria is commonly studied by high-throughput sequencing of size-selected cDNA fragment libraries of reverse-transcribed RNA preparations. However, the depletion of ribosomal RNAs, enzymatic reverse transcription and the fragmentation, size selection and amplification during library preparation lead to inevitable losses of information about the initial composition of the RNA pool. We demonstrate that direct RNA sequencing on the nanopore platform can overcome these limitations. Nanopore sequencing of total RNA yielded novel insights into the Pseudomonas aeruginosa transcriptome that - if replicated in other species - will change our view of the bacterial RNA world. The discovery of sense - antisense pairs of tmRNA, tRNAs and mRNAs indicates a further and unknown level of gene regulation in bacteria.

Composition and Dynamics of H1N1 and H7N9 Influenza A Virus Quasispecies in a Co-infected Patient Analyzed by Single Molecule Sequencing Technology

A human co-infected with H1N1 and H7N9 subtypes influenza A virus (IAV) causes a complex infectious disease. The identification of molecular-level variations in composition and dynamics of IAV quasispecies will help to understand the pathogenesis and provide guidance for precision medicine treatment. In this study, using single-molecule real-time sequencing (SMRT) technology, we successfully acquired full-length IAV genomic sequences and quantified their genotypes abundance in serial samples from an 81-year-old male co-infected with H1N1 and H7N9 subtypes IAV. A total of 26 high diversity nucleotide loci was detected, in which the A-G base transversion was the most abundant substitution type (67 and 64%, in H1N1 and H7N9, respectively). Seven significant amino acid variations were detected, such as NA:H275Y and HA: R222K in H1N1 as well as PB2:E627K and NA: K432E in H7N9, which are related to viral drug-resistance or mammalian adaptation. Furtherly, we retrieved 25 H1N1 and 22 H7N9 genomic segment haplotypes from the eight samples based on combining high-diversity nucleotide loci, which provided a more concise overview of viral quasispecies composition and dynamics. Our approach promotes the popularization of viral quasispecies analysis in a complex infectious disease, which will boost the understanding of viral infections, pathogenesis, evolution, and precision medicine.