VirSorter: mining viral signal from microbial genomic data

Microbial diversity, genomics, and phage-host interactions of cyanobacterial harmful algal blooms

The occurrence of cyanobacterial harmful algal blooms (cyanoHABs) is related to their physical and chemical environment. However, less is known about their associated microbial interactions and processes. In this study, cyanoHABs were analyzed as a microbial ecosystem, using 1 year of 16S rRNA sequencing and 70 metagenomes collected during the bloom season from Lake Okeechobee (Florida, USA). Biogeographical patterns observed in microbial community composition and function reflected ecological zones distinct in their physical and chemical parameters that resulted in bloom "hotspots" near major lake inflows. Changes in relative abundances of taxa within multiple phyla followed increasing bloom severity. Functional pathways that correlated with increasing bloom severity encoded organic nitrogen and phosphorus utilization, storage of nutrients, exchange of genetic material, phage defense, and protection against oxidative stress, suggesting that microbial interactions may promote cyanoHAB resilience. Cyanobacterial communities were highly diverse, with picocyanobacteria ubiquitous and oftentimes most abundant, especially in the absence of blooms. The identification of novel bloom-forming cyanobacteria and genomic comparisons indicated a functionally diverse cyanobacterial community with differences in its capability to store nitrogen using cyanophycin and to defend against phage using CRISPR and restriction-modification systems. Considering blooms in the context of a microbial ecosystem and their interactions in nature, physiologies and interactions supporting the proliferation and stability of cyanoHABs are proposed, including a role for phage infection of picocyanobacteria. This study displayed the power of "-omics" to reveal important biological processes that could support the effective management and prediction of cyanoHABs.

IMPORTANCE

Cyanobacterial harmful algal blooms pose a significant threat to aquatic ecosystems and human health. Although physical and chemical conditions in aquatic systems that facilitate bloom development are well studied, there are fundamental gaps in the biological understanding of the microbial ecosystem that makes a cyanobacterial bloom. High-throughput sequencing was used to determine the drivers of cyanobacteria blooms in nature. Multiple functions and interactions important to consider in cyanobacterial bloom ecology were identified. The microbial biodiversity of blooms revealed microbial functions, genomic characteristics, and interactions between cyanobacterial populations that could be involved in bloom stability and more coherently define cyanobacteria blooms. Our results highlight the importance of considering cyanobacterial blooms as a microbial ecosystem to predict, prevent, and mitigate them.

DeepMicroClass sorts metagenomic contigs into prokaryotes, eukaryotes and viruses

Sequence classification facilitates a fundamental understanding of the structure of microbial communities. Binary metagenomic sequence classifiers are insufficient because environmental metagenomes are typically derived from multiple sequence sources. Here we introduce a deep-learning based sequence classifier, DeepMicroClass, that classifies metagenomic contigs into five sequence classes, i.e. viruses infecting prokaryotic or eukaryotic hosts, eukaryotic or prokaryotic chromosomes, and prokaryotic plasmids. DeepMicroClass achieved high performance for all sequence classes at various tested sequence lengths ranging from 500 bp to 100 kbps. By benchmarking on a synthetic dataset with variable sequence class composition, we showed that DeepMicroClass obtained better performance for eukaryotic, plasmid and viral contig classification than other state-of-the-art predictors. DeepMicroClass achieved comparable performance on viral sequence classification with geNomad and VirSorter2 when benchmarked on the CAMI II marine dataset. Using a coastal daily time-series metagenomic dataset as a case study, we showed that microbial eukaryotes and prokaryotic viruses are integral to microbial communities. By analyzing monthly metagenomes collected at HOT and BATS, we found relatively higher viral read proportions in the subsurface layer in late summer, consistent with the seasonal viral infection patterns prevalent in these areas. We expect DeepMicroClass will promote metagenomic studies of under-appreciated sequence types.

VIBES: a workflow for annotating and visualizing viral sequences integrated into bacterial genomes

Bacteriophages are viruses that infect bacteria. Many bacteriophages integrate their genomes into the bacterial chromosome and become prophages. Prophages may substantially burden or benefit host bacteria fitness, acting in some cases as parasites and in others as mutualists. Some prophages have been demonstrated to increase host virulence. The increasing ease of bacterial genome sequencing provides an opportunity to deeply explore prophage prevalence and insertion sites. Here we present VIBES (Viral Integrations in Bacterial genomES), a workflow intended to automate prophage annotation in complete bacterial genome sequences. VIBES provides additional context to prophage annotations by annotating bacterial genes and viral proteins in user-provided bacterial and viral genomes. The VIBES pipeline is implemented as a Nextflow-driven workflow, providing a simple, unified interface for execution on local, cluster and cloud computing environments. For each step of the pipeline, a container including all necessary software dependencies is provided. VIBES produces results in simple tab-separated format and generates intuitive and interactive visualizations for data exploration. Despite VIBES's primary emphasis on prophage annotation, its generic alignment-based design allows it to be deployed as a general-purpose sequence similarity search manager. We demonstrate the utility of the VIBES prophage annotation workflow by searching for 178 Pf phage genomes across 1072 Pseudomonas spp. genomes.

Phage phylogeny, molecular signaling, and auxiliary antimicrobial resistance in aerobic and anaerobic membrane bioreactors

Phage emit communication signals that inform their lytic and lysogenic life cycles. However, little is known regarding the abundance and diversity of the genes associated with phage communication systems in wastewater treatment microbial communities. This study focused on phage communities within two distinct biochemical wastewater environments, specifically aerobic membrane bioreactors (AeMBRs) and anaerobic membrane bioreactors (AnMBRs) exposed to varying antibiotic concentrations. Metagenomic data from the bench-scale systems were analyzed to explore phage phylogeny, life cycles, and genetic capacity for antimicrobial resistance and quorum sensing. Two dominant phage families, Schitoviridae and Peduoviridae, exhibited redox-dependent dynamics. Schitoviridae prevailed in anaerobic conditions, while Peduoviridae dominated in aerobic conditions. Notably, the abundance of lytic and lysogenic proteins varied across conditions, suggesting the coexistence of both life cycles. Furthermore, the presence of antibiotic resistance genes (ARGs) within viral contigs highlighted the potential for phage to transfer ARGs in AeMBRs. Finally, quorum sensing genes in the virome of AeMBRs indicated possible molecular signaling between phage and bacteria. Overall, this study provides insights into the dynamics of viral communities across varied redox conditions in MBRs. These findings shed light on phage life cycles, and auxiliary genetic capacity such as antibiotic resistance and bacterial quorum sensing within wastewater treatment microbial communities.

Genome sequence of the plant-growth-promoting bacterium Bacillus velezensis EU07

Many Gram-positive spore-forming rhizobacteria of the genus Bacillus show potential as biocontrol biopesticides that promise improved sustainability and ecological safety in agriculture. Here, we present a draft-quality genome sequence for Bacillus velezensis EU07, which shows growth-promotion in tomato plants and biocontrol against Fusarium head blight. We found that the genome of EU07 is almost identical to that of the commercially used strain QST713, but identified 46 single-nucleotide differences that distinguish these strains from each other. The availability of this genome sequence will facilitate future efforts to unravel the genetic and molecular basis for EU07's beneficial properties.

New avenues for potentially seeking microbial responses to climate change beneath Antarctic ice shelves

