Towards quantitative viromics for both double-stranded and single-stranded DNA viruses

Identification and characterization of a novel circovirus in Iberian lynx in Spain

Circoviruses cause severe disease in pigs and birds. Canine circovirus has thus far only been associated with respiratory and gastrointestinal disorders and systemic disease in dogs. The Iberian lynx (Lynx pardinus) is one of the most endangered carnivores in Europe and the most endangered felid worldwide. Exploring the virome of these animals may be important in terms of virus discovery and assessing the interspecies-circulation of viruses from related carnivores. In this study, 162 spleen samples from Iberian lynx were screened for CRESS DNA viruses. Overall, 11 (6.8%) of 162 samples tested positive using a consensus PCR. Partial rep sequences were tightly related to each other (96.6-100%). Specific molecular protocols were designed on the partial rep sequences of the novel virus, Iberian lynx-associated circovirus-1 (ILCV-1). By screening a subset of 45 spleen samples, the infection rate of ILCV-1 in Iberian lynxes was 57.8% (26/45). ILCV-1 strains formed a separate cluster intermingled with bat, rodent, mongoose, and felid circoviruses. The genome of the novel virus displayed the highest nucleotide identity (64.3-65.3%) to mongoose circoviruses, thus representing a novel candidate circovirus species. The detection of these viruses in the spleen tissues could suggest systemic infection in the animal host. Overall, these findings suggest that this novel circovirus is common in the Iberian lynx. Further studies are warranted to assess the possible health implications of ILCV-1 in this endangered species.

CAIM: Coverage-based Analysis for Identification of Microbiome

Accurate taxonomic profiling of microbial taxa in a metagenomic sample is vital to gain insights into microbial ecology. Recent advancements in sequencing technologies have contributed tremendously toward understanding these microbes at species resolution through a whole shotgun metagenomic (WMS) approach. In this study, we developed a new bioinformatics tool, CAIM, for accurate taxonomic classification and quantification within both long- and short-read metagenomic samples using an alignment-based method. CAIM depends on two different containment techniques to identify species in metagenomic samples using their genome coverage information to filter out false positives rather than the traditional approach of relative abundance. In addition, we propose a nucleotide-count based abundance estimation, which yield lesser root mean square error than the traditional read-count approach. We evaluated the performance of CAIM on 28 metagenomic mock communities and 2 synthetic datasets by comparing it with other top-performing tools. CAIM maintained a consitently good performance across datasets in identifying microbial taxa and in estimating relative abundances than other tools. CAIM was then applied to a real dataset sequenced on both Nanopore (with and without amplification) and Illumina sequencing platforms and found high similality of taxonomic profiles between the sequencing platforms. Lastly, CAIM was applied to fecal shotgun metagenomic datasets of 232 colorectal cancer patients and 229 controls obtained from 4 different countries and primary 44 liver cancer patients and 76 controls. The predictive performance of models using the genome-coverage cutoff was better than those using the relative-abundance cutoffs in discriminating colorectal cancer and primary liver cancer patients from healthy controls with a highly confident species markers.

Identification and classification of the genomes of novel microviruses in poultry slaughterhouse

Microviridae is a family of phages with circular ssDNA genomes and they are widely found in various environments and organisms. In this study, virome techniques were employed to explore potential members of Microviridae in a poultry slaughterhouse, leading to the identification of 98 novel and complete microvirus genomes. Using a similarity clustering network classification approach, these viruses were found to belong to at least 6 new subfamilies within Microviridae and 3 higher-level taxonomic units. Genome size, GC content and genome structure of these new taxa showed evident regularities, validating the rationality of our classification method. Our method can divide microviruses into about 45 additional detailed clusters, which may serve as a new standard for classifying Microviridae members. Furthermore, by addressing the scarcity of host information for microviruses, the current study significantly broadened their host range and discovered over 20 possible new hosts, including important pathogenic bacteria such as Helicobacter pylori and Vibrio cholerae, as well as different taxa demonstrated different host specificities. The findings of this study effectively expand the diversity of the Microviridae family, providing new insights for their classification and identification. Additionally, it offers a novel perspective for monitoring and controlling pathogenic microorganisms in poultry slaughterhouse environments.

Comparison of Three Viral Nucleic Acid Preamplification Pipelines for Sewage Viral Metagenomics

Viral metagenomics is a useful tool for detecting multiple human viruses in urban sewage. However, more refined protocols are required for its effective use in disease surveillance. In this study, we investigated the performance of three different preamplification pipelines (specific to RNA viruses, DNA viruses or both) for viral genome sequencing using spiked-in Phosphate Buffered Saline and sewage samples containing known concentrations of viruses. We found that compared to the pipeline targeting all genome types, the RNA pipeline performed better in detecting RNA viruses in both spiked and unspiked sewage samples, allowing the detection of various mammalian viruses including members from the Reoviridae, Picornaviridae, Astroviridae and Caliciviridae. However, the DNA-specific pipeline did not improve the detection of mammalian DNA viruses. We also measured viral recovery by quantitative reverse transcription polymerase chain reaction and assessed the impact of genetic background (non-viral genetic material) on viral coverage. Our results indicate that viral recoveries were generally lower in sewage (average of 11.0%) and higher in Phosphate Buffered Saline (average of 23.4%) for most viruses. Additionally, spiked-in viruses showed lower genome coverage in sewage, demonstrating the negative effect of genetic background on sequencing. Finally, correlation analysis revealed a relationship between virus concentration and genome normalized reads per million, indicating that viral metagenomic sequencing can be semiquantitative.

Breaking the Ice: A Review of Phages in Polar Ecosystems

Bacteriophages, or phages, are viruses that infect and replicate within bacterial hosts, playing a significant role in regulating microbial populations and ecosystem dynamics. However, phages from extreme environments such as polar regions remain relatively understudied due to challenges such as restricted ecosystem access and low biomass. Understanding the diversity, structure, and functions of polar phages is crucial for advancing our knowledge of the microbial ecology and biogeochemistry of these environments. In this review, we will explore the current state of knowledge on phages from the Arctic and Antarctic, focusing on insights gained from -omic studies, phage isolation, and virus-like particle abundance data. Metagenomic studies of polar environments have revealed a high diversity of phages with unique genetic characteristics, providing insights into their evolutionary and ecological roles. Phage isolation studies have identified novel phage-host interactions and contributed to the discovery of new phage species. Virus-like particle abundance and lysis rate data, on the other hand, have highlighted the importance of phages in regulating bacterial populations and nutrient cycling in polar environments. Overall, this review aims to provide a comprehensive overview of the current state of knowledge about polar phages, and by synthesizing these different sources of information, we can better understand the diversity, dynamics, and functions of polar phages in the context of ongoing climate change, which will help to predict how polar ecosystems and residing phages may respond to future environmental perturbations.

Phylogenomic analysis expands the known repertoire of single-stranded DNA viruses in benthic zones of the South Indian Ocean

Single-stranded (ss) DNA viruses are ubiquitous and constitute some of the most diverse entities on Earth. Most studies have focused on ssDNA viruses from terrestrial environments resulting in a significant deficit in benthic ecosystems including aphotic zones of the South Indian Ocean (SIO). Here, we assess the diversity and phylogeny of ssDNA in deep waters of the SIO using a combination of established viral taxonomy tools and a Hidden Markov Model based approach. Replication initiator protein-associated (Rep) phylogenetic reconstruction and sequence similarity networks were used to show that the SIO hosts divergent and as yet unknown circular Rep-encoding ssDNA viruses. Several sequences appear to represent entirely novel families, expanding the repertoire of known ssDNA viruses. Results suggest that a small proportion of these viruses may be circular genetic elements, which may strongly influence the diversity of both eukaryotes and prokaryotes in the SIO. Taken together, our data show that the SIO harbours a diverse assortment of previously unknown ssDNA viruses. Due to their potential to infect a variety of hosts, these viruses may be crucial for marine nutrient recycling through their influence of the biological carbon pump.

