Viral metagenomics

Unveiling microbial worlds: exploring viral metagenomics among waste pickers at Latin America's largest dumpsite

Waste pickers constitute a marginalized demographic engaged in the collection of refuse, facing considerable occupational hazards that heighten their susceptibility to contract infectious diseases. Moreover, waste pickers contend with societal stigmatization and encounter barriers to accessing healthcare services. To explore the viral profile of waste pickers potentially linked to their occupational environment, we conducted a metagenomic analysis on 120 plasma specimens sampled from individuals employed at the Cidade Estrutural dumpsite in Brasilia city, Brazil. In total, 60 blood donors served as a comparative control group. Specimens were pooled and subjected to Illumina NextSeq 2000 sequencing. Viral abundance among waste pickers revealed the presence of significant pathogens, including HIV, HCV, and Chikungunya, which were not detected in the control group. Additionally, elevated levels of anelloviruses and Human pegivirus-1 were noted, with a comparable incidence in the control group. These findings underscore the utility of metagenomics in identifying clinically relevant viral agents within underserved populations. The implications of this study extend to informing public health policies aimed at surveilling infectious diseases among individuals facing socioeconomic disparities and limited access to healthcare resources.

Identification of infectious viruses for risk-based virus testing of CHO unprocessed bulk using next-generation sequencing

It is important to increase manufacturing speed to make medicines more widely available. One bottleneck for CHO-based drug substance release is the in vitro viral (IVV) cell-based assay on unprocessed bulk. To increase process speed, we evaluate the suitability of replacing the IVV cell-based assay with next-generation sequencing (NGS). First, we outline how NGS is currently used in the pharmaceutical industry, and how it may apply to CHO virus testing. Second, we examine CHO virus contamination history. Since prior virus contaminants can replicate in the production bioreactor, we perform a literature search and classify 159 viruses as high, medium, low, or unknown risk based on their ability to infect CHO cells. Overall, the risk of virus contamination during the CHO manufacturing process is low. Only six viruses were reported to have contaminated CHO bioprocesses over the past several decades, and were primarily caused by fetal bovine serum or cell culture components. These virus contamination events can be mitigated through limitation and control of raw materials, combined with virus testing and virus clearance technologies. The list of CHO infectious viruses provides a starting framework for virus safety risk assessment and NGS development. Furthermore, ICH Q5A (R2) includes NGS as a molecular method for adventitious agent testing, paving a path forward for modernizing CHO virus testing.

Bacteria are a major determinant of Orsay virus transmission and infection in Caenorhabditis elegans

The microbiota is a key determinant of the physiology and immunity of animal hosts. The factors governing the transmissibility of viruses between susceptible hosts are incompletely understood. Bacteria serve as food for Caenorhabditis elegans and represent an integral part of the natural environment of C. elegans. We determined the effects of bacteria isolated with C. elegans from its natural environment on the transmission of Orsay virus in C. elegans using quantitative virus transmission and host susceptibility assays. We observed that Ochrobactrum species promoted Orsay virus transmission, whereas Pseudomonas lurida MYb11 attenuated virus transmission relative to the standard laboratory bacterial food Escherichia coli OP50. We found that pathogenic Pseudomonas aeruginosa strains PA01 and PA14 further attenuated virus transmission. We determined that the amount of Orsay virus required to infect 50% of a C. elegans population on P. lurida MYb11 compared with Ochrobactrum vermis MYb71 was dramatically increased, over three orders of magnitude. Host susceptibility was attenuated even further in the presence of P. aeruginosa PA14. Genetic analysis of the determinants of P. aeruginosa required for attenuation of C. elegans susceptibility to Orsay virus infection revealed a role for regulators of quorum sensing. Our data suggest that distinct constituents of the C. elegans microbiota and potential pathogens can have widely divergent effects on Orsay virus transmission, such that associated bacteria can effectively determine host susceptibility versus resistance to viral infection. Our study provides quantitative evidence for a critical role for tripartite host-virus-bacteria interactions in determining the transmissibility of viruses among susceptible hosts.

Utilizing profile hidden Markov model databases for discovering viruses from metagenomic data: a comprehensive review

Profile hidden Markov models (pHMMs) are able to achieve high sensitivity in remote homology search, making them popular choices for detecting novel or highly diverged viruses in metagenomic data. However, many existing pHMM databases have different design focuses, making it difficult for users to decide the proper one to use. In this review, we provide a thorough evaluation and comparison for multiple commonly used profile HMM databases for viral sequence discovery in metagenomic data. We characterized the databases by comparing their sizes, their taxonomic coverage, and the properties of their models using quantitative metrics. Subsequently, we assessed their performance in virus identification across multiple application scenarios, utilizing both simulated and real metagenomic data. We aim to offer researchers a thorough and critical assessment of the strengths and limitations of different databases. Furthermore, based on the experimental results obtained from the simulated and real metagenomic data, we provided practical suggestions for users to optimize their use of pHMM databases, thus enhancing the quality and reliability of their findings in the field of viral metagenomics.

Benchmarking bioinformatic virus identification tools using real-world metagenomic data across biomes

BACKGROUND

As most viruses remain uncultivated, metagenomics is currently the main method for virus discovery. Detecting viruses in metagenomic data is not trivial. In the past few years, many bioinformatic virus identification tools have been developed for this task, making it challenging to choose the right tools, parameters, and cutoffs. As all these tools measure different biological signals, and use different algorithms and training and reference databases, it is imperative to conduct an independent benchmarking to give users objective guidance.

RESULTS

We compare the performance of nine state-of-the-art virus identification tools in thirteen modes on eight paired viral and microbial datasets from three distinct biomes, including a new complex dataset from Antarctic coastal waters. The tools have highly variable true positive rates (0-97%) and false positive rates (0-30%). PPR-Meta best distinguishes viral from microbial contigs, followed by DeepVirFinder, VirSorter2, and VIBRANT. Different tools identify different subsets of the benchmarking data and all tools, except for Sourmash, find unique viral contigs. Performance of tools improved with adjusted parameter cutoffs, indicating that adjustment of parameter cutoffs before usage should be considered.

CONCLUSIONS

Together, our independent benchmarking facilitates selecting choices of bioinformatic virus identification tools and gives suggestions for parameter adjustments to viromics researchers.

Bacteria Are a Major Determinant of Orsay Virus Transmission and Infection in Caenorhabditis elegans

The microbiota is a key determinant of the physiology and immunity of animal hosts. The factors governing the transmissibility of viruses between susceptible hosts are incompletely understood. Bacteria serve as food for Caenorhabditis elegans and represent an integral part of the natural environment of C. elegans. We determined the effects of bacteria isolated with C. elegans from its natural environment on the transmission of Orsay virus in C. elegans using quantitative virus transmission and host susceptibility assays. We observed that Ochrobactrum species promoted Orsay virus transmission, whereas Pseudomonas lurida MYb11 attenuated virus transmission relative to the standard laboratory bacterial food Escherichia coli OP50. We found that pathogenic Pseudomonas aeruginosa strains PA01 and PA14 further attenuated virus transmission. We determined that the amount of Orsay virus required to infect 50% of a C. elegans population on P. lurida MYb11 compared with Ochrobactrum vermis MYb71 was dramatically increased, over three orders of magnitude. Host susceptibility was attenuated even further in presence of P. aeruginosa PA14. Genetic analysis of the determinants of P. aeruginosa required for attenuation of C. elegans susceptibility to Orsay virus infection revealed a role for regulators of quorum sensing. Our data suggest that distinct constituents of the C. elegans microbiota and potential pathogens can have widely divergent effects on Orsay virus transmission, such that associated bacteria can effectively determine host susceptibility versus resistance to viral infection. Our study provides quantitative evidence for a critical role for tripartite host-virus-bacteria interactions in determining the transmissibility of viruses among susceptible hosts.