The signs of climate change are undeniable, and the impact of these changes on ecosystem function heavily depends on the response of microbes that underpin the food web. Antarctic ice shelf is a massive mass of floating ice that extends from the continent into the ocean, exerting a profound influence on global carbon cycles. Beneath Antarctic ice shelves, marine ice stores valuable genetic information, where marine microbial communities before the industrial revolution are archived. Here, in this proof-of-concept, by employing a combination of single-cell technologiesand metagenomics, we have been able to sequence frozen microbial DNA (≈300 years old) stored in the marine ice core B15 collected from the Filchnner-Ronne Ice Shelf. Metagenomic data indicated that Proteobacteria and Thaumarchaeota (e.g., Nitrosopumilus spp.), followed by Actinobacteria (e.g., Actinomarinales), were abundant. Remarkably, our data allow us to "travel to the past" and calibrate genomic and genetic evolutionary changes for ecologically relevant microbes and functions, such as Nitrosopumilus spp., preserved in the marine ice (≈300 years old) with those collected recently in seawater under an ice shelf (year 2017). The evolutionary divergence for the ammonia monooxygenase gene amoA involved in chemolithoautotrophy was about 0.88 amino acid and 2.8 nucleotide substitution rate per 100 sites in a century, while the accumulated rate of genomic SNPs was 2,467 per 1 Mb of genome and 100 years. Whether these evolutionary changes remained constant over the last 300 years or accelerated during post-industrial periods remains an open question that will be further elucidated.

IMPORTANCE

Several efforts have been undertaken to predict the response of microbes under climate change, mainly based on short-term microcosm experiments under forced conditions. A common concern is that manipulative experiments cannot properly simulate the response of microbes to climate change, which is a long-term evolutionary process. In this proof-of-concept study with a limited sample size, we demonstrate a novel approach yet to be fully explored in science for accessing genetic information from putative past marine microbes preserved under Antarctic ice shelves before the industrial revolution. This potentially allows us estimating evolutionary changes as exemplified in our study. We advocate for gathering a more comprehensive Antarctic marine ice core data sets across various periods and sites. Such a data set would enable the establishment of a robust baseline, facilitating a better assessment of the potential effects of climate change on key genetic signatures of microbes.

Long-read powered viral metagenomics in the oligotrophic Sargasso Sea

Dominant microorganisms of the Sargasso Sea are key drivers of the global carbon cycle. However, associated viruses that shape microbial community structure and function are not well characterised. Here, we combined short and long read sequencing to survey Sargasso Sea phage communities in virus- and cellular fractions at viral maximum (80 m) and mesopelagic (200 m) depths. We identified 2,301 Sargasso Sea phage populations from 186 genera. Over half of the phage populations identified here lacked representation in global ocean viral metagenomes, whilst 177 of the 186 identified genera lacked representation in genomic databases of phage isolates. Viral fraction and cell-associated viral communities were decoupled, indicating viral turnover occurred across periods longer than the sampling period of three days. Inclusion of long-read data was critical for capturing the breadth of viral diversity. Phage isolates that infect the dominant bacterial taxa Prochlorococcus and Pelagibacter, usually regarded as cosmopolitan and abundant, were poorly represented.

Another piece of the puzzle for the human microbiome: the gut virome under dietary modulation

The virome is the most abundant and highly variable microbial consortium in the gut. Because of difficulties in isolating and culturing gut viruses and consequently shorting of reference genomes, the virome has remained a relatively elusive aspect of the human microbiome. In recent years, studies on the virome have accumulated growing evidence showing that the virome is diet-modulated and widely involved in regulating health. Here, we review the responses of the gut virome to dietary intake and the potential health implications, presenting changes in the gut viral community, preferences of viral members to particular diets. We further discuss how viral-bacterial interactions and phage lifestyle shifts shape the gut microbiota. We also discuss the specific functions conferred by diet on the gut virome and bacterial community in the context of horizontal gene transfer, as well as the import of new viral members along with the diet. Collating these studies will expand our understanding of the dietary regulation of the gut virome and inspire dietary interventions and health maintenance strategies targeting the gut microbiota.

Phylogenomics and genetic analysis of solvent-producing Clostridium species

The genus Clostridium is a large and diverse group within the Bacillota (formerly Firmicutes), whose members can encode useful complex traits such as solvent production, gas-fermentation, and lignocellulose breakdown. We describe 270 genome sequences of solventogenic clostridia from a comprehensive industrial strain collection assembled by Professor David Jones that includes 194 C. beijerinckii, 57 C. saccharobutylicum, 4 C. saccharoperbutylacetonicum, 5 C. butyricum, 7 C. acetobutylicum, and 3 C. tetanomorphum genomes. We report methods, analyses and characterization for phylogeny, key attributes, core biosynthetic genes, secondary metabolites, plasmids, prophage/CRISPR diversity, cellulosomes and quorum sensing for the 6 species. The expanded genomic data described here will facilitate engineering of solvent-producing clostridia as well as non-model microorganisms with innately desirable traits. Sequences could be applied in conventional platform biocatalysts such as yeast or Escherichia coli for enhanced chemical production. Recently, gene sequences from this collection were used to engineer Clostridium autoethanogenum, a gas-fermenting autotrophic acetogen, for continuous acetone or isopropanol production, as well as butanol, butanoic acid, hexanol and hexanoic acid production.

Elucidating human gut microbiota interactions that robustly inhibit diverse Clostridioides difficile strains across different nutrient landscapes

The human gut pathogen Clostridioides difficile displays extreme genetic variability and confronts a changeable nutrient landscape in the gut. We mapped gut microbiota inter-species interactions impacting the growth and toxin production of diverse C. difficile strains in different nutrient environments. Although negative interactions impacting C. difficile are prevalent in environments promoting resource competition, they are sparse in an environment containing C. difficile-preferred carbohydrates. C. difficile strains display differences in interactions with Clostridium scindens and the ability to compete for proline. C. difficile toxin production displays substantial community-context dependent variation and does not trend with growth-mediated inter-species interactions. C. difficile shows substantial differences in transcriptional profiles in the presence of the closely related species C. hiranonis or C. scindens. In co-culture with C. hiranonis, C. difficile exhibits massive alterations in metabolism and other cellular processes, consistent with their high metabolic overlap. Further, Clostridium hiranonis inhibits the growth and toxin production of diverse C. difficile strains across different nutrient environments and ameliorates the disease severity of a C. difficile challenge in a murine model. In sum, strain-level variability and nutrient environments are major variables shaping gut microbiota interactions with C. difficile.

Pannonibacter anstelovis sp. nov. Isolated from Two Cases of Bloodstream Infections in Paediatric Patients

This study describes two cases of bacteraemia sustained by a new putative Pannonibacter species isolated at the U.O.C. of Microbiology and Virology of the Policlinico of Bari (Bari, Italy) from the blood cultures of two patients admitted to the Paediatric Oncohaematology Unit. Pannonibacter spp. is an environmental Gram-negative bacterium not commonly associated with nosocomial infections. Species identification was performed using Sanger sequencing of the 16S rRNA gene and Whole-Genome Sequencing (WGS) for both strains. Genomic analyses for the two isolates, BLAST similarity search, and phylogeny for the 16S rDNA sequences lead to an assignment to the species Pannonibacter phragmitetus. However, by performing ANIb, ANIm, tetranucleotide correlation, and DNA-DNA digital hybridization, analyses of the two draft genomes showed that they were very different from those of the species P. phragmitetus. MALDI-TOF analysis, assessment of antimicrobial susceptibility by E-test method, and Analytical Profile Index (API) tests were also performed. This result highlights how environmental bacterial species can easily adapt to the human host and, especially in nosocomial environments, also gain pathogenic potential through antimicrobial resistance.

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.