With a pinch of salt: metagenomic insights into Namib Desert salt pan microbial mats and halites reveal functionally adapted and competitive communities

IMPORTANCE

The hyperarid Namib Desert is one of the oldest deserts on Earth. It contains multiple clusters of playas which are saline-rich springs surrounded by halite evaporites. Playas are of great ecological importance, and their indigenous (poly)extremophilic microorganisms are potentially involved in the precipitation of minerals such as carbonates and sulfates and have been of great biotechnological importance. While there has been a considerable amount of microbial ecology research performed on various Namib Desert edaphic microbiomes, little is known about the microbial communities inhabiting its multiple playas. In this work, we provide a comprehensive taxonomic and functional potential characterization of the microbial, including viral, communities of sediment mats and halites from two distant salt pans of the Namib Desert, contributing toward a better understanding of the ecology of this biome.

Viruses under the Antarctic Ice Shelf are active and potentially involved in global nutrient cycles

Viruses play an important role in the marine ecosystem. However, our comprehension of viruses inhabiting the dark ocean, and in particular, under the Antarctic Ice Shelves, remains limited. Here, we mine single-cell genomic, transcriptomic, and metagenomic data to uncover the viral diversity, biogeography, activity, and their role as metabolic facilitators of microbes beneath the Ross Ice Shelf. This is the largest Antarctic ice shelf with a major impact on global carbon cycle. The viral community found in the cavity under the ice shelf mainly comprises endemic viruses adapted to polar and mesopelagic environments. The low abundance of genes related to lysogenic lifestyle (<3%) does not support a predominance of the Piggyback-the-Winner hypothesis, consistent with a low-productivity habitat. Our results indicate a viral community actively infecting key ammonium and sulfur-oxidizing chemolithoautotrophs (e.g. Nitrosopumilus spp, Thioglobus spp.), supporting a "kill-the-winner" dynamic. Based on genome analysis, these viruses carry specific auxiliary metabolic genes potentially involved in nitrogen, sulfur, and phosphorus acquisition. Altogether, the viruses under Antarctic ice shelves are putatively involved in programming the metabolism of ecologically relevant microbes that maintain primary production in these chemosynthetically-driven ecosystems, which have a major role in global nutrient cycles.

Novel RNA viruses from the Atlantic Ocean: Ecogenomics, biogeography, and total virioplankton mass contribution from surface to the deep ocean

Marine viruses play a major role in the energy and nutrient cycle and affect the evolution of their hosts. Despite their importance, there is still little knowledge about RNA viruses. Here, we have explored the Atlantic Ocean, from surface to deep (4.296 m), and used viromics and quantitative methods to unveil the genomics, biogeography, and the mass contribution of RNA viruses to the total viroplankton. A total of 2481 putative RNA viral contigs (>500 bp) and 107 larger bona fide RNA viral genomes (>2.5 kb) were identified; 88 of them representing novel viruses belonging mostly to two clades: Yangshan assemblage (sister clade to the class Alsuviricetes) and Nodaviridae. These viruses were highly endemic and locally abundant, with little or no presence in other oceans since only ≈15% of them were found in at least one of the Tara sampling metatranscriptomes. Quantitative data indicated that the abundance of RNA viruses in the surface and deep chlorophyll maximum zone was within ≈106  VLP/mL representing a potential contribution of 5.2%-24.4% to the total viroplankton community (DNA and RNA viruses), with DNA viruses being the predominant members (≈107  VLP/mL). However, for the deep sample, the observed trend was the opposite, although as further discussed, several biases should be considered. Together these results contribute to our understanding of the diversity, abundance, and distribution of RNA viruses in the oceans and provide a basis for further investigation into their ecological roles and biogeography.

Benchmarking of virome metagenomic analysis approaches using a large, 60+ members, viral synthetic community

IMPORTANCE

We report here efforts to benchmark performance of two widespread approaches for virome analysis, which target either virion-associated nucleic acids (VANA) or highly purified double-stranded RNAs (dsRNAs). This was achieved using synthetic communities of varying complexity levels, up to a highly complex community of 72 viral agents (115 viral molecules) comprising isolates from 21 families and 61 genera of plant viruses. The results obtained confirm that the dsRNA-based approach provides a more complete representation of the RNA virome, in particular, for high complexity ones. However, for viromes of low to medium complexity, VANA appears a reasonable alternative and would be the preferred choice if analysis of DNA viruses is of importance. Several parameters impacting performance were identified as well as a direct relationship between the completeness of virome description and sample sequencing depth. The strategy, results, and tools used here should prove useful in a range of virome analysis efforts.

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.

Decoding the DNA and RNA viromes of a tropical urban lagoon

Human and livestock sewage is one of the major causes of excess nutrients, leading to the eutrophication of aquatic ecosystems and potentially to the emergence or spread of pathogenic viruses. This study aimed to investigate the composition and diversity of aquatic viromes in a highly anthropized lagoon, to identify the presence of pathogenic taxa and to explore their use as possible viral indicators of faecal contamination. For this, water and sediment samples were collected in the Ebrié Lagoon (Ivory Coast) at seven stations with contrasting levels of eutrophication. The DNA viromes of the planktonic and the benthic compartments were highly divergent, but were not influenced by the level of eutrophication. Conversely, the RNA viromes in the water column were comparable to those found in sediment, but showed significant differences between the stations. We detected the presence of viral DNA and RNA sequences we had assigned as indicators of faecal contamination (smacovirus, pecovirus and pepper mild mottle virus) as well as human pathogens (human cyclovirus, coxsackie B virus and picobirnavirus), which were all enriched in the most eutrophicated sites. These findings suggest that the examination of viromes represents a promising tool for assessing the state of human-induced contamination of aquatic ecosystems.

Microviruses: A World Beyond phiX174

Two decades of metagenomic analyses have revealed that in many environments, small (∼5 kb), single-stranded DNA phages of the family Microviridae dominate the virome. Although the emblematic microvirus phiX174 is ubiquitous in the laboratory, most other microviruses, particularly those of the gokushovirus and amoyvirus lineages, have proven to be much more elusive. This puzzling lack of representative isolates has hindered insights into microviral biology. Furthermore, the idiosyncratic size and nature of their genomes have resulted in considerable misjudgments of their actual abundance in nature. Fortunately, recent successes in microvirus isolation and improved metagenomic methodologies can now provide us with more accurate appraisals of their abundance, their hosts, and their interactions. The emerging picture is that phiX174 and its relatives are rather rare and atypical microviruses, and that a tremendous diversity of other microviruses is ready for exploration.

The impact of storage buffer and storage conditions on fecal samples for bacteriophage infectivity and metavirome analyses

BACKGROUND

There is an increasing interest in investigating the human gut virome for its influence on the gut bacterial community and its putative influence on the trajectory towards health or disease. Most gut virome studies are based on sequencing of stored fecal samples. However, relatively little is known about how conventional storage buffers and storage conditions affect the infectivity of bacteriophages and influence the downstream metavirome sequencing.

RESULTS

