IMG/VR v4: an expanded database of uncultivated virus genomes within a framework of extensive functional, taxonomic, and ecological metadata

Visualizing metagenomic and metatranscriptomic data: A comprehensive review

The fields of Metagenomics and Metatranscriptomics involve the examination of complete nucleotide sequences, gene identification, and analysis of potential biological functions within diverse organisms or environmental samples. Despite the vast opportunities for discovery in metagenomics, the sheer volume and complexity of sequence data often present challenges in processing analysis and visualization. This article highlights the critical role of advanced visualization tools in enabling effective exploration, querying, and analysis of these complex datasets. Emphasizing the importance of accessibility, the article categorizes various visualizers based on their intended applications and highlights their utility in empowering bioinformaticians and non-bioinformaticians to interpret and derive insights from meta-omics data effectively.

Exploring the roles of ribosomal peptides in prokaryote-phage interactions through deep learning-enabled metagenome mining

BACKGROUND

Microbial secondary metabolites play a crucial role in the intricate interactions within the natural environment. Among these metabolites, ribosomally synthesized and post-translationally modified peptides (RiPPs) are becoming a promising source of therapeutic agents due to their structural diversity and functional versatility. However, their biosynthetic capacity and ecological functions remain largely underexplored.

RESULTS

Here, we aim to explore the biosynthetic profile of RiPPs and their potential roles in the interactions between microbes and viruses in the ocean, which encompasses a vast diversity of unique biomes that are rich in interactions and remains chemically underexplored. We first developed TrRiPP to identify RiPPs from ocean metagenomes, a deep learning method that detects RiPP precursors in a hallmark gene-independent manner to overcome the limitations of classic methods in processing highly fragmented metagenomic data. Applying this method to metagenomes from the global ocean microbiome, we uncover a diverse array of previously uncharacterized putative RiPP families with great novelty and diversity. Through correlation analysis based on metatranscriptomic data, we observed a high prevalence of antiphage defense-related and phage-related protein families that were co-expressed with RiPP families. Based on this putative association between RiPPs and phage infection, we constructed an Ocean Virus Database (OVD) and established a RiPP-involving host-phage interaction network through host prediction and co-expression analysis, revealing complex connectivities linking RiPP-encoding prokaryotes, RiPP families, viral protein families, and phages. These findings highlight the potential of RiPP families involved in prokaryote-phage interactions and coevolution, providing insights into their ecological functions in the ocean microbiome.

CONCLUSIONS

This study provides a systematic investigation of the biosynthetic potential of RiPPs from the ocean microbiome at a global scale, shedding light on the essential insights into the ecological functions of RiPPs in prokaryote-phage interactions through the integration of deep learning approaches, metatranscriptomic data, and host-phage connectivity. This study serves as a valuable example of exploring the ecological functions of bacterial secondary metabolites, particularly their associations with unexplored microbial interactions. Video Abstract.

Profiles of phage in global hospital wastewater: Association with microbial hosts, antibiotic resistance genes, metal resistance genes, and mobile genetic elements

Hospital wastewater (HWW) is known to host taxonomically diverse microbial communities, yet limited information is available on the phages infecting these microorganisms. To fill this knowledge gap, we conducted an in-depth analysis using 377 publicly available HWW metagenomic datasets from 16 countries across 4 continents in the NCBI SRA database to elucidate phage-host dynamics and phage contributions to resistance gene transmission. We first assembled a metagenomic HWW phage catalog comprising 13,812 phage operational taxonomic units (pOTUs). The majority of these pOTUs belonged to the Caudoviricetes order, representing 75.29 % of this catalog. Based on the lifestyle of phages, we found that potentially virulent phages predominated in HWW. Specifically, 583 pOTUs have been predicted to have the capability to lyse 81 potentially pathogenic bacteria, suggesting the promising role of HWW phages as a viable alternative to antibiotics. Among all pOTUs, 1.56 % of pOTUs carry 108 subtypes of antibiotic resistance genes (ARGs), 0.96 % of pOTUs carry 76 subtypes of metal resistance genes (MRGs), and 0.96 % of pOTUs carry 22 subtypes of non-phage mobile genetic elements (MGEs). Predictions indicate that certain phages carrying ARGs, MRGs, and non-phage MGEs could infect bacteria hosts, even potential pathogens. This suggests that phages in HWW may contribute to the dissemination of resistance-associated genes in the environment. This meta-analysis provides the first global catalog of HWW phages, revealing their correlations with microbial hosts and pahge-associated ARGs, MRG, and non-phage MGEs. The insights gained from this research hold promise for advancing the applications of phages in medical and industrial contexts.

Microbial adaptation to spaceflight is correlated with bacteriophage-encoded functions

Evidence from the International Space Station suggests microbial populations are rapidly adapting to the spacecraft environment; however, the mechanism of this adaptation is not understood. Bacteriophages are prolific mediators of bacterial adaptation on Earth. Here we survey 245 genomes sequenced from bacterial strains isolated on the International Space Station for dormant (lysogenic) bacteriophages. Our analysis indicates phage-associated genes are significantly different between spaceflight strains and their terrestrial counterparts. In addition, we identify 283 complete prophages, those that could initiate bacterial lysis and infect additional hosts, of which 21% are novel. These prophage regions encode functions that correlate with increased persistence in extreme environments, such as spaceflight, to include antimicrobial resistance and virulence, DNA damage repair, and dormancy. Our results correlate microbial adaptation in spaceflight to bacteriophage-encoded functions that may impact human health in spaceflight.

The bacteriome-coupled phage communities continuously contract and shift to orchestrate the traditional rice vinegar fermentation

Amounts of microbiome studies have uncovered the microbial communities of traditional food fermentations, while in which the phageome development with time is poorly understood. Here, we conducted a study to decipher both phageome and bacteriome of the traditional rice vinegar fermentation. The vinegar phageomes showed significant differences in the alpha diversity, network density and clustering coefficient over time. Peduoviridae had the highest relative abundance. Moreover, the phageome negatively correlated to the cognate bacteriome in alpha diversity, and undergone constantly contracting and shifting across the temporal scale. Nevertheless, 257 core virial clusters (VCs) persistently occurred with time whatever the significant impacts imposed by the varied physiochemical properties. Glycoside hydrolase (GH) and glycosyltransferase (GT) families genes displayed the higher abundances across all samples. Intriguingly, diversely structuring of toxin-antitoxin systems (TAs) and CRISPR-Cas arrays were frequently harbored by phage genomes. Their divergent organization and encoding attributes underlie the multiple biological roles in modulation of network and/or contest of phage community as well as bacterial host community. This phageome-wide mapping will fuel the current insights of phage community ecology in other traditional fermented ecosystems that are challenging to decipher.

Systematic genome-wide discovery of host factors governing bacteriophage infectivity

Bacterial host factors regulate the infection cycle of bacteriophages. Except for some well-studied host factors (e.g., receptors or restriction-modification systems), the contribution of the rest of the host genome on phage infection remains poorly understood. We developed PHAGEPACK, a pooled assay that systematically and comprehensively measures each host-gene impact on phage fitness. PHAGEPACK combines CRISPR interference with phage packaging to link host perturbation to phage fitness during active infection. Using PHAGEPACK, we constructed a genome-wide map of genes impacting T7 phage fitness in permissive E. coli, revealing pathways previously unknown to affect phage packaging. When applied to the non-permissive E. coli O121, PHAGEPACK identified pathways leading to host resistance; their removal increased phage susceptibility up to a billion-fold. Bioinformatic analysis indicates phage genomes carry homologs or truncations of key host factors, potentially for fitness advantage. In summary, PHAGEPACK offers valuable insights into phage-host interactions, phage evolution, and bacterial resistance.