Characterization of the gut bacterial and viral microbiota in latent autoimmune diabetes in adults

Latent autoimmune diabetes in adults (LADA) is a heterogeneous disease characterized by autoantibodies against insulin producing pancreatic beta cells and initial lack of need for insulin treatment. The aim of the present study was to investigate if individuals with LADA have an altered gut microbiota relative to non-diabetic control subjects, individuals with type 1 diabetes (T1D), and individuals with type 2 diabetes (T2D). Bacterial community profiling was performed with primers targeting the variable region 4 of the 16S rRNA gene and sequenced. Amplicon sequence variants (ASVs) were generated with DADA2 and annotated to the SILVA database. The gut virome was sequenced, using a viral particle enrichment and metagenomics approach, assembled, and quantified to describe the composition of the viral community. Comparison of the bacterial alpha- and beta-diversity measures revealed that the gut bacteriome of individuals with LADA resembled that of individuals with T2D. Yet, specific genera were found to differ in abundance in individuals with LADA compared with T1D and T2D, indicating that LADA has unique taxonomical features. The virome composition reflected the stability of the most dominant order Caudovirales and the families Siphoviridae, Podoviridae, and Inoviridae, and the dominant family Microviridae. Further studies are needed to confirm these findings.

A compendium of ruminant gastrointestinal phage genomes revealed a higher proportion of lytic phages than in any other environments

BACKGROUND

Ruminants are important livestock animals that have a unique digestive system comprising multiple stomach compartments. Despite significant progress in the study of microbiome in the gastrointestinal tract (GIT) sites of ruminants, we still lack an understanding of the viral community of ruminants. Here, we surveyed its viral ecology using 2333 samples from 10 sites along the GIT of 8 ruminant species.

RESULTS

We present the Unified Ruminant Phage Catalogue (URPC), a comprehensive survey of phages in the GITs of ruminants including 64,922 non-redundant phage genomes. We characterized the distributions of the phage genomes in different ruminants and GIT sites and found that most phages were organism-specific. We revealed that ~ 60% of the ruminant phages were lytic, which was the highest as compared with those in all other environments and certainly will facilitate their applications in microbial interventions. To further facilitate the future applications of the phages, we also constructed a comprehensive virus-bacteria/archaea interaction network and identified dozens of phages that may have lytic effects on methanogenic archaea.

CONCLUSIONS

The URPC dataset represents a useful resource for future microbial interventions to improve ruminant production and ecological environmental qualities. Phages have great potential for controlling pathogenic bacterial/archaeal species and reducing methane emissions. Our findings provide insights into the virome ecology research of the ruminant GIT and offer a starting point for future research on phage therapy in ruminants. Video Abstract.

Phage-inclusive profiling of human gut microbiomes with Phanta

Due to technical limitations, most gut microbiome studies have focused on prokaryotes, overlooking viruses. Phanta, a virome-inclusive gut microbiome profiling tool, overcomes the limitations of assembly-based viral profiling methods by using customized k-mer-based classification tools and incorporating recently published catalogs of gut viral genomes. Phanta's optimizations consider the small genome size of viruses, sequence homology with prokaryotes and interactions with other gut microbes. Extensive testing of Phanta on simulated data demonstrates that it quickly and accurately quantifies prokaryotes and viruses. When applied to 245 fecal metagenomes from healthy adults, Phanta identifies ~200 viral species per sample, ~5× more than standard assembly-based methods. We observe a ~2:1 ratio between DNA viruses and bacteria, with higher interindividual variability of the gut virome compared to the gut bacteriome. In another cohort, we observe that Phanta performs equally well on bulk versus virus-enriched metagenomes, making it possible to study prokaryotes and viruses in a single experiment, with a single analysis.

Airborne antibiotic resistome and microbiome in pharmaceutical factories

Antimicrobial resistance is considered to be one of the biggest public health problems, and airborne transmission is an important but under-appreciated pathway for the spread of antibiotic resistance genes (ARGs) in the environment. Previous research has shown pharmaceutical factories to be a major source of ARGs and antibiotic resistant bacteria (ARB) in the surrounding receiving water and soil environments. Pharmaceutical factories are hotspots of antibiotic resistance, but the atmospheric transmission and its environmental risk remain more concerns. Here, we conducted a metagenomic investigation into the airborne microbiome and resistome in three pharmaceutical factories in China. Soil (average: 38.45%) and wastewater (average: 28.53%) were major contributors of airborne resistome. ARGs (vanR/vanS, blaOXA, and CfxA) conferring resistance to critically important clinically used antibiotics were identified in the air samples. The wastewater treatment area had significantly higher relative abundances of ARGs (average: 0.64 copies/16S rRNA). Approximately 28.2% of the detected airborne ARGs were found to be associated with plasmids, and this increased to about 50% in the wastewater treatment area. We have compiled a list of high-risk airborne ARGs found in pharmaceutical factories. Moreover, A total of 1,043 viral operational taxonomic units were identified and linked to 47 family-group taxa. Different CRISPR-Cas immune systems have been identified in bacterial hosts in response to phage infection. Similarly, higher phage abundance (average: 2451.70 PPM) was found in the air of the wastewater treatment area. Our data provide insights into the antibiotic resistance gene profiles and microbiome (bacterial and non-bacterial) in pharmaceutical factories and reveal the potential role of horizontal transfer in the spread of airborne ARGs, with implications for human and animal health.

Efficient Recovery of Complete Gut Viral Genomes by Combined Short- and Long-Read Sequencing

Current metagenome assembled human gut phage catalogs contained mostly fragmented genomes. Here, comprehensive gut virome detection procedure is developed involving virus-like particle (VLP) enrichment from ≈500 g feces and combined sequencing of short- and long-read. Applied to 135 samples, a Chinese Gut Virome Catalog (CHGV) is assembled consisting of 21,499 non-redundant viral operational taxonomic units (vOTUs) that are significantly longer than those obtained by short-read sequencing and contained ≈35% (7675) complete genomes, which is ≈nine times more than those in the Gut Virome Database (GVD, ≈4%, 1,443). Interestingly, the majority (≈60%, 13,356) of the CHGV vOTUs are obtained by either long-read or hybrid assemblies, with little overlap with those assembled from only the short-read data. With this dataset, vast diversity of the gut virome is elucidated, including the identification of 32% (6,962) novel vOTUs compare to public gut virome databases, dozens of phages that are more prevalent than the crAssphages and/or Gubaphages, and several viral clades that are more diverse than the two. Finally, the functional capacities are also characterized of the CHGV encoded proteins and constructed a viral-host interaction network to facilitate future research and applications.

Biogeographic patterns and drivers of soil viromes

Viruses are crucial in shaping soil microbial functions and ecosystems. However, studies on soil viromes have been limited in both spatial scale and biome coverage. Here we present a comprehensive synthesis of soil virome biogeographic patterns using the Global Soil Virome dataset (GSV) wherein we analysed 1,824 soil metagenomes worldwide, uncovering 80,750 partial genomes of DNA viruses, 96.7% of which are taxonomically unassigned. The biogeography of soil viral diversity and community structure varies across different biomes. Interestingly, the diversity of viruses does not align with microbial diversity and contrasts with it by showing low diversity in forest and shrubland soils. Soil texture and moisture conditions are further corroborated as key factors affecting diversity by our predicted soil viral diversity atlas, revealing higher diversity in humid and subhumid regions. In addition, the binomial degree distribution pattern suggests a random co-occurrence pattern of soil viruses. These findings are essential for elucidating soil viral ecology and for the comprehensive incorporation of viruses into soil ecosystem models.