We demonstrate that the infectivity and genome recovery rate of different spiked bacteriophages (T4, c2 and Phi X174) are variable and highly dependent on storage buffers. Regardless of the storage temperature and timespan, all tested phages immediately lost 100% (DNA/RNA Shield) or more than 90% (StayRNA and RNAlater) of their infectivity. Generally, in SM buffer at 4 °C phage infectivity was preserved for up to 30 days and phage DNA integrity was maintained for up to 100 days. While in CANVAX, the most effective buffer, all spiked phage genomes were preserved for at least 100 days. Prolonged storage time (500 days) at - 80 °C impacted viral diversity differently in the different buffers. Samples stored in CANVAX or DNA/RNA Shield buffer had the least shifts in metavirome composition, after prolonged storage, but they yielded more contigs classified as "uncharacterised". Moreover, in contrast to the SM buffer, these storage buffers yielded a higher fraction of bacterial DNA in metavirome-sequencing libraries. We demonstrated that the latter was due to inactivation of the DNases employed to remove extra-cellular DNA during virome extraction. The latter could be partly avoided by employing additional washing steps prior to virome extraction.

CONCLUSION

Fecal sample storage buffers and storage conditions (time and temperature) strongly influence bacteriophage infectivity and viral composition as determined by plaque assay and metavirome sequencing. The choice of buffer had a larger effect than storage temperature and storage time on the quality of the viral sequences and analyses. Based on these results, we recommend storage of fecal virome samples at in SM buffer at 4 °C for the isolation of viruses and at - 80 °C for metagenomic applications if practically feasible (i.e., access to cold storage). For fecal samples stored in other buffers, samples should be cleared of these buffers before viral extraction and sequencing. Video Abstract.

Investigating the genetic diversity of CRESS DNA viruses in cats identifies a novel feline circovirus and unveils exposure of cats to canine circovirus

Circular replication-associated protein (Rep)-encoding single stranded (CRESS) DNA viruses include Circoviruses which have been found in several animal species and in human specimens. Circoviruses are associated with severe disease in pigs and birds and with respiratory and gastrointestinal disorders and systemic disease in dogs. In cats there are only a few anecdotical studies reporting CRESS DNA viruses. In this study, a total of 530 samples (361 sera, 131 stools, and 38 respiratory swabs) from cats, were screened for the presence of CRESS DNA viruses. Overall, 48 (9.0%) of 530 samples tested positive using a pan-Rep PCR. A total of 30 Rep sequences were obtained. Ten sequences of fecal origin were tightly related to each other (82.4-100% nt identity) and more distantly related to mongoose circoviruses (68.3 to 77.2% nt identity). At genome level these circoviruses displayed the highest nt identity (74.3-78.7%) to mongoose circoviruses thus representing a novel circovirus species. Circoviruses from different animal hosts (n = 12) and from humans (n = 8) were also identified. However, six Rep sequences were obtained from serum samples, including canine circoviruses, a human cyclovirus and human and fish-associated CRESS DNA viruses. The presence of these viruses in the sera would imply, to various extent, virus replication in the animal host, able to sustain viremia. Overall, these findings indicate a wide genetic diversity of CRESS DNA viruses in cats and warrant further investigations.

Chromosome folding and prophage activation reveal specific genomic architecture for intestinal bacteria

BACKGROUND

Bacteria and their viruses, bacteriophages, are the most abundant entities of the gut microbiota, a complex community of microorganisms associated with human health and disease. In this ecosystem, the interactions between these two key components are still largely unknown. In particular, the impact of the gut environment on bacteria and their associated prophages is yet to be deciphered.

RESULTS

To gain insight into the activity of lysogenic bacteriophages within the context of their host genomes, we performed proximity ligation-based sequencing (Hi-C) in both in vitro and in vivo conditions on the 12 bacterial strains of the OMM¹² synthetic bacterial community stably associated within mice gut (gnotobiotic mouse line OMM¹²). High-resolution contact maps of the chromosome 3D organization of the bacterial genomes revealed a wide diversity of architectures, differences between environments, and an overall stability over time in the gut of mice. The DNA contacts pointed at 3D signatures of prophages leading to 16 of them being predicted as functional. We also identified circularization signals and observed different 3D patterns between in vitro and in vivo conditions. Concurrent virome analysis showed that 11 of these prophages produced viral particles and that OMM¹² mice do not carry other intestinal viruses.

CONCLUSIONS

The precise identification by Hi-C of functional and active prophages within bacterial communities will unlock the study of interactions between bacteriophages and bacteria across conditions (healthy vs disease). Video Abstract.

Expanding known viral diversity in the healthy infant gut

The gut microbiome is shaped through infancy and impacts the maturation of the immune system, thus protecting against chronic disease later in life. Phages, or viruses that infect bacteria, modulate bacterial growth by lysis and lysogeny, with the latter being especially prominent in the infant gut. Viral metagenomes (viromes) are difficult to analyse because they span uncharted viral diversity, lacking marker genes and standardized detection methods. Here we systematically resolved the viral diversity in faecal viromes from 647 1-year-olds belonging to Copenhagen Prospective Studies on Asthma in Childhood 2010, an unselected Danish cohort of healthy mother-child pairs. By assembly and curation we uncovered 10,000 viral species from 248 virus family-level clades (VFCs). Most (232 VFCs) were previously unknown, belonging to the Caudoviricetes viral class. Hosts were determined for 79% of phage using clustered regularly interspaced short palindromic repeat spacers within bacterial metagenomes from the same children. Typical Bacteroides-infecting crAssphages were outnumbered by undescribed phage families infecting Clostridiales and Bifidobacterium. Phage lifestyles were conserved at the viral family level, with 33 virulent and 118 temperate phage families. Virulent phages were more abundant, while temperate ones were more prevalent and diverse. Together, the viral families found in this study expand existing phage taxonomy and provide a resource aiding future infant gut virome research.

Capturing variation in metagenomic assembly graphs with MetaCortex

MOTIVATION

The assembly of contiguous sequence from metagenomic samples presents a particular challenge, due to the presence of multiple species, often closely related, at varying levels of abundance. Capturing diversity within species, for example, viral haplotypes, or bacterial strain-level diversity, is even more challenging.

RESULTS

We present MetaCortex, a metagenome assembler that captures intra-species diversity by searching for signatures of local variation along assembled sequences in the underlying assembly graph and outputting these sequences in sequence graph format. We show that MetaCortex produces accurate assemblies with higher genome coverage and contiguity than other popular metagenomic assemblers on mock viral communities with high levels of strain-level diversity and on simulated communities containing simulated strains.

AVAILABILITY AND IMPLEMENTATION

Source code is freely available to download from https://github.com/SR-Martin/metacortex, is implemented in C and supported on MacOS and Linux. The version used for the results presented in this article is available at doi.org/10.5281/zenodo.7273627.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

KMCP: accurate metagenomic profiling of both prokaryotic and viral populations by pseudo-mapping

MOTIVATION

The growing number of microbial reference genomes enables the improvement of metagenomic profiling accuracy but also imposes greater requirements on the indexing efficiency, database size and runtime of taxonomic profilers. Additionally, most profilers focus mainly on bacterial, archaeal and fungal populations, while less attention is paid to viral communities.

RESULTS

We present KMCP (K-mer-based Metagenomic Classification and Profiling), a novel k-mer-based metagenomic profiling tool that utilizes genome coverage information by splitting the reference genomes into chunks and stores k-mers in a modified and optimized Compact Bit-Sliced Signature Index for fast alignment-free sequence searching. KMCP combines k-mer similarity and genome coverage information to reduce the false positive rate of k-mer-based taxonomic classification and profiling methods. Benchmarking results based on simulated and real data demonstrate that KMCP, despite a longer running time than all other methods, not only allows the accurate taxonomic profiling of prokaryotic and viral populations but also provides more confident pathogen detection in clinical samples of low depth.

AVAILABILITY AND IMPLEMENTATION

The software is open-source under the MIT license and available at https://github.com/shenwei356/kmcp.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Impact of antibiotic perturbation on fecal viral communities in mice