Chapter 2: What Is Life?

The question "What is life?" has existed since the beginning of recorded history. However, the scientific and philosophical contexts of this question have changed and been refined as advancements in technology have revealed both fine details and broad connections in the network of life on Earth. Understanding the framework of the question "What is life?" is central to formulating other questions such as "Where else could life be?" and "How do we search for life elsewhere?" While many of these questions are addressed throughout the Astrobiology Primer 3.0, this chapter gives historical context for defining life, highlights conceptual characteristics shared by all life on Earth as well as key features used to describe it, discusses why it matters for astrobiology, and explores both challenges and opportunities for finding an informative operational definition.

Applications of designer phage encoding recombinant gene payloads

Advances in genetic engineering, synthetic biology, and DNA sequencing have transformed the re-emergent therapeutic bacteriophage field. The increasing rate of multidrug resistant (MDR) infections and the speed at which new bacteriophages can be isolated, sequenced, characterized, and engineered has reinvigorated phage therapy and unlocked new applications of phages for modulating bacteria. The methods used to genetically engineer bacteriophages are undergoing significant development, but identification of heterologous gene payloads with desirable activity and determination of their impact on bacteria or human cells in translationally relevant applications remain underexplored areas. Here, we discuss and categorize recombinant gene payloads for their potential outcome on phage-bacteria interactions when genetically engineered into phage genomes for expression in their bacterial hosts.

Aeromicrobiology: A global review of the cycling and relationships of bioaerosols with the atmosphere

Airborne microorganisms and biological matter (bioaerosols) play a key role in global biogeochemical cycling, human and crop health trends, and climate patterns. Their presence in the atmosphere is controlled by three main stages: emission, transport, and deposition. Aerial survival rates of bioaerosols are increased through adaptations such as ultra-violet radiation and desiccation resistance or association with particulate matter. Current research into modern concerns such as climate change, global gene transfer, and pathogenicity often neglects to consider atmospheric involvement. This comprehensive review outlines the transpiring of bioaerosols across taxa in the atmosphere, with significant focus on their interactions with environmental elements including abiotic factors (e.g., atmospheric composition, water cycle, and pollution) and events (e.g., dust storms, hurricanes, and wildfires). The aim of this review is to increase understanding and shed light on needed research regarding the interplay between global atmospheric phenomena and the aeromicrobiome. The abundantly documented bacteria and fungi are discussed in context of their cycling and human health impacts. Gaps in knowledge regarding airborne viral community, the challenges and importance of studying their composition, concentrations and survival in the air are addressed, along with understudied plant pathogenic oomycetes, and archaea cycling. Key methodologies in sampling, collection, and processing are described to provide an up-to-date picture of ameliorations in the field. We propose optimization to microbiological methods, commonly used in soil and water analysis, that adjust them to the context of aerobiology, along with other directions towards novel and necessary advancements. This review offers new perspectives into aeromicrobiology and calls for advancements in global-scale bioremediation, insights into ecology, climate change impacts, and pathogenicity transmittance.

Ecological Interaction between Bacteriophages and Bacteria in Sub-Arctic Kongsfjorden Bay, Svalbard, Norway

Marine virus diversity and their relationships with their hosts in the marine environment remain unclear. This study investigated the co-occurrence of marine DNA bacteriophages (phages) and bacteria in the sub-Arctic area of Kongsfjorden Bay in Svalbard (Norway) in April and June 2018 using metagenomics tools. Of the marine viruses identified, 48-81% were bacteriophages of the families Myoviridae, Siphoviridae, and Podoviridae. Puniceispirillum phage HMO-2011 was dominant (7.61%) in April, and Puniceispirillum phage HMO-2011 (3.32%) and Pelagibacter phage HTVC008M (3.28%) were dominant in June. Gammaproteobacteria (58%), including Eionea flava (14.3%) and Pseudomonas sabulinigri (12.2%), were dominant in April, whereas Alphaproteobacteria (87%), including Sulfitobacter profundi (51.5%) and Loktanella acticola (32.4%), were dominant in June. The alpha diversity of the bacteriophages and bacterial communities exhibited opposite patterns. The diversity of the bacterial community was higher in April and lower in June. Changes in water temperature and light can influence the relationship between bacteria and bacteriophages.

Tailed bacteriophages (Caudoviricetes) dominate the microbiome of a diseased stingless bee

Bacteriophages, viruses that infect bacterial hosts, are known to rule the dynamics and diversity of bacterial populations in a number of ecosystems. Bacterial communities residing in the gut of animals, known as the gut microbiome, have revolutionized our understanding of many diseases. However, the gut phageome, while of apparent importance in this context, remains an underexplored area of research. Here we identify for the first time genomic sequences from tailed viruses (Caudoviricetes) that are associated with the microbiome of stingless bees (Melipona quadrifasciata). Both DNA and RNA were extracted from virus particles isolated from healthy and diseased forager bees, the latter showing symptoms from an annual syndrome that only affects M. quadrifasciata. Viral contigs from previously sequenced metagenomes of healthy and diseased forager bees were used for the analyses. Using conserved proteins deduced from their genomes, we found that Caudoviricetes were only present in the worker bee gut microbiome from diseased stingless bees. The most abundant phages are phylogenetically related to phages that infect Gram-positive bacteria from the order Lactobacillales and Gram-negative bacteria from the genus Gilliamella and Bartonella, that are common honey bee symbionts. The potential implication of these viruses in the M. quadrifasciata syndrome is discussed.

Bioinformatic Approaches for Comparative Analysis of Viruses

The field of viral genomic studies has experienced an unprecedented increase in data volume. New strains of known viruses are constantly being added to the GenBank database and so are completely new species with little or no resemblance to our databases of sequences. In addition to this, metagenomic techniques have the potential to further increase the number and rate of sequenced genomes. Besides, it is important to consider that viruses have a set of unique features that often break down molecular biology dogmas, e.g., the flux of information from RNA to DNA in retroviruses and the use of RNA molecules as genomes. As a result, extracting meaningful information from viral genomes remains a challenge and standard methods for comparing the unknown and our databases of characterized sequences may need adaptations. Thus, several bioinformatic approaches and tools have been created to address the challenge of analyzing viral data. This chapter offers descriptions and protocols of some of the most important bioinformatic techniques for comparative analysis of viruses. The authors also provide comments and discussion on how viruses' unique features can affect standard analyses and how to overcome some of the major sources of problems. Protocols and topics emphasize online tools (which are more accessible to users) and give the real experience of what most bioinformaticians do in day-by-day work with command-line pipelines. The topics discussed include (1) clustering related genomes, (2) whole genome multiple sequence alignments for small RNA viruses, (3) protein alignment for marker genes and species affiliation, (4) variant calling and annotation, and (5) virome analyses and pathogen identification.

Ancestors in the Extreme: A Genomics View of Microbial Diversity in Hypersaline Aquatic Environments