Developmental Dynamics of the Gut Virome in Tibetan Pigs at High Altitude: A Metagenomic Perspective across Age Groups

Tibetan pig is a geographically isolated pig breed that inhabits high-altitude areas of the Qinghai-Tibetan plateau. At present, there is limited research on viral diseases in Tibetan pigs. This study provides a novel metagenomic exploration of the gut virome in Tibetan pigs (altitude ≈ 3000 m) across three critical developmental stages, including lactation, nursery, and fattening. The composition of viral communities in the Tibetan pig intestine, with a dominant presence of Microviridae phages observed across all stages of development, in combination with the previous literature, suggest that it may be associated with geographical locations with high altitude. Functional annotation of viral operational taxonomic units (vOTUs) highlights that, among the constantly increasing vOTUs groups, the adaptability of viruses to environmental stressors such as salt and heat indicates an evolutionary response to high-altitude conditions. It shows that the lactation group has more abundant viral auxiliary metabolic genes (vAMGs) than the nursery and fattening groups. During the nursery and fattening stages, this leaves only DNMT1 at a high level. which may be a contributing factor in promoting gut health. The study found that viruses preferentially adopt lytic lifestyles at all three developmental stages. These findings not only elucidate the dynamic interplay between the gut virome and host development, offering novel insights into the virome ecology of Tibetan pigs and their adaptation to high-altitude environments, but also provide a theoretical basis for further studies on pig production and epidemic prevention under extreme environmental conditions.

Diverse and abundant phages exploit conjugative plasmids

Phages exert profound evolutionary pressure on bacteria by interacting with receptors on the cell surface to initiate infection. While the majority of phages use chromosomally encoded cell surface structures as receptors, plasmid-dependent phages exploit plasmid-encoded conjugation proteins, making their host range dependent on horizontal transfer of the plasmid. Despite their unique biology and biotechnological significance, only a small number of plasmid-dependent phages have been characterized. Here we systematically search for new plasmid-dependent phages targeting IncP and IncF plasmids using a targeted discovery platform, and find that they are common and abundant in wastewater, and largely unexplored in terms of their genetic diversity. Plasmid-dependent phages are enriched in non-canonical types of phages, and all but one of the 65 phages we isolated were non-tailed, and members of the lipid-containing tectiviruses, ssDNA filamentous phages or ssRNA phages. We show that plasmid-dependent tectiviruses exhibit profound differences in their host range which is associated with variation in the phage holin protein. Despite their relatively high abundance in wastewater, plasmid-dependent tectiviruses are missed by metaviromic analyses, underscoring the continued importance of culture-based phage discovery. Finally, we identify a tailed phage dependent on the IncF plasmid, and find related structural genes in phages that use the orthogonal type 4 pilus as a receptor, highlighting the evolutionarily promiscuous use of these distinct contractile structures by multiple groups of phages. Taken together, these results indicate plasmid-dependent phages play an under-appreciated evolutionary role in constraining horizontal gene transfer via conjugative plasmids.

Hot springs viruses at Yellowstone National Park have ancient origins and are adapted to thermophilic hosts

Geothermal springs house unicellular red algae in the class Cyanidiophyceae that dominate the microbial biomass at these sites. Little is known about host-virus interactions in these environments. We analyzed the virus community associated with red algal mats in three neighboring habitats (creek, endolithic, soil) at Lemonade Creek, Yellowstone National Park (YNP), USA. We find that despite proximity, each habitat houses a unique collection of viruses, with the giant viruses, Megaviricetes, dominant in all three. The early branching phylogenetic position of genes encoded on metagenome assembled virus genomes (vMAGs) suggests that the YNP lineages are of ancient origin and not due to multiple invasions from mesophilic habitats. The existence of genomic footprints of adaptation to thermophily in the vMAGs is consistent with this idea. The Cyanidiophyceae at geothermal sites originated ca. 1.5 Bya and are therefore relevant to understanding biotic interactions on the early Earth.

Discovery and description of novel phage genomes from urban microbiomes sampled by the MetaSUB consortium

Bacteriophages are recognized as the most abundant members of microbiomes and have therefore a profound impact on microbial communities through the interactions with their bacterial hosts. The International Metagenomics and Metadesign of Subways and Urban Biomes Consortium (MetaSUB) has sampled mass-transit systems in 60 cities over 3 years using metagenomics, throwing light into these hitherto largely unexplored urban environments. MetaSUB focused primarily on the bacterial community. In this work, we explored MetaSUB metagenomic data in order to recover and analyze bacteriophage genomes. We recovered and analyzed 1714 phage genomes with size at least 40 kbp, from the class Caudoviricetes, the vast majority of which (80%) are novel. The recovered genomes were predicted to belong to temperate (69%) and lytic (31%) phages. Thirty-three of these genomes have more than 200 kbp, and one of them reaches 572 kbp, placing it among the largest phage genomes ever found. In general, the phages tended to be site-specific or nearly so, but 194 genomes could be identified in every city from which phage genomes were retrieved. We predicted hosts for 48% of the phages and observed general agreement between phage abundance and the respective bacterial host abundance, which include the most common nosocomial multidrug-resistant pathogens. A small fraction of the phage genomes are carriers of antibiotic resistance genes, and such genomes tended to be particularly abundant in the sites where they were found. We also detected CRISPR-Cas systems in five phage genomes. This study expands the previously reported MetaSUB results and is a contribution to the knowledge about phage diversity, global distribution, and phage genome content.

TemStaPro: protein thermostability prediction using sequence representations from protein language models

MOTIVATION

Reliable prediction of protein thermostability from its sequence is valuable for both academic and industrial research. This prediction problem can be tackled using machine learning and by taking advantage of the recent blossoming of deep learning methods for sequence analysis. These methods can facilitate training on more data and, possibly, enable the development of more versatile thermostability predictors for multiple ranges of temperatures.

RESULTS

We applied the principle of transfer learning to predict protein thermostability using embeddings generated by protein language models (pLMs) from an input protein sequence. We used large pLMs that were pre-trained on hundreds of millions of known sequences. The embeddings from such models allowed us to efficiently train and validate a high-performing prediction method using over one million sequences that we collected from organisms with annotated growth temperatures. Our method, TemStaPro (Temperatures of Stability for Proteins), was used to predict thermostability of CRISPR-Cas Class II effector proteins (C2EPs). Predictions indicated sharp differences among groups of C2EPs in terms of thermostability and were largely in tune with previously published and our newly obtained experimental data.

AVAILABILITY AND IMPLEMENTATION

TemStaPro software and the related data are freely available from https://github.com/ievapudz/TemStaPro and https://doi.org/10.5281/zenodo.7743637.

City-scale monitoring of antibiotic resistance genes by digital PCR and metagenomics

BACKGROUND

Anthropogenic activities significantly contribute to the dissemination of antibiotic resistance genes (ARGs), posing a substantial threat to humankind. The development of methods that allow robust ARG surveillance is a long-standing challenge. Here, we use city-scale monitoring of ARGs by using two of the most promising cutting-edge technologies, digital PCR (dPCR) and metagenomics.

METHODS

ARG hot-spots were sampled from the urban water and wastewater distribution systems. Metagenomics was used to provide a broad view of ARG relative abundance and richness in the prokaryotic and viral fractions. From the city-core ARGs in all samples, the worldwide dispersed sul2 and tetW conferring resistance to sulfonamide and tetracycline, respectively, were monitored by dPCR and metagenomics.