Viruses and bacteriophages have a strong impact on intestinal barrier function and the composition and functional properties of commensal bacterial communities. Shifts of the fecal virome might be involved in human diseases, including inflammatory bowel disease (IBD). Loss-of-function variants in the nucleotide-binding oligomerization domain-containing protein 2 (NOD2) gene are associated with an increased risk of developing Crohn's disease, a subtype of human chronic IBD, where specific changes in fecal viral communities have also been described. To improve our understanding of the dynamics of the enteric virome, we longitudinally characterized the virome in fecal samples from wild-type C57BL/6J and NOD2 knock-out mice in response to an antibiotic perturbation. Sequencing of virus-like particles demonstrated both a high diversity and high interindividual variation of the murine fecal virome composed of eukaryotic viruses and bacteriophages. Antibiotics had a significant impact on the fecal murine virome. Viral community composition only partially recovered in the observation period (10 weeks after cessation of antibiotics) irrespective of genotype. However, compositional shifts in the virome and bacteriome were highly correlated, suggesting that the loss of specific phages may contribute to prolonged dysregulation of the bacterial community composition. We suggest that therapeutic interference with the fecal virome may represent a novel approach in microbiota-targeted therapies.

Characterization of the human gut virome in metabolic and autoimmune diseases

The intestinal microbiome is dominated by bacteria and plays a pivotal role in the occurrence and development of disease, including several metabolic and autoimmune disorders. While intestinal viral communities, primarily made up of bacteriophages, are also thought to play a role in disease pathogenesis in the gastrointestinal tract, they have received much less attention than intestinal bacteria. Thus, there is limited information about the relationship between bacteriophages and disease. This review explores a potential role for the intestinal viral microbiome in various metabolic and autoimmune diseases.

Female reproduction and viral infection in a long-lived mammal

For energetically limited organisms, life-history theory predicts trade-offs between reproductive effort and somatic maintenance. This is especially true of female mammals, for whom reproduction presents multifarious energetic and physiological demands. Here, we examine longitudinal changes in the gut virome (viral community) with respect to reproductive status in wild mature female chimpanzees Pan troglodytes schweinfurthii from two communities, Kanyawara and Ngogo, in Kibale National Park, Uganda. We used metagenomic methods to characterize viromes of individual chimpanzees while they were cycling, pregnant and lactating. Females from Kanyawara, whose territory abuts the park's boundary, had higher viral richness and loads (relative quantity of viral sequences) than females from Ngogo, whose territory is more energetically rich and located farther from large human settlements. Viral richness (total number of distinct viruses per sample) was higher when females were lactating than when cycling or pregnant. In pregnant females, viral richness increased with estimated day of gestation. Richness did not vary with age, in contrast to prior research showing increased viral abundance in older males from these same communities. Our results provide evidence of short-term physiological trade-offs between reproduction and infection, which are often hypothesized to constrain health in long-lived species.

The Viral Fraction Metatranscriptomes of Lake Baikal

This article characterises viral fraction metatranscriptomes (smaller than 0.2 µm) from the pelagic zone of oligotrophic Lake Baikal (Russia). The study revealed the dominance of transcripts of DNA viruses: bacteriophages and algal viruses. We identified transcripts similar to Pithovirus sibericum, a nucleocytoplasmic large DNA virus (NCLDV) isolated from the permafrost region of Eastern Siberia. Among the families detected were RNA viruses assigned to Retroviridae, Metaviridae, Potyviridae, Astroviridae, and Closteroviridae. Using the PHROG, SEED subsystems databases, and the VOGDB, we indicated that the bulk of transcripts belong to the functional replication of viruses. In a comparative unweighted pair group method with arithmetic mean (UPGMA) analysis, the transcripts from Lake Baikal formed a separate cluster included in the clade with transcripts from other freshwater lakes, as well as marine and oceanic waters, while there was no separation based on the trophic state of the water bodies, the size of the plankton fraction, or salinity.

Novel lineages of single-stranded DNA phages that coevolved with the symbiotic bacteria Rhizobium

This study describes novel single-stranded DNA phages isolated from common bean agriculture soils by infection of the nitrogen-fixing symbiotic bacteria Rhizobium etli and R. phaseoli. A total of 29 phages analyzed have 4.3–6 kb genomes in size and GC 59–60%. They belong to different clades unrelated to other Microviridae subfamilies. Three-dimensional models of the major capsid protein (MCP) showed a conserved β-barrel structural “jelly-roll” fold. A variable-length loop in the MCPs distinguished three Rhizobium microvirus groups. Microviridae subfamilies were consistent with viral clusters determined by the protein-sharing network. All viral clusters, except for Bullavirinae, included mostly microviruses identified in metagenomes from distinct ecosystems. Two Rhizobium microvirus clusters, chaparroviruses, and chicoviruses, were included within large viral unknown clusters with microvirus genomes identified in diverse metagenomes. A third Rhizobium microvirus cluster belonged to the subfamily Amoyvirinae. Phylogenetic analysis of the MCP confirms the divergence of the Rhizobium microviruses into separate clades. The phylogeny of the bacterial hosts matches the microvirus MCP phylogeny, suggesting a coevolutionary history between the phages and their bacterial host. This study provided essential biological information on cultivated microvirus for understanding the evolution and ecological diversification of the Microviridae family in diverse microbial ecosystems.

Optimization and evaluation of viral metagenomic amplification and sequencing procedures toward a genome-level resolution of the human fecal DNA virome

INTRODUCTION

Viruses in the human gut have been linked to health and disease. Deciphering the gut virome is dependent on metagenomic sequencing of the virus-like particles (VLPs) purified from the fecal specimens. A major limitation of conventional viral metagenomic sequencing is the low recoverability of viral genomes from the metagenomic dataset.

OBJECTIVES

To develop an optimal method for viral amplification and metagenomic sequencing for maximizing the recovery of viral genomes.

METHODS

We performed parallel virus enrichment and DNA extraction to generate ∼ 30 viral DNA samples from each of 5 fresh fecal specimens and conducted the experiments including 1) optimizing the cycle number for high-fidelity enzyme-based PCR amplification, 2) evaluating the reproducibility of the optimally whole viral metagenomic experimental process, 3) evaluating the reliability of multiple displacement amplification (MDA), 4) testing the capability of long-read sequencing for improving viral metagenomic assembly, and 5) comparing the differences between viral metagenomic and bulk metagenomic approaches.

RESULTS

Our results revealed that the optimal cycle number for PCR amplification is 15. We verified the reliability of MDA and the effectiveness of long-read sequencing. Based on our optimized results, we generated 151 high-quality viruses using the dataset combined from short-read and long-read sequencing. Genomic analysis of these viruses found that most (60.3%) of them were previously unknown and showed a remarkable diversity of viral functions, especially the existence of 206 viral auxiliary metabolic genes. Finally, we uncovered significant differences in the efficiency and coverage of viral identification between viral metagenomic and bulk metagenomic approaches.

CONCLUSIONS

Our study demonstrates the potential of optimized experiment and sequencing strategies in uncovering viral genomes from fecal specimens, which will facilitate future research about the genome-level characterization of complex viral communities.

Viral biogeography of the mammalian gut and parenchymal organs

The mammalian virome has been linked to health and disease but our understanding of how it is structured along the longitudinal axis of the mammalian gastrointestinal tract (GIT) and other organs is limited. Here, we report a metagenomic analysis of the prokaryotic and eukaryotic virome occupying luminal and mucosa-associated habitats along the GIT, as well as parenchymal organs (liver, lung and spleen), in two representative mammalian species, the domestic pig and rhesus macaque (six animals per species). Luminal samples from the large intestine of both mammals harboured the highest loads and diversity of bacteriophages (class Caudoviricetes, family Microviridae and others). Mucosal samples contained much lower viral loads but a higher proportion of eukaryotic viruses (families Astroviridae, Caliciviridae, Parvoviridae). Parenchymal organs contained bacteriophages of gut origin, in addition to some eukaryotic viruses. Overall, GIT virome composition was specific to anatomical region and host species. Upper GIT and mucosa-specific viruses were greatly under-represented in distal colon samples (a proxy for faeces). Nonetheless, certain viral and phage species were ubiquitous in all samples from the oral cavity to the distal colon. The dataset and its accompanying methodology may provide an important resource for future work investigating the biogeography of the mammalian gut virome.