The origin of eukaryotic cells, and especially naturally occurring syncytial cells, remains debatable. While a majority of our biomedical research focuses on the eukaryotic result of evolution, our data remain limiting on the prokaryotic precursors of these cells. This is particularly evident when considering extremophile biology, especially in how the genomes of organisms in extreme environments must have evolved and adapted to unique habitats. Might these rapidly diversifying organisms have created new genetic tools eventually used to enhance the evolution of the eukaryotic single nuclear or syncytial cells? Many organisms are capable of surviving, or even thriving, in conditions of extreme temperature, acidity, organic composition, and then rapidly adapt to yet new conditions. This study identified organisms found in extremes of salinity. A lake and a nearby pond in the Ethiopian Rift Valley were interrogated for life by sequencing the DNA of populations of organism collected from the water in these sites. Remarkably, a vast diversity of microbes were identified, and even though the two sites were nearby each other, the populations of organisms were distinctly different. Since these microbes are capable of living in what for humans would be inhospitable conditions, the DNA sequences identified should inform the next step in these investigations; what new gene families, or modifications to common genes, do these organisms employ to survive in these extreme conditions. The relationship between organisms and their environment can be revealed by decoding genomes of organisms living in extreme environments. These genomes disclose new biological mechanisms that enable life outside moderate environmental conditions, new gene functions for application in biotechnology, and may even result in identification of new species. In this study, we have collected samples from two hypersaline sites in the Danakil depression, the shorelines of Lake As'ale and an actively mixing salt pond called Muda'ara (MUP), to identify the microbial community by metagenomics. Shotgun sequencing was applied to high density sampling, and the relative abundance of Operational Taxonomic Units (OTUs) was calculated. Despite the broad taxonomic similarities among the salt-saturated metagenomes analyzed, MUP stood out from Lake As'ale samples. In each sample site, Archaea accounted for 95% of the total OTUs, largely to the class Halobacteria. The remaining 5% of organisms were eubacteria, with an unclassified strain of Salinibacter ruber as the dominant OTU in both the Lake and the Pond. More than 40 different genes coding for stress proteins were identified in the three sample sites of Lake As'ale, and more than 50% of the predicted stress-related genes were associated with oxidative stress response proteins. Chaperone proteins (DnaK, DnaJ, GrpE, and ClpB) were predicted, with percentage of query coverage and similarities ranging between 9.5% and 99.2%. Long reads for ClpB homologous protein from Lake As'ale metagenome datasets were modeled, and compact 3D structures were generated. Considering the extreme environmental conditions of the Danakil depression, this metagenomics dataset can add and complement other studies on unique gene functions on stress response mechanisms of thriving bio-communities that could have contributed to cellular changes leading to single and/or multinucleated eukaryotic cells.

Targeted Genome Editing of Virulent Pseudomonas Phages Using CRISPR-Cas3

The vast number of unknown phage-encoded ORFan genes and limited insights into the genome organization of phages illustrate the need for efficient genome engineering tools to study bacteriophage genes in their natural context. In addition, there is an application-driven desire to alter phage properties, which is hampered by time constraints for phage genome engineering in the bacterial host. We here describe an optimized CRISPR-Cas3 system in Pseudomonas for straightforward editing of the genome of virulent bacteriophages. The two-vector system combines a broad host range CRISPR-Cas3 targeting plasmid with a SEVA plasmid for homologous directed repair, which enables the creation of clean deletions, insertions, or substitutions in the phage genome within a week. After creating the two plasmids separately, a co-transformation to P. aeruginosa cells is performed. A subsequent infection with the targeted phage allows the CRISPR-Cas3 system to cut the DNA specifically and facilitate or select for homologous recombination. This system has also been successfully applied for P. aeruginosa and Pseudomonas putida genome engineering. The method is straightforward, efficient, and universal, enabling to extrapolate the system to other phage-host pairs.

Computational host range prediction-The good, the bad, and the ugly

The rapid emergence and spread of antimicrobial resistance across the globe have prompted the usage of bacteriophages (i.e. viruses that infect bacteria) in a variety of applications ranging from agriculture to biotechnology and medicine. In order to effectively guide the application of bacteriophages in these multifaceted areas, information about their host ranges-that is the bacterial strains or species that a bacteriophage can successfully infect and kill-is essential. Utilizing sixteen broad-spectrum (polyvalent) bacteriophages with experimentally validated host ranges, we here benchmark the performance of eleven recently developed computational host range prediction tools that provide a promising and highly scalable supplement to traditional, but laborious, experimental procedures. We show that machine- and deep-learning approaches offer the highest levels of accuracy and precision-however, their predominant predictions at the species- or genus-level render them ill-suited for applications outside of an ecosystems metagenomics framework. In contrast, only moderate sensitivity (<80 per cent) could be reached at the strain-level, albeit at low levels of precision (<40 per cent). Taken together, these limitations demonstrate that there remains room for improvement in the active scientific field of in silico host prediction to combat the challenge of guiding experimental designs to identify the most promising bacteriophage candidates for any given application.

Biophysical basis of filamentous phage tactoid-mediated antibiotic tolerance in P. aeruginosa

Inoviruses are filamentous phages infecting numerous prokaryotic phyla. Inoviruses can self-assemble into mesoscale structures with liquid-crystalline order, termed tactoids, which protect bacterial cells in Pseudomonas aeruginosa biofilms from antibiotics. Here, we investigate the structural, biophysical, and protective properties of tactoids formed by the P. aeruginosa phage Pf4 and Escherichia coli phage fd. A cryo-EM structure of the capsid from fd revealed distinct biochemical properties compared to Pf4. Fd and Pf4 formed tactoids with different morphologies that arise from differing phage geometries and packing densities, which in turn gave rise to different tactoid emergent properties. Finally, we showed that tactoids formed by either phage protect rod-shaped bacteria from antibiotic treatment, and that direct association with a tactoid is required for protection, demonstrating the formation of a diffusion barrier by the tactoid. This study provides insights into how filamentous molecules protect bacteria from extraneous substances in biofilms and in host-associated infections.

Exploring the Diversity of Plant-Associated Viruses and Related Viruses in Riverine Freshwater Samples Collected in Berlin, Germany

Plant-infecting RNA viruses from 30 families and floating genera, as well as a great number of uncultured as yet-unclassified plant-associated viruses have been described. Even so, the plant RNA virosphere is still underexplored. RNA extracted from enriched virus particles of 50 L water samples from the Teltow Canal and the Havel River in Berlin, Germany, was sequenced using Illumina next-generation sequencing. Sequences were searched for plant viruses with BLAST and DIAMOND. Phylogenetic analyses were conducted with IQ-TREE 2. Altogether, 647 virus sequences greater than 1 kb were detected and further analyzed. These data revealed the presence of accepted and novel viruses related to Albetovirus, Alphaflexiviridae, Aspiviridae, Bromoviridae, Endornaviridae, Partitiviridae, Potyviridae, Solemoviridae, Tombusviridae and Virgaviridae. The vast majority of the sequences were novel and could not be taxonomically assigned. Several tombus- and endorna-like viruses make use of alternative translation tables that suggest unicellular green algae, ciliates, or diplomonades as their hosts. The identification of 27 albeto-like satellite viruses increases available sequence data five-fold. Sixteen new poty-like viruses align with other poty-like viruses in a link that combines the Astroviridae and Potyviridae families. Further, the identification of viruses with peptidase A6-like and peptidase A21-like capsid proteins suggests horizontal gene transfer in the evolution of these viruses.

Zea mays genotype influences microbial and viral rhizobiome community structure

Plant genotype is recognized to contribute to variations in microbial community structure in the rhizosphere, soil adherent to roots. However, the extent to which the viral community varies has remained poorly understood and has the potential to contribute to variation in soil microbial communities. Here we cultivated replicates of two Zea mays genotypes, parviglumis and B73, in a greenhouse and harvested the rhizobiome (rhizoplane and rhizosphere) to identify the abundance of cells and viruses as well as rhizobiome microbial and viral community using 16S rRNA gene amplicon sequencing and genome resolved metagenomics. Our results demonstrated that viruses exceeded microbial abundance in the rhizobiome of parviglumis and B73 with a significant variation in both the microbial and viral community between the two genotypes. Of the viral contigs identified only 4.5% (n = 7) of total viral contigs were shared between the two genotypes, demonstrating that plants even at the level of genotype can significantly alter the surrounding soil viral community. An auxiliary metabolic gene associated with glycoside hydrolase (GH5) degradation was identified in one viral metagenome-assembled genome (vOTU) identified in the B73 rhizobiome infecting Propionibacteriaceae (Actinobacteriota) further demonstrating the viral contribution in metabolic potential for carbohydrate degradation and carbon cycling in the rhizosphere. This variation demonstrates the potential of plant genotype to contribute to microbial and viral heterogeneity in soil systems and harbors genes capable of contributing to carbon cycling in the rhizosphere.