RESULTS

The largest relative overall ARG abundance and richness were detected in the hospital wastewater and the WWTP inlet (up to ≈6,000 ARGs/Gb metagenome) with a large fraction of unclassified resistant bacteria. The abundance of ARGs in DNA and RNA contigs classified as viruses was notably lower, demonstrating a reduction of up to three orders of magnitude compared to contigs associated to prokaryotes. By metagenomics and dPCR, a similar abundance tendency of sul2 and tetW was obtained, with higher abundances in hospital wastewater and WWTP input (≈125-225 ARGs/Gb metagenome). dPCR absolute abundances were between 6,000 and 18,600 copies per ng of sewage DNA (≈105-7 copies/mL) and 6.8 copies/mL in seawater near the WWTP discharging point.

CONCLUSIONS

dPCR was more sensitive and accurate, while metagenomics provided broader coverage of ARG detection. While desirable, a reliable correlation of dPCR absolute abundance units into metagenomic relative abundance units was not obtained here (r2 < 0.4) suggesting methodological factors that introduce variability. Evolutionary pressure does not significantly select the targeted ARGs in natural aquatic environments.

Draft genome sequences of Arthrobacter sp. AZCC_0090 and Mycobacterium sp. AZCC_0083 isolated from oligotrophic subsurface forest soil in the Santa Catalina mountains of Southern Arizona

Here, we present the genomes of two soil actinobacteria: Arthrobacter sp. strain AZCC_0090 and Mycobacterium sp. strain AZCC_0083, isolated from oligotrophic subsurface soils in Southern Arizona, USA.

MetaCerberus: distributed highly parallelized HMM-based processing for robust functional annotation across the tree of life

MOTIVATION

MetaCerberus is a massively parallel, fast, low memory, scalable annotation tool for inference gene function across genomes to metacommunities. MetaCerberus provides an elusive HMM/HMMER-based tool at a rapid scale with low memory. It offers scalable gene elucidation to major public databases, including KEGG (KO), COGs, CAZy, FOAM, and specific databases for viruses, including VOGs and PHROGs, from single genomes to metacommunities.

RESULTS

MetaCerberus is 1.3× as fast on a single node than eggNOG-mapper v2 on 5× less memory using an exclusively HMM/HMMER mode. In a direct comparison, MetaCerberus provides better annotation of viruses, phages, and archaeal viruses than DRAM, Prokka, or InterProScan. MetaCerberus annotates more KOs across domains when compared to DRAM, with a 186× smaller database, and with 63× less memory. MetaCerberus is fully integrated for automatic analysis of statistics and pathways using differential statistic tools (i.e. DESeq2 and edgeR), pathway enrichment (GAGE R), and pathview R. MetaCerberus provides a novel tool for unlocking the biosphere across the tree of life at scale.

AVAILABILITY AND IMPLEMENTATION

MetaCerberus is written in Python and distributed under a BSD-3 license. The source code of MetaCerberus is freely available at https://github.com/raw-lab/metacerberus compatible with Python 3 and works on both Mac OS X and Linux. MetaCerberus can also be easily installed using bioconda: mamba create -n metacerberus -c bioconda -c conda-forge metacerberus.

Forty-nine metagenomic-assembled genomes from an aquatic virome expand Caudoviricetes by 45 potential new families and the newly uncovered Gossevirus of Bamfordvirae

Twenty complete genomes (29-63 kb) and 29 genomes with an estimated completeness of over 90 % (30-90 kb) were identified for novel dsDNA viruses in the Yangshan Harbor metavirome. These newly discovered viruses contribute to the expansion of viral taxonomy by introducing 46 potential new families. Except for one virus, all others belong to the class Caudoviricetes. The exception is a novel member of the recently characterized viral group known as Gossevirus. Fifteen viruses were predicted to be temperate. The predicted hosts for the viruses appear to be involved in various aspects of the nitrogen cycle, including nitrogen fixation, oxidation and denitrification. Two viruses were identified to have a host of Flavobacterium and Tepidimonas fonticaldi, respectively, by matching CRISPR spacers with viral protospacers. Our findings provide an overview for characterizing and identifying specific viruses from Yangshan Harbor. The Gossevirus-like virus uncovered emphasizes the need for further comprehensive isolation and investigation of polinton-like viruses.

Viral potential to modulate microbial methane metabolism varies by habitat

Methane is a potent greenhouse gas contributing to global warming. Microorganisms largely drive the biogeochemical cycling of methane, yet little is known about viral contributions to methane metabolism (MM). We analyzed 982 publicly available metagenomes from host-associated and environmental habitats containing microbial MM genes, expanding the known MM auxiliary metabolic genes (AMGs) from three to 24, including seven genes exclusive to MM pathways. These AMGs are recovered on 911 viral contigs predicted to infect 14 prokaryotic phyla including Halobacteriota, Methanobacteriota, and Thermoproteota. Of those 24, most were encoded by viruses from rumen (16/24), with substantially fewer by viruses from environmental habitats (0-7/24). To search for additional MM AMGs from an environmental habitat, we generate metagenomes from methane-rich sediments in Vrana Lake, Croatia. Therein, we find diverse viral communities, with most viruses predicted to infect methanogens and methanotrophs and some encoding 13 AMGs that can modulate host metabolisms. However, none of these AMGs directly participate in MM pathways. Together these findings suggest that the extent to which viruses use AMGs to modulate host metabolic processes (e.g., MM) varies depending on the ecological properties of the habitat in which they dwell and is not always predictable by habitat biogeochemical properties.

Role of vertical and horizontal microbial transmission of antimicrobial resistance genes in early life: insights from maternal-infant dyads

Early life represents a critical window for metabolic, cognitive and immune system development, which is influenced by the maternal microbiome as well as the infant gut microbiome. Antibiotic exposure, mode of delivery and breastfeeding practices modulate the gut microbiome and the reservoir of antibiotic resistance genes (ARGs). Vertical and horizontal microbial gene transfer during early life and the mechanisms behind these transfers are being uncovered. In this review, we aim to provide an overview of the current knowledge on the transfer of antibiotic resistance in the mother-infant dyad through vertical and horizontal transmission and to highlight the main gaps and challenges in this area.

Double-stranded RNA sequencing reveals distinct riboviruses associated with thermoacidophilic bacteria from hot springs in Japan

Metatranscriptome sequencing expanded the known diversity of the bacterial RNA virome, suggesting that additional riboviruses infecting bacterial hosts remain to be discovered. Here we employed double-stranded RNA sequencing to recover complete genome sequences of two ribovirus groups from acidic hot springs in Japan. One group, denoted hot spring riboviruses (HsRV), consists of viruses with distinct RNA-directed RNA polymerases (RdRPs) that seem to be intermediates between typical ribovirus RdRPs and viral reverse transcriptases. This group forms a distinct phylum, Artimaviricota, or even kingdom within the realm Riboviria. We identified viruses encoding HsRV-like RdRPs in marine water, river sediments and salt marshes, indicating that this group is widespread beyond extreme ecosystems. The second group, denoted hot spring partiti-like viruses (HsPV), forms a distinct branch within the family Partitiviridae. The genome architectures of HsRV and HsPV and their identification in bacteria-dominated habitats suggest that these viruses infect thermoacidophilic bacteria.

The oral microbiome: diversity, biogeography and human health