Diversity of CRESS DNA Viruses in Squamates Recapitulates Hosts Dietary and Environmental Sources of Exposure

Replication-associated protein (Rep)-encoding single-stranded (CRESS) DNA viruses comprise viruses with covalently closed, circular, single-stranded DNA (ssDNA) genomes, and are considered the smallest known autonomously replicating, capsid-encoding animal pathogens. CRESS DNA viruses (phylum Cressdnaviricota) encompass several viral families including Circoviridae. Circoviruses are classified into two genera, Circovirus and Cyclovirus, and they are known to cause fatal diseases in birds and pigs. Circoviruses have also been identified in human stools, blood, and cerebrospinal fluid (CSF), as well as in various wild and domestic vertebrates, including reptiles. The synanthropic presence of Squamata reptiles has increased in the last century due to the anthropic pressure, which has shifted forested animal behavior to an urban and peri-urban adaptation. In this paper, we explored the diversity of CRESS DNA viruses in Squamata reptiles from different Italian areas representative of the Mediterranean basin. CRESS DNA viruses were detected in 31.7% (33/104) of sampled lizards and geckoes. Different CRESS DNA viruses likely reflected dietary composition or environmental contamination and included avian-like (n = 3), dog (n = 4), bat-like (n = 1), goat-like (n = 1), rodent-like (n = 4), and insect-like (n = 2) viruses. Rep sequences of at least two types of human-associated cycloviruses (CyV) were identified consistently, regardless of geographic location, namely, TN9-like (n = 11) and TN12-like (n = 6). A third human-associated CyV, TN25-like, was detected in a single sample. The complete genome of human-like CyVs, of a rodent-like, insect-like, and of a bat-like virus were generated. Collectively, the results recapitulate hosts dietary and environmental sources of exposure and may suggest unexpected ecological niches for some CRESS DNA viruses.

IMPORTANCE CRESS DNA viruses are significant pathogens of birds and pigs and have been detected repeatedly in human samples (stools, serum, and cerebrospinal fluid), both from healthy individuals and from patients with neurological disease, eliciting in 2013 a risk assessment by the European Centre for Disease Prevention and Control (ECDC). Sequences of CRESS DNA viruses previously reported in humans (TN9, TN12, and TN25), and detected in different animal species (e.g., birds, dogs, and bats) were herein detected in fecal samples of synanthropic squamates (geckos and lizards). The complete genome sequence of six viruses was generated. This study extends the information on the genetic diversity and ecology of CRESS DNA viruses. Because geckos and lizards are synanthropic animals, a role in sustaining CRESS DNA virus circulation and increasing viral pressure in the environment is postulated.

Unravelling the ecological complexity of soil viromes: Challenges and opportunities

Viruses are extremely abundant and ubiquitous in soil, and significantly contribute to various terrestrial ecosystem processes such as biogeochemical nutrient cycling, microbiome regulation and community assembly, and host evolutionary dynamics. Despite their numerous dominance and functional importance, understanding soil viral ecology is a formidable challenge, because of the technological challenges to characterize the abundance, diversity and community compositions of viruses, and their interactions with other organisms in the complex soil environment. Viruses may engage in a myriad of biological interactions within soil food webs across a broad range of spatiotemporal scales and are exposed to various biotic and abiotic disturbances. Current studies on the soil viromes, however, often describe the complexity of their tremendous diversity, but lack of exploring their potential ecological roles. In this article, we summarized the major methods to decipher the ecology of soil viruses, discussed biotic and abiotic factors and global change factors that shape the diversity and composition of soil viromes, and the ecological roles of soil viruses. We also proposed a new framework to understand the ecological complexity of viruses from micro to macro ecosystem scales and to predict and unravel their activities in terrestrial ecosystems.

Bacteriophages and their potential for treatment of gastrointestinal diseases

Although bacteriophages have been overshadowed as therapeutic agents by antibiotics for decades, the emergence of multidrug-resistant bacteria and a better understanding of the role of the gut microbiota in human health and disease have brought them back into focus. In this Perspective, we briefly introduce basic phage biology and summarize recent discoveries about phages in relation to their role in the gut microbiota and gastrointestinal diseases, such as inflammatory bowel disease and chronic liver disease. In addition, we review preclinical studies and clinical trials of phage therapy for enteric disease and explore current challenges and potential future directions.

Phages in the Gut Ecosystem

Phages, short for bacteriophages, are viruses that specifically infect bacteria and are the most abundant biological entities on earth found in every explored environment, from the deep sea to the Sahara Desert. Phages are abundant within the human biome and are gaining increasing recognition as potential modulators of the gut ecosystem. For example, they have been connected to gastrointestinal diseases and the treatment efficacy of Fecal Microbiota Transplant. The ability of phages to modulate the human gut microbiome has been attributed to the predation of bacteria or the promotion of bacterial survival by the transfer of genes that enhance bacterial fitness upon infection. In addition, phages have been shown to interact with the human immune system with variable outcomes. Despite the increasing evidence supporting the importance of phages in the gut ecosystem, the extent of their influence on the shape of the gut ecosystem is yet to be fully understood. Here, we discuss evidence for phage modulation of the gut microbiome, postulating that phages are pivotal contributors to the gut ecosystem dynamics. We therefore propose novel research questions to further elucidate the role(s) that they have within the human ecosystem and its impact on our health and well-being.

Predicting the capsid architecture of phages from metagenomic data

Tailed phages are viruses that infect bacteria and are the most abundant biological entities on Earth. Their ecological, evolutionary, and biogeochemical roles in the planet stem from their genomic diversity. Known tailed phage genomes range from 10 to 735 kilobase pairs thanks to the size variability of the protective protein capsids that store them. However, the role of tailed phage capsids’ diversity in ecosystems is unclear. A fundamental gap is the difficulty of associating genomic information with viral capsids in the environment. To address this problem, here, we introduce a computational approach to predict the capsid architecture (T-number) of tailed phages using the sequence of a single gene—the major capsid protein. This approach relies on an allometric model that relates the genome length and capsid architecture of tailed phages. This allometric model was applied to isolated phage genomes to generate a library that associated major capsid proteins and putative capsid architectures. This library was used to train machine learning methods, and the most computationally scalable model investigated (random forest) was applied to human gut metagenomes. Compared to isolated phages, the analysis of gut data reveals a large abundance of mid-sized (T = 7) capsids, as expected, followed by a relatively large frequency of jumbo-like tailed phage capsids (T ≥ 25) and small capsids (T = 4) that have been under-sampled. We discussed how to increase the method’s accuracy and how to extend the approach to other viruses. The computational pipeline introduced here opens the doors to monitor the ongoing evolution and selection of viral capsids across ecosystems.

iVirus 2.0: Cyberinfrastructure-supported tools and data to power DNA virus ecology