Correlation between the gut microbiome and neurodegenerative diseases: a review of metagenomics evidence

A growing body of evidence suggests that the gut microbiota contributes to the development of neurodegenerative diseases via the microbiota-gut-brain axis. As a contributing factor, microbiota dysbiosis always occurs in pathological changes of neurodegenerative diseases, such as Alzheimer's disease, Parkinson's disease, and amyotrophic lateral sclerosis. High-throughput sequencing technology has helped to reveal that the bidirectional communication between the central nervous system and the enteric nervous system is facilitated by the microbiota's diverse microorganisms, and for both neuroimmune and neuroendocrine systems. Here, we summarize the bioinformatics analysis and wet-biology validation for the gut metagenomics in neurodegenerative diseases, with an emphasis on multi-omics studies and the gut virome. The pathogen-associated signaling biomarkers for identifying brain disorders and potential therapeutic targets are also elucidated. Finally, we discuss the role of diet, prebiotics, probiotics, postbiotics and exercise interventions in remodeling the microbiome and reducing the symptoms of neurodegenerative diseases.

Isolation of viruses, including mollivirus, with the potential to infect Acanthamoeba from a Japanese warm temperate zone

Acanthamoeba castellanii is infected with diverse nucleocytoplasmic large DNA viruses. Here, we report the co-isolation of 12 viral strains from marine sediments in Uranouchi Inlet, Kochi, Japan. Based on the morphological features revealed by electron microscopy, these isolates were classified into four viral groups including Megamimiviridae, Molliviridae, Pandoraviridae, and Pithoviridae. Genomic analyses indicated that these isolates showed high similarities to the known viral genomes with which they are taxonomically clustered, and their phylogenetic relationships were also supported by core gene similarities. It is noteworthy that Molliviridae was isolated from the marine sediments in the Japanese warm temperate zone because other strains have only been found in the subarctic region. Furthermore, this strain has 19 and 4 strain-specific genes found in Mollivirus sibericum and Mollivirus kamchatka, respectively. This study extends our knowledge about the habitat and genomic diversity of Molliviridae.

Advances in phage-host interaction prediction: in silico method enhances the development of phage therapies

Phages can specifically recognize and kill bacteria, which lead to important application value of bacteriophage in bacterial identification and typing, livestock aquaculture and treatment of human bacterial infection. Considering the variety of human-infected bacteria and the continuous discovery of numerous pathogenic bacteria, screening suitable therapeutic phages that are capable of infecting pathogens from massive phage databases has been a principal step in phage therapy design. Experimental methods to identify phage-host interaction (PHI) are time-consuming and expensive; high-throughput computational method to predict PHI is therefore a potential substitute. Here, we systemically review bioinformatic methods for predicting PHI, introduce reference databases and in silico models applied in these methods and highlight the strengths and challenges of current tools. Finally, we discuss the application scope and future research direction of computational prediction methods, which contribute to the performance improvement of prediction models and the development of personalized phage therapy.

Tracking and characterization of a novel conjugative transposon identified by shotgun transposon mutagenesis

The horizontal transfer of mobile genetic elements (MGEs) is an essential process determining the functional and genomic diversity of bacterial populations. MGEs facilitate the exchange of fitness determinant genes like antibiotic resistance and virulence factors. Various computational methods exist to identify potential MGEs, but confirming their ability to transfer requires additional experimental approaches. Here, we apply a transposon (Tn) mutagenesis technique for confirming mobilization without the need for targeted mutations. Using this method, we identified two MGEs, including a previously known conjugative transposon (CTn) called BoCTn found in Bacteroides ovatus and a novel CTn, PvCTn, identified in Phocaeicola vulgatus. In addition, Tn mutagenesis and subsequent genetic deletion enabled our characterization of a helix-turn-helix motif gene, BVU3433 which negatively regulates the conjugation efficiency of PvCTn in vitro. Furthermore, our transcriptomics data revealed that BVU3433 plays a crucial role in the repression of PvCTn genes, including genes involved in forming complete conjugation machinery [Type IV Secretion System (T4SS)]. Finally, analysis of individual strain genomes and community metagenomes identified the widespread prevalence of PvCTn-like elements with putative BVU3433 homologs among human gut-associated bacteria. In summary, this Tn mutagenesis mobilization method (TMMM) enables observation of transfer events in vitro and can ultimately be applied in vivo to identify a broader diversity of functional MGEs that may underly the transfer of important fitness determinants.

Benchmarking informatics approaches for virus discovery: caution is needed when combining in silico identification methods

Understanding the ecological impacts of viruses on natural and engineered ecosystems relies on the accurate identification of viral sequences from community sequencing data. To maximize viral recovery from metagenomes, researchers frequently combine viral identification tools. However, the effectiveness of this strategy is unknown. Here, we benchmarked combinations of six widely used informatics tools for viral identification and analysis (VirSorter, VirSorter2, VIBRANT, DeepVirFinder, CheckV, and Kaiju), called "rulesets." Rulesets were tested against mock metagenomes composed of taxonomically diverse sequence types and diverse aquatic metagenomes to assess the effects of the degree of viral enrichment and habitat on tool performance. We found that six rulesets achieved equivalent accuracy [Matthews Correlation Coefficient (MCC) = 0.77, Padj ≥ 0.05]. Each contained VirSorter2, and five used our "tuning removal" rule designed to remove non-viral contamination. While DeepVirFinder, VIBRANT, and VirSorter were each found once in these high-accuracy rulesets, they were not found in combination with each other: combining tools does not lead to optimal performance. Our validation suggests that the MCC plateau at 0.77 is partly caused by inaccurate labeling within reference sequence databases. In aquatic metagenomes, our highest MCC ruleset identified more viral sequences in virus-enriched (44%-46%) than in cellular metagenomes (7%-19%). While improved algorithms may lead to more accurate viral identification tools, this should be done in tandem with careful curation of sequence databases. We recommend using the VirSorter2 ruleset and our empirically derived tuning removal rule. Our analysis provides insight into methods for in silico viral identification and will enable more robust viral identification from metagenomic data sets.

IMPORTANCE

The identification of viruses from environmental metagenomes using informatics tools has offered critical insights in microbial ecology. However, it remains difficult for researchers to know which tools optimize viral recovery for their specific study. In an attempt to recover more viruses, studies are increasingly combining the outputs from multiple tools without validating this approach. After benchmarking combinations of six viral identification tools against mock metagenomes and environmental samples, we found that these tools should only be combined cautiously. Two to four tool combinations maximized viral recovery and minimized non-viral contamination compared with either the single-tool or the five- to six-tool ones. By providing a rigorous overview of the behavior of in silico viral identification strategies and a pipeline to replicate our process, our findings guide the use of existing viral identification tools and offer a blueprint for feature engineering of new tools that will lead to higher-confidence viral discovery in microbiome studies.

In Vitro Microevolution and Co-Selection Assessment of Amoxicillin and Cefotaxime Impact on Escherichia coli Resistance Development