The human oral microbiota is highly diverse and has a complex ecology, comprising bacteria, microeukaryotes, archaea and viruses. These communities have elaborate and highly structured biogeography that shapes metabolic exchange on a local scale and results from the diverse microenvironments present in the oral cavity. The oral microbiota also interfaces with the immune system of the human host and has an important role in not only the health of the oral cavity but also systemic health. In this Review, we highlight recent advances including novel insights into the biogeography of several oral niches at the species level, as well as the ecological role of candidate phyla radiation bacteria and non-bacterial members of the oral microbiome. In addition, we summarize the relationship between the oral microbiota and the pathology of oral diseases and systemic diseases. Together, these advances move the field towards a more holistic understanding of the oral microbiota and its role in health, which in turn opens the door to the study of novel preventive and therapeutic strategies.

Hiding in plain sight: The discovery of complete genomes of 11 hypothetical spindle-shaped viruses that putatively infect mesophilic ammonia-oxidizing archaea

The genome of a putative Nitrosopumilaceae virus with a hypothetical spindle-shaped particle morphology was identified in the Yangshan Harbour metavirome from the East China Sea through protein similarity comparison and structure analysis. This discovery was accompanied by a set of 10 geographically dispersed close relatives found in the environmental virus datasets from typical locations of ammonia-oxidizing archaeon distribution. Its host prediction was supported by iPHoP prediction and protein sequence similarity. The structure of the predicted major capsid protein, together with the overall N-glycosylation site, the transmembrane helices prediction, the hydrophilicity profile, and the docking simulation of the major capsid proteins, indicate that these viruses resemble spindle-shaped viruses. It suggests a similarly assembled structure and, consequently, a possibly spindle-shaped morphology of these newly discovered archaeal viruses.

Ultraconserved bacteriophage genome sequence identified in 1300-year-old human palaeofaeces

Bacteriophages are widely recognised as rapidly evolving biological entities. However, knowledge about ancient bacteriophages is limited. Here, we analyse DNA sequence datasets previously generated from ancient palaeofaeces and human gut-content samples, and identify an ancient phage genome nearly identical to present-day Mushuvirus mushu, a virus that infects gut commensal bacteria. The DNA damage patterns of the genome are consistent with its ancient origin and, despite 1300 years of evolution, the ancient Mushuvirus genome shares 97.7% nucleotide identity with its modern counterpart, indicating a long-term relationship between the prophage and its host. In addition, we reconstruct and authenticate 297 other phage genomes from the last 5300 years, including those belonging to unknown families. Our findings demonstrate the feasibility of reconstructing ancient phage genome sequences, thus expanding the known virosphere and offering insights into phage-bacteria interactions spanning several millennia.

ProkBERT family: genomic language models for microbiome applications

Background

In the evolving landscape of microbiology and microbiome analysis, the integration of machine learning is crucial for understanding complex microbial interactions, and predicting and recognizing novel functionalities within extensive datasets. However, the effectiveness of these methods in microbiology faces challenges due to the complex and heterogeneous nature of microbial data, further complicated by low signal-to-noise ratios, context-dependency, and a significant shortage of appropriately labeled datasets. This study introduces the ProkBERT model family, a collection of large language models, designed for genomic tasks. It provides a generalizable sequence representation for nucleotide sequences, learned from unlabeled genome data. This approach helps overcome the above-mentioned limitations in the field, thereby improving our understanding of microbial ecosystems and their impact on health and disease.

Methods

ProkBERT models are based on transfer learning and self-supervised methodologies, enabling them to use the abundant yet complex microbial data effectively. The introduction of the novel Local Context-Aware (LCA) tokenization technique marks a significant advancement, allowing ProkBERT to overcome the contextual limitations of traditional transformer models. This methodology not only retains rich local context but also demonstrates remarkable adaptability across various bioinformatics tasks.

Results

In practical applications such as promoter prediction and phage identification, the ProkBERT models show superior performance. For promoter prediction tasks, the top-performing model achieved a Matthews Correlation Coefficient (MCC) of 0.74 for E. coli and 0.62 in mixed-species contexts. In phage identification, ProkBERT models consistently outperformed established tools like VirSorter2 and DeepVirFinder, achieving an MCC of 0.85. These results underscore the models' exceptional accuracy and generalizability in both supervised and unsupervised tasks.

Conclusions