Microbes drive myriad ecosystem processes, but under strong influence from viruses. Because studying viruses in complex systems requires different tools than those for microbes, they remain underexplored. To combat this, we previously aggregated double-stranded DNA (dsDNA) virus analysis capabilities and resources into 'iVirus' on the CyVerse collaborative cyberinfrastructure. Here we substantially expand iVirus's functionality and accessibility, to iVirus 2.0, as follows. First, core iVirus apps were integrated into the Department of Energy's Systems Biology KnowledgeBase (KBase) to provide an additional analytical platform. Second, at CyVerse, 20 software tools (apps) were upgraded or added as new tools and capabilities. Third, nearly 20-fold more sequence reads were aggregated to capture new data and environments. Finally, documentation, as "live" protocols, was updated to maximize user interaction with and contribution to infrastructure development. Together, iVirus 2.0 serves as a uniquely central and accessible analytical platform for studying how viruses, particularly dsDNA viruses, impact diverse microbial ecosystems.

Phage Annotation Guide: Guidelines for Assembly and High-Quality Annotation

All sequencing projects of bacteriophages (phages) should seek to report an accurate and comprehensive annotation of their genomes. This article defines 14 questions for those new to phage genomics that should be addressed before submitting a genome sequence to the International Nucleotide Sequence Database Collaboration or writing a publication.

Analysis of viromes and microbiomes from pig fecal samples reveals that phages and prophages rarely carry antibiotic resistance genes

Understanding the transmission of antibiotic resistance genes (ARGs) is critical for human health. For this, it is necessary to identify which type of mobile genetic elements is able to spread them from animal reservoirs into human pathogens. Previous research suggests that in pig feces, ARGs may be encoded by bacteriophages. However, convincing proof for phage-encoded ARGs in pig viromes is still lacking, because of bacterial DNA contaminating issues. We collected 14 pig fecal samples and performed deep sequencing on both highly purified viral fractions and total microbiota, in order to investigate phage and prophage-encoded ARGs. We show that ARGs are absent from the genomes of active, virion-forming phages (below 0.02% of viral contigs from viromes), but present in three prophages, representing 0.02% of the viral contigs identified in the microbial dataset. However, the corresponding phages were not detected in the viromes, and their genetic maps suggest they might be defective. We conclude that among pig fecal samples, phages and prophages rarely carry ARG. Furthermore, our dataset allows for the first time a comprehensive view of the interplay between prophages and viral particles, and uncovers two large clades, inoviruses and Oengus-like phages.

Long-read metagenomics of multiple displacement amplified DNA of low-biomass human gut phageomes by SACRA pre-processing chimeric reads

The human gut bacteriophage community (phageome) plays an important role in the host's health and disease; however, the entire structure is poorly understood, partly owing to the generation of many incomplete genomes in conventional short-read metagenomics. Here, we show long-read metagenomics of amplified DNA of low-biomass phageomes with multiple displacement amplification (MDA), involving the development of a novel bioinformatics tool, split amplified chimeric read algorithm (SACRA), that efficiently pre-processed numerous chimeric reads generated through MDA. Using five samples, SACRA markedly reduced the average chimera ratio from 72% to 1.5% in PacBio reads with an average length of 1.8 kb. De novo assembly of chimera-less PacBio long reads reconstructed contigs of ≥5 kb with an average proportion of 27%, which was 1% in contigs from MiSeq short reads, thereby dramatically improving contig length and genome completeness. Comparison of PacBio and MiSeq contigs found MiSeq contig fragmentations frequently near local repeats and hypervariable regions in the phage genomes, and those caused by multiple homologous phage genomes coexisting in the community. We also developed a reference-independent method to assess the completeness of the linear phage genomes. Overall, we established a SACRA-coupled long-read metagenomics robust to highly diverse gut phageomes, identifying high-quality circular and linear phage genomes with adequate sequence quantity.

Adaptation of the polony technique to quantify Gokushovirinae, a diverse group of single-stranded DNA phage

Advances in metagenomics have revealed the ubiquity of single-stranded DNA (ssDNA) phage belonging to the subfamily Gokushovirinae in the oceans; however, the abundance and ecological roles of this group are unknown. Here, we quantify gokushoviruses through adaptation of the polony method, in which viral template DNA is immobilized in a gel, amplified by PCR, and subsequently detected by hybridization. Primers and probes for this assay were designed based on PCR amplicon diversity of gokushovirus major capsid protein gene sequences from a depth profile in the Gulf of Aqaba, Red Sea sampled in September 2015. At ≥95% identity, these 87 gokushovirus sequences formed 14 discrete clusters with the largest clades showing distinct depth distributions. The application of the polony method enabled the first quantification of gokushoviruses in any environment. The gokushoviruses were most abundant in the upper 40 m of the stratified water column, with a subsurface peak in abundance of 1.26 × 10⁵ viruses ml^-1 . These findings suggest that discrete gokushovirus genotypes infect bacterial hosts that differentially partition in the water column. Since the designed primers and probe are conserved across marine ecosystems, this polony method can be applied broadly for the quantification of gokushoviruses throughout the global oceans.

Prokaryotic viruses impact functional microorganisms in nutrient removal and carbon cycle in wastewater treatment plants

As one of the largest biotechnological applications, activated sludge (AS) systems in wastewater treatment plants (WWTPs) harbor enormous viruses, with 10-1,000-fold higher concentrations than in natural environments. However, the compositional variation and host-connections of AS viruses remain poorly explored. Here, we report a catalogue of ~50,000 prokaryotic viruses from six WWTPs, increasing the number of described viral species of AS by 23-fold, and showing the very high viral diversity which is largely unknown (98.4-99.6% of total viral contigs). Most viral genera are represented in more than one AS system with 53 identified across all. Viral infection widely spans 8 archaeal and 58 bacterial phyla, linking viruses with aerobic/anaerobic heterotrophs, and other functional microorganisms controlling nitrogen/phosphorous removal. Notably, Mycobacterium, notorious for causing AS foaming, is associated with 402 viral genera. Our findings expand the current AS virus catalogue and provide reference for the phage treatment to control undesired microorganisms in WWTPs.

Gut dsDNA virome shows diversity and richness alterations associated with childhood obesity and metabolic syndrome

Changes in the human gut microbiome are associated with obesity and metabolic syndrome, but the role of the gut virome in both diseases remains largely unknown. We characterized the gut dsDNA virome of 28 school-aged children with healthy normal-weight (NW, n = 10), obesity (O, n = 10), and obesity with metabolic syndrome (OMS, n = 8), using metagenomic sequencing of virus-like particles (VLPs) from fecal samples. The virome classification confirmed the bacteriophages' dominance, mainly composed of Caudovirales. Notably, phage richness and diversity of individuals with O and OMS tended to increase, while the VLP abundance remained the same among all groups. Of the 4,611 phage contigs composing the phageome, 48 contigs were highly prevalent in ≥80% of individuals, suggesting high inter-individual phage diversity. The abundance of several contigs correlated with gut bacterial taxa; and with anthropometric and biochemical parameters altered in O and OMS. To our knowledge, this gut phageome represents one of the largest datasets and suggests disease-specific phage alterations.

HoPhage: an ab initio tool for identifying hosts of phage fragments from metaviromes

SUMMARY

We present HoPhage (Host of Phage) to identify the host of a given phage fragment from metavirome data at the genus level. HoPhage integrates two modules using a deep learning algorithm and a Markov chain model, respectively. HoPhage achieves 47.90% and 82.47% mean accuracy at the genus and phylum levels for ∼1-kb long artificial phage fragments when predicting host among 50 genera, representing 7.54-20.22% and 13.55-24.31% improvement, respectively. By testing on three real virome samples, HoPhage yields 81.11% mean accuracy at the genus level within a much broader candidate host range.

AVAILABILITY AND IMPLEMENTATION

HoPhage is available at http://cqb.pku.edu.cn/ZhuLab/HoPhage/data/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Glacier ice archives nearly 15,000-year-old microbes and phages

BACKGROUND