The global spread of antimicrobial resistance has become a prominent issue in both veterinary and public health in the 21st century. The extensive use of amoxicillin, a beta-lactam antibiotic, and consequent resistance development are particularly alarming in food-producing animals, with a focus on the swine and poultry sectors. Another beta-lactam, cefotaxime, is widely utilized in human medicine, where the escalating resistance to third- and fourth-generation cephalosporins is a major concern. The aim of this study was to simulate the development of phenotypic and genotypic resistance to beta-lactam antibiotics, focusing on amoxicillin and cefotaxime. The investigation of the minimal inhibitory concentrations (MIC) of antibiotics was performed at 1×, 10×, 100×, and 1000× concentrations using the modified microbial evolution and growth arena (MEGA-plate) method. Our results indicate that amoxicillin significantly increased the MIC values of several tested antibiotics, except for oxytetracycline and florfenicol. In the case of cefotaxime, this increase was observed in all classes. A total of 44 antimicrobial resistance genes were identified in all samples. Chromosomal point mutations, particularly concerning cefotaxime, revealed numerous complex mutations, deletions, insertions, and single nucleotide polymorphisms (SNPs) that were not experienced in the case of amoxicillin. The findings suggest that, regarding amoxicillin, the point mutation of the acrB gene could explain the observed MIC value increases due to the heightened activity of the acrAB-tolC efflux pump system. However, under the influence of cefotaxime, more intricate processes occurred, including complex amino acid substitutions in the ampC gene promoter region, increased enzyme production induced by amino acid substitutions and SNPs, as well as mutations in the acrR and robA repressor genes that heightened the activity of the acrAB-tolC efflux pump system. These changes may contribute to the significant MIC increases observed for all tested antibiotics. The results underscore the importance of understanding cross-resistance development between individual drugs when choosing clinical alternative drugs. The point mutations in the mdtB and emrR genes may also contribute to the increased activity of the mdtABC-tolC and emrAB-tolC pump systems against all tested antibiotics. The exceptionally high mutation rate induced by cephalosporins justifies further investigations to clarify the exact mechanism behind.

Transmission and dynamics of mother-infant gut viruses during pregnancy and early life

Early development of the gut ecosystem is crucial for lifelong health. While infant gut bacterial communities have been studied extensively, the infant gut virome remains under-explored. To study the development of the infant gut virome over time and the factors that shape it, we longitudinally assess the composition of gut viruses and their bacterial hosts in 30 women during and after pregnancy and in their 32 infants during their first year of life. Using shotgun metagenomic sequencing applied to dsDNA extracted from Virus-Like Particles (VLPs) and bacteria, we generate 205 VLP metaviromes and 322 total metagenomes. With this data, we show that while the maternal gut virome composition remains stable during late pregnancy and after birth, the infant gut virome is dynamic in the first year of life. Notably, infant gut viromes contain a higher abundance of active temperate phages compared to maternal gut viromes, which decreases over the first year of life. Moreover, we show that the feeding mode and place of delivery influence the gut virome composition of infants. Lastly, we provide evidence of co-transmission of viral and bacterial strains from mothers to infants, demonstrating that infants acquire some of their virome from their mother's gut.

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

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

Environmental viromes reveal the global distribution signatures of deep-sea DNA viruses

INTRODUCTION

Viruses are abundant and ecologically significant in marine ecosystems. However, the virome of deep-sea sediments is not extensively investigated.

OBJECTIVES

To explore the distribution pattern of deep-sea viruses on a global scale, the viromes of DNA viruses isolated from 138 sediments of 5 deep-sea ecosystems were characterized.

METHODS

The viral particles were purified from each sediment sample. Then the viral DNAs were extracted and subjected to viral metagenomic analysis.

RESULTS

Here, we constructed a global deep-sea environmental virome dataset by analyzing the viral DNA of 138 sediment samples. A total of 347,737 viral operational taxonomic units (vOTUs) were identified, of which 84.94% were hitherto unknown, indicating that deep sea was a reservoir of novel DNA viruses. Furthermore, circular viral genome analysis revealed 98,581 complete genomes. The classified vOTUs included eukaryotic (44.55%) and prokaryotic (25.75%) viruses, and were taxonomically assigned to 63 viral families. The composition and abundance of the deep-sea sediment viromes were dependent on the deep-sea ecosystem as opposed to geographical region. Further analysis revealed that the viral community differentiation in different deep-sea ecosystems was driven by the virus-mediated energy metabolism.

CONCLUSION

Our findings showed that deep-sea ecosystems are a reservoir of novel DNA viruses and the viral community is shaped by the environmental characteristics of deep-sea ecosystems, thus presenting critical information for determining the ecological significance of viruses in global deep-sea ecosystems.

Bacteriophages in nature: recent advances in research tools and diverse environmental and biotechnological applications

Bacteriophages infect and replicate within bacteria and play a key role in the environment, particularly in microbial ecosystems and bacterial population dynamics. The increasing recognition of their significance stems from their wide array of environmental and biotechnological uses, which encompass the mounting issue of antimicrobial resistance (AMR). Beyond their therapeutic potential in combating antibiotic-resistant infections, bacteriophages also find vast applications such as water quality monitoring, bioremediation, and nutrient cycling within environmental sciences. Researchers are actively involved in isolating and characterizing bacteriophages from different natural sources to explore their applications. Gaining insights into key aspects such as the life cycle of bacteriophages, their host range, immune interactions, and physical stability is vital to enhance their application potential. The establishment of diverse phage libraries has become indispensable to facilitate their wide-ranging uses. Consequently, numerous protocols, ranging from traditional to cutting-edge techniques, have been developed for the isolation, detection, purification, and characterization of bacteriophages from diverse environmental sources. This review offers an exploration of tools, delves into the methods of isolation, characterization, and the extensive environmental applications of bacteriophages, particularly in areas like water quality assessment, the food sector, therapeutic interventions, and the phage therapy in various infections and diseases.

Identification and characterization of integrated prophages and CRISPR-Cas system in Bacillus subtilis RS10 genome

Bacteriophages have been extensively investigated due to their prominent role in the virulence and resistance of pathogenic bacteria. However, little attention has been given to the non-pathogenic Bacillus phages, and their role in the ecological bacteria genome is overlooked. In the present study, we characterized two Bacillus phages with a linear DNA genome of 33.6 kb with 44.83% GC contents and 129.3 kb with 34.70% GC contents. A total of 46 and 175 putative coding DNA sequences (CDS) were identified in prophage 1 (P1) and prophage 2 (P2), respectively, with no tRNA genes. Comparative genome sequence analysis revealed that P1 shares eight CDS with phage Jimmer 2 (NC-041976), and phage Osiris (NC-028969), and six with phage phi CT9441A (NC-029022). On the other hand, P2 showed high similarity with Bacill_SPbeta_NC_001884 and Bacillus phage phi 105. Further, genome analysis indicates several horizontal gene transfer events in both phages during the evolution process. In addition, we detected two CRISPR-Cas systems for the first time in B. subtilis. The identified CRISPR system consists of 24 and 25 direct repeats and integrase coding genes, while the cas gene which encodes Cas protein involved in the cleavage of a target sequence is missing. These findings will expand the current knowledge of soil phages as well as help to develop a new perspective for investigating more ecological phages to understand their role in bacterial communities and diversity.

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.

Virus-prokaryote infection pairs associated with prokaryotic production in a freshwater lake