The ProkBERT model family is a compact yet powerful tool in the field of microbiology and bioinformatics. Its capacity for rapid, accurate analyses and its adaptability across a spectrum of tasks marks a significant advancement in machine learning applications in microbiology. The models are available on GitHub (https://github.com/nbrg-ppcu/prokbert) and HuggingFace (https://huggingface.co/nerualbioinfo) providing an accessible tool for the community.

Macroalgal virosphere assists with host-microbiome equilibrium regulation and affects prokaryotes in surrounding marine environments

The microbiomes in macroalgal holobionts play vital roles in regulating macroalgal growth and ocean carbon cycling. However, the virospheres in macroalgal holobionts remain largely underexplored, representing a critical knowledge gap. Here we unveil that the holobiont of kelp (Saccharina japonica) harbors highly specific and unique epiphytic/endophytic viral species, with novelty (99.7% unknown) surpassing even extreme marine habitats (e.g. deep-sea and hadal zones), indicating that macroalgal virospheres, despite being closest to us, are among the least understood. These viruses potentially maintain microbiome equilibrium critical for kelp health via lytic-lysogenic infections and the expression of folate biosynthesis genes. In-situ kelp mesocosm cultivation and metagenomic mining revealed that kelp holobiont profoundly reshaped surrounding seawater and sediment virus-prokaryote pairings through changing surrounding environmental conditions and virus-host migrations. Some kelp epiphytic viruses could even infect sediment autochthonous bacteria after deposition. Moreover, the presence of ample viral auxiliary metabolic genes for kelp polysaccharide (e.g. laminarin) degradation underscores the underappreciated viral metabolic influence on macroalgal carbon cycling. This study provides key insights into understanding the previously overlooked ecological significance of viruses within macroalgal holobionts and the macroalgae-prokaryotes-virus tripartite relationship.

IMG/PR: a database of plasmids from genomes and metagenomes with rich annotations and metadata

Plasmids are mobile genetic elements found in many clades of Archaea and Bacteria. They drive horizontal gene transfer, impacting ecological and evolutionary processes within microbial communities, and hold substantial importance in human health and biotechnology. To support plasmid research and provide scientists with data of an unprecedented diversity of plasmid sequences, we introduce the IMG/PR database, a new resource encompassing 699 973 plasmid sequences derived from genomes, metagenomes and metatranscriptomes. IMG/PR is the first database to provide data of plasmid that were systematically identified from diverse microbiome samples. IMG/PR plasmids are associated with rich metadata that includes geographical and ecosystem information, host taxonomy, similarity to other plasmids, functional annotation, presence of genes involved in conjugation and antibiotic resistance. The database offers diverse methods for exploring its extensive plasmid collection, enabling users to navigate plasmids through metadata-centric queries, plasmid comparisons and BLAST searches. The web interface for IMG/PR is accessible at https://img.jgi.doe.gov/pr. Plasmid metadata and sequences can be downloaded from https://genome.jgi.doe.gov/portal/IMG_PR.

Gut virome in inflammatory bowel disease and beyond

OBJECTIVE

The gut virome is a dense community of viruses inhabiting the gastrointestinal tract and an integral part of the microbiota. The virome coexists with the other components of the microbiota and with the host in a dynamic equilibrium, serving as a key contributor to the maintenance of intestinal homeostasis and functions. However, this equilibrium can be interrupted in certain pathological states, including inflammatory bowel disease, causing dysbiosis that may participate in disease pathogenesis. Nevertheless, whether virome dysbiosis is a causal or bystander event requires further clarification.

DESIGN

This review seeks to summarise the latest advancements in the study of the gut virome, highlighting its cross-talk with the mucosal microenvironment. It explores how cutting-edge technologies may build upon current knowledge to advance research in this field. An overview of virome transplantation in diseased gastrointestinal tracts is provided along with insights into the development of innovative virome-based therapeutics to improve clinical management.

RESULTS

Gut virome dysbiosis, primarily driven by the expansion of Caudovirales, has been shown to impact intestinal immunity and barrier functions, influencing overall intestinal homeostasis. Although emerging innovative technologies still need further implementation, they display the unprecedented potential to better characterise virome composition and delineate its role in intestinal diseases.

CONCLUSIONS

The field of gut virome is progressively expanding, thanks to the advancements of sequencing technologies and bioinformatic pipelines. These have contributed to a better understanding of how virome dysbiosis is linked to intestinal disease pathogenesis and how the modulation of virome composition may help the clinical intervention to ameliorate gut disease management.

PhageScope: a well-annotated bacteriophage database with automatic analyses and visualizations

Bacteriophages are viruses that infect bacteria or archaea. Understanding the diverse and intricate genomic architectures of phages is essential to study microbial ecosystems and develop phage therapy strategies. However, the existing phage databases are short of meticulous annotations. To this end, we propose PhageScope (https://phagescope.deepomics.org), an online phage database with comprehensive annotations. PhageScope harbors a collection of 873 718 phage sequences from various sources. Applying fifteen state-of-the-art tools to perform systematic annotations and analyses, PhageScope provides annotations on genome completeness, host range, lifestyle information, taxonomy classification, nine types of structural and functional genetic elements, and three types of comparative genomic studies for curated phages. Additionally, PhageScope incorporates automatic analyses and visualizations for curated and customized phages, serving as an efficient platform for phage study.

dbAPIS: a database of anti-prokaryotic immune system genes

Anti-prokaryotic immune system (APIS) proteins, typically encoded by phages, prophages, and plasmids, inhibit prokaryotic immune systems (e.g. restriction modification, toxin-antitoxin, CRISPR-Cas). A growing number of APIS genes have been characterized and dispersed in the literature. Here we developed dbAPIS (https://bcb.unl.edu/dbAPIS), as the first literature curated data repository for experimentally verified APIS genes and their associated protein families. The key features of dbAPIS include: (i) experimentally verified APIS genes with their protein sequences, functional annotation, PDB or AlphaFold predicted structures, genomic context, sequence and structural homologs from different microbiome/virome databases; (ii) classification of APIS proteins into sequence-based families and construction of hidden Markov models (HMMs); (iii) user-friendly web interface for data browsing by the inhibited immune system types or by the hosts, and functions for searching and batch downloading of pre-computed data; (iv) Inclusion of all types of APIS proteins (except for anti-CRISPRs) that inhibit a variety of prokaryotic defense systems (e.g. RM, TA, CBASS, Thoeris, Gabija). The current release of dbAPIS contains 41 verified APIS proteins and ∼4400 sequence homologs of 92 families and 38 clans. dbAPIS will facilitate the discovery of novel anti-defense genes and genomic islands in phages, by providing a user-friendly data repository and a web resource for an easy homology search against known APIS proteins.

Metatranscriptomic data mining together with microfluidic card uncovered the potential pathogens and seasonal RNA viral ecology in a drinking water source

AIMS

Despite metatranscriptomics becoming an emerging tool for pathogen surveillance, very little is known about the feasibility of this approach for understanding the fate of human-derived pathogens in drinking water sources.

METHODS AND RESULTS

We conducted multiplexed microfluidic cards and metatranscriptomic sequencing of the drinking water source in a border city of North Korea in four seasons. Microfluidic card detected norovirus, hepatitis B virus (HBV), enterovirus, and Vibrio cholerae in the water. Phylogenetic analyses showed that environmental-derived sequences from norovirus GII.17, genotype C of HBV, and coxsackievirus A6 (CA6) were genetically related to the local clinical isolates. Meanwhile, metatranscriptomic assembly suggested that several bacterial pathogens, including Acinetobacter johnsonii and V. cholerae might be prevalent in the studied region. Metatranscriptomic analysis recovered 349 species-level groups with substantial viral diversity without detection of norovirus, HBV, and CA6. Seasonally distinct virus communities were also found. Specifically, 126, 73, 126, and 457 types of viruses were identified in spring, summer, autumn, and winter, respectively. The viromes were dominated by the Pisuviricota phylum, including members from Marnaviridae, Dicistroviridae, Luteoviridae, Potyviridae, Picornaviridae, Astroviridae, and Picobirnaviridae families. Further phylogenetic analyses of RNA (Ribonucleic Acid)-dependent RNA polymerase (RdRp) sequences showed a diverse set of picorna-like viruses associated with shellfish, of which several novel picorna-like viruses were also identified. Additionally, potential animal pathogens, including infectious bronchitis virus, Bat dicibavirus, Bat nodavirus, Bat picornavirus 2, infectious bursal disease virus, and Macrobrachium rosenbergii nodavirus were also identified.

CONCLUSIONS

Our data illustrate the divergence between microfluidic cards and metatranscriptomics, highlighting that the combination of both methods facilitates the source tracking of human viruses in challenging settings without sufficient clinical surveillance.

IPEV: identification of prokaryotic and eukaryotic virus-derived sequences in virome using deep learning

BACKGROUND

The virome obtained through virus-like particle enrichment contains a mixture of prokaryotic and eukaryotic virus-derived fragments. Accurate identification and classification of these elements are crucial to understanding their roles and functions in microbial communities. However, the rapid mutation rates of viral genomes pose challenges in developing high-performance tools for classification, potentially limiting downstream analyses.

FINDINGS

We present IPEV, a novel method to distinguish prokaryotic and eukaryotic viruses in viromes, with a 2-dimensional convolutional neural network combining trinucleotide pair relative distance and frequency. Cross-validation assessments of IPEV demonstrate its state-of-the-art precision, significantly improving the F1-score by approximately 22% on an independent test set compared to existing methods when query viruses share less than 30% sequence similarity with known viruses. Furthermore, IPEV outperforms other methods in accuracy on marine and gut virome samples based on annotations by sequence alignments. IPEV reduces runtime by at most 1,225 times compared to existing methods under the same computing configuration. We also utilized IPEV to analyze longitudinal samples and found that the gut virome exhibits a higher degree of temporal stability than previously observed in persistent personal viromes, providing novel insights into the resilience of the gut virome in individuals.

CONCLUSIONS

IPEV is a high-performance, user-friendly tool that assists biologists in identifying and classifying prokaryotic and eukaryotic viruses within viromes. The tool is available at https://github.com/basehc/IPEV.

Mixed waste contamination selects for a mobile genetic element population enriched in multiple heavy metal resistance genes

Mobile genetic elements (MGEs) like plasmids, viruses, and transposable elements can provide fitness benefits to their hosts for survival in the presence of environmental stressors. Heavy metal resistance genes (HMRGs) are frequently observed on MGEs, suggesting that MGEs may be an important driver of adaptive evolution in environments contaminated with heavy metals. Here, we report the meta-mobilome of the heavy metal-contaminated regions of the Oak Ridge Reservation subsurface. This meta-mobilome was compared with one derived from samples collected from unimpacted regions of the Oak Ridge Reservation subsurface. We assembled 1615 unique circularized DNA elements that we propose to be MGEs. The circular elements from the highly contaminated subsurface were enriched in HMRG clusters relative to those from the nearby unimpacted regions. Additionally, we found that these HMRGs were associated with Gamma and Betaproteobacteria hosts in the contaminated subsurface and potentially facilitate the persistence and dominance of these taxa in this region. Finally, the HMRGs were associated with conjugative elements, suggesting their potential for future lateral transfer. We demonstrate how our understanding of MGE ecology, evolution, and function can be enhanced through the genomic context provided by completed MGE assemblies.

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.

Comparative Analyses of Bacteriophage Genomes

Bacterial viruses (bacteriophages or phages) are the most abundant and diverse biological entities on Earth. There is a renewed worldwide interest in phage-centered research motivated by their enormous potential as antimicrobials to cope with multidrug-resistant pathogens. An ever-growing number of complete phage genomes are becoming available, derived either from newly isolated phages (cultivated phages) or recovered from metagenomic sequencing data (uncultivated phages). Robust comparative analysis is crucial for a comprehensive understanding of genotypic variations of phages and their related evolutionary processes, and to investigate the interaction mechanisms between phages and their hosts. In this chapter, we present a protocol for phage comparative genomics employing tools selected out of the many currently available, focusing on complete genomes of phages classified in the class Caudoviricetes. This protocol provides accurate identification of similarities, differences, and patterns among new and previously known complete phage genomes as well as phage clustering and taxonomic classification.

Diverse and abundant phages exploit conjugative plasmids

Phages exert profound evolutionary pressure on bacteria by interacting with receptors on the cell surface to initiate infection. While the majority of phages use chromosomally-encoded cell surface structures as receptors, plasmid-dependent phages exploit plasmid-encoded conjugation proteins, making their host range dependent on horizontal transfer of the plasmid. Despite their unique biology and biotechnological significance, only a small number of plasmid-dependent phages have been characterized. Here we systematically search for new plasmid-dependent phages targeting IncP and IncF plasmids using a targeted discovery platform, and find that they are common and abundant in wastewater, and largely unexplored in terms of their genetic diversity. Plasmid-dependent phages are enriched in non-canonical types of phages, and all but one of the 64 phages we isolated were non-tailed, and members of the lipid-containing tectiviruses, ssDNA filamentous phages or ssRNA phages. We show that plasmid-dependent tectiviruses exhibit profound differences in their host range which is associated with variation in the phage holin protein. Despite their relatively high abundance in wastewater, plasmid-dependent tectiviruses are missed by metaviromic analyses, underscoring the continued importance of culture-based phage discovery. Finally, we identify a tailed phage dependent on the IncF plasmid, and find related structural genes in phages that use the orthogonal type 4 pilus as a receptor, highlighting the promiscuous use of these distinct contractile structures by multiple groups of phages. Taken together, these results indicate plasmid-dependent phages play an under-appreciated evolutionary role in constraining horizontal gene transfer via conjugative plasmids.

Distribution and Implications of Haloarchaeal Plasmids Disseminated in Self-Encoded Plasmid Vesicles

Even though viruses and plasmids are both drivers of horizontal gene transfer, they differ fundamentally in their mode of transfer. Virus genomes are enclosed in virus capsids and are not dependent on cell-to-cell contacts for their dissemination. In contrast, the transfer of plasmids most often requires physical contact between cells. However, plasmid pR1SE of Halorubrum lacusprofundi is disseminated between cells, independent of cell-cell contacts, in specialized membrane vesicles that contain plasmid proteins. In this study, we searched for pR1SE-like elements in public databases and a metagenomics dataset from Australian salt lakes and identified 40 additional pR1SE-like elements in hypersaline environments worldwide. Herein, these elements are named apHPVs (archaeal plasmids of haloarchaea potentially transferred in plasmid vesicles). They share two sets of closely related proteins with conserved synteny, strongly indicating an organization into different functional clusters. We find that apHPVs, besides transferring themselves, have the potential to transfer large fragments of DNA between host cells, including virus defense systems. Most interestingly, apHPVs likely play an important role in the evolution of viruses and plasmids in haloarchaea, as they appear to recombine with both of them. This further supports the idea that plasmids and viruses are not distinct but closely related mobile genetic elements.

Viromes vs. mixed community metagenomes: choice of method dictates interpretation of viral community ecology

Background

Viruses, the majority of which are uncultivated, are among the most abundant biological entities on Earth. From altering microbial physiology to driving community dynamics, viruses are fundamental members of microbiomes. While the number of studies leveraging viral metagenomics (viromics) for studying uncultivated viruses is growing, standards for viromics research are lacking. Viromics can utilize computational discovery of viruses from total metagenomes of all community members (hereafter metagenomes) or use physical separation of virus-specific fractions (hereafter viromes). However, differences in the recovery and interpretation of viruses from metagenomes and viromes obtained from the same samples remain understudied.

Results

Here, we compare viral communities from paired viromes and metagenomes obtained from 60 diverse samples across human gut, soil, freshwater, and marine ecosystems. Overall, viral communities obtained from viromes were more abundant and species rich than those obtained from metagenomes, although there were some exceptions. Despite this, metagenomes still contained many viral genomes not detected in viromes. We also found notable differences in the predicted lytic state of viruses detected in viromes vs metagenomes at the time of sequencing. Other forms of variation observed include genome presence/absence, genome quality, and encoded protein content between viromes and metagenomes, but the magnitude of these differences varied by environment.

Conclusions

Overall, our results show that the choice of method can lead to differing interpretations of viral community ecology. We suggest that the choice of whether to target a metagenome or virome to study viral communities should be dependent on the environmental context and ecological questions being asked. However, our overall recommendation to researchers investigating viral ecology and evolution is to pair both approaches to maximize their respective benefits.

Virophages Found in Viromes from Lake Baikal

In this study, a previously little-studied group of viruses-virophages-was searched for and identified in the viromes of the ancient oligotrophic Lake Baikal. Virophages are small dsDNA viruses that parasitize giant viruses (e.g., Mimiviridae), which in turn affect unicellular eukaryotes. We analyzed eight viromes obtained from the deep-water areas of three basins of Lake Baikal and the shallow-water strait Maloye More in different seasons. The sequences of virophages were revealed in all viromes and were dominant after bacteriophages and algal viruses. Sixteen putative complete genomes of virophages were assembled, all of which contained four conserved genes encoding major capsid protein (MCP), minor capsid protein (mCP), maturation cysteine protease (PRO), and FtsK-HerA family DNA-packaging ATPase (ATPase). The MCP-based cluster analysis showed a sequence separation according to seasons, and a dependence on the geographical localization was not detected.

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.

Hi-C metagenome sequencing reveals soil phage-host interactions

Bacteriophages are abundant in soils. However, the majority are uncharacterized, and their hosts are unknown. Here, we apply high-throughput chromosome conformation capture (Hi-C) to directly capture phage-host relationships. Some hosts have high centralities in bacterial community co-occurrence networks, suggesting phage infections have an important impact on the soil bacterial community interactions. We observe increased average viral copies per host (VPH) and decreased viral transcriptional activity following a two-week soil-drying incubation, indicating an increase in lysogenic infections. Soil drying also alters the observed phage host range. A significant negative correlation between VPH and host abundance prior to drying indicates more lytic infections result in more host death and inversely influence host abundance. This study provides empirical evidence of phage-mediated bacterial population dynamics in soil by directly capturing specific phage-host interactions.

Enhanced Bacterium-Phage Symbiosis in Attached Microbial Aggregates on a Membrane Surface Facing Elevated Hydraulic Stress

Phages are increasingly recognized for their importance in microbial aggregates, including their influence on microbial ecosystem services and biotechnology applications. However, the adaptive strategies and ecological functions of phages in different aggregates remain largely unexplored. Herein, we used membrane bioreactors to investigate bacterium-phage interactions and related microbial functions within suspended and attached microbial aggregates (SMA vs AMA). SMA and AMA represent distinct microbial habitats where bacterial communities display distinct patterns in terms of dominant species, keystone species, and bacterial networks. However, bacteria and phages in both aggregates exhibited high lysogenicity, with 60% lysogenic phages in the virome and 70% lysogenic metagenome-assembled genomes of bacteria. Moreover, substantial phages exhibited broad host ranges (34% in SMA and 42% in AMA) and closely interacted with habitat generalist species (43% in SMA and 49% in AMA) as adaptive strategies in stressful operation environments. Following a mutualistic pattern, phage-carried auxiliary metabolic genes (pAMGs; 238 types in total) presumably contributed to the bacterial survival and aggregate stability. The SMA-pAMGs were mainly associated with energy metabolism, while the AMA-pAMGs were mainly associated with antioxidant biosynthesis and the synthesis of extracellular polymeric substances, representing habitat-dependent patterns. Overall, this study advanced our understanding of phage adaptive strategies in microbial aggregate habitats and emphasized the importance of bacterium-phage symbiosis in the stability of microbial aggregates.

Characterization and diversity of CRISPR/Cas systems in Klebsiella oxytoca

CRISPR (clustered regularly interspaced short palindromic repeats)/Cas (CRISPR-associated protein) system is a crucial adaptive immune system for bacteria to resist foreign DNA infection. In this study, we investigated the prevalence and diversity of CRISPR/Cas systems in 175 Klebsiella oxytoca (K. oxytoca) strains. Specifically, 58.86% (103/175) of these strains possessed at least one confirmed CRISPR locus. Two CRISPR/Cas system types, I-F and IV-A3, were identified in 69 strains. Type I-F system was the most prevalent in this species, which correlated well with MLST. Differently, type IV-A3 system was randomly distributed. Moreover, the type IV-A3 system was separated into two subgroups, with subgroup-specific cas genes and repeat sequences. In addition, spacer origin analysis revealed that approximately one-fifth of type I-F spacers and one-third of type IV-A3 spacers had a significant match to MGEs. The phage tail tape measure protein and conjunctive transfer system protein were important targets of type I-F and IV-A3 systems in K. oxytoca, respectively. PAM sequences were inferred to be 5'-NCC-3' for type I-F, 5'-AAG-3' for subgroup IV-A3-a, and 5'-AAN-3' for subgroup IV-A3-b. Collectively, our findings will shed light on the prevalence, diversity, and functional effects of the CRISPR/Cas system in K. oxytoca.

Virus diversity and activity is driven by snowmelt and host dynamics in a high-altitude watershed soil ecosystem

BACKGROUND

Viruses impact nearly all organisms on Earth, including microbial communities and their associated biogeochemical processes. In soils, highly diverse viral communities have been identified, with a global distribution seemingly driven by multiple biotic and abiotic factors, especially soil temperature and moisture. However, our current understanding of the stability of soil viral communities across time and their response to strong seasonal changes in environmental parameters remains limited. Here, we investigated the diversity and activity of environmental soil DNA and RNA viruses, focusing especially on bacteriophages, across dynamics' seasonal changes in a snow-dominated mountainous watershed by examining paired metagenomes and metatranscriptomes.

RESULTS

We identified a large number of DNA and RNA viruses taxonomically divergent from existing environmental viruses, including a significant proportion of fungal RNA viruses, and a large and unsuspected diversity of positive single-stranded RNA phages (Leviviricetes), highlighting the under-characterization of the global soil virosphere. Among these, we were able to distinguish subsets of active DNA and RNA phages that changed across seasons, consistent with a "seed-bank" viral community structure in which new phage activity, for example, replication and host lysis, is sequentially triggered by changes in environmental conditions. At the population level, we further identified virus-host dynamics matching two existing ecological models: "Kill-The-Winner" which proposes that lytic phages are actively infecting abundant bacteria, and "Piggyback-The-Persistent" which argues that when the host is growing slowly, it is more beneficial to remain in a dormant state. The former was associated with summer months of high and rapid microbial activity, and the latter with winter months of limited and slow host growth.

CONCLUSION

Taken together, these results suggest that the high diversity of viruses in soils is likely associated with a broad range of host interaction types each adapted to specific host ecological strategies and environmental conditions. As our understanding of how environmental and host factors drive viral activity in soil ecosystems progresses, integrating these viral impacts in complex natural microbiome models will be key to accurately predict ecosystem biogeochemistry. Video Abstract.

Metavirome-assembled genome sequence of a new aquatic RNA virus expands the genus Locarnavirus

RNA viruses in marine environments have long been recognized as major players in the virosphere. In this study, the complete genome sequence of an RNA virus from Yangshan Harbor, named marine RNA virus Yangshan-LWW (YS-LWW), was obtained based on metavirome assembly. The genome of YS-LWW is 8839 nt in length and contains two open reading frames (ORFs). Both RdRP- and whole-genome-based phylogenetic analysis showed that YS-LWW, together with 45 viral isolates with sequences in public datasets, represents a new species in the genus Locarnavirus of the family Marnaviridae. PCR and public dataset mining indicate that YS-LWW and YS-LWW-like viruses have been widely detected in coastal and freshwater environments, suggesting that they might play a significant role in aquatic ecosystems.

Characterization and genomic analysis of Stutzerimonas stutzeri phage vB_PstS_ZQG1, representing a novel viral genus

Stutzerimonas stutzeri is an opportunistic pathogenic bacterium belonging to the Gammaproteobacteria, exhibiting wide distribution in the environment and playing significant ecological roles such as nitrogen fixation or pollutant degradation. Despite its ecological importance, only two S. stutzeri phages have been isolated to date. Here, a novel S. stutzeri phage, vB_PstS_ZQG1, was isolated from the surface seawater of Qingdao, China. Transmission electron microscopy analysis indicates that vB_PstS_ZQG1 has a morphology characterized by a long non-contractile tail. The genomic sequence of vB_PstS_ZQG1 contains a linear, double-strand 61,790-bp with the G+C content of 53.24% and encodes 90 putative open reading frames. Two auxiliary metabolic genes encoding TolA protein and nucleotide pyrophosphohydrolase were identified, which are likely involved in host adaptation and phage reproduction. Phylogenetic and comparative genomic analyses demonstrated that vB_PstS_ZQG1 exhibits low similarity with previously isolated phages or uncultured viruses (average nucleotide identity values range from 21.7 to 29.4), suggesting that it represents a novel viral genus by itself, here named as Fuevirus. Biogeographic analysis showed that vB_PstS_ZQG1 was only detected in epipelagic and mesopelagic zone with low abundance. In summary, our findings of the phage vB_PstS_ZQG1 will provide helpful insights for further research on the interactions between S. stutzeri phages and their hosts, and contribute to discovering unknown viral sequences in the metagenomic database.

Viruses in deep-sea cold seep sediments harbor diverse survival mechanisms and remain genetically conserved within species

Deep sea cold seep sediments have been discovered to harbor novel, abundant, and diverse bacterial and archaeal viruses. However, little is known about viral genetic features and evolutionary patterns in these environments. Here, we examined the evolutionary ecology of viruses across active and extinct seep stages in the area of Haima cold seeps in the South China Sea. A total of 338 viral operational taxonomic units are identified and linked to 36 bacterial and archaeal phyla. The dynamics of host-virus interactions are informed by diverse antiviral defense systems across 43 families found in 487 microbial genomes. Cold seep viruses are predicted to harbor diverse adaptive strategies to persist in this environment, including counter-defense systems, auxiliary metabolic genes, reverse transcriptases, and alternative genetic code assignments. Extremely low nucleotide diversity is observed in cold seep viral populations, being influenced by factors including microbial host, sediment depth, and cold seep stage. Most cold seep viral genes are under strong purifying selection with trajectories that differ depending on whether cold seeps are active or extinct. This work sheds light on the understanding of environmental adaptation mechanisms and evolutionary patterns of viruses in the sub-seafloor biosphere.

Genomic Analysis and Taxonomic Characterization of Seven Bacteriophage Genomes Metagenomic-Assembled from the Dishui Lake

Viruses in aquatic ecosystems exhibit remarkable abundance and diversity. However, scattered studies have been conducted to mine uncultured viruses and identify them taxonomically in lake water. Here, whole genomes (29-173 kbp) of seven uncultured dsDNA bacteriophages were discovered in Dishui Lake, the largest artificial lake in Shanghai. We analyzed their genomic signatures and found a series of viral auxiliary metabolic genes closely associated with protein synthesis and host metabolism. Dishui Lake phages shared more genes with uncultivated environmental viruses than with reference viruses based on the gene-sharing network classification. Phylogeny of proteomes and comparative genomics delineated three new genera within two known viral families of Kyanoviridae and Autographiviridae, and four new families in Caudoviricetes for these seven novel phages. Their potential hosts appeared to be from the dominant bacterial phyla in Dishui Lake. Altogether, our study provides initial insights into the composition and diversity of bacteriophage communities in Dishui Lake, contributing valuable knowledge to the ongoing research on the roles played by viruses in freshwater ecosystems.

The International Virus Bioinformatics Meeting 2023

The 2023 International Virus Bioinformatics Meeting was held in Valencia, Spain, from 24-26 May 2023, attracting approximately 180 participants worldwide. The primary objective of the conference was to establish a dynamic scientific environment conducive to discussion, collaboration, and the generation of novel research ideas. As the first in-person event following the SARS-CoV-2 pandemic, the meeting facilitated highly interactive exchanges among attendees. It served as a pivotal gathering for gaining insights into the current status of virus bioinformatics research and engaging with leading researchers and emerging scientists. The event comprised eight invited talks, 19 contributed talks, and 74 poster presentations across eleven sessions spanning three days. Topics covered included machine learning, bacteriophages, virus discovery, virus classification, virus visualization, viral infection, viromics, molecular epidemiology, phylodynamic analysis, RNA viruses, viral sequence analysis, viral surveillance, and metagenomics. This report provides rewritten abstracts of the presentations, a summary of the key research findings, and highlights shared during the meeting.

Identification of mobile genetic elements with geNomad

Identifying and characterizing mobile genetic elements in sequencing data is essential for understanding their diversity, ecology, biotechnological applications and impact on public health. Here we introduce geNomad, a classification and annotation framework that combines information from gene content and a deep neural network to identify sequences of plasmids and viruses. geNomad uses a dataset of more than 200,000 marker protein profiles to provide functional gene annotation and taxonomic assignment of viral genomes. Using a conditional random field model, geNomad also detects proviruses integrated into host genomes with high precision. In benchmarks, geNomad achieved high classification performance for diverse plasmids and viruses (Matthews correlation coefficient of 77.8% and 95.3%, respectively), substantially outperforming other tools. Leveraging geNomad's speed and scalability, we processed over 2.7 trillion base pairs of sequencing data, leading to the discovery of millions of viruses and plasmids that are available through the IMG/VR and IMG/PR databases. geNomad is available at https://portal.nersc.gov/genomad .

Identification and genomic analysis of temperate Halomonas bacteriophage vB_HmeY_H4907 from the surface sediment of the Mariana Trench at a depth of 8,900 m

Viruses play crucial roles in the ecosystem by modulating the host community structure, mediating biogeochemical cycles, and compensating for the metabolism of host cells. Mariana Trench, the world's deepest hadal habitat, harbors a variety of unique microorganisms that have adapted to its extreme conditions of low temperatures, high pressure, and nutrient scarcity. However, our knowledge about isolated hadal phage strains in the hadal trench is still limited. This study reported the discovery of a temperate phage, vB_HmeY_H4907, infecting Halomonas meridiana H4907, isolated from surface sediment from the Mariana Trench at a depth of 8,900 m. To our best knowledge, it is the deepest isolated siphovirus from the ocean. Its 40,452 bp linear dsDNA genome has 57.64% GC content and 55 open reading frames, and it is highly homologous to its host. Phylogenetic analysis and average nucleotide sequence identification reveal that vB_HmeY_H4907 is separated from the isolated phages and represents a new family, Suviridae, with eight predicted proviruses and six uncultured viral genomes. They are widely distributed in the ocean, suggesting a prevalence of this viral family in the deep sea. These findings expand our understanding of the phylogenetic diversity and genomic features of hadal lysogenic phages, provide essential information for further studies of phage-host interactions and evolution, and may reveal new insights into the lysogenic lifestyles of viruses inhabiting the hadal ocean. IMPORTANCE Halomonas phage vB_HmeY_H4907 is the deepest isolated siphovirus from the ocean, and it represents a novel abundant viral family in the ocean. This study provides insights into the genomic, phylogenetic, and ecological characteristics of the new viral family, namely, Suviridae.

An extended catalog of integrated prophages in the infant and adult fecal microbiome shows high prevalence of lysogeny

Background and aims

The acquisition and gradual maturation of gut microbial communities during early childhood is central to an individual's healthy development. Bacteriophages have the potential to shape the gut bacterial communities. However, the complex ecological interactions between phages and their bacterial host are still poorly characterized. In this study, we investigated the abundance and diversity of integrated prophages in infant and adult gut bacteria by detecting integrated prophages in metagenome assembled genomes (MAGs) of commensal bacteria.

Methods

Our study included 88 infants sampled at 3 weeks, 3 months, 6 months, and 12 months (n = 323 total samples), and their parents around delivery time (n = 138 total samples). Fecal DNA was extracted and characterized by using shotgun metagenomic sequencing, and a collection of prokaryotic MAGs was generated. The MAG collection was screened for the presence of integrated bacteriophage sequences, allowing their taxonomic and functional characterization.

Results

A large collection of 6,186 MAGs from infant and adult gut microbiota was obtained and screened for integrated prophages, allowing the identification of 7,165 prophage sequences longer than 10 kb. Strikingly, more than 70% of the near-complete MAGs were identified as lysogens. The prevalence of prophages in MAGs varied across bacterial families, with a lower prevalence observed among Coriobacteriaceae, Eggerthellaceae, Veillonellaceae and Burkholderiaceae, while a very high prevalence of lysogen MAGs were observed in Oscillospiraceae, Enterococcaceae, and Enterobacteriaceae. Interestingly for several bacterial families such as Bifidobacteriaceae and Bacteroidaceae, the prevalence of prophages in MAGs was higher in early infant time point (3 weeks and 3 months) than in later sampling points (6 and 12 months) and in adults. The prophage sequences were clustered into 5,616 species-like vOTUs, 77% of which were novel. Finally, we explored the functional repertoire of the potential auxiliary metabolic genes carried by these prophages, encoding functions involved in carbohydrate metabolism and degradation, amino acid metabolism and carbon metabolism.

Conclusion

Our study provides an enhanced understanding of the diversity and prevalence of lysogens in infant and adult gut microbiota and suggests a complex interplay between prophages and their bacterial hosts.