Glacier ice archives information, including microbiology, that helps reveal paleoclimate histories and predict future climate change. Though glacier-ice microbes are studied using culture or amplicon approaches, more challenging metagenomic approaches, which provide access to functional, genome-resolved information and viruses, are under-utilized, partly due to low biomass and potential contamination.

RESULTS

We expand existing clean sampling procedures using controlled artificial ice-core experiments and adapted previously established low-biomass metagenomic approaches to study glacier-ice viruses. Controlled sampling experiments drastically reduced mock contaminants including bacteria, viruses, and free DNA to background levels. Amplicon sequencing from eight depths of two Tibetan Plateau ice cores revealed common glacier-ice lineages including Janthinobacterium, Polaromonas, Herminiimonas, Flavobacterium, Sphingomonas, and Methylobacterium as the dominant genera, while microbial communities were significantly different between two ice cores, associating with different climate conditions during deposition. Separately, ~355- and ~14,400-year-old ice were subject to viral enrichment and low-input quantitative sequencing, yielding genomic sequences for 33 vOTUs. These were virtually all unique to this study, representing 28 novel genera and not a single species shared with 225 environmentally diverse viromes. Further, 42.4% of the vOTUs were identifiable temperate, which is significantly higher than that in gut, soil, and marine viromes, and indicates that temperate phages are possibly favored in glacier-ice environments before being frozen. In silico host predictions linked 18 vOTUs to co-occurring abundant bacteria (Methylobacterium, Sphingomonas, and Janthinobacterium), indicating that these phages infected ice-abundant bacterial groups before being archived. Functional genome annotation revealed four virus-encoded auxiliary metabolic genes, particularly two motility genes suggest viruses potentially facilitate nutrient acquisition for their hosts. Finally, given their possible importance to methane cycling in ice, we focused on Methylobacterium viruses by contextualizing our ice-observed viruses against 123 viromes and prophages extracted from 131 Methylobacterium genomes, revealing that the archived viruses might originate from soil or plants.

CONCLUSIONS

Together, these efforts further microbial and viral sampling procedures for glacier ice and provide a first window into viral communities and functions in ancient glacier environments. Such methods and datasets can potentially enable researchers to contextualize new discoveries and begin to incorporate glacier-ice microbes and their viruses relative to past and present climate change in geographically diverse regions globally. Video Abstract.

Application of machine learning in bacteriophage research

Phages are one of the key components in the structure, dynamics, and interactions of microbial communities in different bins. It has a clear impact on human health and the food industry. Bacteriophage characterization using in vitro approaches are time/cost consuming and laborious tasks. On the other hand, with the advent of new high-throughput sequencing technology, the development of a powerful computational framework to characterize the newly identified bacteriophages is inevitable for future research. Machine learning includes powerful techniques that enable the analysis of complex datasets for knowledge discovery and pattern recognition. In this study, we have conducted a comprehensive review of machine learning methods application using different types of features were applied in various aspects of bacteriophage research including, automated curation, identification, classification, host species recognition, virion protein identification, and life cycle prediction. Moreover, potential limitations and advantages of the developed frameworks were discussed.

Integrating Viral Metagenomics into an Ecological Framework

Viral metagenomics has expanded our knowledge of the ecology of uncultured viruses, within both environmental (e.g., terrestrial and aquatic) and host-associated (e.g., plants and animals, including humans) contexts. Here, we emphasize the implementation of an ecological framework in viral metagenomic studies to address questions in virology rarely considered ecological, which can change our perception of viruses and how they interact with their surroundings. An ecological framework explicitly considers diverse variants of viruses in populations that make up communities of interacting viruses, with ecosystem-level effects. It provides a structure for the study of the diversity, distributions, dynamics, and interactions of viruses with one another, hosts, and the ecosystem, including interactions with abiotic factors. An ecological framework in viral metagenomics stands poised to broadly expand our knowledge in basic and applied virology. We highlight specific fundamental research needs to capitalize on its potential and advance the field. Expected final online publication date for the Annual Review of Virology, Volume 8 is September 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Comparison of ultrafiltration and iron chloride flocculation in the preparation of aquatic viromes from contrasting sample types

Viral metagenomes (viromes) are a valuable untargeted tool for studying viral diversity and the central roles viruses play in host disease, ecology, and evolution. Establishing effective methods to concentrate and purify viral genomes prior to sequencing is essential for high quality viromes. Using virus spike-and-recovery experiments, we stepwise compared two common approaches for virus concentration, ultrafiltration and iron chloride flocculation, across diverse matrices: wastewater influent, wastewater secondary effluent, river water, and seawater. Viral DNA was purified by removing cellular DNA via chloroform cell lysis, filtration, and enzymatic degradation of extra-viral DNA. We found that viral genomes were concentrated 1-2 orders of magnitude more with ultrafiltration than iron chloride flocculation for all matrices and resulted in higher quality DNA suitable for amplification-free and long-read sequencing. Given its widespread use and utility as an inexpensive field method for virome sampling, we nonetheless sought to optimize iron flocculation. We found viruses were best concentrated in seawater with five-fold higher iron concentrations than the standard used, inhibition of DNase activity reduced purification effectiveness, and five-fold more iron was needed to flocculate viruses from freshwater than seawater—critical knowledge for those seeking to apply this broadly used method to freshwater virome samples. Overall, our results demonstrated that ultrafiltration and purification performed better than iron chloride flocculation and purification in the tested matrices. Given that the method performance depended on the solids content and salinity of the samples, we suggest spike-and-recovery experiments be applied when concentrating and purifying sample types that diverge from those tested here.

Functional Restoration of Bacteriomes and Viromes by Fecal Microbiota Transplantation

BACKGROUND & AIMS

Fecal microbiota transplantation (FMT) is an effective therapy for recurrent Clostridioides difficile infection (rCDI). However, the overall mechanisms underlying FMT success await comprehensive elucidation, and the safety of FMT has recently become a serious concern because of the occurrence of drug-resistant bacteremia transmitted by FMT. We investigated whether functional restoration of the bacteriomes and viromes by FMT could be an indicator of successful FMT.

METHODS

The human intestinal bacteriomes and viromes from 9 patients with rCDI who had undergone successful FMT and their donors were analyzed. Prophage-based and CRISPR spacer-based host bacteria-phage associations in samples from recipients before and after FMT and in donor samples were examined. The gene functions of intestinal microorganisms affected by FMT were evaluated.

RESULTS

Metagenomic sequencing of both the viromes and bacteriomes revealed that FMT does change the characteristics of intestinal bacteriomes and viromes in recipients after FMT compared with those before FMT. In particular, many Proteobacteria, the fecal abundance of which was high before FMT, were eliminated, and the proportion of Microviridae increased in recipients. Most temperate phages also behaved in parallel with the host bacteria that were altered by FMT. Furthermore, the identification of bacterial and viral gene functions before and after FMT revealed that some distinctive pathways, including fluorobenzoate degradation and secondary bile acid biosynthesis, were significantly represented.

CONCLUSIONS

The coordinated action of phages and their host bacteria restored the recipients' intestinal flora. These findings show that the restoration of intestinal microflora functions reflects the success of FMT.

Extended Evaluation of Viral Diversity in Lake Baikal through Metagenomics

Lake Baikal is a unique oligotrophic freshwater lake with unusually cold conditions and amazing biological diversity. Studies of the lake's viral communities have begun recently, and their full diversity is not elucidated yet. Here, we performed DNA viral metagenomic analysis on integral samples from four different deep-water and shallow stations of the southern and central basins of the lake. There was a strict distinction of viral communities in areas with different environmental conditions. Comparative analysis with other freshwater lakes revealed the highest similarity of Baikal viromes with those of the Asian lakes Soyang and Biwa. Analysis of new data, together with previously published data allowed us to get a deeper insight into the diversity and functional potential of Baikal viruses; however, the true diversity of Baikal viruses in the lake ecosystem remains still unknown. The new metaviromic data will be useful for future studies of viral composition, distribution, and the dynamics associated with global climatic and anthropogenic impacts on this ecosystem.