Viruses infect and kill prokaryotic populations in a density- or frequency-dependent manner and affect carbon cycling. However, the effects of the stratification transition, including the stratified and de-stratified periods, on the changes in prokaryotic and viral communities and their interactions remain unclear. We conducted a monthly survey of the surface and deep layers of a large and deep freshwater lake (Lake Biwa, Japan) for a year and analyzed the prokaryotic production and prokaryotic and viral community composition. Our analysis revealed that, in the surface layer, 19 prokaryotic species, accounting for approximately 40% of the total prokaryotic abundance, could potentially contribute to the majority of prokaryotic production, which is the highest during the summer and is suppressed by viruses. This suggests that a small fraction of prokaryotes and phages were the key infection pairs during the peak period of prokaryotic activity in the freshwater lake. We also found that approximately 50% of the dominant prokaryotic and viral species in the deep layer were present throughout the study period. This suggests that the "kill the winner" model could explain the viral impact on prokaryotes in the surface layer, but other dynamics may be at play in the deep layer. Furthermore, we found that annual vertical mixing could result in a similar rate of community change between the surface and deep layers. These findings may be valuable in understanding how communities and the interaction among them change when freshwater lake stratification is affected by global warming in the future.IMPORTANCEViral infection associated with prokaryotic production occurs in a density- or frequency-dependent manner and regulates the prokaryotic community. Stratification transition and annual vertical mixing in freshwater lakes are known to affect the prokaryotic community and the interaction between prokaryotes and viruses. By pairing measurements of virome analysis and prokaryotic production of a 1-year survey of the depths of surface and deep layers, we revealed (i) the prokaryotic infection pairs associated with prokaryotic production and (ii) the reset in prokaryotic and viral communities through annual vertical mixing in a freshwater lake. Our results provide a basis for future work into changes in stratification that may impact the biogeochemical cycling in freshwater lakes.

Discovering and exploring the hidden diversity of human gut viruses using highly enriched virome samples

Viruses are an abundant and crucial component of the human microbiome, but accurately discovering them via metagenomics is still challenging. Currently, the available viral reference genomes poorly represent the diversity in microbiome samples, and expanding such a set of viral references is difficult. As a result, many viruses are still undetectable through metagenomics even when considering the power of de novo metagenomic assembly and binning, as viruses lack universal markers. Here, we describe a novel approach to catalog new viral members of the human gut microbiome and show how the resulting resource improves metagenomic analyses. We retrieved >3,000 viral-like particles (VLP) enriched metagenomic samples (viromes), evaluated the efficiency of the enrichment in each sample to leverage the viromes of highest purity, and applied multiple analysis steps involving assembly and comparison with hundreds of thousands of metagenome-assembled genomes to discover new viral genomes. We reported over 162,000 viral sequences passing quality control from thousands of gut metagenomes and viromes. The great majority of the retrieved viral sequences (~94.4%) were of unknown origin, most had a CRISPR spacer matching host bacteria, and four of them could be detected in >50% of a set of 18,756 gut metagenomes we surveyed. We included the obtained collection of sequences in a new MetaPhlAn 4.1 release, which can quantify reads within a metagenome matching the known and newly uncovered viral diversity. Additionally, we released the viral database for further virome and metagenomic studies of the human microbiome.

Large language models improve annotation of prokaryotic viral proteins

Viral genomes are poorly annotated in metagenomic samples, representing an obstacle to understanding viral diversity and function. Current annotation approaches rely on alignment-based sequence homology methods, which are limited by the paucity of characterized viral proteins and divergence among viral sequences. Here we show that protein language models can capture prokaryotic viral protein function, enabling new portions of viral sequence space to be assigned biologically meaningful labels. When applied to global ocean virome data, our classifier expanded the annotated fraction of viral protein families by 29%. Among previously unannotated sequences, we highlight the identification of an integrase defining a mobile element in marine picocyanobacteria and a capsid protein that anchors globally widespread viral elements. Furthermore, improved high-level functional annotation provides a means to characterize similarities in genomic organization among diverse viral sequences. Protein language models thus enhance remote homology detection of viral proteins, serving as a useful complement to existing approaches.

Environmental viromes reveal global virosphere of deep-sea sediment RNA viruses

INTRODUCTION

Viruses are the most abundant and diverse life forms on the earth. Both DNA viruses and RNA viruses play important roles in marine ecosystems via regulating biogeochemical cycles.

OBJECTIVES

However, the virome of marine RNA viruses has been rarely explored so far. In this study, therefore, the environmental viromes of deep-sea sediment RNA viruses were characterized on a global scale to reveal the global virosphere of deep-sea RNA viruses.

METHODS

The viral particles were purified from each of 133 deep-sea sediment samples and then characterized based on metagenomes of RNA viruses.

RESULTS

In this study, we established the global virome dataset of deep-sea RNA viruses purified from 133 sediment samples that were collected from typical deep-sea ecosystems of three oceans. A total of 85,059 viral operational taxonomic units (vOTUs) were identified, of which only 1.72% were hitherto known, indicating that the deep-sea sediment is a repository of novel RNA viruses. These vOTUs were classified into 20 viral families, including prokaryotic (7.09%) and eukaryotic (65.81%) RNA viruses. Furthermore, 1,463 deep-sea RNA viruses with complete genomes were obtained. The differentiation of RNA viral communities was driven by the deep-sea ecosystems as opposed to geographical region. Specifically, the virus-encoded metabolic genes took great effects on the differentiation of RNA viral communities by mediating the energy metabolism in the deep-sea ecosystems.

CONCLUSIONS

Therefore, our findings indicate that the deep sea is a vast reservoir of novel RNA viruses for the first time, and the differentiation of RNA viral communities is driven by the deep-sea ecosystems through energy metabolism.

Interrogating two extensively self-targeting Type I CRISPR-Cas systems in Xanthomonas albilineans reveals distinct anti-CRISPR proteins that block DNA degradation

CRISPR-Cas systems store fragments of invader DNA as spacers to recognize and clear those same invaders in the future. Spacers can also be acquired from the host's genomic DNA, leading to lethal self-targeting. While self-targeting can be circumvented through different mechanisms, natural examples remain poorly explored. Here, we investigate extensive self-targeting by two CRISPR-Cas systems encoding 24 self-targeting spacers in the plant pathogen Xanthomonas albilineans. We show that the native I-C and I-F1 systems are actively expressed and that CRISPR RNAs are properly processed. When expressed in Escherichia coli, each Cascade complex binds its PAM-flanked DNA target to block transcription, while the addition of Cas3 paired with genome targeting induces cell killing. While exploring how X. albilineans survives self-targeting, we predicted putative anti-CRISPR proteins (Acrs) encoded within the bacterium's genome. Screening of identified candidates with cell-free transcription-translation systems and in E. coli revealed two Acrs, which we named AcrIC11 and AcrIF12Xal, that inhibit the activity of Cas3 but not Cascade of the respective system. While AcrF12Xal is homologous to AcrIF12, AcrIC11 shares sequence and structural homology with the anti-restriction protein KlcA. These findings help explain tolerance of self-targeting through two CRISPR-Cas systems and expand the known suite of DNA degradation-inhibiting Acrs.

Metagenomic analysis of antibiotic-resistance genes and viruses released from glaciers into downstream habitats

Glaciers serve as effective reservoirs of antibiotic resistance genes (ARGs) and viruses for millions of years. Climate change and anthropogenic activity have accelerated the melting of glaciers, but the patterns of release of ARGs and viruses from melting glaciers into downstream habitats remain unknown. We analyzed 171 metagenomic samples from glaciers and their downstream habitats and found that the abundance and diversity of ARGs were higher in glaciers (polar and plateau glaciers) than downstream habitats (Arctic Ocean, Qinghai Lake, and Yangtze River Basin), with the diversity of viruses having the opposite pattern. Proteobacteria and Actinobacteria were the main potential hosts of ARGs and viruses, and the richness of ARGs carried by the hosts was positively correlated with viral abundance, suggesting that the transmission of viruses in the hosts could disseminate ARGs. Source tracking indicated that >18 % of the ARGs and >25 % of the viruses detected in downstream habitats originated from glaciers, demonstrating that glaciers could be one of the potential sources of ARGs and viruses in downstream habitats. Increased solar radiation and emission of carbon dioxide mainly influenced the release of the ARGs and viruses from glaciers into downstream habitats. This study provides a systematic insight demonstrating the release of ARGs and viruses from the melting glaciers, potentially increasing ecological pressure.

Diversities and interactions of phages and bacteria in deep-sea sediments as revealed by metagenomics

Phages are found virtually everywhere, even in extreme environments, and are extremely diverse both in their virion structures and in their genomic content. They are thought to shape the taxonomic and functional composition of microbial communities as well as their stability. A number of studies on laboratory culture and viral metagenomic research provide deeper insights into the abundance, diversity, distribution, and interaction with hosts of phages across a wide range of ecosystems. Although most of these studies focus on easily accessible samples, such as soils, lakes, and shallow oceans, little is known about bathypelagic phages. In this study, through analyzing the 16S rRNA sequencing and viral metagenomic sequencing data of 25 samples collected from five different bathypelagic ecosystems, we detected a high diversity of bacteria and phages, particularly in the cold seep and hydrothermal vent ecosystems, which have stable chemical energy. The relative abundance of phages in these ecosystems was higher than in other three abyssal ecosystems. The low phage/host ratios obtained from host prediction were different from shallow ecosystems and indicated the prevalence of prophages, suggesting the complexity of phage-bacteria interactions in abyssal ecosystems. In the correlation analysis, we revealed several phages-bacteria interaction networks of potential ecological relevance. Our study contributes to a better understanding of the interactions between bathypelagic bacteria and their phages.

Evidence of interspecific plasmid uptake by pathogenic strains of Klebsiella isolated from microplastic pollution on public beaches

Microplastic beads are becoming a common feature on beaches, and there is increasing evidence that such microplastics can become colonised by potential human pathogens. However, whether the concentrations and pathogenicity of these pathogens pose a public health risk are still unclear. Therefore, the aim of this study was to determine realistic environmental concentrations of potential pathogens colonising microplastic beads, and quantify the expression of virulence and antimicrobial resistance genes (ARGs). Microplastic beads were collected from beaches and a culture-dependent approach was used to determine the concentrations of seven target bacteria (Campylobacter spp.; E. coli; intestinal enterococci; Klebsiella spp.; Pseudomonas aeruginosa; Salmonella spp.; Vibrio spp.). All seven target bacteria were detected without the need for a pre-enrichment step; urban sites had higher bacterial concentrations, whilst polymer type had no influence on bacterial concentrations. Klebsiella was the most abundant target bacteria and possessed virulence and ARGs, some of which were present on plasmids from other species, and showed pathogenicity in a Galleria melonella infection model. Our findings demonstrate how pathogen colonised microplastic beads can pose a heightened public health risk at the beach, and highlights the urgency for improved monitoring and enforcement of regulations on the release of microplastics into the environment.

Novel strains of Culex flavivirus and Hubei chryso-like virus 1 from the Anopheles mosquito in western Kenya

Surveillance of mosquito vectors is critical for early detection, prevention and control of vector borne diseases. In this study we used advanced molecular tools, such as DNA barcoding in combination with novel sequencing technologies to discover new and already known viruses in genetically identified mosquito species. Mosquitoes were captured using BG sentinel traps in Western Kenya during May and July 2019, and homogenized individually before pooled into groups of ten mosquitoes. The pools and individual samples were then used for molecular analysis and to infect cell cultures. Of a total of fifty-four (54) 10-pools, thirteen (13) showed cytopathic effect (CPE) on VeroB4 cells, eighteen (18) showed CPE on C6/36 cells. Eight (8) 10-pools out of the 31 CPE positive pools showed CPE on both VeroB4 and C6/36 cells. When using reverse transcriptase polymerase chain reaction (RT-PCR), Sanger sequencing and Twist Comprehensive Viral Research Panel (CVRP) (Twist Biosciences), all pools were found negative by RT-PCR when using genus specific primers targeting alphaviruses, orthobunyaviruses and virus specific primers towards o'nyong-nyong virus, chikungunya virus and Sindbis virus (previously reported to circulate in the region). Interestingly, five pools were RT-PCR positive for flavivirus. Two of the RT-PCR positive pools showed CPE on both VeroB4 and C6/36 cells, two pools showed CPE on C6/36 cells alone and one pool on VeroB4 cells only. Fifty individual mosquito homogenates from the five RT-PCR positive 10-pools were analyzed further for flavivirus RNA. Of these, 19 out of the 50 individual mosquito homogenates indicated the presence of flavivirus RNA. Barcoding of the flavivirus positive mosquitoes revealed the mosquito species as Aedes aegypti (1), Mansonia uniformis (6), Anopheles spp (3), Culex pipiens (5), Culex spp (1), Coquilletidia metallica (2) and Culex quinquefasciatus (1). Of the 19 flavivirus positive individual mosquitoes, five (5) virus positive homogenates were sequenced. Genome sequences of two viruses were completed. One was identified as the single-stranded RNA Culex flavivirus and the other as the double-stranded RNA Hubei chryso-like virus 1. Both viruses were found in the same Anopheles spp. homogenate extracted from a sample that showed CPE on both VeroB4 and C6/36 cells. The detection of both viruses in a single mosquito homogenate indicated coinfection. Phylogenetic analyses suggested that the Culex flavivirus sequence detected was closely related to a Culex flavivirus isolated from Uganda in 2008. All four Hubei chryso-like virus 1 segments clusters closely to Hubei chryso-like virus 1 strains isolated in Australia, China and USA. Two novel strains of insect-specific viruses in Anopheles mosquitoes were detected and characterized.

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.

Freshwater Viral Metagenome Analyses Targeting dsDNA Viruses

Viral metagenomics is one of the most widely used approaches to study viral population genomics. With the recent development of bioinformatic tools, the number of molecular biological methods, programs, and software to analyze viral metagenome data have greatly increased. Here, we describe the basic analysis workflow along with bioinformatic tools that can be used to analyze viral metagenome data. Although this chapter assumes that the viral metagenome data are prepared from the freshwater samples and are subjected to dsDNA sequencing, the protocol can be applied and modified for other types of metagenome data collected from a variety of sources.

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.

Forecasting the dynamics of a complex microbial community using integrated meta-omics

Predicting the behaviour of complex microbial communities is challenging. However, this is essential for complex biotechnological processes such as those in biological wastewater treatment plants (BWWTPs), which require sustainable operation. Here we summarize 14 months of longitudinal meta-omics data from a BWWTP anaerobic tank into 17 temporal signals, explaining 91.1% of the temporal variance, and link those signals to ecological events within the community. We forecast the signals over the subsequent five years and use 21 extra samples collected at defined time intervals for testing and validation. Our forecasts are correct for six signals and hint on phenomena such as predation cycles. Using all the 17 forecasts and the environmental variables, we predict gene abundance and expression, with a coefficient of determination ≥0.87 for the subsequent three years. Our study demonstrates the ability to forecast the dynamics of open microbial ecosystems using interactions between community cycles and environmental parameters.

Identifying ARG-carrying bacteriophages in a lake replenished by reclaimed water using deep learning techniques

As important mobile genetic elements, phages support the spread of antibiotic resistance genes (ARGs). Previous analyses of metaviromes or metagenome-assembled genomes (MAGs) failed to assess the extent of ARGs transferred by phages, particularly in the generation of antibiotic pathogens. Therefore, we have developed a bioinformatic pipeline that utilizes deep learning techniques to identify ARG-carrying phages and predict their hosts, with a special focus on pathogens. Using this method, we discovered that the predominant types of ARGs carried by temperate phages in a typical landscape lake, which is fully replenished by reclaimed water, were related to multidrug resistance and β-lactam antibiotics. MAGs containing virulent factors (VFs) were predicted to serve as hosts for these ARG-carrying phages, which suggests that the phages may have the potential to transfer ARGs. In silico analysis showed a significant positive correlation between temperate phages and host pathogens (R = 0.503, p < 0.001), which was later confirmed by qPCR. Interestingly, these MAGs were found to be more abundant than those containing both ARGs and VFs, especially in December and March. Seasonal variations were observed in the abundance of phages harboring ARGs (from 5.62 % to 21.02 %) and chromosomes harboring ARGs (from 18.01 % to 30.94 %). In contrast, the abundance of plasmids harboring ARGs remained unchanged. In summary, this study leverages deep learning to analyze phage-transferred ARGs and demonstrates an alternative method to track the production of potential antibiotic-resistant pathogens by metagenomics that can be extended to microbiological risk assessment.

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

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

Pan-metagenome reveals the abiotic stress resistome of cigar tobacco phyllosphere microbiome

The important role of microbial associations in mediating plant protection and responses to abiotic stresses has been widely recognized. However, there have been limited studies on the functional profile of the phyllosphere microbiota from tobacco (Nicotiana tabacum), hindering our understanding of the mechanisms underlying stress resilience in this representative and easy-to-cultivate model species from the solanaceous family. To address this knowledge gap, our study employed shotgun metagenomic sequencing for the first time to analyze the genetic catalog and identify putative plant growth promoting bacteria (PGPB) candidates that confer abiotic stress resilience throughout the growth period of cigar tobacco in the phyllosphere. We identified abundant genes from specific bacterial lineages, particularly Pseudomonas, within the cigar tobacco phyllospheric microbiome. These genes were found to confer resilience against a wide range of stressors, including osmotic and drought stress, heavy metal toxicity, temperature perturbation, organic pollutants, oxidative stress, and UV light damage. In addition, we conducted a virome mining analysis on the metagenome to explore the potential roles of viruses in driving microbial adaptation to environmental stresses. Our results identified a total of 3,320 scaffolds predicted to be viral from the cigar tobacco phyllosphere metagenome, with various phages infecting Pseudomonas, Burkholderia, Enterobacteria, Ralstonia, and related viruses. Within the virome, we also annotated genes associated with abiotic stress resilience, such as alkaline phosphatase D (phoD) for nutrient solubilization and glutamate-5-semialdehyde dehydrogenase (proA) for osmolyte synthesis. These findings shed light on the unexplored roles of viruses in facilitating and transferring abiotic stress resilience in the phyllospheric microbiome through beneficial interactions with their hosts. The findings from this study have important implications for agricultural practices, as they offer potential strategies for harnessing the capabilities of the phyllosphere microbiome to enhance stress tolerance in crop plants.

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

IMPORTANCE

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

In Vitro Microevolution and Co-Selection Assessment of Florfenicol Impact on Escherichia coli Resistance Development

The issue of antimicrobial resistance is becoming an increasingly serious challenge in both human and veterinary medicine. Prudent antimicrobial use in veterinary medicine is warranted and supported by international guidelines, with the Antimicrobial Advice Ad Hoc Expert Group (AMEG) placing particular emphasis on the critically important group B antimicrobials. These antimicrobials are commonly employed, especially in the poultry and swine industry. The impact of florfenicol, a veterinary antibiotic, was studied on the resistance development of Escherichia coli. The aim of the study was to investigate the effect of the use of florfenicol on the development of phenotypic and genomic resistances, not only to the drug itself but also to other drugs. The minimum inhibitory concentrations (MICs) of the antibiotics were investigated at 1×, 10×, 100× and 1000× concentrations using the adapted Microbial Evolution and Growth Arena (MEGA-plate) method. The results demonstrate that florfenicol can select for resistance to fluoroquinolone antibiotics (167× MIC value increase) and cephalosporins (67× MIC value increase). A total of 44 antimicrobial resistance genes were identified, the majority of which were consistent across the samples. Chromosomal point mutations, including alterations in resistance-associated and regulatory genes (acrB, acrR, emrR and robA), are thought to trigger multiple drug efflux pump activations, leading to phenotypically increased resistance. The study underscores the impact of florfenicol and its role in the development of antimicrobial resistance, particularly concerning fluoroquinolone antibiotics and cephalosporins. This study is the first to report florfenicol's dose-dependent enhancement of other antibiotics' MICs, linked to mutations in SOS-box genes (mdtABC-tolC, emrAB-tolC and acrAB-tolC) and increased multidrug efflux pump genes. Mutations in the regulatory genes acrR, emrR and rpbA support the possibility of increased gene expression. The results are crucial for understanding antimicrobial resistance and its development, highlighting the promising potential of in vitro evolutionary and coselection studies for future research.

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.

Genomic and transcriptomic insights into complex virus-prokaryote interactions in marine biofilms

Marine biofilms are complex communities of microorganisms that play a crucial ecological role in oceans. Although prokaryotes are the dominant members of these biofilms, little is known about their interactions with viruses. By analysing publicly available and newly sequenced metagenomic data, we identified 2446 virus-prokaryote connections in 84 marine biofilms. Most of these connections were between the bacteriophages in the Uroviricota phylum and the bacteria of Proteobacteria, Cyanobacteria and Bacteroidota. The network of virus-host pairs is complex; a single virus can infect multiple prokaryotic populations or a single prokaryote is susceptible to several viral populations. Analysis of genomes of paired prokaryotes and viruses revealed the presence of 425 putative auxiliary metabolic genes (AMGs), 239 viral genes related to restriction-modification (RM) systems and 38,538 prokaryotic anti-viral defence-related genes involved in 15 defence systems. Transcriptomic evidence from newly established biofilms revealed the expression of viral genes, including AMGs and RM, and prokaryotic defence systems, indicating the active interplay between viruses and prokaryotes. A comparison between biofilms and seawater showed that biofilm prokaryotes have more abundant defence genes than seawater prokaryotes, and the defence gene composition differs between biofilms and the surrounding seawater. Overall, our study unveiled active viruses in natural biofilms and their complex interplay with prokaryotes, which may result in the blooming of defence strategists in biofilms. The detachment of bloomed defence strategists may reduce the infectivity of viruses in seawater and result in the emergence of a novel role of marine biofilms.