Prokaryotic virus host prediction with graph contrastive augmentaion

Prokaryotic viruses, also known as bacteriophages, play crucial roles in regulating microbial communities and have the potential for phage therapy applications. Accurate prediction of phage-host interactions is essential for understanding the dynamics of these viruses and their impacts on bacterial populations. Numerous computational methods have been developed to tackle this challenging task. However, most existing prediction models can be constrained due to the substantial number of unknown interactions in comparison to the constrained diversity of available training data. To solve the problem, we introduce a model for prokaryotic virus host prediction with graph contrastive augmentation (PHPGCA). Specifically, we construct a comprehensive heterogeneous graph by integrating virus-virus protein similarity and virus-host DNA sequence similarity information. As the backbone encoder for learning node representations in the virus-prokaryote graph, we employ LGCN, a state-of-the-art graph embedding technique. Additionally, we apply graph contrastive learning to augment the node representations without the need for additional labels. We further conducted two case studies aimed at predicting the host range of multi-species phages, helping to understand the phage ecology and evolution.

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.

Extracellular vesicles are the main contributor to the non-viral protected extracellular sequence space

Environmental virus metagenomes, commonly referred to as "viromes", are typically generated by physically separating virus-like particles (VLPs) from the microbial fraction based on their size and mass. However, most methods used to purify VLPs, enrich extracellular vesicles (EVs) and gene transfer agents (GTAs) simultaneously. Consequently, the sequence space traditionally referred to as a "virome" contains host-associated sequences, transported via EVs or GTAs. We therefore propose to call the genetic material isolated from size-fractionated (0.22 µm) and DNase-treated samples protected environmental DNA (peDNA). This sequence space contains viral genomes, DNA transduced by viruses and DNA transported in EVs and GTAs. Since there is no genetic signature for peDNA transported in EVs, GTAs and virus particles, we rely on the successful removal of contaminating remaining cellular and free DNA when analyzing peDNA. Using marine samples collected from the North Sea, we generated a thoroughly purified peDNA dataset and developed a bioinformatic pipeline to determine the potential origin of the purified DNA. This pipeline was applied to our dataset as well as existing global marine "viromes". Through this pipeline, we identified known GTA and EV producers, as well as organisms with actively transducing proviruses as the source of the peDNA, thus confirming the reliability of our approach. Additionally, we identified novel and widespread EV producers, and found quantitative evidence suggesting that EV-mediated gene transfer plays a significant role in driving horizontal gene transfer (HGT) in the world's oceans.

Phables: from fragmented assemblies to high-quality bacteriophage genomes

MOTIVATION

Microbial communities have a profound impact on both human health and various environments. Viruses infecting bacteria, known as bacteriophages or phages, play a key role in modulating bacterial communities within environments. High-quality phage genome sequences are essential for advancing our understanding of phage biology, enabling comparative genomics studies and developing phage-based diagnostic tools. Most available viral identification tools consider individual sequences to determine whether they are of viral origin. As a result of challenges in viral assembly, fragmentation of genomes can occur, and existing tools may recover incomplete genome fragments. Therefore, the identification and characterization of novel phage genomes remain a challenge, leading to the need of improved approaches for phage genome recovery.

RESULTS

We introduce Phables, a new computational method to resolve phage genomes from fragmented viral metagenome assemblies. Phables identifies phage-like components in the assembly graph, models each component as a flow network, and uses graph algorithms and flow decomposition techniques to identify genomic paths. Experimental results of viral metagenomic samples obtained from different environments show that Phables recovers on average over 49% more high-quality phage genomes compared to existing viral identification tools. Furthermore, Phables can resolve variant phage genomes with over 99% average nucleotide identity, a distinction that existing tools are unable to make.

AVAILABILITY AND IMPLEMENTATION

Phables is available on GitHub at https://github.com/Vini2/phables.

The first arriving virus shapes within-host viral diversity during natural epidemics

Viral diversity has been discovered across scales from host individuals to populations. However, the drivers of viral community assembly are still largely unknown. Within-host viral communities are formed through co-infections, where the interval between the arrival times of viruses may vary. Priority effects describe the timing and order in which species arrive in an environment, and how early colonizers impact subsequent community assembly. To study the effect of the first-arriving virus on subsequent infection patterns of five focal viruses, we set up a field experiment using naïve Plantago lanceolata plants as sentinels during a seasonal virus epidemic. Using joint species distribution modelling, we find both positive and negative effects of early season viral infection on late season viral colonization patterns. The direction of the effect depends on both the host genotype and which virus colonized the host early in the season. It is well established that co-occurring viruses may change the virulence and transmission of viral infections. However, our results show that priority effects may also play an important, previously unquantified role in viral community assembly. The assessment of these temporal dynamics within a community ecological framework will improve our ability to understand and predict viral diversity in natural systems.

Phables: from fragmented assemblies to high-quality bacteriophage genomes

Microbial communities influence both human health and different environments. Viruses infecting bacteria, known as bacteriophages or phages, play a key role in modulating bacterial communities within environments. High-quality phage genome sequences are essential for advancing our understanding of phage biology, enabling comparative genomics studies, and developing phage-based diagnostic tools. Most available viral identification tools consider individual sequences to determine whether they are of viral origin. As a result of the challenges in viral assembly, fragmentation of genomes can occur, leading to the need for new approaches in viral identification. Therefore, the identification and characterisation of novel phages remain a challenge. We introduce Phables, a new computational method to resolve phage genomes from fragmented viral metagenome assemblies. Phables identifies phage-like components in the assembly graph, models each component as a flow network, and uses graph algorithms and flow decomposition techniques to identify genomic paths. Experimental results of viral metagenomic samples obtained from different environments show that Phables recovers on average over 49% more high-quality phage genomes compared to existing viral identification tools. Furthermore, Phables can resolve variant phage genomes with over 99% average nucleotide identity, a distinction that existing tools are unable to make. Phables is available on GitHub at https://github.com/Vini2/phables.

Macromonas nakdongensis sp. nov., Isolated from Freshwater and Characterization of Bacteriophage BK-30P-The First Phage That Infects Genus Macromonas

A Gram-stain-negative, non-motile, non-pigmented, rod-shaped bacterium was isolated from a freshwater sample of Nakdong River in South Korea and designated as strain BK-30T. An analysis of the 16S rRNA gene sequence of strain BK-30T revealed its closest phylogenetic neighbors were members of the genus Macromonas. Specifically, the strain formed a robust clade with Macromonas bipunctata DSM 12705T, sharing 98.4% similarity in their 16S rRNA gene sequences. The average nucleotide identity value between strain BK-30T and M. bipunctata DSM 12705T was 79.8%, and the genome-to-genome distance averaged 21.3%, indicating the representation of a novel genomic species. Strain BK-30T exhibited optimum growth at 30 °C and pH 7.0, in the absence of NaCl. The major respiratory isoprenoid quinone identified was ubiquinone-8 (Q-8). The principal fatty acids detected were C16:1ω7c and/or C16:1ω6c (49.6%), C16:0 (27.5%), and C18:1ω7c and/or C18:1 ω6c (9.2%). The DNA G+C content of the strain was determined to be 67.3 mol%. Based on these data, we propose a novel species within the genus Macromonas, named Macromonas nakdongensis sp. nov., to accommodate the bacterial isolate. Strain BK-30T is designated as the type strain (=KCTC 52161T = JCM 31376T = FBCC-B1). Additionally, we present the isolation and complete genome sequence of a lytic phage infecting strain BK-30T, named BK-30P. This bacteriophage is the first reported to infect Macromonas, leading us to propose the name "Macromonasphage".

TRGdb: a universal resource for the exploration of taxonomically restricted genes in bacteria

The TRGdb database is a resource dedicated to taxonomically restricted genes (TRGs) in bacteria. It provides a comprehensive collection of genes that are specific to different genera and species, according to the latest release of bacterial taxonomy. The user interface allows for easy browsing and searching as well as sequence similarity exploration. The website also provides information on each TRG protein sequence, including its level of disorder, complexity and tendency to aggregate. TRGdb is a valuable resource for gaining a deeper understanding of the TRG-associated, unique features, and characteristics of bacterial organisms. Database URL www.combio.pl/trgdb.

Everything OLD is new again: How structural, functional, and bioinformatic advances have redefined a neglected nuclease family

Overcoming lysogenization defect (OLD) proteins are a conserved family of ATP-powered nucleases that function in anti-phage defense. Recent bioinformatic, genetic, and crystallographic studies have yielded new insights into the structure, function, and evolution of these enzymes. Here we review these developments and propose a new classification scheme to categorize OLD homologs that relies on gene neighborhoods, biochemical properties, domain organization, and catalytic machinery. This taxonomy reveals important similarities and differences between family members and provides a blueprint to contextualize future in vivo and in vitro findings. We also detail how OLD nucleases are related to PARIS and Septu anti-phage defense systems and discuss important mechanistic questions that remain unanswered.

Metavirome Profiling and Dynamics of the DNA Viral Community in Seawater in Chuuk State, Federated States of Micronesia

Despite their abundance and ecological importance, little is known about the diversity of marine viruses, in part because most cannot be cultured in the laboratory. Here, we used high-throughput viral metagenomics of uncultivated viruses to investigate the dynamics of DNA viruses in tropical seawater sampled from Chuuk State, Federated States of Micronesia, in March, June, and December 2014. Among the identified viruses, 71-79% were bacteriophages belonging to the families Myoviridae, Siphoviridae, and Podoviridae (Caudoviriales), listed in order of abundance at all sampling times. Although the measured environmental factors (temperature, salinity, and pH) remained unchanged in the seawater over time, viral dynamics changed. The proportion of cyanophages (34.7%) was highest in June, whereas the proportion of mimiviruses, phycodnaviruses, and other nucleo-cytoplasmic large DNA viruses (NCLDVs) was higher in March and December. Although host species were not analysed, the dramatic viral community change observed in June was likely due to changes in the abundance of cyanophage-infected cyanobacteria, whereas that in NCLDVs was likely due to the abundance of potential eukaryote-infected hosts. These results serve as a basis for comparative analyses of other marine viral communities, and guide policy-making when considering marine life care in Chuuk State.

Metagenomics of gut microbiome for migratory seagulls in Kunming city revealed the potential public risk to human health

BACKGROUND

Seagull as a migratory wild bird has become most popular species in southwest China since 1980s. Previously, we analyzed the gut microbiota and intestinal pathogenic bacteria configuration for this species by using 16S rRNA sequencing and culture methods. To continue in-depth research on the gut microbiome of migratory seagulls, the metagenomics, DNA virome and RNA virome were both investigated for their gut microbial communities of abundance and diversity in this study.

RESULTS

The metagenomics results showed 99.72% of total species was bacteria, followed by viruses, fungi, archaea and eukaryota. In particular, Shigella sonnei, Escherichia albertii, Klebsiella pneumonia, Salmonella enterica and Shigella flexneri were the top distributed taxa at species level. PCoA, NMDS, and statistics indicated some drug resistant genes, such as adeL, evgS, tetA, PmrF, and evgA accumulated as time went by from November to January of the next year, and most of these genes were antibiotic efflux. DNA virome composition demonstrated that Caudovirales was the most abundance virus, followed by Cirlivirales, Geplafuvirales, Petitvirales and Piccovirales. Most of these phages corresponded to Enterobacteriaceae and Campylobacteriaceae bacterial hosts respectively. Caliciviridae, Coronaviridae and Picornaviridae were the top distributed RNA virome at family level of this migratory animal. Phylogenetic analysis indicated the sequences of contigs of Gammacoronavirus and Deltacoronavirus had highly similarity with some coronavirus references.

CONCLUSIONS

In general, the characteristics of gut microbiome of migratory seagulls were closely related to human activities, and multiomics still revealed the potential public risk to human health.

Ecological drivers and potential functions of viral communities in flooded arsenic-contaminated paddy soils

This work revealed the profile of viral communities in paddy soils with different levels of arsenic (As) contamination during the flooded period. The structure of viral communities differed significantly in highly and moderately As-contaminated soils. The diversity of soil viral communities under high As contamination decreased. Siphoviridae, Podoviridae, Myoviridae, and Microviridae were the dominant viral families in all samples, and the relative abundances of five of the top 20 viral genera were significantly different between highly and moderately As-contaminated groups. Seventeen dissimilatory As(V)-reducing bacteria were predicted to host 161 viral operational taxonomic units (vOTUs), mainly affiliated with the genera of Sulfurospirillum, Deferribacter, Bacillus and Fusibacter. Among them, 28 vOTUs were also associated with Fe(III)-reducing bacteria, which belonged to different species of the genus Shewanella. Procrustes analysis showed that the community structure of soil viruses was strongly correlated with both prokaryotic community structure and geochemical properties. Random forest analyses revealed that the Total-Fe, DCB-Fe and oxalate-Fe were the most significant variables on viral community richness, while the total-As concentration was an important factor on the Shannon index. Furthermore, As resistance genes (ArsC, ArsR and ArsD), As methylation genes (arsM) and As transporter genes (Pst and Pit) were identified among the auxiliary metabolic genes (AMGs) of the virome. This work revealed that the viruses might influence microbial adaptation in response to As-induced stress, and provided a perspective on the potential virus-mediated biogeochemical cycling of As.

Occurrences of similar viral diversity in campus wastewater and reclaimed water of a university dormitory

Water reuse from wastewater sources still remain some critical safety concerns associated with treacherous contaminants like pathogenic viruses. In this study, viral diversities in campus wastewater (CWW) and its reclaimed water (RCW) recycled for toilet flushing and garden irrigation of a university dormitory were assessed using metagenomic sequencing for acquisition of more background data. Results suggested majority (>80%) of gene sequences within assembled contigs predicted by open reading frame (ORF) finder were no-hit yet believed to be novel/unrevealed viral genomic information whereas hits matched bacteriophages (i.e., mainly Myoviridae, Podoviridae, and Siphoviridae families) were predominant in both CWW and RCW samples. Moreover, few pathogenic viruses (<1%) related to infections of human skin (e.g., Molluscum contagiosum virus, MCV), digestion system (e.g., hepatitis C virus, HCV), and gastrointestinal tract (e.g., human norovirus, HuNoV) were also noticed raising safety concerns about application of reclaimed waters. Low-affinity interactions of particular viral exterior proteins (e.g., envelope glycoproteins or spike proteins) for disinfectant ligand (e.g., chlorite) elucidated treatment limitations of current sewage processing systems even with membrane bioreactor and disinfectant contactor. Revolutionary disinfection approaches together with routine monitoring and new regulations are prerequisite to secure pathogen-correlated water quality for safer reuse of reclaimed waters.

The Promise and Pitfalls of Prophages

Phages dominate every ecosystem on the planet. While virulent phages sculpt the microbiome by killing their bacterial hosts, temperate phages provide unique growth advantages to their hosts through lysogenic conversion. Many prophages benefit their host, and prophages are responsible for genotypic and phenotypic differences that separate individual microbial strains. However, the microbes also endure a cost to maintain those phages: additional DNA to replicate and proteins to transcribe and translate. We have never quantified those benefits and costs. Here, we analysed over two and a half million prophages from over half a million bacterial genome assemblies. Analysis of the whole dataset and a representative subset of taxonomically diverse bacterial genomes demonstrated that the normalised prophage density was uniform across all bacterial genomes above 2 Mbp. We identified a constant carrying capacity of phage DNA per bacterial DNA. We estimated that each prophage provides cellular services equivalent to approximately 2.4 % of the cell's energy or 0.9 ATP per bp per hour. We demonstrate analytical, taxonomic, geographic, and temporal disparities in identifying prophages in bacterial genomes that provide novel targets for identifying new phages. We anticipate that the benefits bacteria accrue from the presence of prophages balance the energetics involved in supporting prophages. Furthermore, our data will provide a new framework for identifying phages in environmental datasets, diverse bacterial phyla, and from different locations.

Discovery of phage determinants that confer sensitivity to bacterial immune systems

Over the past few years, numerous anti-phage defense systems have been discovered in bacteria. Although the mechanism of defense for some of these systems is understood, a major unanswered question is how these systems sense phage infection. To systematically address this question, we isolated 177 phage mutants that escape 15 different defense systems. In many cases, these escaper phages were mutated in the gene sensed by the defense system, enabling us to map the phage determinants that confer sensitivity to bacterial immunity. Our data identify specificity determinants of diverse retron systems and reveal phage-encoded triggers for multiple abortive infection systems. We find general themes in phage sensing and demonstrate that mechanistically diverse systems have converged to sense either the core replication machinery of the phage, phage structural components, or host takeover mechanisms. Combining our data with previous findings, we formulate key principles on how bacterial immune systems sense phage invaders.

Unlocking the microbial studies through computational approaches: how far have we reached?

The metagenomics approach accelerated the study of genetic information from uncultured microbes and complex microbial communities. In silico research also facilitated an understanding of protein-DNA interactions, protein-protein interactions, docking between proteins and phyto/biochemicals for drug design, and modeling of the 3D structure of proteins. These in silico approaches provided insight into analyzing pathogenic and nonpathogenic strains that helped in the identification of probable genes for vaccines and antimicrobial agents and comparing whole-genome sequences to microbial evolution. Artificial intelligence, more precisely machine learning (ML) and deep learning (DL), has proven to be a promising approach in the field of microbiology to handle, analyze, and utilize large data that are generated through nucleic acid sequencing and proteomics. This enabled the understanding of the functional and taxonomic diversity of microorganisms. ML and DL have been used in the prediction and forecasting of diseases and applied to trace environmental contaminants and environmental quality. This review presents an in-depth analysis of the recent application of silico approaches in microbial genomics, proteomics, functional diversity, vaccine development, and drug design.

First Discovery of Phenuiviruses within Diverse RNA Viromes of Asiatic Toad (Bufo gargarizans) by Metagenomics Sequencing

Most zoonotic pathogens originate from mammals and avians, but viral diversity and related biosafety risk assessment in lower vertebrates also need to be explored. Amphibians are an important group of lower vertebrates that played a momentous role in animal evolution. To elucidate the diversity of RNA viruses in one important species of amphibians, the Asiatic toad (Bufo gargarizans), we obtained 44 samples including lung, gut, liver, and kidney tissues from Asiatic toads in Sichuan and Jilin provinces, China, for viral metagenomics sequencing. More than 20 novel RNA viruses derived from the order Bunyavirales and 7 families of Astroviridae, Dicistroviridae, Leviviridae, Partitiviridae, Picornaviridae, Rhabdoviridae, and Virgaviridae were discovered, which were distinct from previously described viruses and formed new clusters, as revealed by phylogenetic analyses. Notably, a novel bastrovirus, AtBastV/GCCDC11/2022, of the family Astroviridae was identified from the gut library, the genome of which contains three open reading frames, with the RNA-dependent RNA polymerase (RdRp) coded by ORF1 closely related to that of hepeviruses, and ORF2 encoding an astrovirus-related capsid protein. Notably, phenuiviruses were discovered for the first time in amphibians. AtPhenV1/GCCDC12/2022 and AtPhenV2/GCCDC13/2022 clustered together and formed a clade with the group of phenuiviruses identified from rodents. Picornaviruses and several invertebrate RNA viruses were also detected. These findings improve our understanding of the high RNA viral diversity in the Asiatic toad and provide new insights in the evolution of RNA viruses in amphibians.

Assembly and Annotation of Viral Metagenomes from Short-Read Sequencing Data

Viral metagenomics enables the detection, characterization, and quantification of viral sequences present in shotgun-sequenced datasets of purified virus-like particles and whole metagenomes. Next generation sequencing (Illumina) derived short single or paired-end read runs are a principal platform for metagenomics, and assembly of short reads allows for the identification of distinguishing viral signatures and complex genomic features for taxonomy and functional annotation. Here we describe the identification and characterization of viral genome sequences, bacteriophages, and eukaryotic viruses, from a cohort of human stool samples, using multiple methods. Following the purification of virus-like particles, sequencing, quality refinement, and genome assembly, we begin the protocol with raw short reads deposited in an open-source nucleotide archive. We highlight the use of VIBRANT, an automated computational tool for the characterization of microbial viruses and their viral community function. Finally, we also describe an alternative assembly-free option of mapping reads to established databases of reference genomes and previously characterized metagenome-assembled viral genomes.

Virus diversity and interactions with hosts in deep-sea hydrothermal vents

BACKGROUND

The deep sea harbors many viruses, yet their diversity and interactions with hosts in hydrothermal ecosystems are largely unknown. Here, we analyzed the viral composition, distribution, host preference, and metabolic potential in different habitats of global hydrothermal vents, including vent plumes, background seawater, diffuse fluids, and sediments.

RESULTS

From 34 samples collected at eight vent sites, a total of 4662 viral populations (vOTUs) were recovered from the metagenome assemblies, encompassing diverse phylogenetic groups and defining many novel lineages. Apart from the abundant unclassified viruses, tailed phages are most predominant across the global hydrothermal vents, while single-stranded DNA viruses, including Microviridae and small eukaryotic viruses, also constitute a significant part of the viromes. As revealed by protein-sharing network analysis, hydrothermal vent viruses formed many novel genus-level viral clusters and are highly endemic to specific vent sites and habitat types. Only 11% of the vOTUs can be linked to hosts, which are the key microbial taxa of hydrothermal habitats, such as Gammaproteobacteria and Campylobacterota. Intriguingly, vent viromes share some common metabolic features in that they encode auxiliary genes that are extensively involved in the metabolism of carbohydrates, amino acids, cofactors, and vitamins. Specifically, in plume viruses, various auxiliary genes related to methane, nitrogen, and sulfur metabolism were observed, indicating their contribution to host energy conservation. Moreover, the prevalence of sulfur-relay pathway genes indicated the significant role of vent viruses in stabilizing the tRNA structure, which promotes host adaptation to steep environmental gradients.

CONCLUSIONS

The deep-sea hydrothermal systems hold untapped viral diversity with novelty. They may affect both vent prokaryotic and eukaryotic communities and modulate host metabolism related to vent adaptability. More explorations are needed to depict global vent virus diversity and its roles in this unique ecosystem. Video Abstract.

Over two decades of research on the marine RNA virosphere

RNA viruses (realm: Riboviria), including RNA phages and eukaryote-infecting RNA viruses, are essential components of marine ecosystems. A large number of marine RNA viruses have been discovered in the last two decades because of the rapid development of next-generation sequencing (NGS) technology. Indeed, the combination of NGS and state-of-the-art meta-omics methods (viromics, the study of all viruses in a specific environment) has led to a fundamental understanding of the taxonomy and genetic diversity of RNA viruses in the sea, suggesting the complex ecological roles played by RNA viruses in this complex ecosystem. Furthermore, comparisons of viromes in the context of highly variable marine niches reveal the biogeographic patterns and ecological impact of marine RNA viruses, whose role in global ecology is becoming increasingly clearer. In this review, we summarize the characteristics of the global marine RNA virosphere and outline the taxonomic hierarchy of RNA viruses with a specific focus on their ancient evolutionary history. We also review the development of methodology and the major progress resulting from its applications in RNA viromics. The aim of this review is not only to provide an in-depth understanding of multifaceted aspects of marine RNA viruses, but to offer future perspectives on developing a better methodology for discovery, and exploring the evolutionary origin and major ecological significance of marine RNA virosphere.

Diversity in the soil virosphere: to infinity and beyond?

Viruses are key members of Earth's microbiomes, shaping microbial community composition and metabolism. Here, we describe recent advances in 'soil viromics', that is, virus-focused metagenome and metatranscriptome analyses that offer unprecedented windows into the soil virosphere. Given the emerging picture of high soil viral activity, diversity, and dynamics over short spatiotemporal scales, we then outline key eco-evolutionary processes that we hypothesize are the major diversity drivers for soil viruses. We argue that a community effort is needed to establish a 'global soil virosphere atlas' that can be used to address the roles of viruses in soil microbiomes and terrestrial biogeochemical cycles across spatiotemporal scales.

Highly diverse ribonucleic acid viruses in the viromes of eukaryotic host species in Yunnan province, China

Background

The diversity in currently documented viruses and their morphological characteristics indicates the need for understanding the evolutionary characteristics of viruses. Notably, further studies are needed to obtain a comprehensive landscape of virome, the virome of host species in Yunnan province, China.

Materials and methods

We implemented the metagenomic next-generation sequencing strategy to investigate the viral diversity, which involved in 465 specimens collected from bats, pangolins, monkeys, and other species. The diverse RNA viruses were analyzed, especially focusing on the genome organization, genetic divergence and phylogenetic relationships.

Results

In this study, we investigated the viral composition of eight libraries from bats, pangolins, monkeys, and other species, and found several diverse RNA viruses, including the Alphacoronavirus from bat specimens. By characterizing the genome organization, genetic divergence, and phylogenetic relationships, we identified five Alphacoronavirus strains, which shared phylogenetic association with Bat-CoV-HKU8-related strains. The pestivirus-like virus related to recently identified Dongyang pangolin virus (DYPV) strains from dead pangolin specimens, suggesting that these viruses are evolving. Some genomes showed higher divergence from known species (e.g., calicivirus CS9-Cali-YN-CHN-2020), and many showed evidence of recombination events with unknown or known strains (e.g., mamastroviruses BF2-astro-YN-CHN-2020 and EV-A122 AKM5-YN-CHN-2020). The newly identified viruses showed extensive changes and could be assigned as new species, or even genus (e.g., calicivirus CS9-Cali-YN-CHN-2020 and iflavirus Ifla-YN-CHN-2020). Moreover, we identified several highly divergent RNA viruses and estimated their evolutionary characteristics among different hosts, providing data for further examination of their evolutionary dynamics.

Conclusion

Overall, our study emphasizes the close association between emerging viruses and infectious diseases, and the need for more comprehensive surveys.

Insights into the global freshwater virome

Viruses are by far the most abundant life forms on this planet. Yet, the full viral diversity remains mostly unknown, especially in environments like freshwater. Therefore, we aimed to study freshwater viruses in a global context. To this end, we downloaded 380 publicly available viral metagenomes (>1 TB). More than 60% of these metagenomes were discarded based on their levels of cellular contamination assessed by ribosomal DNA content. For the remaining metagenomes, assembled contigs were decontaminated using two consecutive steps, eventually yielding 273,365 viral contigs longer than 1,000 bp. Long enough contigs (≥ 10 kb) were clustered to identify novel genomes/genome fragments. We could recover 549 complete circular and high-quality draft genomes, out of which 10 were recognized as being novel. Functional annotation of these genomes showed that most of the annotated coding sequences are DNA metabolic genes or phage structural genes. On the other hand, taxonomic analysis of viral contigs showed that most of the assigned contigs belonged to the order Caudovirales, particularly the families of Siphoviridae, Myoviridae, and Podoviridae. The recovered viral contigs contained several auxiliary metabolic genes belonging to several metabolic pathways, especially carbohydrate and amino acid metabolism in addition to photosynthesis as well as hydrocarbon degradation and antibiotic resistance. Overall, we present here a set of prudently chosen viral contigs, which should not only help better understanding of freshwater viruses but also be a valuable resource for future virome studies.

A comparative study of flow cytometry-sorted communities and shotgun viral metagenomics in a Singapore municipal wastewater treatment plant

Traditional or "bulk" viral enrichment and amplification methods used in viral metagenomics introduce unavoidable bias in viral diversity. This bias is due to shortcomings in existing viral enrichment methods and overshadowing by the more abundant viral populations. To reduce the complexity and improve the resolution of viral diversity, we developed a strategy coupling fluorescence-activated cell sorting (FACS) with random amplification and compared this to bulk metagenomics. This strategy was validated on both influent and effluent samples from a municipal wastewater treatment plant using the Modified Ludzack-Ettinger (MLE) process as the treatment method. We found that DNA and RNA communities generated using bulk samples were mostly different from those derived following FACS for both treatments before and after MLE. Before MLE treatment, FACS identified five viral families and 512 viral annotated contigs. Up to 43% of mapped reads were not detected in bulk samples. Nucleo-cytoplasmic large DNA viral families were enriched to a greater extent in the FACS-coupled subpopulations compared with bulk samples. FACS-coupled viromes captured a single-contig viral genome associated with Anabaena phage, which was not observed in bulk samples or in FACS-sorted samples after MLE. These short metagenomic reads, which were assembled into a high-quality draft genome of 46 kbp, were found to be highly dominant in one of the pre-MLE FACS annotated virome fractions (57.4%). Using bulk metagenomics, we identified that between Primary Settling Tank and Secondary Settling Tank viromes, Virgaviridae, Astroviridae, Parvoviridae, Picobirnaviridae, Nodaviridae, and Iridoviridae were susceptible to MLE treatment. In all, bulk and FACS-coupled metagenomics are complementary approaches that enable a more thorough understanding of the community structure of DNA and RNA viruses in complex environmental samples, of which the latter is critical for increasing the sensitivity of detection of viral signatures that would otherwise be lost through bulk viral metagenomics.

Viral Infection in Esophageal, Gastric, and Colorectal Cancer

The human gastrointestinal tract, which constitutes the digestive system, contains a large number of virus particles that maintain organizational homeostasis and health. Conversely, viral pathogens have also attracted attention for their involvement in the pathogenesis of certain cancers, including gastrointestinal cancers. To aid prevention and treatment of these cancers, the relevance of gastrointestinal viral factors as potential risk factors needs to be carefully investigated. This review summarizes and discusses the available literature on the relationship between the development of esophageal, gastric, and colorectal cancers and their corresponding viruses. This review reveals that research on the association between colorectal cancer and viruses, in particular, is still in its infancy compared to the association between HPV and esophageal cancer and between EBV and gastric cancer.

The complexity landscape of viral genomes

BACKGROUND

Viruses are among the shortest yet highly abundant species that harbor minimal instructions to infect cells, adapt, multiply, and exist. However, with the current substantial availability of viral genome sequences, the scientific repertory lacks a complexity landscape that automatically enlights viral genomes' organization, relation, and fundamental characteristics.

RESULTS

This work provides a comprehensive landscape of the viral genome's complexity (or quantity of information), identifying the most redundant and complex groups regarding their genome sequence while providing their distribution and characteristics at a large and local scale. Moreover, we identify and quantify inverted repeats abundance in viral genomes. For this purpose, we measure the sequence complexity of each available viral genome using data compression, demonstrating that adequate data compressors can efficiently quantify the complexity of viral genome sequences, including subsequences better represented by algorithmic sources (e.g., inverted repeats). Using a state-of-the-art genomic compressor on an extensive viral genomes database, we show that double-stranded DNA viruses are, on average, the most redundant viruses while single-stranded DNA viruses are the least. Contrarily, double-stranded RNA viruses show a lower redundancy relative to single-stranded RNA. Furthermore, we extend the ability of data compressors to quantify local complexity (or information content) in viral genomes using complexity profiles, unprecedently providing a direct complexity analysis of human herpesviruses. We also conceive a features-based classification methodology that can accurately distinguish viral genomes at different taxonomic levels without direct comparisons between sequences. This methodology combines data compression with simple measures such as GC-content percentage and sequence length, followed by machine learning classifiers.

CONCLUSIONS

This article presents methodologies and findings that are highly relevant for understanding the patterns of similarity and singularity between viral groups, opening new frontiers for studying viral genomes' organization while depicting the complexity trends and classification components of these genomes at different taxonomic levels. The whole study is supported by an extensive website (https://asilab.github.io/canvas/) for comprehending the viral genome characterization using dynamic and interactive approaches.

DeePVP: Identification and classification of phage virion proteins using deep learning

BACKGROUND

Many biological properties of phages are determined by phage virion proteins (PVPs), and the poor annotation of PVPs is a bottleneck for many areas of viral research, such as viral phylogenetic analysis, viral host identification, and antibacterial drug design. Because of the high diversity of PVP sequences, the PVP annotation of a phage genome remains a particularly challenging bioinformatic task.

FINDINGS

Based on deep learning, we developed DeePVP. The main module of DeePVP aims to discriminate PVPs from non-PVPs within a phage genome, while the extended module of DeePVP can further classify predicted PVPs into the 10 major classes of PVPs. Compared with the present state-of-the-art tools, the main module of DeePVP performs better, with a 9.05% higher F1-score in the PVP identification task. Moreover, the overall accuracy of the extended module of DeePVP in the PVP classification task is approximately 3.72% higher than that of PhANNs. Two application cases show that the predictions of DeePVP are more reliable and can better reveal the compact PVP-enriched region than the current state-of-the-art tools. Particularly, in the Escherichia phage phiEC1 genome, a novel PVP-enriched region that is conserved in many other Escherichia phage genomes was identified, indicating that DeePVP will be a useful tool for the analysis of phage genomic structures.

CONCLUSIONS

DeePVP outperforms state-of-the-art tools. The program is optimized in both a virtual machine with graphical user interface and a docker so that the tool can be easily run by noncomputer professionals. DeePVP is freely available at https://github.com/fangzcbio/DeePVP/.

Identification of bacteriophage genome sequences with representation learning

MOTIVATION

Bacteriophages/phages are the viruses that infect and replicate within bacteria and archaea, and rich in human body. To investigate the relationship between phages and microbial communities, the identification of phages from metagenome sequences is the first step. Currently, there are two main methods for identifying phages: database-based (alignment-based) methods and alignment-free methods. Database-based methods typically use a large number of sequences as references; alignment-free methods usually learn the features of the sequences with machine learning and deep learning models.

RESULTS

We propose INHERIT which uses a deep representation learning model to integrate both database-based and alignment-free methods, combining the strengths of both. Pre-training is used as an alternative way of acquiring knowledge representations from existing databases, while the BERT-style deep learning framework retains the advantage of alignment-free methods. We compare INHERIT with four existing methods on a third-party benchmark dataset. Our experiments show that INHERIT achieves a better performance with the F1-score of 0.9932. In addition, we find that pre-training two species separately helps the non-alignment deep learning model make more accurate predictions.

AVAILABILITY AND IMPLEMENTATION

The codes of INHERIT are now available in: https://github.com/Celestial-Bai/INHERIT.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Molecular Characterization of a DNA Polymerase from Thermus thermophilus MAT72 Phage vB_Tt72: A Novel Type-A Family Enzyme with Strong Proofreading Activity

We present a structural and functional analysis of the DNA polymerase of thermophilic Thermus thermophilus MAT72 phage vB_Tt72. The enzyme shows low sequence identity (<30%) to the members of the type-A family of DNA polymerases, except for two yet uncharacterized DNA polymerases of T. thermophilus phages: φYS40 (91%) and φTMA (90%). The Tt72 polA gene does not complement the Escherichia colipolA− mutant in replicating polA-dependent plasmid replicons. It encodes a 703-aa protein with a predicted molecular weight of 80,490 and an isoelectric point of 5.49. The enzyme contains a nucleotidyltransferase domain and a 3′-5′ exonuclease domain that is engaged in proofreading. Recombinant enzyme with His-tag at the N-terminus was overproduced in E. coli, subsequently purified by immobilized metal affinity chromatography, and biochemically characterized. The enzyme exists in solution in monomeric form and shows optimum activity at pH 8.5, 25 mM KCl, and 0.5 mM Mg2+. Site-directed analysis proved that highly-conserved residues D15, E17, D78, D180, and D184 in 3′-5′ exonuclease and D384 and D615 in the nucleotidyltransferase domain are critical for the enzyme’s activity. Despite the source of origin, the Tt72 DNA polymerase has not proven to be highly thermoresistant, with a temperature optimum at 55 °C. Above 60 °C, the rapid loss of function follows with no activity > 75 °C. However, during heat treatment (10 min at 75 °C), trehalose, trimethylamine N-oxide, and betaine protected the enzyme against thermal inactivation. A midpoint of thermal denaturation at Tm = 74.6 °C (ΔHcal = 2.05 × 104 cal mol−1) and circular dichroism spectra > 60 °C indicate the enzyme’s moderate thermal stability.

Ecosystem-specific microbiota and microbiome databases in the era of big data

The rapid development of sequencing methods over the past decades has accelerated both the potential scope and depth of microbiota and microbiome studies. Recent developments in the field have been marked by an expansion away from purely categorical studies towards a greater investigation of community functionality. As in-depth genomic and environmental coverage is often distributed unequally across major taxa and ecosystems, it can be difficult to identify or substantiate relationships within microbial communities. Generic databases containing datasets from diverse ecosystems have opened a new era of data accessibility despite costs in terms of data quality and heterogeneity. This challenge is readily embodied in the integration of meta-omics data alongside habitat-specific standards which help contextualise datasets both in terms of sample processing and background within the ecosystem. A special case of large genomic repositories, ecosystem-specific databases (ES-DB's), have emerged to consolidate and better standardise sample processing and analysis protocols around individual ecosystems under study, allowing independent studies to produce comparable datasets. Here, we provide a comprehensive review of this emerging tool for microbial community analysis in relation to current trends in the field. We focus on the factors leading to the formation of ES-DB's, their comparison to traditional microbial databases, the potential for ES-DB integration with meta-omics platforms, as well as inherent limitations in the applicability of ES-DB's.

The simple emergence of complex molecular function

At odds with a traditional view of molecular evolution that seeks a descent-with-modification relationship between functional sequences, new functions can emerge de novo with relative ease. At early times of molecular evolution, random polymers could have sufficed for the appearance of incipient chemical activity, while the cellular environment harbours a myriad of proto-functional molecules. The emergence of function is facilitated by several mechanisms intrinsic to molecular organization, such as redundant mapping of sequences into structures, phenotypic plasticity, modularity or cooperative associations between genomic sequences. It is the availability of niches in the molecular ecology that filters new potentially functional proposals. New phenotypes and subsequent levels of molecular complexity could be attained through combinatorial explorations of currently available molecular variants. Natural selection does the rest. This article is part of the theme issue 'Emergent phenomena in complex physical and socio-technical systems: from cells to societies'.