A library preparation optimized for metagenomics of RNA viruses

Our understanding of the viral communities associated to animals has not yet reached the level attained on the bacteriome. This situation is due to, among others, technical challenges in adapting metagenomics using high-throughput sequencing to the study of RNA viromes in animals. Although important developments have been achieved in most steps of viral metagenomics, there is yet a key step that has received little attention: the library preparation. This situation differs from bacteriome studies in which developments in library preparation have largely contributed to the democratisation of metagenomics. Here, we present a library preparation optimized for metagenomics of RNA viruses from insect vectors of viral diseases. The library design allows a simple PCR-based preparation, such as those routinely used in bacterial metabarcoding, that is adapted to shotgun sequencing as required in viral metagenomics. We first optimized our library preparation using mock viral communities and then validated a full metagenomic approach incorporating our preparation in two pilot studies with field-caught insect vectors; one including a comparison with a published metagenomic protocol. Our approach provided a fold increase in virus-like sequences compared to other studies, and nearly-full genomes from new virus species. Moreover, our results suggested conserved trends in virome composition within a population of a mosquito species. Finally, the sensitivity of our approach was compared to a commercial diagnostic PCR for the detection of an arbovirus in field-caught insect vectors. Our approach could facilitate studies on viral communities from animals and the democratization of metagenomics in community ecology of viruses.

Potential virus-mediated nitrogen cycling in oxygen-depleted oceanic waters

Viruses play an important role in the ecology and biogeochemistry of marine ecosystems. Beyond mortality and gene transfer, viruses can reprogram microbial metabolism during infection by expressing auxiliary metabolic genes (AMGs) involved in photosynthesis, central carbon metabolism, and nutrient cycling. While previous studies have focused on AMG diversity in the sunlit and dark ocean, less is known about the role of viruses in shaping metabolic networks along redox gradients associated with marine oxygen minimum zones (OMZs). Here, we analyzed relatively quantitative viral metagenomic datasets that profiled the oxygen gradient across Eastern Tropical South Pacific (ETSP) OMZ waters, assessing whether OMZ viruses might impact nitrogen (N) cycling via AMGs. Identified viral genomes encoded six N-cycle AMGs associated with denitrification, nitrification, assimilatory nitrate reduction, and nitrite transport. The majority of these AMGs (80%) were identified in T4-like Myoviridae phages, predicted to infect Cyanobacteria and Proteobacteria, or in unclassified archaeal viruses predicted to infect Thaumarchaeota. Four AMGs were exclusive to anoxic waters and had distributions that paralleled homologous microbial genes. Together, these findings suggest viruses modulate N-cycling processes within the ETSP OMZ and may contribute to nitrogen loss throughout the global oceans thus providing a baseline for their inclusion in the ecosystem and geochemical models.

Viruses in fermented foods: are they good or bad? Two sides of the same coin

The emergence of Coronavirus disease 2019 as a global pandemic has increased popular concerns about diseases caused by viruses. Fermented foods containing high loads of viable fungi and bacteria are potential sources for virus contamination. The most common include viruses that infect bacteria (bacteriophage) and yeasts reported in fermented milks, sausages, vegetables, wine, sourdough, and cocoa beans. Recent molecular studies have also associated fermented foods as vehicles for pathogenic human viruses. Human noroviruses, rotavirus, and hepatitis virus have been identified in different fermented foods through multiple routes. No severe acute respiratory syndrome Coronavirus 2 (SARS-CoV-2) virus or close members were found in fermented foods to date. However, the occurrence/persistence of other pathogenic viruses reveals a potential vulnerability of fermented foods to SARS-CoV-2 contamination. On the other side of the coin, some bacteriophages are being suggested for improving the fermentation process and food safety, as well as owing potential probiotic properties in modern fermented foods. This review will address the diversity and characteristics of viruses associated with fermented foods and what has been changed after a short introduction to the most common next-generation sequencing platforms. Also, the risk of SARS-CoV-2 transmission via fermented foods and preventive measures will be discussed.

Hybrid assembly of an agricultural slurry virome reveals a diverse and stable community with the potential to alter the metabolism and virulence of veterinary pathogens

BACKGROUND

Viruses are the most abundant biological entities on Earth, known to be crucial components of microbial ecosystems. However, there is little information on the viral community within agricultural waste. There are currently ~ 2.7 million dairy cattle in the UK producing 7-8% of their own bodyweight in manure daily, and 28 million tonnes annually. To avoid pollution of UK freshwaters, manure must be stored and spread in accordance with guidelines set by DEFRA. Manures are used as fertiliser, and widely spread over crop fields, yet little is known about their microbial composition. We analysed the virome of agricultural slurry over a 5-month period using short and long-read sequencing.

RESULTS

Hybrid sequencing uncovered more high-quality viral genomes than long or short-reads alone; yielding 7682 vOTUs, 174 of which were complete viral genomes. The slurry virome was highly diverse and dominated by lytic bacteriophage, the majority of which represent novel genera (~ 98%). Despite constant influx and efflux of slurry, the composition and diversity of the slurry virome was extremely stable over time, with 55% of vOTUs detected in all samples over a 5-month period. Functional annotation revealed a diverse and abundant range of auxiliary metabolic genes and novel features present in the community, including the agriculturally relevant virulence factor VapE, which was widely distributed across different phage genera that were predicted to infect several hosts. Furthermore, we identified an abundance of phage-encoded diversity-generating retroelements, which were previously thought to be rare on lytic viral genomes. Additionally, we identified a group of crAssphages, including lineages that were previously thought only to be found in the human gut.

CONCLUSIONS

The cattle slurry virome is complex, diverse and dominated by novel genera, many of which are not recovered using long or short-reads alone. Phages were found to encode a wide range of AMGs that are not constrained to particular groups or predicted hosts, including virulence determinants and putative ARGs. The application of agricultural slurry to land may therefore be a driver of bacterial virulence and antimicrobial resistance in the environment. Video abstract.

Probing the "Dark Matter" of the Human Gut Phageome: Culture Assisted Metagenomics Enables Rapid Discovery and Host-Linking for Novel Bacteriophages

Recent years have been marked by the growing interest towards virulent and temperate bacteriophage populations inhabiting the human lower gastrointestinal tract – the gut phageome. A number of studies demonstrated high levels of specificity and temporal stability of individual gut phageomes, as well as their specific alterations in disease cohorts, in parallel with changes in the bacteriome. It has been speculated that phages might have an active role in shaping the taxonomic composition and functional properties of the human gut bacteriome. An overwhelming majority of gut bacteriophages, however, remain uncultured, unclassified, and their specific hosts and infection strategies are still unknown. They are often referred to as “the viral dark matter”. A possible breakthrough in understanding of the phageome can only become possible when a significant proportion of the “the viral dark matter” is identified and linked to bacterial hosts. Here, we describe a method that enables rapid discovery and host-linking of novel bacteriophages in the gut via a combination of serial enrichment cultures and shotgun metagenomics of viral DNA. Using this approach dozens of novel and previously known bacteriophages were detected, including the ones infecting difficult-to-culture anaerobic bacteria. The majority of phages failed to produce lysis and propagate on host cultures in traditional assays. The newly identified phages include representatives of Siphoviridae, Myoviridae, Podoviridae, and crAss-like viruses, infecting diverse bacterial taxa of Bacteroidetes, Firmicutes, Actinobacteria, Verrucomicrobia and Proteobacteria phyla. The proposed new method has a potential for high-throughput screening applications for mass discovery of new phages in different